Genomic Answers for Kids: Toward more equitable access to genomic testing for rare diseases in rural populations.

Next-generation sequencing has revolutionized the speed of rare disease (RD) diagnoses. While clinical exome and genome sequencing represent an effective tool for many RD diagnoses, there is room to further improve the diagnostic odyssey of many RD patients. One recognizable intervention lies in increasing equitable access to genomic testing. Rural communities represent a significant portion of underserved and underrepresented individuals facing additional barriers to diagnosis and treatment. Primary care providers (PCPs) at local clinics, though sometimes suspicious of a potential benefit of genetic testing for their patients, have significant constraints in pursuing it themselves and rely on referrals to specialists. Yet, these referrals are typically followed by long waitlists and significant delays in clinical assessment, insurance clearance, testing, and initiation of diagnosis-informed care management. Not only is this process time intensive, but it also often requires multiple visits to urban medical centers for which distance may be a significant barrier to rural families. Therefore, providing early, "direct-to-provider" (DTP) local access to unrestrictive genomic testing is likely to help speed up diagnostic times and access to care for RD patients in rural communities. In a pilot study with a PCP clinic in rural Kansas, we observed a minimum 5.5 months shortening of time to diagnosis through the DTP exome sequencing program as compared to rural patients receiving genetic testing through the "traditional" PCP-referral-to-specialist scheme. We share our experience to encourage future partnerships beyond our center. Our efforts represent just one step in fostering greater diversity and equity in genomic studies.

Deleterious, protein-altering variants in the transcriptional coregulator ZMYM3 in 27 individuals with a neurodevelopmental delay phenotype.

Neurodevelopmental disorders (NDDs) result from highly penetrant variation in hundreds of different genes, some of which have not yet been identified. Using the MatchMaker Exchange, we assembled a cohort of 27 individuals with rare, protein-altering variation in the transcriptional coregulator ZMYM3, located on the X chromosome. Most (n = 24) individuals were males, 17 of which have a maternally inherited variant; six individuals (4 male, 2 female) harbor de novo variants. Overlapping features included developmental delay, intellectual disability, behavioral abnormalities, and a specific facial gestalt in a subset of males. Variants in almost all individuals (n = 26) are missense, including six that recurrently affect two residues. Four unrelated probands were identified with inherited variation affecting Arg441, a site at which variation has been previously seen in NDD-affected siblings, and two individuals have de novo variation resulting in p.Arg1294Cys (c.3880C>T). All variants affect evolutionarily conserved sites, and most are predicted to damage protein structure or function. ZMYM3 is relatively intolerant to variation in the general population, is widely expressed across human tissues, and encodes a component of the KDM1A-RCOR1 chromatin-modifying complex. ChIP-seq experiments on one variant, p.Arg1274Trp, indicate dramatically reduced genomic occupancy, supporting a hypomorphic effect. While we are unable to perform statistical evaluations to definitively support a causative role for variation in ZMYM3, the totality of the evidence, including 27 affected individuals, recurrent variation at two codons, overlapping phenotypic features, protein-modeling data, evolutionary constraint, and experimentally confirmed functional effects strongly support ZMYM3 as an NDD-associated gene.

De novo EIF2AK1 and EIF2AK2 Variants Are Associated with Developmental Delay, Leukoencephalopathy, and Neurologic Decompensation.

EIF2AK1 and EIF2AK2 encode members of the eukaryotic translation initiation factor 2 alpha kinase (EIF2AK) family that inhibits protein synthesis in response to physiologic stress conditions. EIF2AK2 is also involved in innate immune response and the regulation of signal transduction, apoptosis, cell proliferation, and differentiation. Despite these findings, human disorders associated with deleterious variants in EIF2AK1 and EIF2AK2 have not been reported. Here, we describe the identification of nine unrelated individuals with heterozygous de novo missense variants in EIF2AK1 (1/9) or EIF2AK2 (8/9). Features seen in these nine individuals include white matter alterations (9/9), developmental delay (9/9), impaired language (9/9), cognitive impairment (8/9), ataxia (6/9), dysarthria in probands with verbal ability (6/9), hypotonia (7/9), hypertonia (6/9), and involuntary movements (3/9). Individuals with EIF2AK2 variants also exhibit neurological regression in the setting of febrile illness or infection. We use mammalian cell lines and proband-derived fibroblasts to further confirm the pathogenicity of variants in these genes and found reduced kinase activity. EIF2AKs phosphorylate eukaryotic translation initiation factor 2 subunit 1 (EIF2S1, also known as EIF2α), which then inhibits EIF2B activity. Deleterious variants in genes encoding EIF2B proteins cause childhood ataxia with central nervous system hypomyelination/vanishing white matter (CACH/VWM), a leukodystrophy characterized by neurologic regression in the setting of febrile illness and other stressors. Our findings indicate that EIF2AK2 missense variants cause a neurodevelopmental syndrome that may share phenotypic and pathogenic mechanisms with CACH/VWM.

Insurance denials and diagnostic rates in a pediatric genomic research cohort.

PURPOSE: This study aimed to assess the amount and types of clinical genetic testing denied by insurance and the rate of diagnostic and candidate genetic findings identified through research in patients who faced insurance denials. METHODS: Analysis consisted of review of insurance denials in 801 patients enrolled in a pediatric genomic research repository with either no previous genetic testing or previous negative genetic testing result identified through cross-referencing with insurance prior-authorizations in patient medical records. Patients and denials were also categorized by type of insurance coverage. Diagnostic findings and candidate genetic findings in these groups were determined through review of our internal variant database and patient charts. RESULTS: Of the 801 patients analyzed, 147 had insurance prior-authorization denials on record (18.3%). Exome sequencing and microarray were the most frequently denied genetic tests. Private insurance was significantly more likely to deny testing than public insurance (odds ratio = 2.03 [95% CI = 1.38-2.99] P = .0003). Of the 147 patients with insurance denials, 53.7% had at least 1 diagnostic or candidate finding and 10.9% specifically had a clinically diagnostic finding. Fifty percent of patients with clinically diagnostic results had immediate medical management changes (5.4% of all patients experiencing denials). CONCLUSION: Many patients face a major barrier to genetic testing in the form of lack of insurance coverage. A number of these patients have clinically diagnostic findings with medical management implications that would not have been identified without access to research testing. These findings support re-evaluation of insurance carriers' coverage policies.

Genetic heterogeneity and enrichment of variants in DNA-repair genes in ameloblastoma.

OBJECTIVE: Ameloblastomas are a group of relatively common odontogenic tumors that frequently originate from the dental epithelium. These tumors are aggressive in nature and present as slow-growing painless cortical expansion of the jaw. Histologically, the follicular and plexiform subtypes constitute two-thirds of solid/multicystic ameloblastomas. The objective of this study was to understand the genetic architecture of follicular and plexiform ameloblastomas using deep whole-exome sequencing. METHODS: Archived formalin-fixed paraffin-embedded tissue blocks of follicular (n = 4) and plexiform (n = 6) ameloblastomas were retrieved and genomic DNAs were isolated from the tumor tissue dissected from the formalin-fixed paraffin-embedded block. The exomes were enriched using the Integrated DNA Technologies Exome Research Panel (IDT, Coralville, IA) and paired-end sequencing was completed on an Illumina NovaSeq 6000 with an average output of 20 GB of data resulting in a mean coverage of 400×. Variant analysis was completed using custom-developed software: Rapid Understanding of Nucleotide variant Effect Software and variant integration and knowledge interpretation in genomes. RESULTS: Our analyses focused on examining somatic variants (gnomAD minor allele frequency ≤1%) in genes found on an Food and Drug Administration -approved clinical cancer sequencing panel (FoundationOne®CDx). In follicular tumors, variants (>20% of the reads) were identified in BRAF, KMT2D, and ABL1 genes. In plexiform tumors, variants (>20% of the reads) were identified in ALK, BRAF, KRAS, KMT2D, SMO, KMT2A, and BRCA2 genes. Enrichment analysis showed a significant role of DNA repair genes in the development of these tumors. CONCLUSION: The variants identified in follicular and plexiform ameloblastomas were enriched in DNA-repair genes. The observed genetic heterogeneity in these ameloblastomas may contribute to the aggressive nature and recurrence risk of these tumors.

Genomic answers for children: Dynamic analyses of >1000 pediatric rare disease genomes.

PURPOSE: This study aimed to provide comprehensive diagnostic and candidate analyses in a pediatric rare disease cohort through the Genomic Answers for Kids program. METHODS: Extensive analyses of 960 families with suspected genetic disorders included short-read exome sequencing and short-read genome sequencing (srGS); PacBio HiFi long-read genome sequencing (HiFi-GS); variant calling for single nucleotide variants (SNV), structural variant (SV), and repeat variants; and machine-learning variant prioritization. Structured phenotypes, prioritized variants, and pedigrees were stored in PhenoTips database, with data sharing through controlled access the database of Genotypes and Phenotypes. RESULTS: Diagnostic rates ranged from 11% in patients with prior negative genetic testing to 34.5% in naive patients. Incorporating SVs from genome sequencing added up to 13% of new diagnoses in previously unsolved cases. HiFi-GS yielded increased discovery rate with >4-fold more rare coding SVs compared with srGS. Variants and genes of unknown significance remain the most common finding (58% of nondiagnostic cases). CONCLUSION: Computational prioritization is efficient for diagnostic SNVs. Thorough identification of non-SNVs remains challenging and is partly mitigated using HiFi-GS sequencing. Importantly, community research is supported by sharing real-time data to accelerate gene validation and by providing HiFi variant (SNV/SV) resources from >1000 human alleles to facilitate implementation of new sequencing platforms for rare disease diagnoses.

Elucidating the clinical spectrum and molecular basis of HYAL2 deficiency.

PURPOSE: We previously defined biallelic HYAL2 variants causing a novel disorder in 2 families, involving orofacial clefting, facial dysmorphism, congenital heart disease, and ocular abnormalities, with Hyal2 knockout mice displaying similar phenotypes. In this study, we better define the phenotype and pathologic disease mechanism. METHODS: Clinical and genomic investigations were undertaken alongside molecular studies, including immunoblotting and immunofluorescence analyses of variant/wild-type human HYAL2 expressed in mouse fibroblasts, and in silico modeling of putative pathogenic variants. RESULTS: Ten newly identified individuals with this condition were investigated, and they were associated with 9 novel pathogenic variants. Clinical studies defined genotype-phenotype correlations and confirmed a recognizable craniofacial phenotype in addition to myopia, cleft lip/palate, and congenital cardiac anomalies as the most consistent manifestations of the condition. In silico modeling of missense variants identified likely deleterious effects on protein folding. Consistent with this, functional studies indicated that these variants cause protein instability and a concomitant cell surface absence of HYAL2 protein. CONCLUSION: These studies confirm an association between HYAL2 alterations and syndromic cleft lip/palate, provide experimental evidence for the pathogenicity of missense alleles, enable further insights into the pathomolecular basis of the disease, and delineate the core and variable clinical outcomes of the condition.

Clinical Validation of Genome Reference Consortium Human Build 38 in a Laboratory Utilizing Next-Generation Sequencing Technologies.

BACKGROUND: Laboratories utilizing next-generation sequencing align sequence data to a standardized human reference genome (HRG). Several updated versions, or builds, have been released since the original HRG in 2001, including the Genome Reference Consortium Human Build 38 (GRCh38) in 2013. However, most clinical laboratories still use GRCh37, which was released in 2009. We report our laboratory's clinical validation of GRCh38. METHODS: Migration to GRCh38 was validated by comparing the coordinates (lifting over) of 9443 internally curated variants from GRCh37 to GRCh38, globally comparing protein coding sequence variants aligned with GRCh37 vs GRCh38 from 917 exomes, assessing genes with known discrepancies, comparing coverage differences, and establishing the analytic sensitivity and specificity of variant detection using Genome in a Bottle data. RESULTS: Eight discrepancies, due to strand swap or reference base, were observed. Three clinically relevant variants had the GRCh37 alternate allele as the reference allele in GRCh38. A comparison of 88 295 calls between builds identified 8 disease-associated genes with sequence differences: ABO, BNC2, KIZ, NEFL, NR2E3, PTPRQ, SHANK2, and SRD5A2. Discrepancies in coding regions in GRCh37 were resolved in GRCh38. CONCLUSIONS: There were a small number of clinically significant changes between the 2 genome builds. GRCh38 provided improved detection of nucleotide changes due to the resolution of discrepancies present in GRCh37. Implementation of GRCh38 results in more accurate and consistent reporting.

Correction: Whole exome sequencing reveals HSPA1L as a genetic risk factor for spontaneous preterm birth.

[This corrects the article DOI: 10.1371/journal.pgen.1007394.].

Whole exome sequencing reveals HSPA1L as a genetic risk factor for spontaneous preterm birth.

Preterm birth is a leading cause of morbidity and mortality in infants. Genetic and environmental factors play a role in the susceptibility to preterm birth, but despite many investigations, the genetic basis for preterm birth remain largely unknown. Our objective was to identify rare, possibly damaging, nucleotide variants in mothers from families with recurrent spontaneous preterm births (SPTB). DNA samples from 17 Finnish mothers who delivered at least one infant preterm were subjected to whole exome sequencing. All mothers were of northern Finnish origin and were from seven multiplex families. Additional replication samples of European origin consisted of 93 Danish sister pairs (and two sister triads), all with a history of a preterm delivery. Rare exonic variants (frequency <1%) were analyzed to identify genes and pathways likely to affect SPTB susceptibility. We identified rare, possibly damaging, variants in genes that were common to multiple affected individuals. The glucocorticoid receptor signaling pathway was the most significant (p<1.7e-8) with genes containing these variants in a subgroup of ten Finnish mothers, each having had 2-4 SPTBs. This pathway was replicated among the Danish sister pairs. A gene in this pathway, heat shock protein family A (Hsp70) member 1 like (HSPA1L), contains two likely damaging missense alleles that were found in four different Finnish families. One of the variants (rs34620296) had a higher frequency in cases compared to controls (0.0025 vs. 0.0010, p = 0.002) in a large preterm birth genome-wide association study (GWAS) consisting of mothers of general European ancestry. Sister pairs in replication samples also shared rare, likely damaging HSPA1L variants. Furthermore, in silico analysis predicted an additional phosphorylation site generated by rs34620296 that could potentially affect chaperone activity or HSPA1L protein stability. Finally, in vitro functional experiment showed a link between HSPA1L activity and decidualization. In conclusion, rare, likely damaging, variants in HSPA1L were observed in multiple families with recurrent SPTB.

GPR37L1 modulates seizure susceptibility: Evidence from mouse studies and analyses of a human GPR37L1 variant.

Progressive myoclonus epilepsies (PMEs) are disorders characterized by myoclonic and generalized seizures with progressive neurological deterioration. While several genetic causes for PMEs have been identified, the underlying causes remain unknown for a substantial portion of cases. Here we describe several affected individuals from a large, consanguineous family presenting with a novel PME in which symptoms begin in adolescence and result in death by early adulthood. Whole exome analyses revealed that affected individuals have a homozygous variant in GPR37L1 (c.1047G>T [Lys349Asn]), an orphan G protein-coupled receptor (GPCR) expressed predominantly in the brain. In vitro studies demonstrated that the K349N substitution in Gpr37L1 did not grossly alter receptor expression, surface trafficking or constitutive signaling in transfected cells. However, in vivo studies revealed that a complete loss of Gpr37L1 function in mice results in increased seizure susceptibility. Mice lacking the related receptor Gpr37 also exhibited an increase in seizure susceptibility, while genetic deletion of both receptors resulted in an even more dramatic increase in vulnerability to seizures. These findings provide evidence linking GPR37L1 and GPR37 to seizure etiology and demonstrate an association between a GPR37L1 variant and a novel progressive myoclonus epilepsy.

Alström Syndrome: Mutation Spectrum of ALMS1.

Alström Syndrome (ALMS), a recessive, monogenic ciliopathy caused by mutations in ALMS1, is typically characterized by multisystem involvement including early cone-rod retinal dystrophy and blindness, hearing loss, childhood obesity, type 2 diabetes mellitus, cardiomyopathy, fibrosis, and multiple organ failure. The precise function of ALMS1 remains elusive, but roles in endosomal and ciliary transport and cell cycle regulation have been shown. The aim of our study was to further define the spectrum of ALMS1 mutations in patients with clinical features of ALMS. Mutational analysis in a world-wide cohort of 204 families identified 109 novel mutations, extending the number of known ALMS1 mutations to 239 and highlighting the allelic heterogeneity of this disorder. This study represents the most comprehensive mutation analysis in patients with ALMS, identifying the largest number of novel mutations in a single study worldwide. Here, we also provide an overview of all ALMS1 mutations identified to date.

Loss of function variants in human PNPLA8 encoding calcium-independent phospholipase A2 γ recapitulate the mitochondriopathy of the homologous null mouse.

Mitochondriopathies are a group of clinically heterogeneous genetic diseases caused by defects in mitochondrial metabolism, bioenergetic efficiency, and/or signaling functions. The large majority of proteins involved in mitochondrial function are encoded by nuclear genes, with many yet to be associated with human disease. We performed exome sequencing on a young girl with a suspected mitochondrial myopathy that manifested as progressive muscle weakness, hypotonia, seizures, poor weight gain, and lactic acidosis. She was compound heterozygous for two frameshift mutations, p.Asn112HisfsX29 and p.Leu659AlafsX4, in the PNPLA8 gene, which encodes mitochondrial calcium-independent phospholipase A2 γ (iPLA2 γ). Western blot analysis of affected muscle displayed the absence of PNPLA8 protein. iPLA2 s are critical mediators of a variety of cellular processes including growth, metabolism, and lipid second messenger generation, exerting their functions through catalyzing the cleavage of the acyl groups in glycerophospholipids. The clinical presentation, muscle histology and the mitochondrial ultrastructural abnormalities of this proband are highly reminiscent of Pnpla8 null mice. Although other iPLA2 -related diseases have been identified, namely, infantile neuroaxonal dystrophy and neutral lipid storage disease with myopathy, this is the first report of PNPLA8-related disease in a human. We suggest PNPLA8 join the increasing list of human genes involved in lipid metabolism associated with neuromuscular diseases due to mitochondrial dysfunction.

Iron deficiency drives an autosomal dominant hypophosphatemic rickets (ADHR) phenotype in fibroblast growth factor-23 (Fgf23) knock-in mice.

Autosomal dominant hypophosphatemic rickets (ADHR) is unique among the disorders involving Fibroblast growth factor 23 (FGF23) because individuals with R176Q/W and R179Q/W mutations in the FGF23 (176)RXXR(179)/S(180) proteolytic cleavage motif can cycle from unaffected status to delayed onset of disease. This onset may occur in physiological states associated with iron deficiency, including puberty and pregnancy. To test the role of iron status in development of the ADHR phenotype, WT and R176Q-Fgf23 knock-in (ADHR) mice were placed on control or low-iron diets. Both the WT and ADHR mice receiving low-iron diet had significantly elevated bone Fgf23 mRNA. WT mice on a low-iron diet maintained normal serum intact Fgf23 and phosphate metabolism, with elevated serum C-terminal Fgf23 fragments. In contrast, the ADHR mice on the low-iron diet had elevated intact and C-terminal Fgf23 with hypophosphatemic osteomalacia. We used in vitro iron chelation to isolate the effects of iron deficiency on Fgf23 expression. We found that iron chelation in vitro resulted in a significant increase in Fgf23 mRNA that was dependent upon Mapk. Thus, unlike other syndromes of elevated FGF23, our findings support the concept that late-onset ADHR is the product of gene-environment interactions whereby the combined presence of an Fgf23-stabilizing mutation and iron deficiency can lead to ADHR.

Diagnosis of mitochondrial disorders by concomitant next-generation sequencing of the exome and mitochondrial genome.

Mitochondrial diseases are notoriously difficult to diagnose due to extreme locus and allelic heterogeneity, with both nuclear and mitochondrial genomes potentially liable. Using exome sequencing we demonstrate the ability to rapidly and cost effectively evaluate both the nuclear and mitochondrial genomes to obtain a molecular diagnosis for four patients with three distinct mitochondrial disorders. One patient was found to have Leigh syndrome due to a mutation in MT-ATP6, two affected siblings were discovered to be compound heterozygous for mutations in the NDUFV1 gene, which causes mitochondrial complex I deficiency, and one patient was found to have coenzyme Q10 deficiency due to compound heterozygous mutations in COQ2. In all cases conventional diagnostic testing failed to identify a molecular diagnosis. We suggest that additional studies should be conducted to evaluate exome sequencing as a primary diagnostic test for mitochondrial diseases, including those due to mtDNA mutations.

Pangenome graphs improve the analysis of structural variants in rare genetic diseases.

Rare DNA alterations that cause heritable diseases are only partially resolvable by clinical next-generation sequencing due to the difficulty of detecting structural variation (SV) in all genomic contexts. Long-read, high fidelity genome sequencing (HiFi-GS) detects SVs with increased sensitivity and enables assembling personal and graph genomes. We leverage standard reference genomes, public assemblies (n = 94) and a large collection of HiFi-GS data from a rare disease program (Genomic Answers for Kids, GA4K, n = 574 assemblies) to build a graph genome representing a unified SV callset in GA4K, identify common variation and prioritize SVs that are more likely to cause genetic disease (MAF < 0.01). Using graphs, we obtain a higher level of reproducibility than the standard reference approach. We observe over 200,000 SV alleles unique to GA4K, including nearly 1000 rare variants that impact coding sequence. With improved specificity for rare SVs, we isolate 30 candidate SVs in phenotypically prioritized genes, including known disease SVs. We isolate a novel diagnostic SV in KMT2E, demonstrating use of personal assemblies coupled with pangenome graphs for rare disease genomics. The community may interrogate our pangenome with additional assemblies to discover new SVs within the allele frequency spectrum relevant to genetic diseases.

Complex trait associations in rare diseases and impacts on Mendelian variant interpretation.

Emerging evidence implicates common genetic variation - aggregated into polygenic scores (PGS) - impacting the onset and phenotypic presentation of rare diseases. In this study, we quantified individual polygenic liability for 1,151 previously published PGS in a cohort of 2,374 probands enrolled in the Genomic Answers for Kids (GA4K) rare disease study, revealing widespread associations between rare disease phenotypes and PGSs for common complex diseases and traits, blood protein levels, and brain and other organ morphological measurements. We observed increased polygenic burden in probands with variants of unknown significance (VUS) compared to unaffected carrier parents. We further observed an enrichment in overlap between diagnostic and candidate rare disease genes and large-effect PGS genes. Overall, our study supports and expands on previous findings of complex trait associations in rare disease phenotypes and provides a framework for identifying novel candidate rare disease genes and in understanding variable penetrance of candidate Mendelian disease variants.

Significance Associated with Phenotype Score Aids in Variant Prioritization for Exome Sequencing Analysis.

Several in silico annotation-based methods have been developed to prioritize variants in exome sequencing analysis. This study introduced a novel metric Significance Associated with Phenotypes (SAP) score, which generates a statistical score by comparing an individual's observed phenotypes against existing gene-phenotype associations. To evaluate the SAP score, a retrospective analysis was performed on 219 exomes. Among them, 82 family-based and 35 singleton exomes had at least one disease-causing variant that explained the patient's clinical features. SAP scores were calculated, and the rank of the disease-causing variant was compared with a known method, Exomiser. Using the SAP score, the known causative variant was ranked in the top 10 retained variants for 94% (77 of 82) of the family-based exomes and in first place for 73% of these cases. For singleton exomes, the SAP score analysis ranked the known pathogenic variants within the top 10 for 80% (28 of 35) of cases. The SAP score, which is independent of detected variants, demonstrates comparable performance with Exomiser, which considers both phenotype and variant-level evidence simultaneously. Among 102 cases with negative results or variants of uncertain significance, SAP score analysis revealed two cases with a potential new diagnosis based on rank. The SAP score, a phenotypic quantitative metric, can be used in conjunction with standard variant filtration and annotation to enhance variant prioritization in exome analysis.