Genomic architecture of autism from comprehensive whole-genome sequence annotation.

Fully understanding autism spectrum disorder (ASD) genetics requires whole-genome sequencing (WGS). We present the latest release of the Autism Speaks MSSNG resource, which includes WGS data from 5,100 individuals with ASD and 6,212 non-ASD parents and siblings (total n = 11,312). Examining a wide variety of genetic variants in MSSNG and the Simons Simplex Collection (SSC; n = 9,205), we identified ASD-associated rare variants in 718/5,100 individuals with ASD from MSSNG (14.1%) and 350/2,419 from SSC (14.5%). Considering genomic architecture, 52% were nuclear sequence-level variants, 46% were nuclear structural variants (including copy-number variants, inversions, large insertions, uniparental isodisomies, and tandem repeat expansions), and 2% were mitochondrial variants. Our study provides a guidebook for exploring genotype-phenotype correlations in families who carry ASD-associated rare variants and serves as an entry point to the expanded studies required to dissect the etiology in the ∼85% of the ASD population that remain idiopathic.

CLP1 founder mutation links tRNA splicing and maturation to cerebellar development and neurodegeneration.

Neurodegenerative diseases can occur so early as to affect neurodevelopment. From a cohort of more than 2,000 consanguineous families with childhood neurological disease, we identified a founder mutation in four independent pedigrees in cleavage and polyadenylation factor I subunit 1 (CLP1). CLP1 is a multifunctional kinase implicated in tRNA, mRNA, and siRNA maturation. Kinase activity of the CLP1 mutant protein was defective, and the tRNA endonuclease complex (TSEN) was destabilized, resulting in impaired pre-tRNA cleavage. Germline clp1 null zebrafish showed cerebellar neurodegeneration that was rescued by wild-type, but not mutant, human CLP1 expression. Patient-derived induced neurons displayed both depletion of mature tRNAs and accumulation of unspliced pre-tRNAs. Transfection of partially processed tRNA fragments into patient cells exacerbated an oxidative stress-induced reduction in cell survival. Our data link tRNA maturation to neuronal development and neurodegeneration through defective CLP1 function in humans.

SCIP: software for efficient clinical interpretation of copy number variants detected by whole-genome sequencing.

Copy number variants (CNVs) represent major etiologic factors in rare genetic diseases. Current clinical CNV interpretation workflows require extensive back-and-forth with multiple tools and databases. This increases complexity and time burden, potentially resulting in missed genetic diagnoses. We present the Suite for CNV Interpretation and Prioritization (SCIP), a software package for the clinical interpretation of CNVs detected by whole-genome sequencing (WGS). The SCIP Visualization Module near-instantaneously displays all information necessary for CNV interpretation (variant quality, population frequency, inheritance pattern, and clinical relevance) on a single page-supported by modules providing variant filtration and prioritization. SCIP was comprehensively evaluated using WGS data from 1027 families with congenital cardiac disease and/or autism spectrum disorder, containing 187 pathogenic or likely pathogenic (P/LP) CNVs identified in previous curations. SCIP was efficient in filtration and prioritization: a median of just two CNVs per case were selected for review, yet it captured all P/LP findings (92.5% of which ranked 1st). SCIP was also able to identify one pathogenic CNV previously missed. SCIP was benchmarked against AnnotSV and a spreadsheet-based manual workflow and performed superiorly than both. In conclusion, SCIP is a novel software package for efficient clinical CNV interpretation, substantially faster and more accurate than previous tools (available at https://github.com/qd29/SCIP , a video tutorial series is available at https://bit.ly/SCIPVideos ).

Pharmacogenetic profiling via genome sequencing in children with medical complexity.

BACKGROUND: Children with medical complexity (CMC) are a priority pediatric population, with high resource use and associated costs. Genome-wide sequencing is increasingly organized for CMC early in life as a diagnostic test. Polypharmacy becomes common as CMC age. Clinically relevant pharmacogenetic (PGx) information can be extracted from existing genome sequencing (GS) data via GS-PGx profiling. The role of GS-PGx profiling in the CMC population is unclear. METHODS: Prescribed medications were extracted from care plans of 802 eligible CMC enrolled in a structured Complex Care Program over a 10-year period. Drug-gene associations were annotated using curated Clinical Pharmacogenetics Implementation Consortium data. GS-PGx profiling was then performed for a subset of 50 CMC. RESULTS: Overall, 546 CMC (68%) were prescribed at least one medication with an established PGx association. In the GS-PGx subgroup, 24 (48%) carried variants in pharmacogenes with drug-gene guidelines for one or more of their current medications. All had findings of potential relevance to some medications, including 32 (64%) with variants in CYP2C19 that could affect their metabolism of proton-pump inhibitors. CONCLUSION: GS-PGx profiling at the time of diagnostics-focused genetic testing could be an efficient way to incorporate precision prescribing practices into the lifelong care of CMC. IMPACT: Polypharmacy and genetic test utilization are both common in children with medical complexity. The role of repurposing genome sequencing data for pharmacogenetic profiling in children with medical complexity was previously unclear. We identified a high rate of medication use with clinically relevant drug-gene associations in this priority pediatric population and demonstrated that relevant pharmacogenetic information can be extracted from their existing genome sequencing data. Pharmacogenetic profiling at the time of diagnostics-focused genetic testing could be an efficient way to incorporate precision prescribing practices into the lifelong care of children with medical complexity.

Trio genome sequencing for developmental delay and pediatric heart conditions: A comparative microcost analysis.

PURPOSE: Genome sequencing (GS) can aid clinical management of multiple pediatric conditions. Insurers require accurate cost information to inform funding and implementation decisions. The objective was to compare the laboratory workflows and microcosts of trio GS testing in children with developmental delay (DD) and in children with cardiac conditions. METHODS: Cost items related to each step in trio GS (child and 2 parents) for both populations were identified and measured. Program costs over 5 years were estimated. Probabilistic and deterministic analyses were conducted. RESULTS: The mean cost per trio GS was CAD$6634.11 (95% CI = 6352.29-6913.40) for DD and CAD$8053.10 (95% CI = 7699.30-8558.10) for cardiac conditions. The 5-year program cost was CAD$28.11 million (95% CI = 26.91-29.29) for DD and CAD$5.63 million (95% CI = 5.38-5.98) for cardiac conditions. Supplies constituted the largest cost component for both populations. The higher cost per sample for the population with cardiac conditions was due to the inclusion of pharmacogenomics, higher bioinformatics labor costs, and a more labor intensive case review. CONCLUSION: This analysis indicated important variation in trio GS workflow and costs between pediatric populations in a single institution. Enhanced understanding of the clinical utility and costs of GS can inform harmonization and implementation decision-making.

Genome sequencing in families with congenital limb malformations.

The extensive clinical and genetic heterogeneity of congenital limb malformation calls for comprehensive genome-wide analysis of genetic variation. Genome sequencing (GS) has the potential to identify all genetic variants. Here we aim to determine the diagnostic potential of GS as a comprehensive one-test-for-all strategy in a cohort of undiagnosed patients with congenital limb malformations. We collected 69 cases (64 trios, 1 duo, 5 singletons) with congenital limb malformations with no molecular diagnosis after standard clinical genetic testing and performed genome sequencing. We also developed a framework to identify potential noncoding pathogenic variants. We identified likely pathogenic/disease-associated variants in 12 cases (17.4%) including four in known disease genes, and one repeat expansion in HOXD13. In three unrelated cases with ectrodactyly, we identified likely pathogenic variants in UBA2, establishing it as a novel disease gene. In addition, we found two complex structural variants (3%). We also identified likely causative variants in three novel high confidence candidate genes. We were not able to identify any noncoding variants. GS is a powerful strategy to identify all types of genomic variants associated with congenital limb malformation, including repeat expansions and complex structural variants missed by standard diagnostic approaches. In this cohort, no causative noncoding SNVs could be identified.

A Comprehensive Workflow for Read Depth-Based Identification of Copy-Number Variation from Whole-Genome Sequence Data.

A remaining hurdle to whole-genome sequencing (WGS) becoming a first-tier genetic test has been accurate detection of copy-number variations (CNVs). Here, we used several datasets to empirically develop a detailed workflow for identifying germline CNVs >1 kb from short-read WGS data using read depth-based algorithms. Our workflow is comprehensive in that it addresses all stages of the CNV-detection process, including DNA library preparation, sequencing, quality control, reference mapping, and computational CNV identification. We used our workflow to detect rare, genic CNVs in individuals with autism spectrum disorder (ASD), and 120/120 such CNVs tested using orthogonal methods were successfully confirmed. We also identified 71 putative genic de novo CNVs in this cohort, which had a confirmation rate of 70%; the remainder were incorrectly identified as de novo due to false positives in the proband (7%) or parental false negatives (23%). In individuals with an ASD diagnosis in which both microarray and WGS experiments were performed, our workflow detected all clinically relevant CNVs identified by microarrays, as well as additional potentially pathogenic CNVs < 20 kb. Thus, CNVs of clinical relevance can be discovered from WGS with a detection rate exceeding microarrays, positioning WGS as a single assay for genetic variation detection.

The Cardiac Genome Clinic: implementing genome sequencing in pediatric heart disease.

PURPOSE: This study investigated the diagnostic utility of nontargeted genomic testing in patients with pediatric heart disease. METHODS: We analyzed genome sequencing data of 111 families with cardiac lesions for rare, disease-associated variation. RESULTS: In 14 families (12.6%), we identified causative variants: seven were de novo (ANKRD11, KMT2D, NR2F2, POGZ, PTPN11, PURA, SALL1) and six were inherited from parents with no or subclinical heart phenotypes (FLT4, DNAH9, MYH11, NEXMIF, NIPBL, PTPN11). Outcome of the testing was associated with the presence of extracardiac features (p = 0.02), but not a positive family history for cardiac lesions (p = 0.67). We also report novel plausible gene-disease associations for tetralogy of Fallot/pulmonary stenosis (CDC42BPA, FGD5), hypoplastic left or right heart (SMARCC1, TLN2, TRPM4, VASP), congenitally corrected transposition of the great arteries (UBXN10), and early-onset cardiomyopathy (TPCN1). The identified candidate genes have critical functions in heart development, such as angiogenesis, mechanotransduction, regulation of heart size, chromatin remodeling, or ciliogenesis. CONCLUSION: This data set demonstrates the diagnostic and scientific value of genome sequencing in pediatric heart disease, anticipating its role as a first-tier diagnostic test. The genetic heterogeneity will necessitate large-scale genomic initiatives for delineating novel gene-disease associations.

7q31.2q31.31 deletion downstream of FOXP2 segregating in a family with speech and language disorder.

Chromosomal 7q31 deletions have been described in individuals with variable neurodevelopmental phenotypes including speech and language impairment. These copy number variants usually encompass FOXP2, haploinsufficiency of which represents a widely acknowledged cause for specific speech and language disorders. By chromosomal microarray analysis we identified a 4.7 Mb microdeletion at 7q31.2q31.31 downstream of FOXP2 in three family members presenting with variable speech, language and neurodevelopmental phenotypes. The index individual showed delayed speech development with impaired speech production, reduced language comprehension, and additionally learning difficulties, microcephaly, and attention deficit. His younger sister had delayed speech development with impaired speech production and partially reduced language comprehension. Their mother had attended a school for children with speech and language deficiencies and presented with impaired articulation. The deletion had occurred de novo in the mother, includes 15 protein-coding genes and is located in close proximity to the 3' end of FOXP2. Though a novel locus at 7q31.2q31.31 associated with mild neurodevelopmental and more prominent speech and language impairment is possible, the close phenotypic overlap with FOXP2-associated speech and language disorder rather suggests a positional effect on FOXP2 expression and function.

A novel intronic variant in UBE3A identified by genome sequencing in a patient with an atypical presentation of Angelman syndrome.

Angelman syndrome (AS) is a genetic neurodevelopmental disorder caused by loss or deficient expression of UBE3A on the maternally inherited allele. In 10-15% of individuals with a clinical diagnosis of AS, a molecular diagnosis cannot be established with conventional testing. We describe a 13-year-old male with an atypical presentation of AS, who was found to have a novel, maternally inherited, intronic variant in UBE3A (c.3-12T>A) using genome sequencing (GS). Targeted sequencing of RNA isolated from blood confirmed the creation of a new acceptor splice site. These GS results ended a six-year diagnostic odyssey and revealed a 50% recurrence risk for the unaffected parents. This case illustrates a previously unreported splicing variant causing AS. Intronic variants identifiable by GS may account for a proportion of individuals who are suspected of having well-known genetic disorders despite negative prior genetic testing.

Autosomal-Recessive Mutations in the tRNA Splicing Endonuclease Subunit TSEN15 Cause Pontocerebellar Hypoplasia and Progressive Microcephaly.

The tRNA splicing endonuclease is a highly evolutionarily conserved protein complex, involved in the cleavage of intron-containing tRNAs. In human it consists of the catalytic subunits TSEN2 and TSEN34, as well as the non-catalytic TSEN54 and TSEN15. Recessive mutations in the corresponding genes of the first three are known to cause pontocerebellar hypoplasia (PCH) types 2A-C, 4, and 5. Here, we report three homozygous TSEN15 variants that cause a milder version of PCH2. The affected individuals showed progressive microcephaly, delayed developmental milestones, intellectual disability, and, in two out of four cases, epilepsy. None, however, displayed the central visual failure seen in PCH case subjects where other subunits of the TSEN are mutated, and only one was affected by the extensive motor defects that are typical in other forms of PCH2. The three amino acid substitutions impacted the protein level of TSEN15 and the stoichiometry of the interacting subunits in different ways, but all resulted in an almost complete loss of in vitro tRNA cleavage activity. Taken together, our results demonstrate that mutations in any known subunit of the TSEN complex can cause PCH and progressive microcephaly, emphasizing the importance of its function during brain development.

Haploinsufficiency of vascular endothelial growth factor related signaling genes is associated with tetralogy of Fallot.

PURPOSE: To determine disease-associated single-gene variants in conotruncal defects, particularly tetralogy of Fallot (TOF). METHODS: We analyzed for rare loss-of-function and deleterious variants in FLT4 (VEGFR3) and other genes in the vascular endothelial growth factor (VEGF) pathway, as part of a genome sequencing study involving 175 adults with TOF from a single site. RESULTS: We identified nine (5.1%) probands with novel FLT4 variants: seven loss-of-function, including an 8-kb deletion, and two predicted damaging. In ten other probands we found likely disruptive variants in VEGF-related genes: KDR (VEGFR2; two stopgain and two nonsynonymous variants), VEGFA, FGD5, BCAR1, IQGAP1, FOXO1, and PRDM1. Detection of VEGF-related variants (19/175, 10.9%) was associated with an increased prevalence of absent pulmonary valve (26.3% vs. 3.4%, p < 0.0001) and right aortic arch (52.6% vs. 29.1%, p = 0.029). Extracardiac anomalies were rare. In an attempt to replicate findings, we identified three loss-of-function or damaging variants in FLT4, KDR, and IQGAP1 in ten independent families with TOF. CONCLUSION: Loss-of-function variants in FLT4 and KDR contribute substantially to the genetic basis of TOF. The findings support dysregulated VEGF signaling as a novel mechanism contributing to the pathogenesis of TOF.

Improved diagnostic yield compared with targeted gene sequencing panels suggests a role for whole-genome sequencing as a first-tier genetic test.

PurposeGenetic testing is an integral diagnostic component of pediatric medicine. Standard of care is often a time-consuming stepwise approach involving chromosomal microarray analysis and targeted gene sequencing panels, which can be costly and inconclusive. Whole-genome sequencing (WGS) provides a comprehensive testing platform that has the potential to streamline genetic assessments, but there are limited comparative data to guide its clinical use.MethodsWe prospectively recruited 103 patients from pediatric non-genetic subspecialty clinics, each with a clinical phenotype suggestive of an underlying genetic disorder, and compared the diagnostic yield and coverage of WGS with those of conventional genetic testing.ResultsWGS identified diagnostic variants in 41% of individuals, representing a significant increase over conventional testing results (24%; P = 0.01). Genes clinically sequenced in the cohort (n = 1,226) were well covered by WGS, with a median exonic coverage of 40 × ±8 × (mean ±SD). All the molecular diagnoses made by conventional methods were captured by WGS. The 18 new diagnoses made with WGS included structural and non-exonic sequence variants not detectable with whole-exome sequencing, and confirmed recent disease associations with the genes PIGG, RNU4ATAC, TRIO, and UNC13A.ConclusionWGS as a primary clinical test provided a higher diagnostic yield than conventional genetic testing in a clinically heterogeneous cohort.

The Personal Genome Project Canada: findings from whole genome sequences of the inaugural 56 participants.

BACKGROUND: The Personal Genome Project Canada is a comprehensive public data resource that integrates whole genome sequencing data and health information. We describe genomic variation identified in the initial recruitment cohort of 56 volunteers. METHODS: Volunteers were screened for eligibility and provided informed consent for open data sharing. Using blood DNA, we performed whole genome sequencing and identified all possible classes of DNA variants. A genetic counsellor explained the implication of the results to each participant. RESULTS: Whole genome sequencing of the first 56 participants identified 207 662 805 sequence variants and 27 494 copy number variations. We analyzed a prioritized disease-associated data set (n = 1606 variants) according to standardized guidelines, and interpreted 19 variants in 14 participants (25%) as having obvious health implications. Six of these variants (e.g., in BRCA1 or mosaic loss of an X chromosome) were pathogenic or likely pathogenic. Seven were risk factors for cancer, cardiovascular or neurobehavioural conditions. Four other variants - associated with cancer, cardiac or neurodegenerative phenotypes - remained of uncertain significance because of discrepancies among databases. We also identified a large structural chromosome aberration and a likely pathogenic mitochondrial variant. There were 172 recessive disease alleles (e.g., 5 individuals carried mutations for cystic fibrosis). Pharmacogenomics analyses revealed another 3.9 potentially relevant genotypes per individual. INTERPRETATION: Our analyses identified a spectrum of genetic variants with potential health impact in 25% of participants. When also considering recessive alleles and variants with potential pharmacologic relevance, all 56 participants had medically relevant findings. Although access is mostly limited to research, whole genome sequencing can provide specific and novel information with the potential of major impact for health care.

FOXP2 variants in 14 individuals with developmental speech and language disorders broaden the mutational and clinical spectrum.

BACKGROUND: Disruptions of the FOXP2 gene, encoding a forkhead transcription factor, are the first known monogenic cause of a speech and language disorder. So far, mainly chromosomal rearrangements such as translocations or larger deletions affecting FOXP2 have been reported. Intragenic deletions or convincingly pathogenic point mutations in FOXP2 have up to date only been reported in three families. We thus aimed at a further characterisation of the mutational and clinical spectrum. METHODS: Chromosomal microarray testing, trio exome sequencing, multigene panel sequencing and targeted sequencing of FOXP2 were performed in individuals with variable developmental disorders, and speech and language deficits. RESULTS: We identified four different truncating mutations, two novel missense mutations within the forkhead domain and an intragenic deletion in FOXP2 in 14 individuals from eight unrelated families. Mutations occurred de novo in four families and were inherited from an affected parent in the other four. All index patients presented with various manifestations of language and speech impairment. Apart from two individuals with normal onset of speech, age of first words was between 4 and 7 years. Articulation difficulties such as slurred speech, dyspraxia, stuttering and poor pronunciation were frequently noted. Motor development was normal or only mildly delayed. Mild cognitive impairment was reported for most individuals. CONCLUSIONS: By identifying intragenic deletions or mutations in 14 individuals from eight unrelated families with variable developmental delay/cognitive impairment and speech and language deficits, we considerably broaden the mutational and clinical spectrum associated with aberrations in FOXP2.

Haploinsufficiency of NR4A2 is associated with a neurodevelopmental phenotype with prominent language impairment.

Non-recurrent deletions in 2q24.1, minimally overlapping two genes, NR4A2 and GPD2, were recently described in individuals with language impairment and behavioral and cognitive symptoms. We herewith report on a female patient with a similar phenotype of severe language and mild cognitive impairment, in whom we identified a de novo deletion covering only NR4A2. NR4A2 encodes a transcription factor highly expressed in brain regions critical for speech and language and implicated in dopaminergic neuronal development. Our findings of a de novo deletion of NR4A2 in an individual with mild intellectual disability and prominent speech and language impairment provides further evidence for NR4A2 haploinsufficiency being causative for neurodevelopmental and particularly language phenotypes.

De novo mutations in beta-catenin (CTNNB1) appear to be a frequent cause of intellectual disability: expanding the mutational and clinical spectrum.

Recently, de novo heterozygous loss-of-function mutations in beta-catenin (CTNNB1) were described for the first time in four individuals with intellectual disability (ID), microcephaly, limited speech and (progressive) spasticity, and functional consequences of CTNNB1 deficiency were characterized in a mouse model. Beta-catenin is a key downstream component of the canonical Wnt signaling pathway. Somatic gain-of-function mutations have already been found in various tumor types, whereas germline loss-of-function mutations in animal models have been shown to influence neuronal development and maturation. We report on 16 additional individuals from 15 families in whom we newly identified de novo loss-of-function CTNNB1 mutations (six nonsense, five frameshift, one missense, two splice mutation, and one whole gene deletion). All patients have ID, motor delay and speech impairment (both mostly severe) and abnormal muscle tone (truncal hypotonia and distal hypertonia/spasticity). The craniofacial phenotype comprised microcephaly (typically -2 to -4 SD) in 12 of 16 and some overlapping facial features in all individuals (broad nasal tip, small alae nasi, long and/or flat philtrum, thin upper lip vermillion). With this detailed phenotypic characterization of 16 additional individuals, we expand and further establish the clinical and mutational spectrum of inactivating CTNNB1 mutations and thereby clinically delineate this new CTNNB1 haploinsufficiency syndrome.

Axonal transport deficits and degeneration can evolve independently in mouse models of amyotrophic lateral sclerosis.

Axonal transport deficits have been reported in many neurodegenerative conditions and are widely assumed to be an immediate causative step of axon and synapse loss. By imaging changes in axonal morphology and organelle transport over time in several animal models of amyotrophic lateral sclerosis (ALS), we now find that deficits in axonal transport of organelles (mitochondria, endosomes) and axon degeneration can evolve independently. This conclusion rests on the following results: (i) Axons can survive despite long-lasting transport deficits: In the SOD(G93A) model of ALS, transport deficits are detected soon after birth, months before the onset of axon degeneration. (ii) Transport deficits are not necessary for axon degeneration: In the SOD(G85R) model of ALS, motor axons degenerate, but transport is unaffected. (iii) Axon transport deficits are not sufficient to cause immediate degeneration: In mice that overexpress wild-type superoxide dismutase-1 (SOD(WT)), axons show chronic transport deficits, but survive. These data suggest that disturbances of organelle transport are not a necessary step in the emergence of motor neuron degeneration.

Expanding the clinical and mutational spectrum of Kaufman oculocerebrofacial syndrome with biallelic UBE3B mutations.

Biallelic mutations of UBE3B have recently been shown to cause Kaufman oculocerebrofacial syndrome (also reported as blepharophimosis-ptosis-intellectual disability syndrome), an autosomal recessive condition characterized by hypotonia, developmental delay, intellectual disability, congenital anomalies, characteristic facial dysmorphic features, and low cholesterol levels. To date, six patients with either missense mutations affecting the UBE3B HECT domain or truncating mutations have been described. Here, we report on the identification of homozygous or compound heterozygous UBE3B mutations in six additional patients from five unrelated families using either targeted UBE3B sequencing in individuals with suggestive facial dysmorphic features, or exome sequencing. Our results expand the clinical and mutational spectrum of the UBE3B-related disorder in several ways. First, we have identified UBE3B mutations in individuals who previously received distinct clinical diagnoses: two sibs with Toriello-Carey syndrome as well as the patient reported to have a "new" syndrome by Buntinx and Majewski in 1990. Second, we describe the adult phenotype and clinical variability of the syndrome. Third, we report on the first instance of homozygous missense alterations outside the HECT domain of UBE3B, observed in a patient with mildly dysmorphic facial features. We conclude that UBE3B mutations cause a clinically recognizable and possibly underdiagnosed syndrome characterized by distinct craniofacial features, hypotonia, failure to thrive, eye abnormalities, other congenital malformations, low cholesterol levels, and severe intellectual disability. We review the UBE3B-associated phenotypes, including forms that can mimick Toriello-Carey syndrome, and suggest the single designation "Kaufman oculocerebrofacial syndrome".