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.

MARK2 variants cause autism spectrum disorder via the downregulation of WNT/ß-catenin signaling pathway.

Microtubule affinity-regulating kinase 2 (MARK2) contributes to establishing neuronal polarity and developing dendritic spines. Although large-scale sequencing studies have associated MARK2 variants with autism spectrum disorder (ASD), the clinical features and variant spectrum in affected individuals with MARK2 variants, early developmental phenotypes in mutant human neurons, and the pathogenic mechanism underlying effects on neuronal development have remained unclear. Here, we report 31 individuals with MARK2 variants and presenting with ASD, other neurodevelopmental disorders, and distinctive facial features. Loss-of-function (LoF) variants predominate (81%) in affected individuals, while computational analysis and in vitro expression assay of missense variants supported the effect of MARK2 loss. Using proband-derived and CRISPR-engineered isogenic induced pluripotent stem cells (iPSCs), we show that MARK2 loss leads to early neuronal developmental and functional deficits, including anomalous polarity and dis-organization in neural rosettes, as well as imbalanced proliferation and differentiation in neural progenitor cells (NPCs). Mark2+/- mice showed abnormal cortical formation and partition and ASD-like behavior. Through the use of RNA sequencing (RNA-seq) and lithium treatment, we link MARK2 loss to downregulation of the WNT/ß-catenin signaling pathway and identify lithium as a potential drug for treating MARK2-associated ASD.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Rare SUZ12 variants commonly cause an overgrowth phenotype.

The Polycomb repressive complex 2 is an epigenetic writer and recruiter with a role in transcriptional silencing. Constitutional pathogenic variants in its component proteins have been found to cause two established overgrowth syndromes: Weaver syndrome (EZH2-related overgrowth) and Cohen-Gibson syndrome (EED-related overgrowth). Imagawa et al. (2017) initially reported a singleton female with a Weaver-like phenotype with a rare coding SUZ12 variant-the same group subsequently reported two additional affected patients. Here we describe a further 10 patients (from nine families) with rare heterozygous SUZ12 variants who present with a Weaver-like phenotype. We report four frameshift, two missense, one nonsense, and two splice site variants. The affected patients demonstrate variable pre- and postnatal overgrowth, dysmorphic features, musculoskeletal abnormalities and developmental delay/intellectual disability. Some patients have genitourinary and structural brain abnormalities, and there may be an association with respiratory issues. The addition of these 10 patients makes a compelling argument that rare pathogenic SUZ12 variants frequently cause overgrowth, physical abnormalities, and abnormal neurodevelopmental outcomes in the heterozygous state. Pathogenic SUZ12 variants may be de novo or inherited, and are sometimes inherited from a mildly-affected parent. Larger samples sizes will be needed to elucidate whether one or more clinically-recognizable syndromes emerge from different variant subtypes.

Weaver Syndrome-Associated EZH2 Protein Variants Show Impaired Histone Methyltransferase Function In Vitro.

Weaver syndrome (WS) is a rare congenital disorder characterized by generalized overgrowth, macrocephaly, specific facial features, accelerated bone age, intellectual disability, and susceptibility to cancers. De novo mutations in the enhancer of zeste homolog 2 (EZH2) have been shown to cause WS. EZH2 is a histone methyltransferase that acts as the catalytic agent of the polycomb-repressive complex 2 (PRC2) to maintain gene repression via methylation of lysine 27 on histone H3 (H3K27). Functional studies investigating histone methyltransferase activity of mutant EZH2 from various cancers have been reported, whereas WS-associated mutations remain poorly characterized. To investigate the role of EZH2 in WS, we performed functional studies using artificially assembled PRC2 complexes containing mutagenized human EZH2 that reflected the codon changes predicted from patients with WS. We found that WS-associated amino acid alterations reduce the histone methyltransferase function of EZH2 in this in vitro assay. Our results support the hypothesis that WS is caused by constitutional mutations in EZH2 that alter the histone methyltransferase function of PRC2. However, histone methyltransferase activities of different EZH2 variants do not appear to correlate directly with the phenotypic variability between WS patients and individuals with a common c.553G>C (p.Asp185His) polymorphism in EZH2.

A novel mutation in EED associated with overgrowth.

In a patient suspected clinically to have Weaver syndrome, we ruled out mutations in EZH2 and NSD1, then identified a previously undescribed de novo mutation in EZH2's partner protein EED. Both proteins are members of the Polycomb Repressive Complex 2 that maintains gene silencing. On the basis of the similarities of the patient's phenotype to Weaver syndrome, which is caused by de novo mutations in EZH2, and on other lines of evidence including mouse Eed hypomorphs, we characterize this mutation as probably pathogenic for a Weaver-like overgrowth syndrome. This is the first report of overgrowth and related phenotypes associated with a constitutional mutation in human EED.

Somatic mosaicism for the p.His1047Arg mutation in PIK3CA in a girl with mesenteric lipomatosis.

We describe a patient who presented with a localized growth of mature fat tissue, which was surgically removed. MRI imaging identified diffuse increase in visceral adipose tissue. Targeted deep sequencing of the resected tissue uncovered a p.H1047R variant in PIK3CA, which was absent in blood. This report expands the phenotypic spectrum of mosaic PIK3CA mutations.