Rare variation in noncoding regions with evolutionary signatures contributes to autism spectrum disorder risk.

Little is known about the role of noncoding regions in the etiology of autism spectrum disorder (ASD). We examined three classes of noncoding regions: Human Accelerated Regions (HARs), which show signatures of positive selection in humans; experimentally validated neural Vista Enhancers (VEs); and conserved regions predicted to act as neural enhancers (CNEs). Targeted and whole genome analysis of >16,600 samples and >4900 ASD probands revealed that likely recessive, rare, inherited variants in HARs, VEs, and CNEs substantially contribute to ASD risk in probands whose parents share ancestry, which enriches for recessive contributions, but modestly, if at all, in simplex family structures. We identified multiple patient variants in HARs near IL1RAPL1 and in a VE near SIM1 and showed that they change enhancer activity. Our results implicate both human-evolved and evolutionarily conserved noncoding regions in ASD risk and suggest potential mechanisms of how changes in regulatory regions can modulate social behavior.

Comparative transcriptomics reveals human-specific cortical features.

The cognitive abilities of humans are distinctive among primates, but their molecular and cellular substrates are poorly understood. We used comparative single-nucleus transcriptomics to analyze samples of the middle temporal gyrus (MTG) from adult humans, chimpanzees, gorillas, rhesus macaques, and common marmosets to understand human-specific features of the neocortex. Human, chimpanzee, and gorilla MTG showed highly similar cell-type composition and laminar organization as well as a large shift in proportions of deep-layer intratelencephalic-projecting neurons compared with macaque and marmoset MTG. Microglia, astrocytes, and oligodendrocytes had more-divergent expression across species compared with neurons or oligodendrocyte precursor cells, and neuronal expression diverged more rapidly on the human lineage. Only a few hundred genes showed human-specific patterning, suggesting that relatively few cellular and molecular changes distinctively define adult human cortical structure.

Exome Sequencing and the Identification of New Genes and Shared Mechanisms in Polymicrogyria.

Importance: Polymicrogyria is the most commonly diagnosed cortical malformation and is associated with neurodevelopmental sequelae including epilepsy, motor abnormalities, and cognitive deficits. Polymicrogyria frequently co-occurs with other brain malformations or as part of syndromic diseases. Past studies of polymicrogyria have defined heterogeneous genetic and nongenetic causes but have explained only a small fraction of cases. Objective: To survey germline genetic causes of polymicrogyria in a large cohort and to consider novel polymicrogyria gene associations. Design, Setting, and Participants: This genetic association study analyzed panel sequencing and exome sequencing of accrued DNA samples from a retrospective cohort of families with members with polymicrogyria. Samples were accrued over more than 20 years (1994 to 2020), and sequencing occurred in 2 stages: panel sequencing (June 2015 to January 2016) and whole-exome sequencing (September 2019 to March 2020). Individuals seen at multiple clinical sites for neurological complaints found to have polymicrogyria on neuroimaging, then referred to the research team by evaluating clinicians, were included in the study. Targeted next-generation sequencing and/or exome sequencing were performed on probands (and available parents and siblings) from 284 families with individuals who had isolated polymicrogyria or polymicrogyria as part of a clinical syndrome and no genetic diagnosis at time of referral from clinic, with sequencing from 275 families passing quality control. Main Outcomes and Measures: The number of families in whom genetic sequencing yielded a molecular diagnosis that explained the polymicrogyria in the family. Secondarily, the relative frequency of different genetic causes of polymicrogyria and whether specific genetic causes were associated with co-occurring head size changes were also analyzed. Results: In 32.7% (90 of 275) of polymicrogyria-affected families, genetic variants were identified that provided satisfactory molecular explanations. Known genes most frequently implicated by polymicrogyria-associated variants in this cohort were PIK3R2, TUBB2B, COL4A1, and SCN3A. Six candidate novel polymicrogyria genes were identified or confirmed: de novo missense variants in PANX1, QRICH1, and SCN2A and compound heterozygous variants in TMEM161B, KIF26A, and MAN2C1, each with consistent genotype-phenotype relationships in multiple families. Conclusions and Relevance: This study's findings reveal a higher than previously recognized rate of identifiable genetic causes, specifically of channelopathies, in individuals with polymicrogyria and support the utility of exome sequencing for families affected with polymicrogyria.

Cbln1 Directs Axon Targeting by Corticospinal Neurons Specifically toward Thoraco-Lumbar Spinal Cord.

Corticospinal neurons (CSN) are centrally required for skilled voluntary movement, which necessitates that they establish precise subcerebral connectivity with the brainstem and spinal cord. However, molecular controls regulating specificity of this projection targeting remain largely unknown. We previously identified that developing CSN subpopulations exhibit striking axon targeting specificity in the spinal white matter. These CSN subpopulations with segmentally distinct spinal projections are also molecularly distinct; a subset of differentially expressed genes between these distinct CSN subpopulations regulate differential axon projection targeting. Rostrolateral CSN extend axons exclusively to bulbar-cervical segments (CSNBC-lat), while caudomedial CSN (CSNmedial) are more heterogeneous, with distinct, intermingled subpopulations extending axons to either bulbar-cervical or thoraco-lumbar segments. Here, we report, in male and female mice, that Cerebellin 1 (Cbln1) is expressed specifically by CSN in medial, but not lateral, sensorimotor cortex. Cbln1 shows highly dynamic temporal expression, with Cbln1 levels in CSN highest during the period of peak axon extension toward thoraco-lumbar segments. Using gain-of-function experiments, we identify that Cbln1 is sufficient to direct thoraco-lumbar axon extension by CSN. Misexpression of Cbln1 in CSNBC-lat either by in utero electroporation, or by postmitotic AAV-mediated gene delivery, redirects these axons past their normal bulbar-cervical targets toward thoracic segments. Further, Cbln1 overexpression in postmitotic CSNBC-lat increases the number of CSNmedial axons that extend past cervical segments into the thoracic cord. Collectively, these results identify that Cbln1 functions as a potent molecular control over thoraco-lumbar CSN axon extension, part of an integrated network of controls over segmentally-specific CSN axon projection targeting.SIGNIFICANCE STATEMENT Corticospinal neurons (CSN) exhibit remarkable diversity and precision of axonal projections to targets in the brainstem and distinct spinal segments; the molecular basis for this targeting diversity is largely unknown. CSN subpopulations projecting to distinct targets are also molecularly distinguishable. Distinct subpopulations degenerate in specific motor neuron diseases, further suggesting that intrinsic molecular differences might underlie differential vulnerability to disease. Here, we identify a novel molecular control, Cbln1, expressed by CSN extending axons to thoraco-lumbar spinal segments. Cbln1 is sufficient, but not required, for CSN axon extension toward distal spinal segments, and Cbln1 expression is controlled by recently identified, CSN-intrinsic regulators of axon extension. Our results identify that Cbln1, together with other regulators, coordinates segmentally precise CSN axon targeting.

TMEM161B regulates cerebral cortical gyration, Sonic Hedgehog signaling, and ciliary structure in the developing central nervous system.

Sonic hedgehog signaling regulates processes of embryonic development across multiple tissues, yet factors regulating context-specific Shh signaling remain poorly understood. Exome sequencing of families with polymicrogyria (disordered cortical folding) revealed multiple individuals with biallelic deleterious variants in TMEM161B, which encodes a multi-pass transmembrane protein of unknown function. Tmem161b null mice demonstrated holoprosencephaly, craniofacial midline defects, eye defects, and spinal cord patterning changes consistent with impaired Shh signaling, but were without limb defects, suggesting a CNS-specific role of Tmem161b. Tmem161b depletion impaired the response to Smoothened activation in vitro and disrupted cortical histogenesis in vivo in both mouse and ferret models, including leading to abnormal gyration in the ferret model. Tmem161b localizes non-exclusively to the primary cilium, and scanning electron microscopy revealed shortened, dysmorphic, and ballooned ventricular zone cilia in the Tmem161b null mouse, suggesting that the Shh-related phenotypes may reflect ciliary dysfunction. Our data identify TMEM161B as a regulator of cerebral cortical gyration, as involved in primary ciliary structure, as a regulator of Shh signaling, and further implicate Shh signaling in human gyral development.

Genetic studies of human-chimpanzee divergence using stem cell fusions.

Complete genome sequencing has identified millions of DNA changes that differ between humans and chimpanzees. Although a subset of these changes likely underlies important phenotypic differences between humans and chimpanzees, it is currently difficult to distinguish causal from incidental changes and to map specific phenotypes to particular genome locations. To facilitate further genetic study of human-chimpanzee divergence, we have generated human and chimpanzee autotetraploids and allotetraploids by fusing induced pluripotent stem cells (iPSCs) of each species. The resulting tetraploid iPSCs can be stably maintained and retain the ability to differentiate along ectoderm, mesoderm, and endoderm lineages. RNA sequencing identifies thousands of genes whose expression differs between humans and chimpanzees when assessed in single-species diploid or autotetraploid iPSCs. Analysis of gene expression patterns in interspecific allotetraploid iPSCs shows that human-chimpanzee expression differences arise from substantial contributions of both cis-acting changes linked to the genes themselves and trans-acting changes elsewhere in the genome. To enable further genetic mapping of species differences, we tested chemical treatments for stimulating genome-wide mitotic recombination between human and chimpanzee chromosomes, and CRISPR methods for inducing species-specific changes on particular chromosomes in allotetraploid cells. We successfully generated derivative cells with nested deletions or interspecific recombination on the X chromosome. These studies confirm an important role for the X chromosome in trans regulation of expression differences between species and illustrate the potential of this system for more detailed cis and trans mapping of the molecular basis of human and chimpanzee evolution.

Corticospinal neuron subpopulation-specific developmental genes prospectively indicate mature segmentally specific axon projection targeting.

For precise motor control, distinct subpopulations of corticospinal neurons (CSN) must extend axons to distinct spinal segments, from proximal targets in the brainstem and cervical cord to distal targets in thoracic and lumbar spinal segments. We find that developing CSN subpopulations exhibit striking axon targeting specificity in spinal white matter, which establishes the foundation for durable specificity of adult corticospinal circuitry. Employing developmental retrograde and anterograde labeling, and their distinct neocortical locations, we purified developing CSN subpopulations using fluorescence-activated cell sorting to identify genes differentially expressed between bulbar-cervical and thoracolumbar-projecting CSN subpopulations at critical developmental times. These segmentally distinct CSN subpopulations are molecularly distinct from the earliest stages of axon extension, enabling prospective identification even before eventual axon targeting decisions are evident in the spinal cord. This molecular delineation extends beyond simple spatial separation of these subpopulations in the cortex. Together, these results identify candidate molecular controls over segmentally specific corticospinal axon projection targeting.

Rewiring of human neurodevelopmental gene regulatory programs by human accelerated regions.

Human accelerated regions (HARs) are the fastest-evolving regions of the human genome, and many are hypothesized to function as regulatory elements that drive human-specific gene regulatory programs. We interrogate the in vitro enhancer activity and in vivo epigenetic landscape of more than 3,100 HARs during human neurodevelopment, demonstrating that many HARs appear to act as neurodevelopmental enhancers and that sequence divergence at HARs has largely augmented their neuronal enhancer activity. Furthermore, we demonstrate PPP1R17 to be a putative HAR-regulated gene that has undergone remarkable rewiring of its cell type and developmental expression patterns between non-primates and primates and between non-human primates and humans. Finally, we show that PPP1R17 slows neural progenitor cell cycle progression, paralleling the cell cycle length increase seen predominantly in primate and especially human neurodevelopment. Our findings establish HARs as key components in rewiring human-specific neurodevelopmental gene regulatory programs and provide an integrated resource to study enhancer activity of specific HARs.

PRG2 and AQPEP are misexpressed in fetal membranes in placenta previa and percreta†.

The obstetrical conditions placenta accreta spectrum (PAS) and placenta previa are a significant source of pregnancy-associated morbidity and mortality, yet the specific molecular and cellular underpinnings of these conditions are not known. In this study, we identified misregulated gene expression patterns in tissues from placenta previa and percreta (the most extreme form of PAS) compared with control cases. By comparing this gene set with existing placental single-cell and bulk RNA-Seq datasets, we show that the upregulated genes predominantly mark extravillous trophoblasts. We performed immunofluorescence on several candidate molecules and found that PRG2 and AQPEP protein levels are upregulated in both the fetal membranes and the placental disk in both conditions. While this increased AQPEP expression remains restricted to trophoblasts, PRG2 is mislocalized and is found throughout the fetal membranes. Using a larger patient cohort with a diverse set of gestationally aged-matched controls, we validated PRG2 as a marker for both previa and PAS and AQPEP as a marker for only previa in the fetal membranes. Our findings suggest that the extraembryonic tissues surrounding the conceptus, including both the fetal membranes and the placental disk, harbor a signature of previa and PAS that is characteristic of EVTs and that may reflect increased trophoblast invasiveness.

Characterization of a Human-Specific Tandem Repeat Associated with Bipolar Disorder and Schizophrenia.

Bipolar disorder (BD) and schizophrenia (SCZ) are highly heritable diseases that affect more than 3% of individuals worldwide. Genome-wide association studies have strongly and repeatedly linked risk for both of these neuropsychiatric diseases to a 100 kb interval in the third intron of the human calcium channel gene CACNA1C. However, the causative mutation is not yet known. We have identified a human-specific tandem repeat in this region that is composed of 30 bp units, often repeated hundreds of times. This large tandem repeat is unstable using standard polymerase chain reaction and bacterial cloning techniques, which may have resulted in its incorrect size in the human reference genome. The large 30-mer repeat region is polymorphic in both size and sequence in human populations. Particular sequence variants of the 30-mer are associated with risk status at several flanking single-nucleotide polymorphisms in the third intron of CACNA1C that have previously been linked to BD and SCZ. The tandem repeat arrays function as enhancers that increase reporter gene expression in a human neural progenitor cell line. Different human arrays vary in the magnitude of enhancer activity, and the 30-mer arrays associated with increased psychiatric disease risk status have decreased enhancer activity. Changes in the structure and sequence of these arrays likely contribute to changes in CACNA1C function during human evolution and may modulate neuropsychiatric disease risk in modern human populations.