RCFGL: Rapid Condition adaptive Fused Graphical Lasso and application to modeling brain region co-expression networks

Inferring gene co-expression networks is a useful process for understanding gene regulation and pathway activity. The networks are usually undirected graphs where genes are represented as nodes and an edge represents a significant co-expression relationship. When expression data of multiple (p) genes in multiple (K) conditions (e.g., treatments, tissues, strains) are available, joint estimation of networks harnessing shared information across them can significantly increase the power of analysis. In addition, examining condition-specific patterns of co-expression can provide insights into the underlying cellular processes activated in a particular condition. Condition adaptive fused graphical lasso (CFGL) is an existing method that incorporates condition specificity in a fused graphical lasso (FGL) model for estimating multiple co-expression networks. However, with computational complexity of O(p²K log K), the current implementation of CFGL is prohibitively slow even for a moderate number of genes and can only be used for a maximum of three conditions. In this paper, we propose a faster alternative of CFGL named rapid condition adaptive fused graphical lasso (RCFGL). In RCFGL, we incorporate the condition specificity into another popular model for joint network estimation, known as fused multiple graphical lasso (FMGL). We use a more efficient algorithm in the iterative steps compared to CFGL, enabling faster computation with complexity of O(p²K) and making it easily generalizable for more than three conditions. We also present a novel screening rule to determine if the full network estimation problem can be broken down into estimation of smaller disjoint sub-networks, thereby reducing the complexity further. We demonstrate the computational advantage and superior performance of our method compared to two non-condition adaptive methods, FGL and FMGL, and one condition adaptive method, CFGL in both simulation study and real data analysis. We used RCFGL to jointly estimate the gene co-expression networks in different brain regions (conditions) using a cohort of heterogeneous stock rats. We also provide an accommodating C and Python based package that implements RCFGL.     

Synaptophysin I selectively specifies the exocytic pathway of synaptobrevin 2/VAMP2

Biogenesis and recycling of synaptic vesicles are accompanied by sorting processes that preserve the molecular composition of the compartments involved. In the present study, we have addressed the targeting of synaptobrevin 2/VAMP2 (vesicle-associated membrane protein 2), a critical component of the synaptic vesicle--fusion machinery, in a heterotypic context where its sorting is not confounded by the presence of other neuron-specific molecules. Ectopically expressed synaptophysin I interacts with VAMP2 and alters its default surface targeting to a prominent vesicular distribution, with no effect on the targeting of other membrane proteins. Protein-protein interaction is not sufficient for the control of VAMP2 sorting, which is mediated by the C-terminal domain of synaptophysin I. Synaptophysin I directs the sorting of VAMP2 to vesicles before surface delivery, without influencing VAMP2 endocytosis. Consistent with this, dynamin and alpha-SNAP (soluble N-ethylmaleimide-sensitive fusion protein-attachment protein) mutants which block trafficking at the plasma membrane do not abrogate the effect of synaptophysin I on VAMP2 sorting. These results indicate that the sorting determinants of synaptic vesicle proteins can operate independently of a neuronal context and implicate the association of VAMP2 with synaptophysin I in the specification of the pathway of synaptic vesicle biogenesis.     

Bacterial colonization of a fumigated alkaline saline soil

After chloroform fumigating an arable soil, the relative abundance of phylotypes belonging to only two phyla (Actinobacteria and Firmicutes) and two orders [Actinomycetales and Bacillales (mostly Bacillus)] increased in a subsequent aerobic incubation, while it decreased for a wide range of bacterial groups. It remained to be seen if similar bacterial groups were affected when an extreme alkaline saline soil was fumigated. Soil with electrolytic conductivity between 139 and 157 dS m^?1, and pH 10.0 and 10.3 was fumigated and the bacterial community structure determined after 0, 1, 5 and 10 days by analysis of the 16S rRNA gene, while an unfumigated soil served as control. The relative abundance of the Firmicutes increased in the fumigated soil (52.8 %) compared to the unfumigated soil (34.2 %), while that of the Bacteroidetes decreased from 16.2 % in the unfumigated soil to 8.8 % in the fumigated soil. Fumigation increased the relative abundance of the genus Bacillus from 14.7 % in the unfumigated soil to 25.7 %. It was found that phylotypes belonging to the Firmicutes, mostly of the genus Bacillus, were dominant in colonizing the fumigated alkaline saline as found in the arable soil, while the relative abundance of a wide range of bacterial groups decreased.     

Rethinking the dispersal of Homo sapiens out of Africa

Current fossil, genetic, and archeological data indicate that Homo sapiens originated in Africa in the late Middle Pleistocene. By the end of the Late Pleistocene, our species was distributed across every continent except Antarctica, setting the foundations for the subsequent demographic and cultural changes of the Holocene. The intervening processes remain intensely debated and a key theme in hominin evolutionary studies. We review archeological, fossil, environmental, and genetic data to evaluate the current state of knowledge on the dispersal of Homo sapiens out of Africa. The emerging picture of the dispersal process suggests dynamic behavioral variability, complex interactions between populations, and an intricate genetic and cultural legacy. This evolutionary and historical complexity challenges simple narratives and suggests that hybrid models and the testing of explicit hypotheses are required to understand the expansion of Homo sapiens into Eurasia.     

Lack of alterations on meiotic chromosomes and morphological characteristics of male germ cells in mice exposed to a 60 Hz and 2.0 mT magnetic field

The effect of in vivo exposure of mice to a 60 Hz sinusoidal magnetic field (MF) at 2.0 mT on male germ cells was studied. The cytological endpoints measured included meiotic chromosome aberrations in spermatocytes and sperm morphology. Three independent experiments were carried out: (a) animals exposed for 72 h, (b) 10 days/8 h daily, and (c) 72 h exposure to MF plus 5 mg/kg of Mitomycin-C. No statistically significant differences indicative of MF effects were observed between MF exposed and control animals. In addition, an opposite effect between MF exposure and Mitomycin-C treatment in terms of chromosomal aberrations and sperm morphology was observed.     

Subtle mutational changes in the SU protein of a natural feline leukemia virus subgroup A isolate alter disease spectrum

FeLV-945 is a representative isolate of the natural feline leukemia virus (FeLV) variant predominant in non-T-cell malignant, proliferative, and degenerative diseases in a geographic cohort. The FeLV-945 surface glycoprotein (SU) is closely related to natural horizontally transmissible FeLV subgroup A (FeLV-A) but was found to differ from a prototype to a larger extent than the members of FeLV-A differ among themselves. The sequence differences included point mutations restricted largely to the functional domains of SU, i.e., VRA, VRB, and PRR. Despite the sequence differences in these critical domains, measurements of receptor utilization, including host range and superinfection interference, confirmed the assignment of FeLV-945 to subgroup A. Other proviruses isolated from the cohort contained similar sequence hallmarks and were assigned to FeLV subgroup A. A provirus from cat 1046 contained a histidine-to-proline change at SU residue 6 within an SPHQ motif that was previously identified as a critical mediator of fusion events during virus entry. The 1046 pseudotype virus entered cells only in the presence of the soluble cofactor FeLIX provided in trans, but it retained an ecotropic host range even in the presence of FeLIX. The mutational changes in FeLV-945 were shown to confer significant functional differences compared to prototype FeLV-A viruses. The substitution of FeLV-945 envelope gene sequences for FeLV-A/61E sequences conferred a small but statistically significant replicative advantage in some feline cells. Moreover, substitution of the unique FeLV-945 long terminal repeat and envelope gene for those of FeLV-A/61E altered the disease spectrum entirely, from a thymic lymphoma of a T-cell origin to an as yet uncharacterized multicentric lymphoma that did not contain T cells.     

Light- and singlet oxygen-mediated antifungal activity of phenylphenalenone phytoalexins.

The light-induced singlet oxygen production and antifungal activity of phenylphenalenone phytoalexins isolated from infected banana plants (Musa acuminata) are reported. Upon absorption of light energy all studied phenylphenalenones sensitise the production of singlet oxygen in polar and non-polar media. Antifungal activity of these compounds towards Fusarium oxysporum is enhanced in the presence of light. These results, together with the correlation of IC50 values under illumination with the quantum yield of singlet oxygen production and the enhancing effect of D2O on the antifungal activity, suggest the intermediacy of singlet oxygen produced by electronic excitation of the phenylphenalenone phytoalexins.     

Inflorescence architecture: A developmental genetics approach

We are characterizing a suiteof Pisum sativum mutants that alter inflorescence architecture to construct a model for the genetic regulation of inflorescence development in a plant with a compound raceme. Such a model, when compared with those created forAntirrhinum majus andArabidopsis thaliana, both of which have simple racemes, should provide insight into the evolution of the development of inflorescence architecture. The highly conserved nature of cloned genes that regulate reproductive development in plants and the morphological similarities among our mutants and those identified inA. majus andA. thaliana enhance the probability that a developmental genetics approach will be fruitful. Here we describe sixP. sativum mutants that affect morphologically and architecturally distinct aspects of the inflorescence, and we analyze interactions among these genes. Both vegetative and inflorescence growth of the primary axis is affected byUNIFOLIA TA, which is necessary for the function ofDETERMINATE (DET).DET maintains indeterminacy in the first-order axis. In its absence, the meristem differentiates as a stub covered with epidermal hairs.DET interacts withVEGETATIVE1 (VEG1).VEG1 appears essential for second-order inflorescence (I₂) development.veg1 mutants fail to flower or differentiate the I₂ meristem into a rudimentary stub,det veg1 double mutants produce true terminal flowers with no stubs, indicating that two genes must be eliminated for terminal flower formation inP. sativum, whereas elimination of a single gene accomplishes this inA. thaliana andA. majus. NEPTUNE also affects I₂ development by limiting to two the number of flowers produced prior to stub formation. Its role is independent ofDET, as indicated by the additive nature of the double mutantdet nep. UNI, BROC, and PIM all play roles in assigning floral meristem identity to the third-order branch.pim mutants continue to produce inflorescence branches, resulting in a highly complex architecture and aberrant flowers.uni mutants initiate a whorl of sepals, but floral organogenesis is aberrant beyond that developmental point, and the double mutantuni pim lacks identifiable floral organs. A wild-type phenotype is observed inbroc plants, butbroc enhancesthe pim phenotype in the double mutant, producing inflorescences that resemble broccoli. Collectively these genes ensure that only the third-order meristem, not higher- or lower-order meristems, generates floral organs, thus precisely regulating the overall architecture of the plant. Gene symbols used in this article: For clarity a common symbolization is used for genes of all species discussed in this article. Genes are symbolized with italicized capital letters. Mutant alleles are represented by lowercase, italicized letters. In both cases, the number immediately following the gene symbol differentiates among genes with the same symbol. If there are multiple alleles, a hyphen followed by a number is used to distinguish alleles. Protein products are represented by capital letters without italics.