Next Article in Journal
Alteration to Dopaminergic Synapses Following Exposure to Perfluorooctane Sulfonate (PFOS), in Vitro and in Vivo
Next Article in Special Issue
Combination of mTOR and MAPK Inhibitors—A Potential Way to Treat Renal Cell Carcinoma
Previous Article in Journal
Relationship between Neutrophil-to-Lymphocyte Ratio and Inflammatory Markers in Sickle Cell Anaemia Patients with Proteinuria
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Medical Sciences
Volume 4
Issue 3
10.3390/medsci4030012
medsci-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Novel Somatic Copy Number Alteration Identified for Cervical Cancer in the Mexican American Population

by
Alireza Torabi
1,
Javier Ordonez
2,
Brenda Bin Su
4,
Laura Palmer
3,
Chunxiang Mao
3,
Katherine E. Lara
3,
Lewis P. Rubin
3 and
Chun Xu
3,*
1
Department of Pathology, TTUHSC, El Paso 79905, TX, USA
2
Department of Biomedical Science, TTUHSC, El Paso 79905, TX, USA
3
Department of Pediatrics, Texas Tech University Health Sciences Center (TTUHSC), El Paso 79905, TX, USA
4
Department of Internal Medicine, College of Medicine and Health Sciences, UAE University, Al-Ain 15551, UAE
*
Author to whom correspondence should be addressed.
Med. Sci. 2016, 4(3), 12; https://doi.org/10.3390/medsci4030012
Submission received: 1 June 2016 / Revised: 15 July 2016 / Accepted: 25 July 2016 / Published: 3 August 2016
(This article belongs to the Special Issue Gene Therapy and Cancer: Current Developments)
Browse Figures
Review Reports Versions Notes

Abstract

:
Cervical cancer affects millions of Americans, but the rate for cervical cancer in the Mexican American is approximately twice that for non-Mexican Americans. The etiologies of cervical cancer are still not fully understood. A number of somatic mutations, including several copy number alterations (CNAs), have been identified in the pathogenesis of cervical carcinomas in non-Mexican Americans. Thus, the purpose of this study was to investigate CNAs in association with cervical cancer in the Mexican American population. We conducted a pilot study of genome-wide CNA analysis using 2.5 million markers in four diagnostic groups: reference (n = 125), low grade dysplasia (cervical intraepithelial neoplasia (CIN)-I, n = 4), high grade dysplasia (CIN-II and -III, n = 5) and invasive carcinoma (squamous cell carcinoma (SCC), n = 5) followed by data analyses using Partek. We observed a statistically-significant difference of CNA burden between case and reference groups of different sizes (>100 kb, 10–100 kb and 1–10 kb) of CNAs that included deletions and amplifications, e.g., a statistically-significant difference of >100 kb deletions was observed between the reference (6.6%) and pre-cancer and cancer (91.3%) groups. Recurrent aberrations of 98 CNA regions were also identified in cases only. However, none of the CNAs have an impact on cancer progression. A total of 32 CNA regions identified contained tumor suppressor genes and oncogenes. Moreover, the pathway analysis revealed endometrial cancer and estrogen signaling pathways associated with this cancer (p < 0.05) using Kyoto Encyclopedia of Genes and Genomes (KEGG). This is the first report of CNAs identified for cervical cancer in the U.S. Latino population using high density markers. We are aware of the small sample size in the study. Thus, additional studies with a larger sample are needed to confirm the current findings.

1. Introduction

Based on the U.S. Centers for Disease Control and Prevention (CDC) 2009 report, 12,357 women were diagnosed with cervical cancer, and 3909 women died from cervical cancer in the United States. This huge number of patients with cervical cancer has been a critical issue and is responsible for 10%–15% of cancer-related deaths in females globally [1]. Cervical cancer is the second most common malignant tumor in women worldwide. Among these tumors, approximately 80% are squamous cell carcinomas (SCCs), and 5%–20% are adenocarcinomas (AdCAs) [2,3]. Cervical SCCs are developed from a premalignant disease known as cervical intraepithelial neoplasia (CIN) Graded 1–3 with increasing atypical features. The five-year overall survival for cervical cancer is only 66%.
In addition, increasing evidence suggests that infection with high-risk subtypes of human papillomavirus (HPV) (e.g., HPV-16 and HPV-18) is the most common cause and is the primary initiator of premalignant lesions [4]. However, only a small proportion of women infected with oncogenic HPV subtypes develop cervical cancer, which suggests that HPV infection alone is insufficient to cause cancer, and there is a possibility of other host factors linked to the development of invasive cervical cancer [5], like genetic variation, including polymorphisms, insertions or deletions in the host genome [6,7]. Increasing evidence demonstrates that there was a consistent relationship between certain genetic variants (such as the tumor protein 53 (TP53) Arg72Pro polymorphism) and cervical cancer, most likely modulated by the presence of high-risk HPV during progression from squamous intraepithelial lesions (SIL) to cervical cancer.
Moreover, Mexican American women in Texas have among the highest rates of cervical cancer incidence and mortality in the country [8] with twice the frequency as compared to their non-Mexican American counterparts. The annual death rate from cervical cancer for Mexican Americans is 24.2 out of 100,000 [9]. However, there is a lack of studies on somatic mutation identification for cervical cancer in this population. Only four studies have reported on genetic basis of cervical cancer in Mexican American women based on a PubMed search (19 November 2015). The number of the Mexican American population was estimated at 50.5 million in the 2010 census, making Mexican American the largest minority group in the U.S., as well as a rapidly growing segment of the U.S. population.
Recent advances in genome studies have led to the discovery of one important type of variation that can be assessed with recent technology: copy number alterations (CNAs), usually for a cancer study, somatic copy number changes and/or copy number variations (CNVs), usually for a non-cancer study, or germline copy number changes, such as CNVs identified for neuropsychiatric disorders in our recent study [10]. These CNAs or CNVs are by definition chromosomal regions with sizes of 1 kb to several Mb, which vary across individuals with regard to the number of copies of a chromosomal segment. CNVs refer to structural variations of the DNA that include insertions, deletions and duplications. Studies have found that CNVs cover as much as 14% of the human genome [11], and there is a much higher de novo rate as many as 10–1000-fold in CNVs as compared to single nucleotide polymorphisms (SNPs) [11,12]. Furthermore, CNVs have been shown to account for more genomic differences between individuals than SNPs [13,14]. Therefore, CNVs may contribute a sizeable amount of disease phenotypic variation in each individual of a population [15].
Analyses integrating mutation information with data on rearrangements and CNAs have revealed a higher-order organization of the seemingly random genetic events that lead to cancer [16]. Interestingly, genes in regions subject to copy number changes appear to be organized along functional ontological terms related to cancer [16]. Studies have also implicated a number of somatic mutations, including TP53, phosphatidylinositol-4,5-bisphosphate 3-kinase catalytic subunit alpha (PIK3CA), phosphatase and tensin homolog (PTEN), serine/threonine kinase 11 (STK11) and V-Ki-ras2 Kirsten rat sarcoma viral oncogene homolog (KRAS) [17,18], and several CNAs in the pathogenesis of cervical carcinomas [19,20] in non-Mexican Americans.
In this study, we carried out a genome-wide survey of potential somatic CNAs, including amplifications and deletions, in apparently normal tissues (n = 2), low grade dysplasia (n = 4), high grade dysplasia (n = 5), invasive carcinoma (n = 5) and blood samples (n = 125, serving as a reference group) of female subjects from the HapMap data. We are aware of the limited number of cases and the lack of a control group. Thus, in the future, a large study with a control sample and more cases as a methodological alternative is needed. We genotyped 2.5 million markers and analyzed somatic CNAs in a total of 14 tissues using the Illumina HumanOmni2.5-8 BeadChip Kit at tissue-level resolution. We mapped genomic changes: (1) between peripheral blood samples of the reference subjects and cervical tissues from the cases (cervical dysplasia and invasive carcinoma); because only two normal cervical tissue samples had insufficient statistical power, we excluded these two samples from the further analysis; (2) we also analyzed genomic changes among four diagnostic groups (normal, low, high grade and invasive carcinoma). We expect that this study: (i) will provide an estimate of the prevalence of somatic CNAs by identifying specific patterns, genes and/or biological pathways associated with different stages of cervical dysplasia in Mexican Americans; (ii) will investigate the genomic context of these somatic CNAs; and (iii) will evaluate whether the burden of somatic mutations predicts tumor progression.

2. Materials and Methods

2.1. Materials

A total of 14 tissues (low, high grade dysplasia and invasive carcinoma) from cases and 125 female subjects, serving as a reference group, were used for this study. The demographic information is shown in Table 1. The cases were categorized into three groups, including low grade dysplasia (CIN-I, n = 4), high grade dysplasia (CIN-II and III, n = 5) and squamous cell carcinoma (SCC, n = 5) groups. For a broader definition of the case group, we also divided the cases into two groups, pre-cancer (CIN-I, -II and -III) and cancer (SCC). All of the case subjects in the current study were from the Mexican American population recruited from the outpatient clinics at the University Medical Center (UMC) and Texas Tech University Health Sciences Center (TTUHSC)-El Paso. All cases of cervical cancer were diagnosed as SCC by histopathological examinations, whereas healthy women had no abnormal cytological findings in the Pap smear tests of the uterine cervix.
The Illumina HumanOmni2.5-8 BeadChip data from a total of 125 blood samples of female subjects from the HapMap were used as a reference in the current study.

2.2. Methods

2.2.1. Tissue Specimens

Sixteen cervical tissue samples were obtained from the Department of Pathology, TTUHSC-El Paso. We excluded two normal tissues due to an insufficient tissue sample size. All case subjects were HPV positive, from the Mexican American population and had signed Institutional Review Board-approved written informed consent forms prior to enrolling in the study. The procedures were approved by the Institutional Review Boards of TTUHSC (IRB #E13107), and the study was performed in accordance with the Helsinki Declaration of 1975.

2.2.2. Microdissection

Paraffin-embedded tissues were first sectioned into 10-μm slices, which were hematoxylin-eosin stained for the selection of the appropriate tissue area. The corresponding selected areas of each tissue sample were then collected under microscopic observation with a 30-gauge needle (Becton-Dickinson, Franklin Lakes, NJ, USA).
Genomic DNA of micro-dissected tissue was extracted by proteinase K digestion followed by standard phenol-chloroform extraction. The QIAamp Formalin-Fixed, Paraffin-Embedded (FFPE) Tissue Kit from Qiagen (Valencia, CA, USA), which is widely used for extracting DNA from FFPE sections, was used. The experiment was performed according to the manufacturer’s handbook. The total amount of DNA was spectrophotometrically determined by measuring the absorbance at 260 nm (A260), and DNA purity was assessed by detecting the A260/A280 ratio using the Varioskan Flash (Thermo Scientific, Rockford, IL, USA) according to the manufacturer’s instructions.

2.2.3. High Density Genotyping

Genomic DNA from the case group was used to obtain genotypes by the Illumina HumanOmni2.5-8 BeadChip (Illumina, San Diego, CA, USA). This DNA chip provided over 2.5 million markers at a median spacing of 1.2 kb and full support of CNV or CNA applications, which was a powerful genotyping tool and allowed us to make more meaningful discoveries. Genomic annotations were based on National Center for Biotechnology Information (NCBI) Human Genome Build 37 (University of California, Santa Cruz Genome Browser Release 19), and genotyping experiments were performed at the Genomics Core at the TTUHSC-El Paso.
Due to our limited number of control samples of cervical tissues, we used publicly available HapMap data with genotypes, including CNV data of the same DNA chip, the Illumina HumanOmni2.5-8 BeadChip Kit, in the 125 female subjects as a reference group. The HapMap data with genotypes were downloaded from [21].
In addition, there is no genotype data available for the Mexican American population using the HapMap data, so we selected these 125 female subjects from the Admix populations (Utah residents with Northern and Western European ancestry, CEU; Yoruba, YRI; Han Chinese, CHB, and Japanese, JPT) since the various combinations of reference panels with the multi-ethnic populations had better accuracy than those containing only single ethnic samples [22] for genetic study and genetic imputation analysis.
The raw genotyping signal data were processed by the Illumina GenomeStudio software Genotyping Module Version 3.2.33 (Illumina) and converted to allele-specific intensity values. The genotype call rate threshold was set at ≥95% for all samples, and a total of 139 including case (n = 14) and reference (n = 125) subjects passed the quality control. A Partek customized report of the normalized genotype data, composed of those 139 samples, including low grade (CIN-I = 4), high grade (CIN-II and CIN-III = 5), SCC (n = 5) and 125 females from the reference group (Table 1), was transferred to Partek® Genomic Suite® software, Version 6.6, Copyright 2014, Partek Inc. (Saint Louis, MO, USA), for downstream analysis.

2.2.4. CNA Detection

Unpaired copy number analysis was performed in the Partek Genomics Suite comparing allele intensities to a reference baseline of 125 female HapMap admixed samples using a similar analysis strategy as a previous publication [23]. The genomic segmentation algorithm was applied to find break points and to detect amplifications (gains) and deletions (losses). The following stringent parameters were used to identify CNAs and CNA regions as a previous study [24]: (1) each segment must contain a minimum of 10 consecutive filtered probe sets; (2) a p value threshold of 0.001 when compared to the neighboring adjacent regions; and (3) a signal-to-noise threshold of 0.5 and diploid copy number 1.7 to 2.3.

2.2.5. Pathway Analysis

We further examined whether these CNAs in the various genes have an impact on gene functions. Gene ontology (GO) analysis was also performed using the Partek® Genomic Suite® software, Version 6.6, Copyright® 2014, Partek Inc., to investigate whether there was enrichment for cancer-associated CNAs (p < 0.05) in genes from any ontology categories.
Disease association of individual CNA frequencies with patients, as a group, compared to a reference group, was assessed using 2 × 2 or 2 × 3 contingency tables, two-tailed χ2 tests or Fisher exact tests. Statistical analyses were performed using the SPSS statistical software package (Version 10, IBM, Chicago, IL, USA). Differences with two-tailed probability values of p ≤ 0.05 were taken as statistically significant. p-values for tests of CNA association were conservatively corrected for multiple testing using the Bonferroni method.

3. Results

3.1. Overall CNA Patterns

We processed genotype and intensity data for all 2.5 million probes of the Illumina HumanOmni2.5-8 BeadChip Kit for the 14 cases and 125 reference subjects. A number of CNAs (deletions or losses) was present in many chromosomal regions (Figure 1). A total of 2220 CNAs >1000 bp was identified, including 725 (32.7%) amplifications (gains) and 1495 (67.3%) deletions (losses) mainly on the 22 autosomal chromosomes.
The chromosomal locations of the copy number amplifications and deletions of the 22 autosomes and X chromosomes were shown by karyograms (Figure 1). Most of the amplifications are found in the cases at the short arms of chromosomes 1, 5, 8, 16, 19 and 20, as well as the centromeres of chromosomes 14, 15 and 21. Most of the deletions in the cases were mainly observed at the short arms of chromosomes 1, 4, 5, 7, 11, 12, 16, 17 and 19, as well as the centromeres of chromosome 19.

3.2. Principal Component Analysis (PCA)

To characterize aberration profiles in different diagnostic groups (low grade dysplasia, high grade dysplasia, SCC and reference), we performed a PCA (Figure 2). There was a great variability in the three diagnostic groups (reference, pre-cancer and cancer) and four diagnostic groups (reference, low grade, high grade and SCC). These results indicate that there is a distinct difference in the patterns of the CNAs between case and reference groups based on the genomic profiles of these subjects.
Overall, there was a clear separation between the reference and the case group when examined by PCA clustering.

3.3. Genome-Wide CNA and Cancer Progression

Increasing evidence demonstrates that more numerical genomic alterations are associated with progression from precursor lesions to invasive cancer [25]. Thus, we conducted CNA genomic analysis in different diagnostic groups to understand the complexities of the genomic architecture of these highly heterogeneous groups of cervical cancer and to assess whether CNA burden impacts the prognosis of cervical cancer in this Mexican American population. We analyzed CNAs among two (pre-cancer and cancer; Table 2B) and three diagnostic groups (low grade, high grade and SCC; Table 2A) in cases only. There was no statistically-significant difference of CNA burden among diagnostic groups in cases in any size of CNA, although there are slightly higher frequencies of amplifications in the 100 kb, 10–100 kb and 1–10 kb categories in the cancer group (11.0%, 23.3% and 13.4%, respectively) as compared to the pre-cancer group (7.5%, 18.5% and 12.3%, respectively) (Table 2) using 205 two-tailed χ2 tests with the SPSS statistical software package.

3.4. Recurrent Aberrations Identified in Cases Only

We are interested in recurrent CNAs since the chromosome regions with recurrent CNAs more likely to harbor disease-critical genes are those that show alterations that are recurrent among individuals with cancer or other diseases [26,27]. In this context, we can define a recurrent CNA region as a set of contiguous genes (a region) that shows a high enough probability (or evidence) of being altered (e.g., deletion or insertion) in two or more samples, as previously described [28]. The recurrent aberrations in a total of 98 CNA regions containing 849 genes/loci were identified only in the case group (pre-cancer and cancer subjects) with at least two cases who carried these CNAs. None of these CNAs were observed in the reference group (Table 3). Among these 98 recurrent CNA regions, there were 545 deletions (84%) and 104 amplifications (16%) observed in cases only using statistical analyses described in the Methods section.
Table 3 shows the top 20 chromosome regions with a high number of recurrent somatic CNAs (>100 kb) in the cases with at least nine subjects carrying these deletions in each CNA region; none of these CNAs were observed in the reference group. Some of the CNAs contain tumor suppressor genes (e.g., axis inhibition protein 1 (AXIN1) and tuberous sclerosis complex 2 (TSC2)).
Eight exceptionally large CNA regions with deletions of between 1591 kb and 517 kb were detected in cases on short arms of chromosomes 1, 16 and 19.
In addition to recurrent deletions identified in the cases (Table 3), two recurrent CNA regions with a 220-kb amplification on 21q11.2 and a 113.5-kb deletion on 7p22.2 were observed in only CIN-III and SCC.
Next, we checked the Cervical Cancer Database (CCDB) [29] and found three genes, haemoglobin alpha 2 (HBA2), mesothelin (MSLN) and STK11, overlapped between 849 genes/loci observed from our current study and 538 genes listed in the CCDB.

3.5. CNA Regions Contained Known Tumor Suppressor Genes and Oncogenes in Cancer Tissues but Not in the Reference Group

We were also interested in identifying CNA regions that were observed only in cases and over 100 kb containing known cancer-related genes as compared to the references (Table 3, Table 4 and Table 5). Thus, we examined whether newly-identified CNAs contained tumor suppressor genes and/or oncogenes. We observed a total of 26 CNA regions containing known tumor suppressor gene regions. These CNAs were deletions and only seen in the case group. Furthermore, none of these regions were observed in the reference group (Table 4). We also observed a total of six CNA regions with deletions (no amplification) containing known oncogenes in the case group, and none of these regions were observed in the reference group (Table 5).

3.6. Pathway Analysis

Recent studies have incorporated protein networks into the results of genome-wide CNA data using networks or GO analysis to discover disease-associated and/or enriched pathways [30,31,32,33]. Therefore, we conducted GO analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis.
The results of the functional enrichment of KEGG pathway analysis in the identified CNAs >100 kb that occurred in two or more cases are shown in Table 6. Fourteen pathways were discovered, including the insulin signaling pathway and the endometrial cancer and estrogen signaling pathway, with enrichment p values lower than 0.05. A GO analysis was performed by selecting CNAs in the genes to see if any GO categories were overrepresented among CNAs identified in cases [34]. Many biological processes were undisturbed at the molecular level, while others were frequently affected across multiple cases, as shown in Table 7, which lists the most commonly-affected function with p-values lower than 0.05 and more than four genes in each function category.

4. Discussion

It has been demonstrated that somatic structural alterations (e.g., amplifications or deletions) of human chromosomes represent a common class of mutations, which may cause gene disruption (e.g., deletion or rearrangement), gene activation (e.g., CNAs, gain or amplification) or the formation of novel oncogenic gene products (gene fusions). Many of these events actively drive carcinogenesis [35,36]. In our initial cervical cancer cohort, some CNA patterns identified were previously found to be correlated with cervical cancer.
We screened whole genomes with 2.5 million markers of an array to discover any recurrent copy number alterations in cases (pre-cancer and SCC). We identified a total of 98 CNA regions >100 kb in the case group, including low grade dysplasia, high grade dysplasia and SCC. These CNAs occurred in two or more cases, including two large deletions (1591 kb and 1123 kb) (Table 3). Of the top 20 CNA regions >100 kb deletions in cases only, six CNAs occurred on 16p13 that have been reported in cervical cancer [37], and mutation analysis of the AXIN1 gene located at 16p13 was reported to be involved in the Wnt pathway in cervical carcinomas [38]. The CNA containing the SET and MYND domain-containing protein 4 (SMYD4) gene on 17p13.3 was also reported previously in this cancer [37]. The SMYD4 gene is demonstrated to be a potential tumor suppressor and plays a critical role in breast carcinogenesis, at least partly through inhibiting the expression of platelet-derived growth factor receptor-alpha (PDGFRA), and this gene could be a novel target for improving the treatment of breast cancer [39].
A number of previous studies demonstrate that numerical chromosomal aberrations were found to progress to invasive cancer from precursor lesions [25,37]; however, we did not observe this feature in our current study sample. This might be due to the small sample size or any other unknown or not yet identified factors, such as ethnicity, since most previous study populations are non-Hispanic populations.
Among two amplifications identified in CIN-III and SCC, one recurrent CNA region with a 220-kb amplification on 21q11.2 was observed in CIN-III and co-occurred with previous findings, where the amplification had been identified in breast cancer subjects with tamoxifen resistance [40].
The recurrent CNA regions we identified contain a number of tumor suppressor genes, oncogenes and cancer-related genes. A CNA that contains the interferon-induced transmembrane protein 1 (IFITM1) gene at 11p15.5 was identified in eight cases. This gene has been reported to be involved in cervical carcinogenesis [41]. A total of nine cases with the CNA containing the fucosyltransferase 3 (FUT3) gene at 19p13.3 were observed in the current study, and this gene is associated with breast cancer [42]. A CNA region that includes the naked cuticle 2 (NKD2) gene at 5p15.35 was detected in six cases. This gene was found to suppress breast cancer proliferation by inhibiting Wnt signaling [43]. The ovarian cancer-associated gene 2 protein (OVCA2) at 17p13.3 was repeatedly reported to be associated with ovarian cancer, and we observed five cases carrying the CNA at this gene region. A cluster of three tumor suppressor genes, cadherin 1 (CDH1), death-associated protein kinase (DAPK) and hypermethylated in cancer 1 (HIC1), displayed a significantly increased frequency of promoter methylation with progressively more severe cervical neoplasia. In addition, the Hes family BHLH Transcription Factor 5 (HES5) gene was reported to be associated with cervical carcinoma cells using immunocytochemistry, Western blot and methyl thiazolyl tetrazolium assays [44], and the CNAs on 1p36.33-1p36.32 containing the HES5 gene were also identified in the current study. Our newly-identified CNA containing this gene on 7p22.2 was also observed in the recurrent pre-cancer and cancer subjects (Table 3).
Moreover, we also examined deletion burden observed to be similar to those for other cancers and found more deletions identified than amplification in most of the cancer studies, which supports our findings. A recent study using the TCGA data identified nine regions of deletion that were unique to estrogen receptor positive (ER+) post menopause tumors in patients with breast cancer [45], including deletion in 7p22.3, where our newly-identified deletion in cases only was located, and it contains a known tumor suppressor gene.
To analyze the possible effect of genomic alterations, to further capture cancer-causing gene information and to see whether any GO categories are overrepresented among CNA regions, we searched the KEGG pathway database and GO categories and identified a number of pathways and functions where CNAs occur in the SCC group, but not in the reference group (Table 6 and Table 7). Furthermore, our results were partly in agreement with previous reports about cervical cancer. For example, using the functional enrichment of KEGG pathway analysis, we discovered the insulin signaling pathway in the case group, and it was evidenced that HPV 16 E6 oncoprotein interferes with the insulin signaling pathway by binding to tuberin [46]. As is already known, high-risk HPV infection is a causal agent for cervical cancer. In addition, two pathways of interest are endometrial cancer and the estrogen signaling pathway observed in the current study. The NF-κB signaling pathway was also identified in our cancer group (SCC), but not in the reference group, which supports a previous study where bisphenol A (BPA) stimulated the cervical cancer migration via IKK-β/NF-κB signals [47]. The estrogen signaling pathway was also shown in our KEGG analysis. A previous study was performed to evaluate the potential of miRNAs as novel markers for the post-therapeutic monitoring of cervical SCC patients. A regulatory network of differentially-expressed serum miRNAs was identified, and a number of target genes was predicted in the estrogen-mediated signal pathways [48]. Ours and others’ findings support that cervical cancer is a hormone-associated gynecological disease. For example, HPV infection has been associated with the deregulation of the PI3K-Akt-mTOR pathway in invasive cervical carcinomas [49]. The PI3K-Akt signaling pathway was also found in our current study, and 14 genes in this pathway were found in two or more SCC cases (Table 6).
Moreover, signal transducer activity with an enrichment score of 4.2 and a p value of 0.01 containing 11 genes was also identified in a previous report on cervical cancer [50] using GO analysis. Signal transduction was identified in our current study in cervical cancer (SCC), but not in the reference group; a recent study also supports that a signal transduction cascade and mitogen-activated protein kinase kinase kinase 3 (MEKK3) serve as key integration points and are important factors in regulating cellular responses to environmental stress. This signal transduction and MEKK3 not only link diverse extracellular stimuli to subsequent signaling molecules, but they also amplify the initiating signals to ultimately activate effector molecules and induce cell proliferation, differentiation and survival of cervical cancer [51].
These pathways and functions were identified as overlapping biological themes, and these data may provide useful information on the molecular mechanisms for cervical cancer, its prognosis and treatment responses.
There are a number of novelties in this study: (1) this is the first report of CNAs in the relatively ethnically homogeneous group of Mexican Americans using high density mapping (2.5 million in number); previous studies using less density markers may result in an underestimation of the genetic changes that take place in cervical cancer; (2) newly-identified CNAs in cases together with results in the in silico analysis using KEGG and GO function provide insight into how multiple CNAs might contribute to cervical cancer development. We also are aware of some limitations in our study. (1) Our small sample size (14 cases and 125 references) is a major limitation for this type of study. Due to a small number of samples, the nature of genome-wide alteration of copy number may not be fully explained in pre-cancer and SCC. Future confirmation studies using an independent sample will provide an opportunity to more accurately dissect the genetic complexity of somatic CNAs for cervical cancer. (2) There might be a bias of the CNA identified in cervical cancer for the Mexican Americans, since we used the 1000 Genome admixed populations, not Mexican Americans; thus, we currently are recruiting more subjects with cervical cancer from the same population and plan to validate the findings in more samples. (3) Other biomarkers, such as RNA-Seq or DNA methylation profile and sequencing data using next generation sequencing technologies (such as target gene sequencing in these CNA regions), will provide an opportunity for in-depth molecular profiling of the fundamental biological processes of cervical cancer.
We expect that with future validation and confirmation, newly-discovered CNAs can be used in cancer classification, diagnosis, prognosis and treatment responses.

5. Conclusions

Using high density SNP array analysis, we have shown extensive genome-wide CNA changes in pre-cancer and cancer groups as compared to the genome CNA profile in the reference group. Our results demonstrated that the recurrent alterations of CNAs occurred in cases of pre-cancer and SCC. Some of the somatic genomic gains and losses in cervical pre-cancer and cancer were in concordance with the results from previous studies. To our knowledge, no previous studies have applied genome-wide copy number analysis using such high density markers for cervical cancers in the Mexican American population; however, validation and confirmation in a larger sample size are needed.

Acknowledgments

We are grateful to all of the subjects who participated in the study and the many dedicated oncologists at TTUHSC and local hospitals in El Paso for help in patient ascertainment. Thanks are given to Arun Joseph for assisting with part of the data analysis. This study was supported, in part, by the TTUHSC Seed grant (PI, Chun Xu) and the TTUHSC Mini grant (PI, Chun Xu).

Author Contributions

Alireza Torabi provided tissue samples and conceived and designed the experiments. Javier Ordonez performed the experiments and offered data analysis (CNV analysis and bioinformatics analysis). Brenda Bin Su was involved in bioinformatics analysis and wrote the paper. Laura Palmer made a contribution to the data analysis. Chuanxiong Mao made a contribution to the data analysis and manuscript revision. Katherine E. Lara made a contribution to the data mining and data analysis. Lewis P. Rubin made a contribution to part of the manuscript and offered critical comments and suggestions on this paper. Chun Xu, as PI, designed the study, supervised the others and wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jemal, A.; Bray, F.; Center, M.M.; Ferlay, J.; Ward, E.; Forman, D. Global cancer statistics. CA Cancer J. Clin. 2011, 61, 69–90. [Google Scholar] [CrossRef] [PubMed]
  2. Green, J.; Berrington de Gonzalez, A.; Sweetland, S.; Beral, V.; Chilvers, C.; Crossley, B.; Deacon, J.; Hermon, C.; Jha, P.; Mant, D.; et al. Risk factors for adenocarcinoma and squamous cell carcinoma of the cervix in women aged 20–44 years: The UK national case-control study of cervical cancer. Br. J. Cancer 2003, 89, 2078–2086. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Pisani, P.; Bray, F.; Parkin, D.M. Estimates of the world-wide prevalence of cancer for 25 sites in the adult population. Int. J. Cancer 2002, 97, 72–81. [Google Scholar] [CrossRef] [PubMed]
  4. Termini, L.; Maciag, P.C.; Soares, F.A.; Nonogaki, S.; Pereira, S.M.; Alves, V.A.; Longatto-Filho, A.; Villa, L.L. Analysis of human kallikrein 7 expression as a potential biomarker in cervical neoplasia. Int. J. Cancer 2010, 127, 485–490. [Google Scholar] [CrossRef] [PubMed]
  5. Lockwood, W.W.; Coe, B.P.; Williams, A.C.; MacAulay, C.; Lam, W.L. Whole genome tiling path array CGH analysis of segmental copy number alterations in cervical cancer cell lines. Int. J. Cancer 2007, 120, 436–443. [Google Scholar] [CrossRef] [PubMed]
  6. Yu, K.J.; Rader, J.S.; Borecki, I.; Zhang, Z.; Hildesheim, A. CD83 polymorphisms and cervical cancer risk. Gynecol. Oncol. 2009, 114, 319–322. [Google Scholar] [CrossRef] [PubMed]
  7. Ma, D.; Hovey, R.L.; Zhang, Z.; Fye, S.; Huettner, P.C.; Borecki, I.B.; Rader, J.S. Genetic variations in EGFR and ERBB4 increase susceptibility to cervical cancer. Gynecol. Oncol. 2013, 131, 445–450. [Google Scholar] [CrossRef] [PubMed]
  8. Cancer Facts and Figures for Hispanics/Latinos 2009–2011; American Cancer Society: Atlanta, GA, USA, 2009; Available online: http://www.cancer.org/acs/groups/content/@nho/documents/document/ffhispanicslatinos20092011.pdf (accessed on 24 July 2016).
  9. McDougall, J.A.; Madeleine, M.M.; Daling, J.R.; Li, C.I. Racial and ethnic disparities in cervical cancer incidence rates in the United States, 1992–2003. Cancer Causes Control 2007, 18, 1175–1186. [Google Scholar] [CrossRef] [PubMed]
  10. Xu, C.; Blackburn, A.N.; Gonzalez, S.; Villa, E.C.; Ramirez, M.; Zavala, J.; Rodriguez, M.; Camarillo, C.; Ordonez, J.; Armas, R.; et al. Rare Copy Number Variants in Schizophrenia and Bipolar Disorder in a Latino Population; American Society of Human Genetics: San Francisco, CA, USA, 2012. [Google Scholar]
  11. Zhang, F.; Gu, W.; Hurles, M.E.; Lupski, J.R. Copy number variation in human health, disease, and evolution. Annu. Rev. Genom. Hum. Genet. 2009, 10, 451–481. [Google Scholar] [CrossRef] [PubMed]
  12. Lupski, J.R. Structural variation in the human genome. New Engl. J. Med. 2007, 356, 1169–1171. [Google Scholar] [CrossRef] [PubMed]
  13. Redon, R.; Ishikawa, S.; Fitch, K.R.; Feuk, L.; Perry, G.H.; Andrews, T.D.; Fiegler, H.; Shapero, M.H.; Carson, A.R.; Chen, W.; et al. Global variation in copy number in the human genome. Nature 2006, 444, 444–454. [Google Scholar] [CrossRef] [PubMed]
  14. McCarroll, S.A.; Kuruvilla, F.G.; Korn, J.M.; Cawley, S.; Nemesh, J.; Wysoker, A.; Shapero, M.H.; de Bakker, P.I.; Maller, J.B.; Kirby, A.; et al. Integrated detection and population-genetic analysis of SNPs and copy number variation. Nat. Genet. 2008, 40, 1166–1174. [Google Scholar] [CrossRef] [PubMed]
  15. Iafrate, A.J.; Feuk, L.; Rivera, M.N.; Listewnik, M.L.; Donahoe, P.K.; Qi, Y.; Scherer, S.W.; Lee, C. Detection of large-scale variation in the human genome. Nat. Genet. 2004, 36, 949–951. [Google Scholar] [CrossRef] [PubMed]
  16. Inaki, K.; Liu, E.T. Structural mutations in cancer: Mechanistic and functional insights. Trends Genet. 2012, 28, 550–559. [Google Scholar] [CrossRef] [PubMed]
  17. Crook, T.; Wrede, D.; Tidy, J.A.; Mason, W.P.; Evans, D.J.; Vousden, K.H. Clonal p53 mutation in primary cervical cancer: Association with human-papillomavirus-negative tumours. Lancet 1992, 339, 1070–1073. [Google Scholar] [CrossRef]
  18. McIntyre, J.B.; Wu, J.S.; Craighead, P.S.; Phan, T.; Kobel, M.; Lees-Miller, S.P.; Ghatage, P.; Magliocco, A.M.; Doll, C.M. PIK3CA mutational status and overall survival in patients with cervical cancer treated with radical chemoradiotherapy. Gynecol. Oncol. 2013, 128, 409–414. [Google Scholar] [CrossRef] [PubMed]
  19. Narayan, G.; Murty, V.V. Integrative genomic approaches in cervical cancer: Implications for molecular pathogenesis. Future Oncol. 2010, 6, 1643–1652. [Google Scholar] [CrossRef] [PubMed]
  20. Vazquez-Mena, O.; Medina-Martinez, I.; Juarez-Torres, E.; Barron, V.; Espinosa, A.; Villegas-Sepulveda, N.; Gomez-Laguna, L.; Nieto-Martinez, K.; Orozco, L.; Roman-Basaure, E.; et al. Amplified genes may be overexpressed, unchanged, or downregulated in cervical cancer cell lines. PLoS ONE 2012, 7, e32667. [Google Scholar] [CrossRef] [PubMed]
  21. http://hapmap.ncbi.nlm.nih.gov/downloads/samplesindividuals/00README.txt (accessed on 23 May 2016).
  22. Chanda, P.; Yuhki, N.; Li, M.; Bader, J.S.; Hartz, A.; Boerwinkle, E.; Kao, W.H.; Arking, D.E. Comprehensive evaluation of imputation performance in African Americans. J. Hum. Genet. 2012, 57, 411–421. [Google Scholar] [CrossRef] [PubMed]
  23. Massink, M.P.; Kooi, I.E.; van Mil, S.E.; Jordanova, E.S.; Ameziane, N.; Dorsman, J.C.; van Beek, D.M.; van der Voorn, J.P.; Sie, D.; Ylstra, B.; et al. Proper genomic profiling of (BRCA1-mutated) basal-like breast carcinomas requires prior removal of tumor infiltrating lymphocytes. Mol. Oncol. 2015, 9, 877–888. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Ali Hassan, N.Z.; Mokhtar, N.M.; Kok Sin, T.; Mohamed Rose, I.; Sagap, I.; Harun, R.; Jamal, R. Integrated analysis of copy number variation and genome-wide expression profiling in colorectal cancer tissues. PLoS ONE 2014, 9, e92553. [Google Scholar] [CrossRef]
  25. Alazawi, W.; Pett, M.; Strauss, S.; Moseley, R.; Gray, J.; Stanley, M.; Coleman, N. Genomic imbalances in 70 snap-frozen cervical squamous intraepithelial lesions: Associations with lesion grade, state of the HPV16 E2 gene and clinical outcome. Br. J. Cancer 2004, 91, 2063–2070. [Google Scholar] [CrossRef] [PubMed]
  26. Pinkel, D.; Albertson, D.G. Array comparative genomic hybridization and its applications in cancer. Nat. Genet. 2005, 37, 11–17. [Google Scholar] [CrossRef] [PubMed]
  27. Misra, A.; Pellarin, M.; Nigro, J.; Smirnov, I.; Moore, D.; Lamborn, K.R.; Pinkel, D.; Albertson, D.G.; Feuerstein, B.G. Array comparative genomic hybridization identifies genetic subgroups in grade 4 human astrocytoma. Clin. Cancer Res. 2005, 11, 2907–2918. [Google Scholar] [CrossRef] [PubMed]
  28. Rueda, O.M.; Diaz-Uriarte, R. Detection of recurrent copy number alterations in the genome: Taking among-subject heterogeneity seriously. BMC Bioinform. 2009, 10, 308. [Google Scholar] [CrossRef] [PubMed]
  29. http://crdd.osdd.net/raghava/ccdb/faq.php (accessed on 16 October 2015).
  30. Emily, M.; Mailund, T.; Hein, J.; Schauser, L.; Schierup, M.H. Using biological networks to search for interacting loci in genome-wide association studies. Eur. J. Hum. Genet. 2009, 17, 1231–1240. [Google Scholar] [CrossRef] [PubMed]
  31. Pan, W. Network-based model weighting to detect multiple loci influencing complex diseases. Hum. Genet. 2008, 124, 225–234. [Google Scholar] [CrossRef] [PubMed]
  32. Baranzini, S.E.; Galwey, N.W.; Wang, J.; Khankhanian, P.; Lindberg, R.; Pelletier, D.; Wu, W.; Uitdehaag, B.M.; Kappos, L.; Polman, C.H.; et al. Pathway and network-based analysis of genome-wide association studies in multiple sclerosis. Hum. Mol. Genet. 2009, 18, 2078–2090. [Google Scholar] [CrossRef] [PubMed]
  33. Correia, C.; Oliveira, G.; Vicente, A.M. Protein interaction networks reveal novel autism risk genes within GWAS statistical noise. PLoS ONE 2014, 9, e112399. [Google Scholar] [CrossRef] [PubMed]
  34. Holmans, P.; Green, E.K.; Pahwa, J.S.; Ferreira, M.A.; Purcell, S.M.; Sklar, P.; Owen, M.J.; O’Donovan, M.C.; Craddock, N. Gene ontology analysis of gwa study data sets provides insights into the biology of bipolar disorder. Am. J. Hum. Genet. 2009, 85, 13–24. [Google Scholar] [CrossRef] [PubMed]
  35. Stephens, P.J.; Greenman, C.D.; Fu, B.; Yang, F.; Bignell, G.R.; Mudie, L.J.; Pleasance, E.D.; Lau, K.W.; Beare, D.; Stebbings, L.A.; et al. Massive genomic rearrangement acquired in a single catastrophic event during cancer development. Cell 2011, 144, 27–40. [Google Scholar] [CrossRef] [PubMed]
  36. Stephens, P.J.; McBride, D.J.; Lin, M.L.; Varela, I.; Pleasance, E.D.; Simpson, J.T.; Stebbings, L.A.; Leroy, C.; Edkins, S.; Mudie, L.J.; et al. Complex landscapes of somatic rearrangement in human breast cancer genomes. Nature 2009, 462, 1005–1010. [Google Scholar] [CrossRef] [PubMed]
  37. Oh, E.K.; Kim, Y.W.; Kim, I.W.; Liu, H.B.; Lee, K.H.; Chun, H.J.; Park, D.C.; Oh, E.J.; Lee, A.W.; Bae, S.M.; et al. Differential DNA copy number aberrations in the progression of cervical lesions to invasive cervical carcinoma. Int. J. Oncol. 2012, 41, 2038–2046. [Google Scholar] [PubMed]
  38. Su, T.H.; Chang, J.G.; Yeh, K.T.; Lin, T.H.; Lee, T.P.; Chen, J.C.; Lin, C.C. Mutation analysis of CTNNB1 (beta-catenin) and AXIN1, the components of Wnt pathway, in cervical carcinomas. Oncol. Rep. 2003, 10, 1195–1200. [Google Scholar] [PubMed]
  39. Hu, L.; Zhu, Y.T.; Qi, C.; Zhu, Y.J. Identification of Smyd4 as a potential tumor suppressor gene involved in breast cancer development. Cancer Res. 2009, 69, 4067–4072. [Google Scholar] [CrossRef] [PubMed]
  40. Achuthan, R.; Bell, S.M.; Roberts, P.; Leek, J.P.; Horgan, K.; Markham, A.F.; MacLennan, K.A.; Speirs, V. Genetic events during the transformation of a tamoxifen-sensitive human breast cancer cell line into a drug-resistant clone. Cancer Genet. Cytogenet. 2001, 130, 166–172. [Google Scholar] [CrossRef]
  41. Pan, Z.; Chen, S.; Pan, X.; Wang, Z.; Han, H.; Zheng, W.; Wang, X.; Li, F.; Qu, S.; Shao, R. Differential gene expression identified in Uigur women cervical squamous cell carcinoma by suppression subtractive hybridization. Neoplasma 2010, 57, 123–128. [Google Scholar] [CrossRef] [PubMed]
  42. Kolbl, A.C.; Hiller, R.A.; Ilmer, M.; Liesche, F.; Heublein, S.; Schroder, L.; Hutter, S.; Friese, K.; Jeschke, U.; Andergassen, U. Glycosyltransferases as marker genes for the quantitative polymerase chain reaction-based detection of circulating tumour cells from blood samples of patients with breast cancer undergoing adjuvant therapy. Mol. Med. Rep. 2015, 12, 2933–2938. [Google Scholar] [PubMed]
  43. Dong, Y.; Cao, B.; Zhang, M.; Han, W.; Herman, J.G.; Fuks, F.; Zhao, Y.; Guo, M. Epigenetic silencing of NKD2, a major component of wnt signaling, promotes breast cancer growth. Oncotarget 2015, 6, 22126–22138. [Google Scholar] [CrossRef] [PubMed]
  44. Liu, J.; Lu, W.G.; Ye, F.; Cheng, X.D.; Hong, D.; Hu, Y.; Chen, H.Z.; Xie, X. Hes1/Hes5 gene inhibits differentiation via down-regulating Hash1 and promotes proliferation in cervical carcinoma cells. Int. J. Gynecol. Cancer 2010, 20, 1109–1116. [Google Scholar] [CrossRef] [PubMed]
  45. Liao, S.; Hartmaier, R.J.; McGuire, K.P.; Puhalla, S.L.; Luthra, S.; Chandran, U.R.; Ma, T.; Bhargava, R.; Modugno, F.; Davidson, N.E.; et al. The molecular landscape of premenopausal breast cancer. Breast Cancer Res. 2015, 17, 104. [Google Scholar] [CrossRef] [PubMed]
  46. Lu, Z.; Hu, X.; Li, Y.; Zheng, L.; Zhou, Y.; Jiang, H.; Ning, T.; Basang, Z.; Zhang, C.; Ke, Y. Human papillomavirus 16 E6 oncoprotein interferences with insulin signaling pathway by binding to tuberin. J. Biol. Chem. 2004, 279, 35664–35670. [Google Scholar] [CrossRef] [PubMed]
  47. Ma, X.F.; Zhang, J.; Shuai, H.L.; Guan, B.Z.; Luo, X.; Yan, R.L. IKKβ/NF-κβ mediated the low doses of bisphenol A induced migration of cervical cancer cells. Arch. Biochem. Biophys. 2015, 573, 52–58. [Google Scholar] [CrossRef] [PubMed]
  48. Wang, W.T.; Zhao, Y.N.; Yan, J.X.; Weng, M.Y.; Wang, Y.; Chen, Y.Q.; Hong, S.J. Differentially expressed microRNAs in the serum of cervical squamous cell carcinoma patients before and after surgery. J. Hematol. Oncol. 2014, 7, 6. [Google Scholar] [CrossRef] [PubMed]
  49. Tinker, A.V.; Ellard, S.; Welch, S.; Moens, F.; Allo, G.; Tsao, M.S.; Squire, J.; Tu, D.; Eisenhauer, E.A.; MacKay, H. Phase II study of temsirolimus (CCI-779) in women with recurrent, unresectable, locally advanced or metastatic carcinoma of the cervix. A trial of the NCIC Clinical Trials Group (NCIC CTG IND 199). Gynecol. Oncol. 2013, 130, 269–274. [Google Scholar] [CrossRef] [PubMed]
  50. Kuo, J.H.; Lin, C.W. Cellular uptake on N- and C-termini conjugated FITC of Rath cell penetrating peptides and its consequences for gene-expression profiling in U-937 human macrophages and HeLa cervical cancer cells. J. Drug Target. 2013, 21, 801–808. [Google Scholar] [CrossRef] [PubMed]
  51. Cao, X.Q.; Lu, H.S.; Zhang, L.; Chen, L.L.; Gan, M.F. MEKK3 and survivin expression in cervical cancer: Association with clinicopathological factors and prognosis. Asian Pac. J. Cancer Prev. APJCP 2014, 15, 5271–5276. [Google Scholar] [CrossRef] [PubMed]
Medsci 04 00012 g001 1024
Figure 1. Karyogram view of detected amplified and deleted regions across autosomes. Amplifications are shown at right side of the chromosomes, and deletions are shown at the left side of the chromosomes. The length of the horizontal bar corresponds to the number of samples observed at the respective cytobands. Most of the amplifications were found in cases at the short arms of chromosomes 1, 5, 8, 16, 19 and 20, as well as the centromeres of chromosomes 14, 15 and 21. Most of the deletions in the cases were mainly observed at the short arms of chromosomes 1, 4, 5, 7, 11, 12, 16, 17 and 19, as well as the centromere of chromosome 19. We observed a statistically-significant difference of CNA burden between case and reference groups (Table 2) for different sizes of CNAs (>100 kb, 10–100 kb and 1–10 kb). For example, statistically-significant differences of >100 kb, 10–100 kb and 1–10 kb deletions were observed between reference (6.7%, 2.5%, 0.8%), pre-cancer (92.5%, 81.5%, 87.7%) and cancer (89.0%, 76.7%, 86.6%) groups, respectively (Table 2).
Figure 1. Karyogram view of detected amplified and deleted regions across autosomes. Amplifications are shown at right side of the chromosomes, and deletions are shown at the left side of the chromosomes. The length of the horizontal bar corresponds to the number of samples observed at the respective cytobands. Most of the amplifications were found in cases at the short arms of chromosomes 1, 5, 8, 16, 19 and 20, as well as the centromeres of chromosomes 14, 15 and 21. Most of the deletions in the cases were mainly observed at the short arms of chromosomes 1, 4, 5, 7, 11, 12, 16, 17 and 19, as well as the centromere of chromosome 19. We observed a statistically-significant difference of CNA burden between case and reference groups (Table 2) for different sizes of CNAs (>100 kb, 10–100 kb and 1–10 kb). For example, statistically-significant differences of >100 kb, 10–100 kb and 1–10 kb deletions were observed between reference (6.7%, 2.5%, 0.8%), pre-cancer (92.5%, 81.5%, 87.7%) and cancer (89.0%, 76.7%, 86.6%) groups, respectively (Table 2).
Medsci 04 00012 g001
Medsci 04 00012 g002 1024
Figure 2. Plot of principal component analysis (PCA) and hierarchical clustering of copy number variation (CNV) or CNA datasets. (A) PCA scatter plot of three diagnostic groups (pre-cancer, cancer and reference). Each point represents a specific sample. Points are colored by group status, with blue representing pre-cancer (CIN-I, -II and -III), red representing invasive cancer and green representing references. (B) PCA scatter plot of four diagnostic groups. Each point represents a specific sample. Points are colored by group status, with blue representing low grade (CIN-I), red representing high grade (CIN- II and -III), purple representing invasive cancer and green representing references.
Figure 2. Plot of principal component analysis (PCA) and hierarchical clustering of copy number variation (CNV) or CNA datasets. (A) PCA scatter plot of three diagnostic groups (pre-cancer, cancer and reference). Each point represents a specific sample. Points are colored by group status, with blue representing pre-cancer (CIN-I, -II and -III), red representing invasive cancer and green representing references. (B) PCA scatter plot of four diagnostic groups. Each point represents a specific sample. Points are colored by group status, with blue representing low grade (CIN-I), red representing high grade (CIN- II and -III), purple representing invasive cancer and green representing references.
Medsci 04 00012 g002
Table
Table 1. Clinical demographics with diagnosis, age for cases and reference subjects.
Table 1. Clinical demographics with diagnosis, age for cases and reference subjects.
Subject IDAge (years)Diagnosis 1Diagnosis 2Diagnosis 3
1442CIN-Ilow grade dysplasiapre-cancer
1624CIN-Ilow grade dysplasiapre-cancer
1753CIN-Ilow grade dysplasiapre-cancer
1531CIN-Ilow grade dysplasiapre-cancer
1227CIN-IIhigh grade dysplasiapre-cancer
839CIN-IIIhigh grade dysplasiapre-cancer
1121CIN-IIIhigh grade dysplasiapre-cancer
630CIN-IIIhigh grade dysplasiapre-cancer
538CIN-IIIhigh grade dysplasiapre-cancer
355SCCSCCSCC
444SCCSCCSCC
251SCCSCCSCC
731SCCSCCSCC
142SCCSCCSCC
n = 125 femalesNAreferencereferencereference
All subjects in the cases were human papillomavirus (HPV) positive and from the Mexican American population. Reference subjects were from an admixed population. NA, not available. CIN, cervical intraepithelial neoplasia; SCC, squamous cell carcinoma.
Table
Table 2. Copy number alteration (CNA) burden (deletion and amplification) among different diagnostic groups.
Table 2. Copy number alteration (CNA) burden (deletion and amplification) among different diagnostic groups.
2 A
Size of CNAsReference (n = 122)Low (n = 4)High (n = 5)SCC (n = 5)
AllAmplification (%)Deletion (%)AllAmplification (%)Deletion (%)AllAmplification (%)Deletion (%)AllAmplification (%)Deletion (%)
>100 kb212198 (93.3)14 (6.7)15615 (9.6)141 (90.4)24415 (6.1)229 (93.9)21824 (11)194 (89)
10–100 kb157153 (97.5)4 (2.5)11924 (20.2)95 (79.8)18332 (17.5)151 (82.5)16338 (23.3)125 (76.7)
1 kb–10 kb128127 (99.2)1 (0.8)12511 (8.8)114 (91.2)19328 (14.5)165 (85.5)17924 (13.4)155 (86.6)
2 B
Size of CNAsReferencePre-Cancer (n = 9)Cancer (n = 5)
AllAmplification (%)Deletion (%)AllAmplification (%)Deletion (%)AllAmplification (%)Deletion (%)
>100 kb212198 (93.3)14 (6.7)40030 (7.5)370 (92.5)21824 (11)194 (89.0)
10–100 kb157153 (97.5)4 (2.5)30256 (18.5)246 (81.5)16338 (23.3)125 (76.7)
1 kb–10 kb128127 (99.2)1 (0.8)31839 (12.3)279 (87.7)17924 (13.4)155 (86.6)
Table
Table 3. Top 20 CNA chromosome regions with recurrent somatic CNAs (deletions >100 kb) in the cases; not observed in the reference group.
Table 3. Top 20 CNA chromosome regions with recurrent somatic CNAs (deletions >100 kb) in the cases; not observed in the reference group.
CytobandStartStop# of Sub. with Del*.Deletions in SubjectsNo. of MarkersLength (bp)Cancer-Related Genes
19p13.3766,4632,357,62510I (2), II (1), III (3), SCC (4)14971,591,163STK11
19p13.32,860,4293,984,02510I (2), II (1), III (3), SCC (4)11541,123,597MATK, DAPK3
16p13.3895,9861,883,51110I (2), II (1), III (3), SCC (4)1022987,526AXIN1
19p13.33,984,0254,738,5369I (2), II (1), III (3), SCC (3)672754,512
1p36.33-1p36.321,825,5102,545,3689I (2), II (0), III (3), SCC (4)732719,859HES5
1p36.322,545,3683,230,8489I (2), II (0), III (3), SCC (4)678685,481PRDM16
16p13.3377,287895,98610I (2), II (1), III (3), SCC (4)525518,700AXIN1, MSLN
1p36.331,282,6541,799,9759I (2), II (0), III (3), SCC (4)242517,322
16p13.32,816,6993,271,9358I (2), II (0), III (3), SCC (3)469455,237BSG, FGF22, AXIN1
7p22.393,965534,5038I (2), II (0), III (3), SCC (3)265440,539
19p13.3366,342762,49710I (2), II (1), III (3), SCC (4)604396,156HBA2
16p13.32,408,9632,706,3798I (2), II (0), III (3), SCC (3)142297,417
16p13.384,130377,28710I (2), II (1), III (3), SCC (4)284293,158AXIN1
16p13.32,132,6512,397,39110I (2), II (1), III (3), SCC (4)141264,741TSC2
1p36.331,021,3801,282,6549I (2), II (0), III (3), SCC (4)233261,275
16p13.31,888,3162,115,23110I (2), II (1), III (3), SCC (4)210226,916
1p36.33810,8361,021,3809I (2), II (0), III (3), SCC (4)221210,545
19p13.32,532,4932,735,49910I (2), II (1), III (3), SCC (4)228203,007
19p13.32,367,8582,532,49310I (2), II (1), III (3), SCC (4)191164,636
19p13.3254,506365,33410I (2), II (1), III (3), SCC (4)108110,829
# of Sub. with Del *: number of subjects with deletion. Serine/threonine kinase 11 (STK11); maturase K (MATK); Death-associated protein kinase 3 (DAPK3); axis inhibition protein 1 (AXIN1); hairy/enhancer of split, drosophila, homolog of, 5 (HES5); PR domain containing 16 (PRDM16); mesothelin (MSLN); basignin (BSG); fibroblast growth factor 22 (FGF22); haemoglobin A2 (HBA2); tuberous sclerosis complex 2 (TSC2).
Table
Table 4. Newly-identified CNAs (deletions) contain known oncogenes in cases only.
Table 4. Newly-identified CNAs (deletions) contain known oncogenes in cases only.
Query PositionGene NameGene IDLength (bp)Deletions in Different Diagnostic Groups
chr4:1694786–1937978WHSC17468243,193ref (0), I (1), II (0), III (3), SCC (1)
chr4:1937978–2141625WHSC17468203,648ref (0), I (1), II (0), III (3), SCC (1)
chr19:6194240–6239714MLLT1429845,475ref (0), I (2), II (0), III (2), SCC (1)
chr19:6239714–6598484MLLT14298358,771ref (0), I (2), II (0), III (2), SCC (1)
chr19:6679537–6796128VAV17409116,592ref (0), I (2), II (0), III (1), SCC (1)
chr19:6796128–6883617VAV1740987,490ref (0), I (2), II (0), III (1), SCC (1)
Wolf-Hirschhorn syndrome candidate 1 (WHSC1); myeloid/lymphoid leukemia; translocated to 1 (MLLT1); vav guanine nucleotide exchange factor 1 (VAV1).
Table
Table 5. Newly-identified CNAs (deletions) contain known tumor suppressor genes in cases only.
Table 5. Newly-identified CNAs (deletions) contain known tumor suppressor genes in cases only.
Query PositionGene NameGene IDLength (bp)Deletions in Different Diagnostic Groups
chr1:3252272–3613519TP737161361,248ref (0), I (1), II (0), III (3), SCC (3)
chr1:3613519–3615102TP7371611584ref (0), I (1), II (0), III (3), SCC (3)
chr1:3615102–3615341TP737161240ref (0), I (1), II (0), III (3), SCC (3)
chr1:3615341–3621322TP7371615982ref (0), I (1), II (0), III (3), SCC (3)
chr1:3621322–4125280TP737161503,959ref (0), I (1), II (0), III (2), SCC (2)
chr5:204517–422812AHRR57491218,296ref (0), I (1), II (0), III (3), SCC (2)
chr5:422812–582058AHRR57491159,247ref (0), I (1), II (0), III (3), SCC (2)
chr7:1216332–1963055MAD1L18379746,724ref (0), I (2), II (0), III (3), SCC (3)
chr7:1963055–1964889MAD1L183791835ref (0), I (2), II (0), III (3), SCC (3)
chr7:1964889–2219927MAD1L18379255,039ref (0), I (2), II (0), III (3), SCC (3)
chr7:2219927–2770642MAD1L18379550,716ref (0), I (2), II (0), III (3), SCC (3)
chr16:1888316–2115231TSC27249226,916ref (0), I (2), II (1), III (3), SCC (4)
chr16:2115231–2116240TSC272491010ref (0), I (2), II (1), III (3), SCC (4)
chr16:2116240–2127184TSC2724910,945ref (0), I (2), II (1), III (3), SCC (4)
chr16:2127184–2132651TSC272495468ref (0), I (2), II (1), III (3), SCC (4)
chr16:2132651–2397391TSC27249264,741ref (0), I (2), II (1), III (3), SCC (4)
chr16:377287–895986AXIN18312518,700ref (0), I (2), II (1), III (3), SCC (4)
chr16:84130–377287AXIN18312293,158ref (0), I (2), II (1), III (3), SCC (4)
chr17:1586256–1692127SMYD4114826105,872ref (0), I (1), II (0), III (3), SCC (1)
chr17:1692127–1737406SMYD411482645,280ref (0), I (1), II (0), III (3), SCC (1)
chr19:5310257–5598484SAFB29667288,228ref (0), I (2), II (0), III (3), SCC (3)
chr19:5598484–5623377SAFB2966724,894ref (0), I (2), II (0), III (3), SCC (3)
chr19:5598484–5623377SAFB629424,894ref (0), I (2), II (0), III (3), SCC (3)
chr19:5623377–5625242SAFB62941,866ref (0), I (2), II (0), III (3), SCC (3)
chr19:5625242–5625729SAFB6294488ref (0), I (2), II (0), III (3), SCC (3)
chr19:5625729–5846816SAFB6294221,088ref (0), I (2), II (0), III (3), SCC (3)
Reference (0), I (1), II (0), III (3), SCC (3) means a 361-kb deletion on chromosome 1 occurred in one subject with CIN-I, 3 subjects with CIN-III and 3 subjects with SCC.
Table
Table 6. Functional enrichment of Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis in the identified CNAs (>100 kb and occurring in two or more patients with SCC).
Table 6. Functional enrichment of Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis in the identified CNAs (>100 kb and occurring in two or more patients with SCC).
Pathway NameEnrichment ScoreEnrichment p Value% Genes in Pathway# Gene in ListKEGG Pathway ID
Insulin signaling pathway4.7490.0097.03916
Endometrial cancer4.6980.00910.64511
Parkinson’s disease4.0540.0177.227215
Glioma3.7130.0248.335182
Huntington’s disease3.6400.0265.84952
NF-κB signaling3.5720.0287.146163
Ribosome3.5220.0307.066275
Alcoholism3.4950.0305.709191
Neurotrophin signaling pathway3.3990.0336.31750
Estrogen signaling pathway3.3260.0366.756154
RNA polymerase3.1950.04110.713238
Dopaminergic synapse3.0420.0485.83720
PI3K-Akt signaling pathway3.0380.0484.5214262
Ubiquitin-mediated proteolysis3.0050.0495.797175
# ≥6 genes in each function were listed here.
Table
Table 7. Gene ontology (GO) categories of CNAs (>100 kb and occurring in two or more than two patients with SCC).
Table 7. Gene ontology (GO) categories of CNAs (>100 kb and occurring in two or more than two patients with SCC).
TypeFunctionEnrichment ScoreEnrichment p Value% Genes in Group# Genes in ListGO ID
Molecular functionsignal transducer activity5.3440.0050.47774871
molecular transducer activity4.7800.0080.430760089
enzyme binding3.0890.0460.376519899
protein kinase binding5.8580.0030.998419901
kinase binding5.4050.0040.879419900
Cellular componentmembrane-bounded vesicle4.4200.0120.338931988
vesicle4.1930.0150.326931982
plasma membrane3.0610.0470.27195886
plasma membrane part3.3880.0340.328744459
plasma membrane protein complex5.5000.0040.903498797
membrane protein complex3.2030.0410.457498796
Biological processsignal transduction4.6480.0100.291127165
cell surface receptor signaling pathway4.5460.0110.34497166
regulation of signal transduction3.4430.0320.33279966
positive regulation of gene expression3.9920.0180.413610628
macromolecular complex subunit organization3.3030.0370.354643933
G-protein coupled receptor signaling pathway4.0890.0170.48957186
negative regulation of response to stimulus3.8210.0220.457548585
defense response3.7830.0230.45256952
positive regulation of protein metabolic process3.5690.0280.428551247
# ≥6 genes in each function were listed here.

Share and Cite

MDPI and ACS Style

Torabi, A.; Ordonez, J.; Su, B.B.; Palmer, L.; Mao, C.; Lara, K.E.; Rubin, L.P.; Xu, C. Novel Somatic Copy Number Alteration Identified for Cervical Cancer in the Mexican American Population. Med. Sci. 2016, 4, 12. https://doi.org/10.3390/medsci4030012

AMA Style

Torabi A, Ordonez J, Su BB, Palmer L, Mao C, Lara KE, Rubin LP, Xu C. Novel Somatic Copy Number Alteration Identified for Cervical Cancer in the Mexican American Population. Medical Sciences. 2016; 4(3):12. https://doi.org/10.3390/medsci4030012

Chicago/Turabian Style

Torabi, Alireza, Javier Ordonez, Brenda Bin Su, Laura Palmer, Chunxiang Mao, Katherine E. Lara, Lewis P. Rubin, and Chun Xu. 2016. "Novel Somatic Copy Number Alteration Identified for Cervical Cancer in the Mexican American Population" Medical Sciences 4, no. 3: 12. https://doi.org/10.3390/medsci4030012

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