People are exposed to phthalates through their wide use as plasticizers and in personal care products. Many phthalates are endocrine disruptors and have been associated with adverse health outcomes. However, knowledge gaps exist in understanding the molecular mechanisms associated with the effects of exposure in early and late pregnancy. In this study, we examined the relationship of eleven urinary phthalate metabolites with isoprostane, an established marker of oxidative stress, among pregnant Mexican-American women from an agricultural cohort. Isoprostane levels were on average 20% higher at 26 weeks than at 13 weeks of pregnancy. Urinary phthalate metabolite concentrations suggested relatively consistent phthalate exposures over pregnancy. The relationship between phthalate metabolite concentrations and isoprostane levels was significant for the sum of di-2-ethylhexyl phthalate and the sum of high molecular weight metabolites with the exception of monobenzyl phthalate, which was not associated with oxidative stress at either time point. In contrast, low molecular weight metabolite concentrations were not associated with isoprostane at 13 weeks, but this relationship became stronger later in pregnancy (p-value = 0.009 for the sum of low molecular weight metabolites). Our findings suggest that prenatal exposure to phthalates may influence oxidative stress, which is consistent with their relationship with obesity and other adverse health outcomes.

Recent genome- and epigenome-wide studies demonstrate that the DNA methylation is controlled in part by genetics, highlighting the importance of integrating genetic and epigenetic data. To better understand molecular mechanisms affecting gene expression, we used the candidate susceptibility gene paraoxonase 1 (PON1) as a model to assess associations of PON1 genetic polymorphisms with DNA methylation and arylesterase activity, a marker of PON1 expression. PON1 has been associated with susceptibility to obesity, cardiovascular disease, and pesticide exposure. In this study, we assessed DNA methylation in 18 CpG sites located along PON1 shores, shelves, and its CpG island in blood specimens collected from newborns and 9-year-old children participating (n = 449) in the CHAMACOS birth cohort study. The promoter polymorphism, PON1-108 , was strongly associated with methylation, particularly for CpG sites located near the CpG island (P << 0.0005). Among newborns, these relationships were even more pronounced after adjusting for blood cell composition. We also observed significant decreases in arylesterase activity with increased methylation at the same nine CpG sites at both ages. Using causal mediation analysis, we found statistically significant indirect effects of methylation (β(95% confidence interval): 6.9(1.5, 12.4)) providing evidence that DNA methylation mediates the relationship between PON1-108 genotype and PON1 expression. Our findings show that integration of genetic, epigenetic, and expression data can shed light on the functional mechanisms involving genetic and epigenetic regulation of candidate susceptibility genes like PON1.

Single-cell transcriptomics allows researchers to investigate complex communities of heterogeneous cells. It can be applied to stem cells and their descendants in order to chart the progression from multipotent progenitors to fully differentiated cells. While a variety of statistical and computational methods have been proposed for inferring cell lineages, the problem of accurately characterizing multiple branching lineages remains difficult to solve.

The oncogenic MUC1-C protein promotes dedifferentiation of castrate-resistant prostate cancer (CRPC) and triple-negative breast cancer (TNBC) cells. Chromatin remodeling is critical for the cancer stem cell (CSC) state; however, there is no definitive evidence that MUC1-C regulates chromatin accessibility and thereby expression of stemness-associated genes. We demonstrate that MUC1-C drives global changes in chromatin architecture in the dedifferentiation of CRPC and TNBC cells. Our results show that MUC1-C induces differentially accessible regions (DAR) across their genomes, which are significantly associated with differentially expressed genes (DEG). Motif and cistrome analysis further demonstrated MUC1-C-induced DARs align with genes regulated by the JUN/AP-1 family of transcription factors. MUC1-C activates the BAF chromatin remodeling complex, which is recruited by JUN in enhancer selection. In studies of the NOTCH1 gene, which is required for CRPC and TNBC cell self-renewal, we demonstrate that MUC1-C is necessary for (i) occupancy of JUN and ARID1A/BAF, (ii) increases in H3K27ac and H3K4me3 signals, and (iii) opening of chromatin accessibility on a proximal enhancer-like signature. Studies of the EGR1 and LY6E stemness-associated genes further demonstrate that MUC1-C-induced JUN/ARID1A complexes regulate chromatin accessibility on proximal and distal enhancer-like signatures. These findings uncover a role for MUC1-C in chromatin remodeling that is mediated at least in part by JUN/AP-1 and ARID1A/BAF in association with driving the CSC state.

Trajectory inference has radically enhanced single-cell RNA-seq research by enabling the study of dynamic changes in gene expression. Downstream of trajectory inference, it is vital to discover genes that are (i) associated with the lineages in the trajectory, or (ii) differentially expressed between lineages, to illuminate the underlying biological processes. Current data analysis procedures, however, either fail to exploit the continuous resolution provided by trajectory inference, or fail to pinpoint the exact types of differential expression. We introduce tradeSeq, a powerful generalized additive model framework based on the negative binomial distribution that allows flexible inference of both within-lineage and between-lineage differential expression. By incorporating observation-level weights, the model additionally allows to account for zero inflation. We evaluate the method on simulated datasets and on real datasets from droplet-based and full-length protocols, and show that it yields biological insights through a clear interpretation of the data.

The olfactory epithelium is one of the few regions of the nervous system that sustains neurogenesis throughout life. Its experimental accessibility makes it especially tractable for studying molecular mechanisms that drive neural regeneration after injury-induced cell death. In this study, we used single cell sequencing to identify major regulatory players in determining olfactory epithelial stem cell fate after acute injury. We combined gene expression and accessible chromatin profiles of individual lineage traced olfactory stem cells to predict transcription factor activity specific to different lineages and stages of recovery. We further identified a discrete stem cell state that appears poised for activation, characterized by accessible chromatin around wound response and lineage specific genes prior to their later expression in response to injury. Together these results provide evidence that a subset of quiescent olfactory epithelial stem cells are epigenetically primed to support injury-induced regeneration.

Tissue homeostasis and regeneration are mediated by programs of adult stem cell renewal and differentiation. However, the mechanisms that regulate stem cell fates under such widely varying conditions are not fully understood. Using single-cell techniques, we assessed the transcriptional changes associated with stem cell self-renewal and differentiation and followed the maturation of stem cell-derived clones using sparse lineage tracing in the regenerating mouse olfactory epithelium. Following injury, quiescent olfactory stem cells rapidly shift to activated, transient states unique to regeneration and tailored to meet the demands of injury-induced repair, including barrier formation and proliferation. Multiple cell fates, including renewed stem cells and committed differentiating progenitors, are specified during this early window of activation. We further show that Sox2 is essential for cells to transition from the activated to neuronal progenitor states. Our study highlights strategies for stem cell-mediated regeneration that may be conserved in other adult stem cell niches.

Altered olfactory function is a common symptom of COVID-19, but its etiology is unknown. A key question is whether SARS-CoV-2 (CoV-2) - the causal agent in COVID-19 - affects olfaction directly, by infecting olfactory sensory neurons or their targets in the olfactory bulb, or indirectly, through perturbation of supporting cells. Here we identify cell types in the olfactory epithelium and olfactory bulb that express SARS-CoV-2 cell entry molecules. Bulk sequencing demonstrated that mouse, non-human primate and human olfactory mucosa expresses two key genes involved in CoV-2 entry, ACE2 and TMPRSS2. However, single cell sequencing revealed that ACE2 is expressed in support cells, stem cells, and perivascular cells, rather than in neurons. Immunostaining confirmed these results and revealed pervasive expression of ACE2 protein in dorsally-located olfactory epithelial sustentacular cells and olfactory bulb pericytes in the mouse. These findings suggest that CoV-2 infection of non-neuronal cell types leads to anosmia and related disturbances in odor perception in COVID-19 patients.

Single-cell transcriptomics can provide quantitative molecular signatures for large, unbiased samples of the diverse cell types in the brain1-3. With the proliferation of multi-omics datasets, a major challenge is to validate and integrate results into a biological understanding of cell-type organization. Here we generated transcriptomes and epigenomes from more than 500,000 individual cells in the mouse primary motor cortex, a structure that has an evolutionarily conserved role in locomotion. We developed computational and statistical methods to integrate multimodal data and quantitatively validate cell-type reproducibility. The resulting reference atlas-containing over 56 neuronal cell types that are highly replicable across analysis methods, sequencing technologies and modalities-is a comprehensive molecular and genomic account of the diverse neuronal and non-neuronal cell types in the mouse primary motor cortex. The atlas includes a population of excitatory neurons that resemble pyramidal cells in layer 4 in other cortical regions4. We further discovered thousands of concordant marker genes and gene regulatory elements for these cell types. Our results highlight the complex molecular regulation of cell types in the brain and will directly enable the design of reagents to target specific cell types in the mouse primary motor cortex for functional analysis.

In single-cell RNA sequencing (scRNA-Seq), gene expression is assessed individually for each cell, allowing the investigation of developmental processes, such as embryogenesis and cellular differentiation and regeneration, at unprecedented resolution. In such dynamic biological systems, cellular states form a continuum, e.g., for the differentiation of stem cells into mature cell types. This process is often represented via a trajectory in a reduced-dimensional representation of the scRNA-Seq dataset. While many methods have been suggested for trajectory inference, it is often unclear how to handle multiple biological groups or conditions, e.g., inferring and comparing the differentiation trajectories of wild-type and knock-out stem cell populations. In this manuscript, we present condiments, a method for the inference and downstream interpretation of cell trajectories across multiple conditions. Our framework allows the interpretation of differences between conditions at the trajectory, cell population, and gene expression levels. We start by integrating datasets from multiple conditions into a single trajectory. By comparing the cell's conditions along the trajectory's path, we can detect large-scale changes, indicative of differential progression or fate selection. We also demonstrate how to detect subtler changes by finding genes that exhibit different behaviors between these conditions along a differentiation path.

The childhood salivary microbiome, which plays an important role in healthy development, may be influenced by breast milk consumption. The composition of the milk microbiome and the role it plays in the establishment of the infant microbiome are not well understood.

The tumor immune microenvironment plays a critical role in cancer progression and response to immunotherapy in clear cell renal cell carcinoma (ccRCC), yet the composition and phenotypic states of immune cells in this tumor are incompletely characterized. We performed single-cell RNA and T cell receptor sequencing on 164,722 individual cells from tumor and adjacent non-tumor tissue in patients with ccRCC across disease stages: early, locally advanced, and advanced/metastatic. Terminally exhausted CD8+ T cells were enriched in metastatic disease and were restricted in T cell receptor diversity. Within the myeloid compartment, pro-inflammatory macrophages were decreased, and suppressive M2-like macrophages were increased in advanced disease. Terminally exhausted CD8+ T cells and M2-like macrophages co-occurred in advanced disease and expressed ligands and receptors that support T cell dysfunction and M2-like polarization. This immune dysfunction circuit is associated with a worse prognosis in external cohorts and identifies potentially targetable immune inhibitory pathways in ccRCC.