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.

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.

A syndromic neurodevelopmental disorder caused by rare variants in PPFIA3.

PPFIA3 encodes the protein-tyrosine phosphatase, receptor-type, F-polypeptide-interacting-protein-alpha-3 (PPFIA3), which is a member of the LAR-protein-tyrosine phosphatase-interacting-protein (liprin) family involved in synapse formation and function, synaptic vesicle transport, and presynaptic active zone assembly. The protein structure and function are evolutionarily well conserved, but human diseases related to PPFIA3 dysfunction are not yet reported in OMIM. Here, we report 20 individuals with rare PPFIA3 variants (19 heterozygous and 1 compound heterozygous) presenting with developmental delay, intellectual disability, hypotonia, dysmorphisms, microcephaly or macrocephaly, autistic features, and epilepsy with reduced penetrance. Seventeen unique PPFIA3 variants were detected in 18 families. To determine the pathogenicity of PPFIA3 variants in vivo, we generated transgenic fruit flies producing either human wild-type (WT) PPFIA3 or five missense variants using GAL4-UAS targeted gene expression systems. In the fly overexpression assays, we found that the PPFIA3 variants in the region encoding the N-terminal coiled-coil domain exhibited stronger phenotypes compared to those affecting the C-terminal region. In the loss-of-function fly assay, we show that the homozygous loss of fly Liprin-α leads to embryonic lethality. This lethality is partially rescued by the expression of human PPFIA3 WT, suggesting human PPFIA3 function is partially conserved in the fly. However, two of the tested variants failed to rescue the lethality at the larval stage and one variant failed to rescue lethality at the adult stage. Altogether, the human and fruit fly data reveal that the rare PPFIA3 variants are dominant-negative loss-of-function alleles that perturb multiple developmental processes and synapse formation.

The landscape of reported VUS in multi-gene panel and genomic testing: Time for a change.

PURPOSE: Variants of uncertain significance (VUS) are a common result of diagnostic genetic testing and can be difficult to manage with potential misinterpretation and downstream costs, including time investment by clinicians. We investigated the rate of VUS reported on diagnostic testing via multi-gene panels (MGPs) and exome and genome sequencing (ES/GS) to measure the magnitude of uncertain results and explore ways to reduce their potentially detrimental impact. METHODS: Rates of inconclusive results due to VUS were collected from over 1.5 million sequencing test results from 19 clinical laboratories in North America from 2020 to 2021. RESULTS: We found a lower rate of inconclusive test results due to VUSs from ES/GS (22.5%) compared with MGPs (32.6%; P < .0001). For MGPs, the rate of inconclusive results correlated with panel size. The use of trios reduced inconclusive rates (18.9% vs 27.6%; P < .0001), whereas the use of GS compared with ES had no impact (22.2% vs 22.6%; P = ns). CONCLUSION: The high rate of VUS observed in diagnostic MGP testing warrants examining current variant reporting practices. We propose several approaches to reduce reported VUS rates, while directing clinician resources toward important VUS follow-up.

Direct haplotype-resolved 5-base HiFi sequencing for genome-wide profiling of hypermethylation outliers in a rare disease cohort.

Long-read HiFi genome sequencing allows for accurate detection and direct phasing of single nucleotide variants, indels, and structural variants. Recent algorithmic development enables simultaneous detection of CpG methylation for analysis of regulatory element activity directly in HiFi reads. We present a comprehensive haplotype resolved 5-base HiFi genome sequencing dataset from a rare disease cohort of 276 samples in 152 families to identify rare (~0.5%) hypermethylation events. We find that 80% of these events are allele-specific and predicted to cause loss of regulatory element activity. We demonstrate heritability of extreme hypermethylation including rare cis variants associated with short (~200 bp) and large hypermethylation events (>1 kb), respectively. We identify repeat expansions in proximal promoters predicting allelic gene silencing via hypermethylation and demonstrate allelic transcriptional events downstream. On average 30-40 rare hypermethylation tiles overlap rare disease genes per patient, providing indications for variation prioritization including a previously undiagnosed pathogenic allele in DIP2B causing global developmental delay. We propose that use of HiFi genome sequencing in unsolved rare disease cases will allow detection of unconventional diseases alleles due to loss of regulatory element activity.

Committing to genomic answers for all kids: Evaluating inequity in genomic research enrollment.

PURPOSE: Persistent inequities in genomic medicine and research contribute to health disparities. This analysis uses a context-specific and equity-focused strategy to evaluate enrollment patterns for Genomic Answers for Kids (GA4K), a large, metropolitan-wide genomic study on children. METHODS: Electronic health records for 2247 GA4K study participants were used to evaluate the distribution of individuals by demographics (race, ethnicity, and payor type) and location (residential address). Addresses were geocoded to produce point density and 3-digit zip code maps showing local and regional enrollment patterns. Health system reports and census data were used to compare participant characteristics with reference populations at different spatial scales. RESULTS: Racial and ethnic minoritized and populations with low-income were underrepresented in the GA4K study cohort. Geographic variation demonstrates inequity in enrollment and participation among children from historically segregated and socially disadvantaged communities. CONCLUSION: Our findings illustrate inequity in enrollment related to both GA4K study design and structural inequalities, which we suspect may exist for similar US-based studies. Our methods provide a scalable framework for continually evaluating and improving study design to ensure equitable participation in and benefits from genomic research and medicine. The use of high-resolution, place-based data represents a novel and practical means of identifying and characterizing inequities and targeting community engagement.

Rare variants in PPFIA3 cause delayed development, intellectual disability, autism, and epilepsy.

PPFIA3 encodes the Protein-Tyrosine Phosphatase, Receptor-Type, F Polypeptide-Interacting Protein Alpha-3 (PPFIA3), which is a member of the LAR protein-tyrosine phosphatase-interacting protein (liprin) family involved in synaptic vesicle transport and presynaptic active zone assembly. The protein structure and function are well conserved in both invertebrates and vertebrates, but human diseases related to PPFIA3 dysfunction are not yet known. Here, we report 14 individuals with rare mono-allelic PPFIA3 variants presenting with features including developmental delay, intellectual disability, hypotonia, autism, and epilepsy. To determine the pathogenicity of PPFIA3 variants in vivo , we generated transgenic fruit flies expressing either human PPFIA3 wildtype (WT) or variant protein using GAL4-UAS targeted gene expression systems. Ubiquitous expression with Actin-GAL4 showed that the PPFIA3 variants had variable penetrance of pupal lethality, eclosion defects, and anatomical leg defects. Neuronal expression with elav-GAL4 showed that the PPFIA3 variants had seizure-like behaviors, motor defects, and bouton loss at the 3 rd instar larval neuromuscular junction (NMJ). Altogether, in the fly overexpression assays, we found that the PPFIA3 variants in the N-terminal coiled coil domain exhibited stronger phenotypes compared to those in the C-terminal region. In the loss-of-function fly assay, we show that the homozygous loss of fly Liprin- α leads to embryonic lethality. This lethality is partially rescued by the expression of human PPFIA3 WT, suggesting human PPFIA3 protein function is partially conserved in the fly. However, the PPFIA3 variants failed to rescue lethality. Altogether, the human and fruit fly data reveal that the rare PPFIA3 variants are dominant negative loss-of-function alleles that perturb multiple developmental processes and synapse formation.

LHX2 haploinsufficiency causes a variable neurodevelopmental disorder.

PURPOSE: LHX2 encodes the LIM homeobox 2 transcription factor (LHX2), which is highly expressed in brain and well conserved across species, but it has not been clearly linked to neurodevelopmental disorders (NDDs) to date. METHODS: Through international collaboration, we identified 19 individuals from 18 families with variable neurodevelopmental phenotypes, carrying a small chromosomal deletion, likely gene-disrupting or missense variants in LHX2. Functional consequences of missense variants were investigated in cellular systems. RESULTS: Affected individuals presented with developmental and/or behavioral abnormalities, autism spectrum disorder, variable intellectual disability, and microcephaly. We observed nucleolar accumulation for 2 missense variants located within the DNA-binding HOX domain, impaired interaction with co-factor LDB1 for another variant located in the protein-protein interaction-mediating LIM domain, and impaired transcriptional activation by luciferase assay for 4 missense variants. CONCLUSION: We implicate LHX2 haploinsufficiency by deletion and likely gene-disrupting variants as causative for a variable NDD. Our findings suggest a loss-of-function mechanism also for likely pathogenic LHX2 missense variants. Together, our observations underscore the importance of LHX2 in the nervous system and for variable neurodevelopmental phenotypes.

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.

Variants in the zinc transporter TMEM163 cause a hypomyelinating leukodystrophy.

Hypomyelinating leukodystrophies comprise a subclass of genetic disorders with deficient myelination of the CNS white matter. Here we report four unrelated families with a hypomyelinating leukodystrophy phenotype harbouring variants in TMEM163 (NM_030923.5). The initial clinical presentation resembled Pelizaeus-Merzbacher disease with congenital nystagmus, hypotonia, delayed global development and neuroimaging findings suggestive of significant and diffuse hypomyelination. Genomic testing identified three distinct heterozygous missense variants in TMEM163 with two unrelated individuals sharing the same de novo variant. TMEM163 is highly expressed in the CNS particularly in newly myelinating oligodendrocytes and was recently revealed to function as a zinc efflux transporter. All the variants identified lie in highly conserved residues in the cytoplasmic domain of the protein, and functional in vitro analysis of the mutant protein demonstrated significant impairment in the ability to efflux zinc out of the cell. Expression of the mutant proteins in an oligodendroglial cell line resulted in substantially reduced mRNA expression of key myelin genes, reduced branching and increased cell death. Our findings indicate that variants in TMEM163 cause a hypomyelinating leukodystrophy and uncover a novel role for zinc homeostasis in oligodendrocyte development and myelin formation.

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.

Phenotypic expansion and variable expressivity in individuals with JARID2-related intellectual disability: A case series.

Loss of function variants in JARID2 were recently reported in 16 patients with a neurodevelopmental disorder characterized by delays, intellectual and learning disability, autism, behavioral abnormalities, and dysmorphic features. Most cases were de novo, with only one variant inherited from an affected parent. Here, we present seven additional individuals from five families with pathogenic or likely pathogenic JARID2 variants, confirming this gene-disease association and highlighting palatal abnormalities and heart defects as part of the phenotype. In addition, we report inheritance of JARID2 variants from mildly affected parents, demonstrating the variable expressivity of the disease. We also note the high prevalence of intragenic JARID2 copy number variants, emphasizing the importance of exon-level analysis.

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.

Heterozygous variants in ZBTB7A cause a neurodevelopmental disorder associated with symptomatic overgrowth of pharyngeal lymphoid tissue, macrocephaly, and elevated fetal hemoglobin.

By clinical whole exome sequencing, we identified 12 individuals with ages 3 to 37 years, including three individuals from the same family, with a consistent phenotype of intellectual disability (ID), macrocephaly, and overgrowth of adenoid tissue. All 12 individuals harbored a rare heterozygous variant in ZBTB7A which encodes the transcription factor Zinc finger and BTB-domain containing protein 7A, known to play a role in lympho- and hematopoiesis. ID was generally mild. Fetal hemoglobin (HbF) fraction was elevated 2.2%-11.2% (reference value <2% in individuals > 6 months) in four of the five individuals for whom results were available. Ten of twelve individuals had undergone surgery at least once for lymphoid hypertrophy limited to the pharynx. In the most severely affected individual (individual 1), airway obstruction resulted in 17 surgical procedures before the age of 13 years. Sleep apnea was present in 8 of 10 individuals. In the nine unrelated individuals, ZBTB7A variants were novel and de novo. The six frameshift/nonsense and four missense variants were spread throughout the gene. This is the first report of a cohort of individuals with this novel syndromic neurodevelopmental disorder.

Deletion of ERF and CIC causes abnormal skull morphology and global developmental delay.

The ETS2 repressor factor (ERF) is a transcription factor in the RAS-MEK-ERK signal transduction cascade that regulates cell proliferation and differentiation, and pathogenic sequence variants in the ERF gene cause variable craniosynostosis inherited in an autosomal dominant pattern. The reported ERF variants are largely loss-of-function, implying haploinsufficiency as a primary disease mechanism; however, ERF gene deletions have not been reported previously. Here we describe three probands with macrocephaly, craniofacial dysmorphology, and global developmental delay. Clinical genetic testing for fragile X and other relevant sequencing panels were negative; however, chromosomal microarray identified heterozygous deletions (63.7-583.2 kb) on Chromosome 19q13.2 in each proband that together included five genes associated with Mendelian diseases (ATP1A3, ERF, CIC, MEGF8, and LIPE). Parental testing indicated that the aberrations were apparently de novo in two of the probands and were inherited in the one proband with the smallest deletion. Deletion of ERF is consistent with the reported loss-of-function ERF variants, prompting clinical copy-number-variant classifications of likely pathogenic. Moreover, the recent characterization of heterozygous loss-of-function CIC sequence variants as a cause of intellectual disability and neurodevelopmental disorders inherited in an autosomal dominant pattern is also consistent with the developmental delays and intellectual disabilities identified among the two probands with CIC deletions. Taken together, this case series adds to the previously reported patients with ERF and/or CIC sequence variants and supports haploinsufficiency of both genes as a mechanism for a variable syndromic cranial phenotype with developmental delays and intellectual disability inherited in an autosomal dominant pattern.

Factors Affecting Migration to GRCh38 in Laboratories Performing Clinical Next-Generation Sequencing.

The most recent build of the human reference genome, GRCh38, was released in 2013. However, many laboratories performing next-generation sequencing (NGS) continue to align to GRCh37. Our aim was to assess the number of clinical diagnostic laboratories that have migrated to GRCh38 and discern factors impeding migration for those still using GRCh37. A brief, five-question survey was electronically administered to 71 clinical laboratories offering constitutional NGS-based testing and analyzed categorically. Twenty-eight responses meeting inclusion criteria were collected from 24 academic and four commercial diagnostic laboratories. Most of these (14; 50%) reported volumes of <500 NGS-based tests in 2019. Only two respondents (7%) had already migrated entirely to GRCh38; most laboratories (15; 54%) had no plans to migrate. The two prevailing reasons for not yet migrating were as follows: laboratories did not feel the benefits outweighed the time and monetary costs (14; 50%); and laboratories had insufficient staff to facilitate the migration (12; 43%). These data, although limited, suggest most clinical molecular laboratories are reluctant to migrate to GRCh38, and there appear to be multiple obstacles to overcome before GRCh38 is widely adopted.

DNA Methylation Signature for EZH2 Functionally Classifies Sequence Variants in Three PRC2 Complex Genes.

Weaver syndrome (WS), an overgrowth/intellectual disability syndrome (OGID), is caused by pathogenic variants in the histone methyltransferase EZH2, which encodes a core component of the Polycomb repressive complex-2 (PRC2). Using genome-wide DNA methylation (DNAm) data for 187 individuals with OGID and 969 control subjects, we show that pathogenic variants in EZH2 generate a highly specific and sensitive DNAm signature reflecting the phenotype of WS. This signature can be used to distinguish loss-of-function from gain-of-function missense variants and to detect somatic mosaicism. We also show that the signature can accurately classify sequence variants in EED and SUZ12, which encode two other core components of PRC2, and predict the presence of pathogenic variants in undiagnosed individuals with OGID. The discovery of a functionally relevant signature with utility for diagnostic classification of sequence variants in EZH2, EED, and SUZ12 supports the emerging paradigm shift for implementation of DNAm signatures into diagnostics and translational research.

Haploinsufficiency of the basic helix-loop-helix transcription factor HAND2 causes congenital heart defects.

Congenital heart defects (CHDs) are caused by a disruption in heart morphogenesis, which is dependent, in part, on a network of transcription factors (TFs) that regulate myocardial development. Heterozygous sequence variants in the basic helix-loop-helix TF gene heart and neural crest derivatives expressed 2 (HAND2) have been reported among some patients with CHDs; however, HAND2 has not yet been established as a Mendelian disease gene. We report a 31-month-old male with unicommissural unicuspid aortic valve, moderate aortic stenosis, and mild pulmonic stenosis. Chromosome analysis revealed a normal 46,XY karyotype, and a CHD sequencing panel was negative for pathogenic variants in NKX2.5, GATA4, TBX5, and CHD7. However, chromosomal microarray (CMA) testing identified a heterozygous 546.0-kb deletion on chromosome 4q34.1 (174364195_174910239[GRCh37/hg19]) that included exons 1 and 2 of SCRG1, HAND2, and HAND2-AS1. Familial CMA testing determined that the deletion was paternally inherited, which supported a likely pathogenic classification as the proband's father had previously undergone surgery for Tetralogy of Fallot. The family history was also notable for a paternal uncle who had previously died from complications related to an unknown heart defect. Taken together, this first report of a HAND2 and HAND2-AS1 deletion in a family with CHDs strongly supports haploinsufficiency of HAND2 as an autosomal dominant cause of CHD.