Screening BRCA1 and BRCA2 unclassified variants for splicing mutations using reverse transcription PCR on patient RNA and an ex vivo assay based on a splicing reporter minigene

A Compound Heterozygous Pathogenic Variant in ZP2 Gene Causes Female Infertility

The oocyte maturation defect 6 is an autosomal recessive hereditary disease caused by a homozygous variant in ZP2 gene. It is characterized by female primary infertility due to an abnormally thin zona pellucida (ZP) and defective sperm binding. Here we identified a compound heterozygous variant (c.1924C > T and c.1695-2A > G) in ZP2 gene in a Chinese Han family. Quantitative real-time PCR showed that the variant c.1924C > T significantly decreased the expression of truncated ZP2 message RNA by the nonsense-mediated decay pathway. Minigene assays showed the c.1695-2A > G variant led to an extra-61-nt preservation of intron 15 at the junction between exons 15 and 16 during transcription. Both variants (c.1924C > T and c.1695-2A > G) resulted in truncated ZP2 proteins (p.R642X and p.C566Hfs*2) that lost the transmembrane domain, which prevented the secretion of the mutant ZP2 proteins and produced a structurally abnormal ZP, thus resulting in female infertility. This study further elucidated the pathogenic mechanism of these two variants and provided new support for the genetic diagnosis of female infertility.

Comprehensive functional splicing analysis of non-canonical CNGB3 variants using in vitro minigene splice assays

Variants in the CNGB3 gene, encoding the B3-subunit of the cone photoreceptor cyclic nucleotide gated channel, are a major cause of autosomal recessive achromatopsia, a rare inherited retinal disease. The mutation spectrum of achromatopsia-associated CNGB3 variants comprises all types of mutations, including those that are straightforward to evaluate in molecular genetic diagnostics, such as frame-shifting, nonsense, and canonical splice site variants. Additionally, variants have been identified within splice regions outside the conserved ±1,2 splice site dinucleotides, making their potential impact on disease association challenging to interpret. This poses a major hurdle for clinical interpretation of causality between the patient's genotype and the proposed clinical diagnosis, but also for the inclusion of such patients into clinical trials for gene augmentation therapy, for which only patients with confirmed (likely) pathogenic CNGB3 variants are eligible. We here performed comprehensive genetic functional analysis of 21 candidate spliceogenic CNGB3 variants-15 reported and 6 novel variants-by means of in vitro minigene splice assays and cDNA analysis, and characterization of spliceogenic events by subcloning, Sanger-sequencing, and capillary fragment analysis. For 16 variants, an impact on splicing was confirmed, supporting the reclassification of 86% of variants of uncertain significance as likely pathogenic or pathogenic according to the ACMG/AMP guidelines. This reclassification enables the confirmation of patients' genotypes, both retrospectively and prospectively. © 2025 The Author(s). The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Spectrum and genotyping strategies of "dark" genetic matter in germline susceptibility genes of tumor syndromes

PURPOSE

Despite the widespread use of high-throughput genotyping strategies, certain mutation types remain understudied. We provide an overview of these often overlooked mutation types, with representative examples from common hereditary cancer syndromes.

METHODS

We conducted a comprehensive review of the literature and locus-specific variant databases to summarize the germline pathogenic variants discovered through non-routine genotyping methods. We evaluated appropriate detection and analysis methods tailored for these specific genetic aberrations. Additionally, we performed in silico splice predictions on deep intronic variants registered in the ClinVar database.

RESULTS

Our study suggests that, aside from founder mutations, most cases are sporadic. However, we anticipate a relatively high likelihood of splice effects for deep intronic variants. The findings underscore the significant clinical utility of genome sequencing techniques and the importance of applying relevant analysis methods.

Defects in meiosis I contribute to the genesis of androgenetic hydatidiform moles

To identify novel genes responsible for recurrent hydatidiform moles (HMs), we performed exome sequencing on 75 unrelated patients who were negative for mutations in the known genes. We identified biallelic deleterious variants in 6 genes, FOXL2, MAJIN, KASH5, SYCP2, MEIOB, and HFM1, in patients with androgenetic HMs, including a familial case of 3 affected members. Five of these genes are essential for meiosis I, and their deficiencies lead to premature ovarian insufficiency. Advanced maternal age is the strongest risk factor for sporadic androgenetic HM, which affects 1 in every 600 pregnancies. We studied Hfm1-/- female mice and found that these mice lost all their oocytes before puberty but retained some at younger ages. Oocytes from Hfm1-/- mice initiated meiotic maturation and extruded the first polar bodies in culture; however, their meiotic spindles were often positioned parallel, instead of perpendicular, to the ooplasmic membrane at telophase I, and some oocytes extruded the entire spindle with all the chromosomes into the polar bodies at metaphase II, a mechanism we previously reported in Mei1-/- oocytes. The occurrence of a common mechanism in two mouse models argues in favor of its plausibility at the origin of androgenetic HM formation in humans.

Identification of a novel splicing variant of thyroid hormone receptor interaction protein 13 (TRIP13) in female infertility characterized by oocyte maturation arrest

PURPOSE

As a cause of primary female infertility, oocyte maturation arrest (OMA) is characterized by failure to obtain mature oocytes due to abnormal meiosis. We aimed to identify pathogenic variants in two sisters with OMA phenotype from a non-consanguineous family.

METHODS

Whole-exome sequencing and Sanger sequencing were conducted to identify and validate the disease-causing gene variant. Additionally, we performed a minigene assay, quantitative reverse transcription PCR, and Western blotting to assess the effects of the variant.

RESULTS

We identified a novel homozygous splicing variant (c.1021-11T>C) in TRIP13, which followed a recessive inheritance pattern. Minigene assay showed that the variant could disrupt the integrity of TRIP13 mRNA, as evidenced by the production of an alternative transcript with intron10 intermediate retention of 79 bp. Compared to normal controls, the expression of TRIP13 mRNA and abundance of TRIP13 protein were also significantly decreased in Epstein-Barr virus-immortalized lymphoblastoid cells derived from affected individuals.

CONCLUSION

Our findings confirm the contribution of genetic factors to OMA and expand the mutation spectrum of TRIP13 in female infertility.

Featuring BRCA1 and BRCA2 germline mutational landscape from Asturias (North Spain)

The singular BRCA1/2 mutational landscape of Asturias is updated 10 years after the first study. We analyzed BRCA1 and BRCA2 pathogenic variants in 1653 index cases. In total, 238 families were identified to carry a pathogenic variant, 163 families in BRCA1 and 75 families in BRCA2. This yielded a prevalence rate of 14.4%. Seven recurrent variants were found accounting for 55% of the cases. Among them, three are widely distributed (BRCA1 c.211A>G, c.470_471del and c.3331_3334del) and four had been reported as novel in Asturias: two in BRCA1 (c.1674del and c.2901_2902dup) and two in BRCA2 (c.2095C>T and c.4040_4035delinsC). A common haplotype was established for all recurrent variants indicating a shared ancestral origin. Three splicing analyses are shown: BRCA1:c.5152+3A>C and BRCA1:c.5333-3T>G that lead to skipping of exon 18, and 22 respectively, and BRCA1:c.5278-1G>T giving rise to two transcripts, one lacking exon 21 (p.Ille1760Glyfs*60) and one lacking the first 8 nucleotides of exon 21 (p.Phe1761Asnfs*14), supporting pathogenicity for these variants.

A novel PKHD1 splicing variant identified in a fetus with autosomal recessive polycystic kidney disease

Objective: Variants of the polycystic kidney and hepatic disease 1 (PKHD1) gene are associated with autosomal recessive polycystic kidney disease (ARPKD). This study aimed to identify the genetic causes in a Chinese pedigree with ARPKD and design a minigene construct of the PKHD1 gene to investigate the impact of its variants on splicing. Methods: Umbilical cord samples from the proband and peripheral blood samples from his parents were collected, and genomic DNA was extracted for whole-exome sequencing (WES). Bioinformatic analysis was used to identify potential genetic causes, and Sanger sequencing confirmed the existence of variants within the pedigree. A minigene assay was performed to validate the effects of an intronic variant on mRNA splicing. Results: Two variants, c.9455del (p.N3152Tfs*10) and c.2408-13C>G, were identified in the PKHD1 gene (NM_138694.4) by WES; the latter has not been previously reported. In silico analysis predicted that this intronic variant is potentially pathogenic. Bioinformatic splice prediction tools revealed that the variant is likely to strongly impact splice site function. An in vitro minigene assay revealed that c.2408-13C>G can cause aberrant splicing, resulting in the retention of 12 bp of intron 23. Conclusion: A novel pathogenic variant of PKHD1, c.2408-13C>G, was found in a fetus with ARPKD, which enriches the variant spectrum of the PKHD1 gene and provides a basis for genetic counseling and the diagnosis of ARPKD. Minigenes are optimal to determine whether intron variants can cause aberrant splicing.

Systematic analysis of CNGA3 splice variants identifies different mechanisms of aberrant splicing

Achromatopsia is an autosomal recessive cone photoreceptor disease that is frequently caused by pathogenic variants in the CNGA3 gene. Here, we present a systematic functional analysis of 20 CNGA3 splice site variants detected in our large cohort of achromatopsia patients and/or listed in common variant databases. All variants were analyzed by functional splice assays based on the pSPL3 exon trapping vector. We demonstrated that ten variants, both at canonical and non-canonical splice sites, induced aberrant splicing, including intronic nucleotide retention, exonic nucleotide deletion and exon skipping, resulting in 21 different aberrant transcripts. Of these, eleven were predicted to introduce a premature termination codon. The pathogenicity of all variants was assessed based on established guidelines for variant classification. Incorporation of the results of our functional analyses enabled re-classification of 75% of variants previously classified as variants of uncertain significance into either likely benign or likely pathogenic. Our study is the first in which a systematic characterization of putative CNGA3 splice variants has been performed. We demonstrated the utility of pSPL3 based minigene assays in the effective assessment of putative splice variants. Our findings improve the diagnosis of achromatopsia patients, who may thus benefit from future gene-based therapeutic strategies.

Analyzing the effects of BRCA1/2 variants on mRNA splicing by minigene assay

As BRCA1/2 gene sequencing become more extensive, a large number VUS (variants of uncertain significance) emerge rapidly. Verifying the splicing effect is an effective means for VUS reclassification. The Minigene Assay platform was established and its reliability was verified in this article. 47 BRCA1 or BRCA2 variants were selected and performed to validate their effect on mRNA splicing. The results showed that, a total of 16 variants were experimentally proved to have effects on mRNA splicing, among which 14 variants were shown to cause truncated proteins by Sanger sequencing. While the other two variants, BRCA2 c.7976 + 3 A > G and BRCA1 c.5152 + 3_5152 + 4insT was analyzed to cause 57 bp and 26 bp base in-frame deletion, respectively. The remaining 31 variants were not shown to cause mRNA splicing abnormity, including several sites at the edge of exons, which were predicted to affect splicing of mRNA by multiple bioinformatic software. Based on our experimental results, 37 variants were reclassified by ACMG rules. Our study showed that experimental splicing analysis was effectual for variants classification, and multiple functional assay or clinical data were also necessary for comprehensive judgment of variants.

Curating the gnomAD database: Report of novel variants in the thyroid peroxidase gene using in silico bioinformatics algorithms and a literature review

Thyroid peroxidase (TPO) is a membrane-bound glycoprotein located at the apical side of the thyroid follicular cells that catalyzes both iodination and coupling of iodotyrosine residues within the thyroglobulin molecule, leading to the synthesis of thyroid hormone. Variants in TPO cause congenital hypothyroidism (CH) by iodide organification defect and are commonly inherited in an autosomal recessive fashion. In the present work, we report a detailed population analysis and bioinformatic prediction of the TPO variants indexed in the Genome Aggregation Database (gnomAD) v2.1.1. The proportion of missense cysteine variants and nonsense, frameshift, and splice acceptor/donor variants were analyzed in each ethnic group (European (Non-Finnish), European (Finnish), African/African Americans, Latino/Admixed American, East Asian, South Asian, Ashkenazi Jewish, Other). The results showed a clear predominance of frameshift variants in the East Asian (82%) and European (Finnish) (75%) population, whereas the splice site variants predominate in African/African Americans (99.46%), Other (96%), Latino/Admixed American (94%), South Asian (86%), European (Non-Finnish) (56%) and Ashkenazi Jewish (56%) populations. The analysis of the distribution of the variants indexed in gnomAD v2.1.1 database revealed that most missense variants identified in the An peroxidase domain map in exon 8, followed by exons 11, 7 and 9, and finally in descending order by exons 10, 6, 12 and 5. In total, 183 novel TPO variants were described (13 missense cysteine's variants, 158 missense variants involving the An peroxidase domain and 12 splicing acceptor or donor sites variants) which were not reported in the literature and that would have deleterious effects on prediction programs. In the gnomAD v2.1.1 population, the estimated prevalence of heterozygous carriers of the potentially damaging variants was 1:77. In conclusion, we provide an updated and curated reference source of new TPO variants for application in clinical diagnosis and genetic counseling. Also, this work contributes to elucidating the molecular basis of CH associated with TPO defects.

Case Review: Whole-Exome Sequencing Analyses Identify Carriers of a Known Likely Pathogenic Intronic BRCA1 Variant in Ovarian Cancer Cases Clinically Negative for Pathogenic BRCA1 and BRCA2 Variants

Background: Detecting pathogenic intronic variants resulting in aberrant splicing remains a challenge in routine genetic testing. We describe germline whole-exome sequencing (WES) analyses and apply in silico predictive tools of familial ovarian cancer (OC) cases reported clinically negative for pathogenic BRCA1 and BRCA2 variants. Methods: WES data from 27 familial OC cases reported clinically negative for pathogenic BRCA1 and BRCA2 variants and 53 sporadic early-onset OC cases were analyzed for pathogenic variants in BRCA1 or BRCA2. WES data from carriers of pathogenic BRCA1 or BRCA2 variants were analyzed for pathogenic variants in 10 other OC predisposing genes. Loss of heterozygosity analysis was performed on tumor DNA from variant carriers. Results: BRCA1 c.5407-25T>A intronic variant, identified in two affected sisters and one sporadic OC case, is predicted to create a new splice effecting transcription of BRCA1. WES data from BRCA1 c.5407-25T>A carriers showed no evidence of pathogenic variants in other OC predisposing genes. Sequencing the tumor DNA from the variant carrier showed complete loss of the wild-type allele. Conclusions: The findings support BRCA1 c.5407-25T>A as a likely pathogenic variant and highlight the importance of investigating intronic sequences as causal variants in OC families where the involvement of BRCA1 is highly suggestive.

Presumed COL4A3/COL4A4 Missense/Synonymous Variants Induce Aberrant Splicing

Background

The incorrect interpretation of missense and synonymous variants can lead to improper molecular diagnosis and subsequent faulty genetic counselling. The aim of this study was to evaluate the pathogenicity of presumed COL4A3/COL4A4 missense and synonymous variants detected by next-generation sequencing to provide evidence for diagnosis and genetic counselling.

Methods

Patients' clinical findings and genetic data were analysed retrospectively. An in vitro minigene assay was conducted to assess the effect of presumed COL4A3/COL4A4 missense and synonymous variants on RNA splicing.

Results

Five unclassified COL4A3/COL4A4 variants, which were detected in five of 343 patients with hereditary kidney diseases, were analysed. All of them were predicted to affect splicing by Human Splicing Finder. The presumed COL4A3 missense variant c.4793T > G [p. (Leu1598Arg)] resulted in a loss of alternative full-length transcript during the splicing process. The COL4A3 transcript carried synonymous variant c.765G > A [p. (Thr255Thr)], led to an in-frame deletion of exon 13. Nevertheless, variants c.3566G > A [p. (Gly1189Glu)] in COL4A3 and c.3990G > A [p. (Pro1330Pro)], c.4766C > T [p. (Pro1589Leu)] in COL4A4 exhibited no deleterious effect on splicing. Among the five patients harbouring the abovementioned COL4A3/COL4A4 variants, three patients were genetically diagnosed with autosomal recessive Alport syndrome, one patient was highly suspected of having thin basement membrane nephropathy, and the other patient was clinically diagnosed with Alport syndrome.

Conclusions

COL4A3 presumed missense variant p. (Leu1598Arg) and synonymous variant p. (Thr255Thr) affect RNA splicing, which highlights the prime importance of transcript analysis of unclassified exonic sequence variants for better molecular diagnosis and genetic counselling. Meanwhile, the reliability of splicing predictions by predictive tools for exonic substitutions needs to be improved.

Sequencing of BRCA1/2-alterations using NGS-based technology: annotation as a challenge

In this study, the molecular profile of different BRCA-associated tumor types was assessed with regard to the classification and annotation of detected BRCA1/2 variants. The aim was to establish guidelines in order to facilitate the interpretation of BRCA1/2 alterations in routine diagnostics. Annotation of detected variants was evaluated compared to background mutations found in normal tissue samples and manually reviewed according to distinct online databases. This retrospective study included 48 samples (45 tumors, three non-tumors), which were sequenced with the GeneReader (QIAGEN). Thereof ten samples were additionally analyzed with the Ion S5™ (Thermo Fisher) and 20 samples with the MiSeq™ (Illumina^®) to compare the different NGS devices, as well as the sequencing results and their quality. The analysis showed that the individual NGS platforms detected different numbers of BRCA1/2 alterations in the respective tumor sample. In addition, the GeneReader revealed variability in the detection and classification of pathogenic alterations within the platform itself as well as in comparison with the other platforms or online databases. The study concluded that the Ion S5™ in combination with the Oncomine™ Comprehensive Assay v3 is most recommendable for current and prospective requirements of molecular analysis in routine diagnostics. In addition to the two BRCA1/2 genes, a broad number of other genes (BRCAness genes and genes involved in the repair pathway) is covered by the panel, which may open up new treatment options for patients depending on the respective eligibility criteria.

Homologous Recombination Deficiency Assays in Epithelial Ovarian Cancer: Current Status and Future Direction

Epithelial ovarian cancer (EOC) patients are generally diagnosed at an advanced stage, usually relapse after initial treatments, which include debulking surgery and adjuvant platinum-based chemotherapy, and eventually have poor 5-year survival of less than 50%. In recent years, promising survival benefits from maintenance therapy with poly(ADP-ribose) polymerase (PARP) inhibitor (PARPi) has changed the management of EOC in newly diagnosed and recurrent disease. Identification of BRCA mutations and/or homologous recombination deficiency (HRD) is critical for selecting patients for PARPi treatment. However, the currently available HRD assays are not perfect predictors of the clinical response to PARPis in EOC patients. In this review, we introduce the concept of synthetic lethality, the rationale of using PARPi when HRD is present in tumor cells, the clinical trials of PARPi incorporating the HRD assays for EOC, the current HRD assays, and other HRD assays in development.

The role of patient-derived ovarian cancer organoids in the study of PARP inhibitors sensitivity and resistance: from genomic analysis to functional testing

Epithelial ovarian cancer (EOC) harbors distinct genetic features such as homologous recombination repair (HRR) deficiency, and therefore may respond to poly ADP-ribose polymerase inhibitors (PARPi). Over the past few years, PARPi have been added to the standard of care for EOC patients in both front-line and recurrent settings. Next-generation sequencing (NGS) genomic analysis provides key information, allowing for the prediction of PARPi response in patients who are PARPi naïve. However, there are indeed some limitations in NGS analyses. A subset of patients can benefit from PARPi, despite the failed detection of the predictive biomarkers such as BRCA1/2 mutations or HRR deficiency. Moreover, in the recurrent setting, the sequencing of initial tumor does not allow for the detection of reversions or secondary mutations restoring proficient HRR and thus leading to PARPi resistance. Therefore, it becomes crucial to better screen patients who will likely benefit from PARPi treatment, especially those with prior receipt of maintenance PARPi therapy. Recently, patient-derived organoids (PDOs) have been regarded as a reliable preclinical platform with clonal heterogeneity and genetic features of original tumors. PDOs are found feasible for functional testing and interrogation of biomarkers for predicting response to PARPi in EOC. Hence, we review the strengths and limitations of various predictive biomarkers and highlight the role of patient-derived ovarian cancer organoids as functional assays in the study of PARPi response. It was found that a combination of NGS and functional assays using PDOs could enhance the efficient screening of EOC patients suitable for PARPi, thus prolonging their survival time.

A DNA Damage Response Gene Panel for Different Histologic Types of Epithelial Ovarian Carcinomas and Their Outcomes

DNA damage response (DDR) is important for maintaining genomic integrity of the cell. Aberrant DDR pathways lead to accumulation of DNA damage, genomic instability and malignant transformations. Gene mutations have been proven to be associated with epithelial ovarian cancer, and the majority of the literature has focused on BRCA. In this study, we investigated the somatic mutation of DNA damage response genes in epithelial ovarian cancer patients using a multiple-gene panel with next-generation sequencing. In all, 69 serous, 39 endometrioid and 64 clear cell carcinoma patients were enrolled. Serous carcinoma patients (69.6%) had higher percentages of DDR gene mutations compared with patients with endometrioid (33.3%) and clear cell carcinoma (26.6%) (p < 0.001, chi-squared test). The percentages of DDR gene mutations in patients with recurrence (53.9 vs. 32.9% p = 0.006, chi-squared test) or cancer-related death (59.2 vs. 34.4% p = 0.001, chi-squared test) were higher than those without recurrence or living patients. In endometrioid carcinoma, patients with ≥2 DDR gene mutations had shorter PFS (p = 0.0035, log-rank test) and OS (p = 0.015, log-rank test) than those with one mutation or none. In clear cell carcinoma, patients with ≥2 DDR gene mutations had significantly shorter PFS (p = 0.0056, log-rank test) and OS (p = 0.0046, log-rank test) than those with 1 DDR mutation or none. In the EOC patients, somatic DDR gene mutations were associated with advanced-stage tumor recurrence and tumor-related death. Type I EOC patients with DDR mutations had an unfavorable prognosis, especially for clear cell carcinoma.

Curating the gnomAD database: Report of novel variants in the thyrogobulin gene using in silico bioinformatics algorithms

Thyroglobulin (TG) is a large glycosylated protein of 2767 amino acids, secreted by the thyrocytes into the follicular lumen. It plays an essential role in the process of thyroid hormone synthesis. TG gene variants lead to permanent congenital hypothyroidism. In the present work, we report a detailed population and bioinformatic prediction analyses of the TG variants indexed in the Genome Aggregation Database (gnomAD). The results showed a clear predominance of nonsense variants in the European (Finnish), European (Non-Finnish) and Ashkenazi Jewish ethnic groups, whereas the splice site variants predominate in South Asian and African/African-American populations. In total, 282 novel TG variants were described (47 missense involving the wild-type cysteine residues, 177 missense located in the ChEL domain and 58 splice site variants) which were not reported in the literature and that would have deleterious effects in prediction programs. In the gnomAD population, the estimated prevalence of heterozygous carriers of the potentially damaging variants was 1:320. In conclusion, we provide an updated and curated reference source for the diagnosis of thyroid disease, mainly to congenital hypothyroidism due to TG deficiency. The identification and characterization of TG variants is undoubtedly a valuable approach to study the TG structure/function relations and an important tool for clinical diagnosis and genetic counseling.

A homozygous variant in TBPL2 was identified in women with oocyte maturation defects and infertility

STUDY QUESTION

What are the genetic causes of oocyte maturation defects?

SUMMARY ANSWER

A homozygous splicing variant (c.788 + 3A>G) in TATA-box binding protein like 2 (TBPL2) was identified as a contributory genetic factor in oocyte maturation defects.

WHAT IS KNOWN ALREADY

TBPL2, a vertebrate oocyte-specific general transcription factor, is essential for oocyte development. TBPL2 variants have not been studied in human oocyte maturation defects.

STUDY DESIGN, SIZE, DURATION

Two infertile families characterized by oocyte maturation defects were recruited for whole-exome sequencing (WES).

PARTICIPANTS/MATERIALS, SETTING, METHODS

Genomic DNA was extracted from peripheral blood for WES analysis. Sanger sequencing was performed for data validation. Pathogenicity of variants was predicted by in silico analysis. Minigene assay and single-oocyte RNA sequencing were performed to investigate the effects of the variant on mRNA integrity and oocyte transcriptome, respectively.

MAIN RESULTS AND THE ROLE OF CHANCE

A homozygous splicing variant (c.788 + 3A>G) in TBPL2 was identified in two unrelated families characterized by oocyte maturation defects. Haplotype analysis indicated that the disease allele of Families 1 and 2 was independent. The variant disrupted the integrity of TBPL2 mRNA. Transcriptome sequencing of affected oocytes showed that vital genes for oocyte maturation and fertilization were widely and markedly downregulated, suggesting that a mutation in the transcriptional factor, TBPL2, led to global gene alterations in oocytes.

LIMITATIONS, REASONS FOR CAUTION

Limitations include the lack of direct functional evidence. Owing to the scarcity of human oocyte samples, only two immature MI oocytes were obtained from the patients, and we could only investigate the effect of the mutation at the transcriptional level by high-throughput sequencing technology. No extra oocytes were obtained to assess the transcriptional activity of the mutant oocytes by immunofluorescence, or investigate the effects on the binding of TBPL2 caused by the mutation.

WIDER IMPLICATIONS OF THE FINDINGS

Our findings highlight a critical role of TBPL2 in female reproduction and identify a homozygous splicing mutation in TBPL2 that might be related to defects in human oocyte maturation. This information will facilitate the genetic diagnosis of infertile individuals with repeated failures of IVF, providing a basis for genetic counseling.

STUDY FUNDING/COMPETING INTEREST(S)

This study was supported by the National Key Research and Development Program of China (2018YFC1004000, 2017YFC1001504 and 2017YFC1001600), the National Natural Science Foundation of China (81871168, 31900409 and 31871509), the Foundation for Distinguished Young Scholars of Shandong Province (JQ201816), the Innovative Research Team of High-Level Local Universities in Shanghai (SSMU-ZLCX20180401) and the Fundamental Research Funds of Shandong University. The authors have no competing interests to declare.

TRIAL REGISTRATION NUMBER

N/A.

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.

SEPTIN12 c.474 G > A polymorphism as a risk factor in teratozoospermic patients

Teratozoospermia is a condition related to poor morphologically normal sperm count below the lower reference limit, which could hinder natural conception. Single nucleotide polymorphisms (SNPs) in the genes involved in sperm production and testicular function are proved to be risk factors, resulting in decreased sperm parameters and defects in sperm morphology. c.474 G > A polymorphism in the SEPTIN12 gene which is one of the testis-specific genes creates a novel splice variant and the resulting truncated protein was previously found to be more prevalent in infertile men. We aimed to investigate the association of SEPTIN12 c.474 G > A polymorphism with male infertility in teratozoospermia patients. Forty-eight teratozoospermic patients, diagnosed according to Kruger's criteria and 164 fertile controls who fathered at least 1 child within 3 years without assisted reproductive technologies were included into our prospective randomized controlled study. PCR-RFLP method was used for genotyping. Although no statistical difference was found between teratozoospermic patients and fertile controls in terms of genotype distributions, significance was identified between the genotypes of all and non-smoking teratozoopermic patients in terms of neck defects. SEPTIN12 c.474 G > A polymorphism was shown to be associated with sperm neck defects in teratozoospermic patients using the dominant statistical model. Smoking was identified as a risk factor for the sperm morphology defects in teratozoospermic A allele carriers.

LAMB2 novel variant c.2885-9 C>A affects RNA splicing in a minigene assay

Background

Both Pierson syndrome (PS) and isolated nephrotic syndrome can be caused by LAMB2 biallelic pathogenic variants. Only 15 causative splicing variants in the LAMB2 gene have been reported. However, the pathogenicity of most of these variants has not been verified, which may lead to incorrect interpretation of the functional consequence of these variants.

Methods

Using high‐throughput DNA sequencing and Sanger sequencing, we detected variants in a female with clinically suspected PS. A minigene splicing assay was performed to assess the effect of LAMB2 intron 20 c.2885‐9C>A on RNA splicing. We also performed the immunohistochemical analysis of laminin beta‐2 in kidney tissues.

Results

Two novel LAMB2 heteroallelic variants were found: a paternally inherited variant c.2885‐9C>A in intron 20 and a maternally inherited variant c. 3658C>T (p. (Gln1220Ter)). In vitro minigene assay showed that the variant c.2885‐9C>A caused erroneous integration of a 7 bp sequence into intron 20. Immunohistochemical analysis revealed the absence of glomerular expression of laminin beta‐2, the protein encoded by LAMB2.

Conclusion

We demonstrated the impact of a novel LAMB2 intronic variant on RNA splicing using the minigene assay firstly. Our results extend the mutational spectrum of LAMB2.

Analysis of BRCA1/2 variants of unknown significance in the prospective Korean Hereditary Breast Cancer study

Genetic testing for BRCA1 and BRCA2 is crucial in diagnosing hereditary breast and ovarian cancer syndromes and has increased with the development of multigene panel tests. However, results classified as variants of uncertain significance (VUS) present challenges to clinicians in attempting to choose an appropriate management plans. We reviewed a total of 676 breast cancer patients included in the Korean Hereditary Breast Cancer (KOHBRA) study with a VUS on BRCA mutation tests between November 2007 and April 2013. These results were compared to the ClinVar database. We calculated the incidence and odds ratios for these variants using the Korean Reference Genome Database. A total of 58 and 91 distinct VUS in BRCA1 and BRCA2 were identified in the KOHBRA study (comprising 278 and 453 patients, respectively). A total of 27 variants in the KOHBRA study were not registered in the Single Nucleotide Polymorphism database. Among BRCA1 VUSs, 20 were reclassified as benign or likely benign, four were reclassified as pathogenic or likely pathogenic, and eight remained as VUSs according to the ClinVar database. Of the BRCA2 VUSs, 25 were reclassified as benign or likely benign, two were reclassified as pathogenic or likely pathogenic, and 33 remained as VUS according to the ClinVar database. There were 12 variants with conflicting interpretations of pathogenicity for BRCA1 and 18 for BRCA2. Among them, p.Leu1780Pro showed a particularly high odds ratio. Six pathogenic variants and one conflicting variant identified using ClinVar could be reclassified as pathogenic variants in this study. Using updated ClinVar information and calculating odds ratios can be helpful when reclassifying VUSs in BRCA1/2.

The BRCA1 c.788G > T (NM_007294.4) variant in a high grade serous ovarian cancer (HGSOC) patient: foods for thought

In this report we described the case of a BRCA1/2 (BRCA) molecular testing performed on tumor sample in a High Grade Serous Ovarian Cancer (HGSOC) patient with two different Next Generation Tumor Sequencing (NGTS) pipelines. The two clinical reports leaded to apparently different BRCA status, providing important foods for thought. After NGTS, the gene sequencing information (i.e., reads) are aligned to the reference gene sequences obtained from public databases, in order to provide an uniform nomenclature for unambiguous variant designation. However, the criteria adopted for variant reporting in tissue test are not always univocal. Particularly, this is the case of rare and unclassified BRCA variants for which the molecular evaluation may be a relevant challenge. Here we described a BRCA1 unclassified variant that may be re-evaluated in the context of alternative BRCA1 transcripts due to its different biological effect. We underlined that an in-depth knowledge of BRCA testing is mandatory for its appropriate use.

A novel mutation in intron 11 donor splice site, responsible of a rare genotype in thyroglobulin gene by altering the pre-mRNA splincing process. Cell expression and bioinformatic analysis

Thyroglobulin (TG) is a homodimeric glycoprotein synthesized by the thyroid gland. To date, two hundred twenty-seven variations of the TG gene have been identified in humans. Thyroid dyshormonogenesis due to TG gene mutations have an estimated incidence of approximately 1 in 100,000 newborns. The clinical spectrum ranges from euthyroid to mild or severe hypothyroidism. The purpose of the present study was to identify and characterize new variants in the TG gene. We report an Argentine patient with congenital hypothyroidism, enlarged thyroid gland and low levels of serum TG. Sequencing of DNA, expression of chimeric minigenes as well as bioinformatics analysis were performed. DNA sequencing identified the presence of compound heterozygous mutations in the TG gene: the maternal mutation consists of a c.3001+5G > A, whereas the paternal mutation consists of p.Arg296*. Minigen analysis of the variant c.3001+5A performed in HeLa, CV1 and Hek293T cell lines, showed a total lack of transcript expression. So, in order to validate that the loss of expression was caused by such variation, site-directed mutagenesis was performed on the mutated clone, which previously had a pSPL3 vector change, to give rise to a wild-type clone c.3001+5G, endorsing that the mutation c.3001+5G > A is the cause of the total lack of expression. In conclusion, we demonstrate that the c.3001+5G > A mutation causes a rare genotype, altering the splicing of the pre-mRNA. This work contributes to elucidating the molecular bases of TG defects associated with congenital hypothyroidism and expands our knowledge in relation to the pathologic roles of the position 5 in the donor splice site.

Novel candidates of pathogenic variants of the BRCA1 and BRCA2 genes from a dataset of 3,552 Japanese whole genomes (3.5KJPNv2)

Identification of the population frequencies of definitely pathogenic germline variants in two major hereditary breast and ovarian cancer syndrome (HBOC) genes, BRCA1/2, is essential to estimate the number of HBOC patients. In addition, the identification of moderately penetrant HBOC gene variants that contribute to increasing the risk of breast and ovarian cancers in a population is critical to establish personalized health care. A prospective cohort subjected to genome analysis can provide both sets of information. Computational scoring and prospective cohort studies may help to identify such likely pathogenic variants in the general population. We annotated the variants in the BRCA1 and BRCA2 genes from a dataset of 3,552 whole-genome sequences obtained from members of a prospective cohorts with genome data in the Tohoku Medical Megabank Project (TMM) with InterVar software. Computational impact scores (CADD_phred and Eigen_raw) and minor allele frequencies (MAFs) of pathogenic (P) and likely pathogenic (LP) variants in ClinVar were used for filtration criteria. Familial predispositions to cancers among the 35,000 TMM genome cohort participants were analyzed to verify the identified pathogenicity. Seven potentially pathogenic variants were newly identified. The sisters of carriers of these moderately deleterious variants and definite P and LP variants among members of the TMM prospective cohort showed a statistically significant preponderance for cancer onset, from the self-reported cancer history. Filtering by computational scoring and MAF is useful to identify potentially pathogenic variants in BRCA genes in the Japanese population. These results should help to follow up the carriers of variants of uncertain significance in the HBOC genes in the longitudinal prospective cohort study.

ESMO recommendations on predictive biomarker testing for homologous recombination deficiency and PARP inhibitor benefit in ovarian cancer

BACKGROUND

Homologous recombination repair deficiency (HRD) is a frequent feature of high-grade serous ovarian, fallopian tube and peritoneal carcinoma (HGSC) and is associated with sensitivity to PARP inhibitor (PARPi) therapy. HRD testing provides an opportunity to optimise PARPi use in HGSC but methodologies are diverse and clinical application remains controversial.

MATERIALS AND METHODS

To define best practice for HRD testing in HGSC the ESMO Translational Research and Precision Medicine Working Group launched a collaborative project that incorporated a systematic review approach. The main aims were to (i) define the term 'HRD test'; (ii) provide an overview of the biological rationale and the level of evidence supporting currently available HRD tests; (iii) provide recommendations on the clinical utility of HRD tests in clinical management of HGSC.

RESULTS

A broad range of repair genes, genomic scars, mutational signatures and functional assays are associated with a history of HRD. Currently, the clinical validity of HRD tests in ovarian cancer is best assessed, not in terms of biological HRD status per se, but in terms of PARPi benefit. Clinical trials evidence supports the use of BRCA mutation testing and two commercially available assays that also incorporate genomic instability for identifying subgroups of HGSCs that derive different magnitudes of benefit from PARPi therapy, albeit with some variation by clinical scenario. These tests can be used to inform treatment selection and scheduling but their use is limited by a failure to consistently identify a subgroup of patients who derive no benefit from PARPis in most studies. Existing tests lack negative predictive value and inadequately address the complex and dynamic nature of the HRD phenotype.

CONCLUSIONS

Currently available HRD tests are useful for predicting likely magnitude of benefit from PARPis but better biomarkers are urgently needed to better identify current homologous recombination proficiency status and stratify HGSC management.

Splicing impact of deep exonic missense variants in CAPN3 explored systematically by minigene functional assay

Improving the accuracy of variant interpretation during diagnostic sequencing is a major goal for genomic medicine. To explore an often-overlooked splicing effect of missense variants, we developed the functional assay ("minigene") for the majority of exons of CAPN3, the gene responsible for limb girdle muscular dystrophy. By systematically screening 21 missense variants distributed along the gene, we found that eight clinically relevant missense variants located at a certain distance from the exon-intron borders (deep exonic missense variants) disrupted normal splicing of CAPN3 exons. Several recent machine learning-based computational tools failed to predict splicing impact for the majority of these deep exonic missense variants, highlighting the importance of including variants of this type in the training sets during the future algorithm development. Overall, 24 variants in CAPN3 gene were explored, leading to the change in the American College of Medical Genetics and Genomics classification of seven of them when results of the "minigene" functional assay were considered. Our findings reveal previously unknown splicing impact of several clinically important variants in CAPN3 and draw attention to the existence of deep exonic variants with a disruptive effect on gene splicing that could be overlooked by the current approaches in clinical genetics.

Calibration of Pathogenicity Due to Variant-Induced Leaky Splicing Defects by Using BRCA2 Exon 3 as a Model System

BRCA2 is a clinically actionable gene implicated in breast and ovarian cancer predisposition that has become a high priority target for improving the classification of variants of unknown significance (VUS). Among all BRCA2 VUS, those causing partial/leaky splicing defects are the most challenging to classify because the minimal level of full-length (FL) transcripts required for normal function remains to be established. Here, we explored BRCA2 exon 3 (BRCA2e3) as a model for calibrating variant-induced spliceogenicity and estimating thresholds for BRCA2 haploinsufficiency. In silico predictions, minigene splicing assays, patients' RNA analyses, a mouse embryonic stem cell (mESC) complementation assay and retrieval of patient-related information were combined to determine the minimal requirement of FL BRCA2 transcripts. Of 100 BRCA2e3 variants tested in the minigene assay, 64 were found to be spliceogenic, causing mild to severe RNA defects. Splicing defects were also confirmed in patients' RNA when available. Analysis of a neutral leaky variant (c.231T>G) showed that a reduction of approximately 60% of FL BRCA2 transcripts from a mutant allele does not cause any increase in cancer risk. Moreover, data obtained from mESCs suggest that variants causing a decline in FL BRCA2 with approximately 30% of wild-type are not pathogenic, given that mESCs are fully viable and resistant to DNA-damaging agents in those conditions. In contrast, mESCs producing lower relative amounts of FL BRCA2 exhibited either null or hypomorphic phenotypes. Overall, our findings are likely to have broader implications on the interpretation of BRCA2 variants affecting the splicing pattern of other essential exons. SIGNIFICANCE: These findings demonstrate that BRCA2 tumor suppressor function tolerates substantial reduction in full-length transcripts, helping to determine the pathogenicity of BRCA2 leaky splicing variants, some of which may not increase cancer risk.

In or Out? New Insights on Exon Recognition through Splice-Site Interdependency

Noncanonical splice-site mutations are an important cause of inherited diseases. Based on in vitro and stem-cell-based studies, some splice-site variants show a stronger splice defect than expected based on their predicted effects, suggesting that other sequence motifs influence the outcome. We investigated whether splice defects due to human-inherited-disease-associated variants in noncanonical splice-site sequences in ABCA4, DMD, and TMC1 could be rescued by strengthening the splice site on the other side of the exon. Noncanonical 5′- and 3′-splice-site variants were selected. Rescue variants were introduced based on an increase in predicted splice-site strength, and the effects of these variants were analyzed using in vitro splice assays in HEK293T cells. Exon skipping due to five variants in noncanonical splice sites of exons in ABCA4, DMD, and TMC1 could be partially or completely rescued by increasing the predicted strengths of the other splice site of the same exon. We named this mechanism “splicing interdependency”, and it is likely based on exon recognition by splicing machinery. Awareness of this interdependency is of importance in the classification of noncanonical splice-site variants associated with disease and may open new opportunities for treatments.

Splicing profile by capture RNA-seq identifies pathogenic germline variants in tumor suppressor genes

Germline variants in tumor suppressor genes (TSGs) can result in RNA mis-splicing and predisposition to cancer. However, identification of variants that impact splicing remains a challenge, contributing to a substantial proportion of patients with suspected hereditary cancer syndromes remaining without a molecular diagnosis. To address this, we used capture RNA-sequencing (RNA-seq) to generate a splicing profile of 18 TSGs (APC, ATM, BRCA1, BRCA2, BRIP1, CDH1, CHEK2, MLH1, MSH2, MSH6, MUTYH, NF1, PALB2, PMS2, PTEN, RAD51C, RAD51D, and TP53) in 345 whole-blood samples from healthy donors. We subsequently demonstrated that this approach can detect mis-splicing by comparing splicing profiles from the control dataset to profiles generated from whole blood of individuals previously identified with pathogenic germline splicing variants in these genes. To assess the utility of our TSG splicing profile to prospectively identify pathogenic splicing variants, we performed concurrent capture DNA and RNA-seq in a cohort of 1000 patients with suspected hereditary cancer syndromes. This approach improved the diagnostic yield in this cohort, resulting in a 9.1% relative increase in the detection of pathogenic variants, demonstrating the utility of performing simultaneous DNA and RNA genetic testing in a clinical context.

Expert consensus guidelines for the genetic diagnosis of Alport syndrome

Recent expert guidelines recommend genetic testing for the diagnosis of Alport syndrome. Here, we describe current best practice and likely future developments. In individuals with suspected Alport syndrome, all three COL4A5, COL4A3 and COL4A4 genes should be examined for pathogenic variants, probably by high throughput-targeted next generation sequencing (NGS) technologies, with a customised panel for simultaneous testing of the three Alport genes. These techniques identify up to 95% of pathogenic COL4A variants. Where causative pathogenic variants cannot be demonstrated, the DNA should be examined for deletions or insertions by re-examining the NGS sequencing data or with multiplex ligation-dependent probe amplification (MLPA). These techniques identify a further 5% of variants, and the remaining few changes include deep intronic splicing variants or cases of somatic mosaicism. Where no pathogenic variants are found, the basis for the clinical diagnosis should be reviewed. Genes in which mutations produce similar clinical features to Alport syndrome (resulting in focal and segmental glomerulosclerosis, complement pathway disorders, MYH9-related disorders, etc.) should be examined. NGS approaches have identified novel combinations of pathogenic variants in Alport syndrome. Two variants, with one in COL4A3 and another in COL4A4, produce a more severe phenotype than an uncomplicated heterozygous change. NGS may also identify further coincidental pathogenic variants in genes for podocyte-expressed proteins that also modify the phenotype. Our understanding of the genetics of Alport syndrome is evolving rapidly, and both genetic and non-genetic factors are likely to contribute to the observed phenotypic variability.

Functional Characterization of Argininosuccinate Lyase Gene Variants by Mini-Gene Splicing Assay

Objective

Argininosuccinate lyase (ASL) gene mutations account for argininosuccinic aciduria (ASA). This study aimed to design a minigene construct of ASL gene in order to investigate the impact of variants on splicing.

Methods

The peripheral blood samples were collected from the family members, and genomic DNA was extracted for gene diagnosis using the total exon sequencing method. The novel mutation gene was cloned into pEGFP-C1 vector, and the pathogenicity of the mutation was examined in cultured cells in vitro.

Results

The clinical diagnosis of the proband as ASA was clear. Two pathogenic mutations, c.281G>T (p.Arg94Leu) and c.208-15 T>A were detected in the ASL gene, and the two mutations had not been reported. The minigene expression in vitro confirmed that c.208-15 T>A could cause aberrant splicing, resulting in the retention of 13 bp in intron 2 to exon 3.

Conclusion

Two new pathogenic mutations of ASL gene, c.208-15 T>A and c.281G>T, were found in an ASA family, which enriches the mutational profile of the ASL gene and provides a basis for genetic diagnosis of ASA. Minigenes are optimal approaches to determine whether the intron mutation can cause aberrant splicing.

BRCA1-No Matter How You Splice It

BRCA1 (breast cancer 1, early onset), a well-known breast cancer susceptibility gene, is a highly alternatively spliced gene. BRCA1 alternative splicing may serve as an alternative regulatory mechanism for the inactivation of the BRCA1 gene in both hereditary and sporadic breast cancers, and other BRCA1-associated cancers. The alternative transcripts of BRCA1 can mimic known functions, possess unique functions compared with the full-length BRCA1 transcript, and in some cases, appear to function in opposition to full-length BRCA1 In this review, we will summarize the functional "naturally occurring" alternative splicing transcripts of BRCA1 and then discuss the latest next-generation sequencing-based detection methods and techniques to detect alternative BRCA1 splicing patterns and their potential use in cancer diagnosis, prognosis, and therapy.

Mis-splicing in breast cancer: identification of pathogenic BRCA2 variants by systematic minigene assays

Splicing disruption is a common mechanism of gene inactivation associated with germline variants of susceptibility genes. To study the role of BRCA2 mis-splicing in hereditary breast/ovarian cancer (HBOC), we performed a comprehensive analysis of variants from BRCA2 exons 2-9, as well as the initial characterization of the regulatory mechanisms of such exons. A pSAD-based minigene with exons 2-9 was constructed and validated in MCF-7 cells, producing the expected transcript (1016-nt/V1-BRCA2_exons_2-9-V2). DNA variants from mutational databases were analyzed by NNSplice and Human Splicing Finder softwares. To refine ESE-variant prediction, we mapped the regulatory regions through a functional strategy whereby 26 exonic microdeletions were introduced into the minigene and tested in MCF-7 cells. Thus, we identified nine spliceogenic ESE-rich intervals where ESE-variants may be located. Combining bioinformatics and microdeletion assays, 83 variants were selected and genetically engineered in the minigene. Fifty-three changes impaired splicing: 28 variants disrupted the canonical sites, four created new ones, 10 abrogated enhancers, eight created silencers and three caused a double-effect. Notably, nine spliceogenic-ESE variants were located within ESE-containing intervals. Capillary electrophoresis and sequencing revealed more than 23 aberrant transcripts, where exon skipping was the most common event. Interestingly, variant c.67G>A triggered the usage of a noncanonical GC-donor 4-nt upstream. Thirty-six variants that induced severe anomalies (>60% aberrant transcripts) were analyzed according to the ACMG guidelines. Thus, 28 variants were classified as pathogenic, five as likely pathogenic and three as variants of uncertain significance. Interestingly, 13 VUS were reclassified as pathogenic or likely pathogenic variants. In conclusion, a large fraction of BRCA2 variants (∼64%) provoked splicing anomalies lending further support to the high prevalence of this disease-mechanism. The low accuracy of ESE-prediction algorithms may be circumvented by functional ESE-mapping that represents an optimal strategy to identify spliceogenic ESE-variants. Finally, systematic functional assays by minigenes depict a valuable tool for the initial characterization of splicing anomalies and the clinical interpretation of variants. © 2019 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Usefulness and Limitations of Comprehensive Characterization of mRNA Splicing Profiles in the Definition of the Clinical Relevance of BRCA1/2 Variants of Uncertain Significance

Highly penetrant variants of BRCA1/2 genes are involved in hereditary predisposition to breast and ovarian cancer. The detection of pathogenic BRCA variants has a considerable clinical impact, allowing appropriate cancer-risk management. However, a major drawback is represented by the identification of variants of uncertain significance (VUS). Many VUS potentially affect mRNA splicing, making transcript analysis an essential step for the definition of their pathogenicity. Here, we characterize the impact on splicing of ten BRCA1/2 variants. Aberrant splicing patterns were demonstrated for eight variants whose alternative transcripts were fully characterized. Different events were observed, including exon skipping, intron retention, and usage of de novo and cryptic splice sites. Transcripts with premature stop codons or in-frame loss of functionally important residues were generated. Partial/complete splicing effect and quantitative contribution of different isoforms were assessed, leading to variant classification according to Evidence-based Network for the Interpretation of Mutant Alleles (ENIGMA) consortium guidelines. Two variants could be classified as pathogenic and two as likely benign, while due to a partial splicing effect, six variants remained of uncertain significance. The association with an undefined tumor risk justifies caution in recommending aggressive risk-reduction treatments, but prevents the possibility of receiving personalized therapies with potential beneficial effect. This indicates the need for applying additional approaches for the analysis of variants resistant to classification by gene transcript analyses.

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.

Pathogenicity and selective constraint on variation near splice sites

Mutations that perturb normal pre-mRNA splicing are significant contributors to human disease. We used exome sequencing data from 7833 probands with developmental disorders (DDs) and their unaffected parents, as well as more than 60,000 aggregated exomes from the Exome Aggregation Consortium, to investigate selection around the splice sites and quantify the contribution of splicing mutations to DDs. Patterns of purifying selection, a deficit of variants in highly constrained genes in healthy subjects, and excess de novo mutations in patients highlighted particular positions within and around the consensus splice site of greater functional relevance. By using mutational burden analyses in this large cohort of proband–parent trios, we could estimate in an unbiased manner the relative contributions of mutations at canonical dinucleotides (73%) and flanking noncanonical positions (27%), and calculate the positive predictive value of pathogenicity for different classes of mutations. We identified 18 patients with likely diagnostic de novo mutations in dominant DD-associated genes at noncanonical positions in splice sites. We estimate 35%–40% of pathogenic variants in noncanonical splice site positions are missing from public databases.

A novel LRAT mutation affecting splicing in a family with early onset retinitis pigmentosa

Background and purpose

Retinitis pigmentosa is an important cause of severe visual dysfunction. This study reports a novel splicing mutation in the lecithin retinol acyltransferase (LRAT) gene associated with early onset retinitis pigmentosa and characterizes the effects of this mutation on mRNA splicing and structure.

Methods

Genome-wide linkage analysis followed by dideoxy sequencing of the linked candidate gene LRAT was performed in a consanguineous Pakistani family with autosomal recessive retinitis pigmentosa. In silico prediction and minigene assays were used to investigate the effects of the presumptive splicing mutation.

Results

ARRP in this family was linked to chromosome 4q31.21-q32.1 with a maximum LOD score of 5.40. A novel homozygous intronic mutation (NM_004744.4: c.541-15T>G) was detected in LRAT. In silico tools predicted that the AG-creating mutation would activate an intronic cryptic acceptor site, but cloning fragments of wild-type and mutant sequences of LRAT into Exontrap Cloning Vector pET01 and Expression Cloning Vector pCMV-(DYKD₄K)-C showed that the primary effect of the sequence change was to weaken the nearby authentic acceptor site and cause exon skipping, with only a small fraction of transcripts utilizing the acceptor site producing the reference transcript.

Conclusions

The c.541-15T>G mutation in LRAT results in aberrant splicing and is therefore predicted to be causal for the early onset retinitis pigmentosa in this family. In addition, this work suggests that minigenes adapted to the specific gene and exon may need to be designed for variants in the first and last exon and intron to mimic the authentic splicing mechanism in vivo.

Electronic supplementary material

The online version of this article (10.1186/s40246-018-0165-3) contains supplementary material, which is available to authorized users.

Hereditary breast and ovarian cancer in Andalusian families: a genetic population study

Background

The BRCA1/2 mutation profile varies in Spain according to the geographical area studied. The mutational profile of BRCA1/2 in families at risk for hereditary breast and ovarian cancer has not so far been reported in Andalusia (southern Spain).

Methods

We analysed BRCA1/2 germline mutations in 562 high-risk cases with breast and/or ovarian cancer from Andalusian families from 2010 to 2015.

Results

Among the 562 cases, 120 (21.4%) carried a germline pathogenic mutation in BRCA1/2; 50 in BRCA1 (41.7%) and 70 in BRCA2 (58.3%). We detected 67 distinct mutations (29 in BRCA1 and 38 in BRCA2), of which 3 in BRCA1 (c.845C > A, c.1222_1223delAC, c.2527delA) and 5 in BRCA2 (c.293 T > G, c.5558_5559delGT, c.6034delT, c.6650_6654delAAGAT, c.6652delG) had not been previously described. The most frequent mutations in BRCA1 were c.5078_5080delCTG (10%) and c.5123C > A (10%), and in BRCA2 they were c.9018C > A (14%) and c.5720_5723delCTCT (8%). We identified 5 variants of unknown significance (VUS), all in BRCA2 (c.5836 T > C, c.6323G > T, c.9501 + 3A > T, c.8022_8030delGATAATGGA, c.10186A > C). We detected 76 polymorphisms (31 in BRCA1, 45 in BRCA2) not associated with breast cancer risk.

Conclusions

This is the first study reporting the mutational profile of BRCA1/2 in Andalusia. We identified 21.4% of patients harbouring BRCA1/2 mutations, 58.3% of them in BRCA2. We also characterized the clinical data, mutational profile, VUS and haplotype profile.

Electronic supplementary material

The online version of this article (10.1186/s12885-018-4537-9) contains supplementary material, which is available to authorized users.

Identification of Eight Spliceogenic Variants in BRCA2 Exon 16 by Minigene Assays

Genetic testing of BRCA1 and BRCA2 identifies a large number of variants of uncertain clinical significance whose functional and clinical interpretations pose a challenge for genetic counseling. Interestingly, a relevant fraction of DNA variants can disrupt the splicing process in cancer susceptibility genes. We have tested more than 200 variants throughout 19 BRCA2 exons mostly by minigene assays, 54% of which displayed aberrant splicing, thus confirming the utility of this assay to check genetic variants in the absence of patient RNA. Our goal was to investigate BRCA2 exon 16 with a view to characterizing spliceogenic variants recorded at the mutational databases. Seventy-two different BIC and UMD variants were analyzed with NNSplice and Human Splicing Finder, 12 of which were selected because they were predicted to disrupt essential splice motifs: canonical splice sites (ss; eight variants) and exonic/intronic splicing enhancers (four variants). These 12 candidate variants were introduced into the BRCA2 minigene with seven exons (14–20) by site-directed mutagenesis and then transfected into MCF-7 cells. Seven variants (six intronic and one missense) induced complete abnormal splicing patterns: c.7618-2A>T, c.7618-2A>G, c.7618-1G>C, c.7618-1G>A, c.7805G>C, c.7805+1G>A, and c.7805+3A>C, as well as a partial anomalous outcome by c.7802A>G. They generated at least 10 different transcripts: Δ16p₄₄ (alternative 3’ss 44-nt downstream; acceptor variants), Δ16 (exon 16-skipping; donor variants), Δ16p₅₅ (alternative 3’ss 55-nt downstream), Δ16q₄ (alternative 5’ss 4-nt upstream), Δ16q₁₀₀ (alternative 5’ss 4-nt upstream), ▾16q²⁰ (alternative 5’ss 20-nt downstream), as well as minor (Δ16p₉₃ and Δ16,17p₆₉) and uncharacterized transcripts of 893 and 954 nucleotides. Isoforms Δ16p₄₄, Δ16, Δ16p₅₅, Δ16q₄, Δ16q₁₀₀, and ▾16q²⁰ introduced premature termination codons which presumably inactivate BRCA2. According to the guidelines the American College of Medical Genetics and Genomics these eight variants could be classified as pathogenic or likely pathogenic whereas the Evidence-based Network for the Interpretation of Germline Mutant Alleles rules suggested seven class 4 and one class 3 variants. In conclusion, our study highlights the relevance of splicing functional assays by hybrid minigenes for the clinical classification of genetic variations. Hence, we provide new data about spliceogenic variants of BRCA2 exon 16 that are directly correlated with breast cancer susceptibility.

The BRCA2 c.68-7T > A variant is not pathogenic: A model for clinical calibration of spliceogenicity

Although the spliceogenic nature of the BRCA2 c.68-7T > A variant has been demonstrated, its association with cancer risk remains controversial. In this study, we accurately quantified by real-time PCR and digital PCR (dPCR), the BRCA2 isoforms retaining or missing exon 3. In addition, the combined odds ratio for causality of the variant was estimated using genetic and clinical data, and its associated cancer risk was estimated by case-control analysis in 83,636 individuals. Co-occurrence in trans with pathogenic BRCA2 variants was assessed in 5,382 families. Exon 3 exclusion rate was 4.5-fold higher in variant carriers (13%) than controls (3%), indicating an exclusion rate for the c.68-7T > A allele of approximately 20%. The posterior probability of pathogenicity was 7.44 × 10-115 . There was neither evidence for increased risk of breast cancer (OR 1.03; 95% CI 0.86-1.24) nor for a deleterious effect of the variant when co-occurring with pathogenic variants. Our data provide for the first time robust evidence of the nonpathogenicity of the BRCA2 c.68-7T > A. Genetic and quantitative transcript analyses together inform the threshold for the ratio between functional and altered BRCA2 isoforms compatible with normal cell function. These findings might be exploited to assess the relevance for cancer risk of other BRCA2 spliceogenic variants.

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*10²⁵. 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.

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.

Splicing Analysis of Exonic OCRL Mutations Causing Lowe Syndrome or Dent-2 Disease

Mutations in the OCRL gene are associated with both Lowe syndrome and Dent-2 disease. Patients with Lowe syndrome present congenital cataracts, mental disabilities and a renal proximal tubulopathy, whereas patients with Dent-2 disease exhibit similar proximal tubule dysfunction but only mild, or no additional clinical defects. It is not yet understood why some OCRL mutations cause the phenotype of Lowe syndrome, while others develop the milder phenotype of Dent-2 disease. Our goal was to gain new insights into the consequences of OCRL exonic mutations on pre-mRNA splicing. Using predictive bioinformatics tools, we selected thirteen missense mutations and one synonymous mutation based on their potential effects on splicing regulatory elements or splice sites. These mutations were analyzed in a minigene splicing assay. Results of the RNA analysis showed that three presumed missense mutations caused alterations in pre-mRNA splicing. Mutation c.741G>T; p.(Trp247Cys) generated splicing silencer sequences and disrupted splicing enhancer motifs that resulted in skipping of exon 9, while mutations c.2581G>A; p.(Ala861Thr) and c.2581G>C; p.(Ala861Pro) abolished a 5′ splice site leading to skipping of exon 23. Mutation c.741G>T represents the first OCRL exonic variant outside the conserved splice site dinucleotides that results in alteration of pre-mRNA splicing. Our results highlight the importance of evaluating the effects of OCRL exonic mutations at the mRNA level.

Evaluation of reported pathogenic variants and their frequencies in a Japanese population based on a whole-genome reference panel of 2049 individuals

Clarifying allele frequencies of disease-related genetic variants in a population is important in genomic medicine; however, such data is not yet available for the Japanese population. To estimate frequencies of actionable pathogenic variants in the Japanese population, we examined the reported pathological variants in genes recommended by the American College of Medical Genetics and Genomics (ACMG) in our reference panel of genomic variations, 2KJPN, which was created by whole-genome sequencing of 2049 individuals of the resident cohort of the Tohoku Medical Megabank Project. We searched for pathogenic variants in 2KJPN for 57 autosomal ACMG-recommended genes responsible for 26 diseases and then examined their frequencies. By referring to public databases of pathogenic variations, we identified 143 reported pathogenic variants in 2KJPN for the 57 ACMG recommended genes based on a classification system. At the individual level, 21% of the individuals were found to have at least one reported pathogenic allele. We then conducted a literature survey to review the variants and to check for evidence of pathogenicity. Our results suggest that a substantial number of people have reported pathogenic alleles for the ACMG genes, and reviewing variants is indispensable for constructing the information infrastructure of genomic medicine for the Japanese population.

Characterization of BRCA1 and BRCA2 variants found in a Norwegian breast or ovarian cancer cohort

Germline mutations in BRCA1 and BRCA2 cause hereditary breast and ovarian cancer. Molecular screening of these two genes in patients with a family history of breast or ovarian cancer has revealed pathogenic variants as well as genetic variants of unknown significance (VUS). These VUS may cause a challenge in the genetic counseling process regarding clinical management of the patient and the family. We investigated 32 variants previously detected in 33 samples from patients with a family history of breast or ovarian cancer. cDNA was analyzed for alternative transcripts and selected missense variants located in the BRCT domains of BRCA1 were assessed for their trans-activation ability. Although an extensive cDNA analysis was done, only three of the 32 variants appeared to affect the splice-process (BRCA1 c.213-5T>A, BRCA1 c.5434C>G and BRCA2 c.68-7T>A). In addition, two variants located in the BRCT domains of BRCA1 (c.5075A>C p.Asp1692Ala and c.5513T>G p.Val1838Gly) were shown to abolish the BRCT domain trans-activation ability, whereas BRCA1 c.5125G>A p.Gly1709Arg exhibited equal trans-activation capability as the WT domain. These functional studies may offer further insights into the pathogenicity of certain identified variants; however, this assay is only applicable for a subset of missense variants.

Alternative Splicing in Genetic Diseases: Improved Diagnosis and Novel Treatment Options

Alternative splicing is an important mechanism to regulate gene expression and to expand the repertoire of gene products in order to accommodate an increase in complexity of multicellular organisms. It needs to be precisely regulated, which is achieved via RNA structure, splicing factors, transcriptional regulation, and chromatin. Changes in any of these factors can lead to disease. These may include the core spliceosome, splicing enhancer/repressor sequences and their interacting proteins, the speed of transcription by RNA polymerase II, and histone modifications. While the basic principle of splicing is well understood, it is still very difficult to predict splicing outcome, due to the multiple levels of regulation. Current molecular diagnostics mainly uses Sanger sequencing of exons, or next-generation sequencing of gene panels or the whole exome. Functional analysis of potential splicing variants is scarce, and intronic variants are often not considered. This likely results in underestimation of the percentage of splicing variants. Understanding how sequence variants may affect splicing is not only crucial for confirmation of diagnosis and for genetic counseling, but also for the development of novel treatment options. These include small molecules, transsplicing, antisense oligonucleotides, and gene therapy. Here we review the current state of molecular mechanisms of splicing regulation and how deregulation can lead to human disease, diagnostics to detect splicing variants, and novel treatment options based on splicing correction.