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.

Biallelic BORCS8 variants cause an infantile-onset neurodegenerative disorder with altered lysosome dynamics.

BLOC-one-related complex (BORC) is a multiprotein complex composed of eight subunits named BORCS1-8. BORC associates with the cytosolic face of lysosomes, where it sequentially recruits the small GTPase ARL8 and kinesin-1 and -3 microtubule motors to promote anterograde transport of lysosomes toward the peripheral cytoplasm in non-neuronal cells and the distal axon in neurons. The physiological and pathological importance of BORC in humans, however, remains to be determined. Here, we report the identification of compound heterozygous variants [missense c.85T>C (p.Ser29Pro) and frameshift c.71-75dupTGGCC (p.Asn26Trpfs*51)] and homozygous variants [missense c.196A>C (p.Thr66Pro) and c.124T>C (p.Ser42Pro)] in BORCS8 in five children with a severe early-infantile neurodegenerative disorder from three unrelated families. The children exhibit global developmental delay, severe-to-profound intellectual disability, hypotonia, limb spasticity, muscle wasting, dysmorphic facies, optic atrophy, leuko-axonopathy with hypomyelination, and neurodegenerative features with prevalent supratentorial involvement. Cellular studies using a heterologous transfection system show that the BORCS8 missense variants p.Ser29Pro, p.Ser42Pro and p.Thr66Pro are expressed at normal levels but exhibit reduced assembly with other BORC subunits and reduced ability to drive lysosome distribution toward the cell periphery. The BORCS8 frameshift variant p.Asn26Trpfs*51, on the other hand, is expressed at lower levels and is completely incapable of assembling with other BORC subunits and promoting lysosome distribution toward the cell periphery. Therefore, all the BORCS8 variants are partial or total loss-of-function alleles and are thus likely pathogenic. Knockout of the orthologous borcs8 in zebrafish causes decreased brain and eye size, neuromuscular anomalies and impaired locomotion, recapitulating some of the key traits of the human disease. These findings thus identify BORCS8 as a novel genetic locus for an early-infantile neurodegenerative disorder and highlight the critical importance of BORC and lysosome dynamics for the development and function of the central nervous system.

Biallelic Loss of Function Variants in SENP7 Cause Immunodeficiency with Neurologic and Muscular Phenotypes.

To evaluate a novel candidate disease gene, we engaged international collaborators and identified rare, biallelic, specifically homozygous, loss of function variants in SENP7 in 4 children from 3 unrelated families presenting with neurodevelopmental abnormalities, dysmorphism, and immunodeficiency. Their clinical presentations were characterized by hypogammaglobulinemia, intermittent neutropenia, and ultimately death in infancy for all 4 patients. SENP7 is a sentrin-specific protease involved in posttranslational modification of proteins essential for cell regulation, via a process referred to as deSUMOylation. We propose that deficiency of deSUMOylation may represent a novel mechanism of primary immunodeficiency.

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.

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.

Characterization and visualization of tandem repeats at genome scale.

Tandem repeat (TR) variation is associated with gene expression changes and numerous rare monogenic diseases. Although long-read sequencing provides accurate full-length sequences and methylation of TRs, there is still a need for computational methods to profile TRs across the genome. Here we introduce the Tandem Repeat Genotyping Tool (TRGT) and an accompanying TR database. TRGT determines the consensus sequences and methylation levels of specified TRs from PacBio HiFi sequencing data. It also reports reads that support each repeat allele. These reads can be subsequently visualized with a companion TR visualization tool. Assessing 937,122 TRs, TRGT showed a Mendelian concordance of 98.38%, allowing a single repeat unit difference. In six samples with known repeat expansions, TRGT detected all expansions while also identifying methylation signals and mosaicism and providing finer repeat length resolution than existing methods. Additionally, we released a database with allele sequences and methylation levels for 937,122 TRs across 100 genomes.

Large-scale mutational analysis identifies UNC93B1 variants that drive TLR-mediated autoimmunity in mice and humans.

Nucleic acid-sensing Toll-like receptors (TLR) 3, 7/8, and 9 are key innate immune sensors whose activities must be tightly regulated to prevent systemic autoimmune or autoinflammatory disease or virus-associated immunopathology. Here, we report a systematic scanning-alanine mutagenesis screen of all cytosolic and luminal residues of the TLR chaperone protein UNC93B1, which identified both negative and positive regulatory regions affecting TLR3, TLR7, and TLR9 responses. We subsequently identified two families harboring heterozygous coding mutations in UNC93B1, UNC93B1+/T93I and UNC93B1+/R336C, both in key negative regulatory regions identified in our screen. These patients presented with cutaneous tumid lupus and juvenile idiopathic arthritis plus neuroinflammatory disease, respectively. Disruption of UNC93B1-mediated regulation by these mutations led to enhanced TLR7/8 responses, and both variants resulted in systemic autoimmune or inflammatory disease when introduced into mice via genome editing. Altogether, our results implicate the UNC93B1-TLR7/8 axis in human monogenic autoimmune diseases and provide a functional resource to assess the impact of yet-to-be-reported UNC93B1 mutations.