Functional Analyses of a Novel Splice Variant in the CHD7 Gene, Found by Next Generation Sequencing, Confirm Its Pathogenicity in a Spanish Patient and Diagnose Him with CHARGE Syndrome

Splicing Analysis of 16 PALB2 ClinVar Variants by Minigene Assays: Identification of Six Likely Pathogenic Variants

PALB2 loss-of-function variants are associated with significant increased risk of breast cancer as well as other types of tumors. Likewise, splicing disruptions are a common mechanism of disease susceptibility. Indeed, we previously showed, by minigene assays, that 35 out of 42 PALB2 variants impaired splicing. Taking advantage of one of these constructs (mgPALB2_ex1-3), we proceeded to analyze other variants at exons 1 to 3 reported at the ClinVar database. Thirty-one variants were bioinformatically analyzed with MaxEntScan and SpliceAI. Then, 16 variants were selected for subsequent RNA assays. We identified a total of 12 spliceogenic variants, 11 of which did not produce any trace of the expected minigene full-length transcript. Interestingly, variant c.49-1G > A mimicked previous outcomes in patient RNA (transcript ∆(E2p6)), supporting the reproducibility of the minigene approach. A total of eight variant-induced transcripts were characterized, three of which (∆(E1q17), ∆(E3p11), and ∆(E3)) were predicted to introduce a premature termination codon and to undergo nonsense-mediated decay, and five (▼(E1q9), ∆(E2p6), ∆(E2), ▼(E3q48)-a, and ▼(E3q48)-b) maintained the reading frame. According to an ACMG/AMP (American College of Medical Genetics and Genomics/Association for Molecular Pathology)-based classification scheme, which integrates mgPALB2 data, six PALB2 variants were classified as pathogenic/likely pathogenic, five as VUS, and five as likely benign. Furthermore, five ±1,2 variants were catalogued as VUS because they produced significant proportions of in-frame transcripts of unknown impact on protein function.

Opitz syndrome: improving clinical interpretation of intronic variants in MID1 gene

BACKGROUND

Loss-of-function variants in MID1 are the most common cause of Opitz G/BBB syndrome (OS). The interpretation of intronic variants affecting the splicing is a rising issue in OS.

METHODS

Exon sequencing of a 2-year-old boy with OS showed that he was a carrier of the de novo c.1286-10G>T variant in MID1. In silico predictions and minigene assays explored the effect of the variant on splicing. The minigene approach was also applied to two previously identified MID1 c.864+1G>T and c.1285+1G>T variants.

RESULTS

Minigene assay demonstrated that the c.1286-10G>T variant generated the inclusion of eight nucleotides that predicted generation of a frameshift. The c.864+1G>T and c.1285+1G>T variants resulted in an in-frame deletion predicted to generate a shorter MID1 protein. In hemizygous males, this allowed reclassification of all the identified variants from "of unknown significance" to "likely pathogenic."

CONCLUSIONS

Minigene assay supports functional effects from MID1 intronic variants. This paves the way to the introduction of similar second-tier investigations in the molecular diagnostics workflow of OS.

IMPACT

Causative intronic variants in MID1 are rarely investigated in Opitz syndrome. MID1 is not expressed in blood and mRNA studies are hardly accessible in routine diagnostics. Minigene assay is an alternative for assessing the effect of intronic variants on splicing. This is the first study characterizing the molecular consequences of three MID1 variants for diagnostic purposes and demonstrating the efficacy of minigene assays in supporting their clinical interpretation. Review of the criteria according to the American College of Medical Genetics reassessed all variants as likely pathogenic.

RNA Analysis and Clinical Characterization of a Novel Splice Variant in the NSD1 Gene Causing Familial Sotos Syndrome

BACKGROUND

Sotos syndrome is an autosomal dominant disorder characterized by overgrowth, macrocephaly, distinctive facial features and learning disabilities. Haploinsufficiency of the nuclear receptor SET domain-containing protein 1 (NSD1) gene located on chromosome 5q35 is the major cause of the syndrome. This syndrome shares characteristics with other overgrowth syndromes, which can complicate the differential diagnosis.

METHODS

Genomic DNA was extracted from peripheral blood samples of members of the same family and targeted exome analysis was performed. In silico study of the variant found by next-generation sequencing was used to predict disruption/creation of splice sites and the identification of potential cryptic splice sites. RNA was extracted from peripheral blood samples of patients and functional analyses were performed to confirm the pathogenicity.

RESULTS

We found a novel c.6463 + 5G>A heterozygous NSD1 gene pathogenic variant in a son and his father. Molecular analyses revealed that part of the intron 22 of NSD1 is retained due to the destruction of the splicing donor site, causing the appearance of a premature stop codon in the NSD1 protein.

CONCLUSIONS

Our findings underline the importance of performing RNA functional assays in order to determine the clinical significance of intronic variants, and contribute to the genetic counseling and clinical management of patients and their relatives. Our work also highlights the relevance of using in silico prediction tools to detect a potential alteration in the splicing process.

Hybrid Minigene Assay: An Efficient Tool to Characterize mRNA Splicing Profiles of NF1 Variants

Neurofibromatosis type 1 (NF1) is caused by heterozygous loss of function mutations in the NF1 gene. Although patients are diagnosed according to clinical criteria and few genotype-phenotype correlations are known, molecular analysis remains important. NF1 displays allelic heterogeneity, with a high proportion of variants affecting splicing, including deep intronic alleles and changes outside the canonical splice sites, making validation problematic. Next Generation Sequencing (NGS) technologies integrated with multiplex ligation-dependent probe amplification (MLPA) have largely overcome RNA-based techniques but do not detect splicing defects. A rapid minigene-based system was set up to test the effects of NF1 variants on splicing. We investigated 29 intronic and exonic NF1 variants identified in patients during the diagnostic process. The minigene assay showed the coexistence of multiple mechanisms of splicing alterations for seven variants. A leaky effect on splicing was documented in one de novo substitution detected in a sporadic patient with a specific phenotype without neurofibromas. Our splicing assay proved to be a reliable and fast method to validate novel NF1 variants potentially affecting splicing and to detect hypomorphic effects that might have phenotypic consequences, avoiding the requirement of patient's RNA.

MDAKRLS: Predicting human microbe-disease association based on Kronecker regularized least squares and similarities

BACKGROUND

Microbes are closely related to human health and diseases. Identification of disease-related microbes is of great significance for revealing the pathological mechanism of human diseases and understanding the interaction mechanisms between microbes and humans, which is also useful for the prevention, diagnosis and treatment of human diseases. Considering the known disease-related microbes are still insufficient, it is necessary to develop effective computational methods and reduce the time and cost of biological experiments.

METHODS

In this work, we developed a novel computational method called MDAKRLS to discover potential microbe-disease associations (MDAs) based on the Kronecker regularized least squares. Specifically, we introduced the Hamming interaction profile similarity to measure the similarities of microbes and diseases besides Gaussian interaction profile kernel similarity. In addition, we introduced the Kronecker product to construct two kinds of Kronecker similarities between microbe-disease pairs. Then, we designed the Kronecker regularized least squares with different Kronecker similarities to obtain prediction scores, respectively, and calculated the final prediction scores by integrating the contributions of different similarities.

RESULTS

The AUCs value of global leave-one-out cross-validation and 5-fold cross-validation achieved by MDAKRLS were 0.9327 and 0.9023 ± 0.0015, which were significantly higher than five state-of-the-art methods used for comparison. Comparison results demonstrate that MDAKRLS has faster computing speed under two kinds of frameworks. In addition, case studies of inflammatory bowel disease (IBD) and asthma further showed 19 (IBD), 19 (asthma) of the top 20 prediction disease-related microbes could be verified by previously published biological or medical literature.

CONCLUSIONS

All the evaluation results adequately demonstrated that MDAKRLS has an effective and reliable prediction performance. It may be a useful tool to seek disease-related new microbes and help biomedical researchers to carry out follow-up studies.

Comprehensive Functional Characterization and Clinical Interpretation of 20 Splice-Site Variants of the RAD51C Gene

Hereditary breast and/or ovarian cancer is a highly heterogeneous disease with more than 10 known disease-associated genes. In the framework of the BRIDGES project (Breast Cancer Risk after Diagnostic Gene Sequencing), the RAD51C gene has been sequenced in 60,466 breast cancer patients and 53,461 controls. We aimed at functionally characterizing all the identified genetic variants that are predicted to disrupt the splicing process. Forty RAD51C variants of the intron-exon boundaries were bioinformatically analyzed, 20 of which were selected for splicing functional assays. To test them, a splicing reporter minigene with exons 2 to 8 was designed and constructed. This minigene generated a full-length transcript of the expected size (1062 nucleotides), sequence, and structure (Vector exon V1- RAD51C exons_2-8- Vector exon V2). The 20 candidate variants were genetically engineered into the wild type minigene and functionally assayed in MCF-7 cells. Nineteen variants (95%) impaired splicing, while 18 of them produced severe splicing anomalies. At least 35 transcripts were generated by the mutant minigenes: 16 protein-truncating, 6 in-frame, and 13 minor uncharacterized isoforms. According to ACMG/AMP-based standards, 15 variants could be classified as pathogenic or likely pathogenic variants: c.404G > A, c.405-6T > A, c.571 + 4A > G, c.571 + 5G > A, c.572-1G > T, c.705G > T, c.706-2A > C, c.706-2A > G, c.837 + 2T > C, c.905-3C > G, c.905-2A > C, c.905-2_905-1del, c.965 + 5G > A, c.1026 + 5_1026 + 7del, and c.1026 + 5G > T.

Variable phenotypes and outcomes associated with the MMACHC c.609G>A homologous mutation: long term follow-up in a large cohort of cases

Background

Cobalamin C deficiency (cblC) caused by the MMACHC mutations is the most common type of the disorders of intracellular cobalamin metabolism. While the c.609G > A mutation is most frequent in Chinese cblC patients, its correlation with phenotype has not been delineated. Here we aim to investigate the factors affecting variable phenotypes and outcomes associated with the MMACHC c.609G > A homologous mutation in 149 Chinese cases to have implications for treatment and prevention.

Methods

We assessed 149 cblC patients caused by MMACHC c.609G > A homozygous mutation. The clinical manifestations, complications, treatment, and outcomes were evaluated; 120 patients were followed-up till December 2019.

Results

Two patients (1.3%) were prenatally diagnosed, treated after birth and consequently showed normal development. In 15 patients (10.1%) detected by newborn screening, 10 were treated at the age of 2 weeks and showed normal development, while the other 5 were treated after onset and showed neurologic disorders. All 132 clinically diagnosed patients (88.6%) developed symptoms at age from few minutes after birth to 72 months. Among them, 101 (76.5%) had early-onset (before the age of 12 months) and 31 (23.5%) had late-onset (after the age of 12 months). Totally 5 patients died and 24 were lost to follow-up. Of the 132 clinical diagnosed patients, 92 (69.7%) presented with developmental delay, 65 (49.2%) had seizures, 37 (28.0%) had anemia, 24 (18.2%) had feeding difficulty, 23 (17.4%) had ocular problems, and 22 (16.7%) had hydrocephalus. Compared with the non-developmental delay group, the onset age, the age at treatment initiation and the time from onset to treatment initiation were later in the developmental delay group. Seizure group showed significantly higher urinary methylmalonic acid concentration. During long-term follow-up, plasma total homocysteine (tHcy) levels were significantly higher in patients in the uncontrolled group than those in the seizure-free group.

Conclusions

Most cblC patients caused by MMACHC c.609G > A homozygous mutation showed early-onset. The clinically diagnosed patients usually showed the presence of irreversible brain disorders. Patients treated from the pre-symptomatic stage showed favorable outcomes. Therefore, newborn screening, prenatal diagnosis and early treatment are crucial and the c.609G > A mutant allele should be listed in the pre-pregnancy carrier screening panel in China.

A Novel Cosegregating DCTN1 Splice Site Variant in a Family with Bipolar Disorder May Hold the Key to Understanding the Etiology

A novel cosegregating splice site variant in the Dynactin-1 (DCTN1) gene was discovered by Next Generation Sequencing (NGS) in a family with a history of bipolar disorder (BD) and major depressive diagnosis (MDD). Psychiatric illness in this family follows an autosomal dominant pattern. DCTN1 codes for the largest dynactin subunit, namely p150^Glued, which plays an essential role in retrograde axonal transport and in neuronal autophagy. A GT→TT transversion in the DCTN1 gene, uncovered in the present work, is predicted to disrupt the invariant canonical splice donor site IVS22 + 1G > T and result in intron retention and a premature termination codon (PTC). Thus, this splice site variant is predicted to trigger RNA nonsense-mediated decay (NMD) and/or result in a C-terminal truncated p150^Glued protein (ct-p150^Glued), thereby negatively impacting retrograde axonal transport and neuronal autophagy. BD prophylactic medications, and most antipsychotics and antidepressants, are known to enhance neuronal autophagy. This variant is analogous to the dominant-negative GLUED Gl¹ mutation in Drosophila, which is responsible for a neurodegenerative phenotype. The newly identified variant may reflect an autosomal dominant cause of psychiatric pathology in this affected family. Factors that affect alternative splicing of the DCTN1 gene, leading to NMD and/or ct-p150^Glued, may be of fundamental importance in contributing to our understanding of the etiology of BD as well as MDD.

UGT1A1 Variants c.864+5G>T and c.996+2_996+5del of a Crigler-Najjar Patient Induce Aberrant Splicing in Minigene Assays

A large fraction of DNA variants impairs pre-mRNA splicing in human hereditary disorders. Crigler-Najjar syndrome (CNS) is characterized by a severe unconjugated hyperbilirubinemia caused by variants in the UGT1A1 gene. We previously reported one CNS-type II patient with two splice-site variants in trans (c.864+5G>T and c.996+2_996+5del). According to MaxEntScan, both disrupt their corresponding donor sites (c.864+5G>T: 6.99 → 2.28; c.996+2_996+5del: 5.96 → -11.02), so they were selected for subsequent functional tests. Given the unavailability of patient RNA, we constructed an UGT1A1 splicing-reporter minigene with exons 1-4 to characterize the underlying splicing anomaly. The variant c.996+2_996+5del generated two aberrant transcripts, Δ(E2) (exon 2 skipping/64%) and ▼(E2q135) (intron retention of 135-nt/36%), which lead to the loss of 18 conserved amino-acids and the gain of 45 new ones of a critical functional domain, respectively. The c.864+5G>T variant mainly produced the aberrant transcript Δ(E1q141) (141-nt deletion/70.4%) and the full-length isoform (29.6%). Δ(E1q141) would provoke the loss of 47 amino-acids of the N-terminal domain that encodes for substrate specificity. Thus, the three anomalous transcripts are likely to inactivate UGT1A1. Moreover, this patient is also homozygous for the promoter variant A(TA)7TAA that decreases UGT1A1 expression by 70%, so the full-length transcript produced by c.864+5G>T would be even more reduced (<9%), thus supporting the diagnosis of CNS-type II. Therefore, minigenes represent valuable tools for the functional and clinical classifications of genetic variants.

A novel CHD7 variant disrupting acceptor splice site in a patient with mild features of CHARGE syndrome: a case report

Background

CHARGE syndrome (MIM# 214800)—which is characterised by a number of congenital anomalies including coloboma, ear anomalies, deafness, facial anomalies, heart defects, atresia choanae, genital hypoplasia, growth retardation, and developmental delay—is caused by a heterozygous variant in the CHD7 (MIM# 608892) gene located on chromosome 8q12. We report the identification of a novel c.5535-1G > A variant in CHD7 and provide the evaluation of its effect on pre-mRNA splicing.

Case presentation

In this study, we report on a female presenting features of CHARGE syndrome. A novel heterozygous CHD7 variant c.5535-1G > A located in the acceptor splice site of intron 26 was identified in the proband’s DNA sample after analysis of whole exome sequencing data. In silico predictions indicating that the variant is probably pathogenic by affecting pre-mRNA splicing were verified by genetic analysis based on reverse transcription of the patient’s RNA followed by PCR amplifications performed on synthesised cDNA and Sanger sequencing. Sanger sequencing of cDNA revealed that the c.5535-1G > A variant disrupts the original acceptor splice site and activates a cryptic splice site only one nucleotide downstream of the pathogenic variant site. This change causes the omission of the first nucleotide of exon 27, leading to a frameshift in the mRNA of the CHD7 gene. Our results suggest that the alteration induces the premature truncation of the CHD7 protein (UniProtKB: Q9P2D1), thus resulting in CHARGE syndrome.

Conclusion

Genetic analysis of novel splice site variant underlines its importance for studying the pathogenic splicing mechanism as well as for confirming a diagnosis.

Electronic supplementary material

The online version of this article (10.1186/s12881-019-0859-y) contains supplementary material, which is available to authorized users.

Minigene Splicing Assays Identify 12 Spliceogenic Variants of BRCA2 Exons 14 and 15

A relevant fraction of BRCA2 variants is associated with splicing alterations and with an increased risk of hereditary breast and ovarian cancer (HBOC). In this work, we have carried out a thorough study of variants from BRCA2 exons 14 and 15 reported at mutation databases. A total of 294 variants from exons 14 and 15 and flanking intronic sequences were analyzed with the online splicing tools NNSplice and Human Splicing Finder. Fifty-three out of these 294 variants were selected as candidate splicing variants. All variants but one, were introduced into the minigene MGBR2_ex14-20 (with exons 14–20) by site-directed mutagenesis and assayed in MCF-7 cells. Twelve of the remaining 52 variants (23.1%) impaired splicing at different degrees, yielding from 5 to 100% of aberrant transcripts. Nine variants affected the natural acceptor or donor sites of both exons and three affected putative enhancers or silencers. Fluorescent capillary electrophoresis revealed at least 10 different anomalous transcripts: (E14q5), Δ (E14p10), Δ(E14p246), Δ(E14q256), Δ(E14), Δ(E15p12), Δ(E15p13), Δ(E15p83), Δ(E15) and a 942-nt fragment of unknown structure. All transcripts, except for Δ(E14q256) and Δ(E15p12), are expected to truncate the BRCA2 protein. Nine variants induced severe splicing aberrations with more than 90% of abnormal transcripts. Thus, according to the guidelines of the American College of Medical Genetics and Genomics, eight variants should be classified as pathogenic (c.7008-2A > T, c.7008-1G > A, c.7435+1G > C, c.7436-2A > T, c.7436-2A > G, c.7617+1G > A, c.7617+1G > T, and c.7617+2T > G), one as likely pathogenic (c.7008-3C > G) and three remain as variants of uncertain clinical significance or VUS (c.7177A > G, c.7447A > G and c.7501C > T). In conclusion, functional assays by minigenes constitute a valuable strategy to primarily check the splicing impact of DNA variants and their clinical interpretation. While bioinformatics predictions of splice site variants were accurate, those of enhancer or silencer variants were poor (only 3/23 spliceogenic variants) which showed weak impacts on splicing (∼5–16% of aberrant isoforms). So, the Exonic Splicing Enhancer and Silencer (ESE and ESS, respectively) prediction algorithms require further improvement.

Mis-splicing in breast cancer: identification of pathogenic BRCA2 variants by systematic minigene assays

Splicing disruption is a common mechanism of gene inactivation associated with germline variants of susceptibility genes. To study the role of BRCA2 mis-splicing in hereditary breast/ovarian cancer (HBOC), we performed a comprehensive analysis of variants from BRCA2 exons 2-9, as well as the initial characterization of the regulatory mechanisms of such exons. A pSAD-based minigene with exons 2-9 was constructed and validated in MCF-7 cells, producing the expected transcript (1016-nt/V1-BRCA2_exons_2-9-V2). DNA variants from mutational databases were analyzed by NNSplice and Human Splicing Finder softwares. To refine ESE-variant prediction, we mapped the regulatory regions through a functional strategy whereby 26 exonic microdeletions were introduced into the minigene and tested in MCF-7 cells. Thus, we identified nine spliceogenic ESE-rich intervals where ESE-variants may be located. Combining bioinformatics and microdeletion assays, 83 variants were selected and genetically engineered in the minigene. Fifty-three changes impaired splicing: 28 variants disrupted the canonical sites, four created new ones, 10 abrogated enhancers, eight created silencers and three caused a double-effect. Notably, nine spliceogenic-ESE variants were located within ESE-containing intervals. Capillary electrophoresis and sequencing revealed more than 23 aberrant transcripts, where exon skipping was the most common event. Interestingly, variant c.67G>A triggered the usage of a noncanonical GC-donor 4-nt upstream. Thirty-six variants that induced severe anomalies (>60% aberrant transcripts) were analyzed according to the ACMG guidelines. Thus, 28 variants were classified as pathogenic, five as likely pathogenic and three as variants of uncertain significance. Interestingly, 13 VUS were reclassified as pathogenic or likely pathogenic variants. In conclusion, a large fraction of BRCA2 variants (∼64%) provoked splicing anomalies lending further support to the high prevalence of this disease-mechanism. The low accuracy of ESE-prediction algorithms may be circumvented by functional ESE-mapping that represents an optimal strategy to identify spliceogenic ESE-variants. Finally, systematic functional assays by minigenes depict a valuable tool for the initial characterization of splicing anomalies and the clinical interpretation of variants. © 2019 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

ALS and CHARGE syndrome: a clinical and genetic study

Amyotrophic Lateral Sclerosis and CHARGE syndrome are complex neurological disorders, which never occurred together in the same family and, to date, no putative correlation between them has been described on PubMed Central. Due to our aim was to evaluate the presence of different genetic variants involved in these pathologies, we reported a clinical and genetic description of two sisters affected by these two different disorders. In the CHARGE patient, molecular analysis of the CHD7 gene revealed the c.8016G >A de novo variant in exon 37. The ALS patient had been screened negative for mutations in SOD1, TARDBP, FUS/TLS, C9orf72 and KIF5A genes. Anyway, targeted next generation sequencing analysis identified known and unknown genetic variations in 39 ALS-related genes: a total of 380 variants were reported, of which 194 in the ALS patient and 186 in the CHARGE patient. To date, although the results suggest that the occurrence of the two syndromes in the same family is co-incidental rather than based on a causative genetic variant, we could hypothesize that other factors might act as modulators in the pathogenesis of these different phenotypes.

Identification of Eight Spliceogenic Variants in BRCA2 Exon 16 by Minigene Assays

Genetic testing of BRCA1 and BRCA2 identifies a large number of variants of uncertain clinical significance whose functional and clinical interpretations pose a challenge for genetic counseling. Interestingly, a relevant fraction of DNA variants can disrupt the splicing process in cancer susceptibility genes. We have tested more than 200 variants throughout 19 BRCA2 exons mostly by minigene assays, 54% of which displayed aberrant splicing, thus confirming the utility of this assay to check genetic variants in the absence of patient RNA. Our goal was to investigate BRCA2 exon 16 with a view to characterizing spliceogenic variants recorded at the mutational databases. Seventy-two different BIC and UMD variants were analyzed with NNSplice and Human Splicing Finder, 12 of which were selected because they were predicted to disrupt essential splice motifs: canonical splice sites (ss; eight variants) and exonic/intronic splicing enhancers (four variants). These 12 candidate variants were introduced into the BRCA2 minigene with seven exons (14–20) by site-directed mutagenesis and then transfected into MCF-7 cells. Seven variants (six intronic and one missense) induced complete abnormal splicing patterns: c.7618-2A>T, c.7618-2A>G, c.7618-1G>C, c.7618-1G>A, c.7805G>C, c.7805+1G>A, and c.7805+3A>C, as well as a partial anomalous outcome by c.7802A>G. They generated at least 10 different transcripts: Δ16p₄₄ (alternative 3’ss 44-nt downstream; acceptor variants), Δ16 (exon 16-skipping; donor variants), Δ16p₅₅ (alternative 3’ss 55-nt downstream), Δ16q₄ (alternative 5’ss 4-nt upstream), Δ16q₁₀₀ (alternative 5’ss 4-nt upstream), ▾16q²⁰ (alternative 5’ss 20-nt downstream), as well as minor (Δ16p₉₃ and Δ16,17p₆₉) and uncharacterized transcripts of 893 and 954 nucleotides. Isoforms Δ16p₄₄, Δ16, Δ16p₅₅, Δ16q₄, Δ16q₁₀₀, and ▾16q²⁰ introduced premature termination codons which presumably inactivate BRCA2. According to the guidelines the American College of Medical Genetics and Genomics these eight variants could be classified as pathogenic or likely pathogenic whereas the Evidence-based Network for the Interpretation of Germline Mutant Alleles rules suggested seven class 4 and one class 3 variants. In conclusion, our study highlights the relevance of splicing functional assays by hybrid minigenes for the clinical classification of genetic variations. Hence, we provide new data about spliceogenic variants of BRCA2 exon 16 that are directly correlated with breast cancer susceptibility.