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.

Rare de novo gain-of-function missense variants in DOT1L are associated with developmental delay and congenital anomalies.

Misregulation of histone lysine methylation is associated with several human cancers and with human developmental disorders. DOT1L is an evolutionarily conserved gene encoding a lysine methyltransferase (KMT) that methylates histone 3 lysine-79 (H3K79) and was not previously associated with a Mendelian disease in OMIM. We have identified nine unrelated individuals with seven different de novo heterozygous missense variants in DOT1L through the Undiagnosed Disease Network (UDN), the SickKids Complex Care genomics project, and GeneMatcher. All probands had some degree of global developmental delay/intellectual disability, and most had one or more major congenital anomalies. To assess the pathogenicity of the DOT1L variants, functional studies were performed in Drosophila and human cells. The fruit fly DOT1L ortholog, grappa, is expressed in most cells including neurons in the central nervous system. The identified DOT1L variants behave as gain-of-function alleles in flies and lead to increased H3K79 methylation levels in flies and human cells. Our results show that human DOT1L and fly grappa are required for proper development and that de novo heterozygous variants in DOT1L are associated with a Mendelian disease.

Monoallelic loss-of-function BMP2 variants result in BMP2-related skeletal dysplasia spectrum.

PURPOSE: Bone morphogenic proteins (BMPs) regulate gene expression that is related to many critical developmental processes, including osteogenesis for which they are named. In addition, BMP2 is widely expressed in cells of mesenchymal origin, including bone, cartilage, skeletal and cardiac muscle, and adipose tissue. It also participates in neurodevelopment by inducing differentiation of neural stem cells. In humans, BMP2 variants result in a multiple congenital anomaly syndrome through a haploinsufficiency mechanism. We sought to expand the phenotypic spectrum and highlight phenotypes of patients harboring monoallelic missense variants in BMP2. METHODS: We used retrospective chart review to examine phenotypes from an international cohort of 18 individuals and compared these with published cases. Patient-derived missense variants were modeled in zebrafish to examine their effect on the ability of bmp2b to promote embryonic ventralization. RESULTS: The presented cases recapitulated existing descriptions of BMP2-related disorders, including craniofacial, cardiac, and skeletal anomalies and exhibit a wide phenotypic spectrum. We also identified patients with neural tube defects, structural brain anomalies, and endocrinopathies. Missense variants modeled in zebrafish resulted in loss of protein function. CONCLUSION: We use this expansion of reported phenotypes to suggest multidisciplinary medical monitoring and management of patients with BMP2-related skeletal dysplasia spectrum.

Finding the sweet spot: a qualitative study exploring patients' acceptability of chatbots in genetic service delivery.

Chatbots, web-based artificial intelligence tools that simulate human conversation, are increasingly in use to support many areas of genomic medicine. However, patient preferences towards using chatbots across the range of clinical settings are unknown. We conducted a qualitative study with individuals who underwent genetic testing for themselves or their child. Participants were asked about their preferences for using a chatbot within the genetic testing journey. Thematic analysis employing interpretive description was used. We interviewed 30 participants (67% female, 50% 50 + years). Participants considered chatbots to be inefficient for very simple tasks (e.g., answering FAQs) or very complex tasks (e.g., explaining results). Chatbots were acceptable for moderately complex tasks where participants perceived a favorable return on their investment of time and energy. In addition to achieving this "sweet spot," participants anticipated that their comfort with chatbots would increase if the chatbot was used as a complement to but not a replacement for usual care. Participants wanted a "safety net" (i.e., access to a clinician) for needs not addressed by the chatbot. This study provides timely insights into patients' comfort with and perceived limitations of chatbots for genomic medicine and can inform their implementation in practice.

Bridging clinical care and research in Ontario, Canada: Maximizing diagnoses from reanalysis of clinical exome sequencing data.

We examined the utility of clinical and research processes in the reanalysis of publicly-funded clinical exome sequencing data in Ontario, Canada. In partnership with eight sites, we recruited 287 families with suspected rare genetic diseases tested between 2014 and 2020. Data from seven laboratories was reanalyzed with the referring clinicians. Reanalysis of clinically relevant genes identified diagnoses in 4% (13/287); four were missed by clinical testing. Translational research methods, including analysis of novel candidate genes, identified candidates in 21% (61/287). Of these, 24 families have additional evidence through data sharing to support likely diagnoses (8% of cohort). This study indicates few diagnoses are missed by clinical laboratories, the incremental gain from reanalysis of clinically-relevant genes is modest, and the highest yield comes from validation of novel disease-gene associations. Future implementation of translational research methods, including continued reporting of compelling genes of uncertain significance by clinical laboratories, should be considered to maximize diagnoses.

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

EEF1A2 pathogenic variant presenting in an infant with failure to thrive and frequent apneas requiring respiratory support.

EEF1A2 is a gene whose protein product, eukaryotic translation elongation factor 1 alpha 2 (eEF1A2), plays an important role in neurodevelopment. Reports of individuals with pathogenic variants in EEF1A2 are rare, with less than 40 individuals reported world-wide, however a common feature is the association of the variant with developmental and epileptic encephalopathy. Thus far, there have been limited reports of other organ systems or body functions affected by variants in this gene. Here, we present a case of a child with EEF1A2-related disorder who presented at 3 months of age with hypotonia, microcephaly, failure to thrive, and respiratory insufficiency with central apneas requiring respiratory support. Our case highlights the notion that the respiratory system may be highly implicated in EEF1A2-related disorder, allowing for better phenotypic characterization of the disorder.

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.

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.

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.

Neurodegenerative VPS41 variants inhibit HOPS function and mTORC1-dependent TFEB/TFE3 regulation.

Vacuolar protein sorting 41 (VPS41) is as part of the Homotypic fusion and Protein Sorting (HOPS) complex required for lysosomal fusion events and, independent of HOPS, for regulated secretion. Here, we report three patients with compound heterozygous mutations in VPS41 (VPS41S285P and VPS41R662* ; VPS41c.1423-2A>G and VPS41R662* ) displaying neurodegeneration with ataxia and dystonia. Cellular consequences were investigated in patient fibroblasts and VPS41-depleted HeLa cells. All mutants prevented formation of a functional HOPS complex, causing delayed lysosomal delivery of endocytic and autophagic cargo. By contrast, VPS41S285P enabled regulated secretion. Strikingly, loss of VPS41 function caused a cytosolic redistribution of mTORC1, continuous nuclear localization of Transcription Factor E3 (TFE3), enhanced levels of LC3II, and a reduced autophagic response to nutrient starvation. Phosphorylation of mTORC1 substrates S6K1 and 4EBP1 was not affected. In a C. elegans model of Parkinson's disease, co-expression of VPS41S285P /VPS41R662* abolished the neuroprotective function of VPS41 against α-synuclein aggregates. We conclude that the VPS41 variants specifically abrogate HOPS function, which interferes with the TFEB/TFE3 axis of mTORC1 signaling, and cause a neurodegenerative disease.

Heterozygous loss of WBP11 function causes multiple congenital defects in humans and mice.

The genetic causes of multiple congenital anomalies are incompletely understood. Here, we report novel heterozygous predicted loss-of-function (LoF) and predicted damaging missense variants in the WW domain binding protein 11 (WBP11) gene in seven unrelated families with a variety of overlapping congenital malformations, including cardiac, vertebral, tracheo-esophageal, renal and limb defects. WBP11 encodes a component of the spliceosome with the ability to activate pre-messenger RNA splicing. We generated a Wbp11 null allele in mouse using CRISPR-Cas9 targeting. Wbp11 homozygous null embryos die prior to E8.5, indicating that Wbp11 is essential for development. Fewer Wbp11 heterozygous null mice are found than expected due to embryonic and postnatal death. Importantly, Wbp11 heterozygous null mice are small and exhibit defects in axial skeleton, kidneys and esophagus, similar to the affected individuals, supporting the role of WBP11 haploinsufficiency in the development of congenital malformations in humans. LoF WBP11 variants should be considered as a possible cause of VACTERL association as well as isolated Klippel-Feil syndrome, renal agenesis or esophageal atresia.

YIF1B mutations cause a post-natal neurodevelopmental syndrome associated with Golgi and primary cilium alterations.

Human post-natal neurodevelopmental delay is often associated with cerebral alterations that can lead, by themselves or associated with peripheral deficits, to premature death. Here, we report the clinical features of 10 patients from six independent families with mutations in the autosomal YIF1B gene encoding a ubiquitous protein involved in anterograde traffic from the endoplasmic reticulum to the cell membrane, and in Golgi apparatus morphology. The patients displayed global developmental delay, motor delay, visual deficits with brain MRI evidence of ventricle enlargement, myelination alterations and cerebellar atrophy. A similar profile was observed in the Yif1b knockout (KO) mouse model developed to identify the cellular alterations involved in the clinical defects. In the CNS, mice lacking Yif1b displayed neuronal reduction, altered myelination of the motor cortex, cerebellar atrophy, enlargement of the ventricles, and subcellular alterations of endoplasmic reticulum and Golgi apparatus compartments. Remarkably, although YIF1B was not detected in primary cilia, biallelic YIF1B mutations caused primary cilia abnormalities in skin fibroblasts from both patients and Yif1b-KO mice, and in ciliary architectural components in the Yif1b-KO brain. Consequently, our findings identify YIF1B as an essential gene in early post-natal development in human, and provide a new genetic target that should be tested in patients developing a neurodevelopmental delay during the first year of life. Thus, our work is the first description of a functional deficit linking Golgipathies and ciliopathies, diseases so far associated exclusively to mutations in genes coding for proteins expressed within the primary cilium or related ultrastructures. We therefore propose that these pathologies should be considered as belonging to a larger class of neurodevelopmental diseases depending on proteins involved in the trafficking of proteins towards specific cell membrane compartments.

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.

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.

Phenotypes and genotypes of mitochondrial aminoacyl-tRNA synthetase deficiencies from a single neurometabolic clinic.

Mitochondrial aminoacyl-tRNA synthetases play a major role in protein translation, synthesis, and oxidative phosphorylation. We reviewed all patients diagnosed with mitochondrial aminoacyl-tRNA synthetase deficiencies diagnosed in a single neurometabolic clinic. We report five patients with mitochondrial aminoacyl-tRNA synthetase deficiencies including DARS2, EARS2, PARS2, and RARS2 deficiencies. Siblings with DARS2 deficiency presented with global developmental delay within the first year of life. DARS2, EARS2, PARS2, and RARS2 deficiencies were identified by whole exome sequencing. We report coagulation factor abnormalities in PARS2 deficiency for the first time. We also report symmetric increased signal intensity in globus pallidi in FLAIR images in brain MRI in EARS2 deficiency for the first time. One patient with RARS2 deficiency had compound heterozygous variants in RARS2. One of those variants was an intronic variant. We confirmed the pathogenicity by mRNA studies. Mitochondrial aminoacyl-tRNA synthetase deficiencies are diagnosed by molecular genetic investigations. Clinically available non-invasive biochemical investigations are non-specific for the diagnosis of mitochondrial aminoacyl-tRNA synthetase deficiencies. A combination of brain MRI features and molecular genetic investigations should be undertaken to confirm the diagnosis of mitochondrial aminoacyl-tRNA synthetase deficiencies.

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.