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

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.

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.

Long-read genome sequencing reveals a novel intronic retroelement insertion in NR5A1 associated with 46,XY differences of sexual development

Despite significant advancements in rare genetic disease diagnostics, many patients with rare genetic disease remain without a molecular diagnosis. Novel tools and methods are needed to improve the detection of disease-associated variants and understand the genetic basis of many rare diseases. Long-read genome sequencing provides improved sequencing in highly repetitive, homologous, and low-complexity regions, and improved assessment of structural variation and complex genomic rearrangements compared to short-read genome sequencing. As such, it is a promising method to explore overlooked genetic variants in rare diseases with a high suspicion of a genetic basis. We therefore applied PacBio HiFi sequencing in a large multi-generational family presenting with autosomal dominant 46,XY differences of sexual development (DSD), for whom extensive molecular testing over multiple decades had failed to identify a molecular diagnosis. This revealed a rare SINE-VNTR-Alu retroelement insertion in intron 4 of NR5A1, a gene in which loss-of-function variants are an established cause of 46,XY DSD. The insertion segregated among affected family members and was associated with loss-of-expression of alleles in cis, demonstrating a functional impact on NR5A1. This case highlights the power of long-read genome sequencing to detect genomic variants that have previously been intractable to detection by standard short-read genomic testing.

Genome and RNA sequencing boost neuromuscular diagnoses to 62% from 34% with exome sequencing alone

OBJECTIVE

Most families with heritable neuromuscular disorders do not receive a molecular diagnosis. Here we evaluate diagnostic utility of exome, genome, RNA sequencing, and protein studies and provide evidence-based recommendations for their integration into practice.

METHODS

In total, 247 families with suspected monogenic neuromuscular disorders who remained without a genetic diagnosis after standard diagnostic investigations underwent research-led massively parallel sequencing: neuromuscular disorder gene panel, exome, genome, and/or RNA sequencing to identify causal variants. Protein and RNA studies were also deployed when required.

RESULTS

Integration of exome sequencing and auxiliary genome, RNA and/or protein studies identified causal or likely causal variants in 62% (152 out of 247) of families. Exome sequencing alone informed 55% (83 out of 152) of diagnoses, with remaining diagnoses (45%; 69 out of 152) requiring genome sequencing, RNA and/or protein studies to identify variants and/or support pathogenicity. Arrestingly, novel disease genes accounted for <4% (6 out of 152) of diagnoses while 36.2% of solved families (55 out of 152) harbored at least one splice-altering or structural variant in a known neuromuscular disorder gene. We posit that contemporary neuromuscular disorder gene-panel sequencing could likely provide 66% (100 out of 152) of our diagnoses today.

INTERPRETATION

Our results emphasize thorough clinical phenotyping to enable deep scrutiny of all rare genetic variation in phenotypically consistent genes. Post-exome auxiliary investigations extended our diagnostic yield by 81% overall (34-62%). We present a diagnostic algorithm that details deployment of genomic and auxiliary investigations to obtain these diagnoses today most effectively. We hope this provides a practical guide for clinicians as they gain greater access to clinical genome and transcriptome sequencing.

Mono-allelic KCNB2 variants lead to a neurodevelopmental syndrome caused by altered channel inactivation

Ion channels mediate voltage fluxes or action potentials that are central to the functioning of excitable cells such as neurons. The KCNB family of voltage-gated potassium channels (Kv) consists of two members (KCNB1 and KCNB2) encoded by KCNB1 and KCNB2, respectively. These channels are major contributors to delayed rectifier potassium currents arising from the neuronal soma which modulate overall excitability of neurons. In this study, we identified several mono-allelic pathogenic missense variants in KCNB2, in individuals with a neurodevelopmental syndrome with epilepsy and autism in some individuals. Recurrent dysmorphisms included a broad forehead, synophrys, and digital anomalies. Additionally, we selected three variants where genetic transmission has not been assessed, from two epilepsy studies, for inclusion in our experiments. We characterized channel properties of these variants by expressing them in oocytes of Xenopus laevis and conducting cut-open oocyte voltage clamp electrophysiology. Our datasets indicate no significant change in absolute conductance and conductance-voltage relationships of most disease variants as compared to wild type (WT), when expressed either alone or co-expressed with WT-KCNB2. However, variants c.1141A>G (p.Thr381Ala) and c.641C>T (p.Thr214Met) show complete abrogation of currents when expressed alone with the former exhibiting a left shift in activation midpoint when expressed alone or with WT-KCNB2. The variants we studied, nevertheless, show collective features of increased inactivation shifted to hyperpolarized potentials. We suggest that the effects of the variants on channel inactivation result in hyper-excitability of neurons, which contributes to disease manifestations.

Severus: accurate detection and characterization of somatic structural variation in tumor genomes using long reads

Most current studies rely on short-read sequencing to detect somatic structural variation (SV) in cancer genomes. Long-read sequencing offers the advantage of better mappability and long-range phasing, which results in substantial improvements in germline SV detection. However, current long-read SV detection methods do not generalize well to the analysis of somatic SVs in tumor genomes with complex rearrangements, heterogeneity, and aneuploidy. Here, we present Severus: a method for the accurate detection of different types of somatic SVs using a phased breakpoint graph approach. To benchmark various short- and long-read SV detection methods, we sequenced five tumor/normal cell line pairs with Illumina, Nanopore, and PacBio sequencing platforms; on this benchmark Severus showed the highest F1 scores (harmonic mean of the precision and recall) as compared to long-read and short-read methods. We then applied Severus to three clinical cases of pediatric cancer, demonstrating concordance with known genetic findings as well as revealing clinically relevant cryptic rearrangements missed by standard genomic panels.

Mapping structural variants to rare disease genes using long-read whole genome sequencing and trait-relevant polygenic scores

Recent studies have revealed the pervasive landscape of rare structural variants (rSVs) present in human genomes. rSVs can have extreme effects on the expression of proximal genes and, in a rare disease context, have been implicated in patient cases where no diagnostic single nucleotide variant (SNV) was found. Approaches for integrating rSVs to date have focused on targeted approaches in known Mendelian rare disease genes. This approach is intractable for rare diseases with many causal loci or patients with complex, multi-phenotype syndromes. We hypothesized that integrating trait-relevant polygenic scores (PGS) would provide a substantial reduction in the number of candidate disease genes in which to assess rSV effects. We further implemented a method for ranking PGS genes to define a set of core/key genes where a rSV has the potential to exert relatively larger effects on disease risk. Among a subset of patients enrolled in the Genomic Answers for Kids (GA4K) rare disease program (N=497), we used PacBio HiFi long-read whole genome sequencing (lrWGS) to identify rSVs intersecting genes in trait-relevant PGSs. Illustrating our approach in Autism (N=54 cases), we identified 1,827 deletions, 158 duplications, 619 insertions, and 14 inversions overlapping putative core/key PGS genes. Additionally, by integrating genomic constraint annotations from gnomAD, we observed that rare duplications overlapping putative core/key PGS genes were frequently in higher constraint regions compared to controls (P = 2×10-04). This difference was not observed in the lowest-ranked gene set (P = 0.18). Overall, our study provides a framework for the annotation of long-read rSVs from lrWGS data and prioritization of disease-linked genomic regions for downstream functional validation of rSV impacts. To enable reuse by other researchers, we have made SV allele frequencies and gene associations freely available.

Exome and genome sequencing in a heterogeneous population of patients with rare disease: Identifying predictors of a diagnosis

PURPOSE

Exome (ES) and genome sequencing (GS) are increasingly being utilized for individuals with rare and undiagnosed diseases; however, guidelines on their use remain limited. This study aimed to identify factors associated with diagnosis by ES and/or GS in a heterogeneous population of patients with rare and undiagnosed diseases.

METHODS

In this case control study, we reviewed data from 400 diagnosed and 400 undiagnosed randomly selected participants in the Undiagnosed Diseases Network, all of whom had undergone ES and/or GS. We analyzed factors associated with receiving a diagnosis by ES and/or GS.

RESULTS

Factors associated with a decreased odds of being diagnosed included adult symptom onset, singleton sequencing, and having undergone ES and/or GS before acceptance to the Undiagnosed Diseases Network (48%, 51%, and 32% lower odds, respectively). Factors that increased the odds of being diagnosed by ES and/or GS included having primarily neurological symptoms and having undergone prior chromosomal microarray testing (44% and 59% higher odds, respectively).

CONCLUSION

We identified several factors that were associated with receiving a diagnosis by ES and/or GS. This will ideally inform the utilization of ES and/or GS and help manage expectations of individuals and families undergoing these tests.

A 25-year odyssey of genomic technology advances and structural variant discovery

This perspective focuses on advances in genome technology over the last 25 years and their impact on germline variant discovery within the field of human genetics. The field has witnessed tremendous technological advances from microarrays to short-read sequencing and now long-read sequencing. Each technology has provided genome-wide access to different classes of human genetic variation. We are now on the verge of comprehensive variant detection of all forms of variation for the first time with a single assay. We predict that this transition will further transform our understanding of human health and biology and, more importantly, provide novel insights into the dynamic mutational processes shaping our genomes.

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.

Referral networks for pediatric patients with genetic conditions: The perspective of occupational therapists

Families of children with developmental delays but no diagnosed genetic condition may benefit from connection to genetic systems of care. This work examines the role of occupational therapy as a space for families of pediatric patients to gain access to genetic services. Between September 2021 and February 2022, we interviewed 20 occupational therapists in New England who work primarily with pediatric patients. We transcribed the interviews and used a grounded theory approach to identify and code recurring themes. The data reveal several barriers to linking pediatric patients to genetic systems of care, including lack of insurance coverage, wait times for appointments and test results, hesitant primary care providers, and familial and cultural stigma of disability. We discuss the unique role of occupational therapists as professionals who spend substantial time with patients, often in their everyday environments, to bridge these barriers. We also address challenges associated with occupational therapists facilitating connections to genetics services, including their lack of specialized knowledge of genetics and barriers fully integrating with others on the medical team.

Panels, Exomes, Genomes, and More-Finding the Best Path Through the Diagnostic Odyssey

Selecting the ideal test to evaluate an individual with a suspected genetic disorder can be challenging. While several clinical testing options are available, no single test yet captures all potentially causative genetic variants. Thus, clinicians may order testing in a stepwise fashion, and what to order after non-diagnostic testing can be challenging to determine. Here, we provide an overview of commonly used clinical genetic tests, guidance on when they are best used, and what they may miss. We conclude with a discussion of how new technologies might be used to identify challenging variants and simplify clinical testing in the future.

Synchronized long-read genome, methylome, epigenome, and transcriptome for resolving a Mendelian condition

Resolving the molecular basis of a Mendelian condition (MC) remains challenging owing to the diverse mechanisms by which genetic variants cause disease. To address this, we developed a synchronized long-read genome, methylome, epigenome, and transcriptome sequencing approach, which enables accurate single-nucleotide, insertion-deletion, and structural variant calling and diploid de novo genome assembly, and permits the simultaneous elucidation of haplotype-resolved CpG methylation, chromatin accessibility, and full-length transcript information in a single long-read sequencing run. Application of this approach to an Undiagnosed Diseases Network (UDN) participant with a chromosome X;13 balanced translocation of uncertain significance revealed that this translocation disrupted the functioning of four separate genes (NBEA, PDK3, MAB21L1, and RB1) previously associated with single-gene MCs. Notably, the function of each gene was disrupted via a distinct mechanism that required integration of the four 'omes' to resolve. These included nonsense-mediated decay, fusion transcript formation, enhancer adoption, transcriptional readthrough silencing, and inappropriate X chromosome inactivation of autosomal genes. Overall, this highlights the utility of synchronized long-read multi-omic profiling for mechanistically resolving complex phenotypes.

Case report: a novel deep intronic splice-altering variant in DMD as a cause of Becker muscular dystrophy

We present the case of a male patient who was ultimately diagnosed with Becker muscular dystrophy (BMD; MIM# 300376) after the onset of muscle weakness in his teens progressively led to significant walking difficulties in his twenties. A genetic diagnosis was pursued but initial investigation revealed no aberrations in the dystrophin gene (DMD), although immunohistochemistry and Western blot analysis suggested the diagnosis of dystrophinopathy. Eventually, after more than 10 years, an RNA analysis captured abnormal splicing where 154 nucleotides from intron 43 were inserted between exon 43 and 44 resulting in a frameshift and a premature stop codon. Normal splicing of the DMD gene was also observed. Additionally, a novel variant c.6291-13537A>G in DMD was confirmed in the genomic DNA of the patient. The predicted function of the variant aligns with the mRNA results. To conclude, we here demonstrate that mRNA analysis can guide the diagnosis of non-coding genetic variants in DMD.

Local read haplotagging enables accurate long-read small variant calling

Long-read sequencing technology has enabled variant detection in difficult-to-map regions of the genome and enabled rapid genetic diagnosis in clinical settings. Rapidly evolving third-generation sequencing platforms like Pacific Biosciences (PacBio) and Oxford nanopore technologies (ONT) are introducing newer platforms and data types. It has been demonstrated that variant calling methods based on deep neural networks can use local haplotyping information with long-reads to improve the genotyping accuracy. However, using local haplotype information creates an overhead as variant calling needs to be performed multiple times which ultimately makes it difficult to extend to new data types and platforms as they get introduced. In this work, we have developed a local haplotype approximate method that enables state-of-the-art variant calling performance with multiple sequencing platforms including PacBio Revio system, ONT R10.4 simplex and duplex data. This addition of local haplotype approximation makes DeepVariant a universal variant calling solution for long-read sequencing platforms.

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.

Long-Read DNA Sequencing: Recent Advances and Remaining Challenges

DNA sequencing has revolutionized medicine over recent decades. However, analysis of large structural variation and repetitive DNA, a hallmark of human genomes, has been limited by short-read technology, with read lengths of 100-300 bp. Long-read sequencing (LRS) permits routine sequencing of human DNA fragments tens to hundreds of kilobase pairs in size, using both real-time sequencing by synthesis and nanopore-based direct electronic sequencing. LRS permits analysis of large structural variation and haplotypic phasing in human genomes and has enabled the discovery and characterization of rare pathogenic structural variants and repeat expansions. It has also recently enabled the assembly of a complete, gapless human genome that includes previously intractable regions, such as highly repetitive centromeres and homologous acrocentric short arms. With the addition of protocols for targeted enrichment, direct epigenetic DNA modification detection, and long-range chromatin profiling, LRS promises to launch a new era of understanding of genetic diversity and pathogenic mutations in human populations.

Case report: Severe nonketotic hyperglycinemia in a neonate without apparent seizures but concomitant cleft palate and cerebral sinovenous thrombosis

Nonketotic hyperglycinemia (NKH) is in most cases a fatal inborn error of metabolism which usually presents during the neonatal period as encephalopathy and refractory seizures. The reported congenital anomalies associated with NKH included corpus callosal agenesis, club foot, cleft palate, and congenital heart disease. Here, we report a newborn who presented with encephalopathy without overt seizures, cerebral venous sinus thrombosis, and cleft palate. Electroencephalography showed a burst suppression pattern, which suggests the etiology could be due to a metabolic or genetic disorder. The amino acid analysis of plasma and cerebrospinal fluid showed elevated glycine. Whole exome sequencing identified a heterozygous c.492C > G; p.Tyr164Ter variant in exon 4 of the GLDC gene inherited from the patient's father. Further long-read whole genome sequencing revealed an exon 1-2 deletion in the GLDC gene inherited from the patient's mother. Additional analyses revealed no pathogenic variants of the cleft palate-related genes. The cleft palate may be an associated congenital anomaly in NKH. Regarding cerebral venous sinus thrombosis, we found a heterozygous variant (p.Arg189Trp) of the PROC gene, which is a common cause of thrombophilia among Thai newborns. A neonate with NKH could present with severe encephalopathy without seizures. A close follow up for clinical changes and further next generation sequencing are crucial for definite diagnosis in neonates with encephalopathy of unclear cause.

Beyond the exome: What's next in diagnostic testing for Mendelian conditions

Despite advances in clinical genetic testing, including the introduction of exome sequencing (ES), more than 50% of individuals with a suspected Mendelian condition lack a precise molecular diagnosis. Clinical evaluation is increasingly undertaken by specialists outside of clinical genetics, often occurring in a tiered fashion and typically ending after ES. The current diagnostic rate reflects multiple factors, including technical limitations, incomplete understanding of variant pathogenicity, missing genotype-phenotype associations, complex gene-environment interactions, and reporting differences between clinical labs. Maintaining a clear understanding of the rapidly evolving landscape of diagnostic tests beyond ES, and their limitations, presents a challenge for non-genetics professionals. Newer tests, such as short-read genome or RNA sequencing, can be challenging to order, and emerging technologies, such as optical genome mapping and long-read DNA sequencing, are not available clinically. Furthermore, there is no clear guidance on the next best steps after inconclusive evaluation. Here, we review why a clinical genetic evaluation may be negative, discuss questions to be asked in this setting, and provide a framework for further investigation, including the advantages and disadvantages of new approaches that are nascent in the clinical sphere. We present a guide for the next best steps after inconclusive molecular testing based upon phenotype and prior evaluation, including when to consider referral to research consortia focused on elucidating the underlying cause of rare unsolved genetic disorders.

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.

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.

Solving inherited white matter disorder etiologies in the neurology clinic: Challenges and lessons learned using next-generation sequencing

Introduction

Rare neurodevelopmental disorders, including inherited white matter disorders or leukodystrophies, often present a diagnostic challenge on a genetic level given the large number of causal genes associated with a range of disease subtypes. This study aims to demonstrate the challenges and lessons learned in the genetic investigations of leukodystrophies through presentation of a series of cases solved using exome or genome sequencing.

Methods

Each of the six patients had a leukodystrophy associated with hypomyelination or delayed myelination on MRI, and inconclusive clinical diagnostic genetic testing results. We performed next generation sequencing (case-based exome or genome sequencing) to further investigate the genetic cause of disease.

Results

Following different lines of investigation, molecular diagnoses were obtained for each case, with patients harboring pathogenic variants in a range of genes including TMEM106B, GJA1, AGA, POLR3A, and TUBB4A. We describe the lessons learned in reaching the genetic diagnosis, including the importance of (a) utilizing proper multi-gene panels in clinical testing, (b) assessing the reliability of biochemical assays in supporting diagnoses, and (c) understanding the limitations of exome sequencing methods in regard to CNV detection and region coverage in GC-rich areas.

Discussion

This study illustrates the importance of applying a collaborative diagnostic approach by combining detailed phenotyping data and metabolic results from the clinical environment with advanced next generation sequencing analysis techniques from the research environment to increase the diagnostic yield in patients with genetically unresolved leukodystrophies.

Exome/Genome Sequencing in Undiagnosed Syndromes

Exome sequencing (ES) and genome sequencing (GS) have radically transformed the diagnostic approach to undiagnosed rare/ultrarare Mendelian diseases. Next-generation sequencing (NGS), the technology integral for ES, GS, and most large (100+) gene panels, has enabled previously unimaginable diagnoses, changes in medical management, new treatments, and accurate reproductive risk assessments for patients, as well as new disease gene discoveries. Yet, challenges remain, as most individuals remain undiagnosed with current NGS. Improved NGS technology has resulted in long-read sequencing, which may resolve diagnoses in some patients who do not obtain a diagnosis with current short-read ES and GS, but its effectiveness is unclear, and it is expensive. Other challenges that persist include the resolution of variants of uncertain significance, the urgent need for patients with ultrarare disorders to have access to therapeutics, the need for equity in patient access to NGS-based testing, and the study of ethical concerns. However, the outlook for undiagnosed disease resolution is bright, due to continual advancements in the field.

Long-read HiFi sequencing of NUDT15: Phased full-gene haplotyping and pharmacogenomic allele discovery

To determine the phase of NUDT15 sequence variants for more comprehensive star (*) allele diplotyping, we developed a novel long-read single-molecule real-time HiFi amplicon sequencing method. A 10.5 kb NUDT15 amplicon assay was validated using reference material positive controls and additional samples for specimen type and blinded accuracy assessment. Triplicate NUDT15 HiFi sequencing of two reference material samples had nonreference genotype concordances of >99.9%, indicating that the assay is robust. Notably, short-read genome sequencing of a subset of samples was unable to determine the phase of star (*) allele-defining NUDT15 variants, resulting in ambiguous diplotype results. In contrast, long-read HiFi sequencing phased all variants across the NUDT15 amplicons, including a *2/*9 diplotype that previously was characterized as *1/*2 in the 1000 Genomes Project v3 data set. Assay throughput was also tested using 8.5 kb amplicons from 100 Ashkenazi Jewish individuals, which identified a novel NUDT15 *1 suballele (c.-121G>A) and a rare likely deleterious coding variant (p.Pro129Arg). Both novel alleles were Sanger confirmed and assigned as *1.007 and *20, respectively, by the PharmVar Consortium. Taken together, NUDT15 HiFi amplicon sequencing is an innovative method for phased full-gene characterization and novel allele discovery, which could improve NUDT15 pharmacogenomic testing and subsequent phenotype prediction.

Long-read sequencing for molecular diagnostics in constitutional genetic disorders

Long-read sequencing (LRS) has been around for more than a decade, but widespread adoption of the technology has been slow due to the perceived high error rates and high sequencing cost. This is changing due to the recent advancements to produce highly accurate sequences and the reducing costs. LRS promises significant improvement over short read sequencing in four major areas: (1) better detection of structural variation (2) better resolution of highly repetitive or nonunique regions (3) accurate long-range haplotype phasing and (4) the detection of base modifications natively from the sequencing data. Several successful applications of LRS have demonstrated its ability to resolve molecular diagnoses where short-read sequencing fails to identify a cause. However, the argument for increased diagnostic yield from LRS remains to be validated. Larger cohort studies may be required to establish the realistic boundaries of LRS's clinical utility and analytical validity, as well as the development of standards for clinical applications. We discuss the limitations of the current standard of care, and contrast with the applications and advantages of two major LRS platforms, PacBio and Oxford Nanopore, for molecular diagnostics of constitutional disorders, and present a critical argument about the potential of LRS in diagnostic settings.

Approaches to long-read sequencing in a clinical setting to improve diagnostic rate

Over the past decade, advances in genetic testing, particularly the advent of next-generation sequencing, have led to a paradigm shift in the diagnosis of molecular diseases and disorders. Despite our present collective ability to interrogate more than 90% of the human genome, portions of the genome have eluded us, resulting in stagnation of diagnostic yield with existing methodologies. Here we show how application of a new technology, long-read sequencing, has the potential to improve molecular diagnostic rates. Whole genome sequencing by long reads was able to cover 98% of next-generation sequencing dead zones, which are areas of the genome that are not interpretable by conventional industry-standard short-read sequencing. Through the ability of long-read sequencing to unambiguously call variants in these regions, we discovered an immunodeficiency due to a variant in IKBKG in a subject who had previously received a negative genome sequencing result. Additionally, we demonstrate the ability of long-read sequencing to detect small variants on par with short-read sequencing, its superior performance in identifying structural variants, and thirdly, its capacity to determine genomic methylation defects in native DNA. Though the latter technical abilities have been demonstrated, we demonstrate the clinical application of this technology to successfully identify multiple types of variants using a single test.

Comprehensive cancer predisposition testing within the prospective MASTER trial identifies hereditary cancer patients and supports treatment decisions for rare cancers

BACKGROUND

Germline variant evaluation in precision oncology opens new paths towards the identification of patients with genetic tumor risk syndromes and the exploration of therapeutic relevance. Here, we present the results of germline variant analysis and their clinical implications in a precision oncology study for patients with predominantly rare cancers.

PATIENTS AND METHODS

Matched tumor and control genome/exome and RNA sequencing was performed for 1,485 patients with rare cancers (79%) and/or young adults (77% younger than 51 years) in the NCT/DKTK MASTER trial, a German multicenter, prospective observational precision oncology study. Clinical and therapeutic relevance of prospective pathogenic germline variant (PGV) evaluation was analyzed and compared to other precision oncology studies.

RESULTS

Ten percent of patients (n=157) harbored PGVs in 35 genes associated with autosomal dominant cancer predisposition, whereof up to 75% were unknown before study participation. Another five percent of patients (n=75) were heterozygous carriers for recessive genetic tumor risk syndromes. Particularly high PGV yields were found in patients with gastrointestinal stromal tumors (GISTs) (28%, 11/40), and more specific in wild-type GISTS (50%, n=10/20), leiomyosarcomas (21%, n=19/89), and hepatopancreaticobiliary cancers (16%, n=16/97). Forty-five percent of PGVs (n=100/221) supported treatment recommendations, and its implementation led to a clinical benefit in 40% of patients (n=10/25). A comparison of different precision oncology studies revealed variable PGV yields and considerable differences in germline variant analysis workflows. We therefore propose a detailed workflow for germline variant evaluation.

CONCLUSIONS

Genetic germline testing in patients with rare cancers can identify the very first patient in a hereditary cancer family and can lead to clinical benefit in a broad range of entities. Its routine implementation in precision oncology accompanied by the harmonization of germline variant evaluation workflows will increase clinical benefit and boost research.

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.

Massively parallel identification of functionally consequential noncoding genetic variants in undiagnosed rare disease patients

Clinical whole genome sequencing has enabled the discovery of potentially pathogenic noncoding variants in the genomes of rare disease patients with a prior history of negative genetic testing. However, interpreting the functional consequences of noncoding variants and distinguishing those that contribute to disease etiology remains a challenge. Here we address this challenge by experimentally profiling the functional consequences of rare noncoding variants detected in a cohort of undiagnosed rare disease patients at scale using a massively parallel reporter assay. We demonstrate that this approach successfully identifies rare noncoding variants that alter the regulatory capacity of genomic sequences. In addition, we describe an integrative analysis that utilizes genomic features alongside patient clinical data to further prioritize candidate variants with an increased likelihood of pathogenicity. This work represents an important step towards establishing a framework for the functional interpretation of clinically detected noncoding variants.