Search

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.

The morphology of cells is dynamic and mediated by genetic and environmental factors. Characterizing how genetic variation impacts cell morphology can provide an important link between disease association and cellular function. Here, we combined genomic and high-content imaging approaches on iPSCs from 297 unique donors to investigate the relationship between genetic variants and cellular morphology to map what we term cell morphological quantitative trait loci (cmQTLs). We identified novel associations between rare protein altering variants inWASF2, TSPAN15, andPRLRwith several morphological traits related to cell shape, nucleic granularity, and mitochondrial distribution. Knockdown of these genes by CRISPRi confirmed their role in cell morphology. Analysis of common variants yielded one significant association and nominated over 300 variants with suggestive evidence (P-6) of association with one or more morphology traits. Our results showed that, similar to other molecular phenotypes, morphological profiling can yield insight about the function of genes and variants.

Within the vertebrate neocortex and other telencephalic structures, molecularly-defined neurons tend to segregate at first order into inhibitory (GABAergic) and excitatory (glutamatergic) types. We used single-nucleus RNA sequencing, analyzing over 2.4 million brain cells sampled from 16 locations in a primate (the common marmoset) to ask whether (1) neurons generally segregate by neurotransmitter status, and (2) neurons expressing the same neurotransmitters share additional molecular features in common, beyond the few genes directly responsible for neurotransmitter synthesis and release. We find the answer to both is “no”: there is a remarkable degree of transcriptional similarity between GABAergic and glutamatergic neurons found in the same brain structure, and there is generally little in common between glutamatergic neurons residing in phylogenetically divergent brain structures. The origin effect is permanent: we find that cell types that cross cephalic boundaries in development retain the transcriptional identities of their birthplaces. GABAergic interneurons, which migrate widely, follow highly specialized and distinct distributions in striatum and neocortex. We use interneuron-restricted AAVs to reveal the morphological diversity of molecularly defined types. Our analyses expose how lineage and functional class sculpt the transcriptional identity and biodistribution of primate neurons.One-Sentence SummaryPrimate neurons are primarily imprinted by their region of origin, more so than by their functional identity.

Many regions in the human genome vary in length among individuals due to variable numbers of tandem repeats (VNTRs). We recently showed that protein-coding VNTRs underlie some of the strongest known genetic associations with diverse phenotypes. Here, we assessed the phenotypic impact of VNTRs genome-wide, 99% of which lie in non-coding regions. We applied a statistical imputation approach to estimate the lengths of 9,561 autosomal VNTR loci in 418,136 unrelated UK Biobank participants. Association and statistical fine-mapping analyses identified 107 VNTR-phenotype associations (involving 58 VNTRs) that were assigned a high probability of VNTR causality (PIP≥0.5). Non-coding VNTRs atTMCO1andEIF3Happeared to generate the largest known contributions of common human genetic variation to risk of glaucoma and colorectal cancer, respectively. Each of these two VNTRs associated with a >2- fold risk range across individuals. These results reveal a substantial and previously unappreciated role of non-coding VNTRs in human health.

Marmosets, especially Callithrix jacchus, are becoming valuable, high-demand biomedical models but their supply outside of their native Brazil is limited. Thus, uncovering the ancestry and biogeographic origins of captive marmosets is an essential task to optimize their biomedical use. Such information facilitates assessing standing levels of captive marmoset genetic diversity, identifying hybrids, and maintaining genetic diversity of captive marmoset populations. Here, we present 16 newly sequenced mitochondrial genomes from marmosets with C. jacchus phenotypes from international biomedical and animal supply facilities. We combine these new data with publicly available sequences to compare usage of mitogenome and the mtDNA control region for determining the ancestry and diversity of captive marmosets. MtDNA control region and mitogenome haplotypes from all 16 newly sequenced marmosets grouped within C. jacchus phylogenetic clades. Overall, phylogenies based on mitogenomes rather than the control region showed better resolved branching patterns with stronger statistical support within C. jacchus clades. However, both mtDNA markers show that C. jacchus is the sister clade to the phylogenetic C. penicillata Caatinga biome clade. Callithrix jacchus mtDNA control region haplotypes showed the highest haplotype and nucleotide diversity, and were intermediate in terms of per sequence theta and polymorphic sites. The larger number of available Callithrix mtDNA control region sequences relative to mitogenomes was used to determine ancestral mtDNA haplotype biogeographic origins. However, biogeographic origins for many C. jacchus haplotypes including those from international captive facilities could not be determined, which we attributed to unknown provenance information from natural C. jacchus populations. Given our results, we encourage increased sampling and genomic sequencing efforts for natural C. jacchus populations in Brazil, and utilizing the mitogenome over other mtDNA markers for phylogenetic-based analyses. Further, combining marmoset mitogenomic data with phenotypic and nuclear genomic data will greatly enhance the biomedical use of marmosets.

Humans have unique cognitive abilities among primates, including language, but their molecular, cellular, and circuit substrates are poorly understood. We used comparative single nucleus transcriptomics in adult humans, chimpanzees, gorillas, rhesus macaques, and common marmosets from the middle temporal gyrus (MTG) to understand human-specific features of cellular and molecular organization. Human, chimpanzee, and gorilla MTG showed highly similar cell type composition and laminar organization, and a large shift in proportions of deep layer intratelencephalic-projecting neurons compared to macaque and marmoset. Species differences in gene expression generally mirrored evolutionary distance and were seen in all cell types, although chimpanzees were more similar to gorillas than humans, consistent with faster divergence along the human lineage. Microglia, astrocytes, and oligodendrocytes showed accelerated gene expression changes compared to neurons or oligodendrocyte precursor cells, indicating either relaxed evolutionary constraints or positive selection in these cell types. Only a few hundred genes showed human-specific patterning in all or specific cell types, and were significantly enriched near human accelerated regions (HARs) and conserved deletions (hCONDELS) and in cell adhesion and intercellular signaling pathways. These results suggest that relatively few cellular and molecular changes uniquely define adult human cortical structure, particularly by affecting circuit connectivity and glial cell function.

Astrocytes play essential roles in normal brain function, with dysfunction implicated in diverse developmental and degenerative disease processes. Emerging evidence of profound species divergent features of astrocytes coupled with the relative inaccessibility of human brain tissue underscore the utility of human pluripotent stem cell (hPSC) technologies for the generation and study of human astrocytes. However, existing approaches for hPSC-astrocyte generation are typically lengthy, incompletely characterized, or require intermediate purification steps, limiting their utility for multi-cell line, adequately powered functional studies. Here, we establish a rapid and highly scalable method for generating functional human induced astrocytes (hiAs) based upon transient Neurogenin 2 (NGN2) induction of neural progenitor-like cells followed by maturation in astrocyte media, which demonstrate remarkable homogeneity within the population and across 11 independent cell lines in the absence of additional purification steps. These hiAs express canonical astrocyte markers, respond to pro-inflammatory stimuli, exhibit ATP-induced calcium transients and support neuronal maturation in vitro. Moreover, single-cell transcriptomic analyses reveal the generation of highly reproducible cell populations across individual donors, most closely resembling human fetal astrocytes, and highly similar to hPSC-derived astrocytes generated using more complex approaches. Finally, the hiAs capture key molecular hallmarks in a trisomy 21 disease model. Thus, hiAs provide a valuable and practical resource well-suited for study of basic human astrocyte function and dysfunction in disease.

Age is the dominant risk factor for infectious diseases, but the mechanisms linking the two are incompletely understood1,2. Age-related mosaic chromosomal alterations (mCAs) detected from blood-derived DNA genotyping, are structural somatic variants associated with aberrant leukocyte cell counts, hematological malignancy, and mortality3-11. Whether mCAs represent independent risk factors for infection is unknown. Here we use genome-wide genotyping of blood DNA to show that mCAs predispose to diverse infectious diseases. We analyzed mCAs from 767,891 individuals without hematological cancer at DNA acquisition across four countries. Expanded mCA (cell fraction10%) prevalence approached 4% by 60 years of age and was associated with diverse incident infections, including sepsis, pneumonia, and coronavirus disease 2019 (COVID-19) hospitalization. A genome-wide association study of expanded mCAs identified 63 significant loci. Germline genetic alleles associated with expanded mCAs were enriched at transcriptional regulatory sites for immune cells. Our results link mCAs with impaired immunity and predisposition to infections. Furthermore, these findings may also have important implications for the ongoing COVID-19 pandemic, particularly in prioritizing individual preventive strategies and evaluating immunization responses.

Rare coding variation has historically provided the most direct connections between gene function and disease pathogenesis. By meta-analysing the whole exomes of 24,248 schizophrenia cases and 97,322 controls, we implicate ultra-rare coding variants (URVs) in 10 genes as conferring substantial risk for schizophrenia (odds ratios of 3–50, P < 2.14 × 10−6) and 32 genes at a false discovery rate of 1, epilepsy and severe neurodevelopmental disorders2, although different mutation types are implicated in some shared genes. Most genes described here, however, are not implicated in neurodevelopment. We demonstrate that genes prioritized from common variant analyses of schizophrenia are enriched in rare variant risk3, suggesting that common and rare genetic risk factors converge at least partially on the same underlying pathogenic biological processes. Even after excluding significantly associated genes, schizophrenia cases still carry a substantial excess of URVs, which indicates that more risk genes await discovery using this approach.

Despite their widespread use in research, there has not yet been a systematic genomic analysis of human embryonic stem cell (hESC) lines at a single-nucleotide resolution. We therefore performed whole-genome sequencing (WGS) of 143 hESC lines and annotated their single-nucleotide and structural genetic variants. We found that while a substantial fraction of hESC lines contained large deleterious structural variants, finer-scale structural and single-nucleotide variants (SNVs) that are ascertainable only through WGS analyses were present in hESC genomes and human blood-derived genomes at similar frequencies. Moreover, WGS allowed us to identify SNVs associated with cancer and other diseases that could alter cellular phenotypes and compromise the safety of hESC-derived cellular products transplanted into humans. As a resource to enable reproducible hESC research and safer translation, we provide a user-friendly WGS data portal and a data-driven scheme for cell line maintenance and selection.

Contains fulltext : 248656.pdf (Publisher’s version ) (Closed access) BACKGROUND: Sex differences in incidence and/or presentation of schizophrenia (SCZ), major depressive disorder (MDD), and bipolar disorder (BIP) are pervasive. Previous evidence for shared genetic risk and sex differences in brain abnormalities across disorders suggest possible shared sex-dependent genetic risk. METHODS: We conducted the largest to date genome-wide genotype-by-sex (G×S) interaction of risk for these disorders using 85,735 cases (33,403 SCZ, 19,924 BIP, and 32,408 MDD) and 109,946 controls from the PGC (Psychiatric Genomics Consortium) and iPSYCH. RESULTS: Across disorders, genome-wide significant single nucleotide polymorphism-by-sex interaction was detected for a locus encompassing NKAIN2 (rs117780815, p = 3.2 × 10(-8)), which interacts with sodium/potassium-transporting ATPase (adenosine triphosphatase) enzymes, implicating neuronal excitability. Three additional loci showed evidence (p < 1 × 10(-6)) for cross-disorder G×S interaction (rs7302529, p = 1.6 × 10(-7); rs73033497, p = 8.8 × 10(-7); rs7914279, p = 6.4 × 10(-7)), implicating various functions. Gene-based analyses identified G×S interaction across disorders (p = 8.97 × 10(-7)) with transcriptional inhibitor SLTM. Most significant in SCZ was a MOCOS gene locus (rs11665282, p = 1.5 × 10(-7)), implicating vascular endothelial cells. Secondary analysis of the PGC-SCZ dataset detected an interaction (rs13265509, p = 1.1 × 10(-7)) in a locus containing IDO2, a kynurenine pathway enzyme with immunoregulatory functions implicated in SCZ, BIP, and MDD. Pathway enrichment analysis detected significant G×S interaction of genes regulating vascular endothelial growth factor receptor signaling in MDD (false discovery rate-corrected p < .05). CONCLUSIONS: In the largest genome-wide G×S analysis of mood and psychotic disorders to date, there was substantial genetic overlap between the sexes. However, significant sex-dependent effects were enriched for genes related to neuronal development and immune and vascular functions across and within SCZ, BIP, and MDD at the variant, gene, and pathway levels.

Primates and rodents, which descended from a common ancestor around 90 million years ago1, exhibit profound differences in behaviour and cognitive capacity; the cellular basis for these differences is unknown. Here we use single-nucleus RNA sequencing to profile RNA expression in 188,776 individual interneurons across homologous brain regions from three primates (human, macaque and marmoset), a rodent (mouse) and a weasel (ferret). Homologous interneuron types—which were readily identified by their RNA-expression patterns—varied in abundance and RNA expression among ferrets, mice and primates, but varied less among primates. Only a modest fraction of the genes identified as ‘markers’ of specific interneuron subtypes in any one species had this property in another species. In the primate neocortex, dozens of genes showed spatial expression gradients among interneurons of the same type, which suggests that regional variation in cortical contexts shapes the RNA expression patterns of adult neocortical interneurons. We found that an interneuron type that was previously associated with the mouse hippocampus—the ‘ivy cell’, which has neurogliaform characteristics—has become abundant across the neocortex of humans, macaques and marmosets but not mice or ferrets. We also found a notable subcortical innovation: an abundant striatal interneuron type in primates that had no molecularly homologous counterpart in mice or ferrets. These interneurons expressed a unique combination of genes that encode transcription factors, receptors and neuropeptides and constituted around 30% of striatal interneurons in marmosets and humans. Single-nucleus RNA-sequencing analyses of brain from humans, macaques, marmosets, mice and ferrets reveal diverse ways that interneuron populations have changed during evolution.

Clonally expanded blood cells with somatic mutations (clonal hematopoiesis, CH) are commonly acquired with age and increase risk of blood cancer1–9. The blood clones identified to date contain diverse large-scale mosaic chromosomal alterations (mCAs: deletions, duplications, and copy-neutral loss of heterozygosity [CN-LOH]) on all chromosomes1,2,5,6,9, but the sources of selective advantage that drive expansion of most clones remain unknown. To identify genes, mutations and biological processes that give selective advantage to mutant clones, we analyzed genotyping data from the blood-derived DNA of 482,789 UK Biobank participants10, identifying 19,632 autosomal mCAs which we analyzed for relationships to inherited genetic variation. Fifty-two inherited, rare, large-effect coding or splice variants in seven genes associated with greatly increased (odds ratios 11 to 758) vulnerability to CH with specific acquired CN-LOH mutations. Acquired mutations systematically replaced the inherited risk alleles (at MPL) or duplicated them to the homologous chromosome (at FH, NBN, MRE11, ATM, SH2B3, and TM2D3). Three of the seven genes (MRE11, NBN, and ATM) encode components of the MRN-ATM pathway, which limits cell division after DNA damage and telomere attrition11−13; another two (MPL, SH2B3) encode proteins that regulate stem cell self-renewal14–16. In addition to these monogenic inherited forms of CH, we found a common and surprisingly polygenic form: CN-LOH mutations across the genome tended to cause chromosomal segments with alleles that promote hematopoietic cell expansions to replace their homologous (allelic) counterparts, increasing polygenic drive for blood-cell proliferation traits. Readily-acquired mutations that replace chromosomal segments with their homologous counterparts appear to interact with pervasive inherited variation to create a challenge for lifelong cytopoiesis.

Meiosis, although essential for reproduction, is also variable and error-prone: rates of chromosome crossover vary among gametes, between the sexes, and among humans of the same sex, and chromosome missegregation leads to abnormal chromosome numbers (aneuploidy)1–8. To study diverse meiotic outcomes and how they covary across chromosomes, gametes and humans, we developed Sperm-seq, a way of simultaneously analysing the genomes of thousands of individual sperm. Here we analyse the genomes of 31,228 human gametes from 20 sperm donors, identifying 813,122 crossovers and 787 aneuploid chromosomes. Sperm donors had aneuploidy rates ranging from 0.01 to 0.05 aneuploidies per gamete; crossovers partially protected chromosomes from nondisjunction at the meiosis I cell division. Some chromosomes and donors underwent more-frequent nondisjunction during meiosis I, and others showed more meiosis II segregation failures. Sperm genomes also manifested many genomic anomalies that could not be explained by simple nondisjunction. Diverse recombination phenotypes—from crossover rates to crossover location and separation, a measure of crossover interference—covaried strongly across individuals and cells. Our results can be incorporated with earlier observations into a unified model in which a core mechanism, the variable physical compaction of meiotic chromosomes, generates interindividual and cell-to-cell variation in diverse meiotic phenotypes. Thousands of sperm genomes have been analysed with a new method called Sperm-seq, revealing interconnected meiotic variation at the single-cell and person-to-person levels, and suggesting chromosome compaction as a way to explain the relationships between diverse recombination phenotypes.

Protein-coding de novo mutations (DNMs) are significant risk factors in many neurodevelopmental disorders, whereas schizophrenia (SCZ) risk associated with DNMs has thus far been shown to be modest. We analyzed DNMs from 1,695 SCZ-affected trios and 1,077 published SCZ-affected trios to better understand the contribution to SCZ risk. Among 2,772 SCZ probands, exome-wide DNM burden remained modest. Gene set analyses revealed that SCZ DNMs were significantly concentrated in genes that were highly expressed in the brain, that were under strong evolutionary constraint and/or overlapped with genes identified in other neurodevelopmental disorders. No single gene surpassed exome-wide significance; however, 16 genes were recurrently hit by protein-truncating DNMs, corresponding to a 3.15-fold higher rate than the mutation model expectation (permuted 95% confidence interval: 1-10 genes; permuted P = 3 × 10-5). Overall, DNMs explain a small fraction of SCZ risk, and larger samples are needed to identify individual risk genes, as coding variation across many genes confers risk for SCZ in the population.

The four voltage-gated sodium channels SCN1/2/3/8A have been associated with heterogeneous types of developmental disorders, each presenting with disease specific temporal and cell type specific gene expression. Using single-cell RNA sequencing transcriptomic data from humans and mice, we observe that SCN1A is predominantly expressed in inhibitory neurons. In contrast, SCN2/3/8A are profoundly expressed in excitatory neurons with SCN2/3A starting prenatally, followed by SCN1/8A neonatally. In contrast to previous observations from low resolution RNA screens, we observe that all four genes are expressed in both excitatory and inhibitory neurons, however, exhibit differential expression strength. These findings provide molecular evidence, at single-cell resolution, to support the hypothesis that the excitatory/inhibitory (E/I) neuronal expression ratios of sodium channels are important regulatory mechanisms underlying brain homeostasis and neurological diseases. Modulating the E/I expression balance within cell types of sodium channels could serve as a potential strategy to develop targeted treatment for NaV-associated neuronal developmental disorders.

Mosaic loss of chromosome Y (LOY) in circulating white blood cells is the most common form of clonal mosaicism1–5, yet our knowledge of the causes and consequences of this is limited. Here, using a computational approach, we estimate that 20% of the male population represented in the UK Biobank study (n = 205,011) has detectable LOY. We identify 156 autosomal genetic determinants of LOY, which we replicate in 757,114 men of European and Japanese ancestry. These loci highlight genes that are involved in cell-cycle regulation and cancer susceptibility, as well as somatic drivers of tumour growth and targets of cancer therapy. We demonstrate that genetic susceptibility to LOY is associated with non-haematological effects on health in both men and women, which supports the hypothesis that clonal haematopoiesis is a biomarker of genomic instability in other tissues. Single-cell RNA sequencing identifies dysregulated expression of autosomal genes in leukocytes with LOY and provides insights into why clonal expansion of these cells may occur. Collectively, these data highlight the value of studying clonal mosaicism to uncover fundamental mechanisms that underlie cancer and other ageing-related diseases. A genome-wide association study of mosaic loss of chromosome Y (LOY) in UK Biobank participants identifies 156 genetic determinants of LOY, showing that LOY is associated with cancer and non-haematological health outcomes.

Schizophrenia is a common, chronic and debilitating neuropsychiatric syndrome affecting tens of millions of individuals worldwide. While rare genetic variants play a role in the etiology of schizophrenia, most of the currently explained liability is within common variation, suggesting that variation predating the human diaspora out of Africa harbors a large fraction of the common variant attributable heritability. However, common variant association studies in schizophrenia have concentrated mainly on cohorts of European descent. We describe genome-wide association studies of 6152 cases and 3918 controls of admixed African ancestry, and of 1234 cases and 3090 controls of Latino ancestry, representing the largest such study in these populations to date. Combining results from the samples with African ancestry with summary statistics from the Psychiatric Genomics Consortium (PGC) study of schizophrenia yielded seven newly genome-wide significant loci, and we identified an additional eight loci by incorporating the results from samples with Latino ancestry. Leveraging population differences in patterns of linkage disequilibrium, we achieve improved fine-mapping resolution at 22 previously reported and 4 newly significant loci. Polygenic risk score profiling revealed improved prediction based on trans-ancestry meta-analysis results for admixed African (Nagelkerke’s R2 = 0.032; liability R2 = 0.017; P −52), Latino (Nagelkerke’s R2 = 0.089; liability R2 = 0.021; P −58), and European individuals (Nagelkerke’s R2 = 0.089; liability R2 = 0.037; P −113), further highlighting the advantages of incorporating data from diverse human populations.

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

The primary motor cortex (M1) is essential for voluntary fine-motor control and is functionally conserved across mammals1. Here, using high-throughput transcriptomic and epigenomic profiling of more than 450,000 single nuclei in humans, marmoset monkeys and mice, we demonstrate a broadly conserved cellular makeup of this region, with similarities that mirror evolutionary distance and are consistent between the transcriptome and epigenome. The core conserved molecular identities of neuronal and non-neuronal cell types allow us to generate a cross-species consensus classification of cell types, and to infer conserved properties of cell types across species. Despite the overall conservation, however, many species-dependent specializations are apparent, including differences in cell-type proportions, gene expression, DNA methylation and chromatin state. Few cell-type marker genes are conserved across species, revealing a short list of candidate genes and regulatory mechanisms that are responsible for conserved features of homologous cell types, such as the GABAergic chandelier cells. This consensus transcriptomic classification allows us to use patch–seq (a combination of whole-cell patch-clamp recordings, RNA sequencing and morphological characterization) to identify corticospinal Betz cells from layer 5 in non-human primates and humans, and to characterize their highly specialized physiology and anatomy. These findings highlight the robust molecular underpinnings of cell-type diversity in M1 across mammals, and point to the genes and regulatory pathways responsible for the functional identity of cell types and their species-specific adaptations., An examination of motor cortex in humans, marmosets and mice reveals a generally conserved cellular makeup that is likely to extend to many mammalian species, but also differences in gene expression, DNA methylation and chromatin state that lead to species-dependent specializations.

To study how the 22q11.2 deletion predisposes to psychiatric disease, we generated induced pluripotent stem cells from deletion carriers and controls, as well as utilized CRISPR/Cas9 to introduce the heterozygous deletion into a control cell line. Upon differentiation into neural progenitor cells, we found the deletion acted in trans to alter the abundance of transcripts associated with risk for neurodevelopmental disorders including Autism Spectrum Disorder. In more differentiated excitatory neurons, altered transcripts encoded presynaptic factors and were associated with genetic risk for schizophrenia, including common (per-SNP heritability p (τc)= 4.2 x 10-6) and rare, loss of function variants (p = 1.29×10-12). These findings suggest a potential relationship between cellular states, developmental windows and susceptibility to psychiatric conditions with different ages of onset. To understand how the deletion contributed to these observed changes in gene expression, we developed and applied PPItools, which identifies the minimal protein-protein interaction network that best explains an observed set of gene expression alterations. We found that many of the genes in the 22q11.2 interval interact in presynaptic, proteasome, and JUN/FOS transcriptional pathways that underlie the broader alterations in psychiatric risk gene expression we identified. Our findings suggest that the 22q11.2 deletion impacts genes and pathways that may converge with risk loci implicated by psychiatric genetic studies to influence disease manifestation in each deletion carrier.

Recent genetic discoveries in schizophrenia have highlighted neuronal genes with functions in the synapse. Although emblematic of neurons, the development of synapses and neuronal maturation relies on interactions with glial cells including astrocytes. To study the role of glia-neuron interactions in schizophrenia, we generated RNA sequence data from human pluripotent stem cell (hPSC) derived neurons that were cocultured with glial cells. We found that expression of genes characteristic of astrocytes induced the expression of post-synaptic genetic programs in neurons, consistent with advanced neuronal maturation. We further found that the astrocyte-induced genes in neurons were associated with risk for schizophrenia. To understand how glial cells promoted neuronal maturation, we studied the association of transcript abundances in glial cells to gene expression in neurons. We found that expression of synaptic cell adhesion molecules in glial cells corresponded to induced synaptic transcripts in neurons and were associated with genetic risk for schizophrenia. These included 11 genes in significant GWAS loci and three with direct genetic evidence for the disorder (MAGI2, NRXN1, LRRC4B, and MSI2). Our results suggest that astrocyte-expressed genes with functions in the synapse are associated with schizophrenia and promote synaptic genetic programs in neurons, and further highlight the potential role for astrocyte-neuron interactions in schizophrenia.

Brain function depends on forming and maintaining connections between neurons of specific types, ensuring neural function while allowing the plasticity necessary for cellular and behavioral dynamics. However, systematic descriptions of how brain cell types organize into synaptic networks and which molecules instruct these relationships are not readily available. Here, we introduce SBARRO (Synaptic Barcode Analysis by Retrograde Rabies ReadOut), a method that uses single-cell RNA sequencing to reveal directional, monosynaptic relationships based on the paths of a barcoded rabies virus from its “starter” postsynaptic cell to that cell’s presynaptic partners1. Thousands of these partner relationships can be ascertained in a single experiment, alongside genome-wide RNA profiles – and thus cell identities and molecular states – of each host cell. We used SBARRO to describe synaptic networks formed by diverse mouse brain cell types in vitro, leveraging a system similar to those used to identify synaptogenic molecules. We found that the molecular identity (cell type/subtype) of the starter cell predicted the number and types of cells that had synapsed onto it. Rabies transmission tended to occur into cells with RNA-expression signatures related to developmental maturation and synaptic transmission. The estimated size of a cell’s presynaptic network, relative to that of other cells of the same type, associated with increased expression of Arpp21 and Cdh13. By tracking individual virions and their clonal progeny as they travel among host cells, single-cell, single-virion genomic technologies offer new opportunities to map the synaptic organization of neural circuits in health and disease.

Bipolar disorder (BD) is a serious mental illness with substantial common variant heritability. However, the role of rare coding variation in BD is not well established. We examined the protein-coding (exonic) sequences of 3,987 unrelated individuals with BD and 5,322 controls of predominantly European ancestry across four cohorts from the Bipolar Sequencing Consortium (BSC). We assessed the burden of rare, protein-altering, single nucleotide variants classified as pathogenic or likely pathogenic (P-LP) both exome-wide and within several groups of genes with phenotypic or biologic plausibility in BD. While we observed an increased burden of rare coding P-LP variants within 165 genes identified as BD GWAS regions in 3,987 BD cases (meta-analysis OR = 1.9, 95% CI = 1.3–2.8, one-sided p = 6.0 × 10−4), this enrichment did not replicate in an additional 9,929 BD cases and 14,018 controls (OR = 0.9, one-side p = 0.70). Although BD shares common variant heritability with schizophrenia, in the BSC sample we did not observe a significant enrichment of P-LP variants in SCZ GWAS genes, in two classes of neuronal synaptic genes (RBFOX2 and FMRP) associated with SCZ or in loss-of-function intolerant genes. In this study, the largest analysis of exonic variation in BD, individuals with BD do not carry a replicable enrichment of rare P-LP variants across the exome or in any of several groups of genes with biologic plausibility. Moreover, despite a strong shared susceptibility between BD and SCZ through common genetic variation, we do not observe an association between BD risk and rare P-LP coding variants in genes known to modulate risk for SCZ., International Bipolar Sequencing Consortium (NIMH) R01 MH 110437, Whole Genome Sequencing for Schizophrenia and Biopolar Disorder in the GPS (NIMH) UO1 MH105653, Whole Genome and Exome Sequencing for Bipolar Disorder (NIMH) R01 MH 094145, Multi-ethnic GWAS of Bipolar I Disorder (NIH) R01 MH 085543, Genetic Epidemiology Research in Adult Health and Aging (GERA) RC2 AG036607, Kaiser Permanente Research Program on Genes, Environment, and Health, Rare Bipolar Loci Identification Through Synaptome Sequencing (NIMH) R01 MH 087979 MH 087992, 2/2 Large Scale Genetic Studies of Schizophrenia in Sweden R01MH095034, 1/3 Genetic Analysis of the International Cohort Collection for Bipolar Disorder R01MH106531, 1/2 Large Scale Genetic Studies of Schizophrenia in Sweden R01MH077139, 2/3 Genetic Analysis of the International Cohort Collection for Bipolar Disorder R01MH106527, Dalio Foundation

Polygenic risk scores (PRS) summarize genetic liability to a disease at the individual level, and the aim is to use them as biomarkers of disease and poor outcomes in real-world clinical practice. To date, few studies have assessed the prognostic value of PRS relative to standards of care. Schizophrenia (SCZ), the archetypal psychotic illness, is an ideal test case for this because the predictive power of the SCZ PRS exceeds that of most other common diseases. Here, we analyzed clinical and genetic data from two multi-ethnic cohorts totaling 8,541 adults with SCZ and related psychotic disorders, to assess whether the SCZ PRS improves the prediction of poor outcomes relative to clinical features captured in a standard psychiatric interview. For all outcomes investigated, the SCZ PRS did not improve the performance of predictive models, an observation that was generally robust to divergent case ascertainment strategies and the ancestral background of the study participants. The inclusion of polygenic risk scores does not improve the performance of standard-of-care predictive models of disease outcomes in patients with psychosis.

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.

Significance By evaluating children with a malformed cerebral cortex, we identified an ATPase pump (ATP1A3) with an early role in brain development. The ATP1A3 pump maintains the physiological concentration of sodium and potassium ions in cells, a process critical for osmotic equilibrium and membrane potential across several developing cell populations. We employed single-cell sequencing approaches to identify key enrichments for ATP1A3 expression during human cortex development. Unravelling this early cell-type–specific pathophysiology in the developing brain offers a potential basis for the treatment of ATP1A3-related diseases affecting prenatal and early childhood development., Osmotic equilibrium and membrane potential in animal cells depend on concentration gradients of sodium (Na+) and potassium (K+) ions across the plasma membrane, a function catalyzed by the Na+,K+-ATPase α-subunit. Here, we describe ATP1A3 variants encoding dysfunctional α3-subunits in children affected by polymicrogyria, a developmental malformation of the cerebral cortex characterized by abnormal folding and laminar organization. To gain cell-biological insights into the spatiotemporal dynamics of prenatal ATP1A3 expression, we built an ATP1A3 transcriptional atlas of fetal cortical development using mRNA in situ hybridization and transcriptomic profiling of ∼125,000 individual cells with single-cell RNA sequencing (Drop-seq) from 11 areas of the midgestational human neocortex. We found that fetal expression of ATP1A3 is most abundant to a subset of excitatory neurons carrying transcriptional signatures of the developing subplate, yet also maintains expression in nonneuronal cell populations. Moving forward a year in human development, we profiled ∼52,000 nuclei from four areas of an infant neocortex and show that ATP1A3 expression persists throughout early postnatal development, most predominantly in inhibitory neurons, including parvalbumin interneurons in the frontal cortex. Finally, we discovered the heteromeric Na+,K+-ATPase pump complex may form nonredundant cell-type–specific α-β isoform combinations, including α3-β1 in excitatory neurons and α3-β2 in inhibitory neurons. Together, the developmental malformation phenotype of affected individuals and single-cell ATP1A3 expression patterns point to a key role for α3 in human cortex development, as well as a cell-type basis for pre- and postnatal ATP1A3-associated diseases.

Osmotic equilibrium and membrane potential in animal cells depend on concentration gradients of sodium (Na+) and potassium (K+) ions across the plasma membrane, a function catalyzed by the Na+,K+-ATPase α-subunit. Here, we describe ATP1A3 variants encoding dysfunctional α3-subunits in children affected by polymicrogyria, a developmental malformation of the cerebral cortex characterized by abnormal folding and laminar organization. To gain cell-biological insights into the spatiotemporal dynamics of prenatal ATP1A3 expression, we built an ATP1A3 transcriptional atlas of fetal cortical development using mRNA in situ hybridization and transcriptomic profiling of ∼125,000 individual cells with single-cell RNA sequencing (Drop-seq) from 11 areas of the midgestational human neocortex. We found that fetal expression of ATP1A3 is most abundant to a subset of excitatory neurons carrying transcriptional signatures of the developing subplate, yet also maintains expression in nonneuronal cell populations. Moving forward a year in human development, we profiled ∼52,000 nuclei from four areas of an infant neocortex and show that ATP1A3 expression persists throughout early postnatal development, most predominantly in inhibitory neurons, including parvalbumin interneurons in the frontal cortex. Finally, we discovered the heteromeric Na+,K+-ATPase pump complex may form nonredundant cell-type-specific α-β isoform combinations, including α3-β1 in excitatory neurons and α3-β2 in inhibitory neurons. Together, the developmental malformation phenotype of affected individuals and single-cell ATP1A3 expression patterns point to a key role for α3 in human cortex development, as well as a cell-type basis for pre- and postnatal ATP1A3-associated diseases.

PURPOSE: This pilot study sought to (1) validate the use of a novel technology for single-sperm-cell genome sequencing (Sperm-seq) in infertile men who may not have optimal quantity or quality of sperm for genomic analysis and (2) compare these results to fertile donors. METHODS: Infertile men undergoing IVF with female partners with a previous history of failed fertilization with ICSI (FF) or poor blastulation of embryos (PB) were recruited from a large IVF center. Sperm-seq was used to analyze thousands of individual sperm and was carried out at an affiliated university research institute. Global aneuploidy rate, crossover locations, and crossover frequencies were assessed in the infertile population, and compared with a control group of 20 fertile donors, which were analyzed previously at the same laboratory. RESULTS: Eight patients were initially included, but 3 samples did not yield high-quality genomic data for analysis. A total of 10,042 sperm were analyzed from 5 patients, 2 in the FF group, and 3 in the PB group. The global aneuploidy rate among the samples was 2–4%, similar to the control group. Likewise, crossover locations and frequencies were similar. CONCLUSION: Sperm-seq provides a robust analysis but may not be applicable to all male infertility cases due to technical limitations. This group of male infertility patients did not have higher rates of aneuploidy or abnormal crossover patterns compared to a fertile donor population. Our data may suggest that FF and PB phenotypes may not be related to sperm aneuploidy or meiotic errors but rather to other intrinsic nuclear anomalies. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s10815-021-02300-3.

Schizophrenia, Schizoaffective, and Bipolar disorders share behavioral and phenomenological traits, intermediate phenotypes, and some associated genetic loci with pleiotropic effects. Volumetric abnormalities in brain structures are among the intermediate phenotypes consistently reported associated with these disorders. In order to examine the genetic underpinnings of these structural brain modifications, we performed genome-wide association analyses (GWAS) on 60 quantitative structural brain MRI phenotypes in a sample of 777 subjects (483 cases and 294 controls pooled together). Genotyping was performed with the Illumina PsychChip microarray, followed by imputation to the 1000 genomes multiethnic reference panel. Enlargement of the Temporal Horns of Lateral Ventricles (THLV) is associated with an intronic SNP of the gene NRXN1 (rs12467877, P = 6.76E–10), which accounts for 4.5% of the variance in size. Enlarged THLV is associated with psychosis in this sample, and with reduction of the hippocampus and enlargement of the choroid plexus and caudate. Eight other suggestively significant associations (P

Hundreds of the proteins encoded in human genomes contain domains that vary in size or copy number due to variable numbers of tandem repeats (VNTRs) in protein-coding exons. VNTRs have eluded analysis by the molecular methods-SNP arrays and high-throughput sequencing-used in large-scale human genetic studies to date; thus, the relationships of VNTRs to most human phenotypes are unknown. We developed ways to estimate VNTR lengths from whole-exome sequencing data, identify the SNP haplotypes on which VNTR alleles reside, and use imputation to project these haplotypes into abundant SNP data. We analyzed 118 protein-altering VNTRs in 415,280 UK Biobank participants for association with 791 phenotypes. Analysis revealed some of the strongest associations of common variants with human phenotypes including height, hair morphology, and biomarkers of human health; for example, a VNTR encoding 13-44 copies of a 19-amino-acid repeat in the chondroitin sulfate domain of aggrecan (ACAN) associated with height variation of 3.4 centimeters (s.e. 0.3 cm). Incorporating large-effect VNTRs into analysis also made it possible to map many additional effects at the same loci: for the blood biomarker lipoprotein(a), for example, analysis of the kringle IV-2 VNTR within the LPA gene revealed that 18 coding SNPs and the VNTR in LPA explained 90% of lipoprotein(a) heritability in Europeans, enabling insights about population differences and epidemiological significance of this clinical biomarker. These results point to strong, cryptic effects of highly polymorphic common structural variants that have largely eluded molecular analyses to date.

BACKGROUND: Polygenic scores (PGSs), which assess the genetic risk of individuals for a disease, are calculated as a weighted count of risk alleles identified in genome-wide association studies. PGS methods differ in which DNA variants are included and the weights assigned to them; some require an independent tuning sample to help inform these choices. PGSs are evaluated in independent target cohorts with known disease status. Variability between target cohorts is observed in applications to real data sets, which could reflect a number of factors, e.g., phenotype definition or technical factors.METHODS: The Psychiatric Genomics Consortium Working Groups for schizophrenia and major depressive disorder bring together many independently collected case-control cohorts. We used these resources (31,328 schizophrenia cases, 41,191 controls; 248,750 major depressive disorder cases, 563,184 controls) in repeated application of leave-one-cohort-out meta-analyses, each used to calculate and evaluate PGS in the left-out (target) cohort. Ten PGS methods (the baseline PC+T method and 9 methods that model genetic architecture more formally: SBLUP, LDpred2-Inf, LDpred-funct, LDpred2, Lassosum, PRS-CS, PRS-CS-auto, SBayesR, MegaPRS) were compared.RESULTS: Compared with PC+T, the other 9 methods gave higher prediction statistics, MegaPRS, LDPred2, and SBayesR significantly so, explaining up to 9.2% variance in liability for schizophrenia across 30 target cohorts, an increase of 44%. For major depressive disorder across 26 target cohorts, these statistics were 3.5% and 59%, respectively.CONCLUSIONS: Although the methods that more formally model genetic architecture have similar performance, MegaPRS, LDpred2, and SBayesR rank highest in most comparisons and are recommended in applications to psychiatric disorders.