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

Mutations in CHD7 have been shown to be a major cause of CHARGE syndrome, which presents many symptoms and features common to other syndromes making its diagnosis difficult. Next generation sequencing (NGS) of a panel of intellectual disability related genes was performed in an adult patient without molecular diagnosis. A splice donor variant in CHD7 (c.5665 + 1G > T) was identified. To study its potential pathogenicity, exons and flanking intronic sequences were amplified from patient DNA and cloned into the pSAD® splicing vector. HeLa cells were transfected with this construct and a wild-type minigene and functional analysis were performed. The construct with the c.5665 + 1G > T variant produced an aberrant transcript with an insert of 63 nucleotides of intron 28 creating a premature termination codon (TAG) 25 nucleotides downstream. This would lead to the insertion of 8 new amino acids and therefore a truncated 1896 amino acid protein. As a result of this, the patient was diagnosed with CHARGE syndrome. Functional analyses underline their usefulness for studying the pathogenicity of variants found by NGS and therefore its application to accurately diagnose patients.This work was funded by Jesús de Gangoiti Barrera Foundation (FJGB15/005). The EAV laboratory is funded by projects of the Spanish Ministry of Economy and Competitiveness, National Plan for R & D 2013–2016, ISCIII (FIS: PI13/01749) co-financed by FEDER from Regional Development European Funds (European Union) and the project CSI090U14 of the Regional ministry of Education (ORDER EDU/122/2014) (Castilla y León, Spain). This study made use of data generated by the UK10K Project. Funding for the UK10K Project was provided by the Wellcome Trust under award WT091310.Peer reviewe

Functional classification of DNA variants by hybrid minigenes: identification of 30 spliceogenic variants of BRCA2 exons 17 and 18

Mutation screening of the breast cancer genes BRCA1 and BRCA2 identifies a large fraction of variants of uncertain clinical significance (VUS) whose functional and clinical interpretations pose a challenge for genomic medicine. Likewise, an increasing amount of evidence indicates that genetic variants can have deleterious effects on pre-mRNA splicing. Our goal was to investigate the impact on splicing of a set of reported variants of BRCA2 exons 17 and 18 to assess their role in hereditary breast cancer and to identify critical regulatory elements that may constitute hotspots for spliceogenic variants. A splicing reporter minigene with BRCA2 exons 14 to-20 (MGBR2_ ex14-20) was constructed in the pSAD vector. Fiftytwo candidate variants were selected with splicing prediction programs, introduced in MGBR2_ ex14-20 by site-directed mutagenesis and assayed in triplicate in MCF-7 cells. Wild type MGBR2_ ex14-20 produced a stable transcript of the expected size (1,802 nucleotides) and structure (V1-[BRCA2_ exons_ 14-20]-V2). Functional mapping by microdeletions revealed essential sequences for exon recognition on the 3 ' end of exon 17 (c. 79447973) and the 5 ' end of exon 18 (c. 7979-7988, c. 7999-8013). Thirty out of the 52 selected variants induced anomalous splicing in minigene assays with > 16 different aberrant transcripts, where exon skipping was the most common event. A wide range of splicing motifs were affected including the canonical splice sites (15 variants), novel alternative sites (3 variants), the polypyrimidine tract (3 variants) and enhancers/silencers (9 variants). According to the guidelines of the American College of Medical Genetics and Genomics (ACMG), 20 variants could be classified as pathogenic (c. 7806-2A > G, c. 7806-1G > A, c. 7806-1G > T, c. 7806-1_ 7806-2dup, c. 7976+ 1G > A, c. 7977-3_ 7978del, c. 7977-2A > T, c. 7977-1G > T, c. 7977-1G > C, c. 8009C > A, c. 8331+ 1G > T and c. 8331+ 2T > C) or likely pathogenic (c. 78069T > G, c. 7976G > C, c. 7976G > A, c. 7977-7C > G, c. 7985C > G, c. 8023A > G, c. 8035G > T and c. 8331G > A), accounting for 30.8% of all pathogenic/likely pathogenic variants of exons 17-18 at the BRCA Share database. The remaining 8 variants (c. 7975A > G, c. 7977-6T > G, c. 7988A > T, c. 7992T > A, c. 8007A > G, c. 8009C > T, c. 8009C > G, and c. 8072C > T) induced partial splicing anomalies with important ratios of the full-length transcript (>= 70%), so that they remained classified as VUS. Aberrant splicing is therefore especially prevalent in BRCA2 exons 17 and 18 due to the presence of active ESEs involved in exon recognition. Splicing functional assays with minigenes are a valuable strategy for the initial characterization of the splicing outcomes and the subsequent clinical interpretation of variants of any disease-gene, although these results should be checked, whenever possible, against patient RNA

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

Germline variants of CYBA and TRPM4 predispose to familial colorectal cancer

Simple Summary Whole-genome sequencing and bioinformatics analysis on unique colorectal cancer families revealed two attractive candidate predisposition genes, CYBA and TRPM4, each with a loss-of-function variant. Supported by our functional studies, we suggest that the two gene defects mechanistically involve intestinal barrier integrity through reactive oxygen species and mucus biology, which converges in chronic bowel inflammation, a known risk factor for colorectal cancer. Abstract Familial colorectal cancer (CRC) is only partially explained by known germline predisposing genes. We performed whole-genome sequencing in 15 Polish families of many affected individuals, without mutations in known CRC predisposing genes. We focused on loss-of-function variants and functionally characterized them. We identified a frameshift variant in the CYBA gene (c.246delC) in one family and a splice site variant in the TRPM4 gene (c.25–1 G > T) in another family. While both variants were absent or extremely rare in gene variant databases, we identified four additional Polish familial CRC cases and two healthy elderly individuals with the CYBA variant (odds ratio 2.46, 95% confidence interval 0.48–12.69). Both variants led to a premature stop codon and to a truncated protein. Functional characterization of the variants showed that knockdown of CYBA or TRPM4 depressed generation of reactive oxygen species (ROS) in LS174T and HT-29 cell lines. Knockdown of TRPM4 resulted in decreased MUC2 protein production. CYBA encodes a component in the NADPH oxidase system which generates ROS and controls, e.g., bacterial colonization in the gut. Germline CYBA variants are associated with early onset inflammatory bowel disease, supported with experimental evidence on loss of intestinal mucus barrier function due to ROS deficiency. TRPM4 encodes a calcium-activated ion channel, which, in a human colonic cancer cell line, controls calcium-mediated secretion of MUC2, a major component of intestinal mucus barrier. We suggest that the gene defects in CYBA and TRPM4 mechanistically involve intestinal barrier integrity through ROS and mucus biology, which converges in chronic bowel inflammation

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

Splicing predictions, minigene analyses, and ACMG-AMP clinical classification of 42 germline PALB2 splice-site variants.

PALB2 loss-of-function variants confer high risk of developing breast cancer. Here we present a systematic functional analysis of PALB2 splice-site variants detected in approximately 113,000 women in the large-scale sequencing project Breast Cancer After Diagnostic Gene Sequencing (BRIDGES; https://bridges-research.eu/). Eighty-two PALB2 variants at the intron-exon boundaries were analyzed with MaxEntScan. Forty-two variants were selected for the subsequent splicing functional assays. For this purpose, three splicing reporter minigenes comprising exons 1-12 were constructed. The 42 potential spliceogenic variants were introduced into the minigenes by site-directed mutagenesis and assayed in MCF-7/MDA-MB-231 cells. Splicing anomalies were observed in 35 variants, 23 of which showed no traces or minimal amounts of the expected full-length transcripts of each minigene. More than 30 different variant-induced transcripts were characterized, 23 of which were predicted to truncate the PALB2 protein. The pathogenicity of all variants was interpreted according to an in-house adaptation of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology (ACMG-AMP) variant classification scheme. Up to 23 variants were classified as pathogenic/likely pathogenic. Remarkably, three ±1,2 variants (c.49-2A>T, c.108+2T>C, and c.211+1G>A) were classified as variants of unknown significance, as they produced significant amounts of either in-frame transcripts of unknown impact on the PALB2 protein function or the minigene full-length transcripts. In conclusion, we have significantly contributed to the ongoing effort of identifying spliceogenic variants in the clinically relevant PALB2 cancer susceptibility gene. Moreover, we suggest some approaches to classify the findings in accordance with the ACMG-AMP rationale. © 2021 The Authors. The Journal of Pathology published by John Wiley & Sons, Ltd on behalf of The Pathological Society of Great Britain and Ireland

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: Δ16p44 (alternative 3'ss 44-nt downstream; acceptor variants), Δ16 (exon 16-skipping; donor variants), Δ16p55 (alternative 3'ss 55-nt downstream), Δ16q4 (alternative 5'ss 4-nt upstream), Δ16q100 (alternative 5'ss 4-nt upstream), ▾16q20 (alternative 5'ss 20-nt downstream), as well as minor (Δ16p93 and Δ16,17p69) and uncharacterized transcripts of 893 and 954 nucleotides. Isoforms Δ16p44, Δ16, Δ16p55, Δ16q4, Δ16q100, and ▾16q20 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.EV’s lab was supported by grants from the Spanish Ministry of Economy and Competitivity, Plan Nacional de I+D+I 2013–2016, ISCIII (Grants: PI13/01749 and PI17/00227) co-funded by FEDER from European Regional Development Funds (European Union), and grant CSI090U14 from the Consejería de Educación (ORDEN EDU/122/2014) and Junta de Castilla y León. EF-B was supported by a predoctoral fellowship from the University of Valladolid and Banco Santander (2015–2019).We acknowledge the support of the publication fee by the CSIC Open Access Publication Support Initiative through its Unit of Information Resources for Research (URICI).Peer reviewe

Table_1_Identification of Eight Spliceogenic Variants in BRCA2 Exon 16 by Minigene Assays.xlsx

<p>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.</p

Association between bat vitamin D receptor 30 haplotypes and vitamin D levels at baseline and a lower response after increased vitamin D supplementation and exposure to sunlight

[Objective]: The aim of this study was to evaluate the relationship between vitamin D levels at baseline and after 12 weeks of supplementation/exposure to sunlight and VDR genotypes (BsmI, TaqI and ApaI) and haplotypes in a homogeneous population of postmenopausal women.[Methods]: We made a prospective study in which 151 women were randomized to two groups: One with 1000 mg of calcium and 800 IU vitamin D supplementation (102 women) and a placebo group with neither calcium or vitamin D supplementation (49 women). The follow-up was from May to September 2012.Vitamin D was determined by chemiluminescent immunoassay. Genotypes were determined using the Sequenomi Plexplatform and haplotypes using PHASE software.[Results]: Baseline (25 ± 10 ng/mlvs.23 ± 9 ng/ml, p > 0.05) and 12-week (32 ± 8 ng/mlvs.29 ± 10 ng/ml, p > 0.05) vitamin D levels were similar between the two groups. The genetic study was made in the total population. There were no differences in baseline and final levels of vitamin D in terms of genotypes and haplotypes, except for the Bat haplotype, whose baseline values were lower (25OHD: 21 ± 10 ng/mlvs. 21 ± 10 ng/ml, p = 0.038). The rate of nonresponders in this group was 15 % (p = 0.001), compared with 9 %, 2 % and 3 % in the other groups.[Conclusions]: The Bat haplotype was associated with lower baseline levels of vitamin D and a worse response to supplementation and, therefore, may be a risk factor for vitamin D deficiency.Peer reviewe