- HPV prevalence in primary cervical lesions:
HPV prevalence was found in 88.8% (24/27) and 98.78% (81/82) of CIN and CACX samples respectively [Fig S2]. Out of HPV positive CIN samples, major samples were HPV 16 positive (87.5%; 21/24) followed by HPV 18 (8.33%; 2/24) and HPV 16/18 mixed positive (4.16%; 1/24). Among the HPV positive CACX samples, 82.71% (67/81) and 11.11% (9/81) of the samples were HPV 16 and HPV 18 positive respectively and the remaining 6.17% (5/81) were both HPV 16/18 mixed positive [Fig S2].
2. Expression profile analysis of FA-BRCA pathway genes:
A. In publicly available GEO dataset:
Four publicly available datasets viz. the GSE6722, GSE39001, GSE6791 and GSE9750 from GEO NCBI database were evaluated for expression microarray analysis. Where the different stages of CACX samples were compared with respect to normal samples, in order to understand the role of this pathway in the progression of this disease. In tumour samples, significant and differential downregulation of BRCA1 and BRCA2 was observed in majority of the datasets [Fig. 1; A-B]. Whereas FANCC and FANCD2 genes were found to be significantly down-regulated in tumour samples of GSE67522/GSE39001 and GSE39001/GSE6791 datasets respectively [Fig. 1; A-B]. This indicates that these key regulatory genes may play an important role in the early steps of carcinogenesis from normal to tumour tissues.
B. mRNA expression profile in the samples:
The mRNA expression of the key regulatory genes of the FA-BRCA pathway was comparable (BRCA1/FANCD2) or slightly decreased (BRCA2/FANCC) in CIN followed by gradual down-regulation with progression of the disease [Fig. 1; C-F]. Interestingly, significant down-regulation of FANCC (p= 0.019) was found in stage-I/II samples, whereas, BRCA1 and FANCD2 expression were significantly (p= 0.036, 0.007) reduced in stage-III/IV compared to stage-I/II. Therefore, under expression of key regulatory genes of FA-BRCA pathway, as found from the GEO datasets, was well in concordance with our sample pool. The expression profile of the genes in the CACX cell lines was also checked and seemed to be comparable with the tumour samples [Fig. 1; C-F].
C. Immunohistochemical expression profile in the samples:
The mRNA expression pattern was next correlated with the protein expression of the genes by immunohistochemistry [Fig. 2; A-D]. The prevalence of high/medium (H/M) predominant nuclear expression of BRCA1, BRCA2 and FANCD2 proteins was high in basal/parabasal layers of cervical epithelium with 81.8%, 86.36% and 68.18% respectively followed by their gradual decrease in spinous layer (63.6%, 68.18% and 40.90%) [Fig. 2; A-B, D]. Whereas, high frequency (86.36%) of FANCC expression (H/M) was seen throughout the entire cervical epithelium [Fig. 2; C].
On contrary to Normal cervical epithelium, the nuclear expression of BRCA1 protein was gradually decreased from CIN (66.60%) to stage-I/II (57.69%) and significantly reduced (p=0.018) in later invasive stage (stage-III/IV) (16.66%) [Fig. 2; A]. Interestingly, significant downregulation (p=0.030) of FANCC was preferably noticed in the early invasive stage (stage-I/II) (34.61%) of cervical lesions and maintained in stage-III/IV as well (25%) [Fig. 2; C], this phenomenon may insight that role of downregulation of the key regulatory genes through the different stages of the disease is needed according to the need of tumorigenesis. Furthermore, similar trend of comparable protein expression of both BRCA2 and FANCD2 was observed in both CIN and stage-I/II (around 50%) and reduced in stage- III/IV (around 17%), [Fig. 2; B, D]. Apart from this, protein expression pattern of all the key regulatory genes showed strong concordance with their respective mRNA expression as well [Table S3A].
Among the invasive lesions, around 59.09% of the samples showed co-downregulation of more than one protein of FA-BRCA pathway, signifying the cumulative role of down-regulation of these genes in the development of CACX [Fig. 2; E].
3. Promoter methylation analysis of FA-BRCA pathway genes:
A. In publicly available TCGA-CACX dataset:
In the TCGA cohort, promoter methylation is represented in the form of β-values which ranges from 0 (unmethylated) to 1 (fully methylated). The promoter methylation frequencies of BRCA1, BRCA2 and FANCC were significantly (p= 0.0004, 0.003, 0.004) high in primary tumour compared to normal tissue. On the other hand, FANCD2 showed very slight change in promoter methylation frequency in tumour tissue compared to normal. [Fig 3; A].
B. Qualitative promoter methylation profile in the samples:
The TCGA dataset could not explain protein expression levels of the genes in different layers of normal cervical epithelium. So, we tried to address this problem by undertaking promoter methylation analysis in composite (total), laser-guided micro-dissected basal-parabasal (stem cell rich, undifferentiated cells) and spinous layers (more terminally differentiated cells) of normal epithelium and in cervical lesions. The strategies of MSRA have been shown in Fig S1.
Infrequent promoter methylation of the FA-BRCA genes in normal cervical epithelium (total) was seen in the following order- FANCD2 (17.14%)>BRCA1-FANCC (14.28%) and > BRCA2 (11.42%) [Fig 3; B-E (i-ii)]. However, in respective micro-dissected samples, the methylation frequencies of BRCA1, BRCA2 and FANCD2 were low in basal-parabasal layer with 11.42%, 8.75% and 11.42% respectively than the spinous layer (17.14%, 11.42% and 20.00%) [Fig 3; B-C, E (i-ii)]. Unlike these genes, uniform promoter methylation of FANCC was found in both the layers (14.28%) of normal cervical epithelium [Fig 3; D (i-ii)].
In cervical lesions, promoter methylation of all the genes was gradually increased from CIN to subsequent clinical stages compared to the basal-parabasal layer of normal tissue [Fig 3; B-D (i-ii)]. Noticeably, the methylation frequencies of FANCC and BRCA2/FANCD2 genes were significantly high (p=0.027-0.039) in early (stage-I/II) and later (stage-III/IV) invasive lesions respectively [Fig 3; C-E (i-ii)]. The promoter methylation pattern of all the key regulatory genes showed strong concordance with their respective protein expression as well [Table S3B]. This event indicates the plausible role of promoter methylation of the genes for down-regulation of respective protein expressions during the progression of the disease.
C. Quantitative promoter methylation profile in the samples:
The methylation pattern of the genes was further validated by quantitative MSRA (q-MSRA), as it gives us an idea of the extent of promoter methylation percentage in each sample [24]. Dot-plot obtained from the q-MSRA data revealed the sample-wise details of the changes in promoter methylation events during tumorigenesis from normal cervical epithelium (basal-parabasal and spinous layer). In q-MSRA analysis, the result showed that hypo-methylation frequencies of FA-BRCA pathway genes in basal-parabasal and spinous layers showed concordance with their qualitative MSRA analysis followed by a gradual decrease in CIN and subsequent clinical stages [Fig 3; F (i-iv)]. [lower the percentage of hypomethylation indicates higher frequency of methylation in the samples].
4. Deletion analysis of FA-BRCA pathway genes:
To evaluate the molecular mechanism behind the under expression of the genes, deletion analysis was performed. Deletion frequency of BRCA1 and BRCA2 was gradually increased from CIN (33.33%, 25.92%) to stage-I/II (46.15%, 38.46%) and became significantly (p= 0.024, 0.030) high at stage-III/IV (73.33%, 63.33%) [Fig 4; A-B (i-ii)]. However, infrequent deletion frequency of FANCC in 14.81% of CIN samples became significantly (p= 0.043) increased in early invasive stages (36.53%) and remained comparable in stage-III/IV (40%) as well. Compare to FANCC, deletion frequency of FANCD2 was comparatively high on pre-malignment CIN samples (29.63%) and gradually increased in stage-I/II (48.07%) and stage-III/IV (56.66%) [Fig 4; C-D (i-ii)]. In addition to deletion, infrequent MA of BRCA1 and BRCA2 was noticed in around 11.11% and 7.40% of CIN samples followed by their gradual increase in stage-I/II (17.37%, 13.46%) and became significant (p= 0.048, 0.032) at stage-III/IV (36.66%, 33.33%) [Fig 4; A-B (i-ii)]. This phenomenon insight the role of genomic instability occurred due to alteration of any of the genes of FA-BRCA pathway actually induces the development of early steps of carcinogenesis.
5. Clinico-pathological association with alterations of FA-BRCA pathway genes:
In totality, overall alterations (either methylation and/or deletion) of the genes of CIN samples were seen in the following order- FANCD2 (55.55%)> BRCA2 (44.44%)> BRCA1/FANCC (37.03%) with gradual increase with progression of the disease [Fig 4; E (i-iv)]. Interestingly, overall alterations in FANCC and BRCA2 were significantly high (p= 0.0009, 0.002) in early and later invasive stages respectively [Fig 4; E (ii-iii)]. This event is signifying that alteration of any of the alleles may actually play an anchoring role in regulating the reduced expression of these genes.
On the other hand, co-alterations of the genes (both methylation and deletion) were also gradually increased with tumour development [Fig 4; E (i-iv)]. The result showed that the co-alteration frequencies of both BRCA1 and FANCD2 were significantly raised (p= 0.030, 0.012) at stage-III/IV samples.
In 5-year follow-up study (Kaplan-Meier analysis) revealed patients with BRCA1-deletion in tumours had poor patient outcomes [Fig 5; A]. In addition, patients having overall alterations (Methylation and/or Deletion) of either one of the key regulatory genes like- BRCA1/BRCA2/FANCC/FANCD2 showed the worst prognosis among patients [Fig 5; B-E]. Similarly, patients with co-alteration of BRCA1-BRCA2-FANCC-FANCD2 together had poor survival, suggesting the prognostic importance of the FA-BRCA pathway genes with disease outcome [Fig 5; F].
6. Prevalence of promoter methylation of FA-BRCA pathway genes in plasma DNA of CACX patients and its association with survival:
To evaluate the usefulness of promoter methylation of the FA-BRCA pathway genes in plasma DNA, 36 primary cervical lesions were analysed. Details of therapy (CCRT therapy) and the survival status of the respective patients were also recorded in their follow-up data. In our study, methylation frequencies of the genes in plasma DNA derived from circulating tumour cells (CTCs) were in the following order- BRCA1-BRCA2 (16.66%)> FANCD2 (13.88%)> FANCC (11.11%) [Fig 6; A-D]. Noticeably, it was evident that the presence of promoter methylation of the respective genes was more prevalent in later invasive lesions (stage-III/IV) compare to early stage (stage- II). Prevalence of promoter methylation of the genes (BRCA1/FANCC/FANCD2) in plasma DNA also showed significant (p= 0.029- 0.042) concordance of their methylation frequency in respective primary tumours. [Fig 6; E].
Interestingly, patients having methylation of any one of the genes in plasma showed poor patient outcome even after the therapeutic intervention [Fig 6; E].
7. Effect of cisplatin-treatment in the molecular profile of FA-BRCA pathway genes in CACX cell lines:
After analysing the poor patient survival due to the non-responsiveness of the patients to the treatment in our follow-up data as discussed previously [Fig 6], our aim was to understand the molecular mechanism of the FA-BRCA pathway genes in disease recurrence among the CACX patients. For this purpose, the CACX cell lines were treated with different concentrations of cisplatin and then the molecular profile (expression/methylation profile) of the target genes was evaluated in our in-vitro study.
A. Cisplatin sensitivity of the CACX cell lines:
The CACX cell lines showed different sensitivity towards cisplatin treatment. Hela cells were seen to be more sensitive to cisplatin than SiHa cells [Unpublished Data].
B. Effect of cisplatin on the expression of the genes in CACX cell lines:
After treating the cell lines with different concentrations of cisplatin for 48 hrs, mRNA expression of FA-BRCA pathway genes was analysed. It was evident that mRNA expression was slightly increased at IC30 and became significantly high at a gradual increment of the drug concentrations of IC50 and IC70 in both of the cell lines [Fig 7; A (i-iv)].
In western blot analysis, the protein expression also gradually increased with the increasing drug concentrations and predominantly raised at IC50 and IC70 of cisplatin in both of the cell lines, showing concordance with their mRNA expression in both the cell lines [Fig 7; B (i-v)], indicating transcriptional activation of the genes after cisplatin treatment in the cell lines.
C. Effect of cisplatin on promoter methylation of the genes in CACX cell lines:
To understand the molecular mechanism behind the upregulation of the FA-BRCA pathway genes in cisplatin treated cells, their promoter methylation status was analysed by q-MSRA. The hypomethylation frequency of the genes after cisplatin treatment in the cell lines was gradually increased (higher the value of hypomethylation percentage indicates lower the promoter methylation and vice versa) with drug concentrations and became significantly high at IC50 and IC70 concentrations, showing concordance with expression of the genes [Fig 7; C (i-v)].
8. Effect of cisplatin on DNMT1 expression in CACX cell lines:
To find out the plausible reason for reduced promoter methylation of the genes in cisplatin-treated cell lines, DNMT1 expression (both mRNA and protein) was analysed, due to its importance in the maintenance of tissue-specific methylation [Klein et al, 2011]. It was evident that m-RNA expression of DNMT1 was prominently reduced at IC50 and IC70 concentration of cisplatin in SiHa and only at IC70 in HeLa cell line respectively [Fig 8; A].
In the western blot analysis of DNMT1, gradual decrease in protein expression of DNMT1 was seen with increasing drug concentrations and showed concordance with its mRNA expression [Fig 8; B-C].
9. Effect of cisplatin on immunocytochemical analysis of ϒ-H2X foci in CACX cells:
For functional evaluation of the up-regulated FA-BRCA pathway genes in cisplatin-treated CACX cell lines, the ϒ-H2X nuclear foci formation was analysed after treatment with different concentrations of cisplatin for 48 hrs in CACX cell lines [Fig S3; A-B]. As ϒ-H2X is an important marker of DNA damage [12], our result revealed that there was an increase in the number of ϒ-H2X foci at IC30 (32-37 foci/cell) of the drug compared to untreated control followed by a gradual decrease at IC50 (21-27 foci/cell) and IC70 (12-15 foci/cell) concentrations [Fig S3; C], indicating activation of the DNA repair mechanism in the residual tumour cells in the system.