Summary Amyotrophic lateral sclerosis (ALS) is a devastating condition with very limited treatment options. It is a heterogeneous disease with complex genetics and unclear etiology, making the discovery of disease-modifying interventions very challenging. To discover novel mechanisms underlying ALS, we leverage a unique platform that combines isogenic, induced pluripotent stem cell (iPSC)-derived models of disease-causing mutations with rich phenotyping via high-content imaging and deep learning models. We introduced eight mutations that cause familial ALS (fALS) into multiple donor iPSC lines, and differentiated them into motor neurons to create multiple isogenic pairs of healthy (wild-type) and sick (mutant) motor neurons. We collected extensive high-content imaging data and used machine learning (ML) to process the images, segment the cells, and learn phenotypes. Self-supervised ML was used to create a concise embedding that captured significant, ALS-relevant biological information in these images. We demonstrate that ML models trained on core cell morphology alone can accurately predict TDP-43 mislocalization, a known phenotypic feature related to ALS. In addition, we were able to impute RNA expression from these image embeddings, in a way that elucidates molecular differences between mutants and wild-type cells. Finally, predictors leveraging these embeddings are able to distinguish between mutant and wild-type both within and across donors, defining cellular, ML-derived disease models for diverse fALS mutations. These disease models are the foundation for a novel screening approach to discover disease-modifying targets for familial ALS.
Amyotrophic Lateral Sclerosis (ALS) is a fatal neurodegenerative disorder characterised by a progressive loss of motor function. The eponymous spinal sclerosis observed at autopsy is the result of the degeneration of extratelencephalic neurons, Betz cells (ETNs, Cortico-Spinal Motor Neuron). It remains unclear why this neuronal subtype is selectively affected. To understand the unique molecular properties that sensitise these cells to ALS, we performed RNA sequencing of 79,169 single nuclei from cortices of patients and controls. In unaffected individuals, we found that expression of ALS risk genes was significantly enriched inTHY1+-ETNs and not in other cell types. In patients, these genetic risk factors, as well as genes involved in protein homeostasis and stress responses, were significantly induced in a wide collection of ETNs, but not in neurons with more superficial identities. Examination of oligodendroglial and microglial nuclei revealed patient-specific changes that were at least in part a response to alterations in neurons: downregulation of myelinating genes in oligodendrocytes and upregulation of a reactive state connected to endo-lysosomal pathways in microglia. Our findings suggest that the selective vulnerability of extratelencephalic neurons is partly connected to their intrinsic molecular properties sensitising them to genetics and mechanisms of degeneration.Graphical abstract and working modelOur study highlights cell type specific changes in premotor/motor cortex of sporadic ALS patients. Specifically, we identify upregulation of synaptic molecules in excitatory neurons of upper cortical layers, interestingly correlating to hyperexcitability phenotypes seen in patients. Moreover, excitatory neurons of the deeper layers of the cortex, that project to the spinal cord and are most affected by the disease, show higher levels of cellular stresses than other neuronal types. Correspondently, oligodendrocytes transition from a highly myelinating state to a more neuronally engaged state, probably to counteract stressed phenotypes seen in excitatory neurons. At the same time, microglia show a reactive state with specific upregulation of endo-lysosomal pathways.
The maturation of neurons and the development of synapses, although emblematic of neurons, also relies on interactions with astrocytes and other glia. Here, to study the role of glia-neuron interactions, we analyze the transcriptomes of human pluripotent stem cell (hPSC)-derived neurons, from 80 human donors, that were cultured with or without contact with glial cells. We find that the presence of astrocytes enhances synaptic gene-expression programs in neurons when in physical contact with astrocytes. These changes in neurons correlate with increased expression, in the cocultured glia, of genes that encode synaptic cell adhesion molecules. Both the neuronal and astrocyte gene-expression programs are enriched for genes associated with schizophrenia risk. Our results suggest that astrocyte-expressed genes with synaptic functions are associated with stronger expression of synaptic genetic programs in neurons, and they suggest a potential role for astrocyte-neuron interactions in schizophrenia.
Human pluripotent stem cells (hPSCs) are a powerful tool for disease modeling of hard-to-access tissues (such as the brain). Current protocols either direct neuronal differentiation with small molecules or use transcription-factor-mediated programming. In this study, we couple overexpression of transcription fac-tor Neurogenin2 (Ngn2) with small molecule patterning to differentiate hPSCs into lower induced motor neurons (liMoNes/liMNs). This approach induces canonical MN markers including MN-specific Hb9/ MNX1 in more than 95% of cells. liMNs resemble bona fide hPSC-derived MN, exhibit spontaneous elec-trical activity, express synaptic markers, and can contact muscle cells in vitro. Pooled, multiplexed sin-gle-cell RNA sequencing on 50 hPSC lines reveals reproducible populations of distinct subtypes of cer-vical and brachial MNs that resemble their in vivo, embryonic counterparts. Combining small molecule patterning with Ngn2 overexpression facilitates high-yield, reproducible production of disease-relevant MN subtypes, which is fundamental in propelling our knowledge of MN biology and its disruption in dis-ease.
Human genome variation contributes to diversity in neurodevelopmental outcomes and vulnerabilities; recognizing the underlying molecular and cellular mechanisms will require scalable approaches. Here, we describe a "cell village" experimental platform we used to analyze genetic, molecular, and phenotypic heterogeneity across neural progenitor cells from 44 human donors cultured in a shared in vitro environment using algorithms (Dropulation and Census-seq) to assign cells and phenotypes to individual donors. Through rapid induction of human stem cell-derived neural progenitor cells, measurements of natural genetic variation, and CRISPR-Cas9 genetic perturbations, we identified a common variant that regulates antiviral IFITM3 expression and explains most inter-individual variation in susceptibility to the Zika virus. We also detected expression QTLs corresponding to GWAS loci for brain traits and discovered novel disease-relevant regulators of progenitor proliferation and differentiation such as CACHD1. This approach provides scalable ways to elucidate the effects of genes and genetic variation on cellular phenotypes.
SUMMARYVariation in the human genome contributes to abundant diversity in human traits and vulnerabilities, but the underlying molecular and cellular mechanisms are not yet known, and will need scalable approaches to accelerate their recognition. Here, we advanced and applied an experimental platform that analyzes genetic, molecular, and phenotypic heterogeneity across cells from very many human donors cultured in a single, shared in vitro environment, with algorithms (Dropulation and Census-seq) for assigning phenotypes to individual donors. We used natural genetic variation and synthetic (CRISPR-Cas9) genetic perturbations to analyze the vulnerability of neural progenitor cells to infection with Zika virus. These analyses identified a common variant in the antiviral IFITM3 gene that regulated IFITM3 expression and explained most inter-individual variation in NPCs’ susceptibility to Zika virus infectivity. These and other approaches could provide scalable ways to recognize the impact of genes and genetic variation on cellular phenotypes.HIGHLIGHTSMeasuring cellular phenotypes in iPSCs and hPSC-derived NPCs from many donorsEffects of donor sex, cell source, genetic and other variables on hPSC RNA expressionNatural genetic variation and synthetic perturbation screens both identify IFITM3 in NPC susceptibility to Zika virusA common genetic variant in IFITM3 explains most inter-individual variation in NPC susceptibility to Zika virus
SummaryTens of thousands of genetic variants shape human phenotypes, mostly by unknown cellular mechanisms. Here we describe Census-seq, a way to measure cellular phenotypes in cells from many people simultaneously. Analogous to pooled CRISPR screens but for natural variation, Census-seq associates cellular phenotypes to donors’ genotypes by quantifying the presence of each donor’s DNA in cell “villages” before and after sorting or selection for cellular traits of interest. Census-seq enables population-scale cell-biological phenotyping with low cost and high internal control. We demonstrate Census-seq through investigation of genetic effects on the SMN protein whose deficiency underlies spinal muscular atrophy (SMA). Census-seq quantified and mapped effects of many common alleles on SMN protein levels and response to SMN-targeted therapeutics, including a common, cryptic non-responder allele. We provide tools enabling population-scale cell experiments and explain how Census-seq can be used to map genetic effects on diverse cell phenotypes.Abstract FigureHighlightsCensus-seq reveals how inherited genetic variation affects cell phenotypesGenetic analysis of cellular traits in cell villages of >100 donorsCharacterizing human alleles that shape SMN protein expression and drug responsesDevelopment of protocols and software to enable cellular population genetics
Neural progenitor cells (NPCs) are essential to brain development and their dysfunction is linked to several disorders, including autism, Zika Virus Congenital Syndrome, and cancer. Understanding of these conditions has been improved by advancements with stem cell-derived NPC models. However, current differentiation methods require many days or weeks to generate NPCs and show variability in efficacy among cell lines. Here, we describe human Stem cell-derived NGN2-accelerated Progenitor cells (SNaPs), which are produced in only 48 hours. SNaPs express canonical forebrain NPC protein markers, are proliferative, multipotent, and like other human NPCs, are susceptible to Zika-mediated death. We further demonstrate SNaPs are valuable for large-scale investigations of genetic and environmental influencers of neurodevelopment by deploying them for genome-wide CRISPR-Cas9 screens. Our studies expand knowledge of NPCs by identifying known and novel Zika host factors, as well as new regulators of NPC proliferation validated by re-identification of the autism spectrum gene PTEN .
The authors surveyed whole-exome and RNA-sequencing data from 252 unique pluripotent stem cell lines, some of which are in the pipeline for clinical use, and found that approximately 5% of cell lines had acquired mutations in the TP53 gene that allow mutant cells to rapidly outcompete non-mutant cells, but do not prevent differentiation. Copy number variants at particular genomic locations have been shown to arise in human pluripotent stem cells (hPSCs) under certain culture conditions, but the extent of acquired mutations in such culture remains to be determined. Kevin Eggan and colleagues surveyed the exomes of 140 human embryonic stem cell (hESC) lines, some of which are in the pipeline for clinical use.They identified mosaic mutations in the TP53 gene in a subset of cells for five unrelated hESC lines and show that the cells carrying the mutations outcompeted the non-mutant cells and could readily differentiate. Similar mutations were also identified by mining published datasets for an additional 14 hESC lines and more than 100 human induced PSC lines. The study highlights the need for in-depth characterization of cells derived from hPSCs before their use in the clinic. Human pluripotent stem cells (hPS cells) can self-renew indefinitely, making them an attractive source for regenerative therapies. This expansion potential has been linked with the acquisition of large copy number variants that provide mutated cells with a growth advantage in culture1,2,3. The nature, extent and functional effects of other acquired genome sequence mutations in cultured hPS cells are not known. Here we sequence the protein-coding genes (exomes) of 140 independent human embryonic stem cell (hES cell) lines, including 26 lines prepared for potential clinical use4. We then apply computational strategies for identifying mutations present in a subset of cells in each hES cell line5. Although such mosaic mutations were generally rare, we identified five unrelated hES cell lines that carried six mutations in the TP53 gene that encodes the tumour suppressor P53. The TP53 mutations we observed are dominant negative and are the mutations most commonly seen in human cancers. We found that the TP53 mutant allelic fraction increased with passage number under standard culture conditions, suggesting that the P53 mutations confer selective advantage. We then mined published RNA sequencing data from 117 hPS cell lines, and observed another nine TP53 mutations, all resulting in coding changes in the DNA-binding domain of P53. In three lines, the allelic fraction exceeded 50%, suggesting additional selective advantage resulting from the loss of heterozygosity at the TP53 locus. As the acquisition and expansion of cancer-associated mutations in hPS cells may go unnoticed during most applications, we suggest that careful genetic characterization of hPS cells and their differentiated derivatives be carried out before clinical use.
The coatomer module of the nuclear pore complex borders the cylinder-like nuclear pore-membrane domain of the nuclear envelope. In evolution, a single coatomer module increases in size from hetero-heptamer ( Saccharomyces cerevisiae ) to hetero-octamer ( Schizosaccharomyces pombe ) to hetero-nonamer (Metazoa). Notably, the heptamer–octamer transition proceeds through the acquisition of the nucleoporin Nup37. How Nup37 contacts the heptamer remained unknown. Using recombinant nucleoporins, we show that Sp-Nup37 specifically binds the Sp-Nup120 member of the hetero-heptamer but does not bind an Sc-Nup120 homolog. To elucidate the Nup37–Nup120 interaction at the atomic level, we carried out crystallographic analyses of Sp-Nup37 alone and in a complex with an N-terminal, ∼110-kDa fragment of Sp-Nup120 comprising residues 1–950. Corroborating structural predictions, we determined that Nup37 folds into a seven-bladed β-propeller. Several disordered surface regions of the Nup37 β-propeller assume structure when bound to Sp-Nup120. The N-terminal domain of Sp-Nup120 1–950 also folds into a seven-bladed propeller with a markedly protruding 6D–7A insert and is followed by a contorted helical domain. Conspicuously, this 6D–7A insert contains an extension of 50 residues which also is highly conserved in Metazoa but is absent in Sc-Nup120. Strikingly, numerous contacts with the Nup37 β-propeller are located on this extension of the 6D–7A insert. Another contact region is situated toward the end of the helical region of Sp-Nup120 1–950 . Our findings provide information about the evolution and the assembly of the coatomer module of the nuclear pore complex.
Nuclear pore complexes (NPCs) control the movement of molecules across the nuclear envelope (NE). We investigated the molecular interactions that exist at the interface between the NPC scaffold and the pore membrane. We show that key players mediating these interactions in mammalian cells are the nucleoporins Nup155 and Nup160. Nup155 depletion massively alters NE structure, causing a dramatic decrease in NPC numbers and the improper targeting of membrane proteins to the inner nuclear membrane. The role of Nup155 in assembly is likely closely linked to events at the membrane as we show that Nup155 interacts with pore membrane proteins Pom121 and NDC1. Furthermore, we demonstrate that the N terminus of Pom121 directly binds the β-propeller regions of Nup155 and Nup160. We propose a model in which the interactions of Pom121 with Nup155 and Nup160 are predicted to assist in the formation of the nuclear pore and the anchoring of the NPC to the pore membrane.
Oxytocin (OT) is a potent uterine agonist. Its receptor (OTR) is a G protein-coupled receptor that is downregulated by prolonged exposure to OT. We hypothesized that activation of PKC mediated this OT-induced decrease in OTR expression. Diminished PKC activity in late pregnancy could underlie the increased expression of uterine OTR preceding labor onset. Using cell cultures of transformed human uterine myocytes, we determined the effects of PKC agonists and antagonists on the expression of OTR. We also explored the effects of overexpression of activator protein-1 (AP-1, a mediator of many PKC- and phorbol ester-induced effects) using adenoviral expression vectors for the AP-1 subunits c-Jun and c-Fos. Stimulation of PKC using the phorbol ester 12-O-tetradecanoylphorbol 13-acetate caused a rapid, significant (P < or = 0.05) increase in c-Jun and c-Fos concentrations but a significant decrease in mRNA for OTR within 6 h followed by a significant decrease in OT binding by 24 h. Adenoviral infection of the cells with expression vectors for c-Jun and c-Fos increased the AP-1 subunits but had no effect on OTR expression. Furthermore, there were no changes in c-Fos or c-Jun levels in human intrauterine tissues around the time of labor onset, as measured by Western analyses. We conclude that phorbol ester treatment decreases OTR expression, likely through a mechanism that does not involve AP-1.
Compromised fetal growth impairs vascular function; however, it is unclear whether chronic hypoxia in utero affects adult endothelial function. We hypothesized that maternal hypoxia (H, 12% O2, n= 9) or nutrient restriction (NR, 40% of control, n= 7) imposed from day 15-21 pregnancy in rats would impair endothelial function in adult male offspring (relative to control, C, n= 10). Using a wire myograph, endothelium-dependent relaxation in response to methacholine was assessed in small mesenteric arteries from 4- and 7-month-old (mo) male offspring. Nitric oxide (NO) mediation of endothelium-dependent relaxation was evaluated using N(omega)-nitro-L-arginine methyl ester (L-NAME; NO synthase inhibitor). Observed differences in the NO pathway at 7 months were investigated using exogenous superoxide dismutase (SOD) to reduce NO scavenging, and sodium nitroprusside (SNP; NO donor) to assess smooth muscle sensitivity to NO. Sensitivity to methacholine-induced endothelium-dependent relaxation was reduced in H offspring at 4 months (P < 0.05), but was not different among groups at 7 months. L-NAME reduced methacholine sensitivity in C (P < 0.01), H (P < 0.01) and NR (P < 0.05) offspring at 4 months, but at 7 months L-NAME reduced sensitivity in C (P < 0.05), tended to in NR (P= 0.055) but had no effect in H offspring. SOD did not alter sensitivity to methacholine in C, but increased sensitivity in H offspring (P < 0.01). SNP responses did not differ among groups. In summary, prenatal hypoxia, but not nutrient restriction impaired endothelium-dependent relaxation at 4 months, and reduced NO mediation of endothelial function at 7 months, in part through reduced NO bio-availability. Distinct effects following reduced maternal oxygen versus nutrition suggest that decreased oxygen supply during fetal life may specifically impact adult vascular function.
Increased concentrations of IL-1beta and TNF-alpha have been associated with parturition. However, the role of these cytokines is unknown. Before parturition, the uterus undergoes a process of activation, during which there are significant changes in expression of genes associated with increased uterine contractility, including the receptors for oxytocin (OT) and prostaglandin (PG)F(2alpha) (FP), PGH(2) synthase isoform 2 (PGHS2), the gap junction protein connexin-43 (Cx-43), and the inducible isoform of nitric oxide synthase (iNOS). To determine whether IL-1beta or TNF-alpha was part of the causal mechanism for increased uterine contractions, we placed osmotic pumps infusing IL-1beta or TNF-alpha into the peritoneal cavity of late pregnant rats (gestation day 19) and measured the effects on uterine contractility and on the uterine concentrations of mRNA for the contraction-associated genes 24 h later. Maternal serum concentrations of IL-1beta and TNF-alpha were increased significantly. By day 21, the control animals had significant increases (P < or = 0.05) in mRNA for OT, FP, PGHS2, and Cx-43, a decrease (P < or = 0.05) in iNOS, and an increase (P < or = 0.05) in uterine sensitivity and responsiveness to OT. Infusion of IL-1beta or TNF-alpha had no effect on uterine contractility or on expression of the activation-associated genes. We conclude that intraperitoneal infusion of IL-1beta or TNF-alpha resulting in significantly increased maternal serum cytokine levels does not cause uterine activation. The role of proinflammatory cytokines in the mechanism of parturition remains unclear.