Cell-type-specific transcriptome architecture underlying the establishment and exacerbation of systemic lupus erythematosus

Systemic lupus erythematosus (SLE) is a complex and heterogeneous autoimmune disease involving multiple immune cells. A major hurdle to the elucidation of SLE pathogenesis is our limited understanding of dysregulated gene expression linked to various clinical statuses with a high cellular resolution. Here, we conducted a large-scale transcriptome study with 6,386 RNA sequencing data covering 27 immune cell types from 159 SLE and 89 healthy donors. We first profiled two distinct cell-type-specific transcriptomic signatures: disease-state and disease-activity signatures, reflecting disease establishment and exacerbation, respectively. We next identified candidate biological processes unique to each signature. This study suggested the clinical value of disease-activity signatures, which were associated with organ involvement and responses to therapeutic agents such as belimumab. However, disease-activity signatures were less enriched around SLE risk variants than disease-state signatures, suggesting that the genetic studies to date may not well capture clinically vital biology in SLE. Together, we identified comprehensive gene signatures of SLE, which will provide essential foundations for future genomic, genetic, and clinical studies.

Gardnerella Vaginalis Alters Cervicovaginal Epithelial Cell Function Through Epithelial Cell-type Specific Immune Responses

Background: The cervicovaginal (CV) microbiome is highly associated with vaginal health and disease in both pregnant and non-pregnant individuals. An overabundance of Gardnerella vaginalis in the CV space is commonly associated with adverse reproductive outcomes including bacterial vaginosis (BV), sexually transmitted diseases and preterm birth while the presence of Lactobacillus spp is often associated with reproductive health. While host-microbial interactions are hypothesized to contribute to CV health and disease, the mechanisms by which these interactions regulate CV epithelial function remain largely unknown. Results: Using an in vitro co-culture model, we assessed the effects of Lactobacillus crispatus and G. vaginalis on the CV epithelial barrier, the immune mediators that could be contributing to decreased barrier integrity and the immune signaling pathways regulating the immune response. G. vaginalis, but not L. crispatus, significantly increased epithelial cell death and decreased epithelial barrier integrity in an epithelial cell-specific manner. A G. vaginalis-mediated epithelial immune response including NFkB activation and proinflammatory cytokine release was initiated partially through TLR2 dependent signaling pathways. Additionally, investigation of the cytokine immune profile in human CV fluid showed distinctive clustering of cytokines by G. vaginalis abundance and birth outcome. Conclusions: The results of this study show both microbe- and epithelial cell-type specific effects on CV epithelial function. Altered epithelial barrier function through cell death and immune mediated mechanisms by G. vaginalis, but not L. crispatus, indicates that host epithelial cells respond to bacteria-associated signals, resulting in altered epithelial function and ultimately CV disease. Additionally, distinct immune signatures associated with G. vaginalis or birth outcome provide further evidence that host-microbial interactions may contribute significantly to the biological mechanisms regulating reproductive outcomes.

In vivo single-cell profiling of lncRNAs during Ebola virus infection

Long non-coding RNAs (lncRNAs) are pivotal mediators of systemic immune response to viral infection, yet most studies concerning their expression and functions upon immune stimulation are limited to in vitro bulk cell populations. This strongly constrains our understanding of how lncRNA expression varies at single-cell resolution, and how their cell-type specific immune regulatory roles may differ compared to protein-coding genes. Here, we perform the first in-depth characterization of lncRNA expression variation at single-cell resolution during Ebola virus (EBOV) infection in vivo. Using bulk RNA-sequencing from 119 samples and 12 tissue types, we significantly expand the current macaque lncRNA annotation. We then profile lncRNA expression variation in immune circulating single-cells during EBOV infection and find that lncRNAs' expression in fewer cells is a major differentiating factor from their protein-coding gene counterparts. Upon EBOV infection, lncRNAs present dynamic and mostly cell-type specific changes in their expression profiles especially in monocytes, the main cell type targeted by EBOV. Such changes are associated with gene regulatory modules related to important innate immune responses such as interferon response and purine metabolism. Within infected cells, several lncRNAs have positively and negatively correlated expression with viral load, suggesting that expression of some of these lncRNAs might be directly hijacked by EBOV to attack host cells. This study provides novel insights into the roles that lncRNAs play in the host response to acute viral infection and paves the way for future lncRNA studies at single-cell resolution.

Cell-Type-Specific Circadian Bioluminescence Rhythms in Dbp Reporter Mice

Circadian rhythms are endogenously generated physiological and molecular rhythms with a cycle length of about 24 h. Bioluminescent reporters have been exceptionally useful for studying circadian rhythms in numerous species. Here, we report development of a reporter mouse generated by modification of a widely expressed and highly rhythmic gene encoding D-site albumin promoter binding protein ( Dbp). In this line of mice, firefly luciferase is expressed from the Dbp locus in a Cre recombinase-dependent manner, allowing assessment of bioluminescence rhythms in specific cellular populations. A mouse line in which luciferase expression was Cre-independent was also generated. The Dbp reporter alleles do not alter Dbp gene expression rhythms in liver or circadian locomotor activity rhythms. In vivo and ex vivo studies show the utility of the reporter alleles for monitoring rhythmicity. Our studies reveal cell-type-specific characteristics of rhythms among neuronal populations within the suprachiasmatic nuclei ex vivo. In vivo studies show Dbp-driven bioluminescence rhythms in the liver of Albumin-Cre;Dbp KI/+ “liver reporter” mice. After a shift of the lighting schedule, locomotor activity achieved the proper phase relationship with the new lighting cycle more rapidly than hepatic bioluminescence did. As previously shown, restricting food access to the daytime altered the phase of hepatic rhythmicity. Our model allowed assessment of the rate of recovery from misalignment once animals were provided with food ad libitum. These studies confirm the previously demonstrated circadian misalignment following environmental perturbations and reveal the utility of this model for minimally invasive, longitudinal monitoring of rhythmicity from specific mouse tissues.

Cell type-specific aging clocks to quantify aging and rejuvenation in regenerative regions of the brain

Aging manifests as progressive dysfunction culminating in death. The diversity of cell types is a challenge to the precise quantification of aging and its reversal. Here we develop a suite of 'aging clocks' based on single cell transcriptomic data to characterize cell type-specific aging and rejuvenation strategies. The subventricular zone (SVZ) neurogenic region contains many cell types and provides an excellent system to study cell-level tissue aging and regeneration. We generated 21,458 single-cell transcriptomes from the neurogenic regions of 28 mice, tiling ages from young to old. With these data, we trained a suite of single cell-based regression models (aging clocks) to predict both chronological age (passage of time) and biological age (fitness, in this case the proliferative capacity of the neurogenic region). Both types of clocks perform well on independent cohorts of mice. Genes underlying the single cell-based aging clocks are mostly cell-type specific, but also include a few shared genes in the interferon and lipid metabolism pathways. We used these single cell-based aging clocks to measure transcriptomic rejuvenation, by generating single cell RNA-seq datasets of SVZ neurogenic regions for two interventions - heterochronic parabiosis (young blood) and exercise. Interestingly, the use of aging clocks reveals that both heterochronic parabiosis and exercise reverse transcriptomic aging in the niche, but in different ways across cell types and genes. This study represents the first development of high-resolution aging clocks from single cell transcriptomic data and demonstrates their application to quantify transcriptomic rejuvenation.

Structural insights in cell-type specific evolution of intra-host diversity by SARS-CoV-2

As the global burden of SARS-CoV-2 infections escalates, so does the evolution of viral variants with increased transmissibility and pathology. In addition to this entrenched diversity, RNA viruses can also display genetic diversity within single infected hosts with co-existing viral variants evolving differently in distinct cell types. The BriSΔ variant, originally identified as a viral subpopulation from SARS-CoV-2 isolate hCoV-19/England/02/2020, comprises in the spike an eight amino-acid deletion encompassing a furin recognition motif and S1/S2 cleavage site. We elucidate the structure, function and molecular dynamics of this spike providing mechanistic insight into how the deletion correlates to viral cell tropism, ACE2 receptor binding and infectivity of this SARS-CoV-2 variant. Our results reveal long-range allosteric communication between functional domains that differ in the wild-type and the deletion variant and support a view of SARS-CoV-2 probing multiple evolutionary trajectories in distinct cell types within the same infected host.

Single-cell network biology characterizes cell-type gene regulation for drug repurposing and phenotype prediction in Alzheimer's disease

Dysregulation of gene expression in Alzheimer's disease (AD) remains elusive, especially at the cell type level. Gene regulatory network, a key molecular mechanism linking transcription factors (TFs) and regulatory elements to govern target gene expression, can change across cell types in the human brain and thus serve as a model for studying gene dysregulation in AD. However, it is still challenging to understand how cell type networks work abnormally under AD. To address this, we integrated single-cell multi-omics data and predicted the gene regulatory networks in AD and control for four major cell types, excitatory and inhibitory neurons, microglia and oligodendrocytes. Importantly, we applied network biology approaches to analyze the changes of network characteristics across these cell types, and between AD and control. For instance, many hub TFs target different genes between AD and control (rewiring). Also, these networks show strong hierarchical structures in which top TFs (master regulators) are largely common across cell types, whereas different TFs operate at the middle levels in some cell types (e.g., microglia). The regulatory logics of enriched network motifs (e.g., feed-forward loops) further uncover cell-type-specific TF-TF cooperativities in gene regulation. The cell type networks are highly modular. Several network modules with cell-type-specific expression changes in AD pathology are enriched with AD-risk genes and putative targets of approved and pending AD drugs, suggesting possible cell-type genomic medicine in AD. Finally, using the cell type gene regulatory networks, we developed machine learning models to classify and prioritize additional AD genes. We found that top prioritized genes predict clinical phenotypes (e.g., cognitive impairment). Overall, this single-cell network biology analysis provides a comprehensive map linking genes, regulatory networks, cell types and drug targets and reveals mechanisms on cell-type gene dyregulation in AD.

The Cell Differentiation of Idioblast Myrosin Cells: Similarities With Vascular and Guard Cells

Idioblasts are defined by abnormal shapes, sizes, and contents that are different from neighboring cells. Myrosin cells are Brassicales-specific idioblasts and accumulate a large amount of thioglucoside glucohydrolases (TGGs, also known as myrosinases) in their vacuoles. Myrosinases convert their substrates, glucosinolates, into toxic compounds when herbivories and pests attack plants. In this review, we highlight the similarities and differences between myrosin cells and vascular cells/guard cells (GCs) because myrosin cells are distributed along vascular cells, especially the phloem parenchyma, and myrosin cells share the master transcription factor FAMA with GCs for their cell differentiation. In addition, we analyzed the overlap of cell type-specific genes between myrosin cells and GCs by using published single-cell transcriptomics (scRNA-seq) data, suggesting significant similarities in the gene expression patterns of these two specialized cells.

Neuronal Yin Yang1 in the prefrontal cortex regulates transcriptional and behavioral responses to chronic stress in mice

Although the synaptic alterations associated with the stress-related mood disorder major depression has been well-documented, the underlying transcriptional mechanisms remain poorly understood. Here, we perform complementary bulk nuclei- and single-nucleus transcriptome profiling and map locus-specific chromatin interactions in mouse neocortex to identify the cell type-specific transcriptional changes associated with stress-induced behavioral maladaptation. We find that cortical excitatory neurons, layer 2/3 neurons in particular, are vulnerable to chronic stress and acquire signatures of gene transcription and chromatin structure associated with reduced neuronal activity and expression of Yin Yang 1 (YY1). Selective ablation of YY1 in cortical excitatory neurons enhances stress sensitivity in both male and female mice and alters the expression of stress-associated genes following an abbreviated stress exposure. These findings demonstrate how chronic stress impacts transcription in cortical excitatory neurons and identify YY1 as a regulator of stress-induced maladaptive behavior in mice.