Miracle-dz/UniFAR-FTUs huggingface dataset:
Effects of treatment with three antibiotics, vancomycin, neomycin, and AVNM on gut microbiome in C57BL/6 mice

Higher organisms, especially mammals, harbor diverse microbiota in the gut that plays a major role in maintaining health and physiological homeostasis. Perturbation of gut flora helps identifying their roles. Antibiotics are potent perturbing agents of microbiome. Select antibiotics like vancomycin, neomycin, and AVNM (an antibiotic cocktail containing ampicillin, vancomycin, neomycin, and metronidazole) were used to perturb the gut microbiota of C57BL/6 male mice to understand their roles in host immunity and metabolism. The current study revealed that the resulting gut microbial composition was different, and diversity (at the phylum and genus level) was reduced differentially following each antibiotic treatment. Vancomycin treatment caused a significant increase in Verrucomicrobia and Proteobacteria phyla. The treatment with neomycin yielded an increase in the Bacteroidetes phylum, while the treatment with AVNM led to an increase in Proteobacteria phylum with lowest diversity of microbiome in the gut. The current results also revealed that the different antibiotics treatment caused variation in the cecal index, expression of immune genes (TNF-α, IL-10, IFN-γ) in the colon, and short-chain fatty acids (SCFA) level in the blood of mice. A strong correlation was observed for antibiotic-induced differential dysbiosis patterns of gut microbiota and the altered immune and SCFA profile of the host. The outcome of the present study could be clinically important.
Co-culturing Hyphomicrobium nitrativorans strain NL23 and Methylophaga nitratireducenticrescens strain JAM1 allows sustainable denitrifying activities under marine conditions

Hyphomicrobium nitrativorans strain NL23 and Methylophaga nitratireducenticrescens strain JAM1 were the principal bacteria involved in the denitrifying activities of a methanol-fed, fluidized marine denitrification reactor. We believe that a tight relationship has developed between these two strains to achieve denitrification in the reactor under marine conditions. To characterize the potential synergy between strain JAM1 and strain NL23, we compared some of their physiological traits, and performed co-cultures. Pure cultures of strain JAM1 had a readiness to reduce nitrate (NO3−) with no lag phase for growth contrary to pure cultures of strain NL23, which has a 2-3 days lag phase before NO3− starts to be consumed and growth to occur. Compared to strain NL23, strain JAM1 has a higher μmax for growth and higher specific NO3− reduction rates. Antagonist assays showed no sign of exclusion by both strains. Planktonic co-cultures could only be performed on low NaCl concentrations for strain NL23 to survive. Denitrification rates were twice higher in the planktonic co-cultures than those measured in strain NL23 pure cultures. Biofilm co-cultures were performed for several months in a 500-mL bioreactor filled with Bioflow supports, and operated under fed-batch mode with increasing concentrations of NaCl for strain NL23 to acclimate to marine conditions. Under these conditions, the biofilm co-cultures showed sustained denitrifying activities and surface colonization by both strains. Increase in ectoine concentrations produced by strain JAM1 was observed in the biofilm with increasing NaCl concentrations. These results illustrate the capacity of both strains to act together in performing denitrification under marine environments. Although strain JAM1 did not contribute in better specific denitrifying activities in the biofilm co-cultures, its presence was essential for strain NL23 to survive in a medium with NaCl concentrations > 1.0%. We believe that ectoine is an important factor for the survival of strain NL23 in these environments.
Single cell transcriptomic analysis reveals cellular diversity of murine esophageal epithelium and age-associated mitochondrial dysfunction

Stratified squamous epithelium of the esophagus is comprised of basal keratinocytes that execute a terminal differentiation program in overlying suprabasal and superficial cell layers. Although morphologic progression coupled with expression of specific molecular markers has been characterized along the esophageal epithelial differentiation gradient, the molecular heterogeneity within the cell types along this trajectory has yet to be classified at the level of single cell resolution. To explore the molecular characteristics of esophageal keratinocytes along the squamous differentiation continuum, we performed single cell RNA-Sequencing transcriptomic profiling of 7,972 cells from murine esophageal epithelial sheets. We identified 8 distinct cell clusters in esophageal epithelium, unveiling an unexpected level of diversity, particularly among basal cells. We further mapped the cellular pathways and lineage trajectories within basal, suprabasal, and superficial clusters as well as within the heterogeneous basal cell populations, providing a comprehensive molecular view of esophageal epithelial cells in the context of squamous differentiation. Finally, we explored the impact of tissue aging upon esophageal epithelial cell clusters and demonstrated that mitochondrial dysfunction is a feature of aging in normal esophageal epithelium. These studies provide an unparalleled molecular perspective on murine esophageal keratinocytes that will serve as a valuable resource for dissecting cell type-specific roles in esophageal biology under conditions of homeostasis, aging, and tissue pathology.
Plasma Aβ ratios in autosomal dominant Alzheimer’s disease: the influence of genotype

In-vitro studies of autosomal dominant Alzheimer’s disease (ADAD) implicate longer Aβ peptides in pathogenesis, however less is known about the behaviour of ADAD mutations in-vivo. In this cross-sectional cohort study, we used liquid chromatography-tandem mass spectrometry to analyse 66 plasma samples from ADAD family members who were at-risk of inheriting a mutation or were already symptomatic. We tested for differences in plasma Aβ42:38, 38:40 and 42:40 ratios between Presenilin1 (PSEN1) and Amyloid Precursor Protein (APP) carriers. We examined the relationship between plasma and in-vitro models of Aβ processing and, among PSEN1 carriers, tested for associations with parental age at onset (AAO). 39 participants were mutation carriers (28 PSEN1 and 11 APP). Age- and sex-adjusted models showed marked differences in plasma Aβ between APP and PSEN1: higher Aβ42:38 in PSEN1 versus APP (p<0.001) and non-carriers (p<0.001); higher Aβ38:40 in APP versus PSEN1 (p<0.001) and non-carriers (p<0.001), while Aβ42:40 was higher in APP and PSEN1 compared to non-carriers (both p<0.001). Aβ profiles were reasonably consistent in plasma and cell lines. Within PSEN1, sex-adjusted models demonstrated negative associations between (i)Aβ42:40 (ii)Aβ42:38 and parental AAO. In-vivo differences in Aβ processing between APP and PSEN1 provide insights into ADAD pathophysiology which can inform therapy development.
Benchmarking Association Analyses of Continuous Exposures with RNA-seq in Observational Studies

Large datasets of hundreds to thousands of individuals measuring RNA-seq in observational studies are becoming available. Many popular software packages for analysis of RNA-seq data were constructed to study differences in expression signatures in an experimental design with well-defined conditions (exposures). In contrast, observational studies may have varying levels of confounding of the transcript-exposure associations; further, exposure measures may vary from discrete (exposed, yes/no) to continuous (levels of exposure), with non-normal distributions of exposure. We compare popular software for gene expression - DESeq2, edgeR, and limma - as well as linear regression-based analyses for studying the association of continuous exposures with RNA-seq. We developed a computation pipeline that includes transformation, filtering, and generation of empirical null distribution of association p-values, and we apply the pipeline to compute empirical p-values with multiple testing correction. We employ a resampling approach that allows for assessment of false positive detection across methods, power comparison, and the computation of quantile empirical p-values. The results suggest that linear regression methods are substantially faster with better control of false detections than other methods, even with the resampling method to compute empirical p-values. We provide the proposed pipeline with fast algorithms in R.
Struo2: efficient metagenome profiling database construction for ever-expanding microbial genome datasets

Mapping metagenome reads to reference databases is the standard approach for assessing microbial taxonomic and functional diversity from metagenomic data. However, public reference databases often lack recently generated genomic data such as metagenome-assembled genomes (MAGs), which can limit the sensitivity of read-mapping approaches. We previously developed the Struo pipeline in order to provide a straight-forward method for constructing custom databases; however, the pipeline does not scale well with the ever-increasing number of publicly available microbial genomes. Moreover, the pipeline does not allow for efficient database updating as new data are generated. To address these issues, we developed Struo2, which is >3.5-fold faster than Struo at database generation and can also efficiently update existing databases. We also provide custom Kraken2, Bracken, and HUMAnN3 databases that can be easily updated with new genomes and/or individual gene sequences. Struo2 enables feasible database generation for continually increasing large-scale genomic datasets. Availability Struo2: https://github.com/leylabmpi/Struo2 Pre-built databases: http://ftp.tue.mpg.de/ebio/projects/struo2/ Utility tools: https://github.com/nick-youngblut/gtdb_to_taxdump
Synaptotagmin 7 is enriched at the plasma membrane through γ-secretase processing to promote vesicle docking and control synaptic plasticity in mouse hippocampal neurons

Synaptotagmin (SYT) 7 has emerged as key regulator of presynaptic function, but its localization and precise function in the synaptic vesicle cycle remain unclear. Here, we used iGluSnFR to optically and directly interrogate glutamate release, at the single bouton level, in SYT7 KO dissociated mouse hippocampal neurons. We analyzed asynchronous release, paired pulse facilitation, and synaptic vesicle replenishment, and found that SYT7 contributes to each of these processes to different degrees. ‘Zap-and-freeze’ electron microscopy revealed that loss of SYT7 impairs the docking of synaptic vesicles after a stimulus and the recovery of depleted synaptic vesicles after a stimulus train. To execute these functions, SYT7 must be targeted to the plasma membrane via γ-secretase-mediated cleavage of the amino terminus, followed by palmitoylation. The complex sorting itinerary of SYT7 endows this Ca2+-sensor with the ability to control crucial forms of synaptic function and plasticity. SYT7 mediated asynchronous release, paired pulse facilitation, and synaptic vesicle replenishment was observed optically at individual hippocampal synapses Localization, trafficking, and stability of SYT7 is dependent on processing by γ-secretase Short term plasticity defects arise in SYT7KOs due to decreased docking of synaptic vesicles after stimulation SYT7 promotes paired-pulse facilitation and asynchronous release via distinct mechanisms
Neuromesodermal progenitor origin of trunk neural crest in vivo

Neural crest (NC) is a vertebrate-specific population of multipotent embryonic cells predisposed to diverse derivatives along the anteroposterior (A-P) axis. Only cranial NC progenitors give rise to ectomesenchymal cell types, whereas trunk NC is biased for neuronal cell fates. By integrating multimodal single-cell analysis, we provide evidence for divergent embryonic origins of cranial vs. trunk NC that explain this dichotomy. We show that the NC regulator foxd3 is heterogeneously expressed across the A-P axis and identify its specific cranial and trunk autoregulatory enhancers. Whereas cranial-specific enhancer is active in the bona fide NC, the trunk foxd3 autoregulatory element surprisingly marked bipotent tailbud neuromesodermal progenitors (NMps). We integrated NMp single cell epigemomics and trasncriptomics data and for the first time reconstructed anamniote NMp gene regulatory network. Moreover, using pseudotime and developmental trajectory analyses of NMps and NC during normal development and in foxd3 mutants, we demonstrate an active role for foxd3 in balancing non-cranial NC and NMp fates during early embryonic development. Strikingly, we show that a portion of posterior NC in the developing zebrafish embryo is derived from the pro-neural NMps. This suggests a common embryonic origin of trunk NC and NM progenitors that is distinct from cranial NC anlage, and elucidates pro-neural bias of trunk NC.
Specification of axial identity by Hoxa2 distinguishes between a phenotypic and molecular ground state in mouse cranial neural crest cells

Hox genes play a key role in head formation by specifying the axial identity of neural crest cells (NCCs) migrating into embryonic pharyngeal arches. In the absence of Hoxa2, NCC derivatives of the second pharyngeal arch (PA2) undergo a homeotic transformation and duplicate structures formed by first arch (PA1) NCCs. Current models postulate that PA1 represents a NCC ‘ground state’ and loss of Hoxa2 causes a reversion of PA2 NCCs to the PA1 ‘ground state’. We use bulk and single-cell RNAseq to investigate the molecular mechanisms driving this phenotypic transformation in the mouse. In Hoxa2-/- mutants, PA2 NCCs generally maintain expression of the PA2 transcriptional signature and fail to strongly upregulate a PA1 transcriptional signature. Our analyses identify putative HOXA2 targets and suggest that subsets of NCCs may respond to HOXA2 activity in distinct manners. This separation of phenotypic and molecular states has significant implications for understanding craniofacial development.
An algal symbiont (Breviolum psygmophilum) responds more strongly to chronic high temperatures than its facultatively symbiotic coral host (Astrangia poculata)

Scleractinian corals form the foundation of coral reefs by secreting skeletons of calcium carbonate. Their intracellular algal symbionts (Symbiodiniaceae) translocate a large proportion of photosynthate to the coral host, which is required to maintain high rates of calcification. Global warming is causing dissociation of coral host and algal symbiont, visibly presented as coral bleaching. Despite decades of study, the precise mechanisms of coral bleaching remain unknown. Separating the thermal stress response of the coral from the algal symbiont is key to understanding bleaching in tropical corals. The facultatively symbiotic northern star coral, Astrangia poculata, naturally occurs as both symbiotic and aposymbiotic (lacking algal symbionts) polyps – sometimes on the same coral colony. Thus, it is possible to separate the heat stress response of the coral host alone from the coral in symbiosis with its symbiont Breviolum psygmophilum. Using replicate symbiotic and aposymbiotic ramets of A. poculata, we conducted a chronic heat stress experiment to increase our understanding of the cellular mechanisms resulting in coral bleaching. Sustained high temperature stress resulted in photosynthetic dysfunction in B. psygmophilum, including a decline in maximum photosynthesis rate, maximum photochemical efficiency, and the absorbance peak of chlorophyll a. Interestingly, the metabolic rates of symbiotic and aposymbiotic corals were differentially impacted. RNAseq analysis revealed more differentially expressed genes between heat-stressed and control aposymbiotic colonies than heat-stressed and control symbiotic colonies. Notably, aposymbiotic colonies increased the expression of inflammation-associated genes such as nitric oxide synthases. Unexpectedly, the largest transcriptional response was observed between heat-stressed and control B. psygmophilum, including genes involved in photosynthesis, response to oxidative stress, and meiosis. Thus, it appears that the algal symbiont suppresses the immune response of the host, potentially increasing the vulnerability of the host to pathogens. The A. poculata-B. psygmophilum symbiosis provides a tractable model system for investigating thermal stress and immune challenge in scleractinian corals.
Extensive N4 Cytosine Methylation is Essential for Marchantia Sperm Function

4-methylcytosine (4mC) is an important DNA modification in prokaryotes, but its relevance, and even presence in eukaryotes have been mysterious. Here we show that spermatogenesis in the liverwort Marchantia polymorpha involves two waves of extensive DNA methylation reprogramming. First, 5-methylcytosine (5mC), a well-known eukaryotic DNA modification, expands from transposons to the entire genome. Notably, the second wave installs 4mC throughout genic regions, covering over 50% of CG sites in sperm. 4mC is catalyzed by two novel methyltransferases (MpDN4MT1a and MpDN4MT1b) specifically expressed during late spermiogenesis. Deletion of MpDN4MT1s eliminates 4mC, alters the sperm transcriptome, and produces sperm with swimming defects. Our results reveal extensive 4mC in a eukaryote and define a new family of eukaryotic methyltransferases, thereby expanding the repertoire of functional eukaryotic DNA modifications.
A macroecological description of alternative stable states reproduces intra- and inter-host variability of gut microbiome

The most fundamental questions in microbial ecology concern the diversity and variability of communities. Their composition varies widely across space and time, as it is determined by a non-trivial combination of stochastic and deterministic processes. The interplay between non-linear community dynamics and environmental fluctuations determines the rich statistical structure of community variability, with both rapid temporal dynamics fluctuations and non-trivial correlations across habitats. Here we analyze long time-series of gut microbiome and compare intra- and inter-community dissimilarity. Under a macroecological framework we characterize their statistical properties. We show that most taxa have large but stationary fluctuations over time, while a minority is characterized by quick changes of average abundance which cluster in time, suggesting the presence of alternative stable states. We disentangle inter-individual variability in a major stochastic component and a deterministic one, the latter recapitulated by differences in the carrying capacities of taxa. Finally, we develop a model which includes environmental fluctuations and alternative stable states. This model quantitatively predicts the statistical properties of both intra- and inter-individual community variability, therefore summarizing variation in a unique macroecological framework.
Argonaute binding within human nuclear RNA and its impact on alternative splicing

Mammalian RNA interference (RNAi) is often linked to the regulation of gene expression in the cytoplasm. Synthetic RNAs, however, can also act through the RNAi pathway to regulate transcription and splicing. While nuclear regulation by synthetic RNAs can be robust, a critical unanswered question is whether endogenous functions for nuclear RNAi exist in mammalian cells. Using enhanced crosslinking immunoprecipitation (eCLIP) in combination with RNA sequencing (RNAseq) and multiple AGO knockout cell lines, we mapped AGO2 protein binding sites within nuclear RNA. The strongest AGO2 binding sites were mapped to micro RNAs (miRNAs). The most abundant miRNAs were distributed similarly between cytoplasm and nucleus, providing no evidence for mechanisms that facilitate localization of miRNAs in one compartment versus the other. Beyond miRNAs, most statistically-significant AGO2 binding was within introns. Splicing changes were confirmed by RT-PCR and were recapitulated by synthetic duplex RNAs and miRNA mimics complementary to the sites of AGO2 binding. These data support the hypothesis that miRNAs can control gene splicing. While nuclear RNAi proteins have the potential to be a natural regulatory mechanism, careful study will be necessary to identify critical RNA drivers of normal physiology and disease.
Disordered information processing dynamics in experimental epilepsy

Neurological disorders share common high-level alterations, such as cognitive deficits, anxiety, and depression. This raises the possibility of fundamental alterations in the way information conveyed by neural firing is maintained and dispatched in the diseased brain. Using experimental epilepsy as a model of neurological disorder we tested the hypothesis of altered information processing, analyzing how neurons in the hippocampus and the entorhinal cortex store and exchange information during slow and theta oscillations. We equate the storage and sharing of information to low level, or primitive, information processing at the algorithmic level, the theoretical intermediate level between structure and function. We find that these low-level processes are organized into substates during brain states marked by theta and slow oscillations. Their internal composition and organization through time are disrupted in epilepsy, loosing brain state-specificity, and shifting towards a regime of disorder in a brain region dependent manner. We propose that the alteration of information processing at an algorithmic level may be a mechanism behind the emergent and widespread co-morbidities associated with epilepsy, and perhaps other disorders.
Quantitative Proteomics Identifies Secreted Diagnostic Biomarkers as well as Tumor-Dependent Prognostic Targets for Clear Cell Renal Cell Carcinoma

Clear cell renal cell carcinoma (ccRCC) is the third most common and most malignant urological cancer, with a 5-year survival rate of 10% for patients with advanced tumors. Here, we identified 10,160 unique proteins by in-depth quantitative proteomics, of which 955 proteins were significantly regulated between tumor and normal adjacent tissues. We verified four putatively secreted biomarker candidates, namely, PLOD2, FERMT3, SPARC, and SIRPα, as highly expressed proteins that are not affected by intratumor and intertumor heterogeneity. Moreover, SPARC displayed a significant increase in urine samples of patients with ccRCC, making it a promising marker for the detection of the disease in body fluids. Furthermore, based on molecular expression profiles, we propose a biomarker panel for the robust classification of ccRCC tumors into two main clusters, which significantly differed in patient outcome with an almost three times higher risk of death for cluster 1 tumors compared with cluster 2 tumors. Moreover, among the most significant clustering proteins, 13 were targets of repurposed inhibitory FDA-approved drugs. Our rigorous proteomics approach identified promising diagnostic and tumor-discriminative biomarker candidates which can serve as therapeutic targets for the treatment of ccRCC. Implications: Our in-depth quantitative proteomics analysis of ccRCC tissues identifies the putatively secreted protein SPARC as a promising urine biomarker and reveals two molecular tumor phenotypes.