Targeted RNA-seq successfully identifies normal and pathogenic splicing events in breast/ovarian cancer susceptibility and Lynch syndrome genes.

A subset of genetic variants found through screening of patients with hereditary breast and ovarian cancer syndrome (HBOC) and Lynch syndrome impact RNA splicing. Through target enrichment of the transcriptome, it is possible to perform deep-sequencing and to identify the different and even rare mRNA isoforms. A targeted RNA-seq approach was used to analyse the naturally-occurring splicing events for a panel of 8 breast and/or ovarian cancer susceptibility genes (BRCA1, BRCA2, RAD51C, RAD51D, PTEN, STK11, CDH1, TP53), 3 Lynch syndrome genes (MLH1, MSH2, MSH6) and the fanconi anaemia SLX4 gene, in which monoallelic mutations were found in non-BRCA families. For BRCA1, BRCA2, RAD51C and RAD51D the results were validated by capillary electrophoresis and were compared to a non-targeted RNA-seq approach. We also compared splicing events from lymphoblastoid cell-lines with those from breast and ovarian fimbriae tissues. The potential of targeted RNA-seq to detect pathogenic changes in RNA-splicing was validated by the inclusion of samples with previously well characterized BRCA1/2 genetic variants. In our study, we update the catalogue of normal splicing events for BRCA1/2, provide an extensive catalogue of normal RAD51C and RAD51D alternative splicing, and list splicing events found for eight other genes. Additionally, we show that our approach allowed the identification of aberrant splicing events due to the presence of BRCA1/2 genetic variants and distinguished between complete and partial splicing events. In conclusion, targeted-RNA-seq can be very useful to classify variants based on their putative pathogenic impact on splicing.

Correlated evolutionary rates across genomic compartments in Annonaceae.

The molecular clock hypothesis is an important concept in biology. Deviations from a constant rate of nucleotide substitution have been found widely among lineages, genomes, genes and individual sites. Phylogenetic research can accommodate for these differences in applying specific models of evolution. Lineage-specific rate heterogeneity however can generate bi- or multimodal distributions of substitution rates across the branches of a tree and this may mislead phylogenetic inferences with currently available models. The plant family Annonaceae is an excellent case to study lineage-specific rate heterogeneity. The two major sister subfamilies, Annonoideae and Malmeoideae, have shown great discrepancies in branch lengths. We used high-throughput sequencing data of 72 genes, 99 spacers and 16 introns from 24 chloroplast genomes and nuclear ribosomal DNA of 23 species to study the molecular rate of evolution in Annonaceae. In all analyses, longer branch lengths and/or higher substitution rates were found for the Annonoideae compared to the Malmeoideae. The Annonaceae had wide variability in chloroplast length, ranging from minimal 175,684bp to 201,723 for Annonoideae and minimal 152,357 to 170,985bp in Malmeoideae, mostly reflecting variation in inverted-repeat length. The Annonoideae showed a higher GC-content in the conserved parts of the chloroplast genome and higher omega (dN/dS)-ratios than the Malmeoideae, which could indicate less stringent purifying selection, a pattern that has been found in groups with small population sizes. This study generates new insights into the processes causing lineage-specific rate heterogeneity, which could lead to improved phylogenetic methods.

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.

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.

Comprehensive annotation of splice junctions supports pervasive alternative splicing at the BRCA1 locus: a report from the ENIGMA consortium.

Loss-of-function germline mutations in BRCA1 (MIM #113705) confer markedly increased risk of breast and ovarian cancer. The full-length transcript codifies for a protein involved in DNA repair pathways and cell-cycle checkpoints. Several BRCA1 splicing isoforms have been described in public domain databases, but the physiological role (if any) of BRCA1 alternative splicing remains to be established. An accurate description of 'naturally occurring' alternative splicing at this locus is a prerequisite to understand its biological significance. However, a systematic analysis of alternative splicing at the BRCA1 locus is yet to be conducted. Here, the Evidence-Based Network for the Interpretation of Germ-Line Mutant Alleles consortium combines RT-PCR, exon scanning, cloning, sequencing and relative semi-quantification to describe naturally occurring BRCA1 alternative splicing with unprecedented resolution. The study has been conducted in blood-related RNA sources, commonly used for clinical splicing assays, as well as in one healthy breast tissue. We have characterized a total of 63 BRCA1 alternative splicing events, including 35 novel findings. A minimum of 10 splicing events (Δ1Aq, Δ5, Δ5q, Δ8p, Δ9, Δ(9,10), Δ9_11, Δ11q, Δ13p and Δ14p) represent a substantial fraction of the full-length expression level (ranging from 5 to 100%). Remarkably, our data indicate that BRCA1 alternative splicing is similar in blood and breast, a finding supporting the clinical relevance of blood-based in vitro splicing assays. Overall, our data suggest an alternative splicing model in which most non-mutually exclusive alternative splicing events are randomly combined into individual mRNA molecules to produce hundreds of different BRCA1 isoforms.

Emerging importance of chemokine receptor CXCR3 and its ligands in cardiovascular diseases.

The CXC chemokines, CXCL4, -9, -10, -11, CXCL4L1, and the CC chemokine CCL21, activate CXC chemokine receptor 3 (CXCR3), a cell-surface G protein-coupled receptor expressed mainly by Th1 cells, cytotoxic T (Tc) cells and NK cells that have a key role in immunity and inflammation. However, CXCR3 is also expressed by vascular smooth muscle and endothelial cells, and appears to be important in controlling physiological vascular function. In the last decade, evidence from pre-clinical and clinical studies has revealed the participation of CXCR3 and its ligands in multiple cardiovascular diseases (CVDs) of different aetiologies including atherosclerosis, hypertension, cardiac hypertrophy and heart failure, as well as in heart transplant rejection and transplant coronary artery disease (CAD). CXCR3 ligands have also proven to be valid biomarkers for the development of heart failure and left ventricular dysfunction, suggesting an underlining pathophysiological relation between levels of these chemokines and the development of adverse cardiac remodelling. The observation that several of the above-mentioned chemokines exert biological actions independent of CXCR3 provides both opportunities and challenges for developing effective drug strategies. In this review, we provide evidence to support our contention that CXCR3 and its ligands actively participate in the development and progression of CVDs, and may additionally have utility as diagnostic and prognostic biomarkers.

A randomised controlled phase II trial of pre-operative celecoxib treatment reveals anti-tumour transcriptional response in primary breast cancer.

INTRODUCTION: Cyclooxygenase-2 (COX-2) is frequently over-expressed in primary breast cancer. In transgenic breast cancer models, over-expression of COX-2 leads to tumour formation while COX-2 inhibition exerts anti-tumour effects in breast cancer cell lines. To further determine the effect of COX-2 inhibition in primary breast cancer, we aimed to identify transcriptional changes in breast cancer tissues of patients treated with the selective COX-2 inhibitor celecoxib. METHODS: In a single-centre double-blind phase II study, thirty-seven breast cancer patients were randomised to receive either pre-operative celecoxib (400 mg) twice daily for two to three weeks (n = 22) or a placebo according to the same schedule (n = 15). Gene expression in fresh-frozen pre-surgical biopsies (before treatment) and surgical excision specimens (after treatment) was profiled by using Affymetrix arrays. Differentially expressed genes and altered pathways were bioinformatically identified. Expression of selected genes was validated by quantitative PCR (qPCR). Immunohistochemical protein expression analyses of the proliferation marker Ki-67, the apoptosis marker cleaved caspase-3 and the neo-angiogenesis marker CD34 served to evaluate biological response. RESULTS: We identified 972 and 586 significantly up- and down-regulated genes, respectively, in celecoxib-treated specimens. Significant expression changes in six out of eight genes could be validated by qPCR. Pathway analyses revealed over-representation of deregulated genes in the networks of proliferation, cell cycle, extracellular matrix biology, and inflammatory immune response. The Ki-67 mean change relative to baseline was -29.1% (P = 0.019) and -8.2% (P = 0.384) in the treatment and control arm, respectively. Between treatment groups, the change in Ki-67 was statistically significant (P = 0.029). Cleaved caspase-3 and CD34 expression were not significantly different between the celecoxib-treated and placebo-treated groups. CONCLUSIONS: Short-term COX-2 inhibition by celecoxib induces transcriptional programs supporting anti-tumour activity in primary breast cancer tissue. The impact on proliferation-associated genes is reflected by a reduction of Ki-67 positive cells. Therefore, COX-2 inhibition should be considered as a treatment strategy for further clinical testing in primary breast cancer. TRIAL REGISTRATION: ClinicalTrials.gov NCT01695226.

A guide for functional analysis of BRCA1 variants of uncertain significance.

Germline mutations in the tumor suppressor gene BRCA1 confer an estimated lifetime risk of 56-80% for breast cancer and 15-60% for ovarian cancer. Since the mid 1990s when BRCA1 was identified, genetic testing has revealed over 1,500 unique germline variants. However, for a significant number of these variants, the effect on protein function is unknown making it difficult to infer the consequences on risks of breast and ovarian cancers. Thus, many individuals undergoing genetic testing for BRCA1 mutations receive test results reporting a variant of uncertain clinical significance (VUS), leading to issues in risk assessment, counseling, and preventive care. Here, we describe functional assays for BRCA1 to directly or indirectly assess the impact of a variant on protein conformation or function and how these results can be used to complement genetic data to classify a VUS as to its clinical significance. Importantly, these methods may provide a framework for genome-wide pathogenicity assignment.

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.

Phylogenomics of the Major Tropical Plant Family Annonaceae Using Targeted Enrichment of Nuclear Genes.

Targeted enrichment and sequencing of hundreds of nuclear loci for phylogenetic reconstruction is becoming an important tool for plant systematics and evolution. Annonaceae is a major pantropical plant family with 110 genera and ca. 2,450 species, occurring across all major and minor tropical forests of the world. Baits were designed by sequencing the transcriptomes of five species from two of the largest Annonaceae subfamilies. Orthologous loci were identified. The resulting baiting kit was used to reconstruct phylogenetic relationships at two different levels using concatenated and gene tree approaches: a family wide Annonaceae analysis sampling 65 genera and a species level analysis of tribe Piptostigmateae sampling 29 species with multiple individuals per species. DNA extraction was undertaken mainly on silicagel dried leaves, with two samples from herbarium dried leaves. Our kit targets 469 exons (364,653 bp of sequence data), successfully capturing sequences from across Annonaceae. Silicagel dried and herbarium DNA worked equally well. We present for the first time a nuclear gene-based phylogenetic tree at the generic level based on 317 supercontigs. Results mainly confirm previous chloroplast based studies. However, several new relationships are found and discussed. We show significant differences in branch lengths between the two large subfamilies Annonoideae and Malmeoideae. A new tribe, Annickieae, is erected containing a single African genus Annickia. We also reconstructed a well-resolved species-level phylogenetic tree of the Piptostigmteae tribe. Our baiting kit is useful for reconstructing well-supported phylogenetic relationships within Annonaceae at different taxonomic levels. The nuclear genome is mainly concordant with plastome information with a few exceptions. Moreover, we find that substitution rate heterogeneity between the two subfamilies is also found within the nuclear compartment, and not just plastomes and ribosomal DNA as previously shown. Our results have implications for understanding the biogeography, molecular dating and evolution of Annonaceae.

Full in-frame exon 3 skipping of BRCA2 confers high risk of breast and/or ovarian cancer.

Germline pathogenic variants in the BRCA2 gene are associated with a cumulative high risk of breast/ovarian cancer. Several BRCA2 variants result in complete loss of the exon-3 at the transcript level. The pathogenicity of these variants and the functional impact of loss of exon 3 have yet to be established. As a collaboration of the COVAR clinical trial group (France), and the ENIGMA consortium for investigating breast cancer gene variants, this study evaluated 8 BRCA2 variants resulting in complete deletion of exon 3. Clinical information for 39 families was gathered from Portugal, France, Denmark and Sweden. Multifactorial likelihood analyses were conducted using information from 293 patients, for 7 out of the 8 variants (including 6 intronic). For all variants combined the likelihood ratio in favor of causality was 4.39*1025. These results provide convincing evidence for the pathogenicity of all examined variants that lead to a total exon 3 skipping, and suggest that other variants that result in complete loss of exon 3 at the molecular level could be associated with a high risk of cancer comparable to that associated with classical pathogenic variants in BRCA1 or BRCA2 gene. In addition, our functional study shows, for the first time, that deletion of exon 3 impairs the ability of cells to survive upon Mitomycin-C treatment, supporting lack of function for the altered BRCA2 protein in these cells. Finally, this study demonstrates that any variant leading to expression of only BRCA2 delta-exon 3 will be associated with an increased risk of breast and ovarian cancer.

Familial breast/ovarian cancer and BRCA1/2 genetic screening: the role of immunohistochemistry as an additional method in the selection of patients.

Only 20-25% of families screened for BRCA1/2 mutations are found positive. Because only a positive result is informative, we studied the role of BRCA1/2 immunohistochemistry as an additional method for patient selection. From 53 high-risk-affected probands, 18 (34%) had available paraffin blocks of their tumors and were selected for this study. Mutation screening was done by conformation-sensitive gel electrophoresis and multiplex ligation-dependent probe amplification. For immunohistochemistry, 21 neoplastic specimens (15 breast carcinomas, 5 ovary neoplasms, and 1 rectal adenocarcinoma) were analyzed with BRCA1 (monoclonal antibody, Ab-1, oncogene) and BRCA2 (polyclonal antibody, Ab-2, oncogene) antibodies. Absence of the BRCA1 protein was confirmed in negative tumors by Western blotting. Seven patients were positive for BRCA1/2 mutations: 5 for BRCA1 and 2 for BRCA2. Four out of five positive patients had tumors negative for BRCA1 immunostaining, and the remaining 13 BRCA1-negative patients had positive BRCA1 immunostaining in all tumor samples. Sensitivity to predict for BRCA1 mutation carriers was 80%, and specificity was 100%, with a positive predictive value of 100% and a negative predictive value of 93%. This correlation was statistically significant (p=0.001). No correlation was observed for BRCA2. If larger studies confirm these results, high-risk patients with BRCA1-negative tumors should be screened first for this gene.

Screening for a BRCA2 rearrangement in high-risk breast/ovarian cancer families: evidence for a founder effect and analysis of the associated phenotypes.

PURPOSE: BRCA2 rearrangements are rare genetic events. A large BRCA2 genomic insertion was recurrently observed in our participants, and we sought to characterize it at the molecular and phenotypic level. PATIENTS AND METHODS: We studied 210 high-risk breast/ovarian cancer families. Fifty-three probands were fully screened for BRCA1/2 mutations, and three of 53 had a large insertion in exon 3 of BRCA2. This finding was analyzed by polymerase chain reaction (PCR), reverse transcriptase PCR (RT-PCR), and sequencing. An additional 157 consecutive families were screened for this mutation by a three-step PCR method. Phenotype and haplotype analysis was also performed. RESULTS: Sixteen BRCA mutations were observed in 19 of 53 patients (36% detection rate). A recurrent Alu motif insertion in position c.156_157 was observed after sequencing of an abnormal fragment obtained after the amplification of BRCA2 exon 3. RT-PCR revealed exon 3 skipping. Screening of this rearrangement identified 14 additional families (out of 157). In total, 17 (8%) of 210 high-risk families ascertained in our clinic were positive for this mutation. Segregation of a common haplotype (from D13S260 to D13S1695) confirmed a common origin, estimated to have occurred 2,400 to 2,600 years ago. The following four cancer phenotypes were observed in the 17 positive families: female breast (n = 9), male breast (n = 4), breast/ovarian (n = 2), and heterogeneous (n = 2). Male breast cancer was more frequently observed in c.156_157insAlu-positive families compared with negative families (23% v 12%, respectively), and 33% of all male breast cancer families with an identified BRCA mutation were c.156_157insAlu positive. CONCLUSION: c.156_157insAlu is a founder mutation of Portuguese origin and is the most frequent BRCA2 rearrangement described to date.