Dissecting mechanisms of chamber-specific cardiac differentiation and its perturbation following retinoic acid exposure.

The specification of distinct cardiac lineages occurs before chamber formation and acquisition of bona fide atrial or ventricular identity. However, the mechanisms underlying these early specification events remain poorly understood. Here, we performed single cell analysis at the murine cardiac crescent, primitive heart tube and heart tube stages to uncover the transcriptional mechanisms underlying formation of atrial and ventricular cells. We find that progression towards differentiated cardiomyocytes occurs primarily based on heart field progenitor identity, and that progenitors contribute to ventricular or atrial identity through distinct differentiation mechanisms. We identify new candidate markers that define such differentiation processes and examine their expression dynamics using computational lineage trajectory methods. We further show that exposure to exogenous retinoic acid causes defects in ventricular chamber size, dysregulation in FGF signaling and a shunt in differentiation towards orthogonal lineages. Retinoic acid also causes defects in cell-cycle exit resulting in formation of hypomorphic ventricles. Collectively, our data identify, at a single cell level, distinct lineage trajectories during cardiac specification and differentiation, and the precise effects of manipulating cardiac progenitor patterning via retinoic acid signaling.

Molecular Profiling of Coronavirus Disease 2019 (COVID-19) Autopsies Uncovers Novel Disease Mechanisms.

Current understanding of coronavirus disease 2019 (COVID-19) pathophysiology is limited by disease heterogeneity, complexity, and a paucity of studies assessing patient tissues with advanced molecular tools. Rapid autopsy tissues were evaluated using multiscale, next-generation RNA-sequencing methods (bulk, single-nuclei, and spatial transcriptomics) to provide unprecedented molecular resolution of COVID-19-induced damage. Comparison of infected/uninfected tissues revealed four major regulatory pathways. Effectors within these pathways could constitute novel therapeutic targets, including the complement receptor C3AR1, calcitonin receptor-like receptor, or decorin. Single-nuclei RNA sequencing of olfactory bulb and prefrontal cortex highlighted remarkable diversity of coronavirus receptors. Angiotensin-converting enzyme 2 was rarely expressed, whereas basigin showed diffuse expression, and alanyl aminopeptidase, membrane, was associated with vascular/mesenchymal cell types. Comparison of lung and lymph node tissues from patients with different symptoms (one had died after a month-long hospitalization with multiorgan involvement, and the other had died after a few days of respiratory symptoms) with digital spatial profiling resulted in distinct molecular phenotypes. Evaluation of COVID-19 rapid autopsy tissues with advanced molecular techniques can identify pathways and effectors, map diverse receptors at the single-cell level, and help dissect differences driving diverging clinical courses among individual patients. Extension of this approach to larger data sets will substantially advance the understanding of the mechanisms behind COVID-19 pathophysiology.

Integration of CRISPR-engineering and hiPSC-based models of psychiatric genomics.

Neuropsychiatric disorders are highly heritable polygenic disorders arising from the complex interplay of highly penetrant rare variants and common variants of small effect. There is a large index of comorbidity and shared genetic risk between disorders, reflecting the pleiotropy of individual variants as well as predicted downstream pathway-level convergence. Importantly, the mechanism(s) through which psychiatric disease-associated variants interact to contribute to disease risk remains unknown. Human induced pluripotent stem cell (hiPSC)-based models are increasingly useful for the systematic study of the complex genetics associated with brain diseases, particularly when combined with CRISPR-mediated genomic engineering, which together facilitate isogenic comparisons of defined neuronal cell types. In this review, we discuss the latest CRISPR technologies and consider how they can be successfully applied to the functional characterization of the growing list genetic variants linked to psychiatric disease.

New considerations for hiPSC-based models of neuropsychiatric disorders.

The development of human-induced pluripotent stem cells (hiPSCs) has made possible patient-specific modeling across the spectrum of human disease. Here, we discuss recent advances in psychiatric genomics and post-mortem studies that provide critical insights concerning cell-type composition and sample size that should be considered when designing hiPSC-based studies of complex genetic disease. We review recent hiPSC-based models of SZ, in light of our new understanding of critical power limitations in the design of hiPSC-based studies of complex genetic disorders. Three possible solutions are a movement towards genetically stratified cohorts of rare variant patients, application of CRISPR technologies to engineer isogenic neural cells to study the impact of common variants, and integration of advanced genetics and hiPSC-based datasets in future studies. Overall, we emphasize that to advance the reproducibility and relevance of hiPSC-based studies, stem cell biologists must contemplate statistical and biological considerations that are already well accepted in the field of genetics. We conclude with a discussion of the hypothesis of biological convergence of disease-through molecular, cellular, circuit, and patient level phenotypes-and how this might emerge through hiPSC-based studies.

AID stabilizes stem-cell phenotype by removing epigenetic memory of pluripotency genes.

The activation-induced cytidine deaminase (AID; also known as AICDA) enzyme is required for somatic hypermutation and class switch recombination at the immunoglobulin locus. In germinal-centre B cells, AID is highly expressed, and has an inherent mutator activity that helps generate antibody diversity. However, AID may also regulate gene expression epigenetically by directly deaminating 5-methylcytosine in concert with base-excision repair to exchange cytosine. This pathway promotes gene demethylation, thereby removing epigenetic memory. For example, AID promotes active demethylation of the genome in primordial germ cells. However, different studies have suggested either a requirement or a lack of function for AID in promoting pluripotency in somatic nuclei after fusion with embryonic stem cells. Here we tested directly whether AID regulates epigenetic memory by comparing the relative ability of cells lacking AID to reprogram from a differentiated murine cell type to an induced pluripotent stem cell. We show that Aid-null cells are transiently hyper-responsive to the reprogramming process. Although they initiate expression of pluripotency genes, they fail to stabilize in the pluripotent state. The genome of Aid-null cells remains hypermethylated in reprogramming cells, and hypermethylated genes associated with pluripotency fail to be stably upregulated, including many MYC target genes. Recent studies identified a late step of reprogramming associated with methylation status, and implicated a secondary set of pluripotency network components. AID regulates this late step, removing epigenetic memory to stabilize the pluripotent state.

A loss-of-function and H2B-Venus transcriptional reporter allele for Gata6 in mice.

BACKGROUND: The GATA-binding factor 6 (Gata6) gene encodes a zinc finger transcription factor that often functions as a key regulator of lineage specification during development. It is the earliest known marker of the primitive endoderm lineage in the mammalian blastocyst. During gastrulation, GATA6 is expressed in early cardiac mesoderm and definitive endoderm progenitors, and is necessary for development of specific mesoderm and endoderm-derived organs including the heart, liver, and pancreas. Furthermore, reactivation or silencing of the Gata6 locus has been associated with certain types of cancer affecting endodermal organs. RESULTS: We have generated a Gata6(H2B-Venus) knock-in reporter mouse allele for the purpose of labeling GATA6-expressing cells with a bright nuclear-localized fluorescent marker that is suitable for live imaging at single-cell resolution. CONCLUSIONS: Expression of the Venus reporter was characterized starting from embryonic stem (ES) cells, through mouse embryos and adult animals. The Venus reporter was not expressed in ES cells, but was activated upon endoderm differentiation. Gata6(H2B-Venus/H2B-Venus) homozygous embryos did not express GATA6 protein and failed to specify the primitive endoderm in the blastocyst. However, null blastocysts continued to express high levels of Venus in the absence of GATA6 protein, suggesting that early Gata6 transcription is independent of GATA6 protein expression. At early post-implantation stages of embryonic development, there was a strong correlation of Venus with endogenous GATA6 protein in endoderm and mesoderm progenitors, then later in the heart, midgut, and hindgut. However, there were discrepancies in reporter versus endogenous protein expression in certain cells, such as the body wall and endocardium. During organogenesis, detection of Venus in specific organs recapitulated known sites of endogenous GATA6 expression, such as in the lung bud epithelium, liver, pancreas, gall bladder, stomach epithelium, and vascular endothelium. In adults, Venus was observed in the lungs, pancreas, liver, gall bladder, ovaries, uterus, bladder, skin, adrenal glands, small intestine and corpus region of the stomach. Overall, Venus fluorescent protein under regulatory control of the Gata6 locus was expressed at levels that were easily visualized directly and could endure live and time-lapse imaging techniques. Venus is co-expressed with endogenous GATA6 throughout development to adulthood, and should provide an invaluable tool for examining the status of the Gata6 locus during development, as well as its silencing or reactivation in cancer or other disease states.

Anatomy of a blastocyst: cell behaviors driving cell fate choice and morphogenesis in the early mouse embryo.

The preimplantation period of mouse early embryonic development is devoted to the specification of two extraembryonic tissues and their spatial segregation from the pluripotent epiblast. During this period two cell fate decisions are made while cells gradually lose their totipotency. The first fate decision involves the segregation of the extraembryonic trophectoderm (TE) lineage from the inner cell mass (ICM); the second occurs within the ICM and involves the segregation of the extraembryonic primitive endoderm (PrE) lineage from the pluripotent epiblast (EPI) lineage, which eventually gives rise to the embryo proper. Multiple determinants, such as differential cellular properties, signaling cues and the activity of transcriptional regulators, influence lineage choice in the early embryo. Here, we provide an overview of our current understanding of the mechanisms governing these cell fate decisions ensuring proper lineage allocation and segregation, while at the same time providing the embryo with an inherent flexibility to adjust when perturbed.

A bright single-cell resolution live imaging reporter of Notch signaling in the mouse.

BACKGROUND: Live imaging provides an essential methodology for understanding complex and dynamic cell behaviors and their underlying molecular mechanisms. Genetically-encoded reporter expressing mouse strains are an important tool for use in live imaging experiments. Such reporter strains can be engineered by placing cis-regulatory elements of interest to direct the expression of desired reporter genes. If these cis-regulatory elements are downstream targets, and thus activated as a consequence of signaling pathway activation, such reporters can provide read-outs of the signaling status of a cell. The Notch signaling pathway is an evolutionary conserved pathway operating in multiple developmental processes as well as being the basis for several congenital diseases. The transcription factor CBF1 is a central evolutionarily conserved component of the Notch signaling pathway. It binds the active form of the Notch receptor (NICD) and subsequently binds to cis-regulatory regions (CBF1 binding sites) in the promoters of Notch responsive genes. In this way, CBF1 binding sites represent a good target for the design of a Notch signaling reporter. RESULTS: To generate a single-cell resolution Notch signaling reporter, we used a CBF responsive element to direct the expression of a nuclear-localized fluorescent protein. To do this, we linked 4 copies of a consensus CBF1 binding site to the basal simian virus 40 (SV40) promoter, placed this cassette in front of a fluorescent protein fusion comprising human histone H2B linked to the yellow fluorescent protein (YFP) Venus, one of the brightest available YFPs. We used the CBF:H2B-Venus construct to generate both transgenic embryonic mouse stem (ES) cell lines and a strain of transgenic mice that would report Notch signaling activity. CONCLUSION: By using multiple CBF1 binding sites together with a subcellular-localized, genetically-encoded fluorescent protein, H2B-Venus, we have generated a transgenic strain of mice that faithfully recapitulates Notch signaling at single-cell resolution. This is the first mouse reporter strain in which individual cells transducing a Notch signal can be visualized. The improved resolution of this reporter makes it ideal for live imaging developmental processes regulated by the Notch signaling pathway as well as a short-term lineage tracer of Notch expressing cells due to the perdurance of the fluorescent reporter. Taken together, the CBF:H2B-Venus mouse strain is a unique tool to study and understand the morphogenetic events regulated by the Notch signaling pathway.

3D-Image analysis platform monitoring relocation of pluripotency genes during reprogramming.

Nuclear organization of chromatin is an important level of genome regulation with positional changes of genes occurring during reprogramming. Inherent variability of biological specimens, wide variety of sample preparation and imaging conditions, though pose significant challenges to data analysis and comparison. Here, we describe the development of a computational image analysis toolbox overcoming biological variability hurdles by a novel single cell randomizing normalization. We performed a comparative analysis of the relationship between spatial positioning of pluripotency genes with their genomic activity and determined the degree of similarity between fibroblasts, induced pluripotent stem cells and embryonic stem cells. Our analysis revealed a preferred positioning of actively transcribed Sox2, Oct4 and Nanog away from the nuclear periphery, but not from pericentric heterochromatin. Moreover, in the silent state, we found no common nuclear localization for any of the genes. Our results suggest that the surrounding gene density hinders relocation from an internal nuclear position. Altogether, our data do not support the hypothesis that the nuclear periphery acts as a general transcriptional silencer, rather suggesting that internal nuclear localization is compatible with expression in pluripotent cells but not sufficient for expression in mouse embryonic fibroblasts. Thus, our computational approach enables comparative analysis of topological relationships in spite of stark morphological variability typical of biological data sets.

Human pluripotent stem cell-derived respiratory airway progenitors generate alveolar epithelial cells and recapitulate features of idiopathic pulmonary fibrosis.

Human lungs contain unique cell populations in distal respiratory airways (RAs). These populations accumulate in patients with lung injury, chronic obstructive pulmonary disease (COPD) and idiopathic pulmonary fibrosis (IPF). Their lineage potentials and roles are unknown, however. As they are absent in rodents, deeper understanding of these cells requires a human in vitro model. Here we report the generation from human pluripotent stem cells (hPSCs) of expandable spheres (induced respiratory airway progenitors (iRAPs)) consisting of all RA-associated cell types. iRAPs could differentiate into type 1 (AT1) and type 2 alveolar (AT2) epithelial cells in defined conditions, showing that alveolar cells can be derived from RAs. iRAPs with deletion of HPS1, which causes pulmonary fibrosis in humans, display defects that are hallmarks of IPF, indicating involvement of intrinsic dysfunction of RA-associated cells in IPF. iRAPs thus provide a model to gain insight into human lung regeneration and into pathogenesis of IPF.

A distal lung organoid model to study interstitial lung disease, viral infection and human lung development.

Organoids have been an exciting advancement in stem cell research. Here we describe a strategy for directed differentiation of human pluripotent stem cells into distal lung organoids. This protocol recapitulates lung development by sequentially specifying human pluripotent stem cells to definitive endoderm, anterior foregut endoderm, ventral anterior foregut endoderm, lung bud organoids and finally lung organoids. The organoids take ~40 d to generate and can be maintained more than 180 d, while progressively maturing up to a stage consistent with the second trimester of human gestation. They are unique because of their branching morphology, the near absence of non-lung endodermal lineages, presence of mesenchyme and capacity to recapitulate interstitial lung diseases. This protocol can be performed by anyone familiar with cell culture techniques, is conducted in serum-free conditions and does not require lineage-specific reporters or enrichment steps. We also provide a protocol for the generation of single-cell suspensions for single-cell RNA sequencing.

HIV-1 activates oxidative phosphorylation in infected CD4 T cells in a human tonsil explant model.

Introduction: Human immunodeficiency virus type 1 (HIV-1) causes a chronic, incurable infection leading to immune activation and chronic inflammation in people with HIV-1 (PWH), even with virologic suppression on antiretroviral therapy (ART). The role of lymphoid structures as reservoirs for viral latency and immune activation has been implicated in chronic inflammation mechanisms. Still, the specific transcriptomic changes induced by HIV-1 infection in different cell types within lymphoid tissue remain unexplored. Methods: In this study, we utilized human tonsil explants from healthy human donors and infected them with HIV-1 ex vivo. We performed single-cell RNA sequencing (scRNA-seq) to analyze the cell types represented in the tissue and to investigate the impact of infection on gene expression profiles and inflammatory signaling pathways. Results: Our analysis revealed that infected CD4+ T cells exhibited upregulation of genes associated with oxidative phosphorylation. Furthermore, macrophages exposed to the virus but uninfected showed increased expression of genes associated with the NLRP3 inflammasome pathway. Discussion: These findings provide valuable insights into the specific transcriptomic changes induced by HIV-1 infection in different cell types within lymphoid tissue. The activation of oxidative phosphorylation in infected CD4+ T cells and the proinflammatory response in macrophages may contribute to the chronic inflammation observed in PWH despite ART. Understanding these mechanisms is crucial for developing targeted therapeutic strategies to eradicate HIV-1 infection in PWH.

Non-additive effects of schizophrenia risk genes reflect convergent downstream function.

Genetic studies of schizophrenia (SCZ) reveal a complex polygenic risk architecture comprised of hundreds of risk variants, the majority of which are common in the population at-large and confer only modest increases in disorder risk. Precisely how genetic variants with individually small predicted effects on gene expression combine to yield substantial clinical impacts in aggregate is unclear. Towards this, we previously reported that the combinatorial perturbation of four SCZ risk genes ("eGenes", whose expression is regulated by common variants) resulted in gene expression changes that were not predicted by individual perturbations, being most non-additive among genes associated with synaptic function and SCZ risk. Now, across fifteen SCZ eGenes, we demonstrate that non-additive effects are greatest within groups of functionally similar eGenes. Individual eGene perturbations reveal common downstream transcriptomic effects ("convergence"), while combinatorial eGene perturbations result in changes that are smaller than predicted by summing individual eGene effects ("sub-additive effects"). Unexpectedly, these convergent and sub-additive downstream transcriptomic effects overlap and constitute a large proportion of the genome-wide polygenic risk score, suggesting that functional redundancy of eGenes may be a major mechanism underlying non-additivity. Single eGene perturbations likewise fail to predict the magnitude or directionality of cellular phenotypes resulting from combinatorial perturbations. Overall, our results indicate that polygenic risk cannot be extrapolated from experiments testing one risk gene at a time and must instead be empirically measured. By unravelling the interactions between complex risk variants, it may be possible to improve the clinical utility of polygenic risk scores through more powerful prediction of symptom onset, clinical trajectory, and treatment response, or to identify novel targets for therapeutic intervention.

Chromatin profiling in human neurons reveals aberrant roles for histone acetylation and BET family proteins in schizophrenia.

Schizophrenia (SZ) is a psychiatric disorder with complex genetic risk dictated by interactions between hundreds of risk variants. Epigenetic factors, such as histone posttranslational modifications (PTMs), have been shown to play critical roles in many neurodevelopmental processes, and when perturbed may also contribute to the precipitation of disease. Here, we apply an unbiased proteomics approach to evaluate combinatorial histone PTMs in human induced pluripotent stem cell (hiPSC)-derived forebrain neurons from individuals with SZ. We observe hyperacetylation of H2A.Z and H4 in neurons derived from SZ cases, results that were confirmed in postmortem human brain. We demonstrate that the bromodomain and extraterminal (BET) protein, BRD4, is a bona fide 'reader' of H2A.Z acetylation, and further provide evidence that BET family protein inhibition ameliorates transcriptional abnormalities in patient-derived neurons. Thus, treatments aimed at alleviating BET protein interactions with hyperacetylated histones may aid in the prevention or treatment of SZ.

Analysis framework and experimental design for evaluating synergy-driving gene expression.

The mechanisms by which genetic risk variants interact with each other, as well as environmental factors, to contribute to complex genetic disorders remain unclear. We describe in detail our recently published approach to resolve distinct additive and synergistic transcriptomic effects after combinatorial manipulation of genetic variants and/or chemical perturbagens. Although first developed for CRISPR-based perturbation studies of isogenic human induced pluripotent stem cell-derived neurons, our methodology can be broadly applied to any RNA sequencing dataset, provided that raw read counts are available. Whereas other differential expression analyses reveal the effect of individual perturbations, here we specifically query interactions between two or more perturbagens, resolving the extent of non-additive (synergistic) interactions between perturbations. We discuss the careful experimental design required to resolve synergistic effects and considerations of statistical power and how to quantify observed synergy between experiments. Additionally, we speculate on potential future applications and explore the obvious limitations of this approach. Overall, by interrogating the effect of independent factors, alone and in combination, our analytic framework and experimental design facilitate the discovery of convergence and synergy downstream of gene and/or treatment perturbations hypothesized to contribute to complex diseases. We think that this protocol can be successfully applied by any scientist with bioinformatic skills and basic proficiency in the R programming language. Our computational pipeline ( https://github.com/nadschro/synergy-analysis ) is straightforward, does not require supercomputing support and can be conducted in a single day upon completion of RNA sequencing experiments.

Using the dCas9-KRAB system to repress gene expression in hiPSC-derived NGN2 neurons.

We describe a CRISPR inhibition (CRISPRi) protocol to repress endogenous gene expression (e.g., ATP6V1A) in human induced pluripotent stem cell-derived NGN2-induced glutamatergic neurons. CRISPRi enables efficient and precise gene repression of one or multiple target genes via delivering gRNA(s) to direct a dCas9-KRAB fusion protein to the gene(s) of interest. This protocol can also be adapted for gene activation and high-throughput gene manipulation, allowing assessment of the transcriptomic and phenotypic impact of candidate gene(s) associated with neurodevelopment or brain disease. For complete details on the use and execution of this protocol, please refer to Ho et al. (2017) and Wang et al. (2021).

Transformative Network Modeling of Multi-omics Data Reveals Detailed Circuits, Key Regulators, and Potential Therapeutics for Alzheimer's Disease.

To identify the molecular mechanisms and novel therapeutic targets of late-onset Alzheimer's Disease (LOAD), we performed an integrative network analysis of multi-omics profiling of four cortical areas across 364 donors with varying cognitive and neuropathological phenotypes. Our analyses revealed thousands of molecular changes and uncovered neuronal gene subnetworks as the most dysregulated in LOAD. ATP6V1A was identified as a key regulator of a top-ranked neuronal subnetwork, and its role in disease-related processes was evaluated through CRISPR-based manipulation in human induced pluripotent stem cell-derived neurons and RNAi-based knockdown in Drosophila models. Neuronal impairment and neurodegeneration caused by ATP6V1A deficit were improved by a repositioned compound, NCH-51. This study provides not only a global landscape but also detailed signaling circuits of complex molecular interactions in key brain regions affected by LOAD, and the resulting network models will serve as a blueprint for developing next-generation therapeutic agents against LOAD.

Cell Type-Specific In Vitro Gene Expression Profiling of Stem Cell-Derived Neural Models.

Genetic and genomic studies of brain disease increasingly demonstrate disease-associated interactions between the cell types of the brain. Increasingly complex and more physiologically relevant human-induced pluripotent stem cell (hiPSC)-based models better explore the molecular mechanisms underlying disease but also challenge our ability to resolve cell type-specific perturbations. Here, we report an extension of the RiboTag system, first developed to achieve cell type-restricted expression of epitope-tagged ribosomal protein (RPL22) in mouse tissue, to a variety of in vitro applications, including immortalized cell lines, primary mouse astrocytes, and hiPSC-derived neurons. RiboTag expression enables depletion of up to 87 percent of off-target RNA in mixed species co-cultures. Nonetheless, depletion efficiency varies across independent experimental replicates, particularly for hiPSC-derived motor neurons. The challenges and potential of implementing RiboTags in complex in vitro cultures are discussed.

Neuronal impact of patient-specific aberrant NRXN1α splicing.

NRXN1 undergoes extensive alternative splicing, and non-recurrent heterozygous deletions in NRXN1 are strongly associated with neuropsychiatric disorders. We establish that human induced pluripotent stem cell (hiPSC)-derived neurons well represent the diversity of NRXN1α alternative splicing observed in the human brain, cataloguing 123 high-confidence in-frame human NRXN1α isoforms. Patient-derived NRXN1+/- hiPSC-neurons show a greater than twofold reduction in half of the wild-type NRXN1α isoforms and express dozens of novel isoforms from the mutant allele. Reduced neuronal activity in patient-derived NRXN1+/- hiPSC-neurons is ameliorated by overexpression of individual control isoforms in a genotype-dependent manner, whereas individual mutant isoforms decrease neuronal activity levels in control hiPSC-neurons. In a genotype-dependent manner, the phenotypic impact of patient-specific NRXN1+/- mutations can occur through a reduction in wild-type NRXN1α isoform levels as well as the presence of mutant NRXN1α isoforms.

Synergistic effects of common schizophrenia risk variants.

The mechanisms by which common risk variants of small effect interact to contribute to complex genetic disorders are unclear. Here, we apply a genetic approach, using isogenic human induced pluripotent stem cells, to evaluate the effects of schizophrenia (SZ)-associated common variants predicted to function as SZ expression quantitative trait loci (eQTLs). By integrating CRISPR-mediated gene editing, activation and repression technologies to study one putative SZ eQTL (FURIN rs4702) and four top-ranked SZ eQTL genes (FURIN, SNAP91, TSNARE1 and CLCN3), our platform resolves pre- and postsynaptic neuronal deficits, recapitulates genotype-dependent gene expression differences and identifies convergence downstream of SZ eQTL gene perturbations. Our observations highlight the cell-type-specific effects of common variants and demonstrate a synergistic effect between SZ eQTL genes that converges on synaptic function. We propose that the links between rare and common variants implicated in psychiatric disease risk constitute a potentially generalizable phenomenon occurring more widely in complex genetic disorders.