Next Article in Journal
The PD-1/PD-L1 Gateway: Peripheral Immune Regulation in the Pathogenesis of Endometriosis
Next Article in Special Issue
Systematic Modeling of Risk-Associated Copy Number Alterations in Cancer
Previous Article in Journal
Serum Caspase-3 Levels as a Predictive Molecular Biomarker for Acute Ischemic Stroke
Previous Article in Special Issue
Genomic and Transcriptomic Research in the Discovery and Application of Colorectal Cancer Circulating Markers
 
 
ijms-logo
Article Menu

Article Menu

Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Molecular Characterization of the Allelic Expression of the BRCA1 Founder Δ9–12 Pathogenic Variant and Its Potential Clinical Relevance in Hereditary Cancer

by
Julieta Dominguez-Ortiz
1,2,3,†,
Rosa M. Álvarez-Gómez
3,†,
Rogelio Montiel-Manríquez
1,
Alberto Cedro-Tanda
4,
Nicolás Alcaraz
5,
Clementina Castro-Hernández
1,
Luis Bautista-Hinojosa
1,
Laura Contreras-Espinosa
1,
Leda Torres-Maldonado
6,
Verónica Fragoso-Ontiveros
3,
Yuliana Sánchez-Contreras
3,
Rodrigo González-Barrios
1,
Marcela Angélica De la Fuente-Hernández
3,
María de la Luz Mejía-Aguayo
3,
Ulises Juárez-Figueroa
6,
Alejandra Padua-Bracho
3,
Rodrigo Sosa-León
3,
Gabriela Obregon-Serrano
3,
Silvia Vidal-Millán
3,
Paulina María Núñez-Martínez
3,
Abraham Pedroza-Torres
3,
Sergio Nicasio-Arzeta
7,
Alfredo Rodríguez
6,8,
Fernando Luna
1,
Fernanda Cisneros-Soberanis
9,
Sara Frías
6,8,
Cristian Arriaga-Canon
1,10,* and
Luis A. Herrera-Montalvo
1,10,*
add Show full author list remove Hide full author list
1
Unidad de Investigación Biomédica en Cáncer, Instituto Nacional de Cancerología-Instituto de Investigaciones Biomédicas, Universidad Nacional Autónoma de México (UNAM), Avenida San Fernando No. 22 Col. Sección XVI, Tlalpan, Mexico City 14080, Mexico
2
Instituto Nacional de Cancerología, Universidad Nacional Autónoma de México (UNAM), Coyoacán, Mexico City 04510, Mexico
3
Clínica de Cáncer Hereditario, Instituto Nacional de Cancerología, Av. San Fernando No. 22 Col. Sección XVI, Tlalpan, Mexico City 14080, Mexico
4
Núcleo B de Innovación en Medicina de Precisión, Instituto Nacional de Medicina Genómica, Periférico Sur 4809, Arenal Tepepan, Tlalpan, Mexico City 14610, Mexico
5
Novo Nordisk Foundation Center for Protein Research, Faculty of Health and Medical Sciences, University of Copenhagen, Blegdamsvej 3A, 2200 Copenhagen, Denmark
6
Instituto Nacional de Pediatría, Insurgentes Sur No. 3700-C. Coyoacán, Mexico City 04530, Mexico
7
Natural Resource Ecology Laboratory, Colorado State University, Fort Collins, CO 80521, USA
8
Departamento de Medicina Genómica y Toxicología Ambiental, Instituto de Investigaciones Biomédicas, Universidad Nacional Autónoma de México (UNAM), México City 04510, Mexico
9
Wellcome Trust Centre for Cell Biology, ICB, University of Edinburgh, Michael Swann Building, King’s Buildings, Max Born Crescent, Edinburgh EH9 3BF, UK
10
Tecnológico de Monterrey, Escuela de Medicina y Ciencias de la Salud, Monterrey 64710, Mexico
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Int. J. Mol. Sci. 2024, 25(12), 6773; https://doi.org/10.3390/ijms25126773
Submission received: 25 April 2024 / Revised: 13 June 2024 / Accepted: 15 June 2024 / Published: 20 June 2024
(This article belongs to the Special Issue Advances in Cancer-Related Transcriptome and Genome Analyses)

Abstract

:
Hereditary breast and ovarian cancer (HBOC) syndrome is a genetic condition that increases the risk of breast cancer by 80% and that of ovarian cancer by 40%. The most common pathogenic variants (PVs) causing HBOC occur in the BRCA1 gene, with more than 3850 reported mutations in the gene sequence. The prevalence of specific PVs in BRCA1 has increased across populations due to the effect of founder mutations. Therefore, when a founder mutation is identified, it becomes key to improving cancer risk characterization and effective screening protocols. The only founder mutation described in the Mexican population is the deletion of exons 9 to 12 of BRCA1 (BRCA1Δ9–12), and its description focuses on the gene sequence, but no transcription profiles have been generated for individuals who carry this gene. In this study, we describe the transcription profiles of cancer patients and healthy individuals who were heterozygous for PV BRCA1Δ9–12 by analyzing the differential expression of both alleles compared with the homozygous BRCA1 control group using RT–qPCR, and we describe the isoforms produced by the BRCA1 wild-type and BRCA1Δ9–12 alleles using nanopore long-sequencing. Using the Kruskal–Wallis test, our results showed a similar transcript expression of the wild-type allele between the healthy heterozygous group and the homozygous BRCA1 control group. An association between the recurrence and increased expression of both alleles in HBOC patients was also observed. An analysis of the sequences indicated four wild-type isoforms with diagnostic potential for discerning individuals who carry the PV BRCA1Δ9–12 and identifying which of them has developed cancer.

1. Introduction

The hereditary component of cancer contributes to 10% of all cancer cases globally. Among these, hereditary breast and ovarian cancer syndrome (HBOC) is the most prevalent, constituting 10% of all breast cancer (BrCa) cases and 20% of all ovarian cancer (OvCa) cases [1,2,3,4]. HBOC is linked to heterozygous and germline pathogenic variants (PVs) in cancer predisposition genes, with BRCA1 and BRCA2 being the most commonly affected genes [5,6]. Identifying individuals with PVs in HBOC has the following significant implications: (I) facilitating timely cancer diagnosis through high-risk screening methods (e.g., magnetic resonance imaging for breast cancer); (II) enabling targeted therapies, such as poly(ADP-ribose) polymerase (PARP) inhibitors, recognized as treatment options for major HBOC-related malignancies; (III) implementing cancer risk reduction strategies, including risk-reducing surgeries (mastectomy and/or salpingo-oophorectomy) or chemoprevention with agents such as tamoxifen or raloxifene for breast cancer; and (IV) identifying asymptomatic individuals at a high risk for cancer through cascade molecular diagnostics, thereby enabling the effective prevention of cancer morbidity and mortality [1,2,3,4,5,6,7]. This comprehensive approach underscores the importance of early detection, targeted interventions, and personalized risk management strategies in mitigating the impact of HBOC.
Hence, systematically studying PVs within HBOC genes is paramount, with a focus on understanding their implications for specific phenotypes. This entails exploring phenotype–genotype correlations in terms of tumor risks, age of onset, differential responses to pharmacological treatments, and drug resistance. Additionally, it is crucial to investigate how these variants differ among populations based on their ancestry. For instance, certain PVs may be prevalent within populations due to founder effects, which result in reduced genetic diversity and an increased frequency of specific genetic variants [8,9]. The identification of founder mutations facilitates the characterization of HBOC patients and enhances screening strategies to identify individuals and families harboring PVs [10,11].
Within the Mexican population, a single BRCA1 founder PV has been documented: the exon 9–12 deletion (BRCA1Δ9–12) [1,11,12]. BRCA1Δ9–12 arises from splicing between the AluSp element in intron 8 and the AluSx element in intron 12, resulting in a significant structural rearrangement and a loss of 15.4 kb of the BRCA1 gene [12]. This deletion induces a premature stop codon in the mRNA, suggesting its likely deleterious nature [12]. Moreover, the exon 9–10 deletion is the most prevalent PV identified in Mexico, associated with a heightened risk of high-grade epithelial ovarian cancer occurring at a younger age. The exon 9–10 deletion also exhibits an extended favorable response to the PARP inhibitor olaparib for recurrent ovarian cancer [1,11,12,13].
However, there is a dearth of studies elucidating the mechanisms underlying these clinical discrepancies compared to observations with other pathogenic variants (PVs). Specifically, there is a notable absence of research on the isoforms generated from both the BRCA1 wild-type (WT) and BRCA1Δ9–12 (Δ9–12) alleles. The evidence suggests that transcripts of BRCA1 PVs may evade nonsense-mediated mRNA decay (NMD) mechanisms via isoform splicing events or the upregulation of WT isoforms, thereby rescuing the homologous recombination (HR) response pathway involving BRCA1. This phenomenon has been associated with resistance to PARP inhibitors and platinum-based therapies [14,15,16,17,18,19]. However, it remains unknown whether the BRCA1Δ9–12 mutation generates a translatable transcript or modulates the expression of WT isoforms.
This study aimed to characterize the transcript expression of both WT and Δ9–12 alleles in cancer patients (CaH) and healthy individuals (HH) heterozygous for the BRCA1Δ9–12 PV, as well as their relation to the homozygous WT status of healthy individuals not carrying the PV (Ctrl). Our preliminary findings suggest that specific expression profiles may hold significance in the early detection of this founder PV in asymptomatic individuals, offering a cost-effective alternative via qPCR for diagnosing and monitoring asymptomatic individuals (Scheme 1).

2. Results

2.1. Clinical Characteristics of the Patient Cohort

The characterization of the patient cohort was an important goal of this study and the starting point for establishing a correlation between transcription allele differences across the groups. The main clinical features were described in relation to tumor pathological characteristics, treatment, disease-free survival time, olaparib maintenance treatment, and recurrence (Table 1). The cohort comprised 11 female patients, all of whom were in the CaH group. One patient (CA11) lacked a description of her clinical characteristics, since her treatment was conducted at another institution and her clinical history was unavailable; her cancer diagnosis and BRCA1 PV confirmation were corroborated. Among the following 10 patients, the age of primary tumor presentation ranged from 30 to 49 years. Four out of ten patients with available clinical histories were diagnosed with ovarian cancer as their first tumor, while the remaining six were diagnosed with breast cancer. The subsequent development of a secondary primary tumor was reported in three of the breast cancer patients and two of the ovarian cancer patients. The disease-free survival period ranged from 27 to 115 months, and six patients were reported to have cancer recurrence, none of which were breast–ovary related. Treatment for four of the six breast cancer patients included radiotherapy, and radioresistance was not reported in any of the patients (Table 1).

2.2. Analysis of BRCA1WT and BRCA1Δ9–12 Transcript Expression in the Δ9–12 Heterozygous Groups and Controls

To describe and establish whether transcriptional differences were present between the allele expression among the patients who presented with HBOC and those who did not develop cancer and to compare it to that of the healthy Ctrl, we designed a specific set of primers to differentiate between the WT and Δ9–12 transcripts according to the size of the amplicons: 210 bp for the WT and 127 bp for the Δ9–12 transcripts (Figure 1A,B).
To ensure primer specificity, a primer mix reaction was utilized where the characteristic band pattern expected for each group could be visualized, one band in the Ctrl group showing the BRCA1 homozygous WT status, and two bands in both the CaH and HH groups exhibiting the BRCA1 heterozygous WT and Δ9–12 alleles (Figure 1C,D). Sanger sequencing confirmed that exons 8–9 were spliced in the 210 bp WT amplicon and that exons 8–13 were spliced in the 127 bp Δ9–12 amplicon (Figure 1E,F). Transcript amplification using a single primer set was also performed to verify a single amplicon for each sample (Supplementary Figure S2), along with the visualization of representative samples of each group in one electrophoresis gel. The samples were amplified by the primer mix reaction and the specific primers for the WT and Δ9–12 transcripts.
Once the specificity of our two sets of primers was verified, we used an RT–qPCR approach to evaluate the transcript expression. We found differences in the WT transcript expression between the groups (χ2 = 12.37; p = 0.002; d.f. = 2), except between the WT transcripts in the Ctrl (0.0094 ± 0.0028) and HH (0.0075 ± 0.0024) groups. The Δ9–12 allele also showed significant differences in gene expression (χ2 = 13.91; p = 0.003; d.f. = 3). The Δ9–12 allele showed a lower expression in HH (0.0043 ± 0.0017) and CaH (0.0026 ± 0.0013) than in their wild-type counterparts (HH: 0.0075 ± 0.0024; CaH: 0.0044 ± 0.0025). A significant difference was also observed between the Δ9–12 allele in CaH and the WT allele in HH (W = 2; p = 0.0003) (Figure 1G). The above patterns were not observed in any of the heterozygous groups, where the Δ9–12 allele expression remained below that of the WT allele; even though a slight increase was observed in the HH group, the difference was not significant. Equal WT transcript expression between HH and Ctrl individuals, as evaluated in lymphocytes, has not been reported, which makes this result an interesting antecedent for healthy PV heterozygous individuals.

2.3. Correlation of the WT and Δ9–12 CaH Transcript Expression Alleles with Clinical Characteristics

After identifying the differences in transcript expression for each allele in the CaH group, different clinical characteristics were evaluated to identify correlations suggesting new information about this PV and the patient’s clinical presentation. We compared previously related clinical characteristics to changes in allele-specific expression patterns for PVs in BRCA1 [20]. Cancer recurrence was the only clinical characteristic showing significant differences for the WT and Δ9–12 alleles (Figure 2A,B). These findings provide a new observation not described for BRCA1 PVs, adding to their genotype–phenotype description and concluding our evaluation of the transcription among HBOC patients.

2.4. Isoform Expression from the BRCA1Δ9–12 Allele

Another main objective of this work was to identify and study the isoforms and sequence heterogeneity in our group’s samples. To achieve this goal, we amplified the full-length BRCA1 WT and Δ9–12 isoform sequences by endpoint PCR using a single primer set that was complementary to the adjacent region of the starting and stop codons (Figure 3A,B). Once the amplicons were obtained from each sample, our first approach to confirm the amplification of the BRCA1 allele isoforms was to select representative samples from each group to be observed in a 1% gel (Figure 3C) and corroborated by Sanger sequencing of the purified bands, which corresponded to the expected weight of these transcripts and the joint of exons 8 and 9 for the WT strain and exons 8 and 13 for the Δ9–12 strain (Figure 3D).
We sequenced all the amplified isoforms from the samples using nanopore technology. The identified isoforms were analyzed using principal component analysis (PCA) of the top 79 expressed isoforms; the samples were separated by differentiation state, mostly along PC1. The expression patterns correlated and grouped the samples according to heterozygous Δ9–12 (CaH and HH) and homozygous WT (Ctrl) characteristics, as the PCA plot exhibited greater variance (35.49%) between these groups (Figure 4A). The CaH and HH groups clustered differently from each other on the PC3 axis, although with less variance (8.3%), as visualized in a 3D PCA plot (Figure 4B). A second analysis was undertaken using a heatmap of dissimilarity. Hierarchical clustering showed two different areas, where the Ctrl samples clustered within themselves with high similarity, and in the other panel, the CaH and HH samples correlated in another area (Figure 4C). This plot suggests that the transcriptional changes in isoform expression can distinguish between the Ctrl group and individuals with the BRCA1Δ9–12 PV (CaH and HH), and confirms the differential expression of WT isoforms when the PV BRCA1Δ9–12 is present.

2.5. Potential Predictors of PV and Cancer

Four isoforms were among the 79 BRCA1 isoforms identified according to changes in expression. Two of them, isoform 1 (ENST00000352993.7) and isoform 2 (ENST00000484087.6), were overexpressed in samples with BRCA1Δ9–12 (CaH and HH) and underexpressed in the Ctrl samples, allowing the differentiation between heterozygotes for PV BRCA1Δ9–12 and healthy homozygotes for BRCA1. Isoforms 3 (ENST00000700082.1) and 4 (ENST00000618469.2) differentiated the samples within the CaH group from the HH group by being overexpressed in the HH samples and subexpressed in the CaH samples (Figure 4D). These results suggest that these isoforms could be used to differentiate the population of PV heterozygotes and identify those who develop cancer, making them potential biomarkers for identifying BRCA1 PV patients.

3. Discussion

It is estimated that 10% of newly diagnosed breast and/or ovarian cancers are hereditary [21]. HBOC is an inherited predisposition to breast and/or ovarian cancers due to specific gene mutations [22]. The most important genes associated with hereditary breast and ovarian cancers are BRCA1 and BRCA2 [23,24,25]. Therefore, testing individuals for specific gene mutations and detecting PVs provides clinical information for genetic counseling, enabling directed oncological management and risk reduction surgeries [26,27]. These advantages expanded when PV detection increased for a founder mutation, such as BRCA1Δ9–12, in the population [9,10,12].
The benefits of identifying PVs in the population are not limited to the gene but also to the isoforms transcribed from it. The presence of specific variants, such as BRCA12288delT and BRCA1185delAG, has been shown to have a substantial impact on patients who are heterozygous for these PVs, as resistance to PARP inhibitor treatment and platinum in the tumor via the rescue of the HR pathway has been explained through isoforms [17,18,28]. However, even though studying the mechanisms by which the isoforms provide molecular advantages in tumor resistance, and thus tumor survival, has been key to understanding their relevance in the context of cancer, their identification and description are the first fundamental steps in scaling this knowledge.

3.1. Group Identification through Specific Allele Transcript Expression

This study describes the changes in the expression of heterozygous allele transcripts of the pathogenic variant 9–12. We aimed to characterize the population that carries this PV and identify differences in its transcription, which may lead to significant clinical benefits from the observed correlations. The correlation between transcription and clinical associations in BRCA1 PVs has been previously explored [20,29,30]. However, this work addresses this issue from an integrative perspective, including heterozygous healthy individuals, patients, healthy individuals, their non-PV counterparts, and allele-specific transcript assessment. The HH group and the individual allele evaluation provided information that allowed us to identify equal WT transcript expression between the two healthy groups, independent of their PV status, and observe the subexpression of the WT transcript only in cancer patients. Considering these results, the prospect of a biomarker for the follow-up of families with PV should be explored.
Given the limited cohort size and the absence of data that continue to inform us of changes in expression during cancer development, we cannot claim that gene compensation is the biological mechanism of this result. However, this approach should be explored as a possible strategy for healthy heterozygote cells to counter the haploinsufficiency of their WT allele [31,32] and avoid carcinogenesis, as has been proposed for biological models such as zebrafish, Drosophila, mice, human breast cancer cell cultures, and colorectal cancer tissues [32,33,34].

3.2. Isoform Pattern and Clinical Correlation

Among the isoforms that allow the study groups to be differentiated, isoforms 1, 2, and 4 have translational potential. Regardless, only isoform 4 was described as 100% similar to the BRCA1 variant 1 protein. This suggests that isoform 4 can potentially rescue the HR repair pathway, as observed in other BRCA1 VPs [15], possibly influencing the similarity observed in the allele-specific transcript analysis. Therefore, functional analyses need to be performed to investigate the effect of this isoform at the transcript and protein levels.
In summary, this work significantly contributes to describing the genotypic characteristics of the BRCA1Δ9–12 PV at the transcriptional level and proposes two main strategies for diagnostic tools: the use of changes in transcript expression and the use of four isoforms with biomarker potential. To translate these results to clinical practice, the further evaluation of patients and their healthy relatives in a larger cohort and a prospective follow-up are needed to delineate when expression levels change between the two groups. Subsequent research into the biological mechanism by which WT allele expression changes when cancer occurs should be conducted. It is fundamental to continue describing the new isoforms identified during this project, not only for their translation potential but also to offer a breakthrough in the characterization of Hispanic populations.

4. Materials and Methods

4.1. Sample Collection

This study included a total of 40 RNA samples, 11 of which came from leukocytes of CaH of the Hereditary Cancer Clinic of the National Cancer Institute (INCan) in Mexico and nine from HH of the patients’ families; all 20 samples were heterozygous for PV Δ9–12. The BRCA1Δ9–12 mutation was confirmed in all the participants, and complete gene sequencing confirmed no other mutation in either allele; these results were published by Fragoso in 2019. Familiar extensions for diagnosis and family consultation were followed at INCan. Informed consent was obtained from all groups. The BRCA1 Ctrl yielded the remaining 20 RNA samples. All the samples were collected under the Ethics in Research Committee of the National Cancer Institute (CEI/1036/16). The sample selection workflow followed throughout the study is described in the supplementary figures (Supplementary Figure S1).

4.2. RNA Extraction

RNA was extracted from the cells using Direct-zolTM RNA MiniPrep (ZYMO, Irvine, CA, USA, ref. R2052) according to the manufacturer’s instructions. Total RNA was quantified using a NanoDrop 3300 (Thermo Scientific, Waltham, MA, USA), and the OD260/280 ratio was determined. The RNA integrity number (RIN) was determined using an Agilent Bioanalyzer 2100 with an RNA 6000 Nano Assay (Agilent Technologies, Santa Clara, CA, USA).

4.3. cDNA Synthesis and Endpoint PCR Assay

Two different protocols were used, considering the purpose of the amplified transcripts. Protocol 1 aimed to obtain the amplicons of a partial region of the BRCAWT and BRCA1Δ9–12 transcripts for expression evaluation by RT–qPCR, and Protocol 2 aimed to amplify the complete sequence of the transcripts and their isoforms for analysis by nanopore sequencing.

4.3.1. Protocol 1

cDNA synthesis was performed using a High-Capacity cDNA Reverse Transcription Kit (Thermo Scientific, Waltham, MA, USA, ref. 4368814) according to the manufacturer’s instructions.
The endpoint PCR was carried out with a total volume of 20 µL with 1X DreamTaq Buffer, 0.06 U of DreamTaq DNA Polymerase (Thermo Scientific, Waltham, MA, USA, ref. EP0701), a concentration of 0.4 µM dNTPs (Thermo Scientific, Waltham, MA, USA, ref. R0181), 0.4 µM of each primer (two sets of primers were used to amplify each transcript allele; detailed information is given in the primers section), and 1 μL of cDNA. The thermal cycling conditions were 94 °C for 5 min, followed by 40 cycles of 94 °C for 30 s, primer annealing at 59 °C for 30 s, primer extension at 72 °C for 45 s, and a final extension at 72 °C for 7 min. PCR products were resolved in a 1% agarose gel using electrophoreses for 45 min, and a 1 Kb ladder was used (Axygen, Central Ave Union City, CA, USA, ref AXY1016). This method corresponds to the results presented in Figure 1 and Figure 2. The band size comparison between alleles in the studied groups can be seen in Supplementary Figure S3.

4.3.2. Protocol 2

cDNA and PCR endpoints were synthesized from 500 ng of RNA using the SuperScript IV One-Step RT–PCR System (Thermo Scientific, Waltham, MA, USA, ref. 12594100) following the manufacturer’s instructions for protocol amplification and the BRCA1 Jong ex1/12 primer set, as detailed in the primers section. The conditions used for the thermal cycler program were as follows: reverse activation, 55 °C for 10 min; RT inactivation/initial denaturation, 98 °C for 2; amplification cycle, 36 cycles of 98 °C for 10 s, 60 °C for 30 s, and 72 °C for 3:30 min; and a final extension at 72 °C for 5 min. The final products were purified using a QIAquick PCR Purification Kit (QIAGEN, Germantown, MD, USA, ref. 28106) and resolved on a 1% agarose gel using electrophoresis for 45 min. A 1 Mb ladder was used (Rockland Immunochemicals, ref MB-204-0500). Representative amplicons from this protocol are presented in Figure 3 and Supplementary Figure S4.

4.4. Primers

All primers were custom-designed by the researchers, except those used for amplifying the complete sequence of BRCA1 (Table 2), which were obtained from the primer set published by Lucy C. de Jong et al. in 2017 [15].

4.5. RT–qPCR Assay

Ten microliters of RT–qPCR was performed using Maxima SYBR Green/ROX qPCR Master Mix (Thermo Scientific, Waltham, MA, USA, ref. K0222) and 0.16 µM of each primer. The thermal cycling conditions were 94 °C for 5 min, followed by 40 cycles of 94 °C for 30 s, primer annealing at 59 °C for 30 s, primer extension at 72 °C for 45 s, and a final extension at 72 °C for 7 min. Relative gene expression levels were assessed using ΔCT for the different sample groups, with normalization to the housekeeping gene GAPDH. This method is shown in the graphic in Figure 3.

4.6. GridION Library Preparation and Sequencing

We selectively included only PCR samples with concentrations exceeding 200 ng, comprising 10 control samples, seven carrier samples, and eight patient samples. The concentrations of these PCR products were determined using a Qubit dsDNA HS Assay Kit (Molecular Probes, Life Technologies, Eugene, OR, USA, ref. Q32851) according to the manufacturer’s guidelines. Each sample was subjected to amplification, yielding several amplicons of varying sizes that ranged from 200 to 6000 bp.
For sample preparation, we used the NEBNext Ultra II End Repair and dA-tailing Kit (New England BioLabs Inc., Ipswich, MA, USA). Native barcoding of the amplicons was conducted in a 20 µL reaction volume, which included 1.5 µL of DNA-repair amplicons, 2.5 µL of Native Barcode EXP-NBD104 and EXP-NBD114 (Oxford Nanopore Technologies, Oxford, Oxfordshire, UK), 10 µL of NEBNext Ultra II Ligation Master Mix (New England Biolabs, Ipswich, MA, USA), 0.5 µL of NEBNext Ligation Enhancer (New England Biolabs, Ipswich, MA, USA), and 5.5 µL of nuclease-free water. The mixture was incubated for 20 min at 20 °C, followed by 10 min at 65 °C. Subsequently, the barcoded amplicons were pooled.
For sequencing, we used the Ligation Sequencing Kit 1D (SQK-LSK110) according to the manufacturer’s protocol for adapter ligation. The library was then purified using AMPure XP beads and quantified using a Qubit dsDNA HS Assay Kit. We loaded 50 ng of the library onto an R.9.4.1 flow cell for sequencing on the GridION platform (Oxford Nanopore Technologies, Oxford, Oxfordshire, UK) over a period of 72 h. ONT MinKNOW software (Oxford Nanopore Technologies, Oxford, Oxfordshire, UK) was used to collect the raw sequencing data. We confirmed the availability of active pores upon receipt and immediately before the sequencing run. Detailed information is given in the primers section. This method is shown in Figure 4.

4.7. Bioinformatic Analysis

Fastq reads were trimmed for adapters and quality with Pychopper. The trimmed reads were then mapped to the genome GRCh38 built using minimap2 in splice mode using the parameters -ax splice --MD. The quantification of BRCA1 transcripts (known and novel) was performed by mapping to the annotated gencode.v44 transcripts with salmon and the long-read parameter --ont --noLengthCorrection. Differential transcript expression was performed with DESeq2. An unsupervised method was applied.

4.8. Statistical Analysis

Considering the sample size, a normal distribution was assessed using the Lilliefors, Kolmogorov–Smirnov, and Shapiro–Wilk tests. The results confirmed that the samples did not show a normal distribution (p > 0.05). Therefore, we analyzed RT–qPCR differences among the three groups (cancer patients, healthy heterozygotes, and controls) using a nonparametric Kruskal–Wallis test. Furthermore, we employed a paired Wilcoxon test with a Bonferroni adjustment as a post hoc analysis to assess pairwise differences. The differential expression of both BRCA1 transcripts in the clinical aspects of CaH was evaluated using the Wilcoxon test in conjunction with the Bonferroni adjustment. The differential statistic was established at p < 0.05, and the range of differences was indicated by the following values: * < 0.05, ** < 0.01, ***< 0.001.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/ijms25126773/s1.

Author Contributions

The authors’ contributions are as follows: Conceptualization, J.D.-O., C.A.-C. and L.A.H.-M.; methodology, J.D.-O., C.A.-C. and L.A.H.-M.; validation, J.D.-O., R.M.-M. and A.C.-T.; formal analysis, J.D.-O., N.A., S.N.-A., A.P.-B., R.S.-L. and U.J.-F.; investigation, J.D.-O., R.M.-M., A.C.-T., L.T.-M., S.F., M.A.D.l.F.-H., M.d.l.L.M.-A., V.F.-O., F.L., A.R., F.C.-S., S.F. and A.P.-T.; resources, C.C.-H., C.A.-C., R.G.-B. and L.A.H.-M.; data curation, J.D.-O., Y.S.-C., G.O.-S., S.V.-M., P.M.N.-M., R.M.Á.-G., L.C.-E. and L.B.-H.; writing—original draft preparation, J.D.-O., C.A.-C., R.M.Á.-G., A.C.-T., N.A., A.R. and L.A.H.-M.; writing—review and editing, J.D.-O., C.A.-C. and L.A.H.-M.; visualization, J.D.-O., C.A.-C. and L.B.-H.; supervision, C.A.-C., C.C.-H. and L.A.H.-M.; project administration, J.D.-O., C.A.-C. and L.A.H.-M.; funding acquisition, L.A.H.-M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by “Programa presupuestario Anexo 13 del Decreto de PEF: 309 Clínica de Cáncer Hereditario” at the National Cancer Institute of Mexico and Consejo Nacional de Humanidades, Ciencias y Tecnologías (CONAHCYT). Julieta Dominguez-Ortiz is a Ph.D. student from the Programa de Doctorado en Ciencias Biomédicas, Universidad Nacional Autónoma de México (UNAM), and she received fellowship 465336 from the Consejo Nacional de Humanidades, Ciencias y Tecnologías (CONAHCYT), Mexico, with CVU 412880.

Institutional Review Board Statement

The ethics committee INCan (Comité de Ética en Investigación, Instituto Nacional de Cancerología, Mexico: 12 CEI0901411) reviews studies involving human participants. Institutional Review Board approval was obtained (CEI/1036/16).

Informed Consent Statement

Informed consent was obtained from all the subjects involved in the study.

Data Availability Statement

Sequencing data are available in Sequence Read Archive (SRA) with BioProject: PRJNA1113661.

Acknowledgments

We are deeply grateful to the patients and their families for their interest and contribution to this work and to the National Cancer Institute of Mexico (INCan) for their support. We would like to sincerely thank the entire team at the hereditary clinic, especially Soco and Lupita, whose participation was key in the execution of the project. We thank the LXII and LXIII Legislatures of the Congress of the Union of Mexico for their support. The design study figure was created using BioRender.com.

Conflicts of Interest

The authors declare no conflicts of interest.

Correction Statement

This article has been republished with a minor correction to the Funding statement. This change does not affect the scientific content of the article.

References

  1. Villarreal-Garza, C.; Alvarez-Gómez, R.M.; Pérez-Plasencia, C.; Herrera, L.A.; Herzog, J.; Castillo, D.; Mohar, A.; Castro, C.; Gallardo, L.N.; Gallardo, D.; et al. Significant Clinical Impact of Recurrent BRCA1 and BRCA2 Mutations in Mexico. Cancer 2015, 121, 372–378. [Google Scholar] [CrossRef] [PubMed]
  2. Narod, S.A. Modifiers of Risk of Hereditary Breast Cancer. Oncogene 2006, 25, 5832–5836. [Google Scholar] [CrossRef] [PubMed]
  3. Toss, A.; Tomasello, C.; Razzaboni, E.; Contu, G.; Grandi, G.; Cagnacci, A.; Schilder, R.J.; Cortesi, L. Hereditary Ovarian Cancer: Not Only BRCA 1 and 2 Genes. BioMed Res. Int. 2015, 2015, 341723. [Google Scholar] [CrossRef] [PubMed]
  4. Daly, M.B.; Pal, T.; Berry, M.P.; Buys, S.S.; Dickson, P.; Domchek, S.M.; Elkhanany, A.; Friedman, S.; Goggins, M.; Hutton, M.L.; et al. Genetic/Familial High-Risk Assessment: Breast, Ovarian, and Pancreatic, Version 2.2021, NCCN Clinical Practice Guidelines in Oncology. J. Natl. Compr. Cancer Netw. 2021, 19, 77–102. [Google Scholar] [CrossRef] [PubMed]
  5. Torres-Mejía, G.; Royer, R.; Llacuachaqui, M.; Akbari, M.R.; Giuliano, A.R.; Martínez-Matsushita, L.; Angeles-Llerenas, A.; Ortega-Olvera, C.; Ziv, E.; Lazcano-Ponce, E.; et al. Recurrent BRCA1 and BRCA2 Mutations in Mexican Women with Breast Cancer. Cancer Epidemiol. Biomark. Prev. 2015, 24, 498–505. [Google Scholar] [CrossRef] [PubMed]
  6. Fairoosa, P.; Witharana, C. Gene Mutations in Hereditary Breast Cancer—A Review. Eur. J. Med. Health Sci. 2020, 2. [Google Scholar] [CrossRef]
  7. Zhao, L.; Lynch, L.; Eiriksson, L. Information Needs of Lynch Syndrome and BRCA 1/2 Mutation Carriers Considering Risk-Reducing Gynecological Surgery: A Qualitative Study of the Decision-Making Process. Hered. Cancer Clin. Pract. 2024, 22, 5. [Google Scholar] [CrossRef]
  8. Kivisild, T. Founder Effect. In Brenner’s Encyclopedia of Genetics, 2nd ed.; Maloy, S., Hughes, K., Eds.; Academic Press: San Diego, CA, USA, 2013; pp. 100–101. ISBN 978-0-08-096156-9. [Google Scholar]
  9. Neuhausen, S.L. Founder Populations and Their Uses for Breast Cancer Genetics. Breast Cancer Res. 2000, 2, 77. [Google Scholar] [CrossRef]
  10. Ossa, C.A.; Torres, D. Founder and Recurrent Mutations in BRCA1 and BRCA2 Genes in Latin American Countries: State of the Art and Literature Review. Oncologist 2016, 21, 832–839. [Google Scholar] [CrossRef] [PubMed]
  11. Fragoso-Ontiveros, V.; Velázquez-Aragón, J.A.; Nuñez-Martínez, P.M.; de la Luz Mejía-Aguayo, M.; Vidal-Millán, S.; Pedroza-Torres, A.; Sánchez-Contreras, Y.; Ramírez-Otero, M.A.; Muñiz-Mendoza, R.; Domínguez-Ortíz, J.; et al. Mexican BRCA1 Founder Mutation: Shortening the Gap in Genetic Assessment for Hereditary Breast and Ovarian Cancer Patients. PLoS ONE 2019, 14, e0222709. [Google Scholar] [CrossRef]
  12. Weitzel, J.N.; Lagos, V.I.; Herzog, J.S.; Judkins, T.; Hendrickson, B.; Ho, J.S.; Ricker, C.N.; Lowstuter, K.J.; Blazer, K.R.; Tomlinson, G.; et al. Evidence for Common Ancestral Origin of a Recurring BRCA1 Genomic Rearrangement Identified in High-Risk Hispanic Families. Cancer Epidemiol. Biomark. Prev. 2007, 16, 1615–1620. [Google Scholar] [CrossRef]
  13. Gallardo-Rincón, D.; Álvarez-Gómez, R.M.; Montes-Servín, E.; Toledo-Leyva, A.; Montes-Servín, E.; Michel-Tello, D.; Alamilla-García, G.; Bahena-González, A.; Hernández-Nava, E.; Fragoso-Ontiveros, V.; et al. Clinical Evaluation of BRCA1/2 Mutation in Mexican Ovarian Cancer Patients. Transl. Oncol. 2020, 13, 212–220. [Google Scholar] [CrossRef] [PubMed]
  14. Nesic, K.; Krais, J.J.; Vandenberg, C.J.; Wang, Y.; Patel, P.; Cai, K.Q.; Kwan, T.; Lieschke, E.; Ho, G.-Y.; Barker, H.E.; et al. BRCA1 Secondary Splice-Site Mutations Drive Exon-Skipping and PARP Inhibitor Resistance. medRxiv 2023. [Google Scholar] [CrossRef] [PubMed]
  15. de Jong, L.C.; Cree, S.; Lattimore, V.; Wiggins, G.A.R.; Spurdle, A.B.; Miller, A.; Kennedy, M.A.; Walker, L.C.; kConFab Investigators. Nanopore Sequencing of Full-Length BRCA1 mRNA Transcripts Reveals Co-Occurrence of Known Exon Skipping Events. Breast Cancer Res. 2017, 19, 127. [Google Scholar] [CrossRef] [PubMed]
  16. Wang, Y.; Bernhardy, A.J.; Cruz, C.; Krais, J.J.; Nacson, J.; Nicolas, E.; Peri, S.; van der Gulden, H.; van der Heijden, I.; O’Brien, S.W.; et al. The BRCA1-Δ11q Alternative Splice Isoform Bypasses Germline Mutations and Promotes Therapeutic Resistance to PARP Inhibition and Cisplatin. Cancer Res. 2016, 76, 2778–2790. [Google Scholar] [CrossRef] [PubMed]
  17. Wang, Y.; Krais, J.J.; Bernhardy, A.J.; Nicolas, E.; Cai, K.Q.; Harrell, M.I.; Kim, H.H.; George, E.; Swisher, E.M.; Simpkins, F.; et al. RING Domain-Deficient BRCA1 Promotes PARP Inhibitor and Platinum Resistance. J. Clin. Investig. 2016, 126, 3145–3157. [Google Scholar] [CrossRef] [PubMed]
  18. Johnson, N.; Johnson, S.F.; Yao, W.; Li, Y.-C.; Choi, Y.-E.; Bernhardy, A.J.; Wang, Y.; Capelletti, M.; Sarosiek, K.A.; Moreau, L.A.; et al. Stabilization of Mutant BRCA1 Protein Confers PARP Inhibitor and Platinum Resistance. Proc. Natl. Acad. Sci. USA 2013, 110, 17041–17046. [Google Scholar] [CrossRef] [PubMed]
  19. Tammaro, C.; Raponi, M.; Wilson, D.I.; Baralle, D. BRCA1 Exon 11 Alternative Splicing, Multiple Functions and the Association with Cancer. Biochem. Soc. Trans. 2012, 40, 768–772. [Google Scholar] [CrossRef]
  20. Kote-Jarai, Z.; Matthews, L.; Osorio, A.; Shanley, S.; Giddings, I.; Moreews, F.; Locke, I.; Evans, D.G.; Eccles, D.; Carrier Clinic Collaborators; et al. Accurate Prediction of BRCA1 and BRCA2 Heterozygous Genotype Using Expression Profiling after Induced DNA Damage. Clin. Cancer Res. 2006, 12, 3896–3901. [Google Scholar] [CrossRef]
  21. Loboda, A.P.; Adonin, L.S.; Zvereva, S.D.; Guschin, D.Y.; Korneenko, T.V.; Telegina, A.V.; Kondratieva, O.K.; Frolova, S.E.; Pestov, N.B.; Barlev, N.A. BRCA Mutations-The Achilles Heel of Breast, Ovarian and Other Epithelial Cancers. Int. J. Mol. Sci. 2023, 24, 4982. [Google Scholar] [CrossRef]
  22. Jasiak, A.; Koczkowska, M.; Stukan, M.; Wydra, D.; Biernat, W.; Izycka-Swieszewska, E.; Buczkowski, K.; Eccles, M.R.; Walker, L.; Wasag, B.; et al. Analysis of BRCA1 and BRCA2 Alternative Splicing in Predisposition to Ovarian Cancer. Exp. Mol. Pathol. 2023, 130, 104856. [Google Scholar] [CrossRef] [PubMed]
  23. Jelovac, D.; Armstrong, D.K. Recent Progress in the Diagnosis and Treatment of Ovarian Cancer. CA Cancer J. Clin. 2011, 61, 183–203. [Google Scholar] [CrossRef] [PubMed]
  24. Karami, F.; Mehdipour, P. A Comprehensive Focus on Global Spectrum of BRCA1 and BRCA2 Mutations in Breast Cancer. BioMed Res. Int. 2013, 2013, 928562. [Google Scholar] [CrossRef] [PubMed]
  25. Mehrgou, A.; Akouchekian, M. The Importance of BRCA1 and BRCA2 Genes Mutations in Breast Cancer Development. Med. J. Islam. Repub. Iran. 2016, 30, 369. [Google Scholar] [PubMed]
  26. Tung, N.M.; Boughey, J.C.; Pierce, L.J.; Robson, M.E.; Bedrosian, I.; Dietz, J.R.; Dragun, A.; Gelpi, J.B.; Hofstatter, E.W.; Isaacs, C.J.; et al. Management of Hereditary Breast Cancer: American Society of Clinical Oncology, American Society for Radiation Oncology, and Society of Surgical Oncology Guideline. J. Clin. Oncol. 2020, 38, 2080–2106. [Google Scholar] [CrossRef] [PubMed]
  27. Richards, S.; Aziz, N.; Bale, S.; Bick, D.; Das, S.; Gastier-Foster, J.; Grody, W.W.; Hegde, M.; Lyon, E.; Spector, E.; et al. Standards and Guidelines for the Interpretation of Sequence Variants: A Joint Consensus Recommendation of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology. Genet. Med. 2015, 17, 405–424. [Google Scholar] [CrossRef] [PubMed]
  28. Seo, A.; Steinberg-Shemer, O.; Unal, S.; Casadei, S.; Walsh, T.; Gumruk, F.; Shalev, S.; Shimamura, A.; Akarsu, N.A.; Tamary, H.; et al. Mechanism for Survival of Homozygous Nonsense Mutations in the Tumor Suppressor Gene BRCA1. Proc. Natl. Acad. Sci. USA 2018, 115, 5241–5246. [Google Scholar] [CrossRef]
  29. Vuillaume, M.-L.; Uhrhammer, N.; Vidal, V.; Vidal, V.S.; Chabaud, V.; Jesson, B.; Kwiatkowski, F.; Bignon, Y.-J. Use of Gene Expression Profiles of Peripheral Blood Lymphocytes to Distinguish BRCA1 Mutation Carriers in High Risk Breast Cancer Families. Cancer Inform. 2009, 7, 41–56. [Google Scholar] [CrossRef] [PubMed]
  30. Jamard, E.; Volard, B.; Dugué, A.E.; Legros, A.; Leconte, A.; Clarisse, B.; Davy, G.; Polycarpe, F.; Dugast, C.; Abadie, C.; et al. BRCA1 Allele-Specific Expression in Genetic Predisposed Breast/Ovarian Cancer. Fam. Cancer 2017, 16, 167–171. [Google Scholar] [CrossRef]
  31. Sedic, M.; Skibinski, A.; Brown, N.; Gallardo, M.; Mulligan, P.; Martinez, P.; Keller, P.J.; Glover, E.; Richardson, A.L.; Cowan, J.; et al. Haploinsufficiency for BRCA1 Leads to Cell-Type-Specific Genomic Instability and Premature Senescence. Nat. Commun. 2015, 6, 7505. [Google Scholar] [CrossRef]
  32. Cousineau, I.; Belmaaza, A. BRCA1 Haploinsufficiency, but Not Heterozygosity for a BRCA1-Truncating Mutation, Deregulates Homologous Recombination. Cell Cycle 2007, 6, 962–971. [Google Scholar] [CrossRef] [PubMed]
  33. Xu, W.; Ou, W.; Feng, Y.; Xu, Q.; Yang, Y.; Cui, L.; Du, P. Genetic Compensation Response Could Exist in Colorectal Cancer: UPF3A Upregulates the Oncogenic Homologue Gene SRSF3 Expression Corresponding to SRSF6 to Promote Colorectal Cancer Metastasis. J. Gastroenterol. Hepatol. 2023, 38, 634–647. [Google Scholar] [CrossRef] [PubMed]
  34. El-Brolosy, M.A.; Stainier, D.Y.R. Genetic Compensation: A Phenomenon in Search of Mechanisms. PLoS Genet. 2017, 13, e1006780. [Google Scholar] [CrossRef] [PubMed]
Scheme 1. Schematic description of experimental design. (A) Total RNA was extracted from peripheral blood lymphocytes obtained from 20 healthy individuals who were homozygous for the BRCA1WT allele, corresponding to the control group (C); 10 healthy individuals who were heterozygous for the BRCA1Δ9–12 variant (HH); and 10 cancer patients who were heterozygous for the BRCA1Δ9–12 variant (CaH), the last two as study groups. (B) A gene expression analysis between alleles in the three groups was performed with a real-time PCR (RT–qPCR) and ΔCT analysis. (C) A transcription variant analysis of the alleles of the three groups was performed using the Oxford Nanopore target sequencing. Each individual was marked with a barcode, allowing for tracking in each sample. (D) Isoform profiling of the expression pattern in the studied groups.
Scheme 1. Schematic description of experimental design. (A) Total RNA was extracted from peripheral blood lymphocytes obtained from 20 healthy individuals who were homozygous for the BRCA1WT allele, corresponding to the control group (C); 10 healthy individuals who were heterozygous for the BRCA1Δ9–12 variant (HH); and 10 cancer patients who were heterozygous for the BRCA1Δ9–12 variant (CaH), the last two as study groups. (B) A gene expression analysis between alleles in the three groups was performed with a real-time PCR (RT–qPCR) and ΔCT analysis. (C) A transcription variant analysis of the alleles of the three groups was performed using the Oxford Nanopore target sequencing. Each individual was marked with a barcode, allowing for tracking in each sample. (D) Isoform profiling of the expression pattern in the studied groups.
Ijms 25 06773 sch001
Figure 1. The expression of the BRCA1WT allele is increased in comparison to that of the BRCA1Δ9–12 allele in healthy and cancer-heterozygous patients. (A) A depiction of the 24 exons composing the BRCA1WT transcript (green). Primers for the amplification of the wild-type product recognize the junction of exons 6 and 7 (forward) and exon 10 (reverse), as depicted with yellow squares. RT–PCR produced a 210 bp amplification product of the wild-type transcript. (B) The deletion of exons 9–12 of the BRCA1 gene (green) leads to a mutant BRCA1Δ9–12 transcript (yellow). The primers used for amplification of the mutant transcript recognized the junction of exons 6 and 7 (forward) and the junction of exons 8 and 13 (reverse). RT–PCR produced a 127 bp amplification product of the BRCA1Δ9–12 transcript. (C) Endpoint PCR amplicons from 10 control samples (C1–C10) amplified with a mixture of one forward primer and two different reverse primers. The presence of a single 210 bp band indicates homozygosity for the BRCA1WT allele. GAPDH was used as a housekeeping control gene. (D) Endpoint PCR amplicons of samples from five patients with cancer (CA1–CA5) and five healthy heterozygotes (P1–P5) amplified with a mixture of one forward primer and two different reverse primers. The presence of two bands at 210 bp and 127 bp indicates heterozygosity for BRCA1 and amplification of the BRCA1WT and BRCA1Δ9–12 transcripts. GAPDH was used as a housekeeping control gene; GAPDH controls from CaH and HH samples were run on different gels. (E) Sanger sequencing electropherogram of a BRCA1WT control sample, corroborating the joining of exons 8 and 9 in the amplification product. (F) Sanger sequencing electropherogram of a BRCA1Δ9–12 sample, corroborating the joining of exons 8 and 13 in the amplification product. (G) A comparison of the differential expression between the BRCA1 alleles among the three groups of samples evaluated by RT–PCR and analyzed using ∆CT. There were differences between all the box data of the groups except for the BRCA1WT allele in the HH group, which showed expression levels similar to those of the homozygous control group. All bar comparisons were performed between alleles, not between study groups. The range of differences is indicated by the following values: * < 0.05, *** < 0.001. In the case of no difference, ns—not significant (p > 0.05) was used.
Figure 1. The expression of the BRCA1WT allele is increased in comparison to that of the BRCA1Δ9–12 allele in healthy and cancer-heterozygous patients. (A) A depiction of the 24 exons composing the BRCA1WT transcript (green). Primers for the amplification of the wild-type product recognize the junction of exons 6 and 7 (forward) and exon 10 (reverse), as depicted with yellow squares. RT–PCR produced a 210 bp amplification product of the wild-type transcript. (B) The deletion of exons 9–12 of the BRCA1 gene (green) leads to a mutant BRCA1Δ9–12 transcript (yellow). The primers used for amplification of the mutant transcript recognized the junction of exons 6 and 7 (forward) and the junction of exons 8 and 13 (reverse). RT–PCR produced a 127 bp amplification product of the BRCA1Δ9–12 transcript. (C) Endpoint PCR amplicons from 10 control samples (C1–C10) amplified with a mixture of one forward primer and two different reverse primers. The presence of a single 210 bp band indicates homozygosity for the BRCA1WT allele. GAPDH was used as a housekeeping control gene. (D) Endpoint PCR amplicons of samples from five patients with cancer (CA1–CA5) and five healthy heterozygotes (P1–P5) amplified with a mixture of one forward primer and two different reverse primers. The presence of two bands at 210 bp and 127 bp indicates heterozygosity for BRCA1 and amplification of the BRCA1WT and BRCA1Δ9–12 transcripts. GAPDH was used as a housekeeping control gene; GAPDH controls from CaH and HH samples were run on different gels. (E) Sanger sequencing electropherogram of a BRCA1WT control sample, corroborating the joining of exons 8 and 9 in the amplification product. (F) Sanger sequencing electropherogram of a BRCA1Δ9–12 sample, corroborating the joining of exons 8 and 13 in the amplification product. (G) A comparison of the differential expression between the BRCA1 alleles among the three groups of samples evaluated by RT–PCR and analyzed using ∆CT. There were differences between all the box data of the groups except for the BRCA1WT allele in the HH group, which showed expression levels similar to those of the homozygous control group. All bar comparisons were performed between alleles, not between study groups. The range of differences is indicated by the following values: * < 0.05, *** < 0.001. In the case of no difference, ns—not significant (p > 0.05) was used.
Ijms 25 06773 g001
Figure 2. Differential expression of both BRCA1 transcripts is associated with recurrence in patients. The relative expression of the samples (n = 9) was assessed using the ΔCT method, and the differences between groups were evaluated using the Wilcoxon test in conjunction with a Bonferroni adjustment. All analyses were performed using R; statistical significance was set at p < 0.05. Statistically significant differences are marked with * p < 0.05. A red square indicates where an overexpression is observed. (A) Association between clinical characteristics and relative expression of BRCA1WT allele. * p = 0.032 is noted in the red square. (B) Association between clinical characteristics and relative expression of BRCA1Δ9–12 allele variant. * p = 0.016 is noted in the red square.
Figure 2. Differential expression of both BRCA1 transcripts is associated with recurrence in patients. The relative expression of the samples (n = 9) was assessed using the ΔCT method, and the differences between groups were evaluated using the Wilcoxon test in conjunction with a Bonferroni adjustment. All analyses were performed using R; statistical significance was set at p < 0.05. Statistically significant differences are marked with * p < 0.05. A red square indicates where an overexpression is observed. (A) Association between clinical characteristics and relative expression of BRCA1WT allele. * p = 0.032 is noted in the red square. (B) Association between clinical characteristics and relative expression of BRCA1Δ9–12 allele variant. * p = 0.016 is noted in the red square.
Ijms 25 06773 g002
Figure 3. Complete amplification of BRCA1 variants for nanopore sequencing. (A) The location of the BRCA1WT transcript (green) and its exons are indicated by primers targeting exons 2 and 24. The BRCA1WT sequence product corresponds to 5802 bp. (B) The BRCA1Δ9–12 transcript (yellow) with the deleted exons shows the same primer localization at exons 2 and 24, with a decreased length of the amplified product of 2164 bp. (C) An endpoint PCR of the control samples (representative samples C3 and C4), healthy heterozygous samples (representative sample P6), and cancer patient samples (representative sample CA9) showing the transcripts and isoforms amplified with the primers described in A and B, with emphasis on the 5802 bp and 2164 bp amplicons. (D) Sanger sequencing corroboration of the BRCA1WT sample C3 sequence (left) from the purified band corresponding to an expected size of 5802 bp. The chromatogram of the BRCA1WT amplification showed the splicing of exons 8 and 9. Sanger sequencing of the sequences of BRCA1WT and BRCA1Δ9–12 in the sample CA9 (right) of the purified bands corresponding to expected sizes of 5802 bp and 2164 bp, confirming in the chromatograms the join of exons 8 and 9 in the BRCA1WT and exons 8 and 13 in BRCA1Δ9–12 amplicons.
Figure 3. Complete amplification of BRCA1 variants for nanopore sequencing. (A) The location of the BRCA1WT transcript (green) and its exons are indicated by primers targeting exons 2 and 24. The BRCA1WT sequence product corresponds to 5802 bp. (B) The BRCA1Δ9–12 transcript (yellow) with the deleted exons shows the same primer localization at exons 2 and 24, with a decreased length of the amplified product of 2164 bp. (C) An endpoint PCR of the control samples (representative samples C3 and C4), healthy heterozygous samples (representative sample P6), and cancer patient samples (representative sample CA9) showing the transcripts and isoforms amplified with the primers described in A and B, with emphasis on the 5802 bp and 2164 bp amplicons. (D) Sanger sequencing corroboration of the BRCA1WT sample C3 sequence (left) from the purified band corresponding to an expected size of 5802 bp. The chromatogram of the BRCA1WT amplification showed the splicing of exons 8 and 9. Sanger sequencing of the sequences of BRCA1WT and BRCA1Δ9–12 in the sample CA9 (right) of the purified bands corresponding to expected sizes of 5802 bp and 2164 bp, confirming in the chromatograms the join of exons 8 and 9 in the BRCA1WT and exons 8 and 13 in BRCA1Δ9–12 amplicons.
Ijms 25 06773 g003
Figure 4. Identification of differential isoform expression. The samples are color-coded as Ctrl (green), HH (yellow), and CaH (red). (A) Principal component analysis (PCA) of WT and 9–12 isoform expression between Ctr, HH, and CaH samples. The plot shows the two principal components. The Ctr samples differentiated from the HH and CaH samples on the PC1 axis with a greater variance (35.49%). (B) Principal component analysis in 3D (PCA 3D), where the samples are visualized in the plot with a third axis, PC3, where HH samples aggregate in the positive values of the axis, separating from CaH samples, which aggregate in the negative values of the axis. The P7 sample can be observed separately from the Ctr group on the PC2 axis. (C) Heatmap of dissimilarity matrix for clustering BRCA1 isoform expression. The color intensity on the map corresponds to the Euclidean distance of samples based on normalized counts of the BRCA1 isoforms, where white indicates high similarity and dark blue indicates low similarity in the isoform’s expression across individuals. (D) Heatmap of isoform z scores of the top differentially expressed isoforms between groups.
Figure 4. Identification of differential isoform expression. The samples are color-coded as Ctrl (green), HH (yellow), and CaH (red). (A) Principal component analysis (PCA) of WT and 9–12 isoform expression between Ctr, HH, and CaH samples. The plot shows the two principal components. The Ctr samples differentiated from the HH and CaH samples on the PC1 axis with a greater variance (35.49%). (B) Principal component analysis in 3D (PCA 3D), where the samples are visualized in the plot with a third axis, PC3, where HH samples aggregate in the positive values of the axis, separating from CaH samples, which aggregate in the negative values of the axis. The P7 sample can be observed separately from the Ctr group on the PC2 axis. (C) Heatmap of dissimilarity matrix for clustering BRCA1 isoform expression. The color intensity on the map corresponds to the Euclidean distance of samples based on normalized counts of the BRCA1 isoforms, where white indicates high similarity and dark blue indicates low similarity in the isoform’s expression across individuals. (D) Heatmap of isoform z scores of the top differentially expressed isoforms between groups.
Ijms 25 06773 g004
Table 1. Clinical characteristics of cancer patients carrying BRCA1Δ9–12 mutation.
Table 1. Clinical characteristics of cancer patients carrying BRCA1Δ9–12 mutation.
1° Primary Tumor2° Primary Tumor
PatientAge of Diagnosis
(Years)
Tumor LocalizationHistology *Clinical StageAge of Diagnosis
(Years)
Tumor LocalizationHistology *Clinical StageDisease-Free Survival (Months)OlaparibRecurrenceRadiotherapy
CA149OvaryHGSCIV55ThyroidPTCIA7227cHepaticNo
CA236Left BreastIDC SBR TN ki67 40%IANA27NANANo
CA335Right BreastGrade 2 IDCIIB51Left BreastIDC SBR 8 TN ki67 30%IA90NANAYes
CA453OvaryHGSCIV56ThyroidPTCIA6060cRetroperitoneumNo
CA543Right BreastIDC SBR 9 TN ki67 80%IIANA27NANAYes
CA630Left BreastIDC SBR 9 RH+ Her2+ ki 67 15%IIANA115NANANo
CA744OvaryHGSCIIICNA3636cHepaticNo
CA835Left BreastILC TNEX39Right BreastEXEXno3cCervicalYes
CA949OvaryHGSCICNA32NoInguinal lymph nodeNo
CA1032Right BreastIDC SBR 9 TNIA40OvaryHGSCIIIB868cRetroperitoneum and BreastYes
HGSC—high-grade serous carcinoma; IDC—infiltrating ductal carcinoma; PTC—papillary thyroid carcinoma; SBR—Scarff–Bloom–Richardson grading; TN—triple negative; HR—hormonal receptor; Her2+—human epidermal growth factor receptor protein positive; ILC—infiltrating lobular carcinoma; NA—not applicable; EX—external diagnosis/treatment. * The availability of clinical information in records has led to limitations in its presentation. No clinical information was available.
Table 2. Primers used to amplify complete and partial BRCA1 cDNA sequences.
Table 2. Primers used to amplify complete and partial BRCA1 cDNA sequences.
Protocol 1
SymbolSequence 5′-3′BasesAmplicon (bp)Transcript amplificated
F BRCA1 Ex.6-7TCC TTG CAG GAA ACC AGT CT20210Partial BRCA1WT
R BRCA1 Ex.10TTC ATC CCT GGT TCC TTG AG20
F BRCA1 Ex.6-7TCC TTG CAG GAA ACC AGT CT20127Partial BRCA1Δ9–12
R BRCA1 Ex.8-13CAT GGT ATC CCT CTG CCA ATT21
Protocol 2
SymbolSequence 5′-3′BasesAmplicon (bp)Transcript amplificated
F BRCA1 Jong ex1GCG CGG GAA TTA CAG ATA AA205802Complete BRCA1WT
R BRCA1 Jong ex24AAG CTC ATT CTT GGG GTC CT202164Complete BRCA1Δ9–12
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Dominguez-Ortiz, J.; Álvarez-Gómez, R.M.; Montiel-Manríquez, R.; Cedro-Tanda, A.; Alcaraz, N.; Castro-Hernández, C.; Bautista-Hinojosa, L.; Contreras-Espinosa, L.; Torres-Maldonado, L.; Fragoso-Ontiveros, V.; et al. A Molecular Characterization of the Allelic Expression of the BRCA1 Founder Δ9–12 Pathogenic Variant and Its Potential Clinical Relevance in Hereditary Cancer. Int. J. Mol. Sci. 2024, 25, 6773. https://doi.org/10.3390/ijms25126773

AMA Style

Dominguez-Ortiz J, Álvarez-Gómez RM, Montiel-Manríquez R, Cedro-Tanda A, Alcaraz N, Castro-Hernández C, Bautista-Hinojosa L, Contreras-Espinosa L, Torres-Maldonado L, Fragoso-Ontiveros V, et al. A Molecular Characterization of the Allelic Expression of the BRCA1 Founder Δ9–12 Pathogenic Variant and Its Potential Clinical Relevance in Hereditary Cancer. International Journal of Molecular Sciences. 2024; 25(12):6773. https://doi.org/10.3390/ijms25126773

Chicago/Turabian Style

Dominguez-Ortiz, Julieta, Rosa M. Álvarez-Gómez, Rogelio Montiel-Manríquez, Alberto Cedro-Tanda, Nicolás Alcaraz, Clementina Castro-Hernández, Luis Bautista-Hinojosa, Laura Contreras-Espinosa, Leda Torres-Maldonado, Verónica Fragoso-Ontiveros, and et al. 2024. "A Molecular Characterization of the Allelic Expression of the BRCA1 Founder Δ9–12 Pathogenic Variant and Its Potential Clinical Relevance in Hereditary Cancer" International Journal of Molecular Sciences 25, no. 12: 6773. https://doi.org/10.3390/ijms25126773

APA Style

Dominguez-Ortiz, J., Álvarez-Gómez, R. M., Montiel-Manríquez, R., Cedro-Tanda, A., Alcaraz, N., Castro-Hernández, C., Bautista-Hinojosa, L., Contreras-Espinosa, L., Torres-Maldonado, L., Fragoso-Ontiveros, V., Sánchez-Contreras, Y., González-Barrios, R., Fuente-Hernández, M. A. D. l., Mejía-Aguayo, M. d. l. L., Juárez-Figueroa, U., Padua-Bracho, A., Sosa-León, R., Obregon-Serrano, G., Vidal-Millán, S., ... Herrera-Montalvo, L. A. (2024). A Molecular Characterization of the Allelic Expression of the BRCA1 Founder Δ9–12 Pathogenic Variant and Its Potential Clinical Relevance in Hereditary Cancer. International Journal of Molecular Sciences, 25(12), 6773. https://doi.org/10.3390/ijms25126773

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop