Expanding the phenotypic spectrum of NOTCH1 variants: clinical manifestations in families with congenital heart disease.

Pathogenic variants in NOTCH1 are associated with non-syndromic congenital heart disease (CHD) and Adams-Oliver syndrome (AOS). The clinical presentation of individuals with damaging NOTCH1 variants is characterized by variable expressivity and incomplete penetrance; however, data on systematic phenotypic characterization are limited. We report the genotype and phenotype of a cohort of 33 individuals (20 females, 13 males; median age 23.4 years, range 2.5-68.3 years) from 11 families with causative NOTCH1 variants (9 inherited, 2 de novo; 9 novel), ascertained from a proband with CHD. We describe the cardiac and extracardiac anomalies identified in these 33 individuals, only four of whom met criteria for AOS. The most common CHD identified was tetralogy of Fallot, though various left- and right-sided lesions and septal defects were also present. Extracardiac anomalies identified include cutis aplasia (5/33), cutaneous vascular anomalies (7/33), vascular anomalies of the central nervous system (2/10), Poland anomaly (1/33), pulmonary hypertension (2/33), and structural brain anomalies (3/14). Identification of these findings in a cardiac proband cohort supports NOTCH1-associated CHD and NOTCH1-associated AOS lying on a phenotypic continuum. Our findings also support (1) Broad indications for NOTCH1 molecular testing (any familial CHD, simplex tetralogy of Fallot or hypoplastic left heart); (2) Cascade testing in all at-risk relatives; and (3) A thorough physical exam, in addition to cardiac, brain (structural and vascular), abdominal, and ophthalmologic imaging, in all gene-positive individuals. This information is important for guiding the medical management of these individuals, particularly given the high prevalence of NOTCH1 variants in the CHD population.

CUX1-related neurodevelopmental disorder: deep insights into phenotype-genotype spectrum and underlying pathology.

Heterozygous, pathogenic CUX1 variants are associated with global developmental delay or intellectual disability. This study delineates the clinical presentation in an extended cohort and investigates the molecular mechanism underlying the disorder in a Cux1+/- mouse model. Through international collaboration, we assembled the phenotypic and molecular information for 34 individuals (23 unpublished individuals). We analyze brain CUX1 expression and susceptibility to epilepsy in Cux1+/- mice. We describe 34 individuals, from which 30 were unrelated, with 26 different null and four missense variants. The leading symptoms were mild to moderate delayed speech and motor development and borderline to moderate intellectual disability. Additional symptoms were muscular hypotonia, seizures, joint laxity, and abnormalities of the forehead. In Cux1+/- mice, we found delayed growth, histologically normal brains, and increased susceptibility to seizures. In Cux1+/- brains, the expression of Cux1 transcripts was half of WT animals. Expression of CUX1 proteins was reduced, although in early postnatal animals significantly more than in adults. In summary, disease-causing CUX1 variants result in a non-syndromic phenotype of developmental delay and intellectual disability. In some individuals, this phenotype ameliorates with age, resulting in a clinical catch-up and normal IQ in adulthood. The post-transcriptional balance of CUX1 expression in the heterozygous brain at late developmental stages appears important for this favorable clinical course.

Decreased left heart flow in fetal lambs causes left heart hypoplasia and pro-fibrotic tissue remodeling.

Low blood flow through the fetal left heart is often conjectured as an etiology for hypoplastic left heart syndrome (HLHS). To investigate if a decrease in left heart flow results in growth failure, we generate left ventricular inflow obstruction (LVIO) in mid-gestation fetal lambs by implanting coils in their left atrium using an ultrasound-guided percutaneous technique. Significant LVIO recapitulates important clinical features of HLHS: decreased antegrade aortic valve flow, compensatory retrograde perfusion of the brain and ascending aorta (AAo) from the arterial duct, severe left heart hypoplasia, a non-apex forming LV, and a thickened endocardial layer. The hypoplastic AAo have miRNA-gene pairs annotating to cell proliferation that are inversely differentially expressed by bulk RNA-seq. Single-nucleus RNA-seq of the hypoplastic LV myocardium shows an increase in fibroblasts with a reciprocal decrease in cardiomyocyte nuclei proportions. Fibroblasts, cardiomyocytes and endothelial cells from hypoplastic myocardium have increased expression of extracellular matrix component or fibrosis genes with dysregulated fibroblast growth factor signaling. Hence, a severe sustained ( ~ 1/3 gestation) reduction in fetal left heart flow is sufficient to cause left heart hypoplasia. This is accompanied by changes in cellular composition and gene expression consistent with a pro-fibrotic environment and aberrant induction of mesenchymal programs.

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.

Clinically Relevant Genetic Considerations for Patients With Tetralogy of Fallot.

Genetic changes affect embryogenesis, cardiac and extracardiac phenotype, development, later onset conditions, and both short- and long-term outcomes and comorbidities in the increasing population of individuals with tetralogy of Fallot (TOF). In this review, we focus on current knowledge about clinically relevant genetics for patients with TOF across the lifespan. The latest findings for TOF genetics that are pertinent to day-to-day practice and lifelong management are highlighted: morbidity/mortality, cardiac/extracardiac features, including neurodevelopmental expression, and recent changes to prenatal screening and diagnostics. Genome-wide microarray is the first-line clinical genetic test for TOF across the lifespan, detecting relevant structural changes including the most common for TOF, the 22q11.2 microdeletion. Accumulating evidence illustrates opportunities for advances in understanding and care that may arise from genetic diagnosis at any age. We also glimpse into the near future when the multigenic nature of TOF will be more fully revealed, further enhancing possibilities for preventive care. Precision medicine is nigh.

Dans la population croissante des personnes atteintes de la tétralogie de Fallot (TF), des modifications génétiques influencent l'embryogenèse, le développement, le phénotype cardiaque et extracardiaque, les complications tardives ainsi que les issues de santé et les états comorbides, à court et à long terme. Notre article de synthèse présente l'état des connaissances sur les renseignements génétiques cliniquement utiles pour les patients atteints de la TF tout au long de leur vie. Nous soulignons les découvertes récentes sur les aspects génétiques de la TF qui sont pertinentes pour la pratique clinique quotidienne et la prise en charge lors des différentes étapes de la vie : la morbidité et la mortalité, les caractéristiques cardiaques et extracardiaques (y compris l'expression neurodéveloppementale) et les changements récents touchant le dépistage et les diagnostics prénataux. La technologie de puce à ADN pour le génome entier constitue le test génétique clinique de première intention pour les personnes de tout âge atteintes de la TF, et elle permet la détection de modifications structurelles pertinentes dont celle le plus fréquemment associée à la TF, la microdélétion 22q11.2. L'utilité d'un diagnostic génétique pour améliorer la compréhension de la situation des patients de tous les âges et les soins qui leur sont offerts est de plus en plus mise en évidence. Nous entrevoyons également un avenir pas si lointain dans lequel la nature multigénique de la TF sera entièrement connue, ce qui ouvrira la voie à des soins préventifs bonifiés. La venue de la médecine de précision est imminente.

Genetic variants in DDX53 contribute to Autism Spectrum Disorder associated with the Xp22.11 locus.

Autism Spectrum Disorder (ASD) exhibits an ~4:1 male-to-female sex bias and is characterized by early-onset impairment of social/communication skills, restricted interests, and stereotyped behaviors. Disruption of the Xp22.11 locus has been associated with ASD in males. This locus includes the three-exon PTCHD1 gene, an adjacent multi-isoform long noncoding RNA (lncRNA) named PTCHD1-AS (spanning ~1Mb), and a poorly characterized single-exon RNA helicase named DDX53 that is intronic to PTCHD1-AS. While the relationship between PTCHD1/PTCHD1-AS and ASD is being studied, the role of DDX53 has not been examined, in part because there is no apparent functional murine orthologue. Through clinical testing, here, we identified 6 males and 1 female with ASD from 6 unrelated families carrying rare, predicted-damaging or loss-of-function variants in DDX53. Then, we examined databases, including the Autism Speaks MSSNG and Simons Foundation Autism Research Initiative, as well as population controls. We identified 24 additional individuals with ASD harboring rare, damaging DDX53 variations, including the same variants detected in two families from the original clinical analysis. In this extended cohort of 31 participants with ASD (28 male, 3 female), we identified 25 mostly maternally-inherited variations in DDX53, including 18 missense changes, 2 truncating variants, 2 in-frame variants, 2 deletions in the 3' UTR and 1 copy number deletion. Our findings in humans support a direct link between DDX53 and ASD, which will be important in clinical genetic testing. These same autism-related findings, coupled with the observation that a functional orthologous gene is not found in mouse, may also influence the design and interpretation of murine-modelling of ASD.

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.

Genome-wide rare variant score associates with morphological subtypes of autism spectrum disorder.

Defining different genetic subtypes of autism spectrum disorder (ASD) can enable the prediction of developmental outcomes. Based on minor physical and major congenital anomalies, we categorize 325 Canadian children with ASD into dysmorphic and nondysmorphic subgroups. We develop a method for calculating a patient-level, genome-wide rare variant score (GRVS) from whole-genome sequencing (WGS) data. GRVS is a sum of the number of variants in morphology-associated coding and non-coding regions, weighted by their effect sizes. Probands with dysmorphic ASD have a significantly higher GRVS compared to those with nondysmorphic ASD (P = 0.03). Using the polygenic transmission disequilibrium test, we observe an over-transmission of ASD-associated common variants in nondysmorphic ASD probands (P = 2.9 × 10-3). These findings replicate using WGS data from 442 ASD probands with accompanying morphology data from the Simons Simplex Collection. Our results provide support for an alternative genomic classification of ASD subgroups using morphology data, which may inform intervention protocols.

Biallelic PAN2 variants in individuals with a syndromic neurodevelopmental disorder and multiple congenital anomalies.

PAN2 encodes a subunit of a deadenylation complex with important functions in mRNA stability and post-transcriptional regulation of gene expression. A homozygous frameshift deletion in PAN2 was reported in a single affected individual with developmental delay and multiple congenital anomalies. Here, we describe five additional individuals from three unrelated families with homozygous predicted loss-of-function variants in PAN2. The affected individuals presented with significant overlap in their clinical features, including mild-moderate intellectual disability, hypotonia, sensorineural hearing loss, EEG abnormalities, congenital heart defects (tetralogy of Fallot, septal defects, dilated aortic root), urinary tract malformations, ophthalmological anomalies, short stature with other skeletal anomalies, and craniofacial features including flat occiput, ptosis, long philtrum, and short neck. Our data confirm that biallelic predicted loss-of-function variants in PAN2 cause a syndrome with multiple congenital anomalies, and suggest an important role of mRNA polyA tail length for proper organ formation.

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.

A recurrent SHANK3 frameshift variant in Autism Spectrum Disorder.

Autism Spectrum Disorder (ASD) is genetically complex with ~100 copy number variants and genes involved. To try to establish more definitive genotype and phenotype correlations in ASD, we searched genome sequence data, and the literature, for recurrent predicted damaging sequence-level variants affecting single genes. We identified 18 individuals from 16 unrelated families carrying a heterozygous guanine duplication (c.3679dup; p.Ala1227Glyfs*69) occurring within a string of 8 guanines (genomic location [hg38]g.50,721,512dup) affecting SHANK3, a prototypical ASD gene (0.08% of ASD-affected individuals carried the predicted p.Ala1227Glyfs*69 frameshift variant). Most probands carried de novo mutations, but five individuals in three families inherited it through somatic mosaicism. We scrutinized the phenotype of p.Ala1227Glyfs*69 carriers, and while everyone (17/17) formally tested for ASD carried a diagnosis, there was the variable expression of core ASD features both within and between families. Defining such recurrent mutational mechanisms underlying an ASD outcome is important for genetic counseling and early intervention.

Clinical Genetic Risk Variants Inform a Functional Protein Interaction Network for Tetralogy of Fallot.

BACKGROUND: Tetralogy of Fallot (TOF)-the most common cyanotic heart defect in newborns-has evidence of multiple genetic contributing factors. Identifying variants that are clinically relevant is essential to understand patient-specific disease susceptibility and outcomes and could contribute to delineating pathomechanisms. METHODS: Using a clinically driven strategy, we reanalyzed exome sequencing data from 811 probands with TOF, to identify rare loss-of-function and other likely pathogenic variants in genes associated with congenital heart disease. RESULTS: We confirmed a major contribution of likely pathogenic variants in FLT4 (VEGFR3 [vascular endothelial growth factor receptor 3]; n=14) and NOTCH1 (n=10) and identified 1 to 3 variants in each of 21 other genes, including ATRX, DLL4, EP300, GATA6, JAG1, NF1, PIK3CA, RAF1, RASA1, SMAD2, and TBX1. In addition, multiple loss-of-function variants provided support for 3 emerging congenital heart disease/TOF candidate genes: KDR (n=4), IQGAP1 (n=3), and GDF1 (n=8). In total, these variants were identified in 63 probands (7.8%). Using the 26 composite genes in a STRING protein interaction enrichment analysis revealed a biologically relevant network (P=3.3×10-16), with VEGFR2 (vascular endothelial growth factor receptor 2; KDR) and NOTCH1 (neurogenic locus notch homolog protein 1) representing central nodes. Variants associated with arrhythmias/sudden death and heart failure indicated factors that could influence long-term outcomes. CONCLUSIONS: The results are relevant to precision medicine for TOF. They suggest considerable clinical yield from genome-wide sequencing, with further evidence for KDR (VEGFR2) as a congenital heart disease/TOF gene and for VEGF (vascular endothelial growth factor) and Notch signaling as mechanisms in human disease. Harnessing the genetic heterogeneity of single gene defects could inform etiopathogenesis and help prioritize novel candidate genes for TOF.

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.

Assessment of the Implementation of Pharmacogenomic Testing in a Pediatric Tertiary Care Setting.

Importance: Pharmacogenomic (PGx) testing provides preemptive pharmacotherapeutic guidance regarding the lack of therapeutic benefit or adverse drug reactions of PGx targeted drugs. Pharmacogenomic information is of particular value among children with complex medical conditions who receive multiple medications and are at higher risk of developing adverse drug reactions. Objectives: To assess the implementation outcomes of a PGx testing program comprising both a point-of-care model that examined targeted drugs and a preemptive model informed by whole-genome sequencing that evaluated a broad range of drugs for potential therapy among children in a pediatric tertiary care setting. Design, Setting, and Participants: This cohort study was conducted at The Hospital for Sick Children in Toronto, Ontario, from January 2017 to September 2020. Pharmacogenomic analyses were performed among 172 children who were categorized into 2 groups: a point-of-care cohort and a preemptive cohort. The point-of-care cohort comprised 57 patients referred to the consultation clinic for planned therapy with PGx targeted drugs and/or for adverse drug reactions, including lack of therapeutic benefit, after the receipt of current or past medications. The preemptive cohort comprised 115 patients who received exploratory whole-genome sequencing-guided PGx testing for their heart conditions from the cardiac genome clinic at the Ted Rogers Centre for Heart Research. Exposures: Patients received PGx analysis of whole-genome sequencing data and/or multiplex genotyping of 6 pharmacogenes (CYP2C19, CYP2C9, CYP2D6, CYP3A5, VKORC1, and TPMT) that have established PGx clinical guidelines. Main Outcomes and Measures: The number of patients for whom PGx test results warranted deviation from standard dosing regimens. Results: A total of 172 children (mean [SD] age, 8.5 [5.6] years; 108 boys [62.8%]) were enrolled in the study. In the point-of-care cohort, a median of 2 target genes (range, 1-5 genes) were investigated per individual, with CYP2C19 being the most frequently examined; genotypes in 21 of 57 children (36.8%) were incompatible with standard treatment regimens. As expected from population allelic frequencies, among the 115 children in the whole-genome sequencing-guided preemptive cohort, 92 children (80.0%) were recommended to receive nonstandard treatment regimens for potential drug therapies based on their 6-gene pharmacogenetic profile. Conclusions and Relevance: In this cohort study, among both the point-of-care and preemptive cohorts, the multiplex PGx testing program provided dosing recommendations that deviated from standard regimens at an overall rate that was similar to the population frequencies of relevant variants.

EIF3F-related neurodevelopmental disorder: refining the phenotypic and expanding the molecular spectrum.

BACKGROUND: An identical homozygous missense variant in EIF3F, identified through a large-scale genome-wide sequencing approach, was reported as causative in nine individuals with a neurodevelopmental disorder, characterized by variable intellectual disability, epilepsy, behavioral problems and sensorineural hearing-loss. To refine the phenotypic and molecular spectrum of EIF3F-related neurodevelopmental disorder, we examined independent patients. RESULTS: 21 patients were homozygous and one compound heterozygous for c.694T>G/p.(Phe232Val) in EIF3F. Haplotype analyses in 15 families suggested that c.694T>G/p.(Phe232Val) was a founder variant. All affected individuals had developmental delays including delayed speech development. About half of the affected individuals had behavioral problems, altered muscular tone, hearing loss, and short stature. Moreover, this study suggests that microcephaly, reduced sensitivity to pain, cleft lip/palate, gastrointestinal symptoms and ophthalmological symptoms are part of the phenotypic spectrum. Minor dysmorphic features were observed, although neither the individuals' facial nor general appearance were obviously distinctive. Symptoms in the compound heterozygous individual with an additional truncating variant were at the severe end of the spectrum in regard to motor milestones, speech delay, organic problems and pre- and postnatal growth of body and head, suggesting some genotype-phenotype correlation. CONCLUSIONS: Our study refines the phenotypic and expands the molecular spectrum of EIF3F-related syndromic neurodevelopmental disorder.

A Distributed Whole Genome Sequencing Benchmark Study.

Population sequencing often requires collaboration across a distributed network of sequencing centers for the timely processing of thousands of samples. In such massive efforts, it is important that participating scientists can be confident that the accuracy of the sequence data produced is not affected by which center generates the data. A study was conducted across three established sequencing centers, located in Montreal, Toronto, and Vancouver, constituting Canada's Genomics Enterprise (www.cgen.ca). Whole genome sequencing was performed at each center, on three genomic DNA replicates from three well-characterized cell lines. Secondary analysis pipelines employed by each site were applied to sequence data from each of the sites, resulting in three datasets for each of four variables (cell line, replicate, sequencing center, and analysis pipeline), for a total of 81 datasets. These datasets were each assessed according to multiple quality metrics including concordance with benchmark variant truth sets to assess consistent quality across all three conditions for each variable. Three-way concordance analysis of variants across conditions for each variable was performed. Our results showed that the variant concordance between datasets differing only by sequencing center was similar to the concordance for datasets differing only by replicate, using the same analysis pipeline. We also showed that the statistically significant differences between datasets result from the analysis pipeline used, which can be unified and updated as new approaches become available. We conclude that genome sequencing projects can rely on the quality and reproducibility of aggregate data generated across a network of distributed sites.

Genes and Pathways Implicated in Tetralogy of Fallot Revealed by Ultra-Rare Variant Burden Analysis in 231 Genome Sequences.

Recent genome-wide studies of rare genetic variants have begun to implicate novel mechanisms for tetralogy of Fallot (TOF), a severe congenital heart defect (CHD). To provide statistical support for case-only data without parental genomes, we re-analyzed genome sequences of 231 individuals with TOF (n = 175) or related CHD. We adapted a burden test originally developed for de novo variants to assess ultra-rare variant burden in individual genes, and in gene-sets corresponding to functional pathways and mouse phenotypes, accounting for highly correlated gene-sets and for multiple testing. For truncating variants, the gene burden test confirmed significant burden in FLT4 (Bonferroni corrected p-value < 0.01). For missense variants, burden in NOTCH1 achieved genome-wide significance only when restricted to constrained genes (i.e., under negative selection, Bonferroni corrected p-value = 0.004), and showed enrichment for variants affecting the extracellular domain, especially those disrupting cysteine residues forming disulfide bonds (OR = 39.8 vs. gnomAD). Individuals with NOTCH1 ultra-rare missense variants, all with TOF, were enriched for positive family history of CHD. Other genes not previously implicated in CHD had more modest statistical support in gene burden tests. Gene-set burden tests for truncating variants identified a cluster of pathways corresponding to VEGF signaling (FDR = 0%), and of mouse phenotypes corresponding to abnormal vasculature (FDR = 0.8%); these suggested additional candidate genes not previously identified (e.g., WNT5A and ZFAND5). Results for the most promising genes were driven by the TOF subset of the cohort. The findings support the importance of ultra-rare variants disrupting genes involved in VEGF and NOTCH signaling in the genetic architecture of TOF, accounting for 11-14% of individuals in the TOF cohort. These proof-of-principle data indicate that this statistical methodology could assist in analyzing case-only sequencing data in which ultra-rare variants, whether de novo or inherited, contribute to the genetic etiopathogenesis of a complex disorder.

Genome Sequencing as a Diagnostic Test in Children With Unexplained Medical Complexity.

Importance: Children with medical complexity (CMC) represent a growing population in the pediatric health care system, with high resource use and associated health care costs. A genetic diagnosis can inform prognosis, anticipatory care, management, and reproductive planning. Conventional genetic testing strategies for CMC are often costly, time consuming, and ultimately unsuccessful. Objective: To evaluate the analytical and clinical validity of genome sequencing as a comprehensive diagnostic genetic test for CMC. Design, Setting, and Participants: In this cohort study of the prospective use of genome sequencing and comparison with standard-of-care genetic testing, CMC were recruited from May 1, 2017, to November 30, 2018, from a structured complex care program based at a tertiary care pediatric hospital in Toronto, Canada. Recruited CMC had at least 1 chronic condition, technology dependence (child is dependent at least part of each day on mechanical ventilators, and/or child requires prolonged intravenous administration of nutritional substances or drugs, and/or child is expected to have prolonged dependence on other device-based support), multiple subspecialist involvement, and substantial health care use. Review of the care plans for 545 CMC identified 143 suspected of having an undiagnosed genetic condition. Fifty-four families met inclusion criteria and were interested in participating, and 49 completed the study. Probands, similarly affected siblings, and biological parents were eligible for genome sequencing. Exposures: Genome sequencing was performed using blood-derived DNA from probands and family members using established methods and a bioinformatics pipeline for clinical genome annotation. Main Outcomes and Measures: The primary study outcome was the diagnostic yield of genome sequencing (proportion of CMC for whom the test result yielded a new diagnosis). Results: Genome sequencing was performed for 138 individuals from 49 families of CMC (29 male and 20 female probands; mean [SD] age, 7.0 [4.5] years). Genome sequencing detected all genomic variation previously identified by conventional genetic testing. A total of 15 probands (30.6%; 95% CI 19.5%-44.6%) received a new primary molecular genetic diagnosis after genome sequencing. Three individuals had novel diseases and an additional 9 had either ultrarare genetic conditions or rare genetic conditions with atypical features. At least 11 families received diagnostic information that had clinical management implications beyond genetic and reproductive counseling. Conclusions and Relevance: This study suggests that genome sequencing has high analytical and clinical validity and can result in new diagnoses in CMC even in the setting of extensive prior investigations. This clinical population may be enriched for ultrarare and novel genetic disorders. Genome sequencing is a potentially first-tier genetic test for CMC.

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.