A pipeline combining multiple strategies for prioritizing heterozygous variants for the identification of candidate genes in exome datasets

Towards a standard benchmark for variant and gene prioritisation algorithms: PhEval - Phenotypic inference Evaluation framework

Background

Computational approaches to support rare disease diagnosis are challenging to build, requiring the integration of complex data types such as ontologies, gene-to-phenotype associations, and cross-species data into variant and gene prioritisation algorithms (VGPAs). However, the performance of VGPAs has been difficult to measure and is impacted by many factors, for example, ontology structure, annotation completeness or changes to the underlying algorithm. Assertions of the capabilities of VGPAs are often not reproducible, in part because there is no standardised, empirical framework and openly available patient data to assess the efficacy of VGPAs - ultimately hindering the development of effective prioritisation tools.

Results

In this paper, we present our benchmarking tool, PhEval, which aims to provide a standardised and empirical framework to evaluate phenotype-driven VGPAs. The inclusion of standardised test corpora and test corpus generation tools in the PhEval suite of tools allows open benchmarking and comparison of methods on standardised data sets.

Conclusions

PhEval and the standardised test corpora solve the issues of patient data availability and experimental tooling configuration when benchmarking and comparing rare disease VGPAs. By providing standardised data on patient cohorts from real-world case-reports and controlling the configuration of evaluated VGPAs, PhEval enables transparent, portable, comparable and reproducible benchmarking of VGPAs. As these tools are often a key component of many rare disease diagnostic pipelines, a thorough and standardised method of assessment is essential for improving patient diagnosis and care.

Variant ranking pipeline for complex familial disorders

Identifying genetic susceptibility factors for complex disorders remains a challenging task. To analyze collections of small and large pedigrees where genetic heterogeneity is likely, but biological commonalities are plausible, we have developed a weights-based pipeline to prioritize variants and genes. The Weights-based vAriant Ranking in Pedigrees (WARP) pipeline prioritizes variants using 5 weights: disease incidence rate, number of cases in a family, genome fraction shared amongst cases in a family, allele frequency and variant deleteriousness. Weights, except for the population allele frequency weight, are normalized between 0 and 1. Weights are combined multiplicatively to produce family-specific-variant weights that are then averaged across all families in which the variant is observed to generate a multifamily weight. Sorting multifamily weights in descending order creates a ranked list of variants and genes for further investigation. WARP was validated using familial melanoma sequence data from the European Genome-phenome Archive. The pipeline identified variation in known germline melanoma genes POT1, MITF and BAP1 in 4 out of 13 families (31%). Analysis of the other 9 families identified several interesting genes, some of which might have a role in melanoma. WARP provides an approach to identify disease predisposing genes in studies with small and large pedigrees.

Evolving fetal phenotypes and clinical impact of progressive prenatal exome sequencing pathways: cohort study

OBJECTIVES

To determine (1) the diagnostic yield and turnaround time (TAT) of two consecutive prenatal exome sequencing (ES) pathways, (2) the evolution of the fetal phenotype and (3) the clinical impact of detecting causative pathogenic variants and incidental findings.

METHODS

This was a retrospective cohort analysis of prospectively collected fetal cases that underwent trio ES in the presence of a structural anomaly and normal chromosomal microarray testing in the West Midlands Regional Genetics Laboratory, Birmingham, UK. The study included two phases: (1) between July 2018 and October 2020, the clinical pathway from the Prenatal Assessment of Genomes and Exomes (PAGE) study was adopted and involved prenatal trio ES based on a panel of 1542 development disorder genes and case selection by a multidisciplinary team; (2) between October 2020 and July 2021, prenatal trio ES investigation was based on the National Health Service (NHS) England R21 pathway, with definitive inclusion criteria and a panel of 1205 prenatally relevant genes. Deep phenotyping was performed throughout pregnancy and postnatally.

RESULTS

A total of 54 cases were included. The diagnostic yield before vs after R21 pathway implementation was 28.0% (7/25) and 55.2% (16/29), respectively (P = 0.04). The respective values for mean TAT were 54.0 days (range, 14-213 days) and 14.2 days (range, 3-29 days). In cases in which a causative pathogenic variant was identified and in which the pregnancy reached the third trimester, additional anomalies were detected between the second and third trimesters in 73.3% (11/15) of cases, predominantly secondary to progressive hydropic features (3/11 (27.3%)), arthrogryposis (3/11 (27.3%)) or brain anomaly (2/11 (18.2%)). In three cases, a variant of uncertain significance was reclassified to likely pathogenic based on postnatal information. Detection of a causative pathogenic variant had a significant clinical impact in 78.3% (18/23) of cases, most frequently affecting decision-making regarding the course of the pregnancy and neonatal management (7/18 (38.9%)).

CONCLUSIONS

Prenatal ES using the NHS England R21 pathway showed great promise when applied to this cohort, allowing a genetic diagnosis to be made in over half of preselected cases with a fetal structural anomaly on ultrasound. Monitoring and real-time updating of fetal phenotype and reclassification of variants based on postnatal findings is vital to increase the clinical impact that is already evident from this emerging genomic technology. © 2021 International Society of Ultrasound in Obstetrics and Gynecology.

Genome sequencing data analysis for rare disease gene discovery

Rare diseases occur in a smaller proportion of the general population, which is variedly defined as less than 200 000 individuals (US) or in less than 1 in 2000 individuals (Europe). Although rare, they collectively make up to approximately 7000 different disorders, with majority having a genetic origin, and affect roughly 300 million people globally. Most of the patients and their families undergo a long and frustrating diagnostic odyssey. However, advances in the field of genomics have started to facilitate the process of diagnosis, though it is hindered by the difficulty in genome data analysis and interpretation. A major impediment in diagnosis is in the understanding of the diverse approaches, tools and datasets available for variant prioritization, the most important step in the analysis of millions of variants to select a few potential variants. Here we present a review of the latest methodological developments and spectrum of tools available for rare disease genetic variant discovery and recommend appropriate data interpretation methods for variant prioritization. We have categorized the resources based on various steps of the variant interpretation workflow, starting from data processing, variant calling, annotation, filtration and finally prioritization, with a special emphasis on the last two steps. The methods discussed here pertain to elucidating the genetic basis of disease in individual patient cases via trio- or family-based analysis of the genome data. We advocate the use of a combination of tools and datasets and to follow multiple iterative approaches to elucidate the potential causative variant.

Integrated Genome and Transcriptome Sequencing to Solve a Neuromuscular Puzzle: Miyoshi Muscular Dystrophy and Early Onset Primary Dystonia in Siblings of the Same Family

BACKGROUND

Neuromuscular disorders (NMD), many of which are hereditary, affect muscular function. Due to advances in high-throughput sequencing technologies, the diagnosis of hereditary NMDs has dramatically improved in recent years.

METHODS AND RESULTS

In this study, we report an family with two siblings exhibiting two different NMD, Miyoshi muscular dystrophy (MMD) and early onset primary dystonia (EOPD). Whole exome sequencing (WES) identified a novel monoallelic frameshift deletion mutation (dysferlin: c.4404delC/p.I1469Sfs∗17) in the Dysferlin gene in the index patient who suffered from MMD. This deletion was inherited from his unaffected father and was carried by his younger sister with EOPD. However, immunostaining staining revealed an absence of dysferlin expression in the proband's muscle tissue and thus suggested the presence of the second underlying mutant allele in dysferlin. Using integrated RNA sequencing (RNA-seq) and whole genome sequencing (WGS) of muscle tissue, a novel deep intronic mutation in dysferlin (dysferlin: c.5341-415A > G) was discovered in the index patient. This mutation caused aberrant mRNA splicing and inclusion of an additional pseudoexon (PE) which we termed PE48.1. This PE was inherited from his unaffected mother. PE48.1 inclusion altered the Dysferlin sequence, causing premature termination of translation.

CONCLUSION

Using integrated genome and transcriptome sequencing, we discovered hereditary MMD and EOPD affecting two siblings of same family. Our results added further weight to the combined use of RNA-seq and WGS as an important method for detection of deep intronic gene mutations, and suggest that integrated sequencing assays are an effective strategy for the diagnosis of hereditary NMDs.

Genetic Inheritance and Its Contribution to Tinnitus

Tinnitus is the abnormal perception of sound that affects more than 15% of adult population around the globe. Severe tinnitus is considered a complex disorder that arises as result of the interaction of genetic and environmental factors, and it is associated with several comorbidities such as hearing loss, anxiety, and insomnia. We begin this review with an introduction to human molecular genetics and the role of genetic variation on the inheritance. There are some genetic reports on tinnitus heritability including concordance studies in twins and adoptees or aggregation in families providing some evidence for familial aggregation in patients with severe tinnitus and high concordance in monozygotic twins with bilateral tinnitus. So, sex differences in familial aggregation and heritability of bilateral tinnitus suggest a potential sexual dimorphism in tinnitus inheritance.Molecular genetic studies have been demonstrated to be a useful tool to understand the role of genetic variation in rare diseases and complex disorders. The reported associations in common variants in neurotrophic factors such as GDNF, BDNF, or potassium channels genes were underpowered, and the lack of replication questions these findings. Although candidate gene approaches have failed in replicating these genetic associations, the development of high throughput sequencing technology and the selection of extreme phenotypes are strategies that will allow the clinicians and researchers to combine genetic information with clinical data to implement a personalized diagnosis and therapy in patients with tinnitus.

TRPM7 as a Candidate Gene for Vestibular Migraine

Objectives: Vestibular migraine (VM) is a common vestibular disorder, and familial aggregation of VM with autosomal-dominant inheritance has been described, which supports a genetic background. This study aimed to describe the clinical phenotype of a family with VM, and identify a candidate gene for VM.

Methods: We recruited six individuals (four affected and two unaffected) from three consecutive generations of a Korean family with VM, and performed whole-exome sequencing to search for candidate genes.

Results: All affected individuals presented with recurrent vertigo, headache, and nausea/vomiting that fulfilled the diagnostic criteria of VM. Two individuals also experienced transient hemiparesis or dysarthria during the episodes. The symptoms were triggered by physical or emotional stress. Interictal examinations showed uni- or bi-directional horizontal gaze-evoked nystagmus in three of the individuals. They had no causative mutations in genes causing familial hemiplegic migraine or episodic ataxia. Through whole-exome sequencing from three affected individuals, we identified a nonsense mutation c.3526C>T in TRPM7 that encodes a cation channel selective to Ca²⁺ and Mg²⁺.

Conclusions: Alterations in intracellular Ca²⁺ and Mg²⁺ homeostasis by TRPM7 mutation may contribute to the development of the VM phenotype. Our result suggest that TRPM7 is a novel candidate gene for VM.

Episodic Vestibular Syndrome with Hyperventilation-Induced Downbeat Nystagmus

Hyperventilation-induced downbeat nystagmus (HV-DBN) has been reported in cerebellar disorders and explained by a loss of the inhibitory cerebellar output via a metabolic effect on cerebellar Ca²⁺ channels. The aim of this study was to determine the clinical characteristics and underlying pathogenesis of episodic vestibular syndrome (EVS) with HV-DBN. Of 667 patients with EVS, we recruited 22 with HV-DBN and assessed their clinical characteristics, video-oculographic findings, and the results of molecular genetic analyses. The age at symptom onset was 47.5 ± 13.0 years (mean ± SD), and there was a female preponderance (n = 15, 68%). The duration of vertigo/dizziness attacks ranged from minutes to a few days, and 11 patients (50%) fulfilled the diagnostic criteria for vestibular migraine. HV-induced new-onset DBN in 8 patients, while the remaining 14 showed augmentation of spontaneous DBN by HV. The maximum slow-phase velocity of HV-DBN ranged from 2.2 to 11.9°/s, which showed a statistical difference with that of spontaneous DBN (median = 4.95, IQR = 3.68-6.55 vs. median = 1.25, IQR = 0.20-2.15, p < 0.001). HV-DBN was either purely downbeat (n = 11) or accompanied with small horizontal components (n = 11). Other neuro-otologic findings included perverted head-shaking nystagmus (n = 11), central positional nystagmus (n = 7), saccadic pursuit (n = 3), and horizontal gaze-evoked nystagmus (n = 1). Gene expression profiling with a bioinformatics analysis identified 43 upregulated and 49 downregulated differentially expressed genes (DEGs) in patients with EVS and HV-DBN and revealed that the downregulated DEGs were significantly enriched in terms related to the ribosome pathway. Our results suggest that the underlying cerebellar dysfunction would be responsible for paroxysmal attacks of vertigo in patients with EVS and HV-DBN.

Specific phenotype semantics facilitate gene prioritization in clinical exome sequencing

Selection and prioritization of phenotype-centric variants remains a challenging part of variant analysis and interpretation in clinical exome sequencing. Phenotype-driven shortlisting of patient-specific gene lists can avoid missed diagnosis. Here, we analyzed the relevance of using primary Human Phenotype Ontology identifiers (HPO IDs) in prioritizing Mendelian disease genes across 30 in-house, 10 previously reported, and 10 recently published cases using three popular web-based gene prioritization tools (OMIMExplorer, VarElect & Phenolyzer). We assessed partial HPO-based gene prioritization using randomly chosen and top 10%, 30%, and 50% HPO IDs based on information content and found high variance within rank ratios across the former vs the latter. This signified that randomly selected less-specific HPO IDs for a given disease phenotype performed poorly by ranking probe gene farther away from the top rank. In contrast, the use of top 10%, 30%, and 50% HPO IDs individually could rank the probe gene among the top 1% in the ranked list of genes that was equivalent to the results when the full list of HPO IDs were used. Hence, we conclude that use of just the top 10% of HPO IDs chosen based on information content is sufficient for ranking the probe gene at top position. Our findings provide practical guidance for utilizing structured phenotype semantics and web-based gene-ranking tools to aid in identifying known as well unknown candidate gene associations in Mendelian disorders.

Identification of Specific Nuclear Genetic Loci and Genes That Interact With the Mitochondrial Genome and Contribute to Fecundity in Caenorhabditis elegans

Previous studies have found that fecundity is a multigenic trait regulated, in part, by mitochondrial-nuclear (mit-n) genetic interactions. However, the identification of specific nuclear genetic loci or genes interacting with the mitochondrial genome and contributing to the quantitative trait fecundity is an unsolved issue. Here, a panel of recombinant inbred advanced intercrossed lines (RIAILs), established from a cross between the N2 and CB4856 strains of C. elegans, were used to characterize the underlying genetic basis of mit-n genetic interactions related to fecundity. Sixty-seven single nucleotide polymorphisms (SNPs) were identified by association mapping to be linked with fecundity among 115 SNPs linked to mitotype. This indicated significant epistatic effects between nuclear and mitochondria genetics on fecundity. In addition, two specific nuclear genetic loci interacting with the mitochondrial genome and contributing to fecundity were identified. A significant reduction in fecundity was observed in the RIAILs that carried CB4856 mitochondria and a N2 genotype at locus 1 or a CB4856 genotype at locus 2 relative to the wild-type strains. Then, a hybrid strain (CNC10) was established, which was bred as homoplasmic for the CB4856 mtDNA genome and N2 genotype at locus 1 in the CB4856 nuclear background. The mean fecundity of CNC10 was half the fecundity of the control strain. Several functional characteristics of the mitochondria in CNC10 were also influenced by mit-n interactions. Overall, experimental evidence was presented that specific nuclear genetic loci or genes have interactions with the mitochondrial genome and are associated with fecundity. In total, 18 genes were identified using integrative approaches to have interactions with the mitochondrial genome and to contribute to fecundity.

Clinical and Functional Characterization of a Missense ELF2 Variant in a CANVAS Family

Cerebellar ataxia with neuropathy and bilateral vestibular areflexia syndrome (CANVAS) is a rare disorder with an unknown etiology. We present a British family with presumed autosomal dominant CANVAS with incomplete penetrance and variable expressivity. Exome sequencing identified a rare missense variant in the ELF2 gene at chr4:g.140058846 C > T, c.10G > A, p.A4T which segregated in all affected patients. By using transduced BE (2)-M17 cells, we found that the mutated ELF2 (mt-ELF2) gene increased ATXN2 and reduced ELOVL5 gene expression, the causal genes of type 2 and type 38 spinocerebellar ataxias. Both, western blot and confocal microscopy confirmed an increase of ataxin-2 in BE(2)-M17 cells transduced with lentivirus expressing mt-ELF2 (CEE-mt-ELF2), which was not observed in cells transduced with lentivirus expressing wt-ELF2 (CEE-wt-ELF2). Moreover, we observed a significant decrease in the number and size of lipid droplets in the CEE-mt-ELF2-transduced BE (2)-M17 cells, but not in the CEE-wt-ELF2-transduced BE (2)-M17. Furthermore, changes in the expression of ELOVL5 could be related with the reduction of lipid droplets in BE (2)-M17 cells. This work supports that ELF2 gene regulates the expression of ATXN2 and ELOVL5 genes, and defines new molecular links in the pathophysiology of cerebellar ataxias.