Roychowdhury, Tanmoy

BACKGROUND: Bicuspid aortic valve (BAV) is a frequent congenital heart defect with a high heritability. Despite this, only a limited number of genes have been associated with the disease, and the molecular mechanisms remain unexplained in most cases. This study aimed to further understand the genetic architecture of BAV. METHODS: A genome-wide association study meta-analysis including 9631 cases among 65 677 participants was performed. Genes were prioritized using transcriptomic analyses based on RNA sequencing in relevant tissues, including human fetal and adult aortic valves. The impact of the knockdown or knockout of 4 candidate genes on cardiac development was verified in zebrafish. A polygenic risk score was developed, its association with BAV was evaluated in an independent cohort, and its association with a wide range of phenotypes (n=976) was evaluated in UK Biobank (n=355 618 individuals). RESULTS: Thirty-six genomic loci were identified, including 32 that were not described previously. Among the prioritized genes, KANK2 and ERBB4 were identified as potentially causal through transcriptomic analyses, colocalization, and Mendelian randomization based on gene expression in human aortic valves (n=484), whereas PRDM6 and STRN were prioritized using similar analyses from aortic (n=326) and left ventricular tissues (n=326), respectively. Targeting 4 candidate genes ( WNT4 , LEF1 , STRN , and KANK2 ) in zebrafish led to disruption in cardiac development. A polygenic risk score was associated with an odds ratio of 2.07 (95% CI, 1.90–2.25; P =5.43×10 -62 ) per SD for BAV and significantly associated with thoracic aortic aneurysm and atrial fibrillation in UK Biobank. CONCLUSIONS: This study supports a significant polygenic contribution to BAV, where the combination of multiple common variants in genes involved in heart morphogenesis disrupts aortic valve development.

We developed a generally applicable method, CRISPR/Cas9-targeted long-read sequencing (CTLR-Seq), to resolve, haplotype-specifically, the large and complex regions in the human genome that had been previously impenetrable to sequencing analysis, such as large segmental duplications (SegDups) and their associated genome rearrangements. CTLR-Seq combines in vitro Cas9-mediated cutting of the genome and pulse-field gel electrophoresis to isolate intact large (i.e., up to 2,000 kb) genomic regions that encompass previously unresolvable genomic sequences. These targets are then sequenced (amplification-free) at high on-target coverage using long-read sequencing, allowing for their complete sequence assembly. We applied CTLR-Seq to the SegDup-mediated rearrangements that constitute the boundaries of, and give rise to, the 22q11.2 Deletion Syndrome (22q11DS), the most common human microdeletion disorder. We then performed de novo assembly to resolve, at base-pair resolution, the full sequence rearrangements and exact chromosomal breakpoints of 22q11.2DS (including all common subtypes). Across multiple patients, we found a high degree of variability for both the rearranged SegDup sequences and the exact chromosomal breakpoint locations, which coincide with various transposons within the 22q11.2 SegDups, suggesting that 22q11DS can be driven by transposon-mediated genome recombination. Guided by CTLR-Seq results from two 22q11DS patients, we performed three-dimensional chromosomal folding analysis for the 22q11.2 SegDups from patient-derived neurons and astrocytes and found chromosome interactions anchored within the SegDups to be both cell type-specific and patient-specific. Lastly, we demonstrated that CTLR-Seq enables cell-type specific analysis of DNA methylation patterns within the deletion haplotype of 22q11DS.

INTRODUCTION The human cerebral cortex has undergone an extraordinary increase in size and complexity during mammalian evolution. Cortical cell lineages are specified in the embryo, and genetic and epidemiological evidence implicates early cortical development in the etiology of neuropsychiatric disorders such as autism spectrum disorder (ASD), intellectual disabilities, and schizophrenia. Most of the disease-implicated genomic variants are located outside of genes, and the interpretation of noncoding mutations is lagging behind owing to limited annotation of functional elements in the noncoding genome. RATIONALE We set out to discover gene-regulatory elements and chart their dynamic activity during prenatal human cortical development, focusing on enhancers, which carry most of the weight upon regulation of gene expression. We longitudinally modeled human brain development using human induced pluripotent stem cell (hiPSC)–derived cortical organoids and compared organoids to isogenic fetal brain tissue. RESULTS Fetal fibroblast–derived hiPSC lines were used to generate cortically patterned organoids and to compare oganoids’ epigenome and transcriptome to that of isogenic fetal brains and external datasets. Organoids model cortical development between 5 and 16 postconception weeks, thus enabling us to study transitions from cortical stem cells to progenitors to early neurons. The greatest changes occur at the transition from stem cells to progenitors. The regulatory landscape encompasses a total set of 96,375 enhancers linked to target genes, with 49,640 enhancers being active in organoids but not in mid-fetal brain, suggesting major roles in cortical neuron specification. Enhancers that gained activity in the human lineage are active in the earliest stages of organoid development, when they target genes that regulate the growth of radial glial cells. Parallel weighted gene coexpression network analysis (WGCNA) of transcriptome and enhancer activities defined a number of modules of coexpressed genes and coactive enhancers, following just six and four global temporal patterns that we refer to as supermodules, likely reflecting fundamental programs in embryonic and fetal brain. Correlations between gene expression and enhancer activity allowed stratifying enhancers into two categories: activating regulators (A-regs) and repressive regulators (R-regs). Several enhancer modules converged with gene modules, suggesting that coexpressed genes are regulated by enhancers with correlated patterns of activity. Furthermore, enhancers active in organoids and fetal brains were enriched for ASD de novo variants that disrupt binding sites of homeodomain, Hes1, NR4A2, Sox3, and NFIX transcription factors. CONCLUSION We validated hiPSC-derived cortical organoids as a suitable model system for studying gene regulation in human embryonic brain development, evolution, and disease. Our results suggest that organoids may reveal how noncoding mutations contribute to ASD etiology. Summary of the study, analyses, and main results. Data were generated for iPSC-derived human telencephalic organoids and isogenic fetal cortex. Organoids modeled embryonic and early fetal cortex and show a larger repertoire of enhancers. Enhancers could be divided into activators and repressors of gene expression. We derived networks of modules and supermodules with correlated gene and enhancer activities, some of which were implicated in autism spectrum disorders (ASD).

Brain mutations, young and old Most neurons that make up the human brain are postmitotic, living and functioning for a very long time without renewal (see the Perspective by Lee). Bae et al. examined the genomes of single neurons from the prenatal developing human brain. Both the type of mutation and the rates of accumulation changed between gastrulation and neurogenesis. These early mutations could be generating useful neuronal diversity or could predispose individuals to later dysfunction. Lodato et al. also found that neurons take on somatic mutations as they age by sequencing single neurons from subjects aged 4 months to 82 years. Somatic mutations accumulated with increasing age and accumulated faster in individuals affected by inborn errors in DNA repair. Postmitotic mutations might only affect one neuron, but the accumulated divergence of genomes across the brain could affect function. Science , this issue p. 550 , p. 555 ; see also p. 521