Diagnostic implications of pitfalls in causal variant identification based on 4577 molecularly characterized families

The evolving genetic landscape of neuromuscular fetal akinesias

Fetal akinesia is a broad term used to describe absent (or reduced, fetal hypokinesia) fetal movements, and it can be detected as early as the first trimester. Depending on the developmental age of onset, anything that interferes or limits the normal in utero movement results in a range of deformations affecting multiple organs and organ systems. Arthrogryposis, also termed arthrogryposis multiplex congenita (AMC), is a definitive terminology for multiple congenital contractures, with two major subgroups; amyoplasia and distal arthrogryposis (DA). The spectrum includes fetal akinesia deformation sequence (FADS), lethal congenital contracture syndrome (LCCS), and multiple pterygium syndrome (MPS). Variants in more than >400 genes are known to cause AMC, and it is increasingly recognized that variants in genes encoding critical components (including ventral horn cell, peripheral nerve, neuromuscular junction, skeletal muscle) of the extended motor unit underlie ∼40% of presentations. With unbiased screening approaches, including sequencing of comprehensive disease gene panels, exomes and genomes, novel genes and phenotypic expansions associated with known human disease genes have been uncovered in the setting of fetal akinesia. Autosomal-recessive titinopathy is the most frequent genetic cause of AMC. Accurate genetic diagnosis is critical to genetic counseling and informing family planning. Around 50% remain undiagnosed following comprehensive prenatal, diagnostic or research screening. Comprehensive phenotyping and periodic reanalysis with appropriate genomic tools are valuable strategies when faced with initial inconclusive results. There are likely many novel causative genes still to identify, which will inform our understanding of the molecular pathways underlying early human development and in utero movement.

Near-complete Middle Eastern genomes refine autozygosity and enhance disease-causing and population-specific variant discovery

Advances in long-read sequencing have enabled routine complete assembly of human genomes, but much remains to be done to represent broader populations and show impact on disease-gene discovery. Here, we report highly accurate, near-complete and phased genomes from six Middle Eastern (ME) family trios (n = 18) with neurodevelopmental conditions, representing ancestries from Sudan, Jordan, Syria, Qatar and Afghanistan. These genomes revealed 42.2 Mb of new sequence (13.8% impacting known genes), 75 new HLA/KIR alleles and strong signals of inbreeding, with ROH covering up to one-third of chromosomes 6 and 12 in one individual. Using assembly-based variant calling, we identified 23 de novo and recessive variants as strong candidates for causing previously unresolved symptoms in the probands. The ME genomes revealed unique variation relative to existing references, showing enhanced mappability and variant calling. These results underscore the value of de novo assembly for disease variant discovery and the need for sampled ME-specific references to better characterize population-relevant variation.

Severe Joubert syndrome in family with homozygous POC1B p.Arg106Pro variant is due to a co-inherited deep-intronic mutation in the neighboring CEP290 gene

"En bloc" inheritance of point mutations in adjacent genes has rarely been described. We have previously reported a family with severe, mostly early-lethal Joubert syndrome (JBTS) with early-onset severe retinal dystrophy (EOSRD) and polycystic kidney disease (PKD), which at that time had been attributed to a homozygous pathogenic missense variant, p.Arg106Pro (c.317G>C), in the ciliary POC1B gene. Because this and other POC1B variants were, in subsequent studies, only reported in patients with non-syndromic childhood or early-adult-onset macular dystrophy, we have now reassessed our index patient by long-read high-fidelity (HiFi) whole-genome sequencing (LR-WGS). We identified a homozygous deep-intronic variant, c.2818-657T>G, in CEP290, a JBTS/Meckel syndrome-associated gene on chromosome 12q21, only 1.28 Mb from the N terminus of POC1B. cDNA analysis revealed aberrant splicing with the frame-shifting inclusion of 37 bp from CEP290 intron 25, predicting the loss of CEP290 function. EOSRD and PKD can fully be ascribed to this CEP290 variant, whose effect outshines the "background" non-syndromic POC1B retinopathy and co-segregates with the severe syndromic phenotype. Our novel findings in this family no longer justify POC1B as a JBTS gene. This co-inheritance of two ciliopathies, with the clinically decisive variant hidden deep in an intron, exemplifies the importance of WGS for achieving the complete diagnosis in challenging cases.

Novel biallelic COL25A1 variants broaden the clinical spectrum from congenital cranial dysinnervation disorders to fetal lethal phenotypes

Biallelic variants in COL25A1 have been associated with isolated congenital cranial dysinnervation disorders (CCDDs) and arthrogryposis multiplex congenital (AMC) with or without CCDD. COL25A1 encodes collagen XXV that belongs to the subfamily of membrane-associated collagens with interrupted triple helices. COL25A1 contains four non-collagenous and three collagenous domains. Three alternatively spliced COL25A1 transcript variants are known. In mice, Col25a1 is required for intramuscular motor innervation and cranial motor neuron development. We report seven subjects with novel biallelic COL25A1 pathogenic variants, including three AMC-affected individuals, one of whom died in infancy, and four unrelated fetuses. We expand the associated phenotypic spectrum as fetuses showed lethal phenotypes including reduced or no movement, contractures, and hydrops in three and growth retardation and skeletal abnormalities in one. The molecular spectrum includes two microdeletions encompassing several 5' or 3' exons, two missense, one nonsense, one frameshift, and one variant affecting splicing. In fibroblasts of the subject who was compound heterozygous for the c.367G > C and c.1198G > T variants, we identified skipping of exon 3 in COL25A1 mRNAs due to the G-to-C change. These aberrantly spliced transcripts were subject to nonsense-mediated mRNA decay. Analysis of transcriptome sequencing data from primary human fibroblasts without COL25A1 pathogenic variants revealed novel COL25A1 exon-exon junctions and 13 not previously annotated alternatively spliced in-frame exons. We hypothesized that interindividual variation in the splicing of COL25A1 exons in different tissues may underlie the variable phenotypes in the affected individuals.

Arab founder variants: Contributions to clinical genomics and precision medicine

BACKGROUND

Founder variants are ancestral variants shared by individuals who are not closely related. The large effect size of some of these variants in the context of Mendelian disorders offers numerous precision medicine opportunities.

METHODS

Using one of the largest datasets on Mendelian disorders in the Middle East, we identified 2,908 medically relevant founder variants derived from 18,360 exomes and genomes and investigated their contribution to the clinical annotation of the human genome.

FINDINGS

Strikingly, ∼34% of Arab founder variants are absent in gnomAD. We found a strong contribution of Arab founder variants to the identification of novel gene-disease links (n = 224) and the support/dispute (n = 81 support, n = 101 dispute) of previously reported candidate gene-disease links. The powerful segregation evidence generated by Arab founder variants allowed many ClinVar and Human Gene Mutation Database variants to be reclassified. Overall, 39.5% of diagnostic reports from our clinical lab are based on founder variants, and 19.41% of tested individuals carry at least one pathogenic founder variant. The presumptive loss-of-function mechanism that typically underlies autosomal recessive diseases means that Arab founder variants also offer unique opportunities in "druggable genome" research. Arab founder variants were also informative of migration patterns in the Middle East consistent with documented historical accounts.

CONCLUSIONS

We highlight the contribution of founder variants from an under-represented population group to precision medicine and inform future prevention programs. Our study also sheds light on the added value of these variants in supplementing other lines of research in tracing population history.

FUNDING

There is no funding for this work.

Genomics of rare diseases in the Greater Middle East

The Greater Middle East (GME) represents a concentrated region of unparalleled genetic diversity, characterized by an abundance of distinct alleles, founder mutations and extensive autozygosity driven by high consanguinity rates. These genetic hallmarks present a unique, yet vastly untapped resource for genomic research on Mendelian diseases. Despite this immense potential, the GME continues to face substantial challenges in comprehensive data collection and analysis. This Perspective highlights the region's unique position as a natural laboratory for genetic discovery and explores the challenges that have stifled progress thus far. Importantly, we propose strategic solutions, advocating for an all-inclusive research approach. With targeted investment and focused efforts, the latent genetic wealth in the GME can be transformed into a global hub for genomic research that will redefine and advance our understanding of the human genome.

Mapping the genetic landscape of treatable inherited metabolic disorders in a large Middle Eastern biobank

PURPOSE

To date, approximately 1400 inherited metabolic disorders (IMDs) have been described, some of which are treatable. It is estimated that 2% to 3% of live births worldwide are affected by treatable IMDs. Roughly 80% of IMDs are autosomal recessive, leading to a potentially higher incidence in regions with high consanguinity.

METHODOLOGY

The study utilized genome sequencing data from 14,060 adult Qatari participants who were recruited by the Qatar Biobank and sequenced by the Qatar Genome Program. The genome sequencing data were analyzed for 125 nuclear genes known to be associated with 115 treatable IMDs.

RESULTS

Our study identified 253 pathogenic/likely pathogenic single-nucleotide variations associated with 69 treatable IMDs, including 211 known and 42 novel predicted loss-of-function variants. We estimated that approximately 1 in 13 unrelated individuals (8%) carry a heterozygous pathogenic variant for at least 1 of 46 treatable IMDs. Notably, phenylketonuria/hyperphenylalaninemia and homocystinuria had among the highest carrier frequencies (1 in 68 and 1 in 85, respectively).

CONCLUSION

Population-based studies of treatable IMDs, particularly in globally under-studied populations, can identify high-frequency alleles segregating in the community and inform public health policies, including carrier and newborn screening.

NanoRanger enables rapid single-base-pair resolution of genomic disorders

BACKGROUND

Delineating base-resolution breakpoints of complex rearrangements is crucial for an accurate clinical understanding of pathogenic variants and for carrier screening within family networks or the broader population. However, despite advances in genetic testing using short-read sequencing (SRS), this task remains costly and challenging.

METHODS

This study addresses the challenges of resolving missing disease-causing breakpoints in complex genomic disorders with suspected homozygous rearrangements by employing multiple long-read sequencing (LRS) strategies, including a novel and efficient strategy named nanopore-based rapid acquisition of neighboring genomic regions (NanoRanger). NanoRanger does not require large amounts of ultrahigh-molecular-weight DNA and stands out for its ease of use and rapid acquisition of large genomic regions of interest with deep coverage.

FINDINGS

We describe a cohort of 16 familial cases, each harboring homozygous rearrangements that defied breakpoint determination by SRS and optical genome mapping (OGM). NanoRanger identified the breakpoints with single-base-pair resolution, enabling accurate determination of the carrier status of unaffected family members as well as the founder nature of these genomic lesions and their frequency in the local population. The resolved breakpoints revealed that repetitive DNA, gene regulatory elements, and transcription activity contribute to genome instability in these novel recessive rearrangements.

CONCLUSIONS

Our data suggest that NanoRanger greatly improves the success rate of resolving base-resolution breakpoints of complex genomic disorders and expands access to LRS for the benefit of patients with Mendelian disorders.

FUNDING

M.L. is supported by KAUST Baseline Award no. BAS/1/1080-01-01 and KAUST Research Translation Fund Award no. REI/1/4742-01.

A founder variant expands the phenotype of WNT7B-related PDAC syndrome

Pulmonary hypoplasia, Diaphragmatic anomalies, Anophthalmia/microphthalmia, and Cardiac defects (PDAC) syndrome is a genetically heterogeneous multiple congenital malformation syndrome. Although pathogenic variants in RARB and STRA6 are established causes of PDAC, many PDAC cases remain unsolved at the molecular level. Recently, we proposed biallelic WNT7B variants as a novel etiology based on several families with typical features of PDAC syndrome albeit with variable expressivity. Here, we report three patients from two families that share a novel founder variant in WNT7B (c.739C > T; Arg247Trp). The phenotypic expression of this variant ranges from typical PDAC features to isolated genitourinary anomalies. Similar to previously reported PDAC-associated WNT7B variants, this variant was found to significantly impair WNT7B signaling activity further corroborating its proposed pathogenicity. This report adds further evidence to WNT7B-related PDAC and expands its variable expressivity.

Human ABL1 deficiency syndrome (HADS) is a recognizable syndrome distinct from ABL1-related congenital heart defects and skeletal malformations syndrome

Germline gain of function variants in the oncogene ABL1 cause congenital heart defects and skeletal malformations (CHDSKM) syndrome. Whether a corresponding ABL1 deficiency disorder exists in humans remains unknown although developmental defects in mice deficient for Abl1 support this notion. Here, we describe two multiplex consanguineous families, each segregating a different homozygous likely loss of function variant in ABL1. The associated phenotype is multiple congenital malformations and distinctive facial dysmorphism that are opposite in many ways to CHDSKM. We suggest that a tight balance of ABL1 activity is required during embryonic development and that both germline gain of function and loss of function variants result in distinctively different allelic congenital malformation disorders.

Further delineation of the phenotypic and metabolomic profile of ALDH1L2-related neurodevelopmental disorder

ALDH1L2, a mitochondrial enzyme in folate metabolism, converts 10-formyl-THF (10-formyltetrahydrofolate) to THF (tetrahydrofolate) and CO2. At the cellular level, deficiency of this NADP+-dependent reaction results in marked reduction in NADPH/NADP+ ratio and reduced mitochondrial ATP. Thus far, a single patient with biallelic ALDH1L2 variants and the phenotype of a neurodevelopmental disorder has been reported. Here, we describe another patient with a neurodevelopmental disorder associated with a novel homozygous missense variant in ALDH1L2, Pro133His. The variant caused marked reduction in the ALDH1L2 enzyme activity in skin fibroblasts derived from the patient as probed by 10-FDDF, a stable synthetic analog of 10-formyl-THF. Additional associated abnormalities in these fibroblasts include reduced NADPH/NADP+ ratio and pool of mitochondrial ATP, upregulated autophagy and dramatically altered metabolomic profile. Overall, our study further supports a link between ALDH1L2 deficiency and abnormal neurodevelopment in humans.

Beyond the exome: utility of long-read whole genome sequencing in exome-negative autosomal recessive diseases

BACKGROUND

Long-read whole genome sequencing (lrWGS) has the potential to address the technical limitations of exome sequencing in ways not possible by short-read WGS. However, its utility in autosomal recessive Mendelian diseases is largely unknown.

METHODS

In a cohort of 34 families in which the suspected autosomal recessive diseases remained undiagnosed by exome sequencing, lrWGS was performed on the Pacific Bioscience Sequel IIe platform.

RESULTS

Likely causal variants were identified in 13 (38%) of the cohort. These include (1) a homozygous splicing SV in TYMS as a novel candidate gene for lethal neonatal lactic acidosis, (2) a homozygous non-coding SV that we propose impacts STK25 expression and causes a novel neurodevelopmental disorder, (3) a compound heterozygous SV in RP1L1 with complex inheritance pattern in a family with inherited retinal disease, (4) homozygous deep intronic variants in LEMD2 and SNAP91 as novel candidate genes for neurodevelopmental disorders in two families, and (5) a promoter SNV in SLC4A4 causing non-syndromic band keratopathy. Surprisingly, we also encountered causal variants that could have been identified by short-read exome sequencing in 7 families. The latter highlight scenarios that are especially challenging at the interpretation level.

CONCLUSIONS

Our data highlight the continued need to address the interpretation challenges in parallel with efforts to improve the sequencing technology itself. We propose a path forward for the implementation of lrWGS sequencing in the setting of autosomal recessive diseases in a way that maximizes its utility.

Evaluating the association of biallelic OGDHL variants with significant phenotypic heterogeneity

BACKGROUND

Biallelic variants in OGDHL, encoding part of the α-ketoglutarate dehydrogenase complex, have been associated with highly heterogeneous neurological and neurodevelopmental disorders. However, the validity of this association remains to be confirmed. A second OGDHL patient cohort was recruited to carefully assess the gene-disease relationship.

METHODS

Using an unbiased genotype-first approach, we screened large, multiethnic aggregated sequencing datasets worldwide for biallelic OGDHL variants. We used CRISPR/Cas9 to generate zebrafish knockouts of ogdhl, ogdh paralogs, and dhtkd1 to investigate functional relationships and impact during development. Functional complementation with patient variant transcripts was conducted to systematically assess protein functionality as a readout for pathogenicity.

RESULTS

A cohort of 14 individuals from 12 unrelated families exhibited highly variable clinical phenotypes, with the majority of them presenting at least one additional variant, potentially accounting for a blended phenotype and complicating phenotypic understanding. We also uncovered extreme clinical heterogeneity and high allele frequencies, occasionally incompatible with a fully penetrant recessive disorder. Human cDNA of previously described and new variants were tested in an ogdhl zebrafish knockout model, adding functional evidence for variant reclassification. We disclosed evidence of hypomorphic alleles as well as a loss-of-function variant without deleterious effects in zebrafish variant testing also showing discordant familial segregation, challenging the relationship of OGDHL as a conventional Mendelian gene. Going further, we uncovered evidence for a complex compensatory relationship among OGDH, OGDHL, and DHTKD1 isoenzymes that are associated with neurodevelopmental disorders and exhibit complex transcriptional compensation patterns with partial functional redundancy.

CONCLUSIONS

Based on the results of genetic, clinical, and functional studies, we formed three hypotheses in which to frame observations: biallelic OGDHL variants lead to a highly variable monogenic disorder, variants in OGDHL are following a complex pattern of inheritance, or they may not be causative at all. Our study further highlights the continuing challenges of assessing the validity of reported disease-gene associations and effects of variants identified in these genes. This is particularly more complicated in making genetic diagnoses based on identification of variants in genes presenting a highly heterogenous phenotype such as "OGDHL-related disorders".