Analysis and interpretation of RNA splicing alterations in genes involved in genetic disorders

Transcript capture and ultradeep long-read RNA sequencing (CAPLRseq) to diagnose HNPCC/Lynch syndrome

PURPOSE

Whereas most human genes encode multiple mRNA isoforms with distinct function, clinical workflows for assessing this heterogeneity are not readily available. This is a substantial shortcoming, considering that up to 25% of disease-causing gene variants are suspected of disrupting mRNA splicing or mRNA abundance. Long-read sequencing can readily portray mRNA isoform diversity, but its sensitivity is relatively low due to insufficient transcriptome penetration.

METHODS

We developed and applied capture-based target enrichment from patient RNA samples combined with Oxford Nanopore long-read sequencing for the analysis of 123 hereditary cancer transcripts (capture and ultradeep long-read RNA sequencing (CAPLRseq)).

RESULTS

Validating CAPLRseq, we confirmed 17 cases of hereditary non-polyposis colorectal cancer/Lynch syndrome based on the demonstration of splicing defects and loss of allele expression of mismatch repair genes MLH1, PMS2, MSH2 and MSH6. Using CAPLRseq, we reclassified two variants of uncertain significance in MSH6 and PMS2 as either likely pathogenic or benign.

CONCLUSION

Our data show that CAPLRseq is an automatable and adaptable workflow for effective transcriptome-based identification of disease variants in a clinical diagnostic setting.

Splicing analyses for variants in MMR genes: best practice recommendations from the European Mismatch Repair Working Group

Over 20% of the DNA mismatch repair (MMR) germline variants in suspected Lynch syndrome patients are classified as variants of uncertain significance (VUS). Well-established functional assays are pivotal for assessing the biological impact of these variants and provide relevant evidence for clinical classification. In our collaborative European Mismatch Repair Working Group (EMMR-WG) we compared three different experimental approaches for evaluating the effect of seven variants on mRNA splicing in MMR genes: (i) RT-PCR of full-length transcripts (FLT), (ii) RT-PCR of targeted transcript sections (TTS), both from patient biological samples and (iii) minigene splicing assays. An overall good concordance was observed between splicing patterns in TTS, FLT and minigene analyses for all variants. The FLT analysis depicted a higher number of different isoforms and mitigated PCR-bias towards shorter isoforms. TTS analyses may miss aberrant isoforms and minigene assays may under/overestimate the severity of certain splicing defects. The interpretation of the experimental findings must be cautious to adequately discriminate abnormal events from physiological complex alternative splicing patterns. A consensus strategy for investigating the impact of MMR variants on splicing was defined. First, RNA should be obtained from patient's cell cultures (such as fresh lymphocyte cultures) incubated with/without a nonsense-mediated decay inhibitor. Second, FLT RT-PCR analysis is recommended to oversee all generated isoforms. Third, TTS analysis and minigene assays are useful independent approaches for verifying and clarifying FLT results. The use of several methodologies is likely to increase the strength of the experimental evidence which contributes to improve variant interpretation.

Systematic Analysis of Splice-Site-Creating Mutations in Cancer

For the past decade, cancer genomic studies have focused on mutations leading to splice-site disruption, overlooking those having splice-creating potential. Here, we applied a bioinformatic tool, MiSplice, for the large-scale discovery of splice-site-creating mutations (SCMs) across 8,656 TCGA tumors. We report 1,964 originally mis-annotated mutations having clear evidence of creating alternative splice junctions. TP53 and GATA3 have 26 and 18 SCMs, respectively, and ATRX has 5 from lower-grade gliomas. Mutations in 11 genes, including PARP1, BRCA1, and BAP1, were experimentally validated for splice-site-creating function. Notably, we found that neoantigens induced by SCMs are likely several folds more immunogenic compared to missense mutations, exemplified by the recurrent GATA3 SCM. Further, high expression of PD-1 and PD-L1 was observed in tumors with SCMs, suggesting candidates for immune blockade therapy. Our work highlights the importance of integrating DNA and RNA data for understanding the functional and the clinical implications of mutations in human diseases.

Families with BAP1-Tumor Predisposition Syndrome in The Netherlands: Path to Identification and a Proposal for Genetic Screening Guidelines

Germline pathogenic variants in the BRCA1-associated protein-1 (BAP1) gene cause the BAP1-tumor predisposition syndrome (BAP1-TPDS, OMIM 614327). BAP1-TPDS is associated with an increased risk of developing uveal melanoma (UM), cutaneous melanoma (CM), malignant mesothelioma (MMe), renal cell carcinoma (RCC), meningioma, cholangiocarcinoma, multiple non-melanoma skin cancers, and BAP1-inactivated nevi. Because of this increased risk, it is important to identify patients with BAP1-TPDS. The associated tumors are treated by different medical disciplines, emphasizing the need for generally applicable guidelines for initiating genetic analysis. In this study, we describe the path to identification of BAP1-TPDS in 21 probands found in the Netherlands and the family history at the time of presentation. We report two cases of de novo BAP1 germline mutations (2/21, 9.5%). Findings of this study combined with previously published literature, led to a proposal of guidelines for genetic referral. We recommend genetic analysis in patients with ≥2 BAP1-TPDS-associated tumors in their medical history and/or family history. We also propose to test germline BAP1 in patients diagnosed with UM <40 years, CM <18 years, MMe <50 years, or RCC <46 years. Furthermore, other candidate susceptibility genes for tumor types associated with BAP1-TPDS are discussed, which can be included in gene panels when testing patients.

Full-length transcript amplification and sequencing as universal method to test mRNA integrity and biallelic expression in mismatch repair genes

In pathogenicity assessment, RNA-based analyses are important for the correct classification of variants, and require gene-specific cut-offs for allelic representation and alternative/aberrant splicing. Beside this, the diagnostic yield of RNA-based techniques capable to detect aberrant splicing or allelic loss due to intronic/regulatory variants has to be elaborated. We established a cDNA analysis for full-length transcripts (FLT) of the four DNA mismatch repair (MMR) genes to investigate the splicing pattern and transcript integrity with active/inhibited nonsense-mediated mRNA-decay (NMD). Validation was based on results from normal controls, samples with premature termination codons (PTC), samples with splice-site defects (SSD), and samples with pathogenic putative missense variants. The method was applied to patients with variants of uncertain significance (VUS) or unexplained immunohistochemical MMR deficiency. We categorized the allelic representation into biallelic (50 ± 10%) or allelic loss (≤10%), and >10% and <40% as unclear. We defined isoforms up to 10% and exon-specific exceptions as alternative splicing, set the cut-off for SSD in cDNA + P to 30-50%, and regard >10% and <30% as unclear. FLT cDNA analyses designated 16% of all putative missense variants and 12% of VUS as SSD, detected MMR-defects in 19% of the unsolved patients, and re-classified >30% of VUS. Our method allows a standardized, systematic cDNA analysis of the MMR FLTs to assess the pathogenicity mechanism of VUS on RNA level, which will gain relevance for precision medicine and gene therapy. Diagnostic accuracy will be enhanced by detecting MMR defects in hitherto unsolved patients. The data generated will help to calibrate a high-throughput NGS-based mRNA-analysis and optimize prediction programs.

Clinical features and mutation spectrum in Chinese patients with CADASIL: A multicenter retrospective study

AIM

To characterize clinical features and mutation spectrum in Chinese patients with CADASIL.

METHODS

We collected 261 clinically suspected Chinese CADASIL patients from three hospitals located in different regions of China. Sanger sequencing is performed to screen the exons 2 to 24 of NOTCH3 gene. Clinical and genetic data were retrospectively studied. Haplotype analyses were performed in patients carrying p.Arg544Cys and p.Arg607Cys, respectively.

RESULTS

A total of 214 patients were finally genetically diagnosed as CADASIL, with 45 known NOTCH3 mutations and a novel c.1817G>T mutation. We found that patients carrying p.Arg607Cys or p.Arg544Cys mutation located in exon 11 occupied nearly 35% in our mutation spectrum. In retrospectively study of clinical data, we found a higher number of patients having cognitive impairment and a lower number of patients having migraine with aura. Furthermore, we identified that patients carrying mutations in exon 11 seemed to experience a later disease onset (p=6.8×10^-5 ). Additionally, a common haplotype was found in patients from eastern China carrying p.Arg607Cys, and the patients from Fujian carrying p.Arg544Cys shared the same haplotype with patients from Taiwan carrying p.Arg544Cys.

CONCLUSIONS

These findings broaden the mutational and clinical spectrum of CADASIL and provide additional evidences for the existence of founder effect in CADASIL patients.

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

BACKGROUND

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Comprehensive Mutation Analysis of PMS2 in a Large Cohort of Probands Suspected of Lynch Syndrome or Constitutional Mismatch Repair Deficiency Syndrome

Monoallelic PMS2 germline mutations cause 5%-15% of Lynch syndrome, a midlife cancer predisposition, whereas biallelic PMS2 mutations cause approximately 60% of constitutional mismatch repair deficiency (CMMRD), a rare childhood cancer syndrome. Recently improved DNA- and RNA-based strategies are applied to overcome problematic PMS2 mutation analysis due to the presence of pseudogenes and frequent gene conversion events. Here, we determined PMS2 mutation detection yield and mutation spectrum in a nationwide cohort of 396 probands. Furthermore, we studied concordance between tumor IHC/MSI (immunohistochemistry/microsatellite instability) profile and mutation carrier state. Overall, we found 52 different pathogenic PMS2 variants explaining 121 Lynch syndrome and nine CMMRD patients. In vitro mismatch repair assays suggested pathogenicity for three missense variants. Ninety-one PMS2 mutation carriers (70%) showed isolated loss of PMS2 in their tumors, for 31 (24%) no or inconclusive IHC was available, and eight carriers (6%) showed discordant IHC (presence of PMS2 or loss of both MLH1 and PMS2). Ten cases with isolated PMS2 loss (10%; 10/97) harbored MLH1 mutations. We confirmed that recently improved mutation analysis provides a high yield of PMS2 mutations in patients with isolated loss of PMS2 expression. Application of universal tumor prescreening methods will however miss some PMS2 germline mutation carriers.

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

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

Clinical significance of mismatch repair gene expression in sporadic colorectal cancer

Mismatch repair (MMR) genes play an important role in the occurrence and development of sporadic colorectal cancer; however, the effect of MMR genes on clinicopathological features and prognosis remains unclear. The aim of the present study was to observe the clinical significance of MMR gene expression in sporadic colorectal cancer. Clinicopathological data and postoperative samples from 404 patients with sporadic colorectal cancer were obtained from the Affiliated Tumor Hospital of Xinjiang Medical University. The immunohistochemistry PV-9000 two-step method was performed to measure the protein expression of human mutL homolog 1 (hMLH1), human mutS homolog (hMSH) 2, human postmeiotic segregation increased 2 (hPSM2) and hMSH6. Differences in clinicopathological features, family history and survival time subsequent to surgery between groups with normal and aberrant MMR protein (MMRP) expression were compared. A total of 27.23% of all patients showed aberrant nuclear staining of MMRP. Among the patients with aberrant MMRP expression, a higher proportion of patients showed aberrant expression of more than one type of MMRP than aberrant expression of only one type of MMRP. Aberrant expression of hMLH1/hPSM2 was most commonly observed (29/404). In addition, aberrant MMRP expression in colorectal cancer was indicated predominantly in the right hemicolon. Histological type primarily showed mucinous adenocarcinoma. In addition, with increasing body mass index (BMI), the MMRP deficiency rate was also shown to increase gradually. There was a close association between MMRP expression deficiency and family history of cancer (P<0.05). For TNM stage III patients, the Kaplan-Meier survival curve showed that the aberrant MMRP expression group had a three-year disease-free survival (DFS) rate of 66.67%, which was longer than the DFS rate of the normal group (55.41%), with no statistical difference (P>0.05). In conclusion, the immunohistochemistry PV-9000 two-step method can be used to measure MMRP expression in colorectal cancer. Aberrant MMRP expression is closely correlated with tumor location, histological type, BMI and tumor family history in sporadic colorectal cancer. Aberrant MMRP expression may have an effect on the prognosis of stage III patients.