Miracle-dz/UniFAR-FTUs huggingface dataset:
A pH-eQTL interaction at the RIT2-SYT4 Parkinson’s disease risk locus in the substantia nigra

Parkinson’s disease (PD) causes severe motor and cognitive disabilities that result from the progressive loss of dopamine neurons in the substantia nigra. The rs12456492 variant in the RIT2 gene has been repeatedly associated with increased risk for Parkinson’s disease. From a transcriptomic perspective, a meta-analysis found that RIT2 gene expression is correlated with pH in the human brain. To assess these pH associations in relation to PD risk, we examined the two datasets that assayed rs12456492, gene expression, and pH in the postmortem human brain. Using the BrainEAC dataset, we replicate the positive correlation between RIT2 gene expression and pH in the human brain (n=100). Furthermore, we found that the relationship between expression and pH is influenced by rs12456492. When tested across ten brain regions, this interaction is specifically found in the substantia nigra. A similar association was found for the co-localized SYT4 gene. In addition, SYT4 associations are stronger in a combined model with both genes, and the SYT4 interaction appears to be specific to males. In the GTEx dataset, the pH associations involving rs12456492 and expression of either SYT4 and RIT2 was not seen. This null finding may be due to the short postmortem intervals (PMI) of the GTEx tissue samples. In the BrainEAC data, we tested the effect of PMI and only observed the interactions in the longer PMI samples. These previously unknown associations suggest novel mechanistic roles for rs12456492, RIT2, and SYT4 in the regulation of pH in the substantia nigra.
Cell patterning by secretion-induced plasma membrane flows

Cells self-organize using reaction-diffusion and fluid-flow principles. Whether bulk membrane flows contribute to cell patterning has not been established. Here, using mathematical modelling, optogenetics and synthetic probes, we show that polarized exocytosis causes lateral membrane flows away from regions of membrane insertion. Plasma membrane-associated proteins with sufficiently low diffusion and/or detachment rates couple to the flows and deplete from areas of exocytosis. In rod-shaped fission yeast cells, zones of Cdc42 GTPase activity driving polarized exocytosis are limited by GTPase activating proteins (GAPs). We show that membrane flows pattern the GAP Rga4 distribution and coupling of a synthetic GAP to membrane flows is sufficient to establish the rod shape. Thus, membrane flows induced by Cdc42-dependent exocytosis form a negative feedback restricting the zone of Cdc42 activity. One Sentence Summary Exocytosis causes bulk membrane flows that drag associated proteins and form a negative feedback restricting the exocytic site.
Longitudinal profiling of the macaque vaginal microbiome reveals similarities to diverse human vaginal communities: implications for use as a pre-clinical model for bacterial vaginosis

The vaginal microbiota plays an important role in women’s reproductive and urogenital health. Disturbances in this microbial community can lead to several adverse outcomes including pelvic inflammatory disease, bacterial vaginosis (BV) as well as increased susceptibility to sexually transmitted infections, miscarriage, and pre-term births. It is now well accepted that while the microbiome of healthy women in the developed world is dominated by Lactobacillus species, vaginal communities in asymptomatic women, especially those in the developing world, can be comprised of a diverse set of micro-organisms. The presence of a diverse vaginal microbiome has been associated with increased susceptibility to HIV infection but their implications for women’s health remain poorly understood. Rhesus macaques are an excellent translational animal model due to significant physiological and genetic homology with humans. In this study, we performed a longitudinal analysis of clinical and microbiome data from 16 reproductive age female rhesus macaques. Many animals showed hallmarks of BV, including Nugent scores above 7 and high vaginal pH. At both the taxonomic and functional level, the rhesus macaque vaginal microbiome was most similar to that of women who harbor a diverse vaginal community associated with asymptomatic/symptomatic bacterial vaginosis. Specifically, rhesus macaque vaginal microbiomes harbored a diverse set of anaerobic gram-negative bacteria, including; Snethia, Prevotella, Porphyromonas, and Mobilluncus. Interestingly, some animals were transiently colonized by Lactobacillus and some with Gardnerella. Our in-depth and comprehensive analysis highlights the importance of the model to test interventions for manipulating the vaginal microbiome. IMPORTANCE It is widely accepted that the “healthy” vaginal microbiome of the majority of women in the developed world is dominated by Lactobacillus species. However, in the developing world, a majority of women are colonized by diverse microbial communities, typically associated with bacterial vaginosis, but remain asymptomatic. Many questions remain about the drivers of this disparity and potential interventions to alter the vaginal microbiome. Rhesus macaques provide an excellent translational model due to significant physiological and genetic homology with humans. In this study, we performed a longitudinal analysis of clinical and microbiome data from a large cohort of reproductive age rhesus macaques. At the taxonomic, genomic, and functional level, the rhesus macaque vaginal microbiome was most similar to that of humans who harbor a diverse vaginal community associated with asymptomatic/symptomatic bacterial vaginosis. Our in-depth and comprehensive analysis highlights the utility of macaques to test interventions for manipulating the vaginal microbiome.
Astrocytes amplify cerebral blood flow elevation to sustained cortical activation in the awake mouse

Brain metabolism requires increased local cerebral blood flow (CBF) during neuronal activation to match O2 supply with demand. This phenomenon, termed functional hyperemia, is also the basis of functional neuroimaging. Ca2+ elevation in astrocytes can drive arteriole dilation to increase CBF, yet affirmative necessity and sufficiency evidence under realistic conditions in vivo for astrocytes in functional hyperemia is lacking. In awake mice we show that functional hyperemia is bi-phasic and arteriole dilation progresses as sensory stimulation is sustained. Reduction of astrocyte Ca2+ signaling in vivo reduces sustained but not brief sensory-evoked arteriole dilation. Chemogenetic activation of astrocytes selectively augments sustained hyperemia. NMDA-receptors and epoxyeicosatrienoic acid explain the astrocyte-mediated effects on the late, escalating phase of functional hyperemia but have no impact on brief blood flow events. We propose that a fundamental role of astrocyte Ca2+ is to amplify functional hyperemia when neuronal activation is prolonged.
Anatomically informed multi-level fiber tractography

Diffusion weighted MR imaging can assist preoperative planning by reconstructing the trajectory of eloquent fiber pathways. A common task is the delineation of the corticospinal tract in its full extent because lesions to this bundle can severely affect the quality of life. However, this is challenging as existing tractography algorithms typically produce either incomplete results or multiple false-positive tracts. In this work, we suggest a novel approach to fiber tractography that reconstructs multi-level structures by progressively taking into account previously unused fiber orientations. Anatomical priors are used in order to minimize the number of false-positive pathways. The devised method was evaluated on synthetic data with different noise levels. Additionally, it was tested on in-vivo data by reconstructing the corticospinal tract and it was compared to conventional deterministic and probabilistic approaches. The corticospinal tract reconstructed by our method includes lateral projections that could not be observed with deterministic methods, while avoiding spurious tracts reconstructed by probabilistic tractography. Furthermore, the proposed algorithm preserves the neuroanatomical topology of the pathways to a larger extent as compared to probabilistic tractography.
Response of benthic macroinvertebrates to dam removal in the restoration of the Boardman River, Michigan, USA

Dam removal is an increasingly important method of stream restoration, but most removal efforts are under-studied in their effects. In order to better understand the effects of such removals on the stream ecosystem, we examined changes in stream macroinvertebrate communities from 2011-2016 above, below, and before and after the October 2012 removal of the Brown Bridge Dam on the Boardman River in Michigan (USA), and to new channel sites created in its former reservoir (2013-2015). Using linear mixed-effect models on the percent abundance of ecologically sensitive taxa (% Ephemeroptera, Plecoptera, Trichoptera (EPT)), total density of all macroinvertebrates, and overall taxa richness, along with multivariate analyses on the community matrix, we examined differences in community composition among sites and years. EPT declined downstream of the dam immediately after dam removal, but recovered in the second year, becoming dominant within 2-4 years. Downstream sites before removal had different community composition than upstream sites and downstream sites after removal (p<0.001), while upstream and downstream sites after removal converged towards similarity. New channel (restored) %EPT, density, and taxa richness were not different from upstream sites in any year following removal, but new channel sites were the most distinct in community composition, possessing multiple indicator taxa characteristic of unique new conditions. The invasive New Zealand mud snail (Potamopyrgus antipodarum) was absent from all sites prior to dam removal, but appeared at low densities in upstream sites in 2013, had spread to all sites by 2015, and showed large increases at all sites by 2016. Managers employing dam removal for stream restoration should include post-removal monitoring for multiple years following removal and conduct risk analysis regarding potential effects on colonization of invasive invertebrate species.
Dynamic methylome modification is associated with mutational signatures in aging and the etiology of disease

Reversible epigenetic changes within the loci of genes that regulate critical cell processes have recently emerged as important biomarkers of disease pathology. It is then natural to consider the consequences for population health risk of such epigenetic changes during the aging process. Specifically, the interplay between dynamic methylation changes that accompany aging and mutations that accrue in an individual’s genome over time need further investigation. The current study investigated the role of dynamic methylation acting together with gene variants in an individual over time to gain insight into the evolving epigenome–genome interplay that affects biochemical pathways controlling physiological processes during aging. We completed a whole-genome methylation and variant analysis in a non-smoking Zoroastrian-Parsi individual, collecting two samples, 12 years apart (at 53 and 65 years respectively) (ZPMetG-Hv2a-1A (old, t0), ZPMetG-Hv2a-1B (recent, t0+12)) and analysing them using a GridION Nanopore sequencer at 13X genome coverage overall. We further identified the single nucleotide variants (SNVs) and indels in known CpG islands by employing the Genome Analysis Tool Kit (GATK) and MuTect2 variant-caller pipeline with the GRCh37 (patch 13) human genome as reference. We found 5258 disease-relevant genes that had been differentially methylated in this individual over 12 years. Employing the GATK pipeline, we found 24,948 genes, corresponding to 4,58,148 variants, specific to ZPMetG-Hv2a-1B, indicating the presence of variants that had accrued over time. A fraction of the gene variants (242/24948) occurred within the CpG regions that were differentially methylated, with 67/247 exactly coincident with a CpG site. Our analysis yielded a critical cluster of 10 genes that were each significantly methylated and had variants at the CpG site or the ±4 bp CpG region window. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment network analysis as well as Reactome and STRING analysis of gene-specific variants indicated an impact on biological processes regulating the immune system, disease networks implicated in cancer and neurodegenerative diseases, and transcriptional control of processes regulating cellular senescence and longevity. Additional analysis of mutational signatures indicated a majority of C>T transitions followed by T>C transitions in the more recent sample, ZPMetG-Hv2a-1B. Our current study provides additional insight into the aging methylome over time and the interplay between different methylation and gene variants in the etiology of disease.
Methane-derived carbon flow through host-virus trophic networks in soil

The concentration of atmospheric methane continues to increase with microbial communities controlling soil-atmosphere fluxes. While there is substantial knowledge of the diversity and function of organisms regulating methane production and consumption, the frequency and impact of interactions with viruses on their activity in soil is unknown. Metagenomic sequencing of soil microbial communities has enabled identification of linkages between viruses and hosts. However, determining host-virus linkages through sequencing does not determine whether a virus or a host are active. In this study, we identified active individual interactions in situ by following the transfer of assimilated carbon from active hosts to viruses. Using DNA stable-isotope probing combined with metagenomic analyses, we characterized methane-fueled microbial networks in acidic and neutral pH soils, specifically primary and secondary utilisers of carbon, together with the recent transfer of methane-derived carbon to viruses. Sixty-three percent of viral contigs from replicated soil incubations contained genes associated with known methanotrophic bacteria. Genomic sequences from 13C-enriched viruses were present in clustered regularly interspaced short palindromic repeats (CRISPR) arrays of multiple, closely-related Methylocystis populations, revealing differences in their history of viral interaction. Viruses infecting non-methanotrophic methylotrophs and heterotrophic predatory bacteria were also identified through the analysis of shared homologous genes, demonstrating that carbon is transferred to a diverse range of viruses associated with methane-fueled microbial food networks.
Proteomic identification of the UDP-GlcNAc: PI α1–6 GlcNAc-transferase subunits of the glycosylphosphatidylinositol biosynthetic pathway of Trypanosoma brucei

The first step of glycosylphosphatidylinositol (GPI) anchor biosynthesis in all eukaryotes is the addition of N-acetylglucosamine (GlcNAc) to phosphatidylinositol (PI) which is catalysed by a UDP-GlcNAc : PI α1-6 GlcNAc-transferase. This enzyme has been shown to be a complex of at least seven subunits in mammalian cells and a similar complex of homologous subunits has been postulated in yeast. Homologs of most of these mammalian and yeast subunits were identified in the Trypanosoma brucei predicted protein database. The putative catalytic subunit of the T. brucei complex, TbGPI3, was epitope tagged with three consecutive c-Myc sequences at its C-terminus. Immunoprecipitation of TbGPI3-3Myc followed by native polyacrylamide gel electrophoresis and anti-Myc Western blot showed that it is present in a ~240 kDa complex. Label-free quantitative proteomics were performed to compare anti-Myc pull-downs from lysates of TbGPI-3Myc expressing and wild type cell lines. TbGPI3-3Myc was the most highly enriched protein in the TbGPI3-3Myc lysate pull-down and partner proteins TbGPI15, TbGPI9, TbGPI2, TbGPI1 and TbERI1 were also identified with significant enrichment. Our proteomics data also suggest that an Arv1-like protein (TbArv1) is a subunit of the T. brucei complex. Yeast and mammalian Arv1 have been previously implicated in GPI biosynthesis, but here we present the first experimental evidence for physical association of Arv1 with GPI biosynthetic machinery. A putative E2-ligase has also been tentatively identified as part of the T. brucei UDP-GlcNAc : PI α1-6 GlcNAc-transferase complex. Graphical abstract First step of GPI anchor biosynthesis pathway in T.brucei BSF is catalysed by TbGPI3 complex.
The molecular framework of heterophylly in Callitriche palustris L. differs from that in other amphibious plants

Heterophylly refers to the development of different leaf forms in a single plant depending on the environmental conditions. It is often observed in amphibious aquatic plants that can grow under aerial and submerged conditions. Although heterophylly is well recognized in aquatic plants, the associated developmental mechanisms and the molecular basis remain unclear. In this study, we analyzed heterophyllous leaf formation in an aquatic plant, Callitriche palustris, to clarify the underlying developmental and molecular mechanisms. Morphological analyses revealed extensive cell elongation and the rearrangement of cortical microtubules in the elongated submerged leaves of C. palustris. Our observations also suggested that gibberellin, ethylene, and abscisic acid regulate the formation of submerged leaves. However, the perturbation of one or more of the hormones was insufficient to induce the formation of submerged leaves under aerial conditions. Finally, we analyzed gene expression changes during aerial and submerged leaf development and narrowed down the candidate genes controlling heterophylly via transcriptomic comparisons, including a comparison with a closely related terrestrial species. We revealed that the molecular mechanism regulating heterophylly in C. palustris is associated with hormonal changes and diverse transcription factor gene expression profiles, suggesting differences from the corresponding mechanisms in previously investigated amphibious plants.
Integration of Cell Growth and Asymmetric Division during Lateral Root Initiation in Arabidopsis thaliana

Lateral root formation determines to a large extent the ability of plants to forage their environment and thus their growth. In Arabidopsis thaliana and other angiosperms, lateral root initiation requires radial cell expansion and several rounds of anticlinal cell divisions that give rise to a central core of small pericycle cells, which express different markers than the larger surrounding cells. These small central cells then switch their plane of divisions to periclinal, and give rise to seemingly morphologically similar daughter cells that have different identities and establish the different cell types of the new root. Although the execution of these two types of divisions is tightly regulated and essential for the correct development of the lateral root, we know little about their geometrical features. Here we analyse a four-dimensional reconstruction of the first stages of lateral root formation and analyze the geometric features of the anticlinal and periclinal divisions. We identify that the periclinal divisions of the small central cells are morphologically dissimilar and asymmetric. We show that mother cell volume is different when looking at anticlinal versus periclinal divisions and the repeated anticlinal divisions do not lead to reduction in cell volume although cells are shorter. Finally, we show that cells undergoing a periclinal division are characterized by a strong cell expansion. Our results indicate that cells integrate growth and division to precisely partition their volume upon division during the first two stages of lateral root formation.
The loss of IL-31 signaling attenuates bleomycin-induced pulmonary fibrosis

Idiopathic Pulmonary Fibrosis (IPF) is a severe fibrotic lung disease characterized by excessive collagen deposition and progressive decline in lung function. Multiple Th2 T cell-derived cytokines including IL-4 and IL-13 have been shown to contribute to inflammation and fibrotic remodeling in multiple tissues. Interleukin-31 (IL-31) is a newly identified cytokine that is predominantly produced by CD4 TH2 T cells, but its signaling receptor called IL-31RA has been shown predominately expressed by non-hematopoietic cells. However, the potential role of the IL-31-IL31RA axis in pulmonary inflammation and fibrosis has remained largely unknown. To determine the role of IL-31 signaling in pulmonary fibrosis, wildtype, and IL-31RA knockout mice were treated with bleomycin and measured changes in collagen deposition and lung function. Notably, the loss of IL-31 signaling attenuated collagen deposition and lung function decline during bleomycin-induced pulmonary fibrosis. However, the loss of IL-31RA signaling did not affect inflammation in the lungs. The total lung transcriptome analysis showed a significant reduction in fibrosis-associated gene transcripts including ECM- and epithelial cell-associated gene networks. Furthermore, the lungs of IPF showed an elevated expression of IL-31 when compared to control subjects. In support, the percentage of IL-31 producing CD4+ T cells was greater in the lungs and PBMCs from IPF patients compared to healthy controls. Our findings suggest a pathogenic role for IL-31/IL-31RA signaling during bleomycin-induced pulmonary fibrosis. In summary, therapeutic targeting of the IL-31-IL-31RA axis could be beneficial in pulmonary fibrosis and has translational benefits specifically by preventing collagen deposition and improving lung function.
CTCF is a Barrier for Totipotent-like Reprogramming

Totipotent cells have the ability of generating embryonic and extra-embryonic tissues1,2. Interestingly, a rare population of cells with totipotent-like potential was identified within ESC cultures3. These cells, known as 2 cell (2C)-like cells, arise from ESC and display similar features to those found in the totipotent 2 cell embryo2-4. However, the molecular determinants of 2C-like conversion have not been completely elucidated. Here, we show that CTCF is a barrier for 2C-like reprogramming. Indeed, forced conversion to a 2C-like state by DUX expression was associated with DNA damage at a subset of CTCF binding sites. Endogenous or DUX-induced 2C-like ESC showed decreased CTCF enrichment at known binding sites, suggesting that acquisition of a totipotent-like state is associated with a highly dynamic chromatin architecture. Accordingly, depletion of CTCF in ESC efficiently promoted spontaneous and asynchronous conversion to a totipotent-like state. This phenotypic reprogramming was reversible upon restoration of CTCF levels. Furthermore, we showed that transcriptional activation of the ZSCAN4 cluster was necessary for successful 2C-like reprogramming. In summary, we revealed the intimate relation between CTCF and totipotent-like reprogramming.
Decision tree models and cell fate choice

Single cell transcriptomics has laid bare the heterogeneity of apparently identical cells at the level of gene expression. For many cell-types we now know that there is variability in the abundance of many transcripts, and that average transcript abun-dance or average gene expression can be a unhelpful concept. A range of clustering and other classification methods have been proposed which use the signal in single cell data to classify, that is assign cell types, to cells based on their transcriptomic states. In many cases, however, we would like to have not just a classifier, but also a set of interpretable rules by which this classification occurs. Here we develop and demonstrate the interpretive power of one such approach, which sets out to establish a biologically interpretable classification scheme. In particular we are interested in capturing the chain of regulatory events that drive cell-fate decision making across a lineage tree or lineage sequence. We find that suitably defined decision trees can help to resolve gene regulatory programs involved in shaping lineage trees. Our approach combines predictive power with interpretabilty and can extract logical rules from single cell data.