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.

Evolution of Dihydropyrimidine Dehydrogenase Diagnostic Testing in a Single Center during an 8-Year Period of Time.

OBJECTIVE: Fluoropyrimidine treatment can be optimized based on dihydropyrimidine dehydrogenase (DPD) activity. DPD dysfunction leads to increased exposure to active metabolites, which can result in severe or even fatal toxicity. METHODS: We provide an overview of 8 years of DPD diagnostic testing (n = 1194). RESULTS: Within the study period, our diagnostic test evolved from a single-enzyme measurement using first a radiochemical and then a nonradiochemical assay by ultra HPLC-MS in peripheral blood mononuclear cells with uracil, to a combined enzymatic and genetic test (ie, polymerase chain reaction) followed by Sanger sequence analysis of 4 variants of the DPYD gene (ie, DPYD*2A, DPYD*13, c.2846A>T, and 1129-5923C>G; allele frequencies 0.58%, 0.03%, 0.29%, and 1.35%, respectively). Patients who have 1 of the 4 variants tested (n = 814) have lower enzyme activity than the overall patient group. The majority of patients with the DPYD*2A variant (83%) consistently showed decreased enzyme activity. Only 24 (25.3%) of 95 patients (tested for 4 variants) with low enzyme activity carried a variant. Complete DPYD sequencing in a subgroup with low enzyme activity and without DPYD*2A variant (n = 47) revealed 10 genetic variants, of which 4 have not been described previously. We did not observe a strong link between DPYD genotype and enzyme activity. CONCLUSIONS: Previous studies have shown that DPD status should be determined before treatment with fluoropyrimidine agents to prevent unnecessary side effects with possible fatal consequences. Our study in combination with literature shows that there is a discrepancy between the DPD enzyme activity and the presence of clinically relevant single nucleotide polymorphisms. At this moment, a combination of a genetic and enzyme test is preferable for diagnostic testing. (Curr Ther Res Clin Exp. 2018; 79:XXX-XXX).

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.

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.

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.

Aging promotes accumulation of senescent and multiciliated cells in human endometrial epithelium.

STUDY QUESTION: What changes occur in the endometrium during aging, and do they impact fertility? SUMMARY ANSWER: Both the transcriptome and cellular composition of endometrial samples from women of advanced maternal age (AMA) are significantly different from that of samples from young women, suggesting specific changes in epithelial cells that may affect endometrial receptivity. WHAT IS KNOWN ALREADY: Aging is associated with the accumulation of senescent cells in aging tissues. Reproductive aging is mostly attributed to the decline in ovarian reserve and oocyte quality, whereas the endometrium is a unique complex tissue that is monthly renewed under hormonal regulation. Several clinical studies have reported lower implantation and pregnancy rates in oocyte recipients of AMA during IVF. Molecular studies have indicated the presence of specific mutations within the epithelial cells of AMA endometrium, along with altered gene expression of bulk endometrial tissue. STUDY DESIGN SIZE DURATION: Endometrial transcriptome profiling was performed for 44 women undergoing HRT during the assessment of endometrial receptivity before IVF. Patients younger than 28 years were considered as the young maternal age (YMA) group (age 23-27 years) and women older than 45 years were considered as the AMA group (age 47-50 years). Endometrial biopsies were obtained on Day 5 of progesterone treatment and RNA was extracted. All endometrial samples were evaluated as being receptive based on the expression of 68 common endometrial receptivity markers. Endometrial samples from another 24 women classified into four age groups (YMA, intermediate age group 1 (IMA1, age 29-35), intermediate age group 2 (IMA2, age 36-44), and AMA) were obtained in the mid-secretory stage of a natural cycle (NC) and used for validation studies across the reproductive lifespan. PARTICIPANTS/MATERIALS SETTING METHODS: A total of 24 HRT samples (12 YMA and 12 AMA) were subject to RNA sequencing (RNA-seq) and differential gene expression analysis, 20 samples (10 YMA and 10 AMA) were used for qPCR validation, and 24 NC samples (6 YMA, 6 IMA1, 6 IMA2 and 6AMA) were used for RNA-seq validation of AMA genes across the woman's reproductive lifespan. Immunohistochemistry (IHC) was used to confirm some expression changes at the protein level. Computational deconvolution using six endometrial cell type-specific transcriptomic profiles was conducted to compare the cellular composition between the groups. MAIN RESULTS AND THE ROLE OF CHANCE: Comparisons between YMA and AMA samples identified a lower proportion of receptive endometria in the AMA group (P = 0.007). Gene expression profiling identified 491 differentially expressed age-sensitive genes (P adj < 0.05) that revealed the effects of age on endometrial epithelial growth and receptivity, likely contributing to decreased reproductive performance. Our results indicate that changes in the expression of the cellular senescence marker p16INK4a and genes associated with metabolism, inflammation, and hormone response are involved in endometrial aging. Importantly, we demonstrate that the proportion of multi-ciliated cells, as discovered based on RNA-seq data deconvolution and tissue IHC results, is affected by endometrial aging, and propose a putative onset of age-related changes. Furthermore, we propose that aging has an impact on the transcriptomic profile of endometrial tissue in the context of endometrial receptivity. LARGE SCALE DATA: The raw sequencing data reported in this article are deposited at the Gene Expression Omnibus under accession code GSE236128. LIMITATIONS REASONS FOR CAUTION: This retrospective study identified changes in the endometrium of patients undergoing hormonal replacement and validated these changes using samples obtained during a NC. However, future studies must clarify the importance of these findings on the clinical outcomes of assisted reproduction. WIDER IMPLICATIONS OF THE FINDINGS: The findings reported in this study have important implications for devising future strategies aimed at improving fertility management in women of advanced reproductive age. STUDY FUNDING/COMPETING INTERESTS: This research was funded by the Estonian Research Council (grant no. PRG1076), Horizon 2020 innovation grant (ERIN, grant no. EU952516), Enterprise Estonia (grant no. EU48695), MSCA-RISE-2020 project TRENDO (grant no. 101008193), EU 874867 project HUTER, the Horizon Europe NESTOR grant (grant no. 101120075) of the European Commission, the EVA specialty program (grant no. KP111513) of the Maastricht University Medical Center (MUMC+), MICIU/AEI/10.13039/501100011033 and FEDER, EU projects Endo-Map (grant no. PID2021-12728OB-100), ROSY (grant no. CNS2022-135999), and the National Science Fund of Bulgaria (grant no. KII-06 H31/2). The authors declare no competing interests.

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.

MLL2 mutation spectrum in 45 patients with Kabuki syndrome.

Kabuki Syndrome (KS) is a rare syndrome characterized by intellectual disability and multiple congenital abnormalities, in particular a distinct dysmorphic facial appearance. KS is caused by mutations in the MLL2 gene, encoding an H3K4 histone methyl transferase which acts as an epigenetic transcriptional activator during growth and development. Direct sequencing of all 54 exons of the MLL2 gene in 45 clinically well-defined KS patients identified 34 (75.6%) different mutations. One mutation has been described previously, all others are novel. Clinically, all KS patients were sporadic, and mutations were de novo for all 27 families for which both parents were available. We detected nonsense (n=11), frameshift (n=17), splice site (n=4) and missense (n=2) mutations, predicting a high frequency of absent or non-functional MLL2 protein. Interestingly, both missense mutations located in the C-terminal conserved functional domains of the protein. Phenotypically our study indicated a statistically significant difference in the presence of a distinct facial appearance (p=0.0143) and growth retardation (p=0.0040) when comparing KS patients with an MLL2 mutation compared to patients without a mutation. Our data double the number of MLL2 mutations in KS reported so far and widen the spectrum of MLL2 mutations and disease mechanisms in KS.

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

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

MECP2 mutations are an infrequent cause of mental retardation associated with neurological problems in male patients.

Mutations in the methyl-CpG-binding protein 2 (MECP2) gene located on Xq28, cause Rett syndrome (RTT) in female patients. Meanwhile, nonmosaic MECP2 mutations unknown in girls have been found in an increasing number of male patients with a normal 46, XY karyotype. They can cause a broad spectrum of neurodevelopmental disorders which often show a combination of mental retardation (MR) with neurological symptoms. We present the results of MECP2 analysis in a group of 72 male patients with an unexplained combination of MR and neurological features, and review the mutational reports published on male patients since the discovery of the MECP2 gene. Analysis included sequencing of exon 1 which thus far was mostly omitted from DNA screening. One pathogenic mutation has been found in a patient with Rett variant, in addition to an unclassified variant and a series of nonpathogenic changes. No changes have been found in exon 1. Criteria for testing of male patients are classic RTT, severe neonatal encephalopathy, and RTT variant which may be clinically underrecognized. Testing can also be considered in males with a combination of unexplained MR and (progressive) neurological manifestations although the yield of MECP2 analysis is probably low in this situation. Based on the literature, MECP2 testing in males with MR only is debatable.

The unfolding clinical spectrum of holoprosencephaly due to mutations in SHH, ZIC2, SIX3 and TGIF genes.

Holoprosencephaly is a severe malformation of the brain characterized by abnormal formation and separation of the developing central nervous system. The prevalence is 1:250 during early embryogenesis, the live-born prevalence is 1:16 000. The etiology of HPE is extremely heterogeneous and can be teratogenic or genetic. We screened four known HPE genes in a Dutch cohort of 86 non-syndromic HPE index cases, including 53 family members. We detected 21 mutations (24.4%), 3 in SHH, 9 in ZIC2 and 9 in SIX3. Eight mutations involved amino-acid substitutions, 7 ins/del mutations, 1 frame-shift, 3 identical poly-alanine tract expansions and 2 gene deletions. Pathogenicity of mutations was presumed based on de novo character, predicted non-functionality of mutated proteins, segregation of mutations with affected family-members or combinations of these features. Two mutations were reported previously. SNP array confirmed detected deletions; one spanning the ZIC2/ZIC5 genes (approx. 100 kb) the other a 1.45 Mb deletion including SIX2/SIX3 genes. The mutation percentage (24%) is comparable with previous reports, but we detected significantly less mutations in SHH: 3.5 vs 10.7% (P=0.043) and significantly more in SIX3: 10.5 vs 4.3% (P=0.018). For TGIF1 and ZIC2 mutation the rate was in conformity with earlier reports. About half of the mutations were de novo, one was a germ line mosaic. The familial mutations displayed extensive heterogeneity in clinical manifestation. Of seven familial index patients only two parental carriers showed minor HPE signs, five were completely asymptomatic. Therefore, each novel mutation should be considered as a risk factor for clinically manifest HPE, with the caveat of reduced clinical penetrance.