An analytical framework for whole-genome sequence association studies and its implications for autism spectrum disorder

Whole-genome sequencing analyses suggest novel genetic factors associated with Alzheimer's disease and a cumulative effects model for risk liability

Genome-wide association studies (GWAS) on Alzheimer's disease (AD) have predominantly focused on identifying common variants in Europeans. Here, we performed whole-genome sequencing (WGS) of 1,559 individuals from a Korean AD cohort to identify various genetic variants and biomarkers associated with AD. Our GWAS analysis identified a previously unreported locus for common variants (APCDD1) associated with AD. Our WGS analysis was extended to explore the less-characterized genetic factors contributing to AD risk. We identified rare noncoding variants located in cis-regulatory elements specific to excitatory neurons associated with cognitive impairment. Moreover, structural variation analysis showed that short tandem repeat expansion was associated with an increased risk of AD, and copy number variant at the HPSE2 locus showed borderline statistical significance. APOE ε4 carriers with high polygenic burden or structural variants exhibited severe cognitive impairment and increased amyloid beta levels, suggesting a cumulative effects model of AD risk.

Neurodevelopmental impact of CNV models in ASD: Recent advances and future directions

Autism spectrum disorder (ASD) is a neurodevelopmental disorder characterized by social communication impairments and restricted, repetitive behaviors. ASD exhibits a strong genetic basis, with rare and common genetic variants contributing to its etiology. Copy number variations (CNVs), deletions or duplications of chromosomal segments, have emerged as key contributors to ASD risk. Rare CNVs often demonstrate large effect sizes and can directly cause ASD, while common variants collectively exert subtle influences. Recent advances have identified numerous ASD-associated CNVs, including recurrent loci such as 1q21.1, 2p16.3, 7q11.23, 15q11.2, 15q11-q13, 16p11.2 and 22q11.2. Mouse models carrying these CNVs have provided profound insights into the underlying neurobiological mechanisms. Recent studies integrating transcriptomic, proteomic, and functional imaging approaches have revealed alterations in synaptic function, neuronal differentiation, myelination, metabolic pathways, and circuit connectivity. Notably, investigations leveraging conditional knockout models, high magnetic field MRI, and single-cell analyses highlight disruptions in excitatory-inhibitory balance, white matter integrity, and dynamic gene regulatory networks. Parallel human-based approaches, including iPSC-derived neurons, cerebral organoids, and large-scale single-nucleus sequencing, are combined with animal model data. These integrative strategies promise to refine our understanding of ASD's genetic architecture, bridging the gap between fundamental discoveries in model organisms and clinically relevant biomarkers, subtypes, and therapeutic targets in humans. This review summarizes key findings from recent CNV mouse model studies and highlights emerging technologies applied to human ASD samples. Finally, we outline prospects for translating findings from mouse studies to humans. By illuminating both unique and convergent genetic mechanisms, these advances offer a critical framework for unraveling etiological complexity in ASD.

Whole genome sequence-based association analysis of African American individuals with bipolar disorder and schizophrenia

In studies of individuals of primarily European genetic ancestry, common and low-frequency variants and rare coding variants have been found to be associated with the risk of bipolar disorder (BD) and schizophrenia (SZ). However, less is known for individuals of other genetic ancestries or the role of rare non-coding variants in BD and SZ risk. We performed whole genome sequencing of African American individuals: 1,598 with BD, 3,295 with SZ, and 2,651 unaffected controls (InPSYght study). We increased power by incorporating 14,812 jointly called psychiatrically unscreened ancestry-matched controls from the Trans-Omics for Precision Medicine (TOPMed) Program for a total of 17,463 controls. To identify variants and sets of variants associated with BD and/or SZ, we performed single-variant tests, gene-based tests for singleton protein truncating variants, and rare and low-frequency variant annotation-based tests with conservation and universal chromatin states and sliding windows. We found suggestive evidence of BD association with single-variants on chromosome 18 and of lower BD risk associated with rare and low-frequency variants on chromosome 11 in a region with multiple BD GWAS loci, using a sliding window approach. We also found that chromatin and conservation state tests can be used to detect differential calling of variants in controls sequenced at different centers and to assess the effectiveness of sequencing metric covariate adjustments. Our findings reinforce the need for continued whole genome sequencing in additional samples of African American individuals and more comprehensive functional annotation of non-coding variants.

The genetic landscape of early-onset Alzheimer's disease in China

INTRODUCTION

Research on somatic and germline mutations in Chinese individuals with early-onset Alzheimer's disease (EOAD) has been limited.

METHODS

We conducted whole-genome sequencing of blood DNA from 108 patients with EOAD and 116 controls. The analysis included somatic and germline mutations across coding and non-coding regions, mutational signature determination, pathway enrichment identification, and predictive model.

RESULTS

The mutational burden was significantly higher in the EOAD group compared to the control group. The prevalence of single-base substitution signature 5, which is strongly associated with aging, was much higher in patients with EOAD than in controls. EOAD-specific somatic mutations were identified in genes such as MIR31HG, TUBB4B, and APP. Germline mutations in DOCK3, PCSK5, and PDE4D were significantly associated with age of dementia onset. Furthermore, a predictive model comprising 15 mutations demonstrated an area under the curve of 0.78.

DISCUSSION

The accumulation of senescence-related somatic mutations may increase the risk of developing EOAD.

HIGHLIGHTS

Whole genome sequencing was used to find somatic and germline mutations in Chinese individuals with early-onset Alzheimer's disease (EOAD). Total number and burden of blood somatic mutations were significantly higher. The prevalence of single-base substitution signature 5 was notably elevated in EOAD. EOAD-specific somatic mutations were identified in MIR31HG, TUBB4B, and APP. DOCK3, PCSK5, and PDE4D germline mutations were associated with the age of EOAD onset.

Diversity and consequences of structural variation in the human genome

The biomedical community is increasingly invested in capturing all genetic variants across human genomes, interpreting their functional consequences and translating these findings to the clinic. A crucial component of this endeavour is the discovery and characterization of structural variants (SVs), which are ubiquitous in the human population, heterogeneous in their mutational processes, key substrates for evolution and adaptation, and profound drivers of human disease. The recent emergence of new technologies and the remarkable scale of sequence-based population studies have begun to crystalize our understanding of SVs as a mutational class and their widespread influence across phenotypes. In this Review, we summarize recent discoveries and new insights into SVs in the human genome in terms of their mutational patterns, population genetics, functional consequences, and impact on human traits and disease. We conclude by outlining three frontiers to be explored by the field over the next decade.

Generation and Characterization of a Knockout Mouse of an Enhancer of EBF3

Genomic studies of autism and other neurodevelopmental disorders have identified several relevant protein-coding and noncoding variants. One gene with an excess of protein-coding de novo variants is EBF3 that also is the gene underlying the Hypotonia, Ataxia, and Delayed Development Syndrome (HADDS). In previous work, we have identified noncoding de novo variants in an enhancer of EBF3 called hs737 and further showed that there was an enrichment of deletions of this enhancer in individuals with neurodevelopmental disorders. In this present study, we generated a novel mouse line that deletes the highly conserved, orthologous mouse region of hs737 within the Rr169617 regulatory region, and characterized the molecular and phenotypic aspects of this mouse model. This line contains a 1,160 bp deletion within Rr169617 and through heterozygous crosses we found a deviation from Mendelian expectation (p = 0.02) with a significant depletion of the deletion allele (p = 5.8 × 10-4). Rr169617 +/- mice had a reduction of Ebf3 expression by 10% and Rr169617 -/- mice had a reduction of Ebf3 expression by 20%. Differential expression analyses in E12.5 forebrain, midbrain, and hindbrain in Rr169617 +/+ versus Rr169617 -/- mice identified dysregulated genes including histone genes (i.e., Hist1h1e, Hist1h2bk, Hist1h3i, Hist1h2ao) and other brain development related genes (e.g., Chd5, Ntng1). A priori phenotyping analysis (open field, hole board and light/dark transition) identified sex-specific differences in behavioral traits when comparing Rr169617 -/- males versus females; whereby, males were observed to be less mobile, move slower, and spend more time in the dark. Furthermore, both sexes when homozygous for the enhancer deletion displayed body composition differences when compared to wild-type mice. Overall, we show that deletion within Rr169617 reduces the expression of Ebf3 and results in phenotypic outcomes consistent with potential sex specific behavioral differences. This enhancer deletion line provides a valuable resource for others interested in noncoding regions in neurodevelopmental disorders and/or those interested in the gene regulatory network downstream of Ebf3.

Rare germline structural variants increase risk for pediatric solid tumors

Pediatric solid tumors are a leading cause of childhood disease mortality. In this work, we examined germline structural variants (SVs) as risk factors for pediatric extracranial solid tumors using germline genome sequencing of 1765 affected children, their 943 unaffected parents, and 6665 adult controls. We discovered a sex-biased association between very large (>1 megabase) germline chromosomal abnormalities and increased risk of solid tumors in male children. The overall impact of germline SVs was greatest in neuroblastoma, where we uncovered burdens of ultrarare SVs that cause loss of function of highly expressed, mutationally constrained genes, as well as noncoding SVs predicted to disrupt chromatin domain boundaries. Collectively, we estimate that rare germline SVs explain 1.1 to 5.6% of pediatric cancer liability, establishing them as an important component of disease predisposition.

The landscape of structural variation in aye-ayes (Daubentonia madagascariensis)

Aye-ayes (Daubentonia madagascariensis) are one of the 25 most critically endangered primate species in the world. Endemic to Madagascar, their small and highly fragmented populations make them particularly vulnerable to both genetic disease and anthropogenic environmental changes. Over the past decade, conservation genomic efforts have largely focused on inferring and monitoring population structure based on single nucleotide variants to identify and protect critical areas of genetic diversity. However, the recent release of a highly contiguous genome assembly allows, for the first time, for the study of structural genomic variation (deletions, duplications, insertions, and inversions) which are likely to impact a substantial proportion of the species' genome. Based on whole-genome, short-read sequencing data from 14 individuals, >1,000 high-confidence autosomal structural variants were detected, affecting ~240 kb of the aye-aye genome. The majority of these variants (>85%) were deletions shorter than 200 bp, consistent with the notion that longer structural mutations are often associated with strongly deleterious fitness effects. For example, two deletions longer than 850 bp located within disease-linked genes were predicted to impose substantial fitness deficits owing to a resulting frameshift and gene fusion, respectively; whereas several other major effect variants outside of coding regions are likely to impact gene regulatory landscapes. Taken together, this first glimpse into the landscape of structural variation in aye-ayes will enable future opportunities to advance our understanding of the traits impacting the fitness of this endangered species, as well as allow for enhanced evolutionary comparisons across the full primate clade.

Key Synaptic Pathology in Autism Spectrum Disorder: Genetic Mechanisms and Recent Advances

Autism spectrum disorder (ASD) is a neurodevelopmental disorder characterized by impaired social interactions and verbal communication, accompanied by symptoms of restricted and repetitive patterns of behavior or interest. Over the past 30 years, the morbidity of ASD has increased in most areas of the world. Although the pathogenesis of ASD is not fully understood, it has been associated with over 1000 genes or genomic loci, indicating the importance and complexity of the genetic mechanisms involved. This review focuses on the synaptic pathology of ASD and particularly on genetic variants involved in synaptic structure and functions. These include SHANK, NLGN, NRXN, FMR1, and MECP2 as well as other potentially novel genes such as CHD8, CHD2, and SYNGAP1 that could be core elements in ASD pathogenesis. Here, we summarize several pathological pathways supporting the hypothesis that synaptic pathology caused by genetic mutations may be the pathogenic basis for ASD.

Whole genome sequencing analysis identifies sex differences of familial pattern contributing to phenotypic diversity in autism

BACKGROUND

Whole-genome sequencing (WGS) analyses have found higher genetic burden in autistic females compared to males, supporting higher liability threshold in females. However, genomic evidence of sex differences has been limited to European ancestry to date and little is known about how genetic variation leads to autism-related traits within families across sex.

METHODS

To address this gap, we present WGS data of Korean autism families (n = 2255) and a Korean general population sample (n = 2500), the largest WGS data of East Asian ancestry. We analyzed sex differences in genetic burden and compared with cohorts of European ancestry (n = 15,839). Further, with extensively collected family-wise Korean autism phenotype data (n = 3730), we investigated sex differences in phenotypic scores and gene-phenotype associations within family.

RESULTS

We observed robust female enrichment of de novo protein-truncating variants in autistic individuals across cohorts. However, sex differences in polygenic burden varied across cohorts and we found that the differential proportion of comorbid intellectual disability and severe autism symptoms mainly drove these variations. In siblings, males of autistic females exhibited the most severe social communication deficits. Female siblings exhibited lower phenotypic severity despite the higher polygenic burden than male siblings. Mothers also showed higher tolerance for polygenic burden than fathers, supporting higher liability threshold in females.

CONCLUSIONS

Our findings indicate that genetic liability in autism is both sex- and phenotype-dependent, expanding the current understanding of autism's genetic complexity. Our work further suggests that family-based assessments of sex differences can help unravel underlying sex-differential liability in autism.

Gene therapy for polygenic or complex diseases

Gene therapy utilizes nucleic acid drugs to treat diseases, encompassing gene supplementation, gene replacement, gene silencing, and gene editing. It represents a distinct therapeutic approach from traditional medications and introduces novel strategies for genetic disorders. Over the past two decades, significant advancements have been made in the field of gene therapy, leading to the approval of various gene therapy drugs. Gene therapy was initially employed for treating genetic diseases and cancers, particularly monogenic conditions classified as orphan diseases due to their low prevalence rates; however, polygenic or complex diseases exhibit higher incidence rates within populations. Extensive research on the etiology of polygenic diseases has unveiled new therapeutic targets that offer fresh opportunities for their treatment. Building upon the progress achieved in gene therapy for monogenic diseases and cancers, extending its application to polygenic or complex diseases would enable targeting a broader range of patient populations. This review aims to discuss the strategies of gene therapy, methods of gene editing (mainly CRISPR-CAS9), and carriers utilized in gene therapy, and highlight the applications of gene therapy in polygenic or complex diseases focused on applications that have either entered clinical stages or are currently undergoing clinical trials.

Genetic modifiers and ascertainment drive variable expressivity of complex disorders

Variable expressivity of disease-associated variants implies a role for secondary variants that modify clinical features. We assessed the effects of modifier variants towards clinical outcomes of 2,252 individuals with primary variants. Among 132 families with the 16p12.1 deletion, distinct rare and common variant classes conferred risk for specific developmental features, including short tandem repeats for neurological defects and SNVs for microcephaly, while additional disease-associated variants conferred multiple genetic diagnoses. Within disease and population cohorts of 773 individuals with the 16p12.1 deletion, we found opposing effects of secondary variants towards clinical features across ascertainments. Additional analysis of 1,479 probands with other primary variants, such as 16p11.2 deletion and CHD8 variants, and 1,084 without primary variants, showed that phenotypic associations differed by primary variant context and were influenced by synergistic interactions between primary and secondary variants. Our study provides a paradigm to dissect the genomic architecture of complex disorders towards personalized treatment.

Decomposition of phenotypic heterogeneity in autism reveals distinct and coherent genetic programs

Unraveling the phenotypic and genetic complexity of autism is extremely challenging yet critical for understanding the biology, inheritance, trajectory, and clinical manifestations of the many forms of the condition. Here, we leveraged broad phenotypic data from a large cohort with matched genetics to characterize classes of autism and their patterns of core, associated, and co-occurring traits, ultimately demonstrating that phenotypic patterns are associated with distinct genetic and molecular programs. We used a generative mixture modeling approach to identify robust, clinically-relevant classes of autism which we validate and replicate in a large independent cohort. We link the phenotypic findings to distinct patterns of de novo and inherited variation which emerge from the deconvolution of these genetic signals, and demonstrate that class-specific common variant scores strongly align with clinical outcomes. We further provide insights into the distinct biological pathways and processes disrupted by the sets of mutations in each class. Remarkably, we discover class-specific differences in the developmental timing of genes that are dysregulated, and these temporal patterns correspond to clinical milestone and outcome differences between the classes. These analyses embrace the phenotypic complexity of children with autism, unraveling genetic and molecular programs underlying their heterogeneity and suggesting specific biological dysregulation patterns and mechanistic hypotheses.

Investigating copy number variants in schizophrenia pedigrees using a new consensus pipeline called PECAN

Copy number variants (CNVs) have been implicated in many human diseases, including psychiatric disorders. Whole genome sequencing offers advantages in CNV calling compared to previous array-based methods. Here we present a robust and transparent CNV calling pipeline, PECAN (PEdigree Copy number vAriaNt calling), for short-read, whole genome sequencing data, comprised of a novel combination of four calling methods and structural variant genotyping. This method is scalable and can incorporate pedigree information to retain lower-confidence CNVs that would otherwise be discarded. We have robustly benchmarked PECAN using gold-standard CNV calls for two well-established evaluation samples, NA12878 and HG002, showing that PECAN performs with high precision and recall on both datasets, outperforming another pedigree-based CNV calling pipeline. As part of this work, we provide a list of high-confidence gold standard CNVs for the NA12878 reference sample, curated from multiple studies. We applied PECAN to a collection of pedigrees multiply affected with schizophrenia and identified a rare deletion that perfectly co-segregates with schizophrenia in one of the pedigrees. The CNV overlaps the gene PITRM1, which has been implicated in a complex phenotype including ataxia, developmental delay, and schizophrenia-like episodes in affected adults.

A Comparison of Structural Variant Calling from Short-Read and Nanopore-Based Whole-Genome Sequencing Using Optical Genome Mapping as a Benchmark

The identification of structural variants (SVs) in genomic data represents an ongoing challenge because of difficulties in reliable SV calling leading to reduced sensitivity and specificity. We prepared high-quality DNA from 9 parent-child trios, who had previously undergone short-read whole-genome sequencing (Illumina platform) as part of the Genomics England 100,000 Genomes Project. We reanalysed the genomes using both Bionano optical genome mapping (OGM; 8 probands and one trio) and Nanopore long-read sequencing (Oxford Nanopore Technologies [ONT] platform; all samples). To establish a "truth" dataset, we asked whether rare proband SV calls (n = 234) made by the Bionano Access (version 1.6.1)/Solve software (version 3.6.1_11162020) could be verified by individual visualisation using the Integrative Genomics Viewer with either or both of the Illumina and ONT raw sequence. Of these, 222 calls were verified, indicating that Bionano OGM calls have high precision (positive predictive value 95%). We then asked what proportion of the 222 true Bionano SVs had been identified by SV callers in the other two datasets. In the Illumina dataset, sensitivity varied according to variant type, being high for deletions (115/134; 86%) but poor for insertions (13/58; 22%). In the ONT dataset, sensitivity was generally poor using the original Sniffles variant caller (48% overall) but improved substantially with use of Sniffles2 (36/40; 90% and 17/23; 74% for deletions and insertions, respectively). In summary, we show that the precision of OGM is very high. In addition, when applying the Sniffles2 caller, the sensitivity of SV calling using ONT long-read sequence data outperforms Illumina sequencing for most SV types.

Short tandem repeat expansions in cortical layer-specific genes implicate in phenotypic severity and adaptability of autism spectrum disorder

AIM

Short tandem repeats (STRs) are repetitive DNA sequences and highly mutable in various human disorders. While the involvement of STRs in various genetic disorders has been extensively studied, their role in autism spectrum disorder (ASD) remains largely unexplored. In this study, we aimed to investigate genetic association of STR expansions with ASD using whole genome sequencing (WGS) and identify risk loci associated with ASD phenotypes.

METHODS

We analyzed WGS data of 634 ASD families and performed genome-wide evaluation for 12,929 STR loci. We found rare STR expansions that exceeded normal repeat lengths in autism cases compared to unaffected controls. By integrating single cell RNA and ATAC sequencing datasets of human postmortem brains, we prioritized STR loci in genes specifically expressed in cortical development stages. A deep learning method was used to predict functionality of ASD-associated STR loci.

RESULTS

In ASD cases, rare STR expansions predominantly occurred in early cortical layer-specific genes involved in neurodevelopment, highlighting the cellular specificity of STR-associated genes in ASD risk. Leveraging deep learning prediction models, we demonstrated that these STR expansions disrupted the regulatory activity of enhancers and promoters, suggesting a potential mechanism through which they contribute to ASD pathogenesis. We found that individuals with ASD-associated STR expansions exhibited more severe ASD phenotypes and diminished adaptability compared to non-carriers.

CONCLUSION

Short tandem repeat expansions in cortical layer-specific genes are associated with ASD and could potentially be a risk genetic factor for ASD. Our study is the first to show evidence of STR expansion associated with ASD in an under-investigated population.

Towards understanding sex differences in autism spectrum disorders

Autism Spectrum Disorder (ASD) is a neurodevelopmental condition characterized by social deficits, repetitive behaviours and lack of empathy. Its significant genetic heritability and potential comorbidities often lead to diagnostic and therapeutic challenges. This review addresses the biological basis of ASD, focusing on the sex differences in gene expression and hormonal influences. ASD is more commonly diagnosed in males at a ratio of 4:1, indicating a potential oversight in female-specific ASD research and a risk of underdiagnosis in females. We consider how ASD manifests differently across sexes by exploring differential gene expression in female and male brains and consider how variations in steroid hormones influence ASD characteristics. Synaptic function, including excitation/inhibition ratio imbalance, is influenced by gene mutations and this is explored as a key factor in the cognitive and behavioural manifestations of ASD. We also discuss the role of micro RNAs (miRNAs) and highlight a novel mutation in miRNA-873, which affects a suite of key synaptic genes, neurexin, neuroligin, SHANK and post-synaptic density proteins, implicated in the pathology of ASD. Our review suggests that genetic predisposition, sex differences in brain gene expression, and hormonal factors significantly contribute to the presentation, identification and severity of ASD, necessitating sex-specific considerations in diagnosis and treatments. These findings advocate for personalized interventions to improve the outcomes for individuals with ASD.

The Importance of Large-Scale Genomic Studies to Unravel Genetic Risk Factors for Autism

Autism spectrum disorder (ASD) is a common and highly heritable neurodevelopmental disorder. During the last 15 years, advances in genomic technologies and the availability of increasingly large patient cohorts have greatly expanded our knowledge of the genetic architecture of ASD and its neurobiological mechanisms. Over two hundred risk regions and genes carrying rare de novo and transmitted high-impact variants have been identified. Additionally, common variants with small individual effect size are also important, and a number of loci are now being uncovered. At the same time, these new insights have highlighted ongoing challenges. In this perspective article, we summarize developments in ASD genetic research and address the enormous impact of large-scale genomic initiatives on ASD gene discovery.

Developmental isoform diversity in the human neocortex informs neuropsychiatric risk mechanisms

RNA splicing is highly prevalent in the brain and has strong links to neuropsychiatric disorders; yet, the role of cell type-specific splicing and transcript-isoform diversity during human brain development has not been systematically investigated. In this work, we leveraged single-molecule long-read sequencing to deeply profile the full-length transcriptome of the germinal zone and cortical plate regions of the developing human neocortex at tissue and single-cell resolution. We identified 214,516 distinct isoforms, of which 72.6% were novel (not previously annotated in Gencode version 33), and uncovered a substantial contribution of transcript-isoform diversity-regulated by RNA binding proteins-in defining cellular identity in the developing neocortex. We leveraged this comprehensive isoform-centric gene annotation to reprioritize thousands of rare de novo risk variants and elucidate genetic risk mechanisms for neuropsychiatric disorders.

CWAS-Plus: estimating category-wide association of rare noncoding variation from whole-genome sequencing data with cell-type-specific functional data

Variants in cis-regulatory elements link the noncoding genome to human pathology; however, detailed analytic tools for understanding the association between cell-level brain pathology and noncoding variants are lacking. CWAS-Plus, adapted from a Python package for category-wide association testing (CWAS), enhances noncoding variant analysis by integrating both whole-genome sequencing (WGS) and user-provided functional data. With simplified parameter settings and an efficient multiple testing correction method, CWAS-Plus conducts the CWAS workflow 50 times faster than CWAS, making it more accessible and user-friendly for researchers. Here, we used a single-nuclei assay for transposase-accessible chromatin with sequencing to facilitate CWAS-guided noncoding variant analysis at cell-type-specific enhancers and promoters. Examining autism spectrum disorder WGS data (n = 7280), CWAS-Plus identified noncoding de novo variant associations in transcription factor binding sites within conserved loci. Independently, in Alzheimer's disease WGS data (n = 1087), CWAS-Plus detected rare noncoding variant associations in microglia-specific regulatory elements. These findings highlight CWAS-Plus's utility in genomic disorders and scalability for processing large-scale WGS data and in multiple-testing corrections. CWAS-Plus and its user manual are available at https://github.com/joonan-lab/cwas/ and https://cwas-plus.readthedocs.io/en/latest/, respectively.

Noncoding de novo mutations in SCN2A are associated with autism spectrum disorders

Coding de novo mutations (DNMs) contribute to the risk for autism spectrum disorders (ASD), but the contribution of noncoding DNMs remains relatively unexplored. Here we use whole genome sequencing (WGS) data of 12,411 individuals (including 3,508 probands and 2,218 unaffected siblings) from 3,357 families collected in Simons Foundation Powering Autism Research for Knowledge (SPARK) to detect DNMs associated with ASD, while examining Simons Simplex Collection (SSC) with 6383 individuals from 2274 families to replicate the results. For coding DNMs, SCN2A reached exome-wide significance (p=2.06×10-11) in SPARK. The 618 known dominant ASD genes as a group are strongly enriched for coding DNMs in cases than sibling controls (fold change=1.51, p =1.13×10-5 for SPARK; fold change=1.86, p =2.06×10-9 for SSC). For noncoding DNMs, we used two methods to assess statistical significance: a point-based test that analyzes sites with a Combined Annotation Dependent Depletion (CADD) score ≥15, and a segment-based test that analyzes 1kb genomic segments with segment-specific background mutation rates (inferred from expected rare mutations in Gnocchi genome constraint scores). The point-based test identified SCN2A as marginally significant (p=6.12×10-4) in SPARK, yet segment-based test identified CSMD1, RBFOX1 and CHD13 as exome-wide significant. We did not identify significant enrichment of noncoding DNMs (in all 1kb segments or those with Gnocchi>4) in the 618 known ASD genes as a group in cases than sibling controls. When combining evidence from both coding and noncoding DNMs, we found that SCN2A with 11 coding and 5 noncoding DNMs exhibited the strongest significance (p=4.15×10-13). In summary, we identified both coding and noncoding DNMs in SCN2A associated with ASD, while nominating additional candidates for further examination in future studies.

Rare germline structural variants increase risk for pediatric solid tumors

Pediatric solid tumors are rare malignancies that represent a leading cause of death by disease among children in developed countries. The early age-of-onset of these tumors suggests that germline genetic factors are involved, yet conventional germline testing for short coding variants in established predisposition genes only identifies pathogenic events in 10-15% of patients. Here, we examined the role of germline structural variants (SVs)-an underexplored form of germline variation-in pediatric extracranial solid tumors using germline genome sequencing of 1,766 affected children, their 943 unaffected relatives, and 6,665 adult controls. We discovered a sex-biased association between very large (>1 megabase) germline chromosomal abnormalities and a four-fold increased risk of solid tumors in male children. The overall impact of germline SVs was greatest in neuroblastoma, where we revealed burdens of ultra-rare SVs that cause loss-of-function of highly expressed, mutationally intolerant, neurodevelopmental genes, as well as noncoding SVs predicted to disrupt three-dimensional chromatin domains in neural crest-derived tissues. Collectively, our results implicate rare germline SVs as a predisposing factor to pediatric solid tumors that may guide future studies and clinical practice.

Contribution of de novo retroelements to birth defects and childhood cancers

Insertion of active retroelements-L1s, Alus, and SVAs-can disrupt proper genome function and lead to various disorders including cancer. However, the role of de novo retroelements (DNRTs) in birth defects and childhood cancers has not been well characterized due to the lack of adequate data and efficient computational tools. Here, we examine whole-genome sequencing data of 3,244 trios from 12 birth defect and childhood cancer cohorts in the Gabriella Miller Kids First Pediatric Research Program. Using an improved version of our tool xTea (x-Transposable element analyzer) that incorporates a deep-learning module, we identified 162 DNRTs, as well as 2 pseudogene insertions. Several variants are likely to be causal, such as a de novo Alu insertion that led to the ablation of a whole exon in the NF1 gene in a proband with brain tumor. We observe a high de novo SVA insertion burden in both high-intolerance loss-of-function genes and exons as well as more frequent de novo Alu insertions of paternal origin. We also identify potential mosaic DNRTs from embryonic stages. Our study reveals the important roles of DNRTs in causing birth defects and predisposition to childhood cancers.

CWAS-Plus: Estimating category-wide association of rare noncoding variation from whole-genome sequencing data with cell-type-specific functional data

Variants in cis-regulatory elements link the noncoding genome to human brain pathology; however, detailed analytic tools for understanding the association between cell-level brain pathology and noncoding variants are lacking. CWAS-Plus, adapted from a Python package for category-wide association testing (CWAS) employs both whole-genome sequencing and user-provided functional data to enhance noncoding variant analysis, with a faster and more efficient execution of the CWAS workflow. Here, we used single-nuclei assay for transposase-accessible chromatin with sequencing to facilitate CWAS-guided noncoding variant analysis at cell-type specific enhancers and promoters. Examining autism spectrum disorder whole-genome sequencing data (n = 7,280), CWAS-Plus identified noncoding de novo variant associations in transcription factor binding sites within conserved loci. Independently, in Alzheimer's disease whole-genome sequencing data (n = 1,087), CWAS-Plus detected rare noncoding variant associations in microglia-specific regulatory elements. These findings highlight CWAS-Plus's utility in genomic disorders and scalability for processing large-scale whole-genome sequencing data and in multiple-testing corrections. CWAS-Plus and its user manual are available at https://github.com/joonan-lab/cwas/ and https://cwas-plus.readthedocs.io/en/latest/, respectively.

Identifying associations of de novo noncoding variants with autism through integration of gene expression, sequence and sex information

Whole-genome sequencing (WGS) data is facilitating genome-wide identification of rare noncoding variants, while elucidating their roles in disease remains challenging. Towards this end, we first revisit a reported significant brain-related association signal of autism spectrum disorder (ASD) detected from de novo noncoding variants attributed to deep-learning and show that local GC content can capture similar association signals. We further show that the association signal appears driven by variants from male proband-female sibling pairs that are upstream of assigned genes. We then develop Expression Neighborhood Sequence Association Study (ENSAS), which utilizes gene expression correlations and sequence information, to more systematically identify phenotype-associated variant sets. Applying ENSAS to the same set of de novo variants, we identify gene expression-based neighborhoods showing significant ASD association signal, enriched for synapse-related gene ontology terms. For these top neighborhoods, we also identify chromatin states annotations of variants that are predictive of the proband-sibling local GC content differences. Our work provides new insights into associations of non-coding de novo mutations in ASD and presents an analytical framework applicable to other phenotypes.

Epistasis between mutator alleles contributes to germline mutation spectrum variability in laboratory mice

Maintaining germline genome integrity is essential and enormously complex. Although many proteins are involved in DNA replication, proofreading, and repair, mutator alleles have largely eluded detection in mammals. DNA replication and repair proteins often recognize sequence motifs or excise lesions at specific nucleotides. Thus, we might expect that the spectrum of de novo mutations - the frequencies of C>T, A>G, etc. - will differ between genomes that harbor either a mutator or wild-type allele. Previously, we used quantitative trait locus mapping to discover candidate mutator alleles in the DNA repair gene Mutyh that increased the C>A germline mutation rate in a family of inbred mice known as the BXDs (Sasani et al., 2022, Ashbrook et al., 2021). In this study we developed a new method to detect alleles associated with mutation spectrum variation and applied it to mutation data from the BXDs. We discovered an additional C>A mutator locus on chromosome 6 that overlaps Ogg1, a DNA glycosylase involved in the same base-excision repair network as Mutyh (David et al., 2007). Its effect depends on the presence of a mutator allele near Mutyh, and BXDs with mutator alleles at both loci have greater numbers of C>A mutations than those with mutator alleles at either locus alone. Our new methods for analyzing mutation spectra reveal evidence of epistasis between germline mutator alleles and may be applicable to mutation data from humans and other model organisms.

Topologically associating domains define the impact of de novo promoter variants on autism spectrum disorder risk

Whole-genome sequencing (WGS) studies of autism spectrum disorder (ASD) have demonstrated the roles of rare promoter de novo variants (DNVs). However, most promoter DNVs in ASD are not located immediately upstream of known ASD genes. In this study analyzing WGS data of 5,044 ASD probands, 4,095 unaffected siblings, and their parents, we show that promoter DNVs within topologically associating domains (TADs) containing ASD genes are significantly and specifically associated with ASD. An analysis considering TADs as functional units identified specific TADs enriched for promoter DNVs in ASD and indicated that common variants in these regions also confer ASD heritability. Experimental validation using human induced pluripotent stem cells (iPSCs) showed that likely deleterious promoter DNVs in ASD can influence multiple genes within the same TAD, resulting in overall dysregulation of ASD-associated genes. These results highlight the importance of TADs and gene-regulatory mechanisms in better understanding the genetic architecture of ASD.

Intellectual Disability and Behavioral Deficits Linked to CYFIP1 Missense Variants Disrupting Actin Polymerization

BACKGROUND

15q11.2 deletions and duplications have been linked to autism spectrum disorder, schizophrenia, and intellectual disability. Recent evidence suggests that dysfunctional CYFIP1 (cytoplasmic FMR1 interacting protein 1) contributes to the clinical phenotypes observed in individuals with 15q11.2 deletion/duplication syndrome. CYFIP1 plays crucial roles in neuronal development and brain connectivity, promoting actin polymerization and regulating local protein synthesis. However, information about the impact of single nucleotide variants in CYFIP1 on neurodevelopmental disorders is limited.

METHODS

Here, we report a family with 2 probands exhibiting intellectual disability, autism spectrum disorder, spastic tetraparesis, and brain morphology defects and who carry biallelic missense point mutations in the CYFIP1 gene. We used skin fibroblasts from one of the probands, the parents, and typically developing individuals to investigate the effect of the variants on the functionality of CYFIP1. In addition, we generated Drosophila knockin mutants to address the effect of the variants in vivo and gain insight into the molecular mechanism that underlies the clinical phenotype.

RESULTS

Our study revealed that the 2 missense variants are in protein domains responsible for maintaining the interaction within the wave regulatory complex. Molecular and cellular analyses in skin fibroblasts from one proband showed deficits in actin polymerization. The fly model for these mutations exhibited abnormal brain morphology and F-actin loss and recapitulated the core behavioral symptoms, such as deficits in social interaction and motor coordination.

CONCLUSIONS

Our findings suggest that the 2 CYFIP1 variants contribute to the clinical phenotype in the probands that reflects deficits in actin-mediated brain development processes.

Finding Needles in the Haystack: Strategies for Uncovering Noncoding Regulatory Variants

Despite accumulating evidence implicating noncoding variants in human diseases, unraveling their functionality remains a significant challenge. Systematic annotations of the regulatory landscape and the growth of sequence variant data sets have fueled the development of tools and methods to identify causal noncoding variants and evaluate their regulatory effects. Here, we review the latest advances in the field and discuss potential future research avenues to gain a more in-depth understanding of noncoding regulatory variants.

A Cre-dependent massively parallel reporter assay allows for cell-type specific assessment of the functional effects of non-coding elements in vivo

The function of regulatory elements is highly dependent on the cellular context, and thus for understanding the function of elements associated with psychiatric diseases these would ideally be studied in neurons in a living brain. Massively Parallel Reporter Assays (MPRAs) are molecular genetic tools that enable functional screening of hundreds of predefined sequences in a single experiment. These assays have not yet been adapted to query specific cell types in vivo in a complex tissue like the mouse brain. Here, using a test-case 3'UTR MPRA library with genomic elements containing variants from autism patients, we developed a method to achieve reproducible measurements of element effects in vivo in a cell type-specific manner, using excitatory cortical neurons and striatal medium spiny neurons as test cases. This targeted technique should enable robust, functional annotation of genetic elements in the cellular contexts most relevant to psychiatric disease.

Developmental isoform diversity in the human neocortex informs neuropsychiatric risk mechanisms

RNA splicing is highly prevalent in the brain and has strong links to neuropsychiatric disorders, yet the role of cell-type-specific splicing or transcript-isoform diversity during human brain development has not been systematically investigated. Here, we leveraged single-molecule long-read sequencing to deeply profile the full-length transcriptome of the germinal zone (GZ) and cortical plate (CP) regions of the developing human neocortex at tissue and single-cell resolution. We identified 214,516 unique isoforms, of which 72.6% are novel (unannotated in Gencode-v33), and uncovered a substantial contribution of transcript-isoform diversity, regulated by RNA binding proteins, in defining cellular identity in the developing neocortex. We leveraged this comprehensive isoform-centric gene annotation to re-prioritize thousands of rare de novo risk variants and elucidate genetic risk mechanisms for neuropsychiatric disorders.

One-Sentence Summary

A cell-specific atlas of gene isoform expression helps shape our understanding of brain development and disease.

Structured Abstract

INTRODUCTION: The development of the human brain is regulated by precise molecular and genetic mechanisms driving spatio-temporal and cell-type-specific transcript expression programs. Alternative splicing, a major mechanism increasing transcript diversity, is highly prevalent in the human brain, influences many aspects of brain development, and has strong links to neuropsychiatric disorders. Despite this, the cell-type-specific transcript-isoform diversity of the developing human brain has not been systematically investigated.RATIONALE: Understanding splicing patterns and isoform diversity across the developing neocortex has translational relevance and can elucidate genetic risk mechanisms in neurodevelopmental disorders. However, short-read sequencing, the prevalent technology for transcriptome profiling, is not well suited to capturing alternative splicing and isoform diversity. To address this, we employed third-generation long-read sequencing, which enables capture and sequencing of complete individual RNA molecules, to deeply profile the full-length transcriptome of the germinal zone (GZ) and cortical plate (CP) regions of the developing human neocortex at tissue and single-cell resolution.RESULTS: We profiled microdissected GZ and CP regions of post-conception week (PCW) 15-17 human neocortex in bulk and at single-cell resolution across six subjects using high-fidelity long-read sequencing (PacBio IsoSeq). We identified 214,516 unique isoforms, of which 72.6% were novel (unannotated in Gencode), and >7,000 novel exons, expanding the proteome by 92,422 putative proteoforms. We uncovered thousands of isoform switches during cortical neurogenesis predicted to impact RNA regulatory domains or protein structure and implicating previously uncharacterized RNA-binding proteins in cellular identity and neuropsychiatric disease. At the single-cell level, early-stage excitatory neurons exhibited the greatest isoform diversity, and isoform-centric single-cell clustering led to the identification of previously uncharacterized cell states. We systematically assessed the contribution of transcriptomic features, and localized cell and spatio-temporal transcript expression signatures across neuropsychiatric disorders, revealing predominant enrichments in dynamic isoform expression and utilization patterns and that the number and complexity of isoforms per gene is strongly predictive of disease. Leveraging this resource, we re-prioritized thousands of rare de novo risk variants associated with autism spectrum disorders (ASD), intellectual disability (ID), and neurodevelopmental disorders (NDDs), more broadly, to potentially more severe consequences and revealed a larger proportion of cryptic splice variants with the expanded transcriptome annotation provided in this study.CONCLUSION: Our study offers a comprehensive landscape of isoform diversity in the human neocortex during development. This extensive cataloging of novel isoforms and splicing events sheds light on the underlying mechanisms of neurodevelopmental disorders and presents an opportunity to explore rare genetic variants linked to these conditions. The implications of our findings extend beyond fundamental neuroscience, as they provide crucial insights into the molecular basis of developmental brain disorders and pave the way for targeted therapeutic interventions. To facilitate exploration of this dataset we developed an online portal ( https://sciso.gandallab.org/ ).

Rare variants found in clinical gene panels illuminate the genetic and allelic architecture of orofacial clefting

PURPOSE

Orofacial clefts (OFCs) are common birth defects including cleft lip, cleft lip and palate, and cleft palate. OFCs have heterogeneous etiologies, complicating clinical diagnostics because it is not always apparent if the cause is Mendelian, environmental, or multifactorial. Sequencing is not currently performed for isolated or sporadic OFCs; therefore, we estimated the diagnostic yield for 418 genes in 841 cases and 294 controls.

METHODS

We evaluated 418 genes using genome sequencing and curated variants to assess their pathogenicity using American College of Medical Genetics criteria.

RESULTS

9.04% of cases and 1.02% of controls had "likely pathogenic" variants (P < .0001), which was almost exclusively driven by heterozygous variants in autosomal genes. Cleft palate (17.6%) and cleft lip and palate (9.09%) cases had the highest yield, whereas cleft lip cases had a 2.80% yield. Out of 39 genes with likely pathogenic variants, 9 genes, including CTNND1 and IRF6, accounted for more than half of the yield (4.64% of cases). Most variants (61.8%) were "variants of uncertain significance", occurring more frequently in cases (P = .004), but no individual gene showed a significant excess of variants of uncertain significance.

CONCLUSION

These results underscore the etiological heterogeneity of OFCs and suggest sequencing could reduce the diagnostic gap in OFCs.

Systematic evaluation of genome sequencing for the diagnostic assessment of autism spectrum disorder and fetal structural anomalies

Short-read genome sequencing (GS) holds the promise of becoming the primary diagnostic approach for the assessment of autism spectrum disorder (ASD) and fetal structural anomalies (FSAs). However, few studies have comprehensively evaluated its performance against current standard-of-care diagnostic tests: karyotype, chromosomal microarray (CMA), and exome sequencing (ES). To assess the clinical utility of GS, we compared its diagnostic yield against these three tests in 1,612 quartet families including an individual with ASD and in 295 prenatal families. Our GS analytic framework identified a diagnostic variant in 7.8% of ASD probands, almost 2-fold more than CMA (4.3%) and 3-fold more than ES (2.7%). However, when we systematically captured copy-number variants (CNVs) from the exome data, the diagnostic yield of ES (7.4%) was brought much closer to, but did not surpass, GS. Similarly, we estimated that GS could achieve an overall diagnostic yield of 46.1% in unselected FSAs, representing a 17.2% increased yield over karyotype, 14.1% over CMA, and 4.1% over ES with CNV calling or 36.1% increase without CNV discovery. Overall, GS provided an added diagnostic yield of 0.4% and 0.8% beyond the combination of all three standard-of-care tests in ASD and FSAs, respectively. This corresponded to nine GS unique diagnostic variants, including sequence variants in exons not captured by ES, structural variants (SVs) inaccessible to existing standard-of-care tests, and SVs where the resolution of GS changed variant classification. Overall, this large-scale evaluation demonstrated that GS significantly outperforms each individual standard-of-care test while also outperforming the combination of all three tests, thus warranting consideration as the first-tier diagnostic approach for the assessment of ASD and FSAs.

The rate of de novo structural variation is increased in in vitro-produced offspring and preferentially affects the paternal genome

Assisted reproductive technologies (ARTs), including in vitro maturation and fertilization (IVF), are increasingly used in human and animal reproduction. Whether these technologies directly affect the rate of de novo mutation (DNM), and to what extent, has been a matter of debate. Here we take advantage of domestic cattle, characterized by complex pedigrees that are ideally suited to detect DNMs and by the systematic use of ART, to study the rate of de novo structural variation (dnSV) in this species and how it is impacted by IVF. By exploiting features of associated de novo point mutations (dnPMs) and dnSVs in clustered DNMs, we provide strong evidence that (1) IVF increases the rate of dnSV approximately fivefold, and (2) the corresponding mutations occur during the very early stages of embryonic development (one- and two-cell stage), yet primarily affect the paternal genome.

GATK-gCNV enables the discovery of rare copy number variants from exome sequencing data

Copy number variants (CNVs) are major contributors to genetic diversity and disease. While standardized methods, such as the genome analysis toolkit (GATK), exist for detecting short variants, technical challenges have confounded uniform large-scale CNV analyses from whole-exome sequencing (WES) data. Given the profound impact of rare and de novo coding CNVs on genome organization and human disease, we developed GATK-gCNV, a flexible algorithm to discover rare CNVs from sequencing read-depth information, complete with open-source distribution via GATK. We benchmarked GATK-gCNV in 7,962 exomes from individuals in quartet families with matched genome sequencing and microarray data, finding up to 95% recall of rare coding CNVs at a resolution of more than two exons. We used GATK-gCNV to generate a reference catalog of rare coding CNVs in WES data from 197,306 individuals in the UK Biobank, and observed strong correlations between per-gene CNV rates and measures of mutational constraint, as well as rare CNV associations with multiple traits. In summary, GATK-gCNV is a tunable approach for sensitive and specific CNV discovery in WES data, with broad applications.

The contributions of rare inherited and polygenic risk to ASD in multiplex families

Autism spectrum disorder (ASD) has a complex genetic architecture involving contributions from both de novo and inherited variation. Few studies have been designed to address the role of rare inherited variation or its interaction with common polygenic risk in ASD. Here, we performed whole-genome sequencing of the largest cohort of multiplex families to date, consisting of 4,551 individuals in 1,004 families having two or more autistic children. Using this study design, we identify seven previously unrecognized ASD risk genes supported by a majority of rare inherited variants, finding support for a total of 74 genes in our cohort and a total of 152 genes after combined analysis with other studies. Autistic children from multiplex families demonstrate an increased burden of rare inherited protein-truncating variants in known ASD risk genes. We also find that ASD polygenic score (PGS) is overtransmitted from nonautistic parents to autistic children who also harbor rare inherited variants, consistent with combinatorial effects in the offspring, which may explain the reduced penetrance of these rare variants in parents. We also observe that in addition to social dysfunction, language delay is associated with ASD PGS overtransmission. These results are consistent with an additive complex genetic risk architecture of ASD involving rare and common variation and further suggest that language delay is a core biological feature of ASD.

MSV: a modular structural variant caller that reveals nested and complex rearrangements by unifying breakends inferred directly from reads

Structural variant (SV) calling belongs to the standard tools of modern bioinformatics for identifying and describing alterations in genomes. Initially, this work presents several complex genomic rearrangements that reveal conceptual ambiguities inherent to the representation via basic SV. We contextualize these ambiguities theoretically as well as practically and propose a graph-based approach for resolving them. For various yeast genomes, we practically compute adjacency matrices of our graph model and demonstrate that they provide highly accurate descriptions of one genome in terms of another. An open-source prototype implementation of our approach is available under the MIT license at https://github.com/ITBE-Lab/MA .

Challenges in screening for de novo noncoding variants contributing to genetically complex phenotypes

Understanding the genetic basis for complex, heterogeneous disorders, such as autism spectrum disorder (ASD), is a persistent challenge in human medicine. Owing to their phenotypic complexity, the genetic mechanisms underlying these disorders may be highly variable across individual patients. Furthermore, much of their heritability is unexplained by known regulatory or coding variants. Indeed, there is evidence that much of the causal genetic variation stems from rare and de novo variants arising from ongoing mutation. These variants occur mostly in noncoding regions, likely affecting regulatory processes for genes linked to the phenotype of interest. However, because there is no uniform code for assessing regulatory function, it is difficult to separate these mutations into likely functional and nonfunctional subsets. This makes finding associations between complex diseases and potentially causal de novo single-nucleotide variants (dnSNVs) a difficult task. To date, most published studies have struggled to find any significant associations between dnSNVs from ASD patients and any class of known regulatory elements. We sought to identify the underlying reasons for this and present strategies for overcoming these challenges. We show that, contrary to previous claims, the main reason for failure to find robust statistical enrichments is not only the number of families sampled, but also the quality and relevance to ASD of the annotations used to prioritize dnSNVs, and the reliability of the set of dnSNVs itself. We present a list of recommendations for designing future studies of this sort that will help researchers avoid common pitfalls.

Detection of trait-associated structural variations using short-read sequencing

Genomic structural variation (SV) affects genetic and phenotypic characteristics in diverse organisms, but the lack of reliable methods to detect SV has hindered genetic analysis. We developed a computational algorithm (MOPline) that includes missing call recovery combined with high-confidence SV call selection and genotyping using short-read whole-genome sequencing (WGS) data. Using 3,672 high-coverage WGS datasets, MOPline stably detected ∼16,000 SVs per individual, which is over ∼1.7-3.3-fold higher than previous large-scale projects while exhibiting a comparable level of statistical quality metrics. We imputed SVs from 181,622 Japanese individuals for 42 diseases and 60 quantitative traits. A genome-wide association study with the imputed SVs revealed 41 top-ranked or nearly top-ranked genome-wide significant SVs, including 8 exonic SVs with 5 novel associations and enriched mobile element insertions. This study demonstrates that short-read WGS data can be used to identify rare and common SVs associated with a variety of traits.

Genome-wide structural variant analysis identifies risk loci for non-Alzheimer's dementias

We characterized the role of structural variants, a largely unexplored type of genetic variation, in two non-Alzheimer's dementias, namely Lewy body dementia (LBD) and frontotemporal dementia (FTD)/amyotrophic lateral sclerosis (ALS). To do this, we applied an advanced structural variant calling pipeline (GATK-SV) to short-read whole-genome sequence data from 5,213 European-ancestry cases and 4,132 controls. We discovered, replicated, and validated a deletion in TPCN1 as a novel risk locus for LBD and detected the known structural variants at the C9orf72 and MAPT loci as associated with FTD/ALS. We also identified rare pathogenic structural variants in both LBD and FTD/ALS. Finally, we assembled a catalog of structural variants that can be mined for new insights into the pathogenesis of these understudied forms of dementia.

Novel correlative analysis identifies multiple genomic variations impacting ASD with macrocephaly

Autism spectrum disorders (ASD) display both phenotypic and genetic heterogeneity, impeding the understanding of ASD and development of effective means of diagnosis and potential treatments. Genes affected by genomic variations for ASD converge in dozens of gene ontologies (GOs), but the relationship between the variations at the GO level have not been well elucidated. In the current study, multiple types of genomic variations were mapped to GOs and correlations among GOs were measured in ASD and control samples. Several ASD-unique GO correlations were found, suggesting the importance of co-occurrence of genomic variations in genes from different functional categories in ASD etiology. Combined with experimental data, several variations related to WNT signaling, neuron development, synapse morphology/function and organ morphogenesis were found to be important for ASD with macrocephaly, and novel co-occurrence patterns of them in ASD patients were found. Furthermore, we applied this gene ontology correlation analysis method to find genomic variations that contribute to ASD etiology in combination with changes in gene expression and transcription factor binding, providing novel insights into ASD with macrocephaly and a new methodology for the analysis of genomic variation.

Regulatory de novo mutations underlying intellectual disability

The genetic aetiology of a major fraction of patients with intellectual disability (ID) remains unknown. De novo mutations (DNMs) in protein-coding genes explain up to 40% of cases, but the potential role of regulatory DNMs is still poorly understood. We sequenced 63 whole genomes from 21 ID probands and their unaffected parents. In addition, we analysed 30 previously sequenced genomes from exome-negative ID probands. We found that regulatory DNMs were selectively enriched in fetal brain-specific enhancers as compared with adult brain enhancers. DNM-containing enhancers were associated with genes that show preferential expression in the prefrontal cortex. Furthermore, we identified recurrently mutated enhancer clusters that regulate genes involved in nervous system development (CSMD1, OLFM1, and POU3F3). Most of the DNMs from ID probands showed allele-specific enhancer activity when tested using luciferase assay. Using CRISPR-mediated mutation and editing of epigenomic marks, we show that DNMs at regulatory elements affect the expression of putative target genes. Our results, therefore, provide new evidence to indicate that DNMs in fetal brain-specific enhancers play an essential role in the aetiology of ID.

Genome-Wide Analysis of Structural Variants in Parkinson Disease

OBJECTIVE

Identification of genetic risk factors for Parkinson disease (PD) has to date been primarily limited to the study of single nucleotide variants, which only represent a small fraction of the genetic variation in the human genome. Consequently, causal variants for most PD risk are not known. Here we focused on structural variants (SVs), which represent a major source of genetic variation in the human genome. We aimed to discover SVs associated with PD risk by performing the first large-scale characterization of SVs in PD.

METHODS

We leveraged a recently developed computational pipeline to detect and genotype SVs from 7,772 Illumina short-read whole genome sequencing samples. Using this set of SV variants, we performed a genome-wide association study using 2,585 cases and 2,779 controls and identified SVs associated with PD risk. Furthermore, to validate the presence of these variants, we generated a subset of matched whole-genome long-read sequencing data.

RESULTS

We genotyped and tested 3,154 common SVs, representing over 412 million nucleotides of previously uncatalogued genetic variation. Using long-read sequencing data, we validated the presence of three novel deletion SVs that are associated with risk of PD from our initial association analysis, including a 2 kb intronic deletion within the gene LRRN4.

INTERPRETATION

We identified three SVs associated with genetic risk of PD. This study represents the most comprehensive assessment of the contribution of SVs to the genetic risk of PD to date. ANN NEUROL 2023;93:1012-1022.

Chromosomal Instability in Genome Evolution: From Cancer to Macroevolution

The integrity of the genome is crucial for the survival of all living organisms. However, genomes need to adapt to survive certain pressures, and for this purpose use several mechanisms to diversify. Chromosomal instability (CIN) is one of the main mechanisms leading to the creation of genomic heterogeneity by altering the number of chromosomes and changing their structures. In this review, we will discuss the different chromosomal patterns and changes observed in speciation, in evolutional biology as well as during tumor progression. By nature, the human genome shows an induction of diversity during gametogenesis but as well during tumorigenesis that can conclude in drastic changes such as the whole genome doubling to more discrete changes as the complex chromosomal rearrangement chromothripsis. More importantly, changes observed during speciation are strikingly similar to the genomic evolution observed during tumor progression and resistance to therapy. The different origins of CIN will be treated as the importance of double-strand breaks (DSBs) or the consequences of micronuclei. We will also explain the mechanisms behind the controlled DSBs, and recombination of homologous chromosomes observed during meiosis, to explain how errors lead to similar patterns observed during tumorigenesis. Then, we will also list several diseases associated with CIN, resulting in fertility issues, miscarriage, rare genetic diseases, and cancer. Understanding better chromosomal instability as a whole is primordial for the understanding of mechanisms leading to tumor progression.

A collection of read depth profiles at structural variant breakpoints

SWaveform, a newly created open genome-wide resource for read depth signal in the vicinity of structural variant (SV) breakpoints, aims to boost development of computational tools and algorithms for discovery of genomic rearrangement events from sequencing data. SVs are a dominant force shaping genomes and substantially contributing to genetic diversity. Still, there are challenges in reliable and efficient genotyping of SVs from whole genome sequencing data, thus delaying translation into clinical applications and wasting valuable resources. SWaveform includes a database containing ~7 M of read depth profiles at SV breakpoints extracted from 911 sequencing samples generated by the Human Genome Diversity Project, generalised patterns of the signal at breakpoints, an interface for navigation and download, as well as a toolbox for local deployment with user's data. The dataset can be of immense value to bioinformatics and engineering communities as it empowers smooth application of intelligent signal processing and machine learning techniques for discovery of genomic rearrangement events and thus opens the floodgates for development of innovative algorithms and software.

Whole Genome Analysis of Dizygotic Twins With Autism Reveals Prevalent Transposon Insertion Within Neuronal Regulatory Elements: Potential Implications for Disease Etiology and Clinical Assessment

Transposable elements (TEs) have been implicated in autism spectrum disorder (ASD). However, our understanding of their roles is far from complete. Herein, we explored de novo TE insertions (dnTEIs) and de novo variants (DNVs) across the genomes of dizygotic twins with ASD and their parents. The neuronal regulatory elements had a tendency to harbor dnTEIs that were shared between twins, but ASD-risk genes had dnTEIs that were unique to each twin. The dnTEIs were 4.6-fold enriched in enhancers that are active in embryonic stem cell (ESC)-neurons (p < 0.001), but DNVs were 1.5-fold enriched in active enhancers of astrocytes (p = 0.0051). Our findings suggest that dnTEIs and DNVs play a role in ASD etiology by disrupting enhancers of neurons and astrocytes.

Single-neuron whole genome sequencing identifies increased somatic mutation burden in Alzheimer's disease related genes

Accumulation of somatic mutations in human neurons is associated with aging and neurodegeneration. To shed light on the somatic mutational burden in Alzheimer's disease (AD) neurons and get more insight into the role of somatic mutations in AD pathogenesis, we performed single-neuron whole genome sequencing to detect genome-wide somatic mutations (single nucleotide variants (SNVs) and Indels) in 96 single prefrontal cortex neurons from 8 AD patients and 8 elderly controls. We found that the mutational burden is ∼3000 somatic mutations per neuron genome in elderly subjects. AD patients have increased somatic mutation burden in AD-related annotation categories, including AD risk genes and differentially expressed genes in AD neurons. Mutational signature analysis showed somatic SNVs (sSNVs) primarily caused by aging and oxidative DNA damage processes but no significant difference was detected between AD and controls. Additionally, functional somatic mutations identified in AD patients showed significant enrichment in several AD-related pathways, including AD pathway, Notch-signaling pathway and Calcium-signaling pathway. These findings provide genetic insights into how somatic mutations may alter the function of single neurons and exert their potential roles in the pathogenesis of AD.

ParseCNV2: efficient sequencing tool for copy number variation genome-wide association studies

Improved copy number variation (CNV) detection remains an area of heavy emphasis for algorithm development; however, both CNV curation and disease association approaches remain in its infancy. The current practice of focusing on candidate CNVs, where researchers study specific CNVs they believe to be pathological while discarding others, refrains from considering the full spectrum of CNVs in a hypothesis-free GWAS. To address this, we present a next-generation approach to CNV association by natively supporting the popular VCF specification for sequencing-derived variants as well as SNP array calls using a PennCNV format. The code is fast and efficient, allowing for the analysis of large (>100,000 sample) cohorts without dividing up the data on a compute cluster. The scripts are condensed into a single tool to promote simplicity and best practices. CNV curation pre and post-association is rigorously supported and emphasized to yield reliable results of highest quality. We benchmarked two large datasets, including the UK Biobank (n > 450,000) and CAG Biobank (n > 350,000) both of which are genotyped at >0.5 M probes, for our input files. ParseCNV has been actively supported and developed since 2008. ParseCNV2 presents a critical addition to formalizing CNV association for inclusion with SNP associations in GWAS Catalog. Clinical CNV prioritization, interactive quality control (QC), and adjustment for covariates are revolutionary new features of ParseCNV2 vs. ParseCNV. The software is freely available at: https://github.com/CAG-CNV/ParseCNV2 .

Characterization of De Novo Promoter Variants in Autism Spectrum Disorder with Massively Parallel Reporter Assays

Autism spectrum disorder (ASD) is a common, complex, and highly heritable condition with contributions from both common and rare genetic variations. While disruptive, rare variants in protein-coding regions clearly contribute to symptoms, the role of rare non-coding remains unclear. Variants in these regions, including promoters, can alter downstream RNA and protein quantity; however, the functional impacts of specific variants observed in ASD cohorts remain largely uncharacterized. Here, we analyzed 3600 de novo mutations in promoter regions previously identified by whole-genome sequencing of autistic probands and neurotypical siblings to test the hypothesis that mutations in cases have a greater functional impact than those in controls. We leveraged massively parallel reporter assays (MPRAs) to detect transcriptional consequences of these variants in neural progenitor cells and identified 165 functionally high confidence de novo variants (HcDNVs). While these HcDNVs are enriched for markers of active transcription, disruption to transcription factor binding sites, and open chromatin, we did not identify differences in functional impact based on ASD diagnostic status.