In silico prediction of splice-altering single nucleotide variants in the human genome

Altered molecular and cellular mechanisms in KIF5A-associated neurodegenerative or neurodevelopmental disorders

Mutations targeting distinct domains of the neuron-specific kinesin KIF5A associate with different neurodegenerative/neurodevelopmental disorders, but the molecular bases of this clinical heterogeneity are unknown. We characterised five key mutants covering the whole spectrum of KIF5A-related phenotypes: spastic paraplegia (SPG, R17Q and R280C), Charcot-Marie-Tooth disease (CMT, R864*), amyotrophic lateral sclerosis (ALS, N999Vfs*40), and neonatal intractable myoclonus (NEIMY, C975Vfs*73) KIF5A mutants. CMT-R864*-KIF5A and ALS-N999Vfs*40-KIF5A showed impaired autoinhibition and peripheral localisation accompanied by altered mitochondrial distribution, suggesting transport competence disruption. ALS-N999Vfs*40-KIF5A formed SQSTM1/p62-positive inclusions sequestering WT-KIF5A, indicating a gain of toxic function. SPG-R17Q-KIF5A and ALS-N999Vfs*40-KIF5A evidenced a shorter half-life compared to WT-KIF5A, and proteasomal blockage determined their accumulation into detergent-insoluble inclusions. Interestingly, SPG-R280C-KIF5A and ALS-N999Vfs*40-KIF5A both competed for degradation with proteasomal substrates. Finally, NEIMY-C975Vfs*73-KIF5A displayed a similar, but more severe aberrant behaviour compared to ALS-N999Vfs*40-KIF5A; these two mutants share an abnormal tail but cause disorders on the opposite end of KIF5A-linked phenotypic spectrum. Thus, our observations support the pathogenicity of novel KIF5A mutants, highlight abnormalities of recurrent variants, and demonstrate that both unique and shared mechanisms underpin KIF5A-related diseases.

Genetic investigations on singleton school aged children reveal novel variants and new candidate genes for hearing loss

Hearing loss affects around 5% of the global population. Two preliminary studies have described genetic variants in sporadic individuals with hearing loss from Pakistan. Here we extend these studies to determine the spectrum of variants in a cohort of individuals with no previous history of hearing loss. Individuals with hearing loss born to consanguineous couples were identified from special schools. Audiograms were assessed. DNA from participants negative for GJB2 pathogenic variants was subjected to exome sequencing. Data were filtered to include variants with frequencies < 0.01 in the public databases. The effects of the missense variants on respective amino acids were analyzed by using PyMol software. Among the 44 participants, hearing loss was moderate for two individuals; 14 exhibited moderately-severe hearing loss while 25 had a severe degree of hearing loss. Hearing loss was reported to have been progressive in four participants and was currently profound in three participants. Variants were unambiguously identified in 17 genes, of which the majority affected SLC26A4. CDH23, MYO15A and OTOF were other significant contributors. Deleterious variants detected in two genes suggest new associations for hearing loss. Molecular characterization of hearing loss in our cohort revealed high genetic heterogeneity with a 75% diagnostic rate.

The recommendation of re-classification of variants of uncertain significance (VUS) in adult genetic disorders patients

Since variants of uncertain significance (VUS) reported in genetic testing cannot be acted upon clinically, this classification may delay or prohibit precise diagnosis and genetic counseling in adult genetic disorders patients. Large-scale analyses about qualitatively distinct lines of evidence used for VUS can make them re-classification more accurately. We analyzed 458 Chinese adult patients WES data, within 15 pathogenic evidence PS1, PS2, PM1, PM6 and PP4 were not used for VUS pathogenic classification, meanwhile the PP3, BP4, PP2 were used much more frequently. The PM2_Supporting was used most widely for all reported variants. There were also 31 null variants (nonsense, frameshift, canonical ±1 or 2 splice sites) which were probably the disease-causing variants of the patients were classified as VUS. By analyzed the evidence used for all VUS we recommend that appropriate genetic counseling, reliable releasing of in-house data, allele frequency comparison between case and control, expanded verification in patient family, co-segregation analysis and functional assays were urgent need to gather more evidence to reclassify VUS. We also found adult patients with nervous system disease were reported the most phenotype-associated VUS and the lower the phenotypic specificity, the more reported VUS. This result emphasized the importance of pretest genetic counseling which would make less reporting of VUS. Our result revealed the characteristics of the pathogenic classification evidence used for VUS in adult genetic disorders patients for the first time, recommend a rules-based process to evaluate the pathogenicity of VUS which could provide a strong basis for accurately evaluating the pathogenicity and clinical grade information of VUS. Meanwhile, we further expanded the genetic spectrum and improve the diagnostic rate of adult genetic disorders.

Dissecting CASK: Novel splice site variant associated with male MICPCH phenotype

CASK (MIM#300172), encoding a calcium/calmodulin-dependent serine protein kinase, is crucial for synaptic transmission and gene regulation during neural development. Pathogenic variants of CASK are known to cause several neurodevelopmental disorders, including X-linked intellectual disability and microcephaly with pontine and cerebellar hypoplasia (MICPCH). This study introduces a novel, de novo synonymous CASK variant (NM_001367721.1: c.1737G>A, p.(Glu579=)), discovered in a male patient diagnosed with MICPCH, characterized by microcephaly, developmental delay, visual impairment, and myoclonic seizures. The variant disrupts a donor splice-site at the end of exon 18. Transcriptomic analysis of blood identified 12 different CASK transcripts secondary to the synonymous variant. Nearly one third of these transcripts were predicted to result in nonsense mediated decay or protein degradation. Protein modeling revealed structural alterations in the PDZ functional domain of CASK, due to exon 18 deletion. Our findings highlight the utility of transcriptomic analysis in demonstrating the underlying disease mechanism in neurodevelopmental disorders.

Variant Impact Predictor database (VIPdb), version 2: trends from three decades of genetic variant impact predictors

BACKGROUND

Variant interpretation is essential for identifying patients' disease-causing genetic variants amongst the millions detected in their genomes. Hundreds of Variant Impact Predictors (VIPs), also known as Variant Effect Predictors (VEPs), have been developed for this purpose, with a variety of methodologies and goals. To facilitate the exploration of available VIP options, we have created the Variant Impact Predictor database (VIPdb).

RESULTS

The Variant Impact Predictor database (VIPdb) version 2 presents a collection of VIPs developed over the past three decades, summarizing their characteristics, ClinGen calibrated scores, CAGI assessment results, publication details, access information, and citation patterns. We previously summarized 217 VIPs and their features in VIPdb in 2019. Building upon this foundation, we identified and categorized an additional 190 VIPs, resulting in a total of 407 VIPs in VIPdb version 2. The majority of the VIPs have the capacity to predict the impacts of single nucleotide variants and nonsynonymous variants. More VIPs tailored to predict the impacts of insertions and deletions have been developed since the 2010s. In contrast, relatively few VIPs are dedicated to the prediction of splicing, structural, synonymous, and regulatory variants. The increasing rate of citations to VIPs reflects the ongoing growth in their use, and the evolving trends in citations reveal development in the field and individual methods.

CONCLUSIONS

VIPdb version 2 summarizes 407 VIPs and their features, potentially facilitating VIP exploration for various variant interpretation applications. VIPdb is available at  https://genomeinterpretation.org/vipdb.

Exome sequencing reveals genetic heterogeneity in consanguineous Pakistani families with neurodevelopmental and neuromuscular disorders

There remains a crucial need to address inequalities in genomic research and include populations from low- and middle-income countries (LMIC). Here we present eight consanguineous families from Pakistan, five with neurodevelopmental disorders (NDDs) and three with neuromuscular disorders (NMDs). Affected individuals were clinically characterized, and genetic variants were identified through exome sequencing (ES), followed by family segregation analysis. Affected individuals in six out of eight families (75%) carried homozygous variants that met ACMG criteria for being pathogenic (in the genes ADGRG1, METTL23, SPG11) or likely pathogenic (in the genes GPAA1, MFN2, SGSH). The remaining two families had homozygous candidate variants in the genes (AP4M1 and FAM126A) associated with phenotypes consistent with their clinical presentations, but the variants did not meet the criteria for pathogenicity and were hence classified as variants of unknown significance. Notably, the variants in ADGRG1, AP4M1, FAM126A, and SGSH did not have prior reports in the literature, demonstrating the importance of including diverse populations in genomic studies. We provide clinical phenotyping along with analyses of ES data that support the utility of ES in making accurate molecular diagnoses in these patients, as well as in unearthing novel variants in known disease-causing genes in underrepresented populations from LMIC.

Identification of deleterious variants in nine polycystic kidney disease affected families

Polycystic kidney disease (PKD) is common genetic renal disorder. In present study, we performed WES to identify pathogenic variant in nine families including 26 patients with PKD and 19 unaffected members. The eight pathogenic variants were identified in known PKD associated genes including PKD1 (n = 6), PKD2 (n = 1), and OFD1 (n = 1) in eight families. There is one missense, one stopgain, two non-frameshifts, two canonical splicing variants, three frameshift variants and one potential non-canonical splicing variant (NCSV) in 8 families. The six variants were novel variants and not reported in ClinVar database. In addition, the compound heterozygous variants in PKHD1 were identified including one frameshift variants (PKHD1: NM_138694.4, c.9841del, p.S3281Lfs*4) and one non-canonical splicing variant (PKHD1: NM_138694.4, c.6332 + 40A > G) which were defined as deleterious variant by four splicing prediction tools (CADD-splice, SpliceAI, Spliceogen, Squirl). We used the minigene method to validate whether the prioritized potential NSCVs disrupt the typical mRNA splicing process and found abnormally larger PCR production of minigene carrying potential NCSV comparing to wild-type minigene. Sanger sequencing confirmed the 39-bp insertion of intron 38 between exon 38 and exon 39, which results in non-frameshift and 13 amino acid insertions. In conclusion, our study expands the variant spectrum and highlight the important role of non-canonical splicing variant in PKD.

Single-Base Substitution Causing Dual-Exon Skipping Event in PKD2 Gene: Unusual Molecular Finding from Exome Sequencing in a Patient with Autosomal Dominant Polycystic Kidney Disease

Background: Pathogenic variants in the Polycystic Kidney Disease 2 (PKD2) gene are associated with Autosomal Dominant Polycystic Kidney Disease (ADPKD) in approximately 30% of cases. In recent years, the high-throughput sequencing techniques have significantly increased the number of variants identified in affected patients. Here, we described the peculiar effect of a PKD2 splicing variant, the c.1717-2A>G, identified in an Italian male patient with ADPKD. This variant led to the unusual and rare skipping of two consecutive exons, causing a large in-frame deletion. Methods: The genetic evaluation of the patient was performed using the Next-Generation Sequencing (NGS) assay Clinical Exome Solution® (SOPHiA Genetics). Bioinformatics analysis was performed using the SOPHiA DDM platform (SOPHiA Genetics). Prediction of pathogenicity was carried out by integrating several in silico tools. RNA evaluation was performed to test the effect of the variant on the PKD2 splicing using a Reverse-Transcription PCR coupled with cDNA sequencing. Results: NGS revealed the presence of the PKD2 c.1717-2A>G variant that lies in the canonical splice site of intron 7. This rare variant was predicted to have a significant impact on the splicing, proved by the RNA-based analysis. We identified the presence of a transcript characterised by the simultaneous skipping of exons 8 and 9, with a retained reading frame and the merging of exons 7-10. Conclusions: We described for the first time a dual-exon skip event related to the presence of a single-base substitution in the PKD2 gene in an ADPKD-affected patient. We assumed that the molecular basis of such a rare mechanism lies in the specific order of intron removal. The finding represents novel evidence of an alternative and unusual splicing mechanism in the PKD2 gene, adding insights to the pathogenesis of the ADPKD.

RNA variant assessment using transactivation and transdifferentiation

Understanding the impact of splicing and nonsense variants on RNA is crucial for the resolution of variant classification as well as their suitability for precision medicine interventions. This is primarily enabled through RNA studies involving transcriptomics followed by targeted assays using RNA isolated from clinically accessible tissues (CATs) such as blood or skin of affected individuals. Insufficient disease gene expression in CATs does however pose a major barrier to RNA based investigations, which we show is relevant to 1,436 Mendelian disease genes. We term these "silent" Mendelian genes (SMGs), the largest portion (36%) of which are associated with neurological disorders. We developed two approaches to induce SMG expression in human dermal fibroblasts (HDFs) to overcome this limitation, including CRISPR-activation-based gene transactivation and fibroblast-to-neuron transdifferentiation. Initial transactivation screens involving 40 SMGs stimulated our development of a highly multiplexed transactivation system culminating in the 6- to 90,000-fold induction of expression of 20/20 (100%) SMGs tested in HDFs. Transdifferentiation of HDFs directly to neurons led to expression of 193/516 (37.4%) of SMGs implicated in neurological disease. The magnitude and isoform diversity of SMG expression following either transactivation or transdifferentiation was comparable to clinically relevant tissues. We apply transdifferentiation and/or gene transactivation combined with short- and long-read RNA sequencing to investigate the impact that variants in USH2A, SCN1A, DMD, and PAK3 have on RNA using HDFs derived from affected individuals. Transactivation and transdifferentiation represent rapid, scalable functional genomic solutions to investigate variants impacting SMGs in the patient cell and genomic context.

Genomic Landscape of Susceptibility to Severe COVID-19 in the Slovenian Population

Determining the genetic contribution of susceptibility to severe SARS-CoV-2 infection outcomes is important for public health measures and individualized treatment. Through intense research on this topic, several hundred genes have been implicated as possibly contributing to the severe infection phenotype(s); however, the findings are complex and appear to be population-dependent. We aimed to determine the contribution of human rare genetic variants associated with a severe outcome of SARS-CoV-2 infections and their burden in the Slovenian population. A panel of 517 genes associated with severe SARS-CoV-2 infection were obtained by combining an extensive review of the literature, target genes identified by the COVID-19 Host Genetic Initiative, and the curated Research COVID-19 associated genes from PanelApp, England Genomics. Whole genome sequencing was performed using PCR-free WGS on DNA from 60 patients hospitalized due to severe COVID-19 disease, and the identified rare genomic variants were analyzed and classified according to the ACMG criteria. Background prevalence in the general Slovenian population was determined by comparison with sequencing data from 8025 individuals included in the Slovenian genomic database (SGDB). Results show that several rare pathogenic/likely pathogenic genomic variants in genes CFTR, MASP2, MEFV, TNFRSF13B, and RNASEL likely contribute to the severe infection outcomes in our patient cohort. These results represent an insight into the Slovenian genomic diversity associated with a severe COVID-19 outcome.

Systematic rare variant analyses identify RAB32 as a susceptibility gene for familial Parkinson's disease

Despite substantial progress, causal variants are identified only for a minority of familial Parkinson's disease (PD) cases, leaving high-risk pathogenic variants unidentified1,2. To identify such variants, we uniformly processed exome sequencing data of 2,184 index familial PD cases and 69,775 controls. Exome-wide analyses converged on RAB32 as a novel PD gene identifying c.213C > G/p.S71R as a high-risk variant presenting in ~0.7% of familial PD cases while observed in only 0.004% of controls (odds ratio of 65.5). This variant was confirmed in all cases via Sanger sequencing and segregated with PD in three families. RAB32 encodes a small GTPase known to interact with LRRK2 (refs. 3,4). Functional analyses showed that RAB32 S71R increases LRRK2 kinase activity, as indicated by increased autophosphorylation of LRRK2 S1292. Here our results implicate mutant RAB32 in a key pathological mechanism in PD-LRRK2 kinase activity5-7-and thus provide novel insights into the mechanistic connections between RAB family biology, LRRK2 and PD risk.

Molecular diagnosis of cystic fibrosis by RNA obtained from nasal epithelial cells

BACKGROUND

The diagnosis of cystic fibrosis (CF) is established when characteristic clinical signs are coupled with biallelic CFTR pathogenic variants. No previously reported non-canonical splice site variants have to be considered as variants of uncertain significance unless their effect on splicing has been validated.

METHODS

Two variants identified by next-generation sequencing were evaluated. We assayed their effects on splicing employing RNA analysis and real-time expression quantification from RNA obtained from the nasal epithelial cells of a patient with clinically suspected CF and of two patients with milder phenotypes (CFTR-related disorders).

RESULTS

The variant c.164+2dup causes skipping of exon 2 (p.(Ser18_Glu54del)) and exon 2 plus 3 (p.(Ser18Argfs*16)) in CFTR mRNA. Exon 2 expression in the patient heterozygous for c.164+2dup was decreased to 7 % of the exon 2 expression in the controls. The synonymous variant c.1584G>A causes a partial skipping of exon 11. The exon 11 expression in the two patients heterozygous for this variant was 22 % and 42 % of that of the controls, respectively.

CONCLUSION

We conclude that variant c.164+2dup affects mRNA processing and can be considered a CF-causing variant. The results of the functional assay also showed that the p.(Glu528=) variant, usually categorized as a neutral variant based on epidemiological data, partially affects mRNA processing in our patients. This finding would allow us to reclassify the variant as a CFTR-related variant with incomplete penetrance. RNA obtained from nasal epithelial cells is an easy and accurate tool for CFTR functional studies in patients with unclassified splice variants.

Variant Impact Predictor database (VIPdb), version 2: Trends from 25 years of genetic variant impact predictors

Background

Variant interpretation is essential for identifying patients' disease-causing genetic variants amongst the millions detected in their genomes. Hundreds of Variant Impact Predictors (VIPs), also known as Variant Effect Predictors (VEPs), have been developed for this purpose, with a variety of methodologies and goals. To facilitate the exploration of available VIP options, we have created the Variant Impact Predictor database (VIPdb).

Results

The Variant Impact Predictor database (VIPdb) version 2 presents a collection of VIPs developed over the past 25 years, summarizing their characteristics, ClinGen calibrated scores, CAGI assessment results, publication details, access information, and citation patterns. We previously summarized 217 VIPs and their features in VIPdb in 2019. Building upon this foundation, we identified and categorized an additional 186 VIPs, resulting in a total of 403 VIPs in VIPdb version 2. The majority of the VIPs have the capacity to predict the impacts of single nucleotide variants and nonsynonymous variants. More VIPs tailored to predict the impacts of insertions and deletions have been developed since the 2010s. In contrast, relatively few VIPs are dedicated to the prediction of splicing, structural, synonymous, and regulatory variants. The increasing rate of citations to VIPs reflects the ongoing growth in their use, and the evolving trends in citations reveal development in the field and individual methods.

Conclusions

VIPdb version 2 summarizes 403 VIPs and their features, potentially facilitating VIP exploration for various variant interpretation applications.

Availability

VIPdb version 2 is available at https://genomeinterpretation.org/vipdb.

Clinical and Genetic Findings in a Cohort of Patients with PRPF31-Associated Retinal Dystrophy

PURPOSE

The purpose of our study was to assess the phenotypic and genotypic spectrum in a large cohort of patients with PRPF31-associated retinal dystrophy.

DESIGN

Retrospective cohort study.

METHODS

In this retrospective chart review study, we collected cross-sectional data on the phenotype and genotype of patients with PRPF31-associated retinal dystrophy from the clinics for inherited retinal dystrophies at the University of Tuebingen and the local RetDis database and biobank. Patients underwent thorough ophthalmological examinations and genetic testing.

RESULTS

Eighty-six patients from 61 families were available for clinical assessment, while genomic DNA was available for 111 individuals (index patients and family members). Fifty-three different disease-associated variants were observed in our cohort. Point mutations were the most common class. All but two patients exhibited features of a typical Retinitis pigmentosa (RP). One patient showed a cone-rod dystrophy pattern. One mutation carrier revealed no signs of a retinal dystrophy. There was a statistically significant better visual acuity for patients with large deletions in the 20-39 age group. Cystoid macular edema was common in those with preserved central retina and showed an association with female sex.

CONCLUSION

Our study confirms high phenotypic variability in disease onset and age at which legal blindness is reached in PRPF31-associated RP. Non-penetrance is commonly documented in family history, although poorly represented in our study, possibly indicating that true asymptomatic mutation carriers are rare if followed-up over lifetime with thorough ophthalmologic workup.

A novel variation in DEPDC5 causing familial focal epilepsy with variable foci

Background

Disheveled, EGL-10, and pleckstrin (DEP) domain-containing protein 5 (DEPDC5) is a component of GTPase-activating protein (GAP) activity toward the RAG complex 1 (GATOR1) protein, which is an inhibitor of the amino acid-sensing branch of the mammalian target of rapamycin complex 1 (mTORC1) pathway. GATOR1 complex variations were reported to correlate with familial focal epilepsy with variable foci (FFEVF). With the wide application of whole exome sequencing (WES), more and more variations in DEPDC5 were uncovered in FFEVF families.

Methods

A family with a proband diagnosed with familial focal epilepsy with variable foci (FFEVF) was involved in this study. Whole exome sequencing (WES) was performed in the proband, and Sanger sequencing was used to confirm the variation carrying status of the family members. Mini-gene splicing assay was performed to validate the effect on the alternative splicing of the variation.

Results

A novel variant, c.1217 + 2T>A, in DEPDC5 was identified by WES in the proband. This splicing variant that occurred at the 5' end of intron 17 was confirmed by mini-gene splicing assays, which impacted alternative splicing and led to the inclusion of an intron fragment. The analysis of the transcribed mRNA sequence indicates that the translation of the protein is terminated prematurely, which is very likely to result in the loss of function of the protein and lead to the occurrence of FFEVF.

Conclusion

The results suggest that c.1217 + 2T>A variations in DEPDC5 might be the genetic etiology for FFEVF in this pedigree. This finding expands the genotype spectrum of FFEVF and provides new etiological information for FFEVF.

Whole-exome sequencing in a cohort of Chinese patients with isolated cervical dystonia

Background

Dystonia is a kind of movement disorder but its pathophysiological mechanisms are still largely unknown. Recent evidence reveals that genetical defects may play important roles in the pathogenesis of dystonia.

Objectives and Methods

-To explore possible causative genes in Chinese dystonia patients, DNA samples from 42 sporadic patients with isolated cervical dystonia were subjected to whole-exome sequencing. Rare deleterious variants associated with dystonia phenotype were screened out and then classified according to the American College of Medical Genetics and Genomics (ACMG) criteria. Phenolyzer was used for analyzing the most probable candidates correlated with dystonia phenotype, and SWISS-MODEL server was for predicting the 3D structures of variant proteins.

Results

Among 42 patients (17 male and 25 female) recruited, a total of 36 potentially deleterious variants of dystonia-associated genes were found in 30 patients (30/42, 71.4 %). Four disease-causing variants including a pathogenic variant in PLA2G6 (c.797G > C) and three likely pathogenic variants in DCTN1 (c.73C > T), SPR (c.1A > C) and TH (c.56C > G) were found in four patients separately. Other 32 variants were classified as uncertain significance in 26 patients. Phenolyzer prioritized genes TH, PLA2G6 and DCTN1 as the most probable candidates correlated with dystonia phenotype. Although 3D prediction of DCTN1 and PLA2G6 variant proteins detected no obvious structural alterations, the mutation in DCTN1 (c.73C > T:p.Arg25Trp) was closely adjacent to its key functional domain.

Conclusion

Our whole-exome sequencing results identified a novel variant in DCTN1 in sporadic Chinese patients with isolated cervical dystonia, which however, needs our further study on its exact role in dystonia pathogenesis.

Marked Response to Nivolumab by a Patient With SMARCA4-Deficient Undifferentiated Urothelial Carcinoma Showing High PD-L1 Expression: A Case Report

BACKGROUND

SMARCA4 is a component gene of the SWI/SNF (SWItch/Sucrose NonFermentable) chromatin remodeling complex; undifferentiated tumors associated with its functional deletion have been described in several organs. However, no established treatment for these tumors currently exists.

CASE

In this study, we report a case of a SMARCA4-deficient undifferentiated urothelial carcinoma with high PD-L1 expression that was effectively treated with nivolumab after early relapse following treatment for non-invasive bladder cancer. The histological morphology of the rhabdoid-like undifferentiated tumor of unknown primary led us to suspect a SWI/SNF-deficient tumor, and subsequent immunostaining led to the diagnosis of a SMARCA4-deficient undifferentiated tumor. This effort also led to the identification of the developmental origin of this SMARCA4-deficient undifferentiated tumor as a non-invasive bladder cancer. We also carried out a detailed immune phenotypic assay on peripheral T cells. In brief, a phenotypic change of CD8+T cells from naive to terminally differentiated effector memory cells was observed.

CONCLUSION

Regardless of the organ of cancer origin or cancer type, SWI/SNF-deficient tumors should be suspected in undifferentiated and dedifferentiated tumors, and immune checkpoint inhibitors may be considered as a promising treatment option for this type of tumor. The pathogenesis of SMARCA4-deficient anaplastic tumors awaits further elucidation for therapeutic development.

Loss-of-function variants in UBAP1L cause autosomal recessive retinal degeneration

PURPOSE

Inherited retinal diseases (IRDs) are a group of monogenic conditions that can lead to progressive blindness. Their missing heritability is still considerable, due in part to the presence of disease genes that await molecular identification. The purpose of this work was to identify novel genetic associations with IRDs.

METHODS

Patients underwent a comprehensive ophthalmological evaluation using standard-of-care tests, such as detailed retinal imaging (macular optical coherence tomography and short-wavelength fundus autofluorescence) and electrophysiological testing. Exome and genome sequencing, as well as computer-assisted data analysis were used for genotyping and detection of DNA variants. A minigene-driven splicing assay was performed to validate the deleterious effects of 1 of such variants.

RESULTS

We identified 8 unrelated families from Hungary, the United States, Israel, and The Netherlands with members presenting with a form of autosomal recessive and nonsyndromic retinal degeneration, predominantly described as rod-cone dystrophy but also including cases of cone/cone-rod dystrophy. Age of disease onset was very variable, with some patients experiencing first symptoms during their fourth decade of life or later. Myopia greater than 5 diopters was present in 5 of 7 cases with available refractive data, and retinal detachment was reported in 2 cases. All ascertained patients carried biallelic loss-of-function variants in UBAP1L (HGNC: 40028), a gene with unknown function and with homologies to UBAP1, encoding a protein involved in ubiquitin metabolism. One of these pathogenic variants, the intronic NM_001163692.2:c.910-7G>A substitution, was identified in 5 unrelated families. Minigene-driven splicing assays in HEK293T cells confirmed that this DNA change is responsible for the creation of a new acceptor splice site, resulting in aberrant splicing.

CONCLUSION

We identified UBAP1L as a novel IRD gene. Although its function is currently unknown, UBAP1L is almost exclusively expressed in photoreceptors and the retinal pigment epithelium, hence possibly explaining the link between pathogenic variants in this gene and an ocular phenotype.

An interpretable model of pre-mRNA splicing for animal and plant genes

Pre-mRNA splicing is a fundamental step in gene expression, conserved across eukaryotes, in which the spliceosome recognizes motifs at the 3' and 5' splice sites (SSs), excises introns, and ligates exons. SS recognition and pairing is often influenced by protein splicing factors (SFs) that bind to splicing regulatory elements (SREs). Here, we describe SMsplice, a fully interpretable model of pre-mRNA splicing that combines models of core SS motifs, SREs, and exonic and intronic length preferences. We learn models that predict SS locations with 83 to 86% accuracy in fish, insects, and plants and about 70% in mammals. Learned SRE motifs include both known SF binding motifs and unfamiliar motifs, and both motif classes are supported by genetic analyses. Our comparisons across species highlight similarities between non-mammals, increased reliance on intronic SREs in plant splicing, and a greater reliance on SREs in mammalian splicing.

Comprehensive Genetic Evaluation in Patients with Special Reference to Late-Onset Sensorineural Hearing Loss

Hearing loss (HL) is a common and multi-complex etiological deficit that can occur at any age and can be caused by genetic variants, aging, toxic drugs, noise, injury, viral infection, and other factors. Recently, a high incidence of genetic etiologies in congenital HL has been reported, and the usefulness of genetic testing has been widely accepted in congenital-onset or early-onset HL. In contrast, there have been few comprehensive reports on the relationship between late-onset HL and genetic causes. In this study, we performed next-generation sequencing analysis for 91 HL patients mainly consisting of late-onset HL patients. As a result, we identified 23 possibly disease-causing variants from 29 probands, affording a diagnostic rate for this study of 31.9%. The highest diagnostic rate was observed in the congenital/early-onset group (42.9%), followed by the juvenile/young adult-onset group (31.7%), and the middle-aged/aged-onset group (21.4%). The diagnostic ratio decreased with age; however, genetic etiologies were involved to a considerable degree even in late-onset HL. In particular, the responsible gene variants were found in 19 (55.9%) of 34 patients with a familial history and progressive HL. Therefore, this phenotype is considered to be a good candidate for genetic evaluation based on this diagnostic panel.

The prevalence and clinical features of MYO7A-related hearing loss including DFNA11, DFNB2 and USH1B

The MYO7A gene is known to be responsible for both syndromic hearing loss (Usher syndrome type1B:USH1B) and non-syndromic hearing loss including autosomal dominant and autosomal recessive inheritance (DFNA11, DFNB2). However, the prevalence and detailed clinical features of MYO7A-associated hearing loss across a large population remain unclear. In this study, we conducted next-generation sequencing analysis for a large cohort of 10,042 Japanese hearing loss patients. As a result, 137 patients were identified with MYO7A-associated hearing loss so that the prevalence among Japanese hearing loss patients was 1.36%. We identified 70 disease-causing candidate variants in this study, with 36 of them being novel variants. All variants identified in autosomal dominant cases were missense or in-frame deletion variants. Among the autosomal recessive cases, all patients had at least one missense variant. On the other hand, in patients with Usher syndrome, almost half of the patients carried biallelic null variants (nonsense, splicing, and frameshift variants). Most of the autosomal dominant cases showed late-onset progressive hearing loss. On the other hand, cases with autosomal recessive inheritance or Usher syndrome showed congenital or early-onset hearing loss. The visual symptoms in the Usher syndrome cases developed between age 5-15, and the condition was diagnosed at about 6-15 years of age.

A Rare Inherited Bone Marrow Failure Syndrome Disclosed by Reanalysis of the Exome Data of a Patient Evaluated for Cytopenia and Dysmorphic Features

BACKGROUND

Multisystemic findings of inherited bone marrow failure syndromes may cause difficulty in diagnosis. Exome sequencing (ES) helps to define the etiology of rare diseases and reanalysis offers a valuable new diagnostic approach. Herein, we present the clinical and molecular characteristics of a girl who was referred for cytopenia and frequent infections.

CASE REPORT

A 5-year-old girl with cytopenia, dysmorphism, short stature, developmental delay, and myopia was referred for genetic counseling. Reanalysis of the ES data revealed a homozygous splice-site variant in the DNAJC21 (NM_001012339.3:c.983+1G>A), causing Shwachman-Diamond Syndrome (SDS). It was shown by the RNA sequencing that exon 7 was skipped, causing an 88-nucleotide deletion.

CONCLUSIONS

Precise genetic diagnosis enables genetic counseling and improves patient management by avoiding inappropriate treatment and unnecessary testing. This report would contribute to the clinical and molecular understanding of this rare type of SDS caused by DNAJC21 variants and expand the phenotypic features of this condition.

A spectrum of TP63-related disorders with eight affected individuals in five unrelated families

TP63-related disdorders broadly involve varying combinations of ectodermal dysplasia (sparse hair, hypohydrosis, tooth abnormalities, nail dysplasia), cleft lip/palate, acromelic malformation, split-hand/foot malformation/syndactyly, ankyloblepharon filiforme adnatum, lacrimal duct obstruction, hypopigmentation, and hypoplastic breasts and/or nipples. TP63-related disorders are associated with heterozygous pathogenic variants in TP63 and include seven overlapping phenotypes; Ankyloblepharon-ectodermal defects-cleft lip/palate syndrome (AEC), Ectrodactyly-ectodermal dysplasia-cleft lip/palate syndrome 3 (EEC3), Limb-mammary syndrome (LMS), Acro-dermo-ungual-lacrimal-tooth syndrome (ADULT), Rapp-Hodgkin syndrome (RHS), Split-hand/foot malformation 4 (SHFM4), and Orofacial cleft 8. We report on five unrelated families with 8 affected individuals in which the probands presented with varying combinations of ectodermal dysplasia, cleft lip/palate, split-hand/foot malformation, lacrimal duct obstruction, and ankyloblepharon filiforme adnatum. The clinical diagnosis involved AEC syndrome (2 patients), EEC3 syndrome (2 patients), and a yet hitherto unclassified TP63-related disorder. Sanger sequence analysis of the TP63 gene was performed revealing five different variants among which four were novel and three were de novo. The identificated TP63 variants co-segregated with the other affected individuals in the families. The abnormalities of ectoderm derived structures including hair, nails, sweat glands, and teeth should alert the physician to the possibility of TP63-related disorders particularly in the presence of orofacial clefting.

Comprehensive whole-genome sequence analyses provide insights into the genomic architecture of cerebral palsy

We performed whole-genome sequencing (WGS) in 327 children with cerebral palsy (CP) and their biological parents. We classified 37 of 327 (11.3%) children as having pathogenic/likely pathogenic (P/LP) variants and 58 of 327 (17.7%) as having variants of uncertain significance. Multiple classes of P/LP variants included single-nucleotide variants (SNVs)/indels (6.7%), copy number variations (3.4%) and mitochondrial mutations (1.5%). The COL4A1 gene had the most P/LP SNVs. We also analyzed two pediatric control cohorts (n = 203 trios and n = 89 sib-pair families) to provide a baseline for de novo mutation rates and genetic burden analyses, the latter of which demonstrated associations between de novo deleterious variants and genes related to the nervous system. An enrichment analysis revealed previously undescribed plausible candidate CP genes (SMOC1, KDM5B, BCL11A and CYP51A1). A multifactorial CP risk profile and substantial presence of P/LP variants combine to support WGS in the diagnostic work-up across all CP and related phenotypes.

Evaluating the efficacy of a long-read sequencing-based approach in the clinical diagnosis of neonatal congenital adrenocortical hyperplasia

Congenital adrenal hyperplasia (CAH) is a group of autosomal recessive disorders predominantly characterized by impaired corticosteroid synthesis. Clinical phenotypes include hypoadrenocorticism, electrolyte disturbances, abnormal gonadal development, and short stature, of which severe hyponadrenocorticism and salt wasting can be life-threatening. Genetic analysis can help in the clinical diagnosis of CAH. However, the 21-OHD-causing gene CYP21A2 is arranged in tandem with the highly homologous CYP21A1P pseudogene, making it difficult to determine the exact genotypes using the traditional method of multiplex ligation-dependent probe amplification (MLPA) plus Sanger sequencing or next-generation sequencing (NGS). We applied a long-read sequencing-based approach termed comprehensive analysis of CAH (CACAH) to 48 newborns with CAH that were diagnosed by clinical features and the traditional MLPA plus Sanger sequencing method for retrospective analysis, to evaluate its efficacy in the clinical diagnosis of neonatal CAH. Compared with the MLPA plus Sanger sequencing method, CACAH showed 100 % consistency in detecting SNV/indel variants located in exons and exon-intron boundary regions of CAH-related genes. It can directly determine the cis-trans relationship without the need to analyze parental genotypes, which reduces the time to diagnosis. Moreover, CACAH was able to distinguish different CYP21A1P/CYP21A2 and TNXA/TNXB chimeras, and detect additional variants (CYP21A2 variants c.-121C > T, c.*13G > A, c.*52C > T, c.*440C > T, c.*443 T > C, and TNXB variants c.12463 + 2 T > C, c.12204 + 5G > A). We also identified the TNXB variant c.11435_11524 + 30del alone instead of as a part of the TNXA/TNXB-CH-1 chimera in two newborns, which might be introduced by gene conversion. All of these characteristics enabled clinicians to better explain the phenotype of subjects and manage them more effectively. CACAH has a great advantage over the traditional MLPA and Sanger sequencing methods, showing substantial potential in the genetic diagnosis and screening of neonatal CAH.

Validation of targeted next-generation sequencing panels in a cohort of Polish patients with epilepsy: assessing variable performance across clinical endophenotypes and uncovering novel genetic variants

Introduction

Targeted Next-Generation Sequencing Panels (TNGSP) have become a standard in global clinical practice. Instead of questioning the necessity of next-generation sequencing in epilepsy patients, contemporary large-scale research focuses on factors such as the size of TNGSP, the comparative advantages of exome or genome-wide sequencing over TNGSP, and the impact of clinical, electrophysiological, and demographic variables on genetic test performance. This study aims to elucidate the demographic and clinical factors influencing the performance of TNGSP in 138 Polish patients with epilepsy, recognizing the pivotal role of genetic testing in guiding patient management and therapy.

Methods

A retrospective analysis was conducted on patients from a genetic clinic in Poznań, Poland, who underwent commercial gene panel studies at Invitae Corporation (USA) between 2020 and 2022. Patient groups were defined based on the age of onset of the first epileptic seizures, seizure type, gender, fever dependence of seizures, presence of intellectual disability or developmental delay, abnormalities in MRI, and the presence of dysmorphic features or congenital malformations. Seizure classification followed the 2017 ILAE criteria.

Results

Among the 138 patients, 30 (21.7%) exhibited a pathogenic or likely pathogenic variant, with a distribution of 20.7% in males and 22.5% in females. Diagnostic performance correlated with the patient's age at the onset of the first seizure and the type of seizure. Predominant variants were identified in the SCN1A, PRRT2, CDKL5, DEPDC5, TSC2, and SLC2A1 genes. Additionally, 12 genes (CACNA1A, SCN2A, GRIN2A, KCNQ2, CHD2, DYNC1H1, NEXMIF, SCN1B, DDX3X, EEF1A2, NPRL3, UBE3A) exhibited single instances of damage. Notably, novel variants were discovered in DEPDC5, SCN1A, TSC2, CDKL5, NPRL3, DYNC1H1, CHD2, and DDX3X.

Discussion

Identified variants were present in genes previously recognized in both European and non-European populations. A thorough examination of Variants of Uncertain Significance (VUSs), specifically focusing on gene copy number changes, may unveil more extensive chromosomal aberrations. The relatively frequent occurrence of pathological variants in X chromosome-linked genes in girls warrants further investigation, challenging the prevailing notion of male predominance in X-linked epilepsy.

Combining a prioritization strategy and functional studies nominates 5'UTR variants underlying inherited retinal disease

BACKGROUND

5' untranslated regions (5'UTRs) are essential modulators of protein translation. Predicting the impact of 5'UTR variants is challenging and rarely performed in routine diagnostics. Here, we present a combined approach of a comprehensive prioritization strategy and functional assays to evaluate 5'UTR variation in two large cohorts of patients with inherited retinal diseases (IRDs).

METHODS

We performed an isoform-level re-analysis of retinal RNA-seq data to identify the protein-coding transcripts of 378 IRD genes with highest expression in retina. We evaluated the coverage of their 5'UTRs by different whole exome sequencing (WES) kits. The selected 5'UTRs were analyzed in whole genome sequencing (WGS) and WES data from IRD sub-cohorts from the 100,000 Genomes Project (n = 2397 WGS) and an in-house database (n = 1682 WES), respectively. Identified variants were annotated for 5'UTR-relevant features and classified into seven categories based on their predicted functional consequence. We developed a variant prioritization strategy by integrating population frequency, specific criteria for each category, and family and phenotypic data. A selection of candidate variants underwent functional validation using diverse approaches.

RESULTS

Isoform-level re-quantification of retinal gene expression revealed 76 IRD genes with a non-canonical retina-enriched isoform, of which 20 display a fully distinct 5'UTR compared to that of their canonical isoform. Depending on the probe design, 3-20% of IRD genes have 5'UTRs fully captured by WES. After analyzing these regions in both cohorts, we prioritized 11 (likely) pathogenic variants in 10 genes (ARL3, MERTK, NDP, NMNAT1, NPHP4, PAX6, PRPF31, PRPF4, RDH12, RD3), of which 7 were novel. Functional analyses further supported the pathogenicity of three variants. Mis-splicing was demonstrated for the PRPF31:c.-9+1G>T variant. The MERTK:c.-125G>A variant, overlapping a transcriptional start site, was shown to significantly reduce both luciferase mRNA levels and activity. The RDH12:c.-123C>T variant was found in cis with the hypomorphic RDH12:c.701G>A (p.Arg234His) variant in 11 patients. This 5'UTR variant, predicted to introduce an upstream open reading frame, was shown to result in reduced RDH12 protein but unaltered mRNA levels.

CONCLUSIONS

This study demonstrates the importance of 5'UTR variants implicated in IRDs and provides a systematic approach for 5'UTR annotation and validation that is applicable to other inherited diseases.

Using multi-scale genomics to associate poorly annotated genes with rare diseases

BACKGROUND

Next-generation sequencing (NGS) has significantly transformed the landscape of identifying disease-causing genes associated with genetic disorders. However, a substantial portion of sequenced patients remains undiagnosed. This may be attributed not only to the challenges posed by harder-to-detect variants, such as non-coding and structural variations but also to the existence of variants in genes not previously associated with the patient's clinical phenotype. This study introduces EvORanker, an algorithm that integrates unbiased data from 1,028 eukaryotic genomes to link mutated genes to clinical phenotypes.

METHODS

EvORanker utilizes clinical data, multi-scale phylogenetic profiling, and other omics data to prioritize disease-associated genes. It was evaluated on solved exomes and simulated genomes, compared with existing methods, and applied to 6260 knockout genes with mouse phenotypes lacking human associations. Additionally, EvORanker was made accessible as a user-friendly web tool.

RESULTS

In the analyzed exomic cohort, EvORanker accurately identified the "true" disease gene as the top candidate in 69% of cases and within the top 5 candidates in 95% of cases, consistent with results from the simulated dataset. Notably, EvORanker outperformed existing methods, particularly for poorly annotated genes. In the case of the 6260 knockout genes with mouse phenotypes, EvORanker linked 41% of these genes to observed human disease phenotypes. Furthermore, in two unsolved cases, EvORanker successfully identified DLGAP2 and LPCAT3 as disease candidates for previously uncharacterized genetic syndromes.

CONCLUSIONS

We highlight clade-based phylogenetic profiling as a powerful systematic approach for prioritizing potential disease genes. Our study showcases the efficacy of EvORanker in associating poorly annotated genes to disease phenotypes observed in patients. The EvORanker server is freely available at https://ccanavati.shinyapps.io/EvORanker/ .

Predicting the impact of rare variants on RNA splicing in CAGI6

Variants which disrupt splicing are a frequent cause of rare disease that have been under-ascertained clinically. Accurate and efficient methods to predict a variant's impact on splicing are needed to interpret the growing number of variants of unknown significance (VUS) identified by exome and genome sequencing. Here, we present the results of the CAGI6 Splicing VUS challenge, which invited predictions of the splicing impact of 56 variants ascertained clinically and functionally validated to determine splicing impact. The performance of 12 prediction methods, along with SpliceAI and CADD, was compared on the 56 functionally validated variants. The maximum accuracy achieved was 82% from two different approaches, one weighting SpliceAI scores by minor allele frequency, and one applying the recently published Splicing Prediction Pipeline (SPiP). SPiP performed optimally in terms of sensitivity, while an ensemble method combining multiple prediction tools and information from databases exceeded all others for specificity. Several challenge methods equalled or exceeded the performance of SpliceAI, with ultimate choice of prediction method likely to depend on experimental or clinical aims. One quarter of the variants were incorrectly predicted by at least 50% of the methods, highlighting the need for further improvements to splicing prediction methods for successful clinical application.

A novel missense COL9A3 variant in a pedigree with multiple lumbar disc herniation

Trp3 allele in COL9A3 gene has been widely studied in populations with intervertebral disc disease. We identified a novel pathogenic variant in COL9A3 gene in a pedigree with multiple lumbar disc herniation (LDH). The proband was a 14-year-old boy who developed LDH at the L4/5 and L5/S1 spinal segments. His father, paternal aunt and grandfather were diagnosed with LDH at an age of 35, 30 and 23, respectively. By applying whole exome sequencing, a heterozygous missense variant (c.1150C > T, p.Arg384Trp) in COL9A3 was identified. According to the ACMG guidelines, this variant is predicted to be pathogenic. In addition, prediction tools found COL9A3 protein of this variant a reduced stability, some changed charge properties, and an altered spatial conformation. Findings expanded the mutational spectrum of LDH and contributed to the understanding of COL9A3 in the pathogenesis of LDH.

Prioritizing de novo potential non-canonical splicing variants in neurodevelopmental disorders

BACKGROUND

Genomic variants outside of the canonical splicing site (±2) may generate abnormal mRNA splicing, which are defined as non-canonical splicing variants (NCSVs). However, the clinical interpretation of NCSVs in neurodevelopmental disorders (NDDs) is largely unknown.

METHODS

We investigated the contribution of NCSVs to NDDs from 345,787 de novo variants (DNVs) in 47,574 patients with NDDs. We performed functional enrichment and protein-protein interaction analysis to assess the association between genes carrying prioritised NCSVs and NDDs. Minigene was used to validate the impact of NCSVs on mRNA splicing.

FINDINGS

We observed significantly more NCSVs (p = 0.02, odds ratio [OR] = 2.05) among patients with NDD than in controls. Both canonical splicing variants (CSVs) and NCSVs contributed to an equal proportion of patients with NDD (0.76% vs. 0.82%). The candidate genes carrying NCSVs were associated with glutamatergic synapse and chromatin remodelling. Minigene successfully validated 59 of 79 (74.68%) NCSVs that led to abnormal splicing in 40 candidate genes, and 9 of the genes (ARID1B, KAT6B, TCF4, SMARCA2, SHANK3, PDHA1, WDR45, SCN2A, SYNGAP1) harboured recurrent NCSVs with the same variant present in more than two unrelated patients with NDD. Moreover, 36 of 59 (61.02%) NCSVs are novel clinically relevant variants, including 34 unreported and 2 clinically conflicting interpretations or of uncertain significance NCSVs in the ClinVar database.

INTERPRETATION

This study highlights the common pathology and clinical importance of NCSVs in unsolved patients with NDD.

FUNDING

The present study was funded by grants from the National Natural Science Foundation of China, China Postdoctoral Science Foundation, the Hunan Youth Science and Technology Innovation Talent Project, the Provincial Natural Science Foundation of Hunan, The Scientific Research Program of FuRong laboratory, and the Natural Science Project of the University of Anhui Province.

Benchmarking splice variant prediction algorithms using massively parallel splicing assays

BACKGROUND

Variants that disrupt mRNA splicing account for a sizable fraction of the pathogenic burden in many genetic disorders, but identifying splice-disruptive variants (SDVs) beyond the essential splice site dinucleotides remains difficult. Computational predictors are often discordant, compounding the challenge of variant interpretation. Because they are primarily validated using clinical variant sets heavily biased to known canonical splice site mutations, it remains unclear how well their performance generalizes.

RESULTS

We benchmark eight widely used splicing effect prediction algorithms, leveraging massively parallel splicing assays (MPSAs) as a source of experimentally determined ground-truth. MPSAs simultaneously assay many variants to nominate candidate SDVs. We compare experimentally measured splicing outcomes with bioinformatic predictions for 3,616 variants in five genes. Algorithms' concordance with MPSA measurements, and with each other, is lower for exonic than intronic variants, underscoring the difficulty of identifying missense or synonymous SDVs. Deep learning-based predictors trained on gene model annotations achieve the best overall performance at distinguishing disruptive and neutral variants, and controlling for overall call rate genome-wide, SpliceAI and Pangolin have superior sensitivity. Finally, our results highlight two practical considerations when scoring variants genome-wide: finding an optimal score cutoff, and the substantial variability introduced by differences in gene model annotation, and we suggest strategies for optimal splice effect prediction in the face of these issues.

CONCLUSION

SpliceAI and Pangolin show the best overall performance among predictors tested, however, improvements in splice effect prediction are still needed especially within exons.

Genome-wide prediction of pathogenic gain- and loss-of-function variants from ensemble learning of a diverse feature set

Gain-of-function (GOF) variants give rise to increased/novel protein functions whereas loss-of-function (LOF) variants lead to diminished protein function. Experimental approaches for identifying GOF and LOF are generally slow and costly, whilst available computational methods have not been optimized to discriminate between GOF and LOF variants. We have developed LoGoFunc, a machine learning method for predicting pathogenic GOF, pathogenic LOF, and neutral genetic variants, trained on a broad range of gene-, protein-, and variant-level features describing diverse biological characteristics. LoGoFunc outperforms other tools trained solely to predict pathogenicity for identifying pathogenic GOF and LOF variants and is available at https://itanlab.shinyapps.io/goflof/ .

Accurate long-read sequencing identified GBA1 as major risk factor in the Luxembourgish Parkinson's study

Heterozygous variants in the glucocerebrosidase GBA1 gene are an increasingly recognized risk factor for Parkinson's disease (PD). Due to the GBAP1 pseudogene, which shares 96% sequence homology with the GBA1 coding region, accurate variant calling by array-based or short-read sequencing methods remains a major challenge in understanding the genetic landscape of GBA1-associated PD. We analyzed 660 patients with PD, 100 patients with Parkinsonism and 808 healthy controls from the Luxembourg Parkinson's study, sequenced using amplicon-based long-read DNA sequencing technology. We found that 12.1% (77/637) of PD patients carried GBA1 variants, with 10.5% (67/637) of them carrying known pathogenic variants (including severe, mild, risk variants). In comparison, 5% (34/675) of the healthy controls carried GBA1 variants, and among them, 4.3% (29/675) were identified as pathogenic variant carriers. We found four GBA1 variants in patients with atypical parkinsonism. Pathogenic GBA1 variants were 2.6-fold more frequently observed in PD patients compared to controls (OR = 2.6; CI = [1.6,4.1]). Three novel variants of unknown significance (VUS) were identified. Using a structure-based approach, we defined a potential risk prediction method for VUS. This study describes the full landscape of GBA1-related parkinsonism in Luxembourg, showing a high prevalence of GBA1 variants as the major genetic risk for PD. Although the long-read DNA sequencing technique used in our study may be limited in its effectiveness to detect potential structural variants, our approach provides an important advancement for highly accurate GBA1 variant calling, which is essential for providing access to emerging causative therapies for GBA1 carriers.

The hospital Israelita Albert Einstein standards for constitutional sequence variants classification: version 2023

BACKGROUND

Next-generation sequencing has had a significant impact on genetic disease diagnosis, but the interpretation of the vast amount of genomic data it generates can be challenging. To address this, the American College of Medical Genetics and Genomics and the Association for Molecular Pathology have established guidelines for standardized variant interpretation. In this manuscript, we present the updated Hospital Israelita Albert Einstein Standards for Constitutional Sequence Variants Classification, incorporating modifications from leading genetics societies and the ClinGen initiative.

RESULTS

First, we standardized the scientific publications, documents, and other reliable sources for this document to ensure an evidence-based approach. Next, we defined the databases that would provide variant information for the classification process, established the terminology for molecular findings, set standards for disease-gene associations, and determined the nomenclature for classification criteria. Subsequently, we defined the general rules for variant classification and the Bayesian statistical reasoning principles to enhance this process. We also defined bioinformatics standards for automated classification. Our workgroup adhered to gene-specific rules and workflows curated by the ClinGen Variant Curation Expert Panels whenever available. Additionally, a distinct set of specifications for criteria modulation was created for cancer genes, recognizing their unique characteristics.

CONCLUSIONS

The development of an internal consensus and standards for constitutional sequence variant classification, specifically adapted to the Brazilian population, further contributes to the continuous refinement of variant classification practices. The aim of these efforts from the workgroup is to enhance the reliability and uniformity of variant classification.

An autosomal-dominant childhood-onset disorder associated with pathogenic variants in VCP

Valosin-containing protein (VCP) is an AAA+ ATPase that plays critical roles in multiple ubiquitin-dependent cellular processes. Dominant pathogenic variants in VCP are associated with adult-onset multisystem proteinopathy (MSP), which manifests as myopathy, bone disease, dementia, and/or motor neuron disease. Through GeneMatcher, we identified 13 unrelated individuals who harbor heterozygous VCP variants (12 de novo and 1 inherited) associated with a childhood-onset disorder characterized by developmental delay, intellectual disability, hypotonia, and macrocephaly. Trio exome sequencing or a multigene panel identified nine missense variants, two in-frame deletions, one frameshift, and one splicing variant. We performed in vitro functional studies and in silico modeling to investigate the impact of these variants on protein function. In contrast to MSP variants, most missense variants had decreased ATPase activity, and one caused hyperactivation. Other variants were predicted to cause haploinsufficiency, suggesting a loss-of-function mechanism. This cohort expands the spectrum of VCP-related disease to include neurodevelopmental disease presenting in childhood.

Whole exome sequencing reveal 83 novel Mendelian disorders carrier P/LP variants in Chinese adult patients

Carrier screening can identify people at risk of conceiving pregnancies affected with inherited genetic disorders or who have a genetic disorder with late or variable onset. Carrier screening based on whole exome sequencing (WES) data can offer more comprehensive assessment than on-target carrier screening tests. A total of 224 Chinese adult patients WES data was analyzed, except positive variants associated with the patients' major complaint, 378 pathogenic (P) and "likely pathogenic" (LP) variants from 175 adult patients were identified. Whole exome-wide frequency of carriers for Mendelian disorders in Chinese adult patients was about 78.13% in this study, which was lower than the previously reported carrier frequency in healthy population. Contrary to expectations, the number of P or LP variants did not increase with larger chromosome size or decrease with smaller chromosome size. Totally 83 novel P or LP variants were identified which could further expand the carrier variants spectrum of the Chinese population. GJB2: NM_004004.6:c.299_300delAT:p.His100fs*14 and C6:NM_000065.4:c.654T>A:p.Cys218* were found in two or more patients, which might be two underestimated carrier variants in Chinese population. We also found 9 late-onset or atypical symptoms autosomal/X-linked dominant Mendelian disorders causative genes, which were easily overlooked during pathogenicity analysis. These results can provide a strong basis for preventing and avoiding the prevalence rates of birth defects and reducing social and family burdens. By comparing with three different expanded carrier screening gene panels, we further confirmed carrier screening based on WES could offer more comprehensive assessment and WES was applicable for carrier screening.

Dominant negative variants in IKZF2 cause ICHAD syndrome, a new disorder characterised by immunodysregulation, craniofacial anomalies, hearing impairment, athelia and developmental delay

BACKGROUND

Helios (encoded by IKZF2), a member of the Ikaros family of transcription factors, is a zinc finger protein involved in embryogenesis and immune function. Although predominantly recognised for its role in the development and function of T lymphocytes, particularly the CD4+ regulatory T cells (Tregs), the expression and function of Helios extends beyond the immune system. During embryogenesis, Helios is expressed in a wide range of tissues, making genetic variants that disrupt the function of Helios strong candidates for causing widespread immune-related and developmental abnormalities in humans.

METHODS

We performed detailed phenotypic, genomic and functional investigations on two unrelated individuals with a phenotype of immune dysregulation combined with syndromic features including craniofacial differences, sensorineural hearing loss and congenital abnormalities.

RESULTS

Genome sequencing revealed de novo heterozygous variants that alter the critical DNA-binding zinc fingers (ZFs) of Helios. Proband 1 had a tandem duplication of ZFs 2 and 3 in the DNA-binding domain of Helios (p.Gly136_Ser191dup) and Proband 2 had a missense variant impacting one of the key residues for specific base recognition and DNA interaction in ZF2 of Helios (p.Gly153Arg). Functional studies confirmed that both these variant proteins are expressed and that they interfere with the ability of the wild-type Helios protein to perform its canonical function-repressing IL2 transcription activity-in a dominant negative manner.

CONCLUSION

This study is the first to describe dominant negative IKZF2 variants. These variants cause a novel genetic syndrome characterised by immunodysregulation, craniofacial anomalies, hearing impairment, athelia and developmental delay.

The burden of splice-disrupting variants in inherited heart disease and unexplained sudden cardiac death

There is an incomplete understanding of the burden of splice-disrupting variants in definitively associated inherited heart disease genes and whether these genes can amplify from blood RNA to support functional confirmation of splicing outcomes. We performed burden testing of rare splice-disrupting variants in people with inherited heart disease and sudden unexplained death compared to 125,748 population controls. ClinGen definitively disease-associated inherited heart disease genes were amplified using RNA extracted from fresh blood, derived cardiomyocytes, and myectomy tissue. Variants were functionally assessed and classified for pathogenicity. We found 88 in silico-predicted splice-disrupting variants in 128 out of 1242 (10.3%) unrelated participants. There was an excess burden of splice-disrupting variants in PKP2 (5.9%), FLNC (2.7%), TTN (2.8%), MYBPC3 (8.2%) and MYH7 (1.3%), in distinct cardiomyopathy subtypes, and KCNQ1 (3.6%) in long QT syndrome. Blood RNA supported the amplification of 21 out of 31 definitive disease-associated inherited heart disease genes. Our functional studies confirmed altered splicing in six variants. Eleven variants of uncertain significance were reclassified as likely pathogenic based on functional studies and six were used for cascade genetic testing in 12 family members. Our study highlights that splice-disrupting variants are a significant cause of inherited heart disease, and that analysis of blood RNA confirms splicing outcomes and supports variant pathogenicity classification.

A Splicing Variant in RDH8 Is Associated with Autosomal Recessive Stargardt Macular Dystrophy

Stargardt macular dystrophy is a genetic disorder, but in many cases, the causative gene remains unrevealed. Through a combined approach (whole-exome sequencing and phenotype/family-driven filtering algorithm) and a multilevel validation (international database searching, prediction scores calculation, splicing analysis assay, segregation analyses), a biallelic mutation in the RDH8 gene was identified to be responsible for Stargardt macular dystrophy in a consanguineous Italian family. This paper is a report on the first family in which a biallelic deleterious mutation in RDH8 is detected. The disease phenotype is consistent with the expected phenotype hypothesized in previous studies on murine models. The application of the combined approach to genetic data and the multilevel validation allowed the identification of a splicing mutation in a gene that has never been reported before in human disorders.

Puzzling Out the Genetic Architecture of Endometriosis: Whole-Exome Sequencing and Novel Candidate Gene Identification in a Deeply Clinically Characterised Cohort

Endometriosis (EM) is a common multifactorial gynaecological disorder. Although Genome-Wide Association Studies have largely been employed, the current knowledge of the genetic mechanisms underlying EM is far from complete, and other approaches are needed. To this purpose, whole-exome sequencing (WES) was performed on a deeply characterised cohort of 80 EM patients aimed at the identification of rare and damaging variants within 46 EM-associated genes and novel candidates. WES analysis detected 63 rare, predicted, and damaging heterozygous variants within 24 genes in 63% of the EM patients. In particular, (1) a total of 43% of patients carried variants within 13 recurrent genes (FCRL3, LAMA5, SYNE1, SYNE2, GREB1, MAP3K4, C3, MMP3, MMP9, TYK2, VEGFA, VEZT, RHOJ); (2) a total of 8.8% carried private variants within eight genes (KAZN, IL18, WT1, CYP19A1, IL1A, IL2RB, LILRB2, ZNF366); (3) a total of 24% carried variants within three novel candidates (ABCA13, NEB, CSMD1). Finally, to deepen the polygenic architecture of EM, a comprehensive evaluation of the analysed genes was performed, revealing a higher burden (p < 0.05) of genes harbouring rare and damaging variants in the EM patients than in the controls. These results highlight new insights into EM genetics, allowing for the definition of novel genotype-phenotype correlations, thereby contributing, in a long-term perspective, to the development of personalised care for EM patients.

Evidence for a Founder Effect of SDHB Exon 1 Deletion in Brazilian Patients With Paraganglioma

CONTEXT

Limited information is available concerning the genetic spectrum of pheochromocytoma and paraganglioma (PPGL) patients in South America. Germline SDHB large deletions are very rare worldwide, but most of the individuals harboring the SDHB exon 1 deletion originated from the Iberian Peninsula.

OBJECTIVE

Our aim was to investigate the spectrum of SDHB genetic defects in a large cohort of Brazilian patients with PPGLs.

METHODS

Genetic investigation of 155 index PPGL patients was performed by Sanger DNA sequencing, multiplex ligation-dependent probe amplification, and/or target next-generation sequencing panel. Common ancestrality was investigated by microsatellite genotyping with haplotype reconstruction, and analysis of deletion breakpoint.

RESULTS

Among 155 index patients, heterozygous germline SDHB pathogenic or likely pathogenic variants were identified in 22 cases (14.2%). The heterozygous SDHB exon 1 complete deletion was the most frequent genetic defect in SDHB, identified in 8 out of 22 (36%) of patients. Haplotype analysis of 5 SDHB flanking microsatellite markers demonstrated a significant difference in haplotype frequencies in a case-control permutation test (P = 0.03). More precisely, 3 closer/informative microsatellites were shared by 6 out of 8 apparently unrelated cases (75%) (SDHB-GATA29A05-D1S2826-D1S2644 | SDHB-186-130-213), which was observed in only 1 chromosome (1/42) without SDHB exon 1 deletion (X2 = 29.43; P < 0.001). Moreover, all cases with SDHB exon 1 deletion had the same gene breakpoint pattern of a 15 678 bp deletion previously described in the Iberian Peninsula, indicating a common origin.

CONCLUSION

The germline heterozygous SDHB exon 1 deletion was the most frequent genetic defect in the Brazilian PPGL cohort. Our findings demonstrated a founder effect for the SDHB exon 1 deletion in Brazilian patients with paragangliomas.

Gene-specific machine learning for pathogenicity prediction of rare BRCA1 and BRCA2 missense variants

Machine learning-based pathogenicity prediction helps interpret rare missense variants of BRCA1 and BRCA2, which are associated with hereditary cancers. Recent studies have shown that classifiers trained using variants of a specific gene or a set of genes related to a particular disease perform better than those trained using all variants, due to their higher specificity, despite the smaller training dataset size. In this study, we further investigated the advantages of "gene-specific" machine learning compared to "disease-specific" machine learning. We used 1068 rare (gnomAD minor allele frequency (MAF) < 0.005) missense variants of 28 genes associated with hereditary cancers for our investigation. Popular machine learning classifiers were employed: regularized logistic regression, extreme gradient boosting, random forests, support vector machines, and deep neural networks. As features, we used MAFs from multiple populations, functional prediction and conservation scores, and positions of variants. The disease-specific training dataset included the gene-specific training dataset and was > 7 × larger. However, we observed that gene-specific training variants were sufficient to produce the optimal pathogenicity predictor if a suitable machine learning classifier was employed. Therefore, we recommend gene-specific over disease-specific machine learning as an efficient and effective method for predicting the pathogenicity of rare BRCA1 and BRCA2 missense variants.

Intronic position +9 and -9 are potentially splicing sites boundary from intronic variants analysis of whole exome sequencing data

Whole exome sequencing (WES) can also detect some intronic variants, which may affect splicing and gene expression, but how to use these intronic variants, and the characteristics about them has not been reported. This study aims to reveal the characteristics of intronic variant in WES data, to further improve the clinical diagnostic value of WES. A total of 269 WES data was analyzed, 688,778 raw variants were called, among these 367,469 intronic variants were in intronic regions flanking exons which was upstream/downstream region of the exon (default is 200 bps). Contrary to expectation, the number of intronic variants with quality control (QC) passed was the lowest at the +2 and -2 positions but not at the +1 and -1 positions. The plausible explanation was that the former had the worst effect on trans-splicing, whereas the latter did not completely abolish splicing. And surprisingly, the number of intronic variants that passed QC was the highest at the +9 and -9 positions, indicating a potential splicing site boundary. The proportion of variants which could not pass QC filtering (false variants) in the intronic regions flanking exons generally accord with "S"-shaped curve. At +5 and -5 positions, the number of variants predicted damaging by software was most. This was also the position at which many pathogenic variants had been reported in recent years. Our study revealed the characteristics of intronic variant in WES data for the first time, we found the +9 and -9 positions might be a potentially splicing sites boundary and +5 and -5 positions were potentially important sites affecting splicing or gene expression, the +2 and -2 positions seem more important splicing site than +1 and -1 positions, and we found variants in intronic regions flanking exons over ± 50 bps may be unreliable. This result can help researchers find more useful variants and demonstrate that WES data is valuable for intronic variants analysis.

PacBio single-molecule long-read sequencing provides new insights into the complexity of full-length transcripts in oriental river prawn, macrobrachium nipponense

BACKGROUND

Oriental river prawn (Macrobrachium nipponense) is one of the most dominant species in shrimp farming in China, which is a rich source of protein and contributes to a significant impact on the quality of human life. Thus, more complete and accurate annotation of gene models are important for the breeding research of oriental river prawn.

RESULTS

A full-length transcriptome of oriental river prawn muscle was obtained using the PacBio Sequel platform. Then, 37.99 Gb of subreads were sequenced, including 584,498 circular consensus sequences, among which 512,216 were full length non-chimeric sequences. After Illumina-based correction of long PacBio reads, 6,599 error-corrected isoforms were identified. Transcriptome structural analysis revealed 2,263 and 2,555 alternative splicing (AS) events and alternative polyadenylation (APA) sites, respectively. In total, 620 novel genes (NGs), 197 putative transcription factors (TFs), and 291 novel long non-coding RNAs (lncRNAs) were identified.

CONCLUSIONS

In summary, this study offers novel insights into the transcriptome complexity and diversity of this prawn species, and provides valuable information for understanding the genomic structure and improving the draft genome annotation of oriental river prawn.

Case report: Autosomal recessive type 3 Stickler syndrome caused by compound heterozygous mutations in COL11A2

Background: Stickler syndrome (SS) is a group of hereditary collagenopathies caused by a variety of collagen and non-collagen genes. Affected patients have characteristic manifestations involving ophthalmic, articular, craniofacial and auditory disorders. SS is classified into several subtypes according to clinical and molecular features. Type 3 SS is an ultra-rare disease, known as non-ocular SS or otospondylomegaepiphyseal dysplasia (OSMED) with only a few pathogenic COL11A2 variants reported to date. Case presentation: A 29-year-old Chinese male was referred to our hospital for hearing loss and multiple joint pain. He presented a phenotype highly suggestive of OSMED, including progressive sensorineural deafness, spondyloepiphyseal dysplasia with large epiphyses, platyspondyly, degenerative osteoarthritis, and sunken nasal bridge. We detected compound heterozygous mutations in COL11A2, both of which were predicted to be splicing mutations. One is synonymous mutation c.3774C>T (p.Gly1258Gly) supposed to be a splice site mutation, the other is a novel intron mutation c.4750 + 5 G>A, which is a highly conservative site across several species. We also present a review of the current known pathogenic mutation spectrum of COL11A2 in patients with type 3 SS. Conclusion: Both synonymous extonic and intronic variants are easily overlooked by whole-exome sequencing. For patients with clinical manifestations suspected of SS syndrome, next-generation whole-genome sequencing is necessary for precision diagnosis and genetic counseling.

In silico methods for predicting functional synonymous variants

Single nucleotide variants (SNVs) contribute to human genomic diversity. Synonymous SNVs are previously considered to be "silent," but mounting evidence has revealed that these variants can cause RNA and protein changes and are implicated in over 85 human diseases and cancers. Recent improvements in computational platforms have led to the development of numerous machine-learning tools, which can be used to advance synonymous SNV research. In this review, we discuss tools that should be used to investigate synonymous variants. We provide supportive examples from seminal studies that demonstrate how these tools have driven new discoveries of functional synonymous SNVs.

Introme accurately predicts the impact of coding and noncoding variants on gene splicing, with clinical applications

Predicting the impact of coding and noncoding variants on splicing is challenging, particularly in non-canonical splice sites, leading to missed diagnoses in patients. Existing splice prediction tools are complementary but knowing which to use for each splicing context remains difficult. Here, we describe Introme, which uses machine learning to integrate predictions from several splice detection tools, additional splicing rules, and gene architecture features to comprehensively evaluate the likelihood of a variant impacting splicing. Through extensive benchmarking across 21,000 splice-altering variants, Introme outperformed all tools (auPRC: 0.98) for the detection of clinically significant splice variants. Introme is available at https://github.com/CCICB/introme .