Genome-wide study of gene-by-sex interactions identifies risks for cleft palate.

Structural birth defects affect 3-4% of all live births and, depending on the type, tend to manifest in a sex-biased manner. Orofacial clefts (OFCs) are the most common craniofacial structural birth defects and are often divided into cleft lip with or without cleft palate (CL/P) and cleft palate only (CP). Previous studies have found sex-specific risks for CL/P, but these risks have yet to be evaluated in CP. CL/P is more common in males and CP is more frequently observed in females, so we hypothesized there would also be sex-specific differences for CP. Using a trio-based cohort, we performed sex-stratified genome-wide association studies (GWAS) based on proband sex followed by a genome-wide gene-by-sex (GxS) interaction testing. There were 13 loci significant for GxS interactions, with the top finding in LTBP1 (RR=3.37 [2.04 - 5.56], p=1.93x10 -6 ). LTBP1 plays a role in regulating TGF-B bioavailability, and knockdown in both mice and zebrafish lead to craniofacial anomalies. Further, there is evidence for differential expression of LTBP1 between males and females in both mice and humans. Therefore, we tested the association between the imputed genetically regulated gene expression of genes with significant GxS interactions and the CP phenotype. We found significant association for LTBP1 in cell cultured fibroblasts in female probands (p=0.0013) but not in males. Taken altogether, we show there are sex-specific risks for CP that are otherwise undetectable in a combined sex cohort, and LTBP1 is a candidate risk gene, particularly in females.

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.

A harmonized public resource of deeply sequenced diverse human genomes.

Underrepresented populations are often excluded from genomic studies due in part to a lack of resources supporting their analyses. The 1000 Genomes Project (1kGP) and Human Genome Diversity Project (HGDP), which have recently been sequenced to high coverage, are valuable genomic resources because of the global diversity they capture and their open data sharing policies. Here, we harmonized a high quality set of 4,094 whole genomes from 80 populations in the HGDP and 1kGP with data from the Genome Aggregation Database (gnomAD) and identified over 153 million high-quality SNVs, indels, and SVs. We performed a detailed ancestry analysis of this cohort, characterizing population structure and patterns of admixture across populations, analyzing site frequency spectra, and measuring variant counts at global and subcontinental levels. We also demonstrate substantial added value from this dataset compared to the prior versions of the component resources, typically combined via liftOver and variant intersection; for example, we catalog millions of new genetic variants, mostly rare, compared to previous releases. In addition to unrestricted individual-level public release, we provide detailed tutorials for conducting many of the most common quality control steps and analyses with these data in a scalable cloud-computing environment and publicly release this new phased joint callset for use as a haplotype resource in phasing and imputation pipelines. This jointly called reference panel will serve as a key resource to support research of diverse ancestry populations.

Expanding the genetics and phenotypes of ocular congenital cranial dysinnervation disorders.

Purpose: To identify genetic etiologies and genotype/phenotype associations for unsolved ocular congenital cranial dysinnervation disorders (oCCDDs). Methods: We coupled phenotyping with exome or genome sequencing of 467 pedigrees with genetically unsolved oCCDDs, integrating analyses of pedigrees, human and animal model phenotypes, and de novo variants to identify rare candidate single nucleotide variants, insertion/deletions, and structural variants disrupting protein-coding regions. Prioritized variants were classified for pathogenicity and evaluated for genotype/phenotype correlations. Results: Analyses elucidated phenotypic subgroups, identified pathogenic/likely pathogenic variant(s) in 43/467 probands (9.2%), and prioritized variants of uncertain significance in 70/467 additional probands (15.0%). These included known and novel variants in established oCCDD genes, genes associated with syndromes that sometimes include oCCDDs (e.g., MYH10, KIF21B, TGFBR2, TUBB6), genes that fit the syndromic component of the phenotype but had no prior oCCDD association (e.g., CDK13, TGFB2), genes with no reported association with oCCDDs or the syndromic phenotypes (e.g., TUBA4A, KIF5C, CTNNA1, KLB, FGF21), and genes associated with oCCDD phenocopies that had resulted in misdiagnoses. Conclusion: This study suggests that unsolved oCCDDs are clinically and genetically heterogeneous disorders often overlapping other Mendelian conditions and nominates many candidates for future replication and functional studies.

Narrowing the diagnostic gap: Genomes, episignatures, long-read sequencing, and health economic analyses in an exome-negative intellectual disability cohort.

PURPOSE: Genome sequencing (GS)-specific diagnostic rates in prospective tightly ascertained exome sequencing (ES)-negative intellectual disability (ID) cohorts have not been reported extensively. METHODS: ES, GS, epigenetic signatures, and long-read sequencing diagnoses were assessed in 74 trios with at least moderate ID. RESULTS: The ES diagnostic yield was 42 of 74 (57%). GS diagnoses were made in 9 of 32 (28%) ES-unresolved families. Repeated ES with a contemporary pipeline on the GS-diagnosed families identified 8 of 9 single-nucleotide variations/copy-number variations undetected in older ES, confirming a GS-unique diagnostic rate of 1 in 32 (3%). Episignatures contributed diagnostic information in 9% with GS corroboration in 1 of 32 (3%) and diagnostic clues in 2 of 32 (6%). A genetic etiology for ID was detected in 51 of 74 (69%) families. Twelve candidate disease genes were identified. Contemporary ES followed by GS cost US$4976 (95% CI: $3704; $6969) per diagnosis and first-line GS at a cost of $7062 (95% CI: $6210; $8475) per diagnosis. CONCLUSION: Performing GS only in ID trios would be cost equivalent to ES if GS were available at $2435, about a 60% reduction from current prices. This study demonstrates that first-line GS achieves higher diagnostic rate than contemporary ES but at a higher cost.

Improving prenatal diagnosis through standards and aggregation.

Advances in sequencing and imaging technologies enable enhanced assessment in the prenatal space, with a goal to diagnose and predict the natural history of disease, to direct targeted therapies, and to implement clinical management, including transfer of care, election of supportive care, and selection of surgical interventions. The current lack of standardization and aggregation stymies variant interpretation and gene discovery, which hinders the provision of prenatal precision medicine, leaving clinicians and patients without an accurate diagnosis. With large amounts of data generated, it is imperative to establish standards for data collection, processing, and aggregation. Aggregated and homogeneously processed genetic and phenotypic data permits dissection of the genomic architecture of prenatal presentations of disease and provides a dataset on which data analysis algorithms can be tuned to the prenatal space. Here we discuss the importance of generating aggregate data sets and how the prenatal space is driving the development of interoperable standards and phenotype-driven tools.

A harmonized public resource of deeply sequenced diverse human genomes.

Underrepresented populations are often excluded from genomic studies due in part to a lack of resources supporting their analyses. The 1000 Genomes Project (1kGP) and Human Genome Diversity Project (HGDP), which have recently been sequenced to high coverage, are valuable genomic resources because of the global diversity they capture and their open data sharing policies. Here, we harmonized a high quality set of 4,094 whole genomes from HGDP and 1kGP with data from the Genome Aggregation Database (gnomAD) and identified over 153 million high-quality SNVs, indels, and SVs. We performed a detailed ancestry analysis of this cohort, characterizing population structure and patterns of admixture across populations, analyzing site frequency spectra, and measuring variant counts at global and subcontinental levels. We also demonstrate substantial added value from this dataset compared to the prior versions of the component resources, typically combined via liftover and variant intersection; for example, we catalog millions of new genetic variants, mostly rare, compared to previous releases. In addition to unrestricted individual-level public release, we provide detailed tutorials for conducting many of the most common quality control steps and analyses with these data in a scalable cloud-computing environment and publicly release this new phased joint callset for use as a haplotype resource in phasing and imputation pipelines. This jointly called reference panel will serve as a key resource to support research of diverse ancestry populations.

PLS3 missense variants affecting the actin-binding domains cause X-linked congenital diaphragmatic hernia and body-wall defects.

Congenital diaphragmatic hernia (CDH) is a relatively common and genetically heterogeneous structural birth defect associated with high mortality and morbidity. We describe eight unrelated families with an X-linked condition characterized by diaphragm defects, variable anterior body-wall anomalies, and/or facial dysmorphism. Using linkage analysis and exome or genome sequencing, we found that missense variants in plastin 3 (PLS3), a gene encoding an actin bundling protein, co-segregate with disease in all families. Loss-of-function variants in PLS3 have been previously associated with X-linked osteoporosis (MIM: 300910), so we used in silico protein modeling and a mouse model to address these seemingly disparate clinical phenotypes. The missense variants in individuals with CDH are located within the actin-binding domains of the protein but are not predicted to affect protein structure, whereas the variants in individuals with osteoporosis are predicted to result in loss of function. A mouse knockin model of a variant identified in one of the CDH-affected families, c.1497G>C (p.Trp499Cys), shows partial perinatal lethality and recapitulates the key findings of the human phenotype, including diaphragm and abdominal-wall defects. Both the mouse model and one adult human male with a CDH-associated PLS3 variant were observed to have increased rather than decreased bone mineral density. Together, these clinical and functional data in humans and mice reveal that specific missense variants affecting the actin-binding domains of PLS3 might have a gain-of-function effect and cause a Mendelian congenital disorder.

Trio-based GWAS identifies novel associations and subtype-specific risk factors for cleft palate.

Cleft palate (CP) is one of the most common craniofacial birth defects; however, there are relatively few established genetic risk factors associated with its occurrence despite high heritability. Historically, CP has been studied as a single phenotype, although it manifests across a spectrum of defects involving the hard and/or soft palate. We performed a genome-wide association study using transmission disequilibrium tests of 435 case-parent trios to evaluate broad risks for any cleft palate (ACP) (n = 435), and subtype-specific risks for any cleft soft palate (CSP), (n = 259) and any cleft hard palate (CHP) (n = 125). We identified a single genome-wide significant locus at 9q33.3 (lead SNP rs7035976, p = 4.24 × 10-8) associated with CHP. One gene at this locus, angiopoietin-like 2 (ANGPTL2), plays a role in osteoblast differentiation. It is expressed both in craniofacial tissue of human embryos and developing mouse palatal shelves. We found 19 additional loci reaching suggestive significance (p < 5 × 10-6), of which only one overlapped between groups (chromosome 17q24.2, ACP and CSP). Odds ratios for the 20 loci were most similar across all 3 groups for SNPs associated with the ACP group, but more distinct when comparing SNPs associated with either subtype. We also found nominal evidence of replication (p < 0.05) for 22 SNPs previously associated with orofacial clefts. Our study to evaluate CP risks in the context of its subtypes and we provide newly reported associations affecting the broad risk for CP as well as evidence of subtype-specific risks.

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.

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.

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.

Phenotype and genetic analysis of data collected within the first year of NeuroDev.

Genetic association studies have made significant contributions to our understanding of the etiology of neurodevelopmental disorders (NDDs). However, these studies rarely focused on the African continent. The NeuroDev Project aims to address this diversity gap through detailed phenotypic and genetic characterization of children with NDDs from Kenya and South Africa. We present results from NeuroDev's first year of data collection, including phenotype data from 206 cases and clinical genetic analyses of 99 parent-child trios. Most cases met criteria for global developmental delay/intellectual disability (GDD/ID, 80.3%). Approximately half of the children with GDD/ID also met criteria for autism. Analysis of exome-sequencing data identified a pathogenic or likely pathogenic variant in 13 (17%) of the 75 cases from South Africa and 9 (38%) of the 24 cases from Kenya. Data from the trio pilot are publicly available, and the NeuroDev Project will continue to develop resources for the global genetics community.

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.

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.

Trio-based GWAS identifies novel associations and subtype-specific risk factors for cleft palate.

Orofacial clefts (OFCs) are the most common craniofacial birth defects and are often categorized into two etiologically distinct groups: cleft lip with or without cleft palate (CL/P) and isolated cleft palate (CP). CP is highly heritable, but there are still relatively few established genetic risk factors associated with its occurrence compared to CL/P. Historically, CP has been studied as a single phenotype despite manifesting across a spectrum of defects involving the hard and/or soft palate. We performed GWAS using transmission disequilibrium tests using 435 case-parent trios to evaluate broad risks for any cleft palate (ACP, n=435), as well as subtype-specific risks for any cleft soft palate (CSP, n=259) and any cleft hard palate (CHP, n=125). We identified a single genome-wide significant locus at 9q33.3 (lead SNP rs7035976, p=4.24×10 -8 ) associated with CHP. One gene at this locus, angiopoietin-like 2 ( ANGPTL2 ), plays a role in osteoblast differentiation. It is expressed in craniofacial tissue of human embryos, as well as in the developing mouse palatal shelves. We found 19 additional loci reaching suggestive significance (p<5×10 -6 ), of which only one overlapped between groups (chromosome 17q24.2, ACP and CSP). Odds ratios (ORs) for each of the 20 loci were most similar across all three groups for SNPs associated with the ACP group, but more distinct when comparing SNPs associated with either the CSP or CHP groups. We also found nominal evidence of replication (p<0.05) for 22 SNPs previously associated with cleft palate (including CL/P). Interestingly, most SNPs associated with CL/P cases were found to convey the opposite effect in those replicated in our dataset for CP only. Ours is the first study to evaluate CP risks in the context of its subtypes and we provide newly reported associations affecting the broad risk for CP as well as evidence of subtype-specific risks.

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.

A cell type-aware framework for nominating non-coding variants in Mendelian regulatory disorders.

Unsolved Mendelian cases often lack obvious pathogenic coding variants, suggesting potential non-coding etiologies. Here, we present a single cell multi-omic framework integrating embryonic mouse chromatin accessibility, histone modification, and gene expression assays to discover cranial motor neuron (cMN) cis-regulatory elements and subsequently nominate candidate non-coding variants in the congenital cranial dysinnervation disorders (CCDDs), a set of Mendelian disorders altering cMN development. We generated single cell epigenomic profiles for ~86,000 cMNs and related cell types, identifying ~250,000 accessible regulatory elements with cognate gene predictions for ~145,000 putative enhancers. Seventy-five percent of elements (44 of 59) validated in an in vivo transgenic reporter assay, demonstrating that single cell accessibility is a strong predictor of enhancer activity. Applying our cMN atlas to 899 whole genome sequences from 270 genetically unsolved CCDD pedigrees, we achieved significant reduction in our variant search space and nominated candidate variants predicted to regulate known CCDD disease genes MAFB, PHOX2A, CHN1, and EBF3 - as well as new candidates in recurrently mutated enhancers through peak- and gene-centric allelic aggregation. This work provides novel non-coding variant discoveries of relevance to CCDDs and a generalizable framework for nominating non-coding variants of potentially high functional impact in other Mendelian disorders.