Comprehensive prediction of mRNA splicing effects of BRCA1 and BRCA2 variants

Solid Pseudopapillary Neoplasm of the Pancreas and Abdominal Desmoid Tumor in a Patient Carrying Two Different BRCA2 Germline Mutations: New Horizons from Tumor Molecular Profiling

This case report describes the history of a 41 year-old woman with a solid pseudopapillary neoplasm (SPN) of the pancreas and a metachronous abdominal desmoid tumor (DT) that occurred two years after the SPN surgical resection. At next-generation sequencing of 174 cancer-related genes, both neoplasms harbored a CTNNB1 somatic mutation which was different in each tumor. Moreover, two BRCA2 pathogenic mutations were found in both tumors, confirmed as germline by the sequencing of normal tissue. The BRCA2 mutations were c.631G>A, resulting in the amino-acid change p.V211I, and c.7008-2A>T, causing a splice acceptor site loss. However, as the two neoplasms showed neither loss of heterozygosity nor somatic mutation in the second BRCA2 allele, they cannot be considered as BRCA-dependent tumors. Nevertheless, this study highlights the important opportunities opened by extensive tumor molecular profiling. In this particular case, it permitted the detection of BRCA2-germline mutations, essential for addressing the necessary BRCA-related genetic counseling, surveillance, and screening for the patient and her family.

Optic Atrophy and Inner Retinal Thinning in CACNA1F-related Congenital Stationary Night Blindness

Hemizygous pathogenic variants in CACNA1F lead to defective signal transmission from retinal photoreceptors to bipolar cells and cause incomplete congenital stationary night blindness in humans. Although the primary defect is at the terminal end of first-order neurons (photoreceptors), there is limited knowledge of higher-order neuronal changes (inner retinal) in this disorder. This study aimed to investigate inner retinal changes in CACNA1F-retinopathy by analyzing macular ganglion cell layer-inner plexiform layer (GCL-IPL) thickness and optic disc pallor in 22 subjects with molecularly confirmed CACNA1F-retinopathy. Detailed ocular phenotypic data including distance and color vision, refraction and electroretinogram (ERG) were collected. Distance vision was universally reduced (mean: 0.42 LogMAR), six had abnormal color vision and myopia was common (n = 15; mean: −6.32 diopters). Mean GCL-IPL thickness was significantly lower in patients (55.00 µm) compared to age-matched controls (n = 87; 84.57 µm; p << 0.001). The GCL-IPL thickness correlated with scotopic standard (p = 0.04) and bright-flash (p = 0.014) ERG b/a ratios and photopic b-wave amplitudes (p = 0.05). Twenty-one patients had some degree of disc pallor (bilateral in 19). Fifteen putative disease-causing, including five novel variants were identified. This study establishes macular inner retinal thinning and optic atrophy as characteristic features of CACNA1F-retinopathy, which are independent of myopia and could impact potential future treatment strategies.

A proposed molecular mechanism for pathogenesis of severe RNA-viral pulmonary infections

Background: Certain riboviruses can cause severe pulmonary complications leading to death in some infected patients. We propose that DNA damage induced-apoptosis accelerates viral release, triggered by depletion of host RNA binding proteins (RBPs) from nuclear RNA bound to replicating viral sequences. Methods: Information theory-based analysis of interactions between RBPs and individual sequences in the Severe Acute Respiratory Syndrome CoronaVirus 2 (SARS-CoV-2), Influenza A (H3N1), HIV-1, and Dengue genomes identifies strong RBP binding sites in these viral genomes. Replication and expression of viral sequences is expected to increasingly sequester RBPs - SRSF1 and RNPS1. Ordinarily, RBPs bound to nascent host transcripts prevents their annealing to complementary DNA. Their depletion induces destabilizing R-loops. Chromosomal breakage occurs when an excess of unresolved R-loops collide with incoming replication forks, overwhelming the DNA repair machinery. We estimated stoichiometry of inhibition of RBPs in host nuclear RNA by counting competing binding sites in replicating viral genomes and host RNA. Results: Host RBP binding sites are frequent and conserved among different strains of RNA viral genomes. Similar binding motifs of SRSF1 and RNPS1 explain why DNA damage resulting from SRSF1 depletion is complemented by expression of RNPS1. Clustering of strong RBP binding sites coincides with the distribution of RNA-DNA hybridization sites across the genome. SARS-CoV-2 replication is estimated to require 32.5-41.8 hours to effectively compete for binding of an equal proportion of SRSF1 binding sites in host encoded nuclear RNAs. Significant changes in expression of transcripts encoding DNA repair and apoptotic proteins were found in an analysis of influenza A and Dengue-infected cells in some individuals. Conclusions: R-loop-induced apoptosis indirectly resulting from viral replication could release significant quantities of membrane-associated virions into neighboring alveoli. These could infect adjacent pneumocytes and other tissues, rapidly compromising lung function, causing multiorgan system failure and other described symptoms.

A proposed molecular mechanism for pathogenesis of severe RNA-viral pulmonary infections

Background: Certain riboviruses can cause severe pulmonary complications leading to death in some infected patients. We propose that DNA damage induced-apoptosis accelerates viral release, triggered by depletion of host RNA binding proteins (RBPs) from nuclear RNA bound to replicating viral sequences. Methods: Information theory-based analysis of interactions between RBPs and individual sequences in the Severe Acute Respiratory Syndrome CoronaVirus 2 (SARS-CoV-2), Influenza A (H3N2), HIV-1, and Dengue genomes identifies strong RBP binding sites in these viral genomes. Replication and expression of viral sequences is expected to increasingly sequester RBPs - SRSF1 and RNPS1. Ordinarily, RBPs bound to nascent host transcripts prevents their annealing to complementary DNA. Their depletion induces destabilizing R-loops. Chromosomal breakage occurs when an excess of unresolved R-loops collide with incoming replication forks, overwhelming the DNA repair machinery. We estimated stoichiometry of inhibition of RBPs in host nuclear RNA by counting competing binding sites in replicating viral genomes and host RNA. Results: Host RBP binding sites are frequent and conserved among different strains of RNA viral genomes. Similar binding motifs of SRSF1 and RNPS1 explain why DNA damage resulting from SRSF1 depletion is complemented by expression of RNPS1. Clustering of strong RBP binding sites coincides with the distribution of RNA-DNA hybridization sites across the genome. SARS-CoV-2 replication is estimated to require 32.5-41.8 hours to effectively compete for binding of an equal proportion of SRSF1 binding sites in host encoded nuclear RNAs. Significant changes in expression of transcripts encoding DNA repair and apoptotic proteins were found in an analysis of influenza A and Dengue-infected cells in some individuals. Conclusions: R-loop-induced apoptosis indirectly resulting from viral replication could release significant quantities of membrane-associated virions into neighboring alveoli. These could infect adjacent pneumocytes and other tissues, rapidly compromising lung function, causing multiorgan system failure and other described symptoms.

Pan-cancer repository of validated natural and cryptic mRNA splicing mutations

We present a major public resource of mRNA splicing mutations validated according to multiple lines of evidence of abnormal gene expression. Likely mutations present in all tumor types reported in the Cancer Genome Atlas (TCGA) and the International Cancer Genome Consortium (ICGC) were identified based on the comparative strengths of splice sites in tumor versus normal genomes, and then validated by respectively comparing counts of splice junction spanning and abundance of transcript reads in RNA-Seq data from matched tissues and tumors lacking these mutations. The comprehensive resource features 341,486 of these validated mutations, the majority of which (69.9%) are not present in the Single Nucleotide Polymorphism Database (dbSNP 150). There are 131,347 unique mutations which weaken or abolish natural splice sites, and 222,071 mutations which strengthen cryptic splice sites (11,932 affect both simultaneously). 28,812 novel or rare flagged variants (with <1% population frequency in dbSNP) were observed in multiple tumor tissue types. An algorithm was developed to classify variants into splicing molecular phenotypes that integrates germline heterozygosity, degree of information change and impact on expression. The classification thresholds were calibrated against the ClinVar clinical database phenotypic assignments. Variants are partitioned into allele-specific alternative splicing, likely aberrant and aberrant splicing phenotypes. Single variants or chromosome ranges can be queried using a Global Alliance for Genomics and Health (GA4GH)-compliant, web-based Beacon "Validated Splicing Mutations" either separately or in aggregate alongside other Beacons through the public Beacon Network, as well as through our website. The website provides additional information, such as a visual representation of supporting RNAseq results, gene expression in the corresponding normal tissues, and splicing molecular phenotypes.

Novel mutations in actionable breast cancer genes by targeted sequencing in an ethnically homogenous cohort

Background

Genetic testing is becoming an essential tool for breast cancer (BC) diagnosis and treatment pathway, and particularly important for early detection and cancer prevention. The purpose of this study was to explore the diagnostic yield of targeted sequencing of the high priority BC genes.

Methods

We have utilized a cost-effective targeted sequencing approach of high priority actionable BC genes (BRCA1, BRCA2, ERBB2 and TP53) in a homogeneous patient cohort from Bangladesh (n = 52) by using tumor and blood samples.

Results

Blood derived targeted sequencing revealed 25.58% (11/43) clinically relevant mutations (both pathogenic and variants of uncertain significance (VUS)), with 13.95% (6/43) of samples carrying a pathogenic mutations. We have identified and validated five novel pathogenic germline mutations in this cohort, comprising of two frameshift deletions in BRCA2, and missense mutations in BRCA1, BRCA2 and ERBB2 gene respectively. Furthermore, we have identified three pathogenic mutations and a VUS within three tumor samples, including a sample carrying pathogenic mutations impacting both TP53 (c.322dupG; a novel frameshift insertion) and BRCA1 genes (c.116G > A). 22% of tissue samples had a clinically relevant TP53 mutation. Although the cohort is small, we have found pathogenic mutations to be enriched in BRCA2 (9.30%, 4/43) compare to BRCA1 (4.65%, 2/43). The frequency of germline VUS mutations found to be similar in both BRCA1 (4.65%; 2/43) and BRCA2 (4.65%; 2/43) compared to ERBB2 (2.32%; 1/43).

Conclusions

This is the first genetic study of BC predisposition genes in this population, implies that genetic screening through targeted sequencing can detect clinically significant and actionable BC-relevant mutations.

Electronic supplementary material

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

Pan-cancer repository of validated natural and cryptic mRNA splicing mutations

We present a major public resource of mRNA splicing mutations validated according to multiple lines of evidence of abnormal gene expression. Likely mutations present in all tumor types reported in the Cancer Genome Atlas (TCGA) and the International Cancer Genome Consortium (ICGC) were identified based on the comparative strengths of splice sites in tumor versus normal genomes, and then validated by respectively comparing counts of splice junction spanning and abundance of transcript reads in RNA-Seq data from matched tissues and tumors lacking these mutations. The comprehensive resource features 341,486 of these validated mutations, the majority of which (69.9%) are not present in the Single Nucleotide Polymorphism Database (dbSNP 150). There are 131,347 unique mutations which weaken or abolish natural splice sites, and 222,071 mutations which strengthen cryptic splice sites (11,932 affect both simultaneously). 28,812 novel or rare flagged variants (with <1% population frequency in dbSNP) were observed in multiple tumor tissue types. An algorithm was developed to classify variants into splicing molecular phenotypes that integrates germline heterozygosity, degree of information change and impact on expression. The classification thresholds were calibrated against the ClinVar clinical database phenotypic assignments. Variants are partitioned into allele-specific alternative splicing, likely aberrant and aberrant splicing phenotypes. Single variants or chromosome ranges can be queried using a Global Alliance for Genomics and Health (GA4GH)-compliant, web-based Beacon "Validated Splicing Mutations" either separately or in aggregate alongside other Beacons through the public Beacon Network, as well as through our website. The website provides additional information, such as a visual representation of supporting RNAseq results, gene expression in the corresponding normal tissues, and splicing molecular phenotypes.

Preliminary molecular evidence associating a novel BRCA1 synonymous variant with hereditary ovarian cancer syndrome

Extensive molecular screening of the BRCA1/2 (BRCA) genes by massively parallel sequencing (MPS) identified variants of uncertain (or unknown) significance (VUS) and novel variants. We performed a molecular characterization of a novel BRCA1 synonymous variant discovered in a family with hereditary ovarian cancer (HOC) syndrome. We showed that the BRCA1 c.5073 A > T variant might play a pathogenic role in HOC syndrome in this family.

Differential expression of alternatively spliced transcripts related to energy metabolism in colorectal cancer

BACKGROUND

Colorectal cancer (CRC) is one of the most common malignant tumors worldwide. CRC molecular pathogenesis is heterogeneous and may be followed by mutations in oncogenes and tumor suppressor genes, chromosomal and microsatellite instability, alternative splicing alterations, hypermethylation of CpG islands, oxidative stress, impairment of different signaling pathways and energy metabolism. In the present work, we have studied the alterations of alternative splicing patterns of genes related to energy metabolism in CRC.

RESULTS

Using CrossHub software, we analyzed The Cancer Genome Atlas (TCGA) RNA-Seq datasets derived from colon tumor and matched normal tissues. The expression of 1014 alternative mRNA isoforms involved in cell energy metabolism was examined. We found 7 genes with differentially expressed alternative transcripts whereas overall expression of these genes was not significantly altered in CRC. A set of 8 differentially expressed transcripts of interest has been validated by qPCR. These eight isoforms encoded by OGDH, COL6A3, ICAM1, PHPT1, PPP2R5D, SLC29A1, and TRIB3 genes were up-regulated in colorectal tumors, and this is in concordance with the bioinformatics data. The alternative transcript NM_057167 of COL6A3 was also strongly up-regulated in breast, lung, prostate, and kidney tumors. Alternative transcript of SLC29A1 (NM_001078177) was up-regulated only in CRC samples, but not in the other tested tumor types.

CONCLUSIONS

We identified tumor-specific expression of alternative spliced transcripts of seven genes involved in energy metabolism in CRC. Our results bring new knowledge on alternative splicing in colorectal cancer and suggest a set of mRNA isoforms that could be used for cancer diagnosis and development of treatment methods.

A unified analytic framework for prioritization of non-coding variants of uncertain significance in heritable breast and ovarian cancer

BACKGROUND

Sequencing of both healthy and disease singletons yields many novel and low frequency variants of uncertain significance (VUS). Complete gene and genome sequencing by next generation sequencing (NGS) significantly increases the number of VUS detected. While prior studies have emphasized protein coding variants, non-coding sequence variants have also been proven to significantly contribute to high penetrance disorders, such as hereditary breast and ovarian cancer (HBOC). We present a strategy for analyzing different functional classes of non-coding variants based on information theory (IT) and prioritizing patients with large intragenic deletions.

METHODS

We captured and enriched for coding and non-coding variants in genes known to harbor mutations that increase HBOC risk. Custom oligonucleotide baits spanning the complete coding, non-coding, and intergenic regions 10 kb up- and downstream of ATM, BRCA1, BRCA2, CDH1, CHEK2, PALB2, and TP53 were synthesized for solution hybridization enrichment. Unique and divergent repetitive sequences were sequenced in 102 high-risk, anonymized patients without identified mutations in BRCA1/2. Aside from protein coding and copy number changes, IT-based sequence analysis was used to identify and prioritize pathogenic non-coding variants that occurred within sequence elements predicted to be recognized by proteins or protein complexes involved in mRNA splicing, transcription, and untranslated region (UTR) binding and structure. This approach was supplemented by in silico and laboratory analysis of UTR structure.

RESULTS

15,311 unique variants were identified, of which 245 occurred in coding regions. With the unified IT-framework, 132 variants were identified and 87 functionally significant VUS were further prioritized. An intragenic 32.1 kb interval in BRCA2 that was likely hemizygous was detected in one patient. We also identified 4 stop-gain variants and 3 reading-frame altering exonic insertions/deletions (indels).

CONCLUSIONS

We have presented a strategy for complete gene sequence analysis followed by a unified framework for interpreting non-coding variants that may affect gene expression. This approach distills large numbers of variants detected by NGS to a limited set of variants prioritized as potential deleterious changes.

Prioritizing Variants in Complete Hereditary Breast and Ovarian Cancer Genes in Patients Lacking Known BRCA Mutations

BRCA1 and BRCA2 testing for hereditary breast and ovarian cancer (HBOC) does not identify all pathogenic variants. Sequencing of 20 complete genes in HBOC patients with uninformative test results (N = 287), including noncoding and flanking sequences of ATM, BARD1, BRCA1, BRCA2, CDH1, CHEK2, EPCAM, MLH1, MRE11A, MSH2, MSH6, MUTYH, NBN, PALB2, PMS2, PTEN, RAD51B, STK11, TP53, and XRCC2, identified 38,372 unique variants. We apply information theory (IT) to predict and prioritize noncoding variants of uncertain significance in regulatory, coding, and intronic regions based on changes in binding sites in these genes. Besides mRNA splicing, IT provides a common framework to evaluate potential affinity changes in transcription factor (TFBSs), splicing regulatory (SRBSs), and RNA-binding protein (RBBSs) binding sites following mutation. We prioritized variants affecting the strengths of 10 splice sites (four natural, six cryptic), 148 SRBS, 36 TFBS, and 31 RBBS. Three variants were also prioritized based on their predicted effects on mRNA secondary (2°) structure and 17 for pseudoexon activation. Additionally, four frameshift, two in-frame deletions, and five stop-gain mutations were identified. When combined with pedigree information, complete gene sequence analysis can focus attention on a limited set of variants in a wide spectrum of functional mutation types for downstream functional and co-segregation analysis.

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.

Interpretation of mRNA splicing mutations in genetic disease: review of the literature and guidelines for information-theoretical analysis

The interpretation of genomic variants has become one of the paramount challenges in the post-genome sequencing era. In this review we summarize nearly 20 years of research on the applications of information theory (IT) to interpret coding and non-coding mutations that alter mRNA splicing in rare and common diseases. We compile and summarize the spectrum of published variants analyzed by IT, to provide a broad perspective of the distribution of deleterious natural and cryptic splice site variants detected, as well as those affecting splicing regulatory sequences. Results for natural splice site mutations can be interrogated dynamically with Splicing Mutation Calculator, a companion software program that computes changes in information content for any splice site substitution, linked to corresponding publications containing these mutations. The accuracy of IT-based analysis was assessed in the context of experimentally validated mutations. Because splice site information quantifies binding affinity, IT-based analyses can discern the differences between variants that account for the observed reduced (leaky) versus abolished mRNA splicing. We extend this principle by comparing predicted mutations in natural, cryptic, and regulatory splice sites with observed deleterious phenotypic and benign effects. Our analysis of 1727 variants revealed a number of general principles useful for ensuring portability of these analyses and accurate input and interpretation of mutations. We offer guidelines for optimal use of IT software for interpretation of mRNA splicing mutations.

The landscape of alternative splicing in cervical squamous cell carcinoma

BACKGROUND

Alternative splicing (AS) is a key regulatory mechanism in protein synthesis and proteome diversity. In this study, we identified alternative splicing events in four pairs of cervical squamous cell carcinoma (CSCC) and adjacent nontumor tissues using RNA sequencing.

METHODS

The transcripts of the four paired samples were thoroughly analyzed by RNA sequencing. SpliceMap software was used to detect the splicing junctions. Kyoto Encyclopedia of Genes and Genomes pathway analysis was conducted to detect the alternative spliced genes-related signal pathways. The alternative spliced genes were validated by reverse transcription-polymerase chain reaction (RT-PCR).

RESULTS

There were 35 common alternative spliced genes in the four CSCC samples; they were novel and CSCC specific. Sixteen pathways were significantly enriched (P<0.05). One novel 5'AS site in the KLHDC7B gene, encoding kelch domain-containing 7B, and an exon-skipping site in the SYCP2 gene, encoding synaptonemal complex 2, were validated by RT-PCR. The KLHDC7B gene with 5'AS was found in 67.5% (27/40) of CSCC samples and was significantly related with cellular differentiation and tumor size. The exon-skipping site of the SYCP2 gene was found in 35.0% (14/40) of CSCC samples and was significantly related with depth of cervical invasion.

CONCLUSION

The KLHDC7B and the SYCP2 genes with alternative spliced events might be involved in the development and progression of CSCC and could be used as biomarkers in the diagnosis and prognosis of CSCC.

Interpretation of mRNA splicing mutations in genetic disease: review of the literature and guidelines for information-theoretical analysis

The interpretation of genomic variants has become one of the paramount challenges in the post-genome sequencing era. In this review we summarize nearly 20 years of research on the applications of information theory (IT) to interpret coding and non-coding mutations that alter mRNA splicing in rare and common diseases. We compile and summarize the spectrum of published variants analyzed by IT, to provide a broad perspective of the distribution of deleterious natural and cryptic splice site variants detected, as well as those affecting splicing regulatory sequences. Results for natural splice site mutations can be interrogated dynamically with Splicing Mutation Calculator, a companion software program that computes changes in information content for any splice site substitution, linked to corresponding publications containing these mutations. The accuracy of IT-based analysis was assessed in the context of experimentally validated mutations. Because splice site information quantifies binding affinity, IT-based analyses can discern the differences between variants that account for the observed reduced (leaky) versus abolished mRNA splicing. We extend this principle by comparing predicted mutations in natural, cryptic, and regulatory splice sites with observed deleterious phenotypic and benign effects. Our analysis of 1727 variants revealed a number of general principles useful for ensuring portability of these analyses and accurate input and interpretation of mutations. We offer guidelines for optimal use of IT software for interpretation of mRNA splicing mutations.

Splicing mutation analysis reveals previously unrecognized pathways in lymph node-invasive breast cancer

Somatic mutations reported in large-scale breast cancer (BC) sequencing studies primarily consist of protein coding mutations. mRNA splicing mutation analyses have been limited in scope, despite their prevalence in Mendelian genetic disorders. We predicted splicing mutations in 442 BC tumour and matched normal exomes from The Cancer Genome Atlas Consortium (TCGA). These splicing defects were validated by abnormal expression changes in these tumours. Of the 5,206 putative mutations identified, exon skipping, leaky or cryptic splicing was confirmed for 988 variants. Pathway enrichment analysis of the mutated genes revealed mutations in 9 NCAM1-related pathways, which were significantly increased in samples with evidence of lymph node metastasis, but not in lymph node-negative tumours. We suggest that comprehensive reporting of DNA sequencing data should include non-trivial splicing analyses to avoid missing clinically-significant deleterious splicing mutations, which may reveal novel mutated pathways present in genetic disorders.

The alternative splicing side of cancer

Alternative splicing emerges as a potent and pervasive mechanism of gene expression regulation that expands the coding capacity of the genome and forms an intermediate layer of regulation between transcriptional and post-translational networks. Indeed, alternative splicing occupies a pivotal position in developmental programs and in the cell response to external and internal stimuli. Not surprisingly, therefore, its deregulation frequently leads to human disease. In this review we provide an updated overview of the impact of alternative splicing on tumorigenesis. Moreover, we discuss the intricacy of the reciprocal interactions between alternative splicing programs and signal transduction pathways, which appear to be crucially linked to cancer progression in response to the tumor microenvironment. Finally, we focus on the recently described interplay between splicing and chromatin organization which is expected to shed new lights into gene expression regulation in normal and cancer cells.

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.

Combined computational-experimental analyses of CFTR exon strength uncover predictability of exon-skipping level

With the increased number of identified nucleotide sequence variations in genes, the current challenge is to classify them as disease causing or neutral. These variants of unknown clinical significance can alter multiple processes, from gene transcription to RNA splicing or protein function. Using an approach combining several in silico tools, we identified some exons presenting weaker splicing motifs than other exons in the Cystic Fibrosis Transmembrane conductance Regulator (CFTR) gene. These exons exhibit higher rates of basal skipping than exons harboring no identifiable weak splicing signals using minigene assays. We then screened 19 described mutations in three different exons, and identified exon-skipping substitutions. These substitutions induced higher skipping levels in exons having one or more weak splicing motifs. Indeed, this level remained under 2% for exons with strong splicing motifs and could reach 40% for exons having at least one weak motif. Further analysis revealed a functional exon splicing enhancer within exon 3 that was associated with the SR protein SF2/ASF and whose disruption induced exon skipping. Exon skipping was confirmed in vivo in two nasal epithelial cell brushing samples. Our approach, which point out exons with some splicing signals weaknesses, will help spot splicing mutations of clinical relevance.

Interpretation, stratification and evidence for sequence variants affecting mRNA splicing in complete human genome sequences

Information theory-based methods have been shown to be sensitive and specific for predicting and quantifying the effects of non-coding mutations in Mendelian diseases. We present the Shannon pipeline software for genome-scale mutation analysis and provide evidence that the software predicts variants affecting mRNA splicing. Individual information contents (in bits) of reference and variant splice sites are compared and significant differences are annotated and prioritized. The software has been implemented for CLC-Bio Genomics platform. Annotation indicates the context of novel mutations as well as common and rare SNPs with splicing effects. Potential natural and cryptic mRNA splicing variants are identified, and null mutations are distinguished from leaky mutations. Mutations and rare SNPs were predicted in genomes of three cancer cell lines (U2OS, U251 and A431), which were supported by expression analyses. After filtering, tractable numbers of potentially deleterious variants are predicted by the software, suitable for further laboratory investigation. In these cell lines, novel functional variants comprised 6–17 inactivating mutations, 1–5 leaky mutations and 6–13 cryptic splicing mutations. Predicted effects were validated by RNA-seq analysis of the three aforementioned cancer cell lines, and expression microarray analysis of SNPs in HapMap cell lines.

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.

Prediction of mutant mRNA splice isoforms by information theory-based exon definition

Mutations that affect mRNA splicing often produce multiple mRNA isoforms, resulting in complex molecular phenotypes. Definition of an exon and its inclusion in mature mRNA relies on joint recognition of both acceptor and donor splice sites. This study predicts cryptic and exon-skipping isoforms in mRNA produced by splicing mutations from the combined information contents (R(i), which measures binding-site strength, in bits) and distribution of the splice sites defining these exons. The total information content of an exon (R(i),total) is the sum of the R(i) values of its acceptor and donor splice sites, adjusted for the self-information of the distance separating these sites, that is, the gap surprisal. Differences between total information contents of an exon (ΔR(i,total)) are predictive of the relative abundance of these exons in distinct processed mRNAs. Constraints on splice site and exon selection are used to eliminate nonconforming and poorly expressed isoforms. Molecular phenotypes are computed by the Automated Splice Site and Exon Definition Analysis (http://splice.uwo.ca) server. Predictions of splicing mutations were highly concordant (85.2%; n = 61) with published expression data. In silico exon definition analysis will contribute to streamlining assessment of abnormal and normal splice isoforms resulting from mutations.

ENIGMA--evidence-based network for the interpretation of germline mutant alleles: an international initiative to evaluate risk and clinical significance associated with sequence variation in BRCA1 and BRCA2 genes

As genetic testing for predisposition to human diseases has become an increasingly common practice in medicine, the need for clear interpretation of the test results is apparent. However, for many disease genes, including the breast cancer susceptibility genes BRCA1 and BRCA2, a significant fraction of tests results in the detection of a genetic variant for which disease association is not known. The finding of an "unclassified" variant (UV)/variant of uncertain significance (VUS) complicates genetic test reporting and counseling. As these variants are individually rare, a large collaboration of researchers and clinicians will facilitate studies to assess their association with cancer predisposition. It was with this in mind that the ENIGMA consortium (www.enigmaconsortium.org) was initiated in 2009. The membership is both international and interdisciplinary, and currently includes more than 100 research scientists and clinicians from 19 countries. Within ENIGMA, there are presently six working groups focused on the following topics: analysis, clinical, database, functional, tumor histopathology, and mRNA splicing. ENIGMA provides a mechanism to pool resources, exchange methods and data, and coordinately develop and apply algorithms for classification of variants in BRCA1 and BRCA2. It is envisaged that the research and clinical application of models developed by ENIGMA will be relevant to the interpretation of sequence variants in other disease genes.