Comprehensive whole-genome sequence analyses provide insights into the genomic architecture of cerebral palsy.

We performed whole-genome sequencing (WGS) in 327 children with cerebral palsy (CP) and their biological parents. We classified 37 of 327 (11.3%) children as having pathogenic/likely pathogenic (P/LP) variants and 58 of 327 (17.7%) as having variants of uncertain significance. Multiple classes of P/LP variants included single-nucleotide variants (SNVs)/indels (6.7%), copy number variations (3.4%) and mitochondrial mutations (1.5%). The COL4A1 gene had the most P/LP SNVs. We also analyzed two pediatric control cohorts (n = 203 trios and n = 89 sib-pair families) to provide a baseline for de novo mutation rates and genetic burden analyses, the latter of which demonstrated associations between de novo deleterious variants and genes related to the nervous system. An enrichment analysis revealed previously undescribed plausible candidate CP genes (SMOC1, KDM5B, BCL11A and CYP51A1). A multifactorial CP risk profile and substantial presence of P/LP variants combine to support WGS in the diagnostic work-up across all CP and related phenotypes.

A call for increased inclusivity and global representation in pharmacogenetic testing.

Commercial pharmacogenetic testing panels capture a fraction of the genetic variation underlying medication metabolism and predisposition to adverse reactions. In this study we compared variation in six pharmacogenes detected by whole genome sequencing (WGS) to a targeted commercial panel in a cohort of 308 individuals with family history of pediatric heart disease. In 1% of the cohort, WGS identified rare variants that altered the interpretation of metabolizer status and would thus prevent potential errors in gene-based dosing.

Gene copy number variation and pediatric mental health/neurodevelopment in a general population.

We assessed the relationship of gene copy number variation (CNV) in mental health/neurodevelopmental traits and diagnoses, physical health and cognition in a community sample of 7100 unrelated children and youth of European or East Asian ancestry (Spit for Science). Clinically significant or susceptibility CNVs were present in 3.9% of participants and were associated with elevated scores on a continuous measure of attention-deficit/hyperactivity disorder (ADHD) traits (P = 5.0 × 10-3), longer response inhibition (a cognitive deficit found in several mental health and neurodevelopmental disorders; P = 1.0 × 10-2) and increased prevalence of mental health diagnoses (P = 1.9 × 10-6, odds ratio: 3.09), specifically ADHD, autism spectrum disorder anxiety and learning problems/learning disorder (P's < 0.01). There was an increased burden of rare deletions in gene-sets related to brain function or expression in brain associated with more ADHD traits. With the current mental health crisis, our data established a baseline for delineating genetic contributors in pediatric-onset conditions.

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 ).

GeneTerpret: a customizable multilayer approach to genomic variant prioritization and interpretation.

BACKGROUND: Variant interpretation is the main bottleneck in medical genomic sequencing efforts. This usually involves genome analysts manually searching through a multitude of independent databases, often with the aid of several, mostly independent, computational tools. To streamline variant interpretation, we developed the GeneTerpret platform which collates data from current interpretation tools and databases, and applies a phenotype-driven query to categorize the variants identified in the genome(s). The platform assigns quantitative validity scores to genes by query and assembly of the genotype-phenotype data, sequence homology, molecular interactions, expression data, and animal models. It also uses the American College of Medical Genetics and Genomics (ACMG) criteria to categorize variants into five tiers of pathogenicity. The final output is a prioritized list of potentially causal variants/genes. RESULTS: We tested GeneTerpret by comparing its performance to expert-curated genes (ClinGen's gene-validity database) and variant pathogenicity reports (DECIPHER database). Output from GeneTerpret was 97.2% and 83.5% concordant with the expert-curated sources, respectively. Additionally, similar concordance was observed when GeneTerpret's performance was compared with our internal expert-interpreted clinical datasets. CONCLUSIONS: GeneTerpret is a flexible platform designed to streamline the genome interpretation process, through a unique interface, with improved ease, speed and accuracy. This modular and customizable system allows the user to tailor the component-programs in the analysis process to their preference. GeneTerpret is available online at https://geneterpret.com .

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.

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.

Genome sequencing broadens the range of contributing variants with clinical implications in schizophrenia.

The range of genetic variation with potential clinical implications in schizophrenia, beyond rare copy number variants (CNVs), remains uncertain. We therefore analyzed genome sequencing data for 259 unrelated adults with schizophrenia from a well-characterized community-based cohort previously examined with chromosomal microarray for CNVs (none with 22q11.2 deletions). We analyzed these genomes for rare high-impact variants considered causal for neurodevelopmental disorders, including single-nucleotide variants (SNVs) and small insertions/deletions (indels), for potential clinical relevance based on findings for neurodevelopmental disorders. Also, we investigated a novel variant type, tandem repeat expansions (TREs), in 45 loci known to be associated with monogenic neurological diseases. We found several of these variants in this schizophrenia population suggesting that these variants have a wider clinical spectrum than previously thought. In addition to known pathogenic CNVs, we identified 11 (4.3%) individuals with clinically relevant SNVs/indels in genes converging on schizophrenia-relevant pathways. Clinical yield was significantly enriched in females and in those with broadly defined learning/intellectual disabilities. Genome analyses also identified variants with potential clinical implications, including TREs (one in DMPK; two in ATXN8OS) and ultra-rare loss-of-function SNVs in ZMYM2 (a novel candidate gene for schizophrenia). Of the 233 individuals with no pathogenic CNVs, we identified rare high-impact variants (i.e., clinically relevant or with potential clinical implications) for 14 individuals (6.0%); some had multiple rare high-impact variants. Mean schizophrenia polygenic risk score was similar between individuals with and without clinically relevant rare genetic variation; common variants were not sufficient for clinical application. These findings broaden the individual and global picture of clinically relevant genetic risk in schizophrenia, and suggest the potential translational value of genome sequencing as a single genetic technology for schizophrenia.

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.

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.

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.

Hybrid-controlled neurofuzzy networks analysis resulting in genetic regulatory networks reconstruction.

Reverse engineering of gene regulatory networks (GRNs) is the process of estimating genetic interactions of a cellular system from gene expression data. In this paper, we propose a novel hybrid systematic algorithm based on neurofuzzy network for reconstructing GRNs from observational gene expression data when only a medium-small number of measurements are available. The approach uses fuzzy logic to transform gene expression values into qualitative descriptors that can be evaluated by using a set of defined rules. The algorithm uses neurofuzzy network to model genes effects on other genes followed by four stages of decision making to extract gene interactions. One of the main features of the proposed algorithm is that an optimal number of fuzzy rules can be easily and rapidly extracted without overparameterizing. Data analysis and simulation are conducted on microarray expression profiles of S. cerevisiae cell cycle and demonstrate that the proposed algorithm not only selects the patterns of the time series gene expression data accurately, but also provides models with better reconstruction accuracy when compared with four published algorithms: DBNs, VBEM, time delay ARACNE, and PF subjected to LASSO. The accuracy of the proposed approach is evaluated in terms of recall and F-score for the network reconstruction task.