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

Splicing Defects and Cell Death Cause SF3B2-Linked Craniofacial Microsomia

Craniofacial microsomia (CFM) is a genetically and phenotypically heterogeneous disorder characterized by hypoplasia of facial tissue that is often asymmetric. Affected tissues typically include the ears (external and internal), mandible, and maxilla, but various extracranial anomalies have also been reported. Loss-of-function variants in the SF3B2 gene have recently been reported in 8 cases of CFM, representing one of the more common genetic causes identified to date. To better define the full phenotypic spectrum associated with variants in SF3B2, we report novel loss-of-function variants in SF3B2 in 5 new families with CFM. Furthermore, to determine the mechanism by which SF3B2 loss-of-function perturbs craniofacial development, we established sf3b2-null mutant zebrafish, which exhibited severe deficiencies in craniofacial cartilage and bone progenitors due to elevated apoptosis and reduced proliferation of cranial neural crest cells. In addition, we generated a heterozygous truncating variant of SF3B2 in human induced pluripotent stem cells using CRISPR/Cas9 gene editing. Differentiation of these cells into neural crest cells was accompanied by increased cell death and reduced proliferation. RNA sequencing of sf3b2 mutant zebrafish revealed widespread disruption of mRNA splicing, including mdm2, a key regulator of Tp53-mediated apoptosis. Genetic inhibition of tp53 in sf3b2 mutants demonstrated that tp53 inhibition reduces early cell death but does not improve proliferation or craniofacial cartilage development. Therefore, our functional studies indicate that widespread mRNA splicing disruption, in addition to Tp53-dependent cell death, contributes to the craniofacial features observed in SF3B2-related CFM.

Expanded methylome and quantitative trait loci detection by long-read profiling of personal DNA

Structural variants (SVs) are omnipresent in human DNA, yet their genotype and methylation statuses are rarely characterized due to previous limitations in genome assembly and detection of modified nucleotides. Also, the extent to which SVs act as methylation quantitative trait loci (SV-mQTLs) is largely unknown. Here, we generated a pangenome graph summarizing SVs in 782 de novo assemblies obtained from Genomic Answers for Kids, capturing 14.6 million CpG dinucleotides that are absent from the CHM13v2 reference (SV-CpGs), thus expanding their number by 43.6%. Using 435 methylomes, we genotyped 4.06 million SV-CpGs, of which 3.93 million (96.8%) are methylated at least once. Nonrepeat sequences contribute 1.59 × 106 novel SV-CpGs, followed by centromeric satellites (6.57 × 105), simple repeats (5.40 × 105), Alu elements (5.07 × 105), satellites (2.17 × 105), LINE-1s (1.83 × 105), and SVA (SINE-VNTR-Alu) elements (1.50 × 105). Centromeric satellites, simple repeats, and SVAs are overrepresented in SV-CpGs versus reference CpGs. Similarly, methylation levels in SV-CpGs are more variable than in reference CpGs. To explore if SVs are potentially causal for functional variation, we measured SV-mQTLs. This revealed over 230,464 methylation bins where the methylation is associated with common SVs within 100 kbp. Finally, we identified 65,659 methylation bins (28.5%) where the leading QTL variant is an SV. In conclusion, we demonstrate that graph pangenomes provide full SV structures, the associated methylation variation, and reveal tens of thousands of SV-mQTLs, underscoring the importance of assembly based analyses of human traits.

Bi-allelic variants in BRF2 are associated with perinatal death and craniofacial anomalies

BACKGROUND

Variants in genes encoding multiple subunits of the RNA Polymerase III complex which synthesizes rRNAs, tRNAs, and other small RNAs were previously associated with neurological disorders, such as syndromic hypomyelination leukodystrophies, pontocerebellar hypoplasia, and cerebellofaciodental syndrome. One new such candidate is BRF2, which encodes a TFIIB-like factor that recruits the RNA polymerase III complex to type 3 promoters to initiate transcription of U6, RnaseP, and 7SK RNAs.

METHODS

We combined sequencing with functional analyses to investigate the effects of BRF2 variants.

RESULTS

We observe that a previously reported significant underrepresentation of double transmission of a splice variant results in recessive lethality in three large Icelandic families with multiple perinatal losses. Using data aggregation, we identified an additional seven individuals worldwide from three unrelated families carrying biallelic variants in BRF2. Affected individuals present a variable phenotype ranging from severe craniofacial anomalies with early death to intellectual disability with motor and speech development. In silico 3D modelling and functional analyses showed functional impairment of the identified variants, e.g., differences in target loci occupancy. Zebrafish knocked down for the orthologous brf2 presented with abnormal escape response, reduced swimming velocity and head size, and craniofacial malformations. These defects were complemented by the human wild-type but not mutated BRF2 mRNA further demonstrating their deleteriousness.

CONCLUSIONS

Overall, our results support the association of biallelic BRF2 variants with a novel neurodevelopmental disease and provide an additional link between RNA polymerase III, its targets and craniofacial anomalies.

Integration of transcriptomics and long-read genomics prioritizes structural variants in rare disease

Rare structural variants (SVs)-insertions, deletions, and complex rearrangements-can cause Mendelian disease, yet they remain difficult to accurately detect and interpret. We sequenced and analyzed Oxford Nanopore Technologies long-read genomes of 68 individuals from the undiagnosed disease network (UDN) with no previously identified diagnostic mutations from short-read sequencing. Using our optimized SV detection pipelines and 571 control long-read genomes, we detected 716 long-read rare (MAF < 0.01) SV alleles per genome on average, achieving a 2.4× increase from short reads. To characterize the functional effects of rare SVs, we assessed their relationship with gene expression from blood or fibroblasts from the same individuals and found that rare SVs overlapping enhancers were enriched (LOR = 0.46) near expression outliers. We also evaluated tandem repeat expansions (TREs) and found 14 rare TREs per genome; notably, these TREs were also enriched near overexpression outliers. To prioritize candidate functional SVs, we developed Watershed-SV, a probabilistic model that integrates expression data with SV-specific genomic annotations, which significantly outperforms baseline models that do not incorporate expression data. Watershed-SV identified a median of eight high-confidence functional SVs per UDN genome. Notably, this included compound heterozygous deletions in FAM177A1 shared by two siblings, which were likely causal for a rare neurodevelopmental disorder. Our observations demonstrate the promise of integrating long-read sequencing with gene expression toward improving the prioritization of functional SVs and TREs in rare disease patients.

The additional diagnostic yield of long-read sequencing in undiagnosed rare diseases

Long-read sequencing (LRS) is a promising technology positioned to study the significant proportion of rare diseases (RDs) that remain undiagnosed as it addresses many of the limitations of short-read sequencing, detecting and clarifying additional disease-associated variants that may be missed by the current standard diagnostic workflow for RDs. Some key areas where additional diagnostic yields may be realized include: (1) detection and resolution of structural variants (SVs); (2) detection and characterization of tandem repeat expansions; (3) coverage of regions of high sequence similarity; (4) variant phasing; (5) the use of de novo genome assemblies for reference-based or graph genome variant detection; and (6) epigenetic and transcriptomic evaluations. Examples from over 50 studies support that the main areas of added diagnostic yield currently lie in SV detection and characterization, repeat expansion assessment, and phasing (with or without DNA methylation information). Several emerging studies applying LRS in cohorts of undiagnosed RDs also demonstrate that LRS can boost diagnostic yields following negative standard-of-care clinical testing and provide an added yield of 7%-17% following negative short-read genome sequencing. With this evidence of improved diagnostic yield, we discuss the incorporation of LRS into the diagnostic care pathway for undiagnosed RDs, including current challenges and considerations, with the ultimate goal of ending the diagnostic odyssey for countless individuals with RDs.

Severus detects somatic structural variation and complex rearrangements in cancer genomes using long-read sequencing

For the detection of somatic structural variation (SV) in cancer genomes, long-read sequencing is advantageous over short-read sequencing with respect to mappability and variant phasing. However, most current long-read SV detection methods are not developed for the analysis of tumor genomes characterized by complex rearrangements and heterogeneity. Here, we present Severus, a breakpoint graph-based algorithm for somatic SV calling from long-read cancer sequencing. Severus works with matching normal samples, supports unbalanced cancer karyotypes, can characterize complex multibreak SV patterns and produces haplotype-specific calls. On a comprehensive multitechnology cell line panel, Severus consistently outperforms other long-read and short-read methods in terms of SV detection F1 score (harmonic mean of the precision and recall). We also illustrate that compared to long-read methods, short-read sequencing systematically misses certain classes of somatic SVs, such as insertions or clustered rearrangements. We apply Severus to several clinical cases of pediatric leukemia/lymphoma, revealing clinically relevant cryptic rearrangements missed by standard genomic panels.

Identification of a novel non-coding deletion in Allan-Herndon-Dudley syndrome by long-read HiFi genome sequencing

BACKGROUND

Allan-Herndon-Dudley syndrome (AHDS) is an X-linked disorder caused by pathogenic variants in the SLC16A2 gene. Although most reported variants are found in protein-coding regions or adjacent junctions, structural variations (SVs) within non-coding regions have not been previously reported.

METHODS

We investigated two male siblings with severe neurodevelopmental disorders and spasticity, who had remained undiagnosed for over a decade and were negative from exome sequencing, utilizing long-read HiFi genome sequencing. We conducted a comprehensive analysis including short-tandem repeats (STRs) and SVs to identify the genetic cause in this familial case.

RESULTS

While coding variant and STR analyses yielded negative results, SV analysis revealed a novel hemizygous deletion in intron 1 of the SLC16A2 gene (chrX:74,460,691 - 74,463,566; 2,876 bp), inherited from their carrier mother and shared by the siblings. Determination of the breakpoints indicates that the deletion probably resulted from Alu/Alu-mediated rearrangements between homologous AluY pairs. The deleted region is predicted to include multiple transcription factor binding sites, such as Stat2, Zic1, Zic2, and FOXD3, which are crucial for the neurodevelopmental process, as well as a regulatory element including an eQTL (rs1263181) that is implicated in the tissue-specific regulation of SLC16A2 expression, notably in skeletal muscle and thyroid tissues.

CONCLUSIONS

This report, to our knowledge, is the first to describe a non-coding deletion associated with AHDS, demonstrating the potential utility of long-read sequencing for undiagnosed patients. Although interpreting variants in non-coding regions remains challenging, our study highlights this region as a high priority for future investigation and functional studies.

Navigating Genetic Testing in Nephrology: Options and Decision-Making Strategies

Technological advances such as next-generation sequencing (NGS) have enabled high-throughput assessment of the human genome, supporting the usage of genetic testing as a first-line tool across clinical medicine. Although individually rare, genetic causes account for end-stage renal disease in 10% to 15% of adults and 70% of children, and in many of these individuals, genetic testing can identify a specific etiology and meaningfully impact management. However, with numerous options for genetic testing available, nephrologists may feel uncomfortable integrating genetics into their clinical practice. Here, we aim to demystify the process of genetic test selection and highlight the opportunities for interdisciplinary collaboration between nephrologists and genetics professionals, thereby supporting precision medicine for patients with kidney disease. We first detail the various clinical genetic testing modalities, highlighting their technical advantages and limitations, and then discuss indications for their usage. Next, we provide a generalized workflow for genetic test selection among individuals with kidney disease and illustrate how this workflow can be applied to genetic test selection across diverse clinical contexts. We then discuss key areas related to the usage of genetic testing in clinical nephrology that merit further research and approaches to investigate them.

Advancing long-read nanopore genome assembly and accurate variant calling for rare disease detection

More than 50% of families with suspected rare monogenic diseases remain unsolved after whole-genome analysis by short-read sequencing (SRS). Long-read sequencing (LRS) could help bridge this diagnostic gap by capturing variants inaccessible to SRS, facilitating long-range mapping and phasing and providing haplotype-resolved methylation profiling. To evaluate LRS's additional diagnostic yield, we sequenced a rare-disease cohort of 98 samples from 41 families, using nanopore sequencing, achieving per sample ∼36× average coverage and 32-kb read N50 from a single flow cell. Our Napu pipeline generated assemblies, phased variants, and methylation calls. LRS covered, on average, coding exons in ∼280 genes and ∼5 known Mendelian disease-associated genes that were not covered by SRS. In comparison to SRS, LRS detected additional rare, functionally annotated variants, including structural variants (SVs) and tandem repeats, and completely phased 87% of protein-coding genes. LRS detected additional de novo variants and could be used to distinguish postzygotic mosaic variants from prezygotic de novos. Diagnostic variants were established by LRS in 11 probands, with diverse underlying genetic causes including de novo and compound heterozygous variants, large-scale SVs, and epigenetic modifications. Our study demonstrates LRS's potential to enhance diagnostic yield for rare monogenic diseases, implying utility in future clinical genomics workflows.

Genomic reanalysis of a pan-European rare-disease resource yields new diagnoses

Genetic diagnosis of rare diseases requires accurate identification and interpretation of genomic variants. Clinical and molecular scientists from 37 expert centers across Europe created the Solve-Rare Diseases Consortium (Solve-RD) resource, encompassing clinical, pedigree and genomic rare-disease data (94.5% exomes, 5.5% genomes), and performed systematic reanalysis for 6,447 individuals (3,592 male, 2,855 female) with previously undiagnosed rare diseases from 6,004 families. We established a collaborative, two-level expert review infrastructure that allowed a genetic diagnosis in 506 (8.4%) families. Of 552 disease-causing variants identified, 464 (84.1%) were single-nucleotide variants or short insertions/deletions. These variants were either located in recently published novel disease genes (n = 67), recently reclassified in ClinVar (n = 187) or reclassified by consensus expert decision within Solve-RD (n = 210). Bespoke bioinformatics analyses identified the remaining 15.9% of causative variants (n = 88). Ad hoc expert review, parallel to the systematic reanalysis, diagnosed 249 (4.1%) additional families for an overall diagnostic yield of 12.6%. The infrastructure and collaborative networks set up by Solve-RD can serve as a blueprint for future further scalable international efforts. The resource is open to the global rare-disease community, allowing phenotype, variant and gene queries, as well as genome-wide discoveries.

Rapid and scalable personalized ASO screening in patient-derived organoids

Personalized antisense oligonucleotides (ASOs) have achieved positive results in the treatment of rare genetic disease1. As clinical sequencing technologies continue to advance, the ability to identify patients with rare disease harbouring pathogenic genetic variants amenable to this therapeutic strategy will probably improve. Here we describe a scalable platform for generating patient-derived cellular models and demonstrate that these personalized models can be used for preclinical evaluation of patient-specific ASOs. We describe protocols for delivery of ASOs to patient-derived organoid models and confirm reversal of disease-associated phenotypes in cardiac organoids derived from a patient with Duchenne muscular dystrophy (DMD) with a structural deletion in the gene encoding dystrophin (DMD) that is amenable to treatment with existing ASO therapeutics. Furthermore, we designed novel patient-specific ASOs for two additional patients with DMD (siblings) with a deep intronic variant in the DMD gene that gives rise to a novel splice acceptor site, incorporation of a cryptic exon and premature transcript termination. We showed that treatment of patient-derived cardiac organoids with patient-specific ASOs results in restoration of DMD expression and reversal of disease-associated phenotypes. The approach outlined here provides the foundation for an expedited path towards the design and preclinical evaluation of personalized ASO therapeutics for a broad range of rare diseases.

Synchronized long-read genome, methylome, epigenome and transcriptome profiling resolve a Mendelian condition

Resolving the molecular basis of a Mendelian condition 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. This 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 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 Mendelian conditions. Notably, the function of each gene was disrupted via a distinct mechanism that required integration of the four 'omes' to resolve. These included 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.

A corpus of GA4GH phenopackets: Case-level phenotyping for genomic diagnostics and discovery

The Global Alliance for Genomics and Health (GA4GH) Phenopacket Schema was released in 2022 and approved by ISO as a standard for sharing clinical and genomic information about an individual, including phenotypic descriptions, numerical measurements, genetic information, diagnoses, and treatments. A phenopacket can be used as an input file for software that supports phenotype-driven genomic diagnostics and for algorithms that facilitate patient classification and stratification for identifying new diseases and treatments. There has been a great need for a collection of phenopackets to test software pipelines and algorithms. Here, we present Phenopacket Store. Phenopacket Store v.0.1.19 includes 6,668 phenopackets representing 475 Mendelian and chromosomal diseases associated with 423 genes and 3,834 unique pathogenic alleles curated from 959 different publications. This represents the first large-scale collection of case-level, standardized phenotypic information derived from case reports in the literature with detailed descriptions of the clinical data and will be useful for many purposes, including the development and testing of software for prioritizing genes and diseases in diagnostic genomics, machine learning analysis of clinical phenotype data, patient stratification, and genotype-phenotype correlations. This corpus also provides best-practice examples for curating literature-derived data using the GA4GH Phenopacket Schema.

Hospital-wide access to genomic data advanced pediatric rare disease research and clinical outcomes

Boston Children's Hospital has established a genomic sequencing and analysis research initiative to improve clinical care for pediatric rare disease patients. Through the Children's Rare Disease Collaborative (CRDC), the hospital offers CLIA-grade exome and genome sequencing, along with other sequencing types, to patients enrolled in specialized rare disease research studies. The data, consented for broad research use, are harmonized and analyzed with CRDC-supported variant interpretation tools. Since its launch, 66 investigators representing 26 divisions and 45 phenotype-based cohorts have joined the CRDC. These studies enrolled 4653 families, with 35% of analyzed cases having a finding either confirmed or under further investigation. This accessible and harmonized genomics platform also supports additional institutional data collections, research and clinical, and now encompasses 13,800+ patients and their families. This has fostered new research projects and collaborations, increased genetic diagnoses and accelerated innovative research via integration of genomics research with clinical care.

Prevalence of Intersex/Differences in Sex Development and Primary Gonadal Insufficiency in a Pediatric Transgender Population

Purpose

This study aims to assess the prevalence of intersex variations/differences in sex development (I/DSDs), associated adrenal conditions, and primary gonadal insufficiency in children with gender dysphoria.

Methods

We performed a comprehensive review of the medical records for individuals who carried the diagnostic codes for gender dysphoria in addition to intersex and/or other conditions associated with sex steroid variations among patients evaluated by pediatric endocrinologists from 2013 to 2022.

Results

We found that 9 of 612 (1.5%) transmasculine (TM) and 4 of 215 (1.9%) transfeminine patients had detectable I/DSDs. Although most patients were diagnosed with I/DSDs before evaluation of gender dysphoria, 4 of 13 (30.7%) were diagnosed with I/DSDs after being referred to endocrinology for gender dysphoria. In all cases, diagnoses were made by the endocrinologists evaluating for gender dysphoria. An additional 0.7% of TM patients were diagnosed with distinct hyperandrogenic adrenal conditions, and 1% of TM patients were diagnosed with primary ovarian insufficiency.

Conclusion

The low, but clinically relevant, prevalence of I/DSDs, distinct adrenal conditions, and primary gonadal insufficiency in this transgender population supports the need for access to individualized expert medical care. Specifically, multidisciplinary clinics with experience in endocrinology may provide specialized support for the transgender community.

Long-read genome sequencing and variant reanalysis increase diagnostic yield in neurodevelopmental disorders

Variant detection from long-read genome sequencing (lrGS) has proven to be more accurate and comprehensive than variant detection from short-read genome sequencing (srGS). However, the rate at which lrGS can increase molecular diagnostic yield for rare disease is not yet precisely characterized. We performed lrGS using Pacific Biosciences "HiFi" technology on 96 short-read-negative probands with rare diseases that were suspected to be genetic. We generated hg38-aligned variants and de novo phased genome assemblies, and subsequently annotated, filtered, and curated variants using clinical standards. New disease-relevant or potentially relevant genetic findings were identified in 16/96 (16.7%) probands, nine of which (8/96, ∼9.4%) harbored pathogenic or likely pathogenic variants. Nine probands (∼9.4%) had variants that were accurately called in both srGS and lrGS and represent changes to clinical interpretation, mostly from recently published gene-disease associations. Seven cases included variants that were only correctly interpreted in lrGS, including copy-number variants (CNVs), an inversion, a mobile element insertion, two low-complexity repeat expansions, and a 1 bp deletion. While evidence for each of these variants is, in retrospect, visible in srGS, they were either not called within srGS data, were represented by calls with incorrect sizes or structures, or failed quality control and filtration. Thus, while reanalysis of older srGS data clearly increases diagnostic yield, we find that lrGS allows for substantial additional yield (7/96, 7.3%) beyond srGS. We anticipate that as lrGS analysis improves, and as lrGS data sets grow allowing for better variant-frequency annotation, the additional lrGS-only rare disease yield will grow over time.

High-coverage nanopore sequencing of samples from the 1000 Genomes Project to build a comprehensive catalog of human genetic variation

Fewer than half of individuals with a suspected Mendelian or monogenic condition receive a precise molecular diagnosis after comprehensive clinical genetic testing. Improvements in data quality and costs have heightened interest in using long-read sequencing (LRS) to streamline clinical genomic testing, but the absence of control data sets for variant filtering and prioritization has made tertiary analysis of LRS data challenging. To address this, the 1000 Genomes Project (1KGP) Oxford Nanopore Technologies Sequencing Consortium aims to generate LRS data from at least 800 of the 1KGP samples. Our goal is to use LRS to identify a broader spectrum of variation so we may improve our understanding of normal patterns of human variation. Here, we present data from analysis of the first 100 samples, representing all 5 superpopulations and 19 subpopulations. These samples, sequenced to an average depth of coverage of 37× and sequence read N50 of 54 kbp, have high concordance with previous studies for identifying single nucleotide and indel variants outside of homopolymer regions. Using multiple structural variant (SV) callers, we identify an average of 24,543 high-confidence SVs per genome, including shared and private SVs likely to disrupt gene function as well as pathogenic expansions within disease-associated repeats that were not detected using short reads. Evaluation of methylation signatures revealed expected patterns at known imprinted loci, samples with skewed X-inactivation patterns, and novel differentially methylated regions. All raw sequencing data, processed data, and summary statistics are publicly available, providing a valuable resource for the clinical genetics community to discover pathogenic SVs.

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.

Leaving no patient behind! Expert recommendation in the use of innovative technologies for diagnosing rare diseases

Genetic diagnosis plays a crucial role in rare diseases, particularly with the increasing availability of emerging and accessible treatments. The International Rare Diseases Research Consortium (IRDiRC) has set its primary goal as: "Ensuring that all patients who present with a suspected rare disease receive a diagnosis within one year if their disorder is documented in the medical literature". Despite significant advances in genomic sequencing technologies, more than half of the patients with suspected Mendelian disorders remain undiagnosed. In response, IRDiRC proposes the establishment of "a globally coordinated diagnostic and research pipeline". To help facilitate this, IRDiRC formed the Task Force on Integrating New Technologies for Rare Disease Diagnosis. This multi-stakeholder Task Force aims to provide an overview of the current state of innovative diagnostic technologies for clinicians and researchers, focusing on the patient's diagnostic journey. Herein, we provide an overview of a broad spectrum of emerging diagnostic technologies involving genomics, epigenomics and multi-omics, functional testing and model systems, data sharing, bioinformatics, and Artificial Intelligence (AI), highlighting their advantages, limitations, and the current state of clinical adaption. We provide expert recommendations outlining the stepwise application of these innovative technologies in the diagnostic pathways while considering global differences in accessibility. The importance of FAIR (Findability, Accessibility, Interoperability, and Reusability) and CARE (Collective benefit, Authority to control, Responsibility, and Ethics) data management is emphasized, along with the need for enhanced and continuing education in medical genomics. We provide a perspective on future technological developments in genome diagnostics and their integration into clinical practice. Lastly, we summarize the challenges related to genomic diversity and accessibility, highlighting the significance of innovative diagnostic technologies, global collaboration, and equitable access to diagnosis and treatment for people living with rare disease.

Complex trait associations in rare diseases and impacts on Mendelian variant interpretation

Emerging evidence implicates common genetic variation - aggregated into polygenic scores (PGS) - in the onset and phenotypic presentation of rare diseases. Here, we comprehensively map individual polygenic liability for 1102 open-source PGS in a cohort of 3059 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. Using this resource, we demonstrate increased polygenic liability in probands with an inherited candidate disease variant (VUS) compared to unaffected carrier parents. Further, we show an enrichment for large-effect rare variants in putative core PGS genes for associated complex traits. Overall, our study supports and expands on previous findings of complex trait associations in rare diseases, implicates polygenic liability as a potential mechanism underlying variable penetrance of candidate causal variants, and provides a framework for identifying novel candidate rare disease genes.

Advancing long-read nanopore genome assembly and accurate variant calling for rare disease detection

More than 50% of families with suspected rare monogenic diseases remain unsolved after whole genome analysis by short read sequencing (SRS). Long-read sequencing (LRS) could help bridge this diagnostic gap by capturing variants inaccessible to SRS, facilitating long-range mapping and phasing, and providing haplotype-resolved methylation profiling. To evaluate LRS's additional diagnostic yield, we sequenced a rare disease cohort of 98 samples, including 41 probands and some family members, using nanopore sequencing, achieving per sample ∼36x average coverage and 32 kilobase (kb) read N50 from a single flow cell. Our Napu pipeline generated assemblies, phased variants, and methylation calls. LRS covered, on average, coding exons in ∼280 genes and ∼5 known Mendelian disease genes that were not covered by SRS. In comparison to SRS, LRS detected additional rare, functionally annotated variants, including SVs and tandem repeats, and completely phased 87% of protein-coding genes. LRS detected additional de novo variants, and could be used to distinguish postzygotic mosaic variants from prezygotic de novos . Eleven probands were solved, with diverse underlying genetic causes including de novo and compound heterozygous variants, large-scale SVs, and epigenetic modifications. Our study demonstrates LRS's potential to enhance diagnostic yield for rare monogenic diseases, implying utility in future clinical genomics workflows.

Beyond the Human Genome Project: The Age of Complete Human Genome Sequences and Pangenome References

The Human Genome Project was an enormous accomplishment, providing a foundation for countless explorations into the genetics and genomics of the human species. Yet for many years, the human genome reference sequence remained incomplete and lacked representation of human genetic diversity. Recently, two major advances have emerged to address these shortcomings: complete gap-free human genome sequences, such as the one developed by the Telomere-to-Telomere Consortium, and high-quality pangenomes, such as the one developed by the Human Pangenome Reference Consortium. Facilitated by advances in long-read DNA sequencing and genome assembly algorithms, complete human genome sequences resolve regions that have been historically difficult to sequence, including centromeres, telomeres, and segmental duplications. In parallel, pangenomes capture the extensive genetic diversity across populations worldwide. Together, these advances usher in a new era of genomics research, enhancing the accuracy of genomic analysis, paving the path for precision medicine, and contributing to deeper insights into human biology.

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 simplifies long-read variant calling with DeepVariant.

Prevalence and Characterization of NOTCH2NLC GGC Repeat Expansions in Koreans: From a Hospital Cohort Analysis to a Population-Wide Study

Background and Objectives

GGC repeat expansions in the NOTCH2NLC gene are associated with a broad spectrum of progressive neurologic disorders, notably, neuronal intranuclear inclusion disease (NIID). We aimed to investigate the population-wide prevalence and clinical manifestations of NOTCH2NLC-related disorders in Koreans.

Methods

We conducted a study using 2 different cohorts from the Korean population. Patients with available brain MRI scans from Seoul National University Hospital (SNUH) were thoroughly reviewed, and NIID-suspected patients presenting the zigzag edging signs underwent genetic evaluation for NOTCH2NLC repeats by Cas9-mediated nanopore sequencing. In addition, we analyzed whole-genome sequencing data from 3,887 individuals in the Korea Biobank cohort to estimate the distribution of the repeat counts in Koreans and to identify putative patients with expanded alleles and neurologic phenotypes.

Results

In the SNUH cohort, among 90 adult-onset leukoencephalopathy patients with unknown etiologies, we found 20 patients with zigzag edging signs. Except for 2 diagnosed with fragile X-associated tremor/ataxia syndrome and 2 with unavailable samples, all 16 patients (17.8%) were diagnosed with NIID (repeat range: 87-217). By analyzing the Korea Biobank cohort, we estimated the distribution of repeat counts and threshold (>64) for Koreans, identifying 6 potential patients with NIID. Furthermore, long-read sequencing enabled the elucidation of transmission and epigenetic patterns of NOTCH2NLC repeats within a family affected by pediatric-onset NIID.

Discussion

This study presents the population-wide distribution of NOTCH2NLC repeats and the estimated prevalence of NIID in Koreans, providing valuable insights into the association between repeat counts and disease manifestations in diverse neurologic disorders.

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 corpus of GA4GH Phenopackets: case-level phenotyping for genomic diagnostics and discovery

The Global Alliance for Genomics and Health (GA4GH) Phenopacket Schema was released in 2022 and approved by ISO as a standard for sharing clinical and genomic information about an individual, including phenotypic descriptions, numerical measurements, genetic information, diagnoses, and treatments. A phenopacket can be used as an input file for software that supports phenotype-driven genomic diagnostics and for algorithms that facilitate patient classification and stratification for identifying new diseases and treatments. There has been a great need for a collection of phenopackets to test software pipelines and algorithms. Here, we present phenopacket-store. Version 0.1.12 of phenopacket-store includes 4916 phenopackets representing 277 Mendelian and chromosomal diseases associated with 236 genes, and 2872 unique pathogenic alleles curated from 605 different publications. This represents the first large-scale collection of case-level, standardized phenotypic information derived from case reports in the literature with detailed descriptions of the clinical data and will be useful for many purposes, including the development and testing of software for prioritizing genes and diseases in diagnostic genomics, machine learning analysis of clinical phenotype data, patient stratification, and genotype-phenotype correlations. This corpus also provides best-practice examples for curating literature-derived data using the GA4GH Phenopacket Schema.

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 22, 019 deletions, 2,041 duplications, 87,826 insertions, and 214 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 = 1×10-03). This difference was not observed in the lowest-ranked gene set (P = 0.15). 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.

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.

Applications of genome sequencing as a single platform for clinical constitutional genetic testing

The number of human disease genes has dramatically increased over the past decade, largely fueled by ongoing advances in sequencing technologies. In parallel, the number of available clinical genetic tests has also increased, including the utilization of exome sequencing for undiagnosed diseases. Although most clinical sequencing tests have been centered on enrichment-based multigene panels and exome sequencing, the continued improvements in performance and throughput of genome sequencing suggest that this technology is emerging as a potential platform for routine clinical genetic testing. A notable advantage is a single workflow with the opportunity to reflexively interrogate content as clinically indicated; however, challenges with implementing routine clinical genome sequencing still remain. This review is centered on evaluating the applications of genome sequencing as a single platform for clinical constitutional genetic testing, including its potential utility for diagnostic testing, carrier screening, cytogenomic molecular karyotyping, prenatal testing, mitochondrial genome interrogation, and pharmacogenomic and polygenic risk score testing.

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.

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.