Cell subtype-specific effects of genetic variation in the Alzheimer's disease brain.

The relationship between genetic variation and gene expression in brain cell types and subtypes remains understudied. Here, we generated single-nucleus RNA sequencing data from the neocortex of 424 individuals of advanced age; we assessed the effect of genetic variants on RNA expression in cis (cis-expression quantitative trait loci) for seven cell types and 64 cell subtypes using 1.5 million transcriptomes. This effort identified 10,004 eGenes at the cell type level and 8,099 eGenes at the cell subtype level. Many eGenes are only detected within cell subtypes. A new variant influences APOE expression only in microglia and is associated with greater cerebral amyloid angiopathy but not Alzheimer's disease pathology, after adjusting for APOEε4, providing mechanistic insights into both pathologies. Furthermore, only a TMEM106B variant affects the proportion of cell subtypes. Integration of these results with genome-wide association studies highlighted the targeted cell type and probable causal gene within Alzheimer's disease, schizophrenia, educational attainment and Parkinson's disease loci.

Systematically characterizing the roles of E3-ligase family members in inflammatory responses with massively parallel Perturb-seq.

E3 ligases regulate key processes, but many of their roles remain unknown. Using Perturb-seq, we interrogated the function of 1,130 E3 ligases, partners and substrates in the inflammatory response in primary dendritic cells (DCs). Dozens impacted the balance of DC1, DC2, migratory DC and macrophage states and a gradient of DC maturation. Family members grouped into co-functional modules that were enriched for physical interactions and impacted specific programs through substrate transcription factors. E3s and their adaptors co-regulated the same processes, but partnered with different substrate recognition adaptors to impact distinct aspects of the DC life cycle. Genetic interactions were more prevalent within than between modules, and a deep learning model, comßVAE, predicts the outcome of new combinations by leveraging modularity. The E3 regulatory network was associated with heritable variation and aberrant gene expression in immune cells in human inflammatory diseases. Our study provides a general approach to dissect gene function.

Functional analyses and single cell immunoprofiling uncover sex-specific differences in SARS-CoV2 immune memory development.

SARS-CoV-2 infection leads to a broad range of outcomes and immune responses, with the development of neutralizing antibodies generally correlated with protection against reinfection. Here, we have characterized both neutralizing activity and T cell responses in a cluster of subjects with mild disease linked to a single spreading event. Surprisingly, we observed sex-specific associations between spike- and particularly nucleoprotein-specific T cell responses and neutralization, with pro-inflammatory cytokines being linked to higher titers only in males. Using single cell immunoprofiling, which provided matched transcriptome and T-cell receptor (TCR) profiles in restimulated CD4 + and CD8 + cells from these subjects, we identified differences in type I IFN signaling that may underlie this difference in antibody generation. Finally, we also identified several TCRs associated with cytokine producing T cells. Altogether, our work maps the breadth of immunological outcomes of SARS-CoV2 infections and highlight the potential role of sex-specific feedback loops during the generation of neutralizing antibodies.

Cycling cancer persister cells arise from lineages with distinct programs.

Non-genetic mechanisms have recently emerged as important drivers of cancer therapy failure1, where some cancer cells can enter a reversible drug-tolerant persister state in response to treatment2. Although most cancer persisters remain arrested in the presence of the drug, a rare subset can re-enter the cell cycle under constitutive drug treatment. Little is known about the non-genetic mechanisms that enable cancer persisters to maintain proliferative capacity in the presence of drugs. To study this rare, transiently resistant, proliferative persister population, we developed Watermelon, a high-complexity expressed barcode lentiviral library for simultaneous tracing of each cell's clonal origin and proliferative and transcriptional states. Here we show that cycling and non-cycling persisters arise from different cell lineages with distinct transcriptional and metabolic programs. Upregulation of antioxidant gene programs and a metabolic shift to fatty acid oxidation are associated with persister proliferative capacity across multiple cancer types. Impeding oxidative stress or metabolic reprogramming alters the fraction of cycling persisters. In human tumours, programs associated with cycling persisters are induced in minimal residual disease in response to multiple targeted therapies. The Watermelon system enabled the identification of rare persister lineages that are preferentially poised to proliferate under drug pressure, thus exposing new vulnerabilities that can be targeted to delay or even prevent disease recurrence.

COVID-19 tissue atlases reveal SARS-CoV-2 pathology and cellular targets.

COVID-19, which is caused by SARS-CoV-2, can result in acute respiratory distress syndrome and multiple organ failure1-4, but little is known about its pathophysiology. Here we generated single-cell atlases of 24 lung, 16 kidney, 16 liver and 19 heart autopsy tissue samples and spatial atlases of 14 lung samples from donors who died of COVID-19. Integrated computational analysis uncovered substantial remodelling in the lung epithelial, immune and stromal compartments, with evidence of multiple paths of failed tissue regeneration, including defective alveolar type 2 differentiation and expansion of fibroblasts and putative TP63+ intrapulmonary basal-like progenitor cells. Viral RNAs were enriched in mononuclear phagocytic and endothelial lung cells, which induced specific host programs. Spatial analysis in lung distinguished inflammatory host responses in lung regions with and without viral RNA. Analysis of the other tissue atlases showed transcriptional alterations in multiple cell types in heart tissue from donors with COVID-19, and mapped cell types and genes implicated with disease severity based on COVID-19 genome-wide association studies. Our foundational dataset elucidates the biological effect of severe SARS-CoV-2 infection across the body, a key step towards new treatments.

A single-cell and spatial atlas of autopsy tissues reveals pathology and cellular targets of SARS-CoV-2.

The SARS-CoV-2 pandemic has caused over 1 million deaths globally, mostly due to acute lung injury and acute respiratory distress syndrome, or direct complications resulting in multiple-organ failures. Little is known about the host tissue immune and cellular responses associated with COVID-19 infection, symptoms, and lethality. To address this, we collected tissues from 11 organs during the clinical autopsy of 17 individuals who succumbed to COVID-19, resulting in a tissue bank of approximately 420 specimens. We generated comprehensive cellular maps capturing COVID-19 biology related to patients' demise through single-cell and single-nucleus RNA-Seq of lung, kidney, liver and heart tissues, and further contextualized our findings through spatial RNA profiling of distinct lung regions. We developed a computational framework that incorporates removal of ambient RNA and automated cell type annotation to facilitate comparison with other healthy and diseased tissue atlases. In the lung, we uncovered significantly altered transcriptional programs within the epithelial, immune, and stromal compartments and cell intrinsic changes in multiple cell types relative to lung tissue from healthy controls. We observed evidence of: alveolar type 2 (AT2) differentiation replacing depleted alveolar type 1 (AT1) lung epithelial cells, as previously seen in fibrosis; a concomitant increase in myofibroblasts reflective of defective tissue repair; and, putative TP63+ intrapulmonary basal-like progenitor (IPBLP) cells, similar to cells identified in H1N1 influenza, that may serve as an emergency cellular reserve for severely damaged alveoli. Together, these findings suggest the activation and failure of multiple avenues for regeneration of the epithelium in these terminal lungs. SARS-CoV-2 RNA reads were enriched in lung mononuclear phagocytic cells and endothelial cells, and these cells expressed distinct host response transcriptional programs. We corroborated the compositional and transcriptional changes in lung tissue through spatial analysis of RNA profiles in situ and distinguished unique tissue host responses between regions with and without viral RNA, and in COVID-19 donor tissues relative to healthy lung. Finally, we analyzed genetic regions implicated in COVID-19 GWAS with transcriptomic data to implicate specific cell types and genes associated with disease severity. Overall, our COVID-19 cell atlas is a foundational dataset to better understand the biological impact of SARS-CoV-2 infection across the human body and empowers the identification of new therapeutic interventions and prevention strategies.

Mechanical forces drive ordered patterning of hair cells in the mammalian inner ear.

Periodic organization of cells is required for the function of many organs and tissues. The development of such periodic patterns is typically associated with mechanisms based on intercellular signaling such as lateral inhibition and Turing patterning. Here we show that the transition from disordered to ordered checkerboard-like pattern of hair cells and supporting cells in the mammalian hearing organ, the organ of Corti, is likely based on mechanical forces rather than signaling events. Using time-lapse imaging of mouse cochlear explants, we show that hair cells rearrange gradually into a checkerboard-like pattern through a tissue-wide shear motion that coordinates intercalation and delamination events. Using mechanical models of the tissue, we show that global shear and local repulsion forces on hair cells are sufficient to drive the transition from disordered to ordered cellular pattern. Our findings suggest that mechanical forces drive ordered hair cell patterning in a process strikingly analogous to the process of shear-induced crystallization in polymer and granular physics.

A synthetic planar cell polarity system reveals localized feedback on Fat4-Ds1 complexes.

The atypical cadherins Fat and Dachsous (Ds) have been found to underlie planar cell polarity (PCP) in many tissues. Theoretical models suggest that polarity can arise from localized feedbacks on Fat-Ds complexes at the cell boundary. However, there is currently no direct evidence for the existence or mechanism of such feedbacks. To directly test the localized feedback model, we developed a synthetic biology platform based on mammalian cells expressing the human Fat4 and Ds1. We show that Fat4-Ds1 complexes accumulate on cell boundaries in a threshold-like manner and exhibit dramatically slower dynamics than unbound Fat4 and Ds1. This suggests a localized feedback mechanism based on enhanced stability of Fat4-Ds1 complexes. We also show that co-expression of Fat4 and Ds1 in the same cells is sufficient to induce polarization of Fat4-Ds1 complexes. Together, these results provide direct evidence that localized feedbacks on Fat4-Ds1 complexes can give rise to PCP.

Cell-Cell Contact Area Affects Notch Signaling and Notch-Dependent Patterning.

During development, cells undergo dramatic changes in their morphology. By affecting contact geometry, these morphological changes could influence cellular communication. However, it has remained unclear whether and how signaling depends on contact geometry. This question is particularly relevant for Notch signaling, which coordinates neighboring cell fates through direct cell-cell signaling. Using micropatterning with a receptor trans-endocytosis assay, we show that signaling between pairs of cells correlates with their contact area. This relationship extends across contact diameters ranging from micrometers to tens of micrometers. Mathematical modeling predicts that dependence of signaling on contact area can bias cellular differentiation in Notch-mediated lateral inhibition processes, such that smaller cells are more likely to differentiate into signal-producing cells. Consistent with this prediction, analysis of developing chick inner ear revealed that ligand-producing hair cell precursors have smaller apical footprints than non-hair cells. Together, these results highlight the influence of cell morphology on fate determination processes.

Quantitative Analysis of Delta-like 1 Membrane Dynamics Elucidates the Role of Contact Geometry on Notch Signaling.

Notch signaling is ubiquitously used to coordinate differentiation between adjacent cells across metazoans. Whereas Notch pathway components have been studied extensively, the effect of membrane distribution and dynamics of Notch receptors and ligands remains poorly understood. It is also unclear how cellular morphology affects these distributions and, ultimately, the signaling between cells. Here, we combine live-cell imaging and mathematical modeling to address these questions. We use a FRAP-TIRF assay to measure the diffusion and endocytosis rates of Delta-like 1 (Dll1) in mammalian cells. We find large cell-to-cell variability in the diffusion coefficients of Dll1 measured in single cells within the same population. Using a simple reaction-diffusion model, we show how membrane dynamics and cell morphology affect cell-cell signaling. We find that differences in the diffusion coefficients, as observed experimentally, can dramatically affect signaling between cells. Together, these results elucidate how membrane dynamics and cellular geometry can affect cell-cell signaling.

The GPSM2/LGN GoLoco motifs are essential for hearing.

The planar cell polarity (PCP) pathway is responsible for polarizing and orienting cochlear hair cells during development through movement of a primary cilium, the kinocilium. GPSM2/LGN, a mitotic spindle-orienting protein associated with deafness in humans, is a PCP effector involved in kinocilium migration. Here, we link human and mouse truncating mutations in the GPSM2/LGN gene, both leading to hearing loss. The human variant, p.(Trp326*), was identified by targeted genomic enrichment of genes associated with deafness, followed by massively parallel sequencing. Lgn (ΔC) mice, with a targeted deletion truncating the C-terminal GoLoco motifs, are profoundly deaf and show misorientation of the hair bundle and severe malformations in stereocilia shape that deteriorates over time. Full-length protein levels are greatly reduced in mutant mice, with upregulated mRNA levels. The truncated Lgn (ΔC) allele is translated in vitro, suggesting that mutant mice may have partially functioning Lgn. Gαi and aPKC, known to function in the same pathway as Lgn, are dependent on Lgn for proper localization. The polarization of core PCP proteins is not affected in Lgn mutants; however, Lgn and Gαi are misoriented in a PCP mutant, supporting the role of Lgn as a PCP effector. The kinocilium, previously shown to be dependent on Lgn for robust localization, is essential for proper localization of Lgn, as well as Gαi and aPKC, suggesting that cilium function plays a role in positioning of apical proteins. Taken together, our data provide a mechanism for the loss of hearing found in human patients with GPSM2/LGN variants.

DSIF restricts NF-κB signaling by coordinating elongation with mRNA processing of negative feedback genes.

NF-κB is central for immune response and cell survival, and its deregulation is linked to chronic inflammation and cancer through poorly defined mechanisms. IκBα and A20 are NF-κB target genes and negative feedback regulators. Upon their activation by NF-κB, DSIF is recruited, P-TEFb is released, and their elongating polymerase II (Pol II) C-terminal domain (CTD) remains hypophosphorylated. We show that upon DSIF knockdown, mRNA levels of a subset of NF-κB targets are not diminished; yet much less IκBα and A20 protein are synthesized, and NF-κB activation is abnormally prolonged. Further analysis of IκBα and A20 mRNA revealed that a significant portion is uncapped, unspliced, and retained in the nucleus. Interestingly, the Spt5 C-terminal repeat (CTR) domain involved in elongation stimulation through P-TEFb is dispensable for IκBα and A20 regulation. These findings assign a function for DSIF in cotranscriptional mRNA processing when elongating Pol II is hypophosphorylated and define DSIF as part of the negative feedback regulation of NF-κB.

Homeodomain protein otp and activity-dependent splicing modulate neuronal adaptation to stress.

Regulation of corticotropin-releasing hormone (CRH) activity is critical for the animal's adaptation to stressful challenges, and its dysregulation is associated with psychiatric disorders in humans. However, the molecular mechanism underlying this transcriptional response to stress is not well understood. Using various stress paradigms in mouse and zebrafish, we show that the hypothalamic transcription factor Orthopedia modulates the expression of CRH as well as the splicing factor Ataxin 2-Binding Protein-1 (A2BP1/Rbfox-1). We further show that the G protein coupled receptor PAC1, which is a known A2BP1/Rbfox-1 splicing target and an important mediator of CRH activity, is alternatively spliced in response to a stressful challenge. The generation of PAC1-hop messenger RNA isoform by alternative splicing is required for termination of CRH transcription, normal activation of the hypothalamic-pituitary-adrenal axis and adaptive anxiety-like behavior. Our study identifies an evolutionarily conserved biochemical pathway that modulates the neuronal adaptation to stress through transcriptional activation and alternative splicing.

The metabolic regulator PGC-1α directly controls the expression of the hypothalamic neuropeptide oxytocin.

The transcriptional coactivator PGC-1α is a key regulator of cellular energy expenditure in peripheral tissues. Recent studies report that PGC-1α-null mice develop late-onset obesity and that the neuronal inactivation of PGC-1α causes increased food intake. However, the exact role of PGC-1α in the CNS remains unclear. Here we show that PGC-1α directly regulates the expression of the hypothalamic neuropeptide oxytocin, a known central regulator of appetite. We developed a unique genetic approach in the zebrafish, allowing us to monitor and manipulate PGC-1α activity in oxytocinergic neurons. We found that PGC-1α is coexpressed with oxytocin in the zebrafish hypothalamus. Targeted knockdown of the zebrafish PGC-1α gene activity caused a marked decrease in oxytocin mRNA levels and inhibited the expression of a transgenic GFP reporter driven by the oxytocin promoter. The effect of PGC-1α loss of function on oxytocin gene activity was rescued by tissue-specific re-expression of either PGC-1α or oxytocin precursor in zebrafish oxytocinergic neurons. PGC-1α activated the oxytocin promoter in a heterologous cell culture system, and overexpression of PGC-1α induced ectopic expression of oxytocin in muscles and neurons. Finally, PGC-1α forms an in vivo complex with the oxytocin promoter in fed but not fasted animals. These findings demonstrate that PGC-1α is both necessary and sufficient for the production of oxytocin, implicating hypothalamic PGC-1α in the direct activation of a hypothalamic hormone known to control energy intake.

Interplay between E-box and NF-kappaB in regulation of A20 gene by DRB sensitivity-inducing factor (DSIF).

The NF-kappaB target gene A20 serves as a paradigm for gene-specific control of transcription elongation. This gene is regulated by the elongation factor DSIF (DRB sensitivity-inducing factor) under basal and NF-kappaB-activated states by two distinct mechanisms. Prior to NF-kappaB stimulation, the A20 gene is occupied by polymerase II, and elongation is inhibited by DSIF. This inhibition is mediated by an upstream promoter element termed ELIE (elongation inhibitory element). Upon NF-kappaB activation, inhibition of the A20 gene by DSIF persists, but now NF-kappaB and the core promoter regulate DSIF instead of ELIE. Here we investigated the regulation of DSIF by ELIE and the regulatory switch from ELIE to NF-kappaB following NF-kappaB induction. Electrophoretic mobility shift assays revealed two distinct protein complexes that specifically interact with ELIE, one of which is the E-box protein USF1. Interestingly, USF1 is displaced from the A20 promoter upon induction of NF-kappaB. A mutation in the E-box section of ELIE diminished the binding of USF1 and DSIF recruitment. Consistent with these findings, the E-box is crucial for DSIF inhibition in resting, but not NF-kappaB-stimulated, cells. These findings reveal a dynamic regulation of DSIF involving either E-box or NF-kappaB depending on the physiological circumstances.

Differential regulation of NF-kappaB by elongation factors is determined by core promoter type.

NF-kappaB transcription factors activate genes important for immune response, inflammation, and cell survival. P-TEFb and DSIF, which are positive and negative transcription elongation factors, respectively, both regulate NF-kappaB-induced transcription, but the mechanism underlying their recruitment to NF-kappaB target genes is unknown. We show here that upon induction of NF-kappaB, a subset of target genes is regulated differentially by either P-TEFb or DSIF. The regulation of these genes and their occupancy by these elongation factors are dependent on the NF-kappaB enhancer and the core promoter type. Converting a TATA-less promoter to a TATA promoter switches the regulation of NF-kappaB from DSIF to P-TEFb. Accumulation or displacement of DSIF and P-TEFb is dictated by the formation of distinct initiation complexes (TFIID dependent or independent) on the two types of core promoter. The underlying mechanism for the dissociation of DSIF from TATA promoters upon NF-kappaB activation involves the phosphorylation of RNA polymerase II by P-TEFb. The results highlight a regulatory link between the initiation and the elongation phases of the transcription reaction and broaden our comprehension of the NF-kappaB pathway.

Elongation inhibition by DRB sensitivity-inducing factor is regulated by the A20 promoter via a novel negative element and NF-kappaB.

A20 is an immediate-early NF-kappaB target gene. Prior to NF-kappaB stimulation, the A20 promoter is bound by the polymerase II machinery to allow rapid transcription activation. Here we show that the basal A20 transcription is repressed at the level of elongation in a promoter-specific fashion. Immunodepletion in vitro and RNA interference in cultured cells suggest that the basal elongation inhibition is conferred by DRB sensitivity-inducing factor (DSIF). We have identified a negative upstream promoter element called ELIE that controls DSIF activity. Remarkably, following NF-kappaB stimulation, inhibition of the A20 promoter by DSIF persists, but it is now regulated by NF-kappaB rather than ELIE. Similar regulation by DSIF is shown for another NF-kappaB-responsive gene, the IkappaBalpha gene. These findings reveal an intimate and dynamic relationship between DSIF inhibition of elongation and promoter-bound transcription factors. The potential significance of the differential regulation of DSIF activity by cis-acting elements is discussed.