An Integral Approach to the Molecular Diagnosis of Tuberous Sclerosis Complex: The Role of Mosaicism and Splicing Variants.

Tuberous sclerosis complex (TSC) is an autosomal dominant genetic disorder characterized by the presence of hamartomas in multiple organs. At the molecular level, the disease is caused by pathogenic variants in the TSC1 and TSC2 genes, and only 10% to 25% of clinically diagnosed patients remain negative after multiplex ligation-dependent probe amplification and exon sequencing of both genes. Here, to improve the molecular diagnosis of TSC, we developed an integral approach that includes multiplex ligation-dependent probe amplification and deep-coverage next-generation sequencing of the entire TSC1 and TSC2 genes, along with an adapted bioinformatic pipeline to detect variants at low allele frequencies (>1%). Using this workflow, the molecular cause was identified in 29 of 42 patients with TSC, describing here, for the first time, 12 novel pathogenic variants in TSC genes. These variants included seven splicing variants, five of which were studied at the cDNA level, determining their effect on splicing. In addition, 8 of the 29 pathogenic variants were detected in mosaicism, including four patients with previous negative study results who presented extremely low mosaic variants (allele frequency, <16%). We demonstrate that this integral approach allows the molecular diagnosis of patients with TSC and improves the conventional one by adapting the technology to the detection of low-frequency mosaics.

A Novel Intragenic Duplication in the HDAC8 Gene Underlying a Case of Cornelia de Lange Syndrome.

Cornelia de Lange syndrome (CdLS) is a multisystemic genetic disorder characterized by distinctive facial features, growth retardation, and intellectual disability, as well as various systemic conditions. It is caused by genetic variants in genes related to the cohesin complex. Single-nucleotide variations are the best-known genetic cause of CdLS; however, copy number variants (CNVs) clearly underlie a substantial proportion of cases of the syndrome. The NIPBL gene was thought to be the locus within which clinically relevant CNVs contributed to CdLS. However, in the last few years, pathogenic CNVs have been identified in other genes such as HDAC8, RAD21, and SMC1A. Here, we studied an affected girl presenting with a classic CdLS phenotype heterozygous for a de novo ~32 kbp intragenic duplication affecting exon 10 of HDAC8. Molecular analyses revealed an alteration in the physiological splicing that included a 96 bp insertion between exons 9 and 10 of the main transcript of HDAC8. The aberrant transcript was predicted to generate a truncated protein whose accessibility to the active center was restricted, showing reduced ease of substrate entry into the mutated enzyme. Lastly, we conclude that the duplication is responsible for the patient's phenotype, highlighting the contribution of CNVs as a molecular cause underlying CdLS.

Role of Splicing Regulatory Elements and In Silico Tools Usage in the Identification of Deep Intronic Splicing Variants in Hereditary Breast/Ovarian Cancer Genes.

The contribution of deep intronic splice-altering variants to hereditary breast and ovarian cancer (HBOC) is unknown. Current computational in silico tools to predict spliceogenic variants leading to pseudoexons have limited efficiency. We assessed the performance of the SpliceAI tool combined with ESRseq scores to identify spliceogenic deep intronic variants by affecting cryptic sites or splicing regulatory elements (SREs) using literature and experimental datasets. Our results with 233 published deep intronic variants showed that SpliceAI, with a 0.05 threshold, predicts spliceogenic deep intronic variants affecting cryptic splice sites, but is less effective in detecting those affecting SREs. Next, we characterized the SRE profiles using ESRseq, showing that pseudoexons are significantly enriched in SRE-enhancers compared to adjacent intronic regions. Although the combination of SpliceAI with ESRseq scores (considering ∆ESRseq and SRE landscape) showed higher sensitivity, the global performance did not improve because of the higher number of false positives. The combination of both tools was tested in a tumor RNA dataset with 207 intronic variants disrupting splicing, showing a sensitivity of 86%. Following the pipeline, five spliceogenic deep intronic variants were experimentally identified from 33 variants in HBOC genes. Overall, our results provide a framework to detect deep intronic variants disrupting splicing.

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

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

Incorporation of semi-quantitative analysis of splicing alterations for the clinical interpretation of variants in BRCA1 and BRCA2 genes.

BRCA1 and BRCA2 (BRCA1/2) genetic variants that disrupt messenger RNA splicing are commonly associated with increased risks of developing breast/ovarian cancer. The majority of splicing studies published to date rely on qualitative methodologies (i.e., Sanger sequencing), but it is necessary to incorporate semi-quantitative or quantitative approaches to accurately interpret the clinical significance of spliceogenic variants. Here, we characterize the splicing impact of 31 BRCA1/2 variants using semi-quantitative capillary electrophoresis of fluorescent amplicons (CE), Sanger sequencing and allele-specific assays. A total of 14 variants were found to disrupt splicing. Allelic-specific assays could be performed for BRCA1 c.302-1G>A and BRCA2 c.516+2T>A, c.1909+1G>A, c.8332-13T>G, c.8332-2A>G, c.8954-2A>T variants, showing a monoallelic contribution to full-length transcript expression that was concordant with semi-quantitative data. The splicing fraction of alternative and aberrant transcripts was also measured by CE, facilitating variant interpretation. Following Evidence-based Network for the Interpretation of Germline Mutant Alleles criteria, we successfully classified eight variants as pathogenic (Class 5), five variants as likely pathogenic (Class 4), and 14 variants as benign (Class 1). We also provide splicing data for four variants classified as uncertain (Class 3), which produced a "leaky" splicing effect or introduced a missense change in the protein sequence, that will require further assessment to determine their clinical significance.

Screening of BRCA1/2 deep intronic regions by targeted gene sequencing identifies the first germline BRCA1 variant causing pseudoexon activation in a patient with breast/ovarian cancer.

BACKGROUND: Genetic analysis of BRCA1 and BRCA2 for the diagnosis of hereditary breast and ovarian cancer (HBOC) is commonly restricted to coding regions and exon-intron boundaries. Although germline pathogenic variants in these regions explain about ~20% of HBOC cases, there is still an important fraction that remains undiagnosed. We have screened BRCA1/2 deep intronic regions to identify potential spliceogenic variants that could explain part of the missing HBOC susceptibility. METHODS: We analysed BRCA1/2 deep intronic regions by targeted gene sequencing in 192 high-risk HBOC families testing negative for BRCA1/2 during conventional analysis. Rare variants (MAF <0.005) predicted to create/activate splice sites were selected for further characterisation in patient RNA. The splicing outcome was analysed by RT-PCR and Sanger sequencing, and allelic imbalance was also determined when heterozygous exonic loci were present. RESULTS: A novel transcript was detected in BRCA1 c.4185+4105C>T variant carrier. This variant promotes the inclusion of a pseudoexon in mature mRNA, generating an aberrant transcript predicted to encode for a non-functional protein. Quantitative and allele-specific assays determined haploinsufficiency in the variant carrier, supporting a pathogenic effect for this variant. Genotyping of 1030 HBOC cases and 327 controls did not identify additional carriers in Spanish population. CONCLUSION: Screening of BRCA1/2 intronic regions has identified the first BRCA1 deep intronic variant associated with HBOC by pseudoexon activation. Although the frequency of deleterious variants in these regions appears to be low, our study highlights the importance of studying non-coding regions and performing comprehensive RNA assays to complement genetic diagnosis.

Alternative transcript imbalance underlying breast cancer susceptibility in a family carrying PALB2 c.3201+5G>T.

PURPOSE: Disruption of splicing motifs by genetic variants can affect the correct generation of mature mRNA molecules leading to aberrant transcripts. In some cases, variants may alter the physiological transcription profile composed of several transcripts, and an accurate in vitro evaluation is crucial to establish their pathogenicity. In this study, we have characterized a novel PALB2 variant c.3201+5G>T identified in a breast cancer family. METHODS: Peripheral blood RNA was analyzed in two carriers and ten controls by RT-PCR and Sanger sequencing. The splicing profile was also characterized by semi-quantitative capillary electrophoresis and quantitative PCR. RAD51 foci formation and PALB2 LOH status were evaluated in primary breast tumor samples from the carriers. RESULTS: PALB2 c.3201+5G>T disrupts intron 11 donor splice site and modifies the abundance of several alternative transcripts (∆11, ∆12, and ∆11,12), also present in control samples. All transcripts are predicted to encode for non-functional proteins. Semi-quantitative and quantitative analysis of PALB2 full-length transcript indicated haploinsufficiency in carriers. One tumor exhibited PALB2 LOH and RAD51 assay indicated homologous recombination deficiency in both tumors. CONCLUSIONS: Our results support a pathogenic classification for PALB2 c.3201+5G>T, highlighting the impact of variants causing an imbalanced expression of natural RNA isoforms in cancer susceptibility.

Multigene panel testing beyond BRCA1/2 in breast/ovarian cancer Spanish families and clinical actionability of findings.

PURPOSE: Few and small studies have been reported about multigene testing usage by massively parallel sequencing in European cancer families. There is an open debate about what genes should be tested, and the actionability of some included genes is under research. METHODS: We investigated a panel of 34 known high/moderate-risk cancer genes, including 16 related to breast or ovarian cancer (BC/OC) genes, and 63 candidate genes to BC/OC in 192 clinically suspicious of hereditary breast/ovarian cancer (HBOC) Spanish families without pathogenic variants in BRCA1 or BRCA2 (BRCA1/2). RESULTS: We identified 16 patients who carried a high- or moderate-risk pathogenic variant in eight genes: 4 PALB2, 3 ATM, 2 RAD51D, 2 TP53, 2 APC, 1 BRIP1, 1 PTEN and 1 PMS2. These findings led to increased surveillance or prevention options in 12 patients and predictive testing in their family members. We detected 383 unique variants of uncertain significance in known cancer genes, of which 35 were prioritized in silico. Eighteen loss-of-function variants were detected in candidate BC/OC genes in 17 patients (1 BARD1, 1 ERCC3, 1 ERCC5, 2 FANCE, 1 FANCI, 2 FANCL, 1 FANCM, 1 MCPH1, 1 PPM1D, 2 RBBP8, 3 RECQL4 and 1 with SLX4 and XRCC2), three of which also carry pathogenic variants in known cancer genes. CONCLUSIONS: Eight percent of the BRCA1/2 negative patients carry pathogenic variants in other actionable genes. The multigene panel usage improves the diagnostic yield in HBOC testing and it is an effective tool to identify potentially new candidate genes.

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.

About 1% of the breast and ovarian Spanish families testing negative for BRCA1 and BRCA2 are carriers of RAD51D pathogenic variants.

RAD51D mutations have been recently identified in breast (BC) and ovarian cancer (OC) families. Although an etiological role in OC appears to be present, the association of RAD51D mutations and BC risk is more unclear. We aimed to determine the prevalence of germline RAD51D mutations in Spanish BC/OC families negative for BRCA1/BRCA2 mutations. We analyzed 842 index patients: 491 from BC/OC families, 171 BC families, 51 OC families and 129 patients without family history but with early-onset BC or OC or metachronous BC and OC. Mutation detection was performed with high-resolution melting, denaturing high-performance liquid chromatography or Sanger sequencing. Three mutations were found in four families with BC and OC cases (0.82%). Two were novel: c.1A>T (p.Met1?) and c.667+2_667+23del, leading to the exon 7 skipping and one previously described: c.674C>T (p.Arg232*). All were present in BC/OC families with only one OC. The c.667+2_667+23del cosegregated in the family with one early-onset BC and two bilateral BC cases. We also identified the c.629C>T (p.Ala210Val) variant, which was predicted in silico to be potentially pathogenic. About 1% of the BC and OC Spanish families negative for BRCA1/BRCA2 are carriers of RAD51D mutations. The presence of several BC mutation carriers, albeit in the context of familial OC, suggests an increased risk for BC, which should be taken into account in the follow-up and early detection measures. RAD51D testing should be considered in clinical setting for families with BC and OC, irrespective of the number of OC cases in the family.

Mutation analysis of the BCCIP gene for breast cancer susceptibility in breast/ovarian cancer families.

OBJECTIVE: About 5%-10% of breast cancer is due to inherited disease predisposition. Currently, mutations in the BRCA1 and BRCA2 genes explain less than 25% of the familial clustering of breast cancer, and additional susceptibility genes are suspected. The BCCIP gene plays an important role in the regulation of gene transcription and cell proliferation and could be involved in the maintenance of genomic integrity. The BCCIP protein binds in mammalian cells to the longest conserved region of the BRCA2 protein and is required for BRCA2 stability and function, making a critical contribution to the function of BRCA2 in mediating homologous recombination. Variants in the BCCIP gene could affect the BRCA2 functionality and be associated to the familial breast/ovarian carcinogenesis. Therefore, BCCIP gene is a potential candidate for being involved in heritable cancer susceptibility. METHODS: We have screened the entire coding region and splice junctions of BCCIP in affected index cases from 215 Spanish breast/ovarian cancer families for germ line defects, using direct sequencing. RESULTS: Mutation analysis revealed 3 different intronic sequence changes. CONCLUSIONS: Based on the in silico and in vitro RNA analyses of these sequence alterations, none of them were predicted to be pathogenic or associated with cancer susceptibility. Our results indicate that BCCIP germ line mutations are unlikely to be a major contributor to familial breast/ovarian cancer risk in our population.

Mutation analysis of the SHFM1 gene in breast/ovarian cancer families.

PURPOSE: About 5-10 % of breast cancer is due to inherited disease predisposition. Currently known susceptibility genes such as BRCA1 and BRCA2 explain less than 25 % of familial aggregation of breast cancer, which suggests the involvement of additional genetic susceptibility. The SHFM1 [split hand/foot malformation (ectrodactyly) type 1] gene plays an important role in the regulation of gene transcription and cell proliferation and may be involved in the maintenance of genomic integrity. It is a potential candidate for being involved in heritable cancer susceptibility due to its biological function. The SHFM1 protein binds in mammalian cells to the longest conserved region of the BRCA2 protein and is required for BRCA2 stability and function, making a critical contribution to the BRCA2 function in mediating homologous recombination. Therefore, variants in the SHFM1 gene could affect the BRCA2 functionality and be associated with the familial breast/ovarian carcinogenesis. METHODS: We have screened the entire coding region and splice junctions of SHFM1 in affected index cases from 369 Spanish breast/ovarian cancer families for germ line defects, using direct sequencing. RESULTS: Mutation analysis revealed seven different sequence changes. Based on the in silico analyses of these sequence alterations, as well as their occurrence in cases and controls, none of them, however, were predicted to be pathogenic or associated with cancer susceptibility. CONCLUSIONS: To our knowledge, this is the most comprehensive study reporting the mutation screening of the SHFM1 gene in familial breast/ovarian cancer cases. No evidence for the association with breast/ovarian cancer was observed.

Novel BRCA1 deleterious mutation (c.1949_1950delTA) in a woman of Senegalese descent with triple-negative early-onset breast cancer.

Limited information exists regarding BRCA1 and BRCA2 genetic testing and genetic diversity in BRCA1 and BRCA2 in sub-Saharan African populations. We report a novel mutation that consists of a deletion of 2 bp (c.1949_1950delTA) in the exon 11 of the BRCA1 gene. This is a frameshift mutation that causes the disruption of the translational reading frame resulting in a premature stop codon downstream in the BRCA1 protein. The mutation was present in a Senegalese woman with a triple-negative breast tumor and a family history of breast cancer.

Identification of a new complex deleterious mutation in exon 18 of the BRCA2 gene in a hereditary male/female breast cancer family.

We report a novel complex mutation that consists of a deletion of 12 bp and an insertion of 2 bp (c.8402_8413del12ins2bp) in the exon 18 of the BRCA2 gene. This is a frameshift mutation that causes a disruption of the translational reading frame resulting in a stop codon downstream in the 2729 position of the BRCA2 protein. The mutation was present in a Spanish hereditary male/female breast cancer family.