Search tips
Search criteria 


Logo of bmcgenoBioMed Centralsearchsubmit a manuscriptregisterthis articleBMC Genomics
BMC Genomics. 2008; 9: 64.
Published online 2008 February 4. doi:  10.1186/1471-2164-9-64
PMCID: PMC2277413

Up regulation in gene expression of chromatin remodelling factors in cervical intraepithelial neoplasia



The highest rates of cervical cancer are found in developing countries. Frontline monitoring has reduced these rates in developed countries and present day screening programs primarily identify precancerous lesions termed cervical intraepithelial neoplasias (CIN). CIN lesions described as mild dysplasia (CIN I) are likely to spontaneously regress while CIN III lesions (severe dysplasia) are likely to progress if untreated. Thoughtful consideration of gene expression changes paralleling the progressive pre invasive neoplastic development will yield insight into the key casual events involved in cervical cancer development.


In this study, we have identified gene expression changes across 16 cervical cases (CIN I, CIN II, CIN III and normal cervical epithelium) using the unbiased long serial analysis of gene expression (L-SAGE) method. The 16 L-SAGE libraries were sequenced to the level of 2,481,387 tags, creating the largest SAGE data collection for cervical tissue worldwide. We have identified 222 genes differentially expressed between normal cervical tissue and CIN III. Many of these genes influence biological functions characteristic of cancer, such as cell death, cell growth/proliferation and cellular movement. Evaluation of these genes through network interactions identified multiple candidates that influence regulation of cellular transcription through chromatin remodelling (SMARCC1, NCOR1, MRFAP1 and MORF4L2). Further, these expression events are focused at the critical junction in disease development of moderate dysplasia (CIN II) indicating a role for chromatin remodelling as part of cervical cancer development.


We have created a valuable publically available resource for the study of gene expression in precancerous cervical lesions. Our results indicate deregulation of the chromatin remodelling complex components and its influencing factors occur in the development of CIN lesions. The increase in SWI/SNF stabilizing molecule SMARCC1 and other novel genes has not been previously illustrated as events in the early stages of dysplasia development and thus not only provides novel candidate markers for screening but a biological function for targeting treatment.


Cervical cancer affects approximately 500,000 women worldwide each year with highest rates in developing countries [1]. Cervical intraepithelial neoplasia (CIN) is a precursor lesion to cervical cancer and can be further subdivided into CIN I, CIN II and CIN III (mild, moderate and severe dysplasia, respectively). Most CIN I lesions spontaneously regress to normal however CIN III lesions are much more likely to progress to cervical cancer if left untreated [1]. The CIN I to CIN II junction may be critical in disease development.

Human Papillomavirus (HPV) is the recognized etiologic agent for cervical cancer however, alone it is not sufficient for invasive disease. HPV is detected in nearly all cervical cancers, 94% of CIN lesions and up to 46% of normal cervical epithelium [1]. Over 100 strains of HPV exist however HPV 16 and HPV 18 are considered highly virulent strains and account for the majority of cervical cancers [1,2].

The study of cervical cancer prevention has progressed impressively in the recent past. Widely implemented screening programs have resulted in 80% reduction of cervical cancer rates in North America within the past fifty years [3]. Cytological assessment is currently the frontline method for identifying precancerous cervical lesions however repeat evaluations can frequently be required due to low sensitivity [3]. Although vaccines against the most virulent strains of HPV have recently become available, vaccination is not an easily utilized resource for those countries most inflicted with the highest cervical cancer rates [3]. This is largely due to cost and that the administration of the vaccine occurs in three doses over six months to a pre-adolescent female population. Few countries have established programmes targeting healthcare to this population [4]. In addition, perceived social implications in developed countries regarding vaccinating girls at an early age are hindering widespread administration. Together, these social, clinical and genetic factors indicate that frontline monitoring will continue to play an important role in cervical cancer prevention and that improved methods and markers for detection are needed.

It is unclear that HPV alone is responsible for disease progression. A thorough understanding of genetic events in precancerous cervical intraepithelial neoplasia is required to both delineate important causal events in cervical cancer and to identify informative candidate biological markers. Gene expression of cervical tissue and changes in expression pattern have been the focus of several recent publications. Studies by Pérez-Plasencia et al and Shadeo et al both characterized the transcriptome of normal cervical epithelium using serial analysis of gene expression (SAGE) [5,6], Additionally, Gius et al reported changes in proproliferative/immunosuppression gene expression in CIN I lesions, as well as proangiogenic and proinvasive expression signatures that coincide with CIN II and CIN III, respectively [7].

In this study, we build upon our previous work in defining the normal cervical epithelial transcriptome and aim to identify genes differentially expressed between normal cervical epithelium and those precancerous lesions which are more apt to progress to cervical cancer if left without treatment (CIN III). In this study we have distinguished gene expression aberrations across mild/moderate dysplasia (CIN I, CIN II) in addition to CIN III and non cancerous (NC) cervical epithelium using an unbiased long serial analysis of gene expression (L-SAGE) method that simultaneously allows for the discovery of tags which map to HPV 16. In total, sixteen L-SAGE libraries were sequenced for a total of 2,481,387 tags, establishing the largest SAGE data collection for cervical tissue worldwide. Upon evaluation of expression differences between NC cervical epithelium and CIN lesions, we have identified two gene networks directly or indirectly involved in chromatin remodelling altered in expression in CIN III.


L-SAGE Libraries from Cervical Tissue

In this study 16 L-SAGE libraries were constructed and analyzed (Figure (Figure1).1). Libraries N1 to N4 were made from NC cervical tissue samples, M1 to M3 from CIN I samples, M4, to M6 from CIN II samples and C1, to C6 from CIN III samples. N1, N2, C1 and C2 were mined in a preliminary study characterizing normal cervical tissue [6]. Collectively, 2,481,387 useful tags were sequenced (Figure (Figure1).1). This data collection has been made publicly available at Gene Expression Omnibus, series accession number GSE7433 [8].

Figure 1
Summary of L-SAGE library analysis. Sixteen L-SAGE libraries were created and a total of 2,481,387 useful tags were sequenced. An average of 130,560 tags were sequenced per NC library (4), 155,807 per CIN I/CIN II library (6) and 170,718 per CIN III library ...

Early Stage Changes

The mean tag counts were compared between NC tissue samples (libraries N1-4) and mild/moderate stage (CIN I/II) samples (libraries M1-6) in order to identify tags differentially expressed in the early stages of neoplasia. One-hundred sixty-nine tags were identified to be differentially expressed according to our selection criteria as described in Methods [see Additional file 1]. Both increased and decreased expression is observed at comparable frequencies (75 tags increased, 94 tags decreased in CIN I/II) and 138 of these tags mapped to known genes including loci encoding hypothetical proteins.

The most commonly affected biological process in the early events of disease is DNA dependent regulation of transcription and transcription (11% and 8% of differentially expressed tags, as assessed by Onto-Express, respectively) (Figure (Figure2A)2A) [9,10].

Figure 2
Functional groupings of tags differentially expressed in SAGE libraries. (A) Genes differentially expressed between NC and CIN I/CINII. (B) Genes differentially expressed between CIN I/CIN II and CIN III. Categories are as described. The Unknown function ...

Late Stage Changes

The mean tag counts were compared between CINI/II tissue samples (libraries M1-6) and CIN III libraries (C1-6). This identified 128 tags differentially expressed in later stages of neoplasia [see Additional file 2]. Seventy-three tags were increased in CIN III while 55 were decreased. Overall 97 tags mapped to known genes. The major pathways affected by these genes were similar to those observed in early stage changes, DNA dependent regulation of transcription and transcription (8% and 7% as assessed by Onto-Express. respectively) (Figure (Figure2B)2B) [9,10]. Twelve tags mapped to the HPV 16 genome that could not be mapped to the human genome or transcripts using methods previously described [see Additional file 3]. [6]. In addition, 42 tags showed concurrent expression between all sixteen libraries [see Additional file 4].

Pathway Analysis of Overall Changes

The mean score of tags in the four NC libraries were compared to the mean score of tags in the six severe neoplasia libraries (CIN III). We identified 108 tags increased in frequency in CIN III and 138 tags decreased and overall observed 246 tags differentially expressed between NC tissue and CIN III [see Additional file 5]. Two hundred and forty tags mapped to 222 unique genes. Genes differentially expressed between normal and CIN III were evaluated for gene network associations using Ingenuity Pathway Analysis [11]. Biological functions most influenced by these genes include cell death, cell growth/proliferation and cellular movement (Figure (Figure3).3). The gene network with the most number of differentially expressed genes influenced cellular processes such as cell cycle and cell morphology and involved 17 of the differentially expressed genes (Table (Table1,1, Figure Figure4).4). The network with the higher magnitude gene expression changes influenced cell cycle and gene expression (Figure (Figure55).

Table 1
Genes differentially expressed in Networks A and B.
Figure 3
Biological functions targeted by gene expression changes between NC and CIN III Functional categories deemed as significant by Ingenuity Pathway Analysis software, are displayed along the x-axis. The y-axis displays the -(log) significance and the orange ...
Figure 4
Gene Network A. Green denotes loss in expression in CIN III while red indicates gain. Intensity of colour signifies magnitude of change. Solid lines indicated direct interactions while dashed lines indicate indirect interactions as described by Ingenuity ...
Figure 5
Gene Network B. Green denotes loss in expression in CIN III while red indicates gain. Intensity of colour signifies magnitude of change. Solid lines indicated direct interactions while dashed lines indicate indirect interactions as described by Ingenuity ...

Eight genes from these networks were selected for validation in an independent panel of cervical samples (Table (Table2).2). Cyclin-dependent kinase inhibitor 2B (CDKN2B), retinoblastoma-like 2(p130 or RBL2) and phosphatase and tensin homolog (PTEN) had lower counts in the CIN III derived libraries. Conversely, Mof4 family associated protein 1 (MRFAP1), mortality factor 4 like 2(MORF4L2), nuclear receptor co-repressor 1 (NCOR1), dihydrofolate reductase (DHFR) and SWI/SNF related, matrix associated, actin dependent regulator of chromatin, subfamily c, member 1(SMARCC1) showed higher tag counts in the CIN III relative to NC.

Table 2
Candidate genes in Networks A and B.

The validation panel consists of 22 samples representing NC, CIN I, CIN II and CIN III in addition to a tumor sample with paired normal tissue. Quantitative analysis of expression was determined for each gene on this panel by real-time PCR using TaqMan® Gene Expression Assays (Figure (Figure6).6). NCOR1 expression increased in nine of ten CIN I/II samples ranging from 1.3–6.0 fold increase relative to NC with a statistically significant overall all trend (p < 0.05). SMARCC1 expression increase was confirmed in six CIN I/II samples and four CIN III cases (Figure (Figure6).6). Similarly, we confirmed the increased expression of DHFR and MRFAP1 in the validation panel including in nine of ten CIN I/II cases and six of seven CIN III cases and the overall trend of increase for this gene was statistically significant (p < .05). MORF4L2 expression increase was confirmed in eight of ten CIN I/II and six of seven CIN III with statistically significant overall trend of increase (p < .05) in the validation panel. Although we were able to confirm decreased CDKN2B expression in the tumour, we were unable to confirm this in the earlier staged cases in the validation panel. Reduced expression of phosphatase and tensin homolog (PTEN) was confirmed in two CIN II and two CIN III cases [see Additional file 6].

Figure 6
Summary of validation panel quantitative PCR results of genes with altered expression in CIN III L-SAGE libraries. A panel of 22 new cases was investigated in triplicate for expression changes for eight genes. Samples are indicated as follows: CIN I (A-C), ...


The study of precancer lesions is essential in understanding the initiating events in cancer development and is of particular importance in cervical cancer as moderate and high grade lesions (CIN II and CINIII) are more likely to progress to cancer than those of low grade (CIN I). In this study, we have comprehensively evaluated gene expression changes in the precancerous stages of cervical cancer by comparing 16 cervical intraepithelial neoplasia and NC cervical specimens using the L-SAGE method of gene expression profile determination. We have sequenced a total of 2,481,387 tags, establishing the largest SAGE data collection for cervical tissue worldwide. CIN I, CIN II, CINIII as well as normal tissues are represented in this collection. In a preliminary study, we had identified Keratin 6A (KRT6A), carcinoembryonic antigen-related cell adhesion molecule 7 (CEACAM7), S100 calcium-binding protein A7 (S100A7) and small proline-rich protein 3 (SPRR33) to be highly expressed (>500 TPM) in NC cervical tissue in the N1 and N2 libraries [6]. Tags mapping to these genes are present at high levels in the additional normal libraries (N3 and N4) unique to this study with the exception of CECEAM7 which shows slightly reduced expression of 363 and 397 TPM (N3 and N4, respectively). The decreasing trend in Galectin 7 (LGALS7) and Gap junction protein A1 (GJA1) previously observed in CIN III when compared to normal mean values, continues to be true in this study (3.0 and 1.7 fold decrease in CIN III, respectively, p value not significant).

Early and Late Events

Of the tags differentially expressed between NC and mild/moderate dysplasia (169 tags), and between mild/moderate and severe dysplasia (128 tags); 25% fewer tags were altered in later stage disease. These differentially expressed tags most frequently mapped to genes involved in DNA dependent regulation of transcription. Interestingly, tags mapping to genes with biological functions unique to late stage changes include cell cycle (5% of tags), cell division (3%), immune (4%) and inflammatory (3%) response.

Biological Characteristics of Differentially Expressed Tags

Of the 246 tags that are altered in expression in CIN III relative to NC, cell death, cell growth and proliferation, cell movement, cell cycle and DNA replication, recombination and repair are the top five biological functions. These functions encompass characteristics frequently described as the hallmarks of cancer [12].

Tags mapping to seven unique genes showed a greater than ten fold increase in expression in CIN III libraries. A tag mapping to SEC13 homolog (SEC13L1) showed the greatest change as it was not present in any one of the four NC libraries. However, an average of 20 TPM were observed in CIN III libraries. Additional tags with greater than ten fold change mapped to MORF4L2, MRFAP1, WD repeat domain 18(WDR18), SMARCC1, eukaryotic translation elongation factor 1 gamma (EEF1G) and G protein-coupled receptor 180 (GPR180). SMARCC1 is a component of the chromatin remodeling complex SWI/SNF while both MRFAP1 and MORF4L2 belong to the MGF/MORF family of transcription factors involved in growth, cell senescence and are implicated indirectly in chromatin remodeling [13-18]. WDR18 is a member of the WD repeat protein family. Members of this gene family are involved in a variety of cellular processes including cell cycle progression, signal transduction, apoptosis, and gene regulation [19]. Tags mapping to this gene varied in counts between NC and CIN III by greater than twelve fold increase in the latter subgroup. EEF1G is a subunit of the elongation factor-1 complex involved in translation and GPR180 contains transmembrane domains and may play a role in vascular remodelling [20]. It is interesting to note that both homeobox B7 (HOXB7) and BH3 interacting domain death agonist (BID) are also increased in CIN III although to a lesser magnitude (seven and five fold, respectively). Many of the described genes may have the potential to influence processes characteristic of epithelial cancers such as, apoptosis (BID), angiogenesis (GPR180), proliferation (WDR18), and transcription influencing events (SMARCC1, MRFAP1, MORFL2). Functional assays are required to delineate their biological function in CIN and cervical cancer.

Network Analysis

Upon further investigation of the 246 differentially expressed tags between NC and CIN III, 240 mapped to 222 unique genes. When analyzed for gene network relationships, many of the genes targeted pathways influencing properties such as cell cycle, cell morphology, cancer related events and gene expression. In this study, we have identified a group of 28 genes that fall into two gene networks encompassing these properties. The first network (A) contains 15 of the differentially expressed genes while the second network (B) includes 13 of the 222 genes. Only cyclin-dependent kinase inhibitor 2b (CDKN2B) and kruppel-like factor 6 (KLF6) overlap between Network A and B. Network B includes three of the seven tags showing greater than 11 fold change and one of the five tags showing more than 20 fold decrease in NC when compared to CIN III. Genes from both networks were selected for validation in a new cervical tissue panel.

Network A

Cell cycle, cancer and cell morphology are the top functions influenced by genes in Network A. Sixteen genes differentially expressed between NC and CIN III are involved in this network including RBL2/p130, SMARCC1, NCOR1, PTEN, DHFR and CDKN2B.

RBL2/p130 is one of three members of the Retinoblastoma (Rb) gene family, others include RBL1/p107 and RB1 [21]. This family of proteins regulate cell cycle through the G1/S phase by sequestering E2F transcriptional regulators. Loss and mutations in RBL2/p130 have been linked to various cancers [22,23]. E7 from high risk HPV strains targets all member of the Rb family including RBL2/p130, for degradation thus releasing sequestered E2F and allowing for passage through the G1/S phase [24,25]. Zhang et al found that low risk HPV 6 E7 also has the capacity to bind RBL2/p130 although not with as high affinity and did not bind to the other member of the Rb family [25]. We observed a progressive reduction in RBL2/p130 transcripts from NC to mild/moderate dysplasia and severe dysplasia suggesting that a reduction in RBL2/p130 may also be regulated at the transcriptional level in cervical preneoplasia and may be the first Rb family gene event in the development of cervical intraepithelial neoplasia.

SMARCC1 is a member of the SWI/SNF family of genes. Members of this evolutionarily conserved gene family are proposed to regulate gene specific transcription with downstream effects in cell cycle progression [26]. SMARCC1 has recently been reported to directly interact with components of the ATP-dependent SWI/SNF chromatin remodelling complex including, SWI/SNF-related, matrix-associated, actin-dependent regulator of chromatin A4 (BRG1), BRG1-associated factor A (BAF60a) and SWI/SNF-related, matrix-associated, actin-dependent regulator of chromatin B1 (SNF5) [15]. The SWI/SNF chromatin remodelling complex consists of approximately ten components and is thought to regulate gene transcription by altering the surrounding chromatin structure. This complex can recruit both histone acetyltransferases (HATs) and histone deacetyltransferases (HDACs) in a gene specific manner thus influencing gene expression [27]. The SWI/SNF chromatin remodelling complex also has been implicated in hormone receptor signalling and growth [16]. SMARCC1 was characterized as a key regulator of core complex components (BRG1, BAF60a and SNF5) by positively influencing the stabilization of the complex proteins as opposed to regulation of their transcription [15]. In our study, there was no expression change in these core components. In contrast, we observed an 11 fold increase in tags mapping to SMARCC1 in CIN III libraries when compared to normal libraries suggesting that SWI/SNF chromatin remodelling complex stability may have a role in severe dysplasia development.

DHFR is involved in folate metabolism in eukaryotes and is essential for purine and thymidylate biosynthesis and thus DNA replication [28]. We observed over a 250% increase in both mild/moderate and severe dysplasia as compared to normal, indicating that this event is present in very early stage lesions. It has recently been reported that DHFR is subjected to RB mediated repression via SWI/SNF chromatin remodelling activity [29]. The increase in DHFR observed is in concordance with the expected downstream effects of decreased RB gene family repression.

NCOR1 is involved in transcription repression via chromatin condensation. It has been found to physically interact with members of the SWI/SNF complex, including SMARCC1 and core component BRG1 [30]. We observed greater than 300% increase in NCOR1 in severe and mild/moderate dysplasia relative to normal suggesting an increase in NCOR1 occurs very early in disease development (CIN II/CIN I).

Phosphatase and tensin homolog (PTEN), an established tumour suppressor gene, functions through the AKT/PKB signalling pathway [31]. We found PTEN to decrease in expression in CIN II/CIN I libraries by >3 fold compared to normal indicating that this may be an early event. CDKN2B negatively regulates CDK4 and CDK6 [32]. We found CDKN2B decreased in CIN III libraries by nearly 3.5 fold when compared to NC. Decreased expression of CDKN2B has not been found due to copy number loss or mutation in cervical cancer [32].

Genes from Network A were selected for validation in a new 14 sample set consisting NC, CIN I, CIN II and CIN III cervix tissue and one tumour and normal pair. A similar pattern of fold change is present in this panel for DHFR, SMARCC1 and NCOR1 in all three subgroups. Specifically, increased expression is observed in all three CIN II and two CIN III in the new cervical tissue panel, validating the disruption of genes involved in chromatin remodelling.

For PTEN, we observed a decrease in one of the three CIN III cases investigated which is consistent with Lee et al who found PTEN to be reduced in protein expression in only 10% of CIN III cervical cases and 18% of cervical cancers [33]. We also found CDKN2B to be over expressed in all CIN I, two CIN II and one CIN III however, greater than four fold decrease in cervical cancer (T1) was observed. It has been reported that cervical cancers and CIN III lesion have elevated levels of CDKN2B (84% and 79%, respectively) [34].

Network B

Fourteen genes differentially expressed between normal and CIN III targeted Network B including MRFAP1, MORFL2 and CDKN2B. Top functions influenced by Network B include cancer related factors, cell cycle and gene expression. The previously discussed gene CDKN2B overlaps with this network.

MRFAP1 binds with strong affinity to mortality factor 4 like 1 (MRG15) [17]. MRG15 is a member of the MRG/MORF family of genes. All members contain a MRG domain which has the capacity to bind multiple transcriptional regulators [18]. Members of this gene family are suggested to be involved in embryonic development, cell proliferation and senescence [35]. MRG15 is a component of the MRG associated factors 1 and 2 complexes (MAF1, MAF2, respectively) and thus has a putative role in chromatin remodelling [36]. MRG15 has been shown to specifically associate with MRFAP1 in MAF1, tumour suppressor RB and the mammalian co repressor complex, mSin3A [14,35]. We observed a 12.9 fold increase of MRFAP1 in CIN III libraries when compared to NC libraries. Notably, the majority of this increase occurs between mild/moderate and severe dysplasia. Although the biological mechanisms of MRFAP1 and MRG15 interaction has yet to be described, it is plausible that the relationship may influence the essential role MRG15 plays in chromatin remodelling.

We observed a 20 fold increase of MORF4L2 in CIN III libraries. Like MRG15, MORF4L2 also belongs to the MRG/MORF transcriptional regulator family of genes involved in senescence and cell growth [37]. Although the two genes share nearly 90% similarity, MORF4L2 is unique to vertebrates whereas MRG15 is evolutionarily conserved [14,37]. Similarly to MRG15, MORF4L2 directly interacts with MRFAP1 and the mSin3A co repressor complex, components of which include HDAC1 and HDCA2 and are involved in chromatin remodelling [13,14]. MORF4L2 is reported to have both repressive and stimulatory activity in transcription regulation of well known cell cycle regulator, v-myb myeloblastosis viral oncogene homolog (avian)-like 2 (MYBL2) [14,38].

MRFAP1 and MORFL2 were selected from Network B for validation in 22 new samples as described above. In both cases, six of seven CIN III samples showed the greatest increase in expression (1.8 – 6.6 fold increase) and the overall increasing trend was found to be statistically significant (p value .033 and .017, respectively). We observed a 2.0–2.5 fold increase in expression in cervical cancer (T1). Together, these results validate the differential expression again supporting a potential role of chromatin remodelling in cervical cancer progression.

This differential expression in the chromatin remodelling genes could result in changes in the organization of the DNA. If these changes are large enough they could be detectable. Previously, we quantified changes in nuclear texture with increasing grade of CIN [39]. These measured phenotype changes, as seen in Figure Figure77 (an updated graph of the data), begin with CIN I and are more prominent in the high grades of CIN (II, III) and cancer [39]. This appears to correlate CIS with the novel observations presented in the present study with respect to the increased expression of chromatin remodelling complex components indicating that these expression changes result in an alteration of the cell/nuclear phenotype.

Figure 7
An updated box plot from Guillaud et al. showing the correlation between a linear discriminant function score (Texture Score) based entirely on only texture phenotype features measured of the nuclei in formalin fixed quantitatively stained sections of ...


Events in expression change involving genes in Network A and Network B (DHFR, MORFL2, MRFAP1, NCOR1, and SMARCC1) occur at or before the stage of moderate dysplasia (CIN II). These events are maintained, although at a reduced frequency, in CIN III suggesting a role in intermediate events prior to severe dysplasia that may be an important stepping stone in disease progression. Genomic instability is characteristic of CIN III lesions and cervical cancer, as such the non-maintenance of such aberrations in later stages may be due to the masking the of initial events by additional changes acquired in severe dysplasia [3,40,41].

Deregulation of chromatin remodelling functions has analogous effects to mutation in DNA repair components in that the repercussions can be genome wide. It is interesting that evidentially two deregulated genes from both networks investigated are associated with chromatin remodelling (SMARCC1, NCOR1, MRFAP1, MORF4L2) and such disruptions are targeted to the critical stage of moderate dysplasia. The increase in SWI/SNF stabilizing molecule SMARCC1 and other novel genes has not been previously illustrated as events in the early stages of dysplasia development thus providing not only novel candidate markers for screening but a biological function for targeting treatment. Together, our results suggest altered expression events in chromatin remodelling complex components and influencing factors occur in precancerous cervical intraepithelial neoplasia. Future investigation on protein DNA interaction will be necessary to further elucidate the precise role of chromatin remodelling in cervical cancer progression.


Sample Selection and Collection

Specimens were collected as previously described [6], immediately prior to the LEEP (Loop electrosurgical excision procedure) using bite biopsy targeting a small portion of the affected epithelium. These specimens were collected with patient consent and University of British Columbia (UBC) British Columbia Cancer Agency (BCCA) Research Ethics Board approval at the Vancouver General Hospital Women's Clinic at the Vancouver Hospital & Health Science Centre (VHHSC). Anonymized cases were assessed by cervical cancer pathologists at BCCA and were selected without prior knowledge of HPV status. Specimens N1to N4 in this study were observed to be non cancerous squamous epithelia whereas C1toC6 were identified as high grade dysplasia or CIN III. Specimens M1to M3 were assessed as CIN I while M4to M6 as CIN II. Detail description of LEEP cone specimens are listed in Additional file 7. All samples were stored in a RNA preservation solution (RNA later, Ambion) within 2 minutes of retrieval and then placed in storage at -80°C.


The biopsies were homogenised in Lysis Binding buffer (100 mM Tris-HCl, pH7.5, 500 mM LiCl, 10 mM EDTA, pH 8.0, 1% LiDS, 5 mM dithiothreitol). Long SAGE libraries were constructed using the L-SAGE kit (Invitrogen, Ontario, Canada) and sequenced to an average depth of 155,000 tags allowing for deep sampling of the transcriptome without an initial restriction to known genes [42-44].


Raw tag count data for each SAGE library was normalized to tags per million (TPM) to facilitate comparison between libraries. Each sequence tag was then mapped to genes using the August 1 2006 version of SAGEGenie [45]. CIN I and CIN II libraries were grouped in analysis (CIN I/CIN II or mild/moderate dysplasia). Subsequently, those tags detected at a level of 20 TPM in at least one of the NC libraries, CINI/CINII, and CINIII groups were retained for differential expression analysis. To determine differential expression between the three groups, a combination of a two-fold difference in means as well as a permutation score ≥ 1.96 (corresponding to an unadjusted p-value of 0.05) was used [46]. Briefly, the permutation test used in this analysis examines the statistical significance of the observed differences in means between two groups in comparison to the average and standard deviation of difference of 10,000 random permutations of the data.

cDNA synthesis

A validation panel of 12 cervical specimens (independent from those used in the SAGE analysis) were assembled: normal (Na, Nb, Nc), CIN I (CIa, CIb, CIc), CIN II (CIIa, CIIb, CIIc) and CIN III (CIIIa, CIIIb, CIIIc). Specimens were collected and stored as described above. RNA was isolated as previously described [6]. In addition, cervix squamous cell carcinoma (T1) and adjacent normal tissue (N1) total RNA sample set was purchased (Ambion, #7276). All cDNA was generated from 200 ng of total RNA using the High Capacity TaqMan Reverse Transcription Reagents (Applied Biosystems, Foster City, CA).

Realtime PCR Analysis

The expression levels of individual genes were analyzed by real-time PCR using TaqMan® Gene Expression Assays in triplicate in a new panel of 22 cases using the ABI 7500 Real-Time PCR System (Applied Biosystems). Taqman probes (Assay IDs) include 18S (Hs99999901_s1), MORF4L2 (Hs0020211_m1), MRFAP1 (Hs00738144_g1), SMARCC1 (Hs00268265_m1), NCOR1 (Hs00196920_m1), DHFR (Hs0075822_s1), PTEN (Hs00829813_s1), RBL2 (Hs00180562_m1) and CDKN2B (Hs00793225_m1). The relative quantification of the target genes in CIN I, CIN II and CIN III samples compared to the average Ct of NC samples was performed using the established 2-ΔΔCt method (Applied Biosystems, Relative Quantitation Of Gene Expression, ABI PRISM 7700 Sequence Detection System: User Bulletin #2). Gene expression is normalized to 18S. The relative quantification values were then plotted, a ratio of one indicating no change with respect to NC cervical tissue. A single squamous cervical carcinoma sample and adjacent normal tissue were also compared using this method (T1 and N1, respectively).


Cervical Intraepithelial Neoplasia (CIN), Human Papillomavirus (HPV), Serial Analysis of Gene Expression (SAGE), Tags per million (TPM)

Authors' contributions

AS designed, organized and performed this study and wrote the manuscript. RC performed data analysis. AP performed gene specific expression assays. KL contributed to the project design, directed library construction and edited the manuscript. DM, TE, DN, JM, RR, MF and MG were responsible for clinical diagnosis, sample acquisition and pathology assessment. WL and CM are shared principle investigators of this study. All authors read and approved the final draft of this manuscript.

Supplementary Material

Additional file 1:

Genes differentially expressed in the early stages of neoplasia. Scales tags differentially expressed between normal and CIN I/II.

Additional file 2:

Genes differentially expressed in the later stages of neoplasia. Scales tags differentially expressed between CIN I/II and CIN III.

Additional file 3:

Tags mapping to the HPV 16 genome. Tags and scaled tag counts described.

Additional file 4:

Genes concurrently expressed in all 16 libraries. Scales tags expressed at similar levels in all 16 libraries.

Additional file 5:

Genes differentially expressed between normal cervical epithelium and severe dysplasia. Scales tags differentially expressed between normal and CIN III.

Additional file 6:

Validation panel quantitative PCR results of genes with altered expression in CIN III L-SAGE libraries. A panel of 22 new cases was investigated in triplicate for expression changes for eight genes. Samples are indicated as follows: CIN I (A-C), CIN I/II (D), CIN II (E-J), CIN III (K-Q) and tumour (R). Zero on the Y-axis denotes mean expression levels of the respective genes in three NC cervical tissues. Expression in tumour is relative to matched normal.

Additional file 7:

Pathology on all cases. Pathology on all cases used as described by pathologist at the British Columbia Cancer Agency.


We are grateful to Anita Carraro for sample coordination and retrieval and Blair Gervan, and Baljit Kamoh for their technical assistance. DNA sequencing was performed at the BC Cancer Agency Michael Smith Genome Sciences Centre. This work was made possible by support from the National Institute of Health (3P0l-CA8270-07 and 7R0l-CA103830-03), the Canadian Institutes for Health Research and the Michael Smith Foundation for Health Research.


  • Scheurer ME, Tortolero-Luna G, Adler-Storthz K. Human papillomavirus infection: biology, epidemiology, and prevention. Int J Gynecol Cancer. 2005;15:727–746. doi: 10.1111/j.1525-1438.2005.00246.x. [PubMed] [Cross Ref]
  • Woodman CB, Collins SI, Young LS. The natural history of cervical HPV infection: unresolved issues. Nat Rev Cancer. 2007;7:11–22. doi: 10.1038/nrc2050. [PubMed] [Cross Ref]
  • Wheeler CM. Advances in primary and secondary interventions for cervical cancer: human papillomavirus prophylactic vaccines and testing. Nat Clin Pract Oncol. 2007;4:224–235. doi: 10.1038/ncponc0770. [PubMed] [Cross Ref]
  • Cutts FT, Franceschi S, Goldie S, Castellsague X, de Sanjose S, Garnett G, Edmunds WJ, Claeys P, Goldenthal KL, Harper DM, Markowitz L. Human papillomavirus and HPV vaccines: a review. Bull World Health Organ. 2007;85:719–726. doi: 10.2471/BLT.06.038414. [PubMed] [Cross Ref]
  • Perez-Plasencia C, Riggins G, Vazquez-Ortiz G, Moreno J, Arreola H, Hidalgo A, Pina-Sanchez P, Salcedo M. Characterization of the global profile of genes expressed in cervical epithelium by Serial Analysis of Gene Expression (SAGE) BMC Genomics. 2005;6:130. doi: 10.1186/1471-2164-6-130. [PMC free article] [PubMed] [Cross Ref]
  • Shadeo A, Chari R, Vatcher G, Campbell J, Lonergan KM, Matisic J, van Niekerk D, Ehlen T, Miller D, Follen M, Lam WL, Macaulay C. Comprehensive serial analysis of gene expression of the cervical transcriptome. BMC Genomics. 2007;8:142. doi: 10.1186/1471-2164-8-142. [PMC free article] [PubMed] [Cross Ref]
  • Gius D, Funk MC, Chuang EY, Feng S, Huettner PC, Nguyen L, Bradbury CM, Mishra M, Gao S, Buttin BM, Cohn DE, Powell MA, Horowitz NS, Whitcomb BP, Rader JS. Profiling microdissected epithelium and stroma to model genomic signatures for cervical carcinogenesis accommodating for covariates. Cancer Res. 2007;67:7113–7123. doi: 10.1158/0008-5472.CAN-07-0260. [PubMed] [Cross Ref]
  • Edgar R, Domrachev M, Lash AE. Gene Expression Omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002;30:207–210. doi: 10.1093/nar/30.1.207. [PMC free article] [PubMed] [Cross Ref]
  • Draghici S, Khatri P, Martins RP, Ostermeier GC, Krawetz SA. Global functional profiling of gene expression. Genomics. 2003;81:98–104. doi: 10.1016/S0888-7543(02)00021-6. [PubMed] [Cross Ref]
  • Draghici S, Khatri P, Bhavsar P, Shah A, Krawetz SA, Tainsky MA. Onto-Tools, the toolkit of the modern biologist: Onto-Express, Onto-Compare, Onto-Design and Onto-Translate. Nucleic Acids Res. 2003;31:3775–3781. doi: 10.1093/nar/gkg624. [PMC free article] [PubMed] [Cross Ref]
  • Ingenuity® Systems, Ingenuity Pathways Analysis
  • Hanahan D, Weinberg RA. The hallmarks of cancer. Cell. 2000;100:57–70. doi: 10.1016/S0092-8674(00)81683-9. [PubMed] [Cross Ref]
  • Yochum GS, Ayer DE. Role for the mortality factors MORF4, MRGX, and MRG15 in transcriptional repression via associations with Pf1, mSin3A, and Transducin-Like Enhancer of Split. Mol Cell Biol. 2002;22:7868–7876. doi: 10.1128/MCB.22.22.7868-7876.2002. [PMC free article] [PubMed] [Cross Ref]
  • Tominaga K, Leung JK, Rookard P, Echigo J, Smith JR, Pereira-Smith OM. MRGX is a novel transcriptional regulator that exhibits activation or repression of the B-myb promoter in a cell type-dependent manner. J Biol Chem. 2003;278:49618–49624. doi: 10.1074/jbc.M309192200. [PubMed] [Cross Ref]
  • Sohn DH, Lee KY, Lee C, Oh J, Chung H, Jeon SH, Seong RH. SRG3 interacts directly with the major components of the SWI/SNF chromatin remodeling complex and protects them from proteasomal degradation. J Biol Chem. 2007;282:10614–10624. doi: 10.1074/jbc.M610563200. [PubMed] [Cross Ref]
  • Shanahan F, Seghezzi W, Parry D, Mahony D, Lees E. Cyclin E associates with BAF155 and BRG1, components of the mammalian SWI-SNF complex, and alters the ability of BRG1 to induce growth arrest. Mol Cell Biol. 1999;19:1460–1469. [PMC free article] [PubMed]
  • Leung JK, Berube N, Venable S, Ahmed S, Timchenko N, Pereira-Smith OM. MRG15 activates the B-myb promoter through formation of a nuclear complex with the retinoblastoma protein and the novel protein PAM14. J Biol Chem. 2001;276:39171–39178. doi: 10.1074/jbc.M103435200. [PubMed] [Cross Ref]
  • Bowman BR, Moure CM, Kirtane BM, Welschhans RL, Tominaga K, Pereira-Smith OM, Quiocho FA. Multipurpose MRG domain involved in cell senescence and proliferation exhibits structural homology to a DNA-interacting domain. Structure. 2006;14:151–158. doi: 10.1016/j.str.2005.08.019. [PubMed] [Cross Ref]
  • Li D, Roberts R. WD-repeat proteins: structure characteristics, biological function, and their involvement in human diseases. Cell Mol Life Sci. 2001;58:2085–2097. doi: 10.1007/PL00000838. [PubMed] [Cross Ref]
  • Tsukada S, Iwai M, Nishiu J, Itoh M, Tomoike H, Horiuchi M, Nakamura Y, Tanaka T. Inhibition of Experimental Intimal Thickening in Mice Lacking a Novel G-Protein-Coupled Receptor. Circulation. 2003;107:313–319. doi: 10.1161/01.CIR.0000043804.29963.B4. [PubMed] [Cross Ref]
  • Giacinti C, Giordano A. RB and cell cycle progression. Oncogene. 2006;25:5220–5227. doi: 10.1038/sj.onc.1209615. [PubMed] [Cross Ref]
  • Claudio PP, Howard CM, Fu Y, Cinti C, Califano L, Micheli P, Mercer EW, Caputi M, Giordano A. Mutations in the retinoblastoma-related gene RB2/p130 in primary nasopharyngeal carcinoma. Cancer Res. 2000;60:8–12. [PubMed]
  • Zamparelli A, Masciullo V, Bovicelli A, Santini D, Ferrandina G, Minimo C, Terzano P, Costa S, Cinti C, Ceccarelli C, Mancuso S, Scambia G, Bovicelli L, Giordano A. Expression of cell-cycle-associated proteins pRB2/p130 and p27kip in vulvar squamous cell carcinomas. Hum Pathol. 2001;32:4–9. doi: 10.1053/hupa.2001.20371. [PubMed] [Cross Ref]
  • Cinti C, Claudio PP, Howard CM, Neri LM, Fu Y, Leoncini L, Tosi GM, Maraldi NM, Giordano A. Genetic alterations disrupting the nuclear localization of the retinoblastoma-related gene RB2/p130 in human tumor cell lines and primary tumors. Cancer Res. 2000;60:383–389. [PubMed]
  • Zhang B, Chen W, Roman A. The E7 proteins of low- and high-risk human papillomaviruses share the ability to target the pRB family member p130 for degradation. Proc Natl Acad Sci U S A. 2006;103:437–442. doi: 10.1073/pnas.0510012103. [PubMed] [Cross Ref]
  • Cristiano S. SWI/SNF: The crossroads where extracellular signaling pathways meet chromatin. Journal of Cellular Physiology. 2006;207:309–314. doi: 10.1002/jcp.20514. [PubMed] [Cross Ref]
  • Cho KS, Elizondo LI, Boerkoel CF. Advances in chromatin remodeling and human disease. Curr Opin Genet Dev. 2004;14:308–315. doi: 10.1016/j.gde.2004.04.015. [PubMed] [Cross Ref]
  • Assaraf YG. Molecular basis of antifolate resistance. Cancer Metastasis Rev. 2007;26:153–181. doi: 10.1007/s10555-007-9049-z. [PubMed] [Cross Ref]
  • Gunawardena RW, Fox SR, Siddiqui H, Knudsen ES. SWI/SNF activity is required for the repression of deoxyribonucleotide triphosphate metabolic enzymes via the recruitment of mSin3B. J Biol Chem. 2007;282:20116–20123. doi: 10.1074/jbc.M701406200. [PubMed] [Cross Ref] [PubMed] [Cross Ref]
  • Underhill C, Qutob MS, Yee SP, Torchia J. A novel nuclear receptor corepressor complex, N-CoR, contains components of the mammalian SWI/SNF complex and the corepressor KAP-1. J Biol Chem. 2000;275:40463–40470. doi: 10.1074/jbc.M007864200. [PubMed] [Cross Ref]
  • Dillon RL, White DE, Muller WJ. The phosphatidyl inositol 3-kinase signaling network: implications for human breast cancer. Oncogene. 26:1338–1345. doi: 10.1038/sj.onc.1210202. [PubMed] [Cross Ref]
  • Kim JW, Namkoong SE, Ryu SW, Kim HS, Shin JW, Lee JM, Kim DH, Kim IK. Absence of p15INK4B and p16INK4A gene alterations in primary cervical carcinoma tissues and cell lines with human papillomavirus infection. Gynecol Oncol. 1998;70:75–79. doi: 10.1006/gyno.1998.5041. [PubMed] [Cross Ref]
  • Lee JS, Choi YD, Lee JH, Nam JH, Choi C, Lee MC, Park CS, Kim HS, Min KW. Expression of PTEN in the progression of cervical neoplasia and its relation to tumor behavior and angiogenesis in invasive squamous cell carcinoma. J Surg Oncol. 2006;93:233–240. doi: 10.1002/jso.20493. [PubMed] [Cross Ref]
  • Brooks LA, Sullivan A, O'Nions J, Bell A, Dunne B, Tidy JA, Evans DJ, Osin P, Vousden KH, Gusterson B, Farrell PJ, Storey A, Gasco M, Sakai T, Crook T. E7 proteins from oncogenic human papillomavirus types transactivate p73: role in cervical intraepithelial neoplasia. Br J Cancer. 2002;86:263–268. doi: 10.1038/sj.bjc.6600033. [PMC free article] [PubMed] [Cross Ref]
  • Zhang P, Zhao J, Wang B, Du J, Lu Y, Chen J, Ding J. The MRG domain of human MRG15 uses a shallow hydrophobic pocket to interact with the N-terminal region of PAM14. Protein Sci. 2006;15:2423–2434. doi: 10.1110/ps.062397806. [PubMed] [Cross Ref]
  • Tominaga K, Kirtane B, Jackson JG, Ikeno Y, Ikeda T, Hawks C, Smith JR, Matzuk MM, Pereira-Smith OM. MRG15 regulates embryonic development and cell proliferation. Mol Cell Biol. 2005;25:2924–2937. doi: 10.1128/MCB.25.8.2924-2937.2005. [PMC free article] [PubMed] [Cross Ref]
  • Tominaga K, Matzuk MM, Pereira-Smith OM. MrgX is not essential for cell growth and development in the mouse. Mol Cell Biol. 2005;25:4873–4880. doi: 10.1128/MCB.25.12.4873-4880.2005. [PMC free article] [PubMed] [Cross Ref]
  • Joaquin M, Watson RJ. Cell cycle regulation by the B-Myb transcription factor. Cell Mol Life Sci. 2003;60:2389–2401. doi: 10.1007/s00018-003-3037-4. [PubMed] [Cross Ref]
  • Guillaud M, Cox D, Adler-Storthz K, Malpica A, Staerkel G, Matisic J, Van Niekerk D, Poulin N, Follen M, MacAulay C. Exploratory analysis of quantitative histopathology of cervical intraepithelial neoplasia: objectivity, reproducibility, malignancy-associated changes, and human papillomavirus. Cytometry A. 2004;60:81–89. doi: 10.1002/cyto.a.20034. [PubMed] [Cross Ref]
  • Lockwood WW, Coe BP, Williams AC, MacAulay C, Lam WL. Whole genome tiling path array CGH analysis of segmental copy number alterations in cervical cancer cell lines. Int J Cancer. 2007;120:436–443. doi: 10.1002/ijc.22335. [PubMed] [Cross Ref]
  • Narayan G, Bourdon V, Chaganti S, Arias-Pulido H, Nandula SV, Rao PH, Gissmann L, Durst M, Schneider A, Pothuri B, Mansukhani M, Basso K, Chaganti RS, Murty VV. Gene dosage alterations revealed by cDNA microarray analysis in cervical cancer: identification of candidate amplified and overexpressed genes. Genes Chromosomes Cancer. 2007;46:373–384. doi: 10.1002/gcc.20418. [PubMed] [Cross Ref]
  • Velculescu VE, Zhang L, Vogelstein B, Kinzler KW. Serial analysis of gene expression. Science. 1995;270:484–487. doi: 10.1126/science.270.5235.484. [PubMed] [Cross Ref]
  • Li YJ, Xu P, Qin X, Schmechel DE, Hulette CM, Haines JL, Pericak-Vance MA, Gilbert JR. A comparative analysis of the information content in long and short SAGE libraries. BMC Bioinformatics. 2006;7:504. doi: 10.1186/1471-2105-7-504. [PMC free article] [PubMed] [Cross Ref]
  • Chen J, Sun M, Lee S, Zhou G, Rowley JD, Wang SM. Identifying novel transcripts and novel genes in the human genome by using novel SAGE tags. Proc Natl Acad Sci U S A. 2002;99:12257–12262. doi: 10.1073/pnas.192436499. [PubMed] [Cross Ref]
  • Boon K, Osorio EC, Greenhut SF, Schaefer CF, Shoemaker J, Polyak K, Morin PJ, Buetow KH, Strausberg RL, De Souza SJ, Riggins GJ. An anatomy of normal and malignant gene expression. Proc Natl Acad Sci U S A. 2002;99:11287–11292. doi: 10.1073/pnas.152324199. [PubMed] [Cross Ref]
  • Good P. Permutation Tests: A Practical Guide to Resampling Methods for Testing Hypotheses. 2nd ed. Springer-Verlag New York, Inc. ; 2000.

Articles from BMC Genomics are provided here courtesy of BioMed Central