Screening BRCA1 and BRCA2 unclassified variants for splicing mutations using reverse transcription PCR on patient RNA and an ex vivo assay based on a splicing reporter minigene

Detecting splicing patterns in genes involved in hereditary breast and ovarian cancer

Interpretation of variants of unknown significance (VUS) is a major challenge for laboratories performing molecular diagnosis of hereditary breast and ovarian cancer (HBOC), especially considering that many genes are now known to be involved in this syndrome. One important way these VUS can have a functional impact is through their effects on RNA splicing. Here we present a custom RNA-Seq assay plus bioinformatics and biostatistics pipeline to analyse specifically alternative and abnormal splicing junctions in 11 targeted HBOC genes. Our pipeline identified 14 new alternative splices in BRCA1 and BRCA2 in addition to detecting the majority of known alternative spliced transcripts therein. We provide here the first global splicing pattern analysis for the other nine genes, which will enable a comprehensive interpretation of splicing defects caused by VUS in HBOC. Previously known splicing alterations were consistently detected, occasionally with a more complex splicing pattern than expected. We also found that splicing in the 11 genes is similar in blood and breast tissue, supporting the utility and simplicity of blood splicing assays. Our pipeline is ready to be integrated into standard molecular diagnosis for HBOC, but it could equally be adapted for an integrative analysis of any multigene disorder.

Functional assessment of a novel COL4A5 splice region variant and immunostaining of plucked hair follicles as an alternative method of diagnosis in X-linked Alport syndrome

BACKGROUND

Many COL4A5 splice region variants have been described in patients with X-linked Alport syndrome, but few have been confirmed by functional analysis to actually cause defective splicing. We sought to demonstrate that a novel COL4A5 splice region variant in a family with Alport syndrome is pathogenic using functional studies. We also describe an alternative method of diagnosis.

METHODS

Targeted next-generation sequencing results of an individual with Alport syndrome were analyzed and the results confirmed by Sanger sequencing in family members. A splicing reporter minigene assay was used to examine the variant's effect on splicing in transfected cells. Plucked hair follicles from patients and controls were examined for collagen IV proteins using immunofluorescence microscopy.

RESULTS

A novel splice region mutation in COL4A5, c.1780-6T>G, was identified and segregated with disease in this family. This variant caused frequent skipping of exon 25, resulting in a frameshift and truncation of collagen α5(IV) protein. We also developed and validated a new approach to characterize the expression of collagen α5(IV) protein in the basement membranes of plucked hair follicles. Using this approach we demonstrated reduced collagen α5(IV) protein in affected male and female individuals in this family, supporting frequent failure of normal splicing.

CONCLUSIONS

Differing normal to abnormal transcript ratios in affected individuals carrying splice region variants may contribute to variable disease severity observed in Alport families. Examination of plucked hair follicles in suspected X-linked Alport syndrome patients may offer a less invasive alternative method of diagnosis and serve as a pathogenicity test for COL4A5 variants of uncertain significance.

Estimating the prevalence of functional exonic splice regulatory information

In addition to coding information, human exons contain sequences necessary for correct splicing. These elements are known to be under purifying selection and their disruption can cause disease. However, the density of functional exonic splicing information remains profoundly uncertain. Several groups have experimentally investigated how mutations at different exonic positions affect splicing. They have found splice information to be distributed widely in exons, with one estimate putting the proportion of splicing-relevant nucleotides at >90%. These results suggest that splicing could place a major pressure on exon evolution. However, analyses of sequence conservation have concluded that the need to preserve splice regulatory signals only slightly constrains exon evolution, with a resulting decrease in the average human rate of synonymous evolution of only 1–4%. Why do these two lines of research come to such different conclusions? Among other reasons, we suggest that the methods are measuring different things: one assays the density of sites that affect splicing, the other the density of sites whose effects on splicing are visible to selection. In addition, the experimental methods typically consider short exons, thereby enriching for nucleotides close to the splice junction, such sites being enriched for splice-control elements. By contrast, in part owing to correction for nucleotide composition biases and to the assumption that constraint only operates on exon ends, the conservation-based methods can be overly conservative.

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. Fifty-two 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,806 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.7944-7973) 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.7806-9T>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.

A significant proportion of disease-causing mutations of inherited disorders impair splicing. Massive sequencing projects of genetic diseases generate thousands of sequence variations that require functional and clinical interpretations. We have shown that splicing reporter minigenes of the breast cancer genes BRCA1 and BRCA2 are useful tools to functionally test DNA variants. In this work, we have constructed a 7-exon BRCA2 minigene (exons 14 to 20) where we mapped critical splicing regulatory sequences and tested 52 selected variants of exons 17 and 18 detected in breast cancer patients. We finely located three DNA segments on both exons that presumably contain splicing enhancer sequences. We observed that a total of 30 variants of any type disrupted the splicing patterns and, given the severity of their outcomes, we classified 20 of them as pathogenic or likely pathogenic. We also showed that a wide range of splicing elements were affected including canonical and novel 5’ and 3’ splice sites, the polypyrimidine tract and enhancer and silencer sequences. We concluded that splicing aberrations are frequent in Hereditary Breast and Ovarian Cancer and that minigenes are valuable tools to functionally classify DNA variants of any human disease gene under the splicing viewpoint.

Next-Generation Sequencing in Oncology: Genetic Diagnosis, Risk Prediction and Cancer Classification

Next-generation sequencing (NGS) technology has expanded in the last decades with significant improvements in the reliability, sequencing chemistry, pipeline analyses, data interpretation and costs. Such advances make the use of NGS feasible in clinical practice today. This review describes the recent technological developments in NGS applied to the field of oncology. A number of clinical applications are reviewed, i.e., mutation detection in inherited cancer syndromes based on DNA-sequencing, detection of spliceogenic variants based on RNA-sequencing, DNA-sequencing to identify risk modifiers and application for pre-implantation genetic diagnosis, cancer somatic mutation analysis, pharmacogenetics and liquid biopsy. Conclusive remarks, clinical limitations, implications and ethical considerations that relate to the different applications are provided.

A Novel Pathogenic BRCA1 Splicing Variant Produces Partial Intron Retention in the Mature Messenger RNA

About 10% of all breast cancers arise from hereditary mutations that increase the risk of breast and ovarian cancers; and about 25% of these are associated with the BRCA1 or BRCA2 genes. The identification of BRCA1/BRCA2 mutations can enable physicians to better tailor the clinical management of patients; and to initiate preventive measures in healthy carriers. The pathophysiological significance of newly identified variants poses challenges for genetic counseling. We characterized a new BRCA1 variant discovered in a breast cancer patient during BRCA1/2 screening by next-generation sequencing. Bioinformatic predictions; indicating that the variant is probably pathogenetic; were verified using retro-transcription of the patient’s RNA followed by PCR amplifications performed on the resulting cDNA. The variant causes the loss of a canonic donor splice site at position +2 in BRCA1 intron 21; and consequently the partial retention of 156 bp of intron 21 in the patient’s transcript; which demonstrates that this novel BRCA1 mutation plays a pathogenetic role in breast cancer. These findings enabled us to initiate appropriate counseling and to tailor the clinical management of this family. Lastly; these data reinforce the importance of studying the effects of sequence variants at the RNA level to verify their potential role in disease onset.

Functional Studies and In Silico Analyses to Evaluate Non-Coding Variants in Inherited Cardiomyopathies

Point mutations are the most common cause of inherited diseases. Bioinformatics tools can help to predict the pathogenicity of mutations found during genetic screening, but they may work less well in determining the effect of point mutations in non-coding regions. In silico analysis of intronic variants can reveal their impact on the splicing process, but the consequence of a given substitution is generally not predictable. The aim of this study was to functionally test five intronic variants (MYBPC3-c.506-2A>C, MYBPC3-c.906-7G>T, MYBPC3-c.2308+3G>C, SCN5A-c.393-5C>A, and ACTC1-c.617-7T>C) found in five patients affected by inherited cardiomyopathies in the attempt to verify their pathogenic role. Analysis of the MYBPC3-c.506-2A>C mutation in mRNA from the peripheral blood of one of the patients affected by hypertrophic cardiac myopathy revealed the loss of the canonical splice site and the use of an alternative splicing site, which caused the loss of the first seven nucleotides of exon 5 (MYBPC3-G169AfsX14). In the other four patients, we generated minigene constructs and transfected them in HEK-293 cells. This minigene approach showed that MYBPC3-c.2308+3G>C and SCN5A-c.393-5C>A altered pre-mRNA processing, thus resulting in the skipping of one exon. No alterations were found in either MYBPC3-c.906-7G>T or ACTC1-c.617-7T>C. In conclusion, functional in vitro analysis of the effects of potential splicing mutations can confirm or otherwise the putative pathogenicity of non-coding mutations, and thus help to guide the patient's clinical management and improve genetic counseling in affected families.

In silico, in vitro and case-control analyses as an effective combination for analyzing BRCA1 and BRCA2 unclassified variants in a population-based sample

Ascertaining the clinical consequences of BRCA1 and BRCA2 variants of uncertain significance (VUS) is currently indispensable for providing effective genetic counseling and preventive actions for families with hereditary breast and ovarian cancer (HBOC). To this end, we conducted a combination of in silico prediction and cDNA splicing analyses of 13 BRCA1 and 10 BRCA2 VUS identified in our cohort as well as a case-control analysis in a population-based sample of 10 recurrent VUS. We observed consistent results between the in silico predictions and sequencing analyses for all analyzed VUS. An abnormal cDNA pattern was observed for variants c.212+1G>A and c.5278-1G>A in BRCA1 and c.516+2T>A and c.8168A>G in BRCA2 according to in silico splicing prediction. A case-control study of VUS confirmed the polymorphisms of the c.67+62A>G, c.7008-62A>G and c.8851G>A BRCA2 variants previously published. c.4068G>A in the BRCA2 gene can also be considered a polymorphism due to its occurrence at a frequency greater than 1% in our population. Our study shows that employing population-based analysis and a combination of several in silico methods yields highly accurate information, resulting in a reliable tool for selecting variants for cDNA sequencing analysis in routine cancer genetic counseling units.

Detection of Germline Mutation in Hereditary Breast and/or Ovarian Cancers by Next-Generation Sequencing on a Four-Gene Panel

Mutation in BRCA1/BRCA2 genes accounts for 20% of familial breast cancers, 5% to 10% of which may be due to other less penetrant genes which are still incompletely studied. Herein, a four-gene panel was used to examine the prevalence of BRCA1, BRCA2, TP53, and PTEN in hereditary breast and ovarian cancers in Southern Chinese population. In this cohort, 948 high-risk breast and/or ovarian patients were recruited for genetic screening by next-generation sequencing (NGS). The performance of our NGS pipeline was evaluated with 80 Sanger-validated known mutations and eight negative cases. With appropriate bioinformatics analysis pipeline, the detection sensitivity of NGS is comparable with Sanger sequencing. The prevalence of BRCA1/BRCA2 germline mutations was 9.4% in our Chinese cohort, of which 48.8% of the mutations arose from hotspot mutations. With the use of a tailor-made algorithm, HomopolymerQZ, more mutations were detected compared with single mutation detection algorithm. The frequencies of PTEN and TP53 were 0.21% and 0.53%, respectively, in the Southern Chinese patients with breast and/or ovarian cancers. High-throughput NGS approach allows the incorporation of control cohort that provides an ethnicity-specific data for polymorphic variants. Our data suggest that hotspot mutations screening such as SNaPshot could be an effective preliminary screening alternative adopted in a standard clinical laboratory without NGS setup.

Adding In Silico Assessment of Potential Splice Aberration to the Integrated Evaluation of BRCA Gene Unclassified Variants

Clinical mutation screening of the cancer susceptibility genes BRCA1 and BRCA2 generates many unclassified variants (UVs). Most of these UVs are either rare missense substitutions or nucleotide substitutions near the splice junctions of the protein coding exons. Previously, we developed a quantitative method for evaluation of BRCA gene UVs—the “integrated evaluation”—that combines a sequence analysis‐based prior probability of pathogenicity with patient and/or tumor observational data to arrive at a posterior probability of pathogenicity. One limitation of the sequence analysis‐based prior has been that it evaluates UVs from the perspective of missense substitution severity but not probability to disrupt normal mRNA splicing. Here, we calibrated output from the splice‐site fitness program MaxEntScan to generate spliceogenicity‐based prior probabilities of pathogenicity for BRCA gene variants; these range from 0.97 for variants with high probability to damage a donor or acceptor to 0.02 for exonic variants that do not impact a splice junction and are unlikely to create a de novo donor. We created a database http://priors.hci.utah.edu/PRIORS/ that provides the combined missense substitution severity and spliceogenicity‐based probability of pathogenicity for BRCA gene single‐nucleotide substitutions. We also updated the BRCA gene Ex‐UV LOVD, available at http://hci‐exlovd.hci.utah.edu, with 77 re‐evaluable variants.

Naturally occurring BRCA2 alternative mRNA splicing events in clinically relevant samples

BACKGROUND

BRCA1 and BRCA2 are the two principal tumour suppressor genes associated with inherited high risk of breast and ovarian cancer. Genetic testing of BRCA1/2 will often reveal one or more sequence variants of uncertain clinical significance, some of which may affect normal splicing patterns and thereby disrupt gene function. mRNA analyses are therefore among the tests used to interpret the clinical significance of some genetic variants. However, these could be confounded by the appearance of naturally occurring alternative transcripts unrelated to germline sequence variation or defects in gene function. To understand which novel splicing events are associated with splicing mutations and which are part of the normal BRCA2 splicing repertoire, a study was undertaken by members of the Evidence-based Network for the Interpretation of Germline Mutant Alleles (ENIGMA) consortium to characterise the spectrum of naturally occurring BRCA2 mRNA alternate-splicing events.

METHODS

mRNA was prepared from several blood and breast tissue-derived cells and cell lines by contributing ENIGMA laboratories. cDNA representing BRCA2 alternate splice sites was amplified and visualised using capillary or agarose gel electrophoresis, followed by sequencing.

RESULTS

We demonstrate the existence of 24 different BRCA2 mRNA alternate-splicing events in lymphoblastoid cell lines and both breast cancer and non-cancerous breast cell lines.

CONCLUSIONS

These naturally occurring alternate-splicing events contribute to the array of cDNA fragments that may be seen in assays for mutation-associated splicing defects. Caution must be observed in assigning alternate-splicing events to potential splicing mutations.

Compound heterozygous DUOX2 gene mutations (c.2335-1G>C/c.3264_3267delCAGC) associated with congenital hypothyroidism. Characterization of complex cryptic splice sites by minigene analysis

Iodide Organification defects (IOD) represent 10% of cases of congenital hypothyroidism (CH) being the main genes affected that of TPO (thyroid peroxidase) and DUOX2 (dual oxidasa 2). From a patient with clinical and biochemical criteria suggestive with CH associated with IOD, TPO and DUOX2 genes were analyzed by means of PCR-Single Strand Conformation Polymorphism analysis and sequencing. A novel heterozygous compound to the mutations c.2335-1G>C (paternal mutation, intron 17) and c.3264_3267delCAGC (maternal mutation, exon 24) was identified in the DUOX2 gene. Ex-vivo splicing assays and subsequent RT-PCR and sequencing analyses were performed on mRNA isolated from the HeLa cells transfected with wild-type and mutant pSPL3 expression vectors. The wild-type and c.2335-1G>C mutant alleles result in the complete inclusion or exclusion of exon 18, or in the activation of an exonic cryptic 5' ss with the consequent deletion of 169 bp at the end of this exon. However, we observed only a band of the expected size in normal thyroid tissue by RT-PCR. Additionally, the c.2335-1G>C mutation activates an unusual cryptic donor splice site in intron 17, located at position -14 of the authentic intron 17/exon 18 junction site, with an insertion of the last 14 nucleotides of the intron 17 in mutant transcripts with complete and partial inclusion of exon 18. The theoretical consequences of splice site mutation, predicted with the bioinformatics NNSplice, Fsplice, SPL, SPLM and MaxEntScan programs were investigated and evaluated in relation with the experimental evidence. These analyses confirm that c.2335-1G>C mutant allele would result in the abolition of the authentic splice acceptor site. The results suggest the coexistence in our patient of four putative truncated proteins of 786, 805, 806 and 1105 amino acids, with conservation of peroxidase-like domain and loss of gp91(phox)/NOX2-like domain. In conclusion a novel heterozygous compound was identified being responsible of IOD. Cryptic splicing sites have been characterized in DUOX2 gene for the first time. The use of molecular biology techniques is a valuable tool for understanding the molecular pathophysiology of this type of thyroid defects.

Splicing defects caused by exonic mutations in PKD1 as a new mechanism of pathogenesis in autosomal dominant polycystic kidney disease

The correct splicing of precursor-mRNA depends on the actual splice sites plus exonic and intronic regulatory elements recognized by the splicing machinery. Surprisingly, an increasing number of examples reveal that exonic mutations disrupt the binding of splicing factors to these sequences or generate new splice sites or regulatory elements, causing disease. This contradicts the general assumption that missense mutations disrupt protein function and that synonymous mutations are merely polymorphisms. Autosomal dominant polycystic kidney disease (ADPKD) is a common inherited disorder caused mainly by mutations in the PKD1 gene. Recently, we analyzed a substantial number of PKD1 missense or synonymous mutations to further characterize their consequences on pre-mRNA splicing. Our results showed that one missense and 2 synonymous mutations induce significant defects in pre-mRNA splicing. Thus, it appears that aberrant splicing as a result of exonic mutations is a previously unrecognized cause of ADPKD.

Functional classification of BRCA2 DNA variants by splicing assays in a large minigene with 9 exons

Numerous pathogenic DNA variants impair the splicing mechanism in human genetic diseases. Minigenes are optimal approaches to test variants under the splicing viewpoint without the need of patient samples. We aimed to design a robust minigene construct of the breast cancer gene BRCA2 in order to investigate the impact of variants on splicing. BRCA2 exons 19-27 (MGBR2_ex19-27) were cloned in the new vector pSAD. It produced a large transcript of the expected size (2,174 nucleotides) and exon structure (V1-ex19-27-V2). Splicing assays showed that 18 (17 splice-site and 1 silencer variants) out of 40 candidate DNA variants induced aberrant patterns. Twenty-four anomalous transcripts were accurately detected by fluorescent-RT-PCR that were generated by exon-skipping, alternative site usage, and intron-retention events. Fourteen variants induced major anomalies and were predicted to disrupt protein function so they could be classified as pathogenic. Furthermore, minigene mimicked previously reported patient RNA outcomes of seven variants supporting the reproducibility of minigene assays. Therefore, a relevant fraction of variants are involved in breast cancer through splicing alterations. MGBR2_ex19-27 is the largest reported BRCA2 minigene and constitutes a valuable tool for the functional and clinical classification of sequence variations.

Splicing analysis for exonic and intronic mismatch repair gene variants associated with Lynch syndrome confirms high concordance between minigene assays and patient RNA analyses

A subset of DNA variants causes genetic disease through aberrant splicing. Experimental splicing assays, either RT-PCR analyses of patient RNA or functional splicing reporter minigene assays, are required to evaluate the molecular nature of the splice defect. Here, we present minigene assays performed for 17 variants in the consensus splice site regions, 14 exonic variants outside these regions, and two deep intronic variants, all in the DNA mismatch-repair (MMR) genes MLH1, MSH2, MSH6, and PMS2, associated with Lynch syndrome. We also included two deep intronic variants in APC and PKD2. For one variant (MLH1 c.122A>G), our minigene assay and patient RNA analysis could not confirm the previously reported aberrant splicing. The aim of our study was to further investigate the concordance between minigene splicing assays and patient RNA analyses. For 30 variants results from patient RNA analyses were available, either performed by our laboratory or presented in literature. Some variants were deliberately included in this study because they resulted in multiple aberrant transcripts in patient RNA analysis, or caused a splice effect other than the prevalent exon skip. While both methods were completely concordant in the assessment of splice effects, four variants exhibited major differences in aberrant splice patterns. Based on the present and earlier studies, together showing an almost 100% concordance of minigene assays with patient RNA analyses, we discuss the weight given to minigene splicing assays in the current criteria proposed by InSiGHT for clinical classification of MMR variants.

Novel compound heterozygous Thyroglobulin mutations c.745+1G>A/c.7036+2T>A associated with congenital goiter and hypothyroidism in a Vietnamese family. Identification of a new cryptic 5' splice site in the exon 6

Several patients were identified with dyshormonogenesis caused by mutations in the thyroglobulin (TG) gene. These defects are inherited in an autosomal recessive manner and affected individuals are either homozygous or compound heterozygous for the mutations. The aim of the present study was to identify new TG mutations in a patient of Vietnamese origin affected by congenital hypothyroidism, goiter and low levels of serum TG. DNA sequencing identified the presence of compound heterozygous mutations in the TG gene: the maternal mutation consists of a novel c.745+1G>A (g.IVS6 + 1G>A), whereas the hypothetical paternal mutation consists of a novel c.7036+2T>A (g.IVS40 + 2T>A). The father was not available for segregation analysis. Ex-vivo splicing assays and subsequent RT-PCR analyses were performed on mRNA isolated from the eukaryotic-cells transfected with normal and mutant expression vectors. Minigene analysis of the c.745+1G>A mutant showed that the exon 6 is skipped during pre-mRNA splicing or partially included by use of a cryptic 5' splice site located to 55 nucleotides upstream of the authentic exon 6/intron 6 junction site. The functional analysis of c.7036+2T>A mutation showed a complete skipping of exon 40. The theoretical consequences of splice site mutations, predicted with the bioinformatics tool NNSplice, Fsplice, SPL, SPLM and MaxEntScan programs were investigated and evaluated in relation with the experimental evidence. These analyses predicted that both mutant alleles would result in the abolition of the authentic splice donor sites. The c.745+1G>A mutation originates two putative truncated proteins of 200 and 1142 amino acids, whereas c.7036+2T>A mutation results in a putative truncated protein of 2277 amino acids. In conclusion, we show that the c.745+1G>A mutation promotes the activation of a new cryptic donor splice site in the exon 6 of the TG gene. The functional consequences of these mutations could be structural changes in the protein molecule that alter the biosynthesis of thyroid hormones.

Splicing analysis of 14 BRCA1 missense variants classifies nine variants as pathogenic

Pathogenic germline mutations in the BRCA1 gene predispose carriers to early onset breast and ovarian cancer. Clinical genetic screening of BRCA1 often reveals variants with uncertain clinical significance, complicating patient and family management. Therefore, functional examinations are urgently needed to classify whether these uncertain variants are pathogenic or benign. In this study, we investigated 14 BRCA1 variants by in silico splicing analysis and mini-gene splicing assay. All 14 alterations were missense variants located within the BRCT domain of BRCA1 and had previously been examined by functional analysis at the protein level. Results from a validated mini-gene splicing assay indicated that nine BRCA1 variants resulted in splicing aberrations leading to truncated transcripts and thus can be considered pathogenic (c.4987A>T/p.Met1663Leu, c.4988T>A/p.Met1663Lys, c.5072C>T/p.Thr1691Ile, c.5074G>C/p.Asp1692His, c.5074G>A/p.Asp1692Asn, c.5074G>T/p.Asp1692Tyr, c.5332G>A/p.Asp1778Asn, c.5332G>T/p.Asp1778Tyr, and c.5408G>C/p.Gly1803Ala), whereas five BRCA1 variants had no effect on splicing (c.4985T>C/p.Phe1662Ser, c.5072C>A/p.Thr1691Lys, c.5153G>C/p.Trp1718Ser, c.5154G>T/p.Trp1718Cys, and c.5333A>G/p.Asp1778Gly). Eight of the variants having an effect on splicing (c.4987A>T/p.Met1663Leu, c.4988T>A/p.Met1663Lys, c.5074G>C/p.Asp1692His, c.5074G>A/p.Asp1692Asn, c.5074G>T/p.Asp1692Tyr, c.5332G>A/p.Asp1778Asn, c.5332G>T/p.Asp1778Tyr, and c.5408G>C/p.Gly1803Ala) were previously determined to have no or an uncertain effect on the protein level, whereas one variant (c.5072C>T/p.Thr1691Ile) were shown to have a strong effect on the protein level as well. In conclusion, our study emphasizes that in silico splicing prediction and mini-gene splicing analysis are important for the classification of BRCA1 missense variants located close to exon/intron boundaries.

Characterization of three alternative transcripts of the BRCA1 gene in patients with breast cancer and a family history of breast and/or ovarian cancer who tested negative for pathogenic mutations

The study of BRCA1 and BRCA2 genes and their alterations has been essential to the understanding of the development of familial breast and ovarian cancers. Many of the variants identified have an unknown pathogenic significance. These include variants which determine alternative mRNA splicing, identified in the intronic regions and those are capable of destroying the splicing ability. The aim of this study was to detect BRCA1/BRCA2 aberrant transcripts resulting from alternative splicing, in women with a known family history and/or early onset of breast and/or ovarian cancer, tested wild-type for BRCA1 and BRCA2. The identification and characterization of aberrant transcripts through the analysis of mRNA levels in blood lymphocytes may help us to recognize families otherwise misclassified as wild-type BRCA1 and BRCA2. Blood samples were collected from 13 women that had a family history of breast and/or ovarian cancer and tested negative for pathogenic mutations in the BRCA1 and BRCA2 genes. Total RNA was analyzed for the presence of BRCA1 and BRCA2 naturally occuring and pathological transcripts using RT-PCR. In 2 out of the 13 samples, 2 alternative transcripts of the BRCA1 gene were identified. These were probably pathogenic as they lacked exon 17 and exon 15, respectively, giving rise to a truncated protein. In addition to these, we identified the Δ17–19 transcript in 1 patient, which gives rise to a protein with an in-frame deletion of 69 amino acids. In conclusion, this study on alternative transcripts of the BRCA1 and BRCA2 genes revealed the presence of isoforms (prevalence of 15%) in blood samples from women with breast and ovarian cancer that were probably pathogenic, that were not detected by conventional methods of mutation screening based on direct sequencing of all coding regions, intron-exons junctions and MLPA analysis.

Detection of BRCA1 and BRCA2 germline mutations in Japanese population using next-generation sequencing

Tumor suppressor genes BRCA1 and BRCA2 are the two main breast and ovarian cancer susceptibility genes, and their genetic testing has been used to evaluate the risk of hereditary breast and ovarian cancer (HBOC). While several studies have reported the prevalence of BRCA1 and BRCA2 mutations in Japanese populations, there is insufficient information about deleterious mutations compared with western countries. Moreover, because many rare variants are found in BRCA1 and BRCA2, both of which encode large proteins, it is difficult to sequence all coding regions using the Sanger method for mutation detection. In this study, therefore, we performed next-generation sequencing (NGS) analysis of the entire coding regions of BRCA1 and BRCA2 in 135 breast and/or ovarian cancer patients. Deleterious BRCA1 and BRCA2 mutations were detected in 10 patients (7.4%) by NGS analysis. Of these, one mutation in BRCA1 and two in BRCA2 had not been reported previously. Furthermore, a BRCA2 mutation found in a proband was also identified in two unaffected relatives. These data suggest the utility of screening BRCA1 and BRCA2 mutations by NGS in clinical diagnosis.

Experimental assessment of splicing variants using expression minigenes and comparison with in silico predictions

Assessment of the functional consequences of variants near splice sites is a major challenge in the diagnostic laboratory. To address this issue, we created expression minigenes (EMGs) to determine the RNA and protein products generated by splice site variants (n = 10) implicated in cystic fibrosis (CF). Experimental results were compared with the splicing predictions of eight in silico tools. EMGs containing the full-length Cystic Fibrosis Transmembrane Conductance Regulator (CFTR) coding sequence and flanking intron sequences generated wild-type transcript and fully processed protein in Human Embryonic Kidney (HEK293) and CF bronchial epithelial (CFBE41o-) cells. Quantification of variant induced aberrant mRNA isoforms was concordant using fragment analysis and pyrosequencing. The splicing patterns of c.1585-1G>A and c.2657+5G>A were comparable to those reported in primary cells from individuals bearing these variants. Bioinformatics predictions were consistent with experimental results for 9/10 variants (MES), 8/10 variants (NNSplice), and 7/10 variants (SSAT and Sroogle). Programs that estimate the consequences of mis-splicing predicted 11/16 (HSF and ASSEDA) and 10/16 (Fsplice and SplicePort) experimentally observed mRNA isoforms. EMGs provide a robust experimental approach for clinical interpretation of splice site variants and refinement of in silico tools.

Predicting survival in head and neck squamous cell carcinoma from TP53 mutation

For TP53-mutated head and neck squamous cell carcinomas (HNSCCs), the codon and specific amino acid sequence change resulting from a patient's mutation can be prognostic. Thus, developing a framework to predict patient survival for specific mutations in TP53 would be valuable. There are many bioinformatics and functional methods for predicting the phenotypic impact of genetic variation, but their overall clinical value remains unclear. Here, we assess the ability of 15 different methods to predict HNSCC patient survival from TP53 mutation, using TP53 mutation and clinical data from patients enrolled in E4393 by the Eastern Cooperative Oncology Group (ECOG), which investigated whether TP53 mutations in surgical margins were predictive of disease recurrence. These methods include: server-based computational tools SIFT, PolyPhen-2, and Align-GVGD; our in-house POSE and VEST algorithms; the rules devised in Poeta et al. with and without considerations for splice-site mutations; location of mutation in the DNA-bound TP53 protein structure; and a functional assay measuring WAF1 transactivation in TP53-mutated yeast. We assessed method performance using overall survival (OS) and progression-free survival (PFS) from 420 HNSCC patients, of whom 224 had TP53 mutations. Each mutation was categorized as "disruptive" or "non-disruptive". For each method, we compared the outcome between the disruptive group vs. the non-disruptive group. The rules devised by Poeta et al. with or without our splice-site modification were observed to be superior to others. While the differences in OS (disruptive vs. non-disruptive) appear to be marginally significant (Poeta rules + splice rules, P = 0.089; Poeta rules, P = 0.053), both algorithms identified the disruptive group as having significantly worse PFS outcome (Poeta rules + splice rules, P = 0.011; Poeta rules, P = 0.027). In general, prognostic performance was low among assessed methods. Further studies are required to develop and validate methods that can predict functional and clinical significance of TP53 mutations in HNSCC patients.

Functional characterization of BRCA1 gene variants by mini-gene splicing assay

Mutational screening of the breast cancer susceptibility gene BRCA1 leads to the identification of numerous pathogenic variants such as frameshift and nonsense variants, as well as large genomic rearrangements. The screening moreover identifies a large number of variants, for example, missense, silent, and intron variants, which are classified as variants of unknown clinical significance owing to the lack of causal evidence. Variants of unknown clinical significance can potentially have an impact on splicing and therefore functional examinations are warranted to classify whether these variants are pathogenic or benign. Here we validate a mini-gene splicing assay by comparing the results of 24 variants with previously published data from RT-PCR analysis on RNA from blood samples/lymphoblastoid cell lines. The analysis showed an overall concordance of 100%. In addition, we investigated 13 BRCA1 variants of unknown clinical significance or putative variants affecting splicing by in silico analysis and mini-gene splicing assay. Both the in silico analysis and mini-gene splicing assay classified six BRCA1 variants as pathogenic (c.80+1G>A, c.132C>T (p.=), c.213-1G>A, c.670+1delG, c.4185+1G>A, and c.5075-1G>C), whereas six BRCA1 variants were classified as neutral (c.-19-22_-19-21dupAT, c.302-15C>G, c.547+14delG, c.4676-20A>G, c.4987-21G>T, and c.5278-14C>G) and one BRCA1 variant remained unclassified (c.670+16G>A). In conclusion, our study emphasizes that in silico analysis and mini-gene splicing assays are important for the classification of variants, especially if no RNA is available from the patient. This knowledge is crucial for proper genetic counseling of patients and their family members.

In vitro analysis of splice site mutations in the CLCN1 gene using the minigene assay

Mutations in the chloride channel gene CLCN1 cause the allelic disorders Thomsen (dominant) and Becker (recessive) myotonia congenita (MC). The encoded protein, ClC-1, is the primary channel that mediates chloride (Cl-) conductance in skeletal muscle. Mutations in CLCN1 lower the channel's threshold voltage, leading to spontaneous action potentials that are not coupled to neuromuscular transmission and resulting in myotonia. Over 120 mutations in CLCN1 have been described, 10% of which are splicing defects. Biological specimens suitable for RNA extraction are not always available, but obtaining genomic DNA for analysis is easy and non-invasive. This is the first study to evaluate the pathogenic potential of novel splicing mutations using the minigene approach, which is based on genomic DNA analysis. Splicing mutations accounted for 23% of all pathogenic variants in our cohort of MC patients. Four were heterozygous mutations in four unrelated individuals, belonging to this cohort: c.563G>T in exon 5; c.1169-5T>G in intron 10; c.1251+1G>A in intron 11, and c.1931-2A>G in intron 16. These variants were expressed in HEK 293 cells, and aberrant splicing was verified by in vitro transcription and sequencing of the cDNA. Our findings confirm the need to further investigate the nature of rearrangements associated with this class of mutations and their effects on mature transcripts. In particular, splicing mutations predicted to generate in-frame transcripts may generate out-of-frame mRNA transcripts that do not produce functional ClC-1. Clinically, incomplete molecular evaluation could lead to delayed or faulty diagnosis.

Comparison of mRNA splicing assay protocols across multiple laboratories: recommendations for best practice in standardized clinical testing

BACKGROUND

Accurate evaluation of unclassified sequence variants in cancer predisposition genes is essential for clinical management and depends on a multifactorial analysis of clinical, genetic, pathologic, and bioinformatic variables and assays of transcript length and abundance. The integrity of assay data in turn relies on appropriate assay design, interpretation, and reporting.

METHODS

We conducted a multicenter investigation to compare mRNA splicing assay protocols used by members of the ENIGMA (Evidence-Based Network for the Interpretation of Germline Mutant Alleles) consortium. We compared similarities and differences in results derived from analysis of a panel of breast cancer 1, early onset (BRCA1) and breast cancer 2, early onset (BRCA2) gene variants known to alter splicing (BRCA1: c.135-1G>T, c.591C>T, c.594-2A>C, c.671-2A>G, and c.5467+5G>C and BRCA2: c.426-12_8delGTTTT, c.7988A>T, c.8632+1G>A, and c.9501+3A>T). Differences in protocols were then assessed to determine which elements were critical in reliable assay design.

RESULTS

PCR primer design strategies, PCR conditions, and product detection methods, combined with a prior knowledge of expected alternative transcripts, were the key factors for accurate splicing assay results. For example, because of the position of primers and PCR extension times, several isoforms associated with BRCA1, c.594-2A>C and c.671-2A>G, were not detected by many sites. Variation was most evident for the detection of low-abundance transcripts (e.g., BRCA2 c.8632+1G>A Δ19,20 and BRCA1 c.135-1G>T Δ5q and Δ3). Detection of low-abundance transcripts was sometimes addressed by using more analytically sensitive detection methods (e.g., BRCA2 c.426-12_8delGTTTT ins18bp).

CONCLUSIONS

We provide recommendations for best practice and raise key issues to consider when designing mRNA assays for evaluation of unclassified sequence variants.

Capillary electrophoresis analysis of conventional splicing assays: IARC analytical and clinical classification of 31 BRCA2 genetic variants

Rare sequence variants in "high-risk" disease genes, often referred as unclassified variants (UVs), pose a serious challenge to genetic testing. However, UVs resulting in splicing alterations can be readily assessed by in vitro assays. Unfortunately, analytical and clinical interpretation of these assays is often challenging. Here, we explore this issue by conducting splicing assays in 31 BRCA2 genetic variants. All variants were assessed by RT-PCR followed by capillary electrophoresis and direct sequencing. If assays did not produce clear-cut outputs (Class-2 or Class-5 according to analytical International Agency for Research on Cancer guidelines), we performed qPCR and/or minigene assays. The latter were performed with a new splicing vector (pSAD) developed by authors of the present manuscript (patent #P201231427 CSIC). We have identified three clinically relevant Class-5 variants (c.682-2A>G, c.7617+1G>A, and c.8954-5A>G), and 27 analytical Class-2 variants (not inducing splicing alterations). In addition, we demonstrate that rs9534262 (c.7806-14T>C) is a BRCA2 splicing quantitative trait locus.

SASD: the Synthetic Alternative Splicing Database for identifying novel isoform from proteomics

BACKGROUND

Alternative splicing is an important and widespread mechanism for generating protein diversity and regulating protein expression. High-throughput identification and analysis of alternative splicing in the protein level has more advantages than in the mRNA level. The combination of alternative splicing database and tandem mass spectrometry provides a powerful technique for identification, analysis and characterization of potential novel alternative splicing protein isoforms from proteomics.

RESULTS

We used a three-step pipeline to create a synthetic alternative splicing database (SASD) to identify novel alternative splicing isoforms and interpret them at the context of pathway, disease, drug and organ specificity or custom gene set with maximum coverage and exclusive focus on alternative splicing. First, we extracted information on gene structures of all genes in the Ensembl Genes 71 database and incorporated the Integrated Pathway Analysis Database. Then, we compiled artificial splicing transcripts. Lastly, we translated the artificial transcripts into alternative splicing peptides.

CONCLUSIONS

The SASD provides the scientific community with an efficient means to identify, analyze, and characterize novel Exon Skipping and Intron Retention protein isoforms from mass spectrometry and interpret them at the context of pathway, disease, drug and organ specificity or custom gene set with maximum coverage and exclusive focus on alternative splicing.

Functional examination of MLH1, MSH2, and MSH6 intronic mutations identified in Danish colorectal cancer patients

Background

Germ-line mutations in the DNA mismatch repair genes MLH1, MSH2, and MSH6 predispose to the development of colorectal cancer (Lynch syndrome or hereditary nonpolyposis colorectal cancer). These mutations include disease-causing frame-shift, nonsense, and splicing mutations as well as large genomic rearrangements. However, a large number of mutations, including missense, silent, and intronic variants, are classified as variants of unknown clinical significance.

Methods

Intronic MLH1, MSH2, or MSH6 variants were investigated using in silico prediction tools and mini-gene assay to asses the effect on splicing.

Results

We describe in silico and in vitro characterization of nine intronic MLH1, MSH2, or MSH6 mutations identified in Danish colorectal cancer patients, of which four mutations are novel. The analysis revealed aberrant splicing of five mutations (MLH1 c.588 + 5G > A, MLH1 c.677 + 3A > T, MLH1 c.1732-2A > T, MSH2 c.1276 + 1G > T, and MSH2 c.1662-2A > C), while four mutations had no effect on splicing compared to wild type (MLH1 c.117-34A > T, MLH1 c.1039-8 T > A, MSH2 c.2459-18delT, and MSH6 c.3439-16C > T).

Conclusions

In conclusion, we classify five MLH1/MSH2 mutations as pathogenic, whereas four MLH1/MSH2/MSH6 mutations are classified as neutral. This study supports the notion that in silico prediction tools and mini-gene assays are important for the classification of intronic variants, and thereby crucial for the genetic counseling of patients and their family members.

Biallelic DICER1 mutations occur in Wilms tumours

DICER1 is an endoribonuclease central to the generation of microRNAs (miRNAs) and short interfering RNAs (siRNAs). Germline mutations in DICER1 have been associated with a pleiotropic tumour predisposition syndrome and Wilms tumour (WT) is a rare manifestation of this syndrome. Three WTs, each in a child with a deleterious germline DICER1 mutation, were screened for somatic DICER1 mutations and were found to bear specific mutations in either the RNase IIIa (n = 1) or the RNase IIIb domain (n = 2). In the two latter cases, we demonstrate that the germline and somatic DICER1 mutations were in trans, suggesting that the two-hit hypothesis of tumour formation applies for these examples of WT. Among 191 apparently sporadic WTs, we identified five different missense or deletion somatic DICER1 mutations (2.6%) in four individual WTs; one tumour had two very likely deleterious somatic mutations in trans in the RNase IIIb domain (c.5438A>G and c.5452G>A). In vitro studies of two somatic single-base substitutions (c.5429A>G and c.5438A>G) demonstrated exon 25 skipping from the transcript, a phenomenon not previously reported in DICER1. Further we show that DICER1 transcripts lacking exon 25 can be translated in vitro. This study has demonstrated that a subset of WTs exhibits two 'hits' in DICER1, suggesting that these mutations could be key events in the pathogenesis of these tumours.

Comparative in vitro and in silico analyses of variants in splicing regions of BRCA1 and BRCA2 genes and characterization of novel pathogenic mutations

Several unclassified variants (UVs) have been identified in splicing regions of disease-associated genes and their characterization as pathogenic mutations or benign polymorphisms is crucial for the understanding of their role in disease development. In this study, 24 UVs located at BRCA1 and BRCA2 splice sites were characterized by transcripts analysis. These results were used to evaluate the ability of nine bioinformatics programs in predicting genetic variants causing aberrant splicing (spliceogenic variants) and the nature of aberrant transcripts. Eleven variants in BRCA1 and 8 in BRCA2, including 8 not previously characterized at transcript level, were ascertained to affect mRNA splicing. Of these, 16 led to the synthesis of aberrant transcripts containing premature termination codons (PTCs), 2 to the up-regulation of naturally occurring alternative transcripts containing PTCs, and one to an in-frame deletion within the region coding for the DNA binding domain of BRCA2, causing the loss of the ability to bind the partner protein DSS1 and ssDNA. For each computational program, we evaluated the rate of non-informative analyses, i.e. those that did not recognize the natural splice sites in the wild-type sequence, and the rate of false positive predictions, i.e., variants incorrectly classified as spliceogenic, as a measure of their specificity, under conditions setting sensitivity of predictions to 100%. The programs that performed better were Human Splicing Finder and Automated Splice Site Analyses, both exhibiting 100% informativeness and specificity. For 10 mutations the activation of cryptic splice sites was observed, but we were unable to derive simple criteria to select, among the different cryptic sites predicted by the bioinformatics analyses, those actually used. Consistent with previous reports, our study provides evidences that in silico tools can be used for selecting splice site variants for in vitro analyses. However, the latter remain mandatory for the characterization of the nature of aberrant transcripts.

Analysis of 30 putative BRCA1 splicing mutations in hereditary breast and ovarian cancer families identifies exonic splice site mutations that escape in silico prediction

Screening for pathogenic mutations in breast and ovarian cancer genes such as BRCA1/2, CHEK2 and RAD51C is common practice for individuals from high-risk families. However, test results may be ambiguous due to the presence of unclassified variants (UCV) in the concurrent absence of clearly cancer-predisposing mutations. Especially the presence of intronic or exonic variants within these genes that possibly affect proper pre-mRNA processing poses a challenge as their functional implications are not immediately apparent. Therefore, it appears necessary to characterize potential splicing UCV and to develop appropriate classification tools. We investigated 30 distinct BRCA1 variants, both intronic and exonic, regarding their spliceogenic potential by commonly used in silico prediction algorithms (HSF, MaxEntScan) along with in vitro transcript analyses. A total of 25 variants were identified spliceogenic, either causing/enhancing exon skipping or activation of cryptic splice sites, or both. Except from a single intronic variant causing minor effects on BRCA1 pre-mRNA processing in our analyses, 23 out of 24 intronic variants were correctly predicted by MaxEntScan, while HSF was less accurate in this cohort. Among the 6 exonic variants analyzed, 4 severely impair correct pre-mRNA processing, while the remaining two have partial effects. In contrast to the intronic alterations investigated, only half of the spliceogenic exonic variants were correctly predicted by HSF and/or MaxEntScan. These data support the idea that exonic splicing mutations are commonly disease-causing and concurrently prone to escape in silico prediction, hence necessitating experimental in vitro splicing analysis.

Analysis of the effects of rare variants on splicing identifies alterations in GABAA receptor genes in autism spectrum disorder individuals

A large-scale sequencing screen of X-linked synaptic genes in individuals with autism spectrum disorder (ASD) or schizophrenia (SCZ), two common neurodevelopmental disorders, identified many variants most of which have no easily predictable effect on gene function. In this report, we evaluated the impact of these rare missense and silent variants on gene splicing. For this purpose, we used complementary in silico analyses, in vitro minigene-based assays and RNA prepared from lymphoblastoid cells derived from patients with these mutations. Our goal was to identify the variants which might either create or disrupt an acceptor splice site, a donor splice site or an exonic splicing enhancer, thus leading to aberrant splicing that could be involved in the pathogenesis of ASD or SCZ. We identified truncating mutations in distinct X-linked gamma-aminobutyric acid A (GABAA) receptor subunit-encoding genes, GABRQ and GABRA3, in two different families. Furthermore, missense and silent variants in nuclear RNA export factor 5 and histone deacetylase 6 were shown to partially disrupt the protein. While genes from the GABAergic pathway have previously been thought to be involved in the pathophysiology of ASD, this is the first report of ASD patients with truncating mutations in GABA receptors genes.

Whole-exome sequencing studies of nonhereditary (sporadic) parathyroid adenomas

CONTEXT

Genetic abnormalities, such as those of multiple endocrine neoplasia type 1 (MEN1) and Cyclin D1 (CCND1) genes, occur in <50% of nonhereditary (sporadic) parathyroid adenomas.

OBJECTIVE

To identify genetic abnormalities in nonhereditary parathyroid adenomas by whole-exome sequence analysis.

DESIGN

Whole-exome sequence analysis was performed on parathyroid adenomas and leukocyte DNA samples from 16 postmenopausal women without a family history of parathyroid tumors or MEN1 and in whom primary hyperparathyroidism due to single-gland disease was cured by surgery. Somatic variants confirmed in this discovery set were assessed in 24 other parathyroid adenomas.

RESULTS

Over 90% of targeted exons were captured and represented by more than 10 base reads. Analysis identified 212 somatic variants (median eight per tumor; range, 2-110), with the majority being heterozygous nonsynonymous single-nucleotide variants that predicted missense amino acid substitutions. Somatic MEN1 mutations occurred in six of 16 (∼35%) parathyroid adenomas, in association with loss of heterozygosity on chromosome 11. However, no other gene was mutated in more than one tumor. Mutations in several genes that may represent low-frequency driver mutations were identified, including a protection of telomeres 1 (POT1) mutation that resulted in exon skipping and disruption to the single-stranded DNA-binding domain, which may contribute to increased genomic instability and the observed high mutation rate in one tumor.

CONCLUSIONS

Parathyroid adenomas typically harbor few somatic variants, consistent with their low proliferation rates. MEN1 mutation represents the major driver in sporadic parathyroid tumorigenesis although multiple low-frequency driver mutations likely account for tumors not harboring somatic MEN1 mutations.

Comprehensive splicing functional analysis of DNA variants of the BRCA2 gene by hybrid minigenes

Introduction

The underlying pathogenic mechanism of a large fraction of DNA variants of disease-causing genes is the disruption of the splicing process. We aimed to investigate the effect on splicing of the BRCA2 variants c.8488-1G > A (exon 20) and c.9026_9030del (exon 23), as well as 41 BRCA2 variants reported in the Breast Cancer Information Core (BIC) mutation database.

Methods

DNA variants were analyzed with the splicing prediction programs NNSPLICE and Human Splicing Finder. Functional analyses of candidate variants were performed by lymphocyte RT-PCR and/or hybrid minigene assays. Forty-one BIC variants of exons 19, 20, 23 and 24 were bioinformatically selected and generated by PCR-mutagenesis of the wild type minigenes.

Results

Lymphocyte RT-PCR of c.8488-1G > A showed intron 19 retention and a 12-nucleotide deletion in exon 20, whereas c.9026_9030del did not show any splicing anomaly. Minigene analysis of c.8488-1G > A displayed the aforementioned aberrant isoforms but also exon 20 skipping. We further evaluated the splicing outcomes of 41 variants of four BRCA2 exons by minigene analysis. Eighteen variants presented splicing aberrations. Most variants (78.9%) disrupted the natural splice sites, whereas four altered putative enhancers/silencers and had a weak effect. Fluorescent RT-PCR of minigenes accurately detected 14 RNA isoforms generated by cryptic site usage, exon skipping and intron retention events. Fourteen variants showed total splicing disruptions and were predicted to truncate or eliminate essential domains of BRCA2.

Conclusions

A relevant proportion of BRCA2 variants are correlated with splicing disruptions, indicating that RNA analysis is a valuable tool to assess the pathogenicity of a particular DNA change. The minigene system is a straightforward and robust approach to detect variants with an impact on splicing and contributes to a better knowledge of this gene expression step.

SEPTIN12 genetic variants confer susceptibility to teratozoospermia

It is estimated that 10-15% of couples are infertile and male factors account for about half of these cases. With the advent of intracytoplasmic sperm injection (ICSI), many infertile men have been able to father offspring. However, teratozoospermia still remains a big challenge to tackle. Septins belong to a family of cytoskeletal proteins with GTPase activity and are involved in various biological processes e.g. morphogenesis, compartmentalization, apoptosis and cytokinesis. SEPTIN12, identified by c-DNA microarray analysis of infertile men, is exclusively expressed in the post meiotic male germ cells. Septin12(+/+)/Septin12(+/-) chimeric mice have multiple reproductive defects including the presence of immature sperm in the semen, and sperm with bent neck (defect of the annulus) and nuclear DNA damage. These facts make SEPTIN12 a potential sterile gene in humans. In this study, we sequenced the entire coding region of SEPTIN12 in infertile men (n = 160) and fertile controls (n = 200) and identified ten variants. Among them is the c.474 G>A variant within exon 5 that encodes part of the GTP binding domain. The variant creates a novel splice donor site that causes skipping of a portion of exon 5, resulting in a truncated protein lacking the C-terminal half of SEPTIN12. Most individuals homozygous for the c.474 A allele had teratozoospermia (abnormal sperm <14%) and their sperm showed bent tail and de-condensed nucleus with significant DNA damage. Ex vivo experiment showed truncated SEPT12 inhibits filament formation in a dose-dependent manner. This study provides the first causal link between SEPTIN12 genetic variant and male infertility with distinctive sperm pathology. Our finding also suggests vital roles of SEPT12 in sperm nuclear integrity and tail development.

Detection of exon skipping events in BRCA1 RNA using MLPA kit P002

A rapid and easy method to screen for aberrant cDNA would be a very useful diagnostic tool in genetics since a fraction of the DNA variants found affect RNA splicing. The currently used RT-PCR methods require new primer combinations to study each variant that might affect splicing. Since MLPA is routinely used to detect large genomic deletions and successfully detected exon skipping events in Duchenne muscular dystrophy in cDNA, we performed a pilot study to evaluate its value for BRCA1 cDNA. The effect of puromycin, DNase I and two different DNA cleaning protocols were tested in the RNA analysis of lymphocyte cultures. We used two samples from unrelated families with two different BRCA1 exon deletion events, two healthy unrelated controls and six samples from hereditary breast/ovarian cancer syndrome (HBOC) patients without BRCA1/2 mutations. Using RNA treated with DNase I and cleaned in a column system from puromycin-treated fractions, we were able to identify the two BRCA1 deletions. Additional HBOC patients did not show additional splice events. However, we were not able to get reproducible results. Therefore, the cDNA-MLPA technique using kit BRCA1 P002 is in our hands currently not reliable enough for routine RNA analysis and needs further optimization.

An entire exon 3 germ-line rearrangement in the BRCA2 gene: pathogenic relevance of exon 3 deletion in breast cancer predisposition

Background

Germ-line mutations in the BRCA1 and BRCA2 genes are major contributors to hereditary breast/ovarian cancer. Large rearrangements are less frequent in the BRCA2 gene than in BRCA1. We report, here, the first total deletion of exon 3 in the BRCA2 gene that was detected during screening of 2058 index cases from breast/ovarian cancer families for BRCA2 large rearrangements. Deletion of exon 3, which is in phase, does not alter the reading frame. Low levels of alternative transcripts lacking exon 3 (Δ3 delta3 transcript) have been reported in normal tissues, which raises the question whether deletion of exon 3 is pathogenic.

Methods

Large BRCA2 rearrangements were analysed by QMPSF (Quantitative Multiplex PCR of Short Fluorescent Fragments) or MLPA (Multiplex Ligation-Dependent Probe Amplification). The exon 3 deletion was characterized with a "zoom-in" dedicated CGH array to the BRCA2 gene and sequencing. To determine the effect of exon 3 deletion and assess its pathogenic effect, three methods of transcript quantification were used: fragment analysis of FAM-labelled PCR products, specific allelic expression using an intron 2 polymorphism and competitive quantitative RT-PCR.

Results

Large rearrangements of BRCA2 were detected in six index cases out of 2058 tested (3% of all deleterious BRCA2 mutations). This study reports the first large rearrangement of the BRCA2 gene that includes all of exon 3 and leads to an in frame deletion of exon 3 at the transcriptional level. Thirty five variants in exon 3 and junction regions of BRCA2 are also reported, that contribute to the interpretation of the pathogenicity of the deletion. The quantitative approaches showed that there are three classes of delta3 BRCA2 transcripts (low, moderate and exclusive). Exclusive expression of the delta3 transcript by the mutant allele and segregation data provide evidence for a causal effect of the exon 3 deletion.

Conclusion

This paper highlights that large rearrangements and total deletion of exon 3 in the BRCA2 gene could contribute to hereditary breast and/or ovarian cancer. In addition, our findings suggest that, to interpret the pathogenic effect of any variants of exon 3, both accurate transcript quantification and co-segregation analysis are required.

Characterization of BRCA1 and BRCA2 splicing variants: a collaborative report by ENIGMA consortium members

Mutations in BRCA1 and BRCA2 predispose carriers to early onset breast and ovarian cancer. A common problem in clinical genetic testing is interpretation of variants with unknown clinical significance. The Evidence-based Network for the Interpretation of Germline Mutant Alleles (ENIGMA) consortium was initiated to evaluate and implement strategies to characterize the clinical significance of BRCA1 and BRCA2 variants. As an initial project of the ENIGMA Splicing Working Group, we report splicing and multifactorial likelihood analysis of 25 BRCA1 and BRCA2 variants from seven different laboratories. Splicing analysis was performed by reverse transcriptase PCR or mini gene assay, and sequencing to identify aberrant transcripts. The findings were compared to bioinformatic predictions using four programs. The posterior probability of pathogenicity was estimated using multifactorial likelihood analysis, including co-occurrence with a deleterious mutation, segregation and/or report of family history. Abnormal splicing patterns expected to lead to a non-functional protein were observed for 7 variants (BRCA1 c.441+2T>A, c.4184_4185+2del, c.4357+1G>A, c.4987-2A>G, c.5074G>C, BRCA2 c.316+5G>A, and c.8754+3G>C). Combined interpretation of splicing and multifactorial analysis classified an initiation codon variant (BRCA2 c.3G>A) as likely pathogenic, uncertain clinical significance for 7 variants, and indicated low clinical significance or unlikely pathogenicity for another 10 variants. Bioinformatic tools predicted disruption of consensus donor or acceptor sites with high sensitivity, but cryptic site usage was predicted with low specificity, supporting the value of RNA-based assays. The findings also provide further evidence that clinical RNA-based assays should be extended from analysis of invariant dinucleotides to routinely include all variants located within the donor and acceptor consensus splicing sites. Importantly, this study demonstrates the added value of collaboration between laboratories, and across disciplines, to collate and interpret information from clinical testing laboratories to consolidate patient management.

Contribution of bioinformatics predictions and functional splicing assays to the interpretation of unclassified variants of the BRCA genes

A large fraction of sequence variants of unknown significance (VUS) of the breast and ovarian cancer susceptibility genes BRCA1 and BRCA2 may induce splicing defects. We analyzed 53 VUSs of BRCA1 or BRCA2, detected in consecutive molecular screenings, by using five splicing prediction programs, and we classified them into two groups according to the strength of the predictions. In parallel, we tested them by using functional splicing assays. A total of 10 VUSs were predicted by two or more programs to induce a significant reduction of splice site strength or activation of cryptic splice sites or generation of new splice sites. Minigene-based splicing assays confirmed four of these predictions. Five additional VUSs, all at internal exon positions, were not predicted to induce alterations of splice sites, but revealed variable levels of exon skipping, most likely induced by the modification of exonic splicing regulatory elements. We provide new data in favor of the pathogenic nature of the variants BRCA1 c.212+3A>G and BRCA1 c.5194-12G>A, which induced aberrant out-of-frame mRNA forms. Moreover, the novel variant BRCA2 c.7977-7C>G induced in frame inclusion of 6 nt from the 3' end of intron 17. The novel variants BRCA2 c.520C>T and BRCA2 c.7992T>A induced incomplete skipping of exons 7 and 18, respectively. This work highlights the contribution of splicing minigene assays to the assessment of pathogenicity, not only when patient RNA is not available, but also as a tool to improve the accuracy of bioinformatics predictions.

Characterisation of unclassified variants in the BRCA1/2 genes with a putative effect on splicing

A subset of the unclassified variants (UVs) identified during genetic screening of BRCA1/2 genes may affect splicing. We assessed at RNA level the effect of four BRCA1 and ten BRCA2 UVs with a putative splice effect, as predicted in silico. The variants selected for this study were beyond the positions -1, -2 or +1, +2 from the exon, and were not previously described (n = 8) or their effect on splicing was not assessed previously (n = 6). Lymphocytes from UV carriers and healthy controls were cultured and treated with puromycin to prevent nonsense-mediated mRNA decay. The relative contribution of each allele to the various transcripts was assessed using combinations of allele-specific and transcript-specific primers. BRCA2 c.425G>T, c.7976+3_7976+4del and c.8754+3G>C give rise to aberrant transcripts BRCA2Δ4, BRCA2Δ17 and retention of 46nt of intron 21, respectively, and were considered pathogenic. BRCA1 c.4987-3C>G gives rise to BRCA1Δ17 that is likely pathogenic; however, residual expression of the full-length transcript from the variant allele could not be excluded. BRCA1 c.692C>T, c.693G>A and BRCA2 c.6935A>T, besides expressing the full-length transcript, increased expression of BRCA1Δ11 and BRCA2Δ12, respectively. As these are naturally occurring isoforms, also observed in controls, the clinical relevance is unclear. The seven remaining UVs did not affect splicing and three intronic variants were therefore classified as neutral. In conclusion, the RNA analysis results clarified the clinical relevance of 6 of the 14 studied UVs and thereby greatly improve the genetic counselling of high-risk breast/ovarian cancer patients carrying these classified variants.

Unclassified variants in BRCA genes: guidelines for interpretation

In the last few years, several studies have focused on the interpretation of unclassified variants (UVs) of BRCA1 and BRCA2 genes. Analysis of UVs through a unique approach is not sufficient to understand their role in the development of tumors. Thus, it is clear that assembling results from different sources (genetic and epidemiological data, histopathological features, and in vitro and in silico analyses) represents a powerful way to classify such variants. Building reliable integrated models for UV classification requires the joining of many working groups to collaborative consortia, allowing data exchange and improvements of methods. This will lead to improvement in the predictivity of gene testing in BRCA1 and BRCA2 and, consequently, to an increase in the number of families that can be correctly classified as linked or unlinked to these genes, allowing more accurate genetic counseling and clinical management.

Prediction of single-nucleotide substitutions that result in exon skipping: identification of a splicing silencer in BRCA1 exon 6

Missense, nonsense, and translationally silent mutations can inactivate genes by altering the inclusion of mutant exons in mRNA, but their overall frequency among disease-causing exonic substitutions is unknown. Here, we have tested missense and silent mutations deposited in the BRCA1 mutation databases of unclassified variants for their effects on exon inclusion. Analysis of 21 BRCA1 variants using minigene assays revealed a single exon-skipping mutation c.231G>T. Comprehensive mutagenesis of an adjacent 12-nt segment showed that this silent mutation resulted in a higher level of exon skipping than the 35 other single-nucleotide substitutions. Exon inclusion levels of mutant constructs correlated significantly with predicted splicing enhancers/silencers, prompting the development of two online utilities freely available at http://www.dbass.org.uk. EX-SKIP quickly estimates which allele is more susceptible to exon skipping, whereas HOT-SKIP examines all possible mutations at each exon position and identifies candidate exon-skipping positions/substitutions. We demonstrate that the distribution of exon-skipping and disease-associated substitutions previously identified in coding regions was biased toward top-ranking HOT-SKIP mutations. Finally, we show that proteins 9G8, SC35, SF2/ASF, Tra2, and hnRNP A1 were associated with significant alterations of BRCA1 exon 6 inclusion in the mRNA. Together, these results facilitate prediction of exonic substitutions that reduce exon inclusion in mature transcripts.

Ex vivo splicing assays of mutations at noncanonical positions of splice sites in USHER genes

Molecular diagnosis in Usher syndrome type 1 and 2 patients led to the identification of 21 sequence variations located in noncanonical positions of splice sites in MYO7A, CDH23, USH1C, and USH2A genes. To establish experimentally the splicing pattern of these substitutions, whose impact on splicing is not always predictable by available softwares, ex vivo splicing assays were performed. The branch-point mapping strategy was also used to investigate further a putative branch-point mutation in USH2A intron 43. Aberrant splicing was demonstrated for 16 of the 21 (76.2%) tested sequence variations. The mutations resulted more frequently in activation of a nearby cryptic splice site or use of a de novo splice site than exon skipping (37.5%). This study allowed the reclassification as splicing mutations of one silent (c.7872G>A (p.Glu2624Glu) in CDH23) and four missense mutations (c.2993G>A (p.Arg998Lys) in USH2A, c.592G>A (p.Ala198Thr), c.3503G>C [p.Arg1168Pro], c.5944G>A (p.Gly1982Arg) in MYO7A), whereas it provided clues about a role in structure/function in four other cases: c.802G>A (p.Gly268Arg), c.653T>A (p.Val218Glu) (USH2A), and c.397C>T (p.His133Tyr), c.3502C>T (p.Arg1168Trp) (MYO7A). Our data provide insights into the contribution of splicing mutations in Usher genes and illustrate the need to define accurately their splicing outcome for diagnostic purposes.

Effect of BRCA2 sequence variants predicted to disrupt exonic splice enhancers on BRCA2 transcripts

Background

Genetic screening of breast cancer patients and their families have identified a number of variants of unknown clinical significance in the breast cancer susceptibility genes, BRCA1 and BRCA2. Evaluation of such unclassified variants may be assisted by web-based bioinformatic prediction tools, although accurate prediction of aberrant splicing by unclassified variants affecting exonic splice enhancers (ESEs) remains a challenge.

Methods

This study used a combination of RT-PCR analysis and splicing reporter minigene assays to assess five unclassified variants in the BRCA2 gene that we had previously predicted to disrupt an ESE using bioinformatic approaches.

Results

Analysis of BRCA2 c.8308 G > A (p.Ala2770Thr) by mRNA analysis, and BRCA2 c.8962A > G (p.Ser2988Gly), BRCA2 c.8972G > A (p.Arg2991His), BRCA2 c.9172A > G (p.Ser3058Gly), and BRCA2 c.9213G > T (p.Glu3071Asp) by a minigene assay, revealed no evidence for aberrant splicing.

Conclusions

These results illustrate the need for improved methods for predicting functional ESEs and the potential consequences of sequence variants contained therein.

Clinical relevance of rare germline sequence variants in cancer genes: evolution and application of classification models

Multifactorial models developed for BRCA1/2 variant classification have proved very useful for delineating BRCA1/2 variants associated with very high risk of cancer, or with little clinical significance. Recent linkage of this quantitative assessment of risk to clinical management guidelines has provided a basis to standardize variant reporting, variant classification and management of families with such variants, and can theoretically be applied to any disease gene. As proof of principle, the multifactorial approach already shows great promise for application to the evaluation of mismatch repair gene variants identified in families with suspected Lynch syndrome. However there is need to be cautious of the noted limitations and caveats of the current model, some of which may be exacerbated by differences in ascertainment and biological pathways to disease for different cancer syndromes.

A high proportion of DNA variants of BRCA1 and BRCA2 is associated with aberrant splicing in breast/ovarian cancer patients

PURPOSE

Most BRCA1/2 mutations are of unknown clinical relevance. An increasing amount of evidence indicates that there can be deleterious effects through the disruption of the splicing process. We have investigated the effect of aberrant splicing of BRCA1/2 on hereditary breast/ovarian cancer (HBOC).

EXPERIMENTAL DESIGN

DNA variants were analyzed with splicing prediction programs to select putative splicing mutations. Splicing assays of 57 genetic variants were done by lymphocyte reverse transcription-PCR and/or hybrid minigenes in HeLa and nontumor breast epithelial cells.

RESULTS

Twenty-four BRCA1/2 variants of Spanish HBOC patients were bioinformatically preselected. Functional assays showed that 12 variants induced anomalous splicing patterns, 6 of which accounted for 58.5% of BRCA1 families. To further evaluate the defective splicing of BRCA1/2, we analyzed 31 Breast Cancer Information Core Database (BIC) and two artificial variants that were generated by mutagenesis. Sixteen variants induced different degrees of aberrant splicing. Altogether, anomalous splicing was caused by 28 BRCA1/2 variants of all types, indicating that any DNA change can disrupt pre-mRNA processing. We show that a wide range of regulatory elements can be involved, including the canonical and cryptic splice sites, the polypyrimidine tract, and splicing enhancers/silencers. Twenty mutations were predicted to truncate the BRCA proteins and/or to delete essential domains, thus supporting a role in HBOC.

CONCLUSIONS

An important fraction of DNA variants of BRCA1/2 presents splicing aberrations that may represent a relevant disease-causing mechanism in HBOC. The identification of splicing disruptions by functional assays is a valuable tool to discriminate between benign polymorphisms and pathogenic mutations.