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

CLN3 transcript complexity revealed by long-read RNA sequencing analysis

BACKGROUND

Batten disease is a group of rare inherited neurodegenerative diseases. Juvenile CLN3 disease is the most prevalent type, and the most common pathogenic variant shared by most patients is the "1-kb" deletion which removes two internal coding exons (7 and 8) in CLN3. Previously, we identified two transcripts in patient fibroblasts homozygous for the 1-kb deletion: the 'major' and 'minor' transcripts. To understand the full variety of disease transcripts and their role in disease pathogenesis, it is necessary to first investigate CLN3 transcription in "healthy" samples without juvenile CLN3 disease.

METHODS

We leveraged PacBio long-read RNA sequencing datasets from ENCODE to investigate the full range of CLN3 transcripts across various tissues and cell types in human control samples. Then we sought to validate their existence using data from different sources.

RESULTS

We found that a readthrough gene affects the quantification and annotation of CLN3. After taking this into account, we detected over 100 novel CLN3 transcripts, with no dominantly expressed CLN3 transcript. The most abundant transcript has median usage of 42.9%. Surprisingly, the known disease-associated 'major' transcripts are detected. Together, they have median usage of 1.5% across 22 samples. Furthermore, we identified 48 CLN3 ORFs, of which 26 are novel. The predominant ORF that encodes the canonical CLN3 protein isoform has median usage of 66.7%, meaning around one-third of CLN3 transcripts encode protein isoforms with different stretches of amino acids. The same ORFs could be found with alternative UTRs. Moreover, we were able to validate the translational potential of certain transcripts using public mass spectrometry data.

CONCLUSION

Overall, these findings provide valuable insights into the complexity of CLN3 transcription, highlighting the importance of studying both canonical and non-canonical CLN3 protein isoforms as well as the regulatory role of UTRs to fully comprehend the regulation and function(s) of CLN3. This knowledge is essential for investigating the impact of the 1-kb deletion and rare pathogenic variants on CLN3 transcription and disease pathogenesis.

Untapped Potential of Poly(ADP-Ribose) Polymerase Inhibitors: Lessons Learned From the Real-World Clinical Homologous Recombination Repair Mutation Testing

Background

Testing for homologous recombination deficiency (HRD) mutations is pivotal to assess individual risk, to proact preventive measures in healthy carriers and to tailor treatments for cancer patients. Increasing prominence of poly(ADP-ribose) polymerase (PARP) inhibitors with remarkable impact on molecular-selected patient survival across diverse nosologies, ingrains testing for BRCA genes and beyond in clinical practice. Nevertheless, testing strategies remain a question of debate. While several pathogenic BRCA1/2 gene variants have been described as founder pathogenic mutations frequently found in patients from Russia, other homologous recombination repair (HRR) genes have not been sufficiently explored. In this study, we present real-world data of routine HRR gene testing in Russia.

Methods

We evaluated clinical and sequencing data from cancer patients who had germline/somatic next-generation sequencing (NGS) HRR gene testing in Russia (BRCA1/2/ATM/CHEK2, or 15 HRR genes). The primary objectives of this study were to evaluate the frequency of BRCA1/2 and non-BRCA gene mutations in real-world unselected patients from Russia, and to determine whether testing beyond BRCA1/2 is feasible.

Results

Data of 2,032 patients were collected from February 2021 to February 2023. Most had breast (n = 715, 35.2%), ovarian (n = 259, 12.7%), pancreatic (n = 85, 4.2%), or prostate cancer (n = 58, 2.9%). We observed 586 variants of uncertain significance (VUS) and 372 deleterious variants (DVs) across 487 patients, with 17.6% HRR-mutation positivity. HRR testing identified 120 (11.8%) BRCA1/2-positive, and 172 (16.9%) HRR-positive patients. With 51 DVs identified in 242 formalin-fixed paraffin-embedded (FFPE), testing for variant origin clarification was required in one case (0.4%). Most BRCA1/2 germline variants were DV (121 DVs, 26 VUS); in non-BRCA1/2 genes, VUS were ubiquitous (53 DVs, 132 VUS). In silico prediction identified additional 4.9% HRR and 1.2% BRCA1/2/ATM/CHEK2 mutation patients.

Conclusions

Our study represents one of the first reports about the incidence of DV and VUS in HRR genes, including genes beyond BRCA1/2, identified in cancer patients from Russia, assessed by NGS. In silico predictions of the observed HRR gene variants suggest that non-BRCA gene testing is likely to result in higher frequency of patients who are candidates for PARP inhibitor therapy. Continuing sequencing efforts should clarify interpretation of frequently observed non-BRCA VUS.

RNA Panel Sequencing Is an Effective Tool to Help Classify Splice Variants for Clinical Oncogenetic Diagnosis

Routine gene panel analysis identifies pathogenic variants in clinically relevant genes. However, variants of unknown significance (VUSs) are commonly observed, many of which potentially have an impact on mRNA transcription and splicing. Several software programs attempt to predict the impact of variants on splicing and thus make it possible to select the variants for which it is important to study the effect on the transcripts. Transcript analysis is also necessary to show the tandem character of large duplications, and it can be useful for the search for deep intronic variants that are difficult to identify in a DNA panel. We analyzed 53 variants of unknown significance by targeted sequencing of 48 genes using RNA extracted from patient blood samples. RT-PCR and Sanger sequencing of patient mRNA or minigene monoallelic analysis was also carried out when necessary. For the 53 VUSs, 21 could be classified as likely neutral and 10 as pathogenic or likely pathogenic. Data are comprehensively presented for four variants: PTEN c.206+6T>G, MLH1 c.791-489_791-20del, BRCA2 c.68-8_68-7delinsAA, and MSH2 c.(1076+1_1077-1)_(1276+1_1277-1)dup. These four examples illustrate the usefulness of blood RNA panel sequencing in clinical oncogenetics to help classify VUSs with predicted splice effects. It could also be useful for characterizing large duplications and for detecting deep intronic variants with an impact on expressed transcripts.

Profiling of the genetic features of patients with breast, ovarian, colorectal and extracolonic cancers: Association to CHEK2 and PALB2 germline mutations

BACKGROUND AND AIMS

Cancer predisposition goes beyond BRCA and DNA Mismatch Repair (MMR) genes since multi-gene panel testing has become the routine diagnostic tool for hereditary cancer suspicion (HCS) cases. CHEK2 and PALB2 are some of the foremost-mutated non-BRCA/MMR actionable genes in families with a significant familial aggregation. Therefore, the purpose of this work is to unravel which tumours other than breast, ovary or colorectal display the patients.

MATERIALS AND METHODS

We have analysed 528 probands that meet the inclusion criteria for Hereditary Breast and Ovarian Cancer and Lynch Syndrome established by our Hereditary Cancer Regional Program with a customized 35 genes-panel by using Ion Torrent™ Technology.

RESULTS

We have identified pathogenic variants (PVs) in 61 families (1.55%), of which more than half (31 probands) harboured PVs in CHEK2 and PALB2 genes. Ours results reveal that not only were PVs CHEK2 and PALB2 carriers more likely to have family history of cancer not limited to breast, ovarian or colorectal cancers, but also they are prone to other extracolonic cancers, noteworthy endometrial and gastric cancers.

CONCLUSIONS

Multigene panel testing improves the chance of finding PVs in actionable genes in families with HCS. In addition, the coexistence of variants should be recorded to implement a polygenic risk algorithm that might explain the missing heritability in the aforementioned families.

Parallel DNA/RNA NGS Using an Identical Target Enrichment Panel in the Analysis of Hereditary Cancer Predisposition

Germline DNA testing using the next-gene-ration sequencing (NGS) technology has become the analytical standard for the diagnostics of hereditary diseases, including cancer. Its increasing use places high demands on correct sample identification, independent confirmation of prioritized variants, and their functional and clinical interpretation. To streamline these processes, we introduced parallel DNA and RNA capture-based NGS using identical capture panel CZECANCA, which is routinely used for DNA analysis of hereditary cancer predisposition. Here, we present the analytical workflow for RNA sample processing and its analytical and diagnostic performance. Parallel DNA/RNA analysis allowed credible sample identification by calculating the kinship coefficient. The RNA capture-based approach enriched transcriptional targets for the majority of clinically relevant cancer predisposition genes to a degree that allowed analysis of the effect of identified DNA variants on mRNA processing. By comparing the panel and whole-exome RNA enrichment, we demonstrated that the tissue-specific gene expression pattern is independent of the capture panel. Moreover, technical replicates confirmed high reproducibility of the tested RNA analysis. We concluded that parallel DNA/RNA NGS using the identical gene panel is a robust and cost-effective diagnostic strategy. In our setting, it allows routine analysis of 48 DNA/RNA pairs using NextSeq 500/550 Mid Output Kit v2.5 (150 cycles) in a single run with sufficient coverage to analyse 226 cancer predisposition and candidate ge-nes. This approach can replace laborious Sanger confirmatory sequencing, increase testing turnaround, reduce analysis costs, and improve interpretation of the impact of variants by analysing their effect on mRNA processing.

Using the ACMG/AMP framework to capture evidence related to predicted and observed impact on splicing: Recommendations from the ClinGen SVI Splicing Subgroup

The American College of Medical Genetics and Genomics (ACMG)/Association for Molecular Pathology (AMP) framework for classifying variants uses six evidence categories related to the splicing potential of variants: PVS1, PS3, PP3, BS3, BP4, and BP7. However, the lack of guidance on how to apply such codes has contributed to variation in the specifications developed by different Clinical Genome Resource (ClinGen) Variant Curation Expert Panels. The ClinGen Sequence Variant Interpretation Splicing Subgroup was established to refine recommendations for applying ACMG/AMP codes relating to splicing data and computational predictions. We utilized empirically derived splicing evidence to (1) determine the evidence weighting of splicing-related data and appropriate criteria code selection for general use, (2) outline a process for integrating splicing-related considerations when developing a gene-specific PVS1 decision tree, and (3) exemplify methodology to calibrate splice prediction tools. We propose repurposing the PVS1_Strength code to capture splicing assay data that provide experimental evidence for variants resulting in RNA transcript(s) with loss of function. Conversely, BP7 may be used to capture RNA results demonstrating no splicing impact for intronic and synonymous variants. We propose that the PS3/BS3 codes are applied only for well-established assays that measure functional impact not directly captured by RNA-splicing assays. We recommend the application of PS1 based on similarity of predicted RNA-splicing effects for a variant under assessment in comparison with a known pathogenic variant. The recommendations and approaches for consideration and evaluation of RNA-assay evidence described aim to help standardize variant pathogenicity classification processes when interpreting splicing-based evidence.

APPLICATION OF THE ACMG/AMP FRAMEWORK TO CAPTURE EVIDENCE RELEVANT TO PREDICTED AND OBSERVED IMPACT ON SPLICING: RECOMMENDATIONS FROM THE CLINGEN SVI SPLICING SUBGROUP

The American College of Medical Genetics and Genomics (ACMG) and the Association for Molecular Pathology (AMP) framework for classifying variants uses six evidence categories related to the splicing potential of variants: PVS1 (null variant in a gene where loss-of-function is the mechanism of disease), PS3 (functional assays show damaging effect on splicing), PP3 (computational evidence supports a splicing effect), BS3 (functional assays show no damaging effect on splicing), BP4 (computational evidence suggests no splicing impact), and BP7 (silent change with no predicted impact on splicing). However, the lack of guidance on how to apply such codes has contributed to variation in the specifications developed by different Clinical Genome Resource (ClinGen) Variant Curation Expert Panels. The ClinGen Sequence Variant Interpretation (SVI) Splicing Subgroup was established to refine recommendations for applying ACMG/AMP codes relating to splicing data and computational predictions. Our study utilised empirically derived splicing evidence to: 1) determine the evidence weighting of splicing-related data and appropriate criteria code selection for general use, 2) outline a process for integrating splicing-related considerations when developing a gene-specific PVS1 decision tree, and 3) exemplify methodology to calibrate bioinformatic splice prediction tools. We propose repurposing of the PVS1_Strength code to capture splicing assay data that provide experimental evidence for variants resulting in RNA transcript(s) with loss of function. Conversely BP7 may be used to capture RNA results demonstrating no impact on splicing for both intronic and synonymous variants, and for missense variants if protein functional impact has been excluded. Furthermore, we propose that the PS3 and BS3 codes are applied only for well-established assays that measure functional impact that is not directly captured by RNA splicing assays. We recommend the application of PS1 based on similarity of predicted RNA splicing effects for a variant under assessment in comparison to a known Pathogenic variant. The recommendations and approaches for consideration and evaluation of RNA assay evidence described aim to help standardise variant pathogenicity classification processes and result in greater consistency when interpreting splicing-based evidence.

SpliceVault predicts the precise nature of variant-associated mis-splicing

Even for essential splice-site variants that are almost guaranteed to alter mRNA splicing, no current method can reliably predict whether exon-skipping, cryptic activation or multiple events will result, greatly complicating clinical interpretation of pathogenicity. Strikingly, ranking the four most common unannotated splicing events across 335,663 reference RNA-sequencing (RNA-seq) samples (300K-RNA Top-4) predicts the nature of variant-associated mis-splicing with 92% sensitivity. The 300K-RNA Top-4 events correctly identify 96% of exon-skipping events and 86% of cryptic splice sites for 140 clinical cases subject to RNA testing, showing higher sensitivity and positive predictive value than SpliceAI. Notably, RNA re-analyses showed we had missed 300K-RNA Top-4 events for several clinical cases tested before the development of this empirical predictive method. Simply, mis-splicing events that happen around a splice site in RNA-seq data are those most likely to be activated by a splice-site variant. The SpliceVault web portal allows users easy access to 300K-RNA for informed splice-site variant interpretation and classification.

Mutational and splicing landscape in a cohort of 43,000 patients tested for hereditary cancer

DNA germline genetic testing can identify individuals with cancer susceptibility. However, DNA sequencing alone is limited in its detection and classification of mRNA splicing variants, particularly those located far from coding sequences. Here we address the limitations of splicing variant identification and interpretation by pairing DNA and RNA sequencing and describe the mutational and splicing landscape in a clinical cohort of 43,524 individuals undergoing genetic testing for hereditary cancer predisposition.

Effects of Laparoscopic Hyperthermic Perfusion Therapy Combined with Adjuvant Treatment of Compound Yew Capsule on Ovarian Blood Flow Parameters and Immune Function in Patients with Ovarian Cancer

Objective

To determine the effects of laparoscopic hyperthermic perfusion therapy combined with adjuvant compound yew capsules on ovarian blood flow parameters and immune function in patients with ovarian cancer (OC).

Methods

A total of 90 OC patients enrolled in our hospital between January 2019 and January 2020 were randomly distributed into the control (Con group) and experimental group (Exp group) based on the sealed envelope method. The Con group was administered laparoscopic hyperthermic perfusion therapy. On this basis, the Exp group was subjected to compound yew capsules; the ovarian blood flow parameters and immune function indexes were compared between the two groups.

Results

The Exp group was reported to perform better than the Con group regarding ovarian blood flow parameters and immune indexes after treatment (p < 0.001).

Conclusion

Laparoscopic hyperthermic perfusion therapy combined with adjuvant compound yew capsules for patients with OC can substantially improve the clinical indexes and immune function. Furthermore, research and adequate promotion are needed to elicit the evidence beyond preclinical studies to understand the intricacies of its implementation.

Identification of Spliceogenic Variants beyond Canonical GT-AG Splice Sites in Hereditary Cancer Genes

Pathogenic/likely pathogenic variants in susceptibility genes that interrupt RNA splicing are a well-documented mechanism of hereditary cancer syndromes development. However, if RNA studies are not performed, most of the variants beyond the canonical GT-AG splice site are characterized as variants of uncertain significance (VUS). To decrease the VUS burden, we have bioinformatically evaluated all novel VUS detected in 732 consecutive patients tested in the routine genetic counseling process. Twelve VUS that were predicted to cause splicing defects were selected for mRNA analysis. Here, we report a functional characterization of 12 variants located beyond the first two intronic nucleotides using RNAseq in APC, ATM, FH, LZTR1, MSH6, PALB2, RAD51C, and TP53 genes. Based on the analysis of mRNA, we have successfully reclassified 50% of investigated variants. 25% of variants were downgraded to likely benign, whereas 25% were upgraded to likely pathogenic leading to improved clinical management of the patient and the family members.

Minigene Splicing Assays Identify 20 Spliceogenic Variants of the Breast/Ovarian Cancer Susceptibility Gene RAD51C

RAD51C loss-of-function variants are associated with an increased risk of breast and ovarian cancers. Likewise, splicing disruptions are a frequent mechanism of gene inactivation. Taking advantage of a previous splicing-reporter minigene with exons 2-8 (mgR51C_ex2-8), we proceeded to check its impact on the splicing of candidate ClinVar variants. A total of 141 RAD51C variants at the intron/exon boundaries were analyzed with MaxEntScan. Twenty variants were selected and genetically engineered into the wild-type minigene. All the variants disrupted splicing, and 18 induced major splicing anomalies without any trace or minimal amounts (<2.4%) of the minigene full-length (FL) transcript. Twenty-seven transcripts (including the wild-type and r.904A FL transcripts) were identified by fluorescent fragment electrophoresis; of these, 14 were predicted to truncate the RAD51C protein, 3 kept the reading frame, and 8 minor isoforms (1.1−4.7% of the overall expression) could not be characterized. Finally, we performed a tentative interpretation of the variants according to an ACMG/AMP (American College of Medical Genetics and Genomics/Association for Molecular Pathology)-based classification scheme, classifying 16 variants as likely pathogenic. Minigene assays have been proven as valuable tools for the initial characterization of potential spliceogenic variants. Hence, minigene mgR51C_ex2-8 provided useful splicing data for 40 RAD51C variants.

Empirical prediction of variant-activated cryptic splice donors using population-based RNA-Seq data

Predicting which cryptic-donors may be activated by a splicing variant in patient DNA is notoriously difficult. Through analysis of 5145 cryptic-donors (versus 86,963 decoy-donors not used; any GT or GC), we define an empirical method predicting cryptic-donor activation with 87% sensitivity and 95% specificity. Strength (according to four algorithms) and proximity to the annotated-donor appear important determinants of cryptic-donor activation. However, other factors such as splicing regulatory elements, which are difficult to identify, play an important role and are likely responsible for current prediction inaccuracies. We find that the most frequently recurring natural mis-splicing events at each exon-intron junction, summarised over 40,233 RNA-sequencing samples (40K-RNA), predict with accuracy which cryptic-donor will be activated in rare disease. 40K-RNA provides an accurate, evidence-based method to predict variant-activated cryptic-donors in genetic disorders, assisting pathology consideration of possible consequences of a variant for the encoded protein and RNA diagnostic testing strategies.

Genetic variants affecting the consensus splicing motifs can alter binding of spliceosomal components and induce mis-splicing. Here, the authors develop a method, showing that ranking the most common recurring mis-splicing events in public RNA-Seq data can predict the activation of cryptic-donors.

BRCA1 and Breast Cancer: Molecular Mechanisms and Therapeutic Strategies

Breast cancer susceptibility gene 1 (BRCA1) is a tumor suppressor gene, which is mainly involved in the repair of DNA damage, cell cycle regulation, maintenance of genome stability, and other important physiological processes. Mutations or defects in the BRCA1 gene significantly increase the risk of breast, ovarian, prostate, and other cancers in carriers. In this review, we summarized the molecular functions and regulation of BRCA1 and discussed recent insights into the detection and treatment of BRCA1 mutated breast cancer.

Pathogenic neurofibromatosis type 1 (NF1) RNA splicing resolved by targeted RNAseq

Neurofibromatosis type 1 (NF1) is caused by loss-of-function variants in the NF1 gene. Approximately 10% of these variants affect RNA splicing and are either missed by conventional DNA diagnostics or are misinterpreted by in silico splicing predictions. Therefore, a targeted RNAseq-based approach was designed to detect pathogenic RNA splicing and associated pathogenic DNA variants. For this method RNA was extracted from lymphocytes, followed by targeted RNAseq. Next, an in-house developed tool (QURNAs) was used to calculate the enrichment score (ERS) for each splicing event. This method was thoroughly tested using two different patient cohorts with known pathogenic splice-variants in NF1. In both cohorts all 56 normal reference transcript exon splice junctions, 24 previously described and 45 novel non-reference splicing events were detected. Additionally, all expected pathogenic splice-variants were detected. Eleven patients with NF1 symptoms were subsequently tested, three of which have a known NF1 DNA variant with a putative effect on RNA splicing. This effect could be confirmed for all 3. The other eight patients were previously without any molecular confirmation of their NF1-diagnosis. A deep-intronic pathogenic splice variant could now be identified for two of them (25%). These results suggest that targeted RNAseq can be successfully used to detect pathogenic RNA splicing variants in NF1.

Classification of MSH6 Variants of Uncertain Significance Using Functional Assays

Lynch syndrome (LS) is one of the most common hereditary cancer predisposition syndromes worldwide. Individuals with LS have a high risk of developing colorectal or endometrial cancer, as well as several other cancers. LS is caused by autosomal dominant pathogenic variants in one of the DNA mismatch repair (MMR) genes MLH1, MSH2, PMS2 or MSH6, and typically include truncating variants, such as frameshift, nonsense or splicing variants. However, a significant number of missense, intronic, or silent variants, or small in-frame insertions/deletions, are detected during genetic screening of the MMR genes. The clinical effects of these variants are often more difficult to predict, and a large fraction of these variants are classified as variants of uncertain significance (VUS). It is pivotal for the clinical management of LS patients to have a clear genetic diagnosis, since patients benefit widely from screening, preventive and personal therapeutic measures. Moreover, in families where a pathogenic variant is identified, testing can be offered to family members, where non-carriers can be spared frequent surveillance, while carriers can be included in cancer surveillance programs. It is therefore important to reclassify VUSs, and, in this regard, functional assays can provide insight into the effect of a variant on the protein or mRNA level. Here, we briefly describe the disorders that are related to MMR deficiency, as well as the structure and function of MSH6. Moreover, we review the functional assays that are used to examine VUS identified in MSH6 and discuss the results obtained in relation to the ACMG/AMP PS3/BS3 criterion. We also provide a compiled list of the MSH6 variants examined by these assays. Finally, we provide a future perspective on high-throughput functional analyses with specific emphasis on the MMR genes.

New Approach for Detection of Normal Alternative Splicing Events and Aberrant Spliceogenic Transcripts with Long-Range PCR and Deep RNA Sequencing

RNA sequencing is a promising technique for detecting normal and aberrant RNA isoforms. Here, we present a new single-gene, straightforward 1-day hands-on protocol for detection of splicing alterations with deep RNA sequencing from blood. We have validated our method's accuracy by detecting previously published normal splicing isoforms of STK11 gene. Additionally, the same technique was used to provide the first comprehensive catalogue of naturally occurring alternative splicing events of the NBN gene in blood. Furthermore, we demonstrate that our approach can be used for detection of splicing impairment caused by genetic variants. Therefore, we were able to reclassify three variants of uncertain significance: NBN:c.584G>A, STK11:c.863-5_863-3delCTC and STK11:c.615G>A. Due to the simplicity of our approach, it can be incorporated into any molecular diagnostics laboratory for determination of variant's impact on splicing.

BRCA1 and BRCA2 whole cDNA analysis in unsolved hereditary breast/ovarian cancer patients

Germline pathogenic variants in BRCA1 and BRCA2 genes (BRCA1/2) explain an important fraction of hereditary breast/ovarian cancer (HBOC) cases. Genetic testing generally involves examining coding regions and exon/intron boundaries, thus the frequency of deleterious variants in non-coding regions is unknown. Here we analysed BRCA1/2 whole cDNA in a large cohort of 320 unsolved high-risk HBOC cases in order to identify potential splicing alterations explained by variants in BRCA1/2 deep intronic regions. Whole RNA splicing profiles were analysed by RT-PCR using Sanger sequencing or high-resolution electrophoresis in a QIAxcel instrument. Known predominant BRCA1/2 alternative splicing events were detected, together with two novel events BRCA1 ▼21 and BRCA2 Δ18q_27p. BRCA2 exon 3 skipping was detected in one patient (male) affected with breast cancer, caused by a known Portuguese founder mutation (c.156_157insAluYa5). An altered BRCA2 splicing pattern was detected in three patients, consisting in the up-regulation of ▼20A, Δ22 and ▼20A+Δ22 transcripts. In silico analysis and semi-quantitative data identified the polymorphism BRCA2 c.8755-66T>C as a potential modifier of Δ22 levels. Our findings suggest that mRNA alterations in BRCA1/2 caused by deep intronic variants are rare in Spanish population. However, RNA analysis complements DNA-based strategies allowing the identification of alterations that could go undetected by conventional testing.

RAD51D Aberrant Splicing in Breast Cancer: Identification of Splicing Regulatory Elements and Minigene-Based Evaluation of 53 DNA Variants

RAD51D loss-of-function variants increase lifetime risk of breast and ovarian cancer. Splicing disruption is a frequent pathogenic mechanism associated with variants in susceptibility genes. Herein, we have assessed the splicing and clinical impact of splice-site and exonic splicing enhancer (ESE) variants identified through the study of ~113,000 women of the BRIDGES cohort. A RAD51D minigene with exons 2-9 was constructed in splicing vector pSAD. Eleven BRIDGES splice-site variants (selected by MaxEntScan) were introduced into the minigene by site-directed mutagenesis and tested in MCF-7 cells. The 11 variants disrupted splicing, collectively generating 25 different aberrant transcripts. All variants but one produced negligible levels (<3.4%) of the full-length (FL) transcript. In addition, ESE elements of the alternative exon 3 were mapped by testing four overlapping exonic microdeletions (≥30-bp), revealing an ESE-rich interval (c.202_235del) with critical sequences for exon 3 recognition that might have been affected by germline variants. Next, 26 BRIDGES variants and 16 artificial exon 3 single-nucleotide substitutions were also assayed. Thirty variants impaired splicing with variable amounts (0-65.1%) of the FL transcript, although only c.202G>A demonstrated a complete aberrant splicing pattern without the FL transcript. On the other hand, c.214T>C increased efficiency of exon 3 recognition, so only the FL transcript was detected (100%). In conclusion, 41 RAD51D spliceogenic variants (28 of which were from the BRIDGES cohort) were identified by minigene assays. We show that minigene-based mapping of ESEs is a powerful approach for identifying ESE hotspots and ESE-disrupting variants. Finally, we have classified nine variants as likely pathogenic according to ACMG/AMP-based guidelines, highlighting the complex relationship between splicing alterations and variant interpretation.

Spectrum of splicing variants in disease genes and the ability of RNA analysis to reduce uncertainty in clinical interpretation

The complexities of gene expression pose challenges for the clinical interpretation of splicing variants. To better understand splicing variants and their contribution to hereditary disease, we evaluated their prevalence, clinical classifications, and associations with diseases, inheritance, and functional characteristics in a 689,321-person clinical cohort and two large public datasets. In the clinical cohort, splicing variants represented 13% of all variants classified as pathogenic (P), likely pathogenic (LP), or variants of uncertain significance (VUSs). Most splicing variants were outside essential splice sites and were classified as VUSs. Among all individuals tested, 5.4% had a splicing VUS. If RNA analysis were to contribute supporting evidence to variant interpretation, we estimated that splicing VUSs would be reclassified in 1.7% of individuals in our cohort. This would result in a clinically significant result (i.e., P/LP) in 0.1% of individuals overall because most reclassifications would change VUSs to likely benign. In ClinVar, splicing VUSs were 4.8% of reported variants and could benefit from RNA analysis. In the Genome Aggregation Database (gnomAD), splicing variants comprised 9.4% of variants in protein-coding genes; most were rare, precluding unambiguous classification as benign. Splicing variants were depleted in genes associated with dominant inheritance and haploinsufficiency, although some genes had rare variants at essential splice sites or had common splicing variants that were most likely compatible with normal gene function. Overall, we describe the contribution of splicing variants to hereditary disease, the potential utility of RNA analysis for reclassifying splicing VUSs, and how natural variation may confound clinical interpretation of splicing variants.

FOXI1 inhibits gastric cancer cell proliferation by activating miR-590/ATF3 axis via integrating ChIP-seq and RNA-seq data

FOXI1 plays a key role in the development of gastric cancer. However, the whole genome FOXI1 binding sites and its target genes are unclear. In the present study, we used ChIP-seq and RNA-seq technologies to identify the target gene of FOXI1. Firstly, ChIP-seq data showed that, 4476 unique peaks in the genome region were captured. Most of these binding peaks are located in introns or intergenic regions. We annotated all the peaks to the nearest gene and identified 404 genes as FOXI1 binding genes. KEGG and GO analysis showed that FOXI1 binding gene to be correlated with the cellular process, cell part, cell, binding, single-organism process. Further, we performed FOXI1-overexpressed RNA-seq experiment. We comprehensively analyzed the ChIP-seq and RNA-seq data and take the intersection of two databases, several genes were identified. ATF3 was selected from the intersection since ATF3 was the most enriched mRNA after FOXI1 overexpressed. ChIP-qPCR and luciferase report gene were used to validate that ATF3 was target gene of FOXI1. Intriguely, ATF3 protein was significantly downregulated after FOXI1 overexpressed. We found FOXI1 can also bind to the promoter of miR-590 and active it which directly target ATF3. The binding site between FOXI1 and miR-590 was verified by ChIP-qPCR and luciferase report gene, and the target relationship between miR-590 and ATF3 was confirmed by dual-luciferase reporter gene. In conclusion, our data identified the genome binding sites of FOXI1, and provide evidence that FOXI1 inhibits gastric cancer cell proliferation by activating miR-590/ATF3 axis.

Interpretation of BRCA2 Splicing Variants: A Case Series of Challenging Variant Interpretations and the Importance of Functional RNA Analysis

A substantial proportion of pathogenic variants associated with an increased risk of hereditary cancer are sequence variants affecting RNA splicing. The classification of these variants can be complex when both non-functional and functional transcripts are produced from the variant allele. We present four BRCA2 splice site variants with complex variant interpretations (BRCA2 c.68-3T>G, c.68-2A>G, c.425G>T, c.8331+2T>C). Evidence supporting a pathogenic classification is available for each variant, including in silico models, absence in population databases, and published functional data. However, comprehensive RNA analysis showed that some functional transcript may be produced by each variant. BRCA2 c.68-3T>G results in a partial splice defect. For BRCA2 c.68-2A>G and c.425G>T, aberrant splicing was shown to produce a potentially functional, in-frame transcript. BRCA2 c.8331+2T>C may utilize a functional GC donor in place of the wild-type GT donor. The severity of cancer history for carriers of these variants was also assessed using a history weighting algorithm and was not consistent with pathogenic controls (carriers of known pathogenic variants in BRCA2). Due to the conflicting evidence, our laboratory classifies these BRCA2 variants as variants of uncertain significance. This highlights the importance of evaluating new and existing evidence to ensure accurate variant classification and appropriate patient care.

Anti-FIRΔexon2 autoantibody as a novel indicator for better overall survival in gastric cancer

There is no clinically available biomarker for efficiently indicating the overall survival or therapy response of gastric cancer (GC). The autoantibodies (Abs) in the sera of anti‐far‐upstream element‐binding protein‐interacting repressor‐lacking exon2 (FIRΔexon2), anti‐sorting nexin 15, and anti‐spermatogenesis and oogenesis–specific basic helix–loop–helix 1 were markedly higher in GC patients than in healthy donors (HDs). These Abs were identified by large‐scale serological identification of antigens by recombinant cDNA expression cloning screenings and their expression levels were evaluated by amplified luminescence proximity homogeneous assay. In particular, compared with age‐matched HDs, the level of anti‐FIRΔexon2 Abs in GC patients was significantly higher (P < .001). The Spearman's rank correlation analysis between anti‐FIRΔexon2 Abs and clinically available tumor markers such as carcinoembryonic antigen (CEA) was statistically insignificant, indicating that FIRΔexon2 Abs is an independent biomarker. We performed receiver‐operating curve analysis to evaluate the anti‐FIRΔexon2 Ab as a candidate biomarker with CEA and carbohydrate antigen 19‐9 (CA19‐9). The overall survival of GC patients with high anti‐FIRΔexon2 Abs titer was significantly favorable (P = .04) than that of GC patients who were below detection level of anti‐FIRΔexon2 Abs. However, clinical stages were not apparently correlated with the levels of anti‐FIRΔexon2 Ab, CEA, and CA19‐9. In conclusion, anti‐FIRΔexon2 Abs detected in GC patients is a potential biomarker for monitoring a better prognosis. Hence, anti‐FIRΔexon2 Abs is a promising biomarker for indicating better overall survival of gastric cancer patients.

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

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

SVA retrotransposon insertion in exon of MMR genes results in aberrant RNA splicing and causes Lynch syndrome

Lynch syndrome is an autosomal dominant hereditary cancer syndrome in which many cancers develop, the main one being colorectal cancer. Germline pathogenic variants in one of four mismatch repair (MMR) genes are known to be causative of this disease. Accurate diagnosis using genetic testing can greatly benefit the health of those affected. Recently, owing to the improvement of sequence techniques, complicated variants affecting the functions of MMR genes were discovered. In this study, we analyzed insertions of a retrotransposon-like sequence in exon 5 of the MSH6 gene and exon 3 of the MSH2 gene found in Japanese families suspected of having Lynch syndrome. Both of these insertions induced aberrant splicing, and these variants were successfully identified by mRNA sequencing or visual observation of mapping results, although a standard DNA-seq analysis pipeline failed to detect them. The insertion sequences were ~2.5 kbp in length and were found to have the structure of an SVA retrotransposon (SVA). One SVA sequence was not present in the hg19 or hg38 reference genome, but was in a Japanese-specific reference sequence (JRGv2). Our study illustrates the difficulties of identifying SVA insertions in disease genes, and that the possibility of polymorphic insertions should be considered when analyzing mobile elements.

Lessons learned from understanding chemotherapy resistance in epithelial tubo-ovarian carcinoma from BRCA1and BRCA2mutation carriers

BRCA1 and BRCA2 are multi-functional proteins and key factors for maintaining genomic stability through their roles in DNA double strand break repair by homologous recombination, rescuing stalled or damaged DNA replication forks, and regulation of cell cycle DNA damage checkpoints. Impairment of any of these critical roles results in genomic instability, a phenotypic hallmark of many cancers including breast and epithelial ovarian carcinomas (EOC). Damaging, usually loss of function germline and somatic variants in BRCA1 and BRCA2, are important drivers of the development, progression, and management of high-grade serous tubo-ovarian carcinoma (HGSOC). However, mutations in these genes render patients particularly sensitive to platinum-based chemotherapy, and to the more innovative targeted therapies with poly-(ADP-ribose) polymerase inhibitors (PARPis) that are targeted to BRCA1/BRCA2 mutation carriers. Here, we reviewed the literature on the responsiveness of BRCA1/2-associated HGSOC to platinum-based chemotherapy and PARPis, and propose mechanisms underlying the frequent development of resistance to these therapeutic agents.

Comparison of suspected Lynch syndrome patients carrying BRCA and BRCA-like variants with Lynch syndrome probands: Phenotypic characteristics and pedigree analyses

Background

Colorectal cancer (CRC) patients diagnosed with Lynch syndrome (LS) are recommended genetic testing. Increasing numbers of germline variants involved in homologous recombination have been identified in suspected LS patients. This study compared phenotypic the characteristics of suspected LS patients carrying BRCA and BRCA‐like variants with those of LS patients.

Methods

Forty‐two patients carrying pathogenic variants of DNA mismatch repair (MMR) genes (MMR group), 9 carrying BRCA variants, and 11 carrying BRCA‐like variants (BRCA/BRCA‐like group) who met LS clinical criteria were enrolled in this study. Clinical characteristics, pedigrees, and survival rates were compared and BRCA variants were analyzed.

Results

The earliest CRC‐onset age and tumor differentiation were higher in the BRCA/BRCA‐like group than in the MMR group. Metachronous CRCs were more numerous in the MMR group, resulting in a higher progression‐free survival rate in the BRCA/BRCA‐like group. Extra‐colorectal cancers were more frequently observed in the BRCA/BRCA‐like group. BRCA2 and BRCA1 variants were clustered in exons 11 and 4/7, respectively.

Conclusion

BRCA and BRCA‐like variants in CRC patients with LS showed moderate penetrance. BRCA/BRCA‐like variant carriers had a higher risk for extra‐colorectal cancers. Surveillance of susceptible organs other than the intestine should be performed for probands and affected family members.