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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

CYP2C8, CYP2C9, and CYP2C19 Characterization Using Next-Generation Sequencing and Haplotype Analysis: A GeT-RM Collaborative Project.

Pharmacogenetic tests typically target selected sequence variants to identify haplotypes that are often defined by star (∗) allele nomenclature. Due to their design, these targeted genotyping assays are unable to detect novel variants that may change the function of the gene product and thereby affect phenotype prediction and patient care. In the current study, 137 DNA samples that were previously characterized by the Genetic Testing Reference Material (GeT-RM) program using a variety of targeted genotyping methods were recharacterized using targeted and whole genome sequencing analysis. Sequence data were analyzed using three genotype calling tools to identify star allele diplotypes for CYP2C8, CYP2C9, and CYP2C19. The genotype calls from next-generation sequencing (NGS) correlated well to those previously reported, except when novel alleles were present in a sample. Six novel alleles and 38 novel suballeles were identified in the three genes due to identification of variants not covered by targeted genotyping assays. In addition, several ambiguous genotype calls from a previous study were resolved using the NGS and/or long-read NGS data. Diplotype calls were mostly consistent between the calling algorithms, although several discrepancies were noted. This study highlights the utility of NGS for pharmacogenetic testing and demonstrates that there are many novel alleles that are yet to be discovered, even in highly characterized genes such as CYP2C9 and CYP2C19.

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.

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.

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.

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.

The NSIGHT1-randomized controlled trial: rapid whole-genome sequencing for accelerated etiologic diagnosis in critically ill infants.

Genetic disorders are a leading cause of morbidity and mortality in infants in neonatal and pediatric intensive care units (NICU/PICU). While genomic sequencing is useful for genetic disease diagnosis, results are usually reported too late to guide inpatient management. We performed an investigator-initiated, partially blinded, pragmatic, randomized, controlled trial to test the hypothesis that rapid whole-genome sequencing (rWGS) increased the proportion of NICU/PICU infants receiving a genetic diagnosis within 28 days. The participants were families with infants aged <4 months in a regional NICU and PICU, with illnesses of unknown etiology. The intervention was trio rWGS. Enrollment from October 2014 to June 2016, and follow-up until November 2016. Of all, 26 female infants, 37 male infants, and 2 infants of undetermined sex were randomized to receive rWGS plus standard genetic tests (n = 32, cases) or standard genetic tests alone (n = 33, controls). The study was terminated early due to loss of equipoise: 73% (24) controls received genomic sequencing as standard tests, and 15% (five) controls underwent compassionate cross-over to receive rWGS. Nevertheless, intention to treat analysis showed the rate of genetic diagnosis within 28 days of enrollment (the primary end-point) to be higher in cases (31%, 10 of 32) than controls (3%, 1 of 33; difference, 28% [95% CI, 10-46%]; p = 0.003). Among infants enrolled in the first 25 days of life, the rate of neonatal diagnosis was higher in cases (32%, 7 of 22) than controls (0%, 0 of 23; difference, 32% [95% CI, 11-53%];p = 0.004). Median age at diagnosis (25 days [range 14-90] in cases vs. 130 days [range 37-451] in controls) and median time to diagnosis (13 days [range 1-84] in cases, vs. 107 days [range 21-429] in controls) were significantly less in cases than controls (p = 0.04). In conclusion, rWGS increased the proportion of NICU/PICU infants who received timely diagnoses of genetic diseases.

Erratum: Constellation: a tool for rapid, automated phenotype assignment of a highly polymorphic pharmacogene, CYP2D6, from whole-genome sequences.

[This corrects the article DOI: 10.1038/npjgenmed.2015.7.].