Differential DNA copy number aberrations in the progression of cervical lesions to invasive cervical carcinoma

Immune landscape and heterogeneity of cervical squamous cell carcinoma and adenocarcinoma

Despite the differences in disease outcomes and pathological features between cervical squamous cell carcinoma (CSCC) and adenocarcinoma (ADC), the molecular characteristics in immune heterogeneity of the tumor microenvironment remain unclear. Here, we explored the immune landscape and heterogeneity between CSCC and ADC. Gene expression and clinical characteristics of cervical carcinoma from The Cancer Genome Atlas (TCGA) were downloaded. Differentially expressed genes (DEGs), immune cell infiltration, and pathway enrichment analyses were used to explore the immune landscape and heterogeneity between CSCC and ADC. Furthermore, distinct immune signatures between CSCC and ADC were validated based on clinical samples. In total, 4,132 upregulated DEGs and 2,307 down-regulated DEGs were identified between CSCC and ADC, with enrichments in immune related-pathways in CSCC. In addition, 54 hub DEGs correlated with patients' prognosis and immunocytes infiltration were identified. The CSCC patients had a higher ImmuneScore and more abundant immunocytes infiltration compared to ADC patients, as validated by immunohistochemistry (IHC) and multicolor immunofluorescence (mIF) analyses of collected samples. Furthermore, CSCC displayed higher inhibitory immune checkpoints expression, tumor mutation burden (TMB), and microsatellite instability (MSI) compared to ADC, which indicated CSCC patients were more likely to benefit from immunotherapy. In summary, our results revealed the huge immune heterogeneity between CSCC and ADC, and provided guidance for immunotherapy selection for different pathological types of cervical cancer.

Comprehensive profile and contrastive analysis of circular RNA expression in cervical squamous carcinoma and adenocarcinoma

Background

Numerous studies have shown circular RNA (circRNA) dysregulation is associated with the pathogenesis of cervical cancer,particularly in individual carcinoma variants. The aim of this study is to investigate and contrastively analyze the expression pattern of circRNAs in cervical squamous carcinoma and adenocarcinoma mediated by human papillomavirus type 16 (HPV-16).

Methods

The expression of circRNAs in cervical squamous carcinoma (SCC), adenocarcinoma (ADC) and adenosquamous carcinoma (ASC) tissues, together with the adjacent normal tissues (ANT), was profiled by high-throughput RNA sequencing (RNA-seq). Bioinformatics analysis and quantitative real time polymerase chain reaction (qRT-PCR) validation of the sequencing data were performed. A network of circRNA-miRNA (microRNA)-mRNA was then constructed according to predicted targets and function of candidate circRNAs.

Results

A total of 11,685 annotated circRNAs were identified in six cervical samples. There were 42 up-regulated and 98 down-regulated circRNAs. 215 circRNAs were up-regulated in SCC but down-regulated circRNAs in ADC, while 50 circRNAs displayed the opposite trend. Function enrichment analysis based on different expressions of circRNAs found that the most enriched pathway in all the three pathologic variants of cervical cancer was the "ubiquitin mediated proteolysis" pathway. Eight key candidate circRNAs derived from this pathway were further validated, and we noticed that several target miRNAs of candidate circRNAs could target the source genes. Based on this we constructed a related competing endogenous RNA (ceRNA) network.

Conclusion

Through a comprehensive interpretation of differentially expressed circRNAs in different pathologic variants of cervical cancer, this study provides new insights into the process of tumor differentiation mediated by HPV. Our results may help to complement the molecular typing and stem cell theory of cervical cancer.

Folate Deficiency Triggers the Abnormal Segregation of a Region With Large Cluster of CG-Rich Trinucleotide Repeats on Human Chromosome 2

Folate deficiency is associated with a broad range of human disorders, including anemia, fetal neural tube defects, age-associated dementia and several types of cancer. It is well established that a subgroup of rare fragile sites (RFSs) containing expanded CGG trinucleotide repeat (TNR) sequences display instability when cells are deprived of folate. However, given that folate sensitive RFSs exist in a very small percentage of the population, they are unlikely to be the cause of the widespread health problems associated with folate deficiency. We hypothesized that folate deficiency could specifically affect DNA replication at regions containing CG-rich repeat sequences. For this, we identified a region on human chromosome 2 (Chr2) comprising more than 300 CG-rich TNRs (termed “FOLD1”) by examining the human genome database. Via the analysis of chromosome shape and segregation in mitosis, we demonstrate that, when human cells are cultured under folate stress conditions, Chr2 is prone to display a “kink” or “bending” at FOLD1 in metaphase and nondisjunction in anaphase. Furthermore, long-term folate deprivation causes Chr2 aneuploidy. Our results provide new evidence on the abnormalities folate deficiency could cause in human cells. This could facilitate future studies on the deleterious health conditions associated with folate deficiency.

The Role of Predictive Biomarkers in Endocervical Adenocarcinoma: Recommendations From the International Society of Gynecological Pathologists

To review the scientific evidence related to predictive biomarkers in cervical adenocarcinoma (ADC). The authors reviewed the literature regarding predictive biomarkers in cervical ADC. There were several limitations: (1) there is an overlap between predictive and prognostic biomarkers, as the vast majority of patients are treated with anticancer strategies; (2) in many studies and clinical trials, cervical ADC patients are included in a large series of patients predominantly composed of cervical squamous cell carcinomas; and (3) in most of the studies, and clinical trials, there is no distinction between human papillomavirus (HPV)-associated and HPV-independent cervical ADCs, or between various histologic subtypes. Results obtained from a small group of studies confirm that cervical ADCs exhibit distinct molecular features as compared with squamous carcinomas, and that there are different molecular features between different types of cervical ADCs. Promising areas of interest include ERBB2 (HER2) mutations and PD-L1 expression as predictive biomarkers for anti-HER2 treatment and immunotherapy, respectively. To date, no definitive data can be obtained from the literature regarding predictive biomarkers for cervical ADC. Clinical trials specifically designed for endocervical ADC patients are required to elucidate the predictive value of HER2 mutations and PD-L1 expression. The distinction between HPV-associated and HPV-independent cervical ADCs as well as early involvement of pathologists in the design of future clinical trials are needed to identify new predictive biomarkers in cervical ADC.

Overexpression of SEZ6L2 predicts poor prognosis in patients with cholangiocarcinoma

Background

With the feature of destructive and biliary malignancy, intrahepatic cholangiocarcinoma (ICC), presents unclear molecular mechanisms which contributes to typically poor prognosis for patients. Seizure-related 6 homolog-like 2 (SEZ6L2) is a gene that encodes for a seizure-associated protein localized on the cell surface. Thus far, the function of SEZ6L2 in ICC has not been reported.

Methods

We used data from The Cancer Genome Atlas and the Gene Expression Omnibus to analyze dynamics behind and levels of expression of SEZ6L2 in ICC. Then we used qRT-PCR and Immunohistochemical staining to detect levels of expression of SEZ6L2 and thereby determined the potential clinical significance of this protein in ICC.

Results

According to qRT-PCR and immunohistochemical analysis results, SEZ6L2 was overexpressed in ICC. Kaplan-Meier and Cox proportional hazard analyses indicated that patients afflicted by ICC with high levels of relative expression of SEZ6L2 have a poorer prognosis and that SEZ6L2 may be an independent prognostic factor which enables to the accurate prediction of overall survival (OS) and disease-free survivals’ (DFS) expected rates. Subcutaneous xenograft models used to explore the role of SEZ6L2 in tumor formation in vivo. The dynamics of the SEZ6L2 gene being promote angiogenesis in cholangiocarcinoma are related to increasing expressive growth factors which include EGF, VEGF, PDGF and the activation of the P38-MAPK pathway.

Conclusions

Our findings suggest that SEZ6L2 can serve as an advanced biomarker that can be used to accurately predict a patient prognosis and be used as a target for ICC treatment.

Genetic variants of DNAH11 and LRFN2 genes and their association with ovarian and breast cancer

OBJECTIVE

To investigate the association of newly identified genetic variants G>A (rs2285947) of the DNAH11 gene and G>A (rs2494938) of the LRFN2 gene with ovarian and breast cancers in women belonging to Jammu and Kashmir state, where the prevalence of ovarian and breast cancers is remarkably high in the population.

METHODS

A candidate gene prospective case-control association study design was adopted, in which 354 cases (219 cases of ovarian cancer and 135 cases of breast cancer) were histopathologically confirmed and 330 healthy controls matched for age and ethnicity were recruited. The details of cases and controls were also recorded in a predesigned pro forma after their written informed consent. Both variants were genotyped by TaqMan allele discrimination assay using real-time polymerase chain reaction. Logistic regression analysis was performed to estimate the corrected odds ratio (OR), confidence interval (CI), and level of significance (P value) for potential confounding factors.

RESULTS

The rs2285947 variant of DNAH11 was found to be significantly associated with both ovarian and breast cancers with adjusted ORs of 1.7 (95% CI 1.2-2.4; P=0.004) and 1.70 (95% CI 1.13-2.54; P=0.0009), respectively. However, no significant association of variant rs2494938 of LRFN2 was observed with ovarian cancer (estimated OR 0.9, 95% CI 0.6-1.4; P=0.919) or breast cancer (estimated OR 1.27, 95% CI 0.8-1.9; P=0.216).

CONCLUSIONS

The collected data proposed that the variant rs2285947 of DNAH11 gene is a potential risk factor for ovarian and breast cancers in the studied population.

Copy number variation in archival melanoma biopsies versus benign melanocytic lesions

BACKGROUND

Skin melanocytes can give rise to benign and malignant neoplasms. Discrimination of an early melanoma from an unusual/atypical benign nevus can represent a significant challenge. However, previous studies have shown that in contrast to benign nevi, melanoma demonstrates pervasive chromosomal aberrations.

OBJECTIVE

This substantial difference between melanoma and benign nevi can be exploited to discriminate between melanoma and benign nevi.

METHODS

Array-comparative genomic hybridization (aCGH) is an approach that can be used on DNA extracted from formalin-fixed paraffin-embedded (FFPE) tissues to assess the entire genome for the presence of changes in DNA copy number. In this study, high resolution, genome-wide single-nucleotide polymorphism (SNP) arrays were utilized to perform comprehensive and detailed analyses of recurrent copy number aberrations in 41 melanoma samples in comparison with 21 benign nevi.

RESULTS

We found statistically significant copy number gains and losses within melanoma samples. Some of the identified aberrations are previously implicated in melanoma. Moreover, novel regions of copy number alterations were identified, revealing new candidate genes potentially involved in melanoma pathogenesis.

CONCLUSIONS

Taken together, these findings can help improve melanoma diagnosis and introduce novel melanoma therapeutic targets.

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

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.

Serine/threonine kinases 31(STK31) may be a novel cellular target gene for the HPV16 oncogene E7 with potential as a DNA hypomethylation biomarker in cervical cancer

Background

Cervical cancer (CC) is a leading cause of mortality in females, especially in developing countries. The two viral oncoproteins E6 and E7 mediate the oncogenic activities of high-risk human papillomavirus (hrHPV), and hrHPV, especially HPV16 or/and HPV18 (HPV16/18) play critical roles in CC through different pathways. STK31 gene of which the expression has been proven to be regulated by the methylation status of its promoter, is one of the novel cancer/testis (CT) genes and plays important roles in human cancers. Reasearches have indicated that viral infection is correlated to the methylation statuses of some genes. Herein, we detected methylation status of the STK31 gene in cervical tumors and explored its interaction with HPV16 or/and HPV18 (HPV16/18) infection.

Methods

Bisulfite genomic sequencing PCR (BGS) combined with TA clone, methylation-specific PCR (MSP) were used to analyze methylation statuses of the STK31 gene promoter/exon 1 region in HPV16/18-positive, HPV-negative CC cell lines; ectopically expressed HPV16 E6, -E7, and -E6/E7 CC cells; normal cervical tissues and cervical tumor tissues of different stages. The mRNA and protein expressions of STK31 were detected by RT-PCR and western blotting.

Results

The STK31 gene promoter/exon 1 was hypomethylated in the HPV16/18-positive cell lines HeLa, SiHa and CaSki, and the mRNA and protein expression were detected. In contrast, the STK31 gene exhibited hypermethylation and silenced expression in the HPV-negative CC cells C33A and HT-3. Compared with the primary HPV-negative CC cell lines, the STK31 methylation was downregulated, and STK31 expression was induced in the HPV16E7/E67 transfected cells. The methylation statuses and expressions of STK31 were verified in the cervical tumor samples at different stages. Additionally, chemotherapy treatment may influence STK31 expression by regulating its methylation status.

Conclusions

STK31 may be a novel cellular target gene for the HPV16 oncogeneE7. The HPV16 oncogene E7 may affect STK31 expression through a methylation-mediated mechanism. The aberrant methylation of the STK31 promoter/exon 1 region may be a precursor of human cervical carcinogenesis and a potential DNA aberrant methylation biomarker of conditions ranging from precancerous disease to invasive cancer.

Electronic supplementary material

The online version of this article (doi:10.1186/s12985-016-0515-5) contains supplementary material, which is available to authorized users.

Progression of naive intraepithelial neoplasia genome to aggressive squamous cell carcinoma genome of uterine cervix

Although cervical intraepithelial neoplasia (CIN) is considered a neoplasia, its genomic alterations remain unknown. For this, we performed whole-exome sequencing and copy number profiling of three CINs, a microinvasive carcinoma (MIC) and four cervical squamous cell carcinomas (CSCC). Both total mutation and driver mutation numbers of the CINs were significantly fewer than those of the MIC/CSCCs (P = 0.036 and P = 0.018, respectively). Importantly, PIK3CA was altered in all MIC/CSCCs by either mutation or amplification, but not in CINs. The CINs harbored significantly lower numbers of copy number alterations (CNAs) than the MIC/CSCCs as well (P = 0.036). Pathway analysis predicted that the MIC/CSCCs were enriched with cancer-related signalings such as cell adhesion, mTOR signaling pathway and cell migration that were depleted in the CINs. The mutation-based estimation of evolutionary ages identified that CIN genomes were younger than MIC/CSCC genomes. The data indicate that CIN genomes harbor unfixed mutations in addition to human papilloma virus infection but require additional driver hits such as PIK3CA, TP53, STK11 and MAPK1 mutations for CSCC progression. Taken together, our data may explain the long latency from CIN to CSCC progression and provide useful information for molecular diagnosis of CIN and CSCC.

Copy Number Gains at 8q24 and 20q11-q13 in Gastric Cancer Are More Common in Intestinal-Type than Diffuse-Type

The present study was aimed at discovering DNA copy number alterations (CNAs) involved in the carcinogenesis of stomach and at understanding their clinicopathological significances in the Korean population. DNA copy numbers were analyzed using Agilent 244K or 400K array comparative genomic hybridization (aCGH) in fresh-frozen tumor and matched normal tissues from 40 gastric cancer patients. Some of the detected CNA regions were validated using multiplex ligation-dependent probe amplification (MLPA) in six of the 40 patients and customized Agilent 60K aCGH in an independent set of 48 gastric cancers. The mRNA levels of genes at common CNA regions were analyzed using quantitative real-time PCR. Copy number gains were more common than losses across the entire genome in tumor tissues compared to matched normal tissues. The mean number of alterations per case was 64 for gains and 40 for losses, and the median aberration length was 44016 bp for gains and 4732 bp for losses. Copy number gains were frequently detected at 7p22.1 (20%), 8q24.21 (27%-30%), 8q24.3 (22%-48%), 13q34 (20%-31%), and 20q11-q13 (25%-30%), and losses at 3p14.2 (43%), 4q35.2 (27%), 6q26 (23%), and 17p13.3 (20%-23%). CNAs at 7p22.1, 13q34, and 17p13.3 have not been reported in other populations. Most of the copy number losses were associated with down-regulation of mRNA levels, but the correlation between copy number gains and mRNA expression levels varied in a gene-dependent manner. In addition, copy number gains tended to occur more commonly in intestinal-type cancers than in diffuse-type cancers. In conclusion, the present study suggests that copy number gains at 8q24 and 20q11-q13 and losses at 3p14.2 may be common events in gastric cancer but CNAs at 7p22.1, 13q34, and 17p13.3 may be Korean-specific.

Genetic variants at 6p21.1 and 7p15.3 Identified by GWASs of multiple cancers and ovarian cancer risk: a case-control study in Han Chinese women

A recent study summarized several published genome-wide association studies (GWASs) of cancer and reported two pleiotropic loci at 6p21.1 and 7p15.3 contributing to multiple cancers including lung cancer, noncardia gastric cancer (NCGC), and esophageal squamous-cell carcinoma (ESCC) in Han Chinese. However, it is not known whether such genetic variants have similar effects on the risk of gynecologic cancers, such as ovarian cancer. Hence, we explored associations between genetic variants in 6p21.1 and 7p15.3 and ovarian cancer risk in Han Chinese women. We performed an independent case-control study by genotyping the two loci (rs2494938 A > G at 6p21.1 and rs2285947 A > G at 7p15.3) in a total of 377 ovarian cancer cases and 1,034 cancer-free controls using TaqMan allelic discrimination assay. We found that rs2285947 at 7p15.3 was significantly associated with risk of ovarian cancer with per allele odds ratio (OR) of 1.33 [95% confidence interval (CI): 1.08-1.64, P=0.008]. However, no significant association was observed between rs2494938 and ovarian cancer risk. Our results showed that rs2285947 at 7p15.3 may also contribute to the development of ovarian cancer in Han Chinese women, further suggesting pleiotropy of 7p15.3 in multiple cancers.

Chromosomal gains and losses in human papillomavirus-associated neoplasia of the lower genital tract - a systematic review and meta-analysis

BACKGROUND

Overexpression of the human papillomavirus (HPV) oncogenes E6 and E7 is necessary for the development of distinct lower genital tract cancers. However, secondary cellular genomic alterations are mandatory to promote progression of HPV-induced premalignant stages. We aimed at identifying the chromosomal regions most frequently gained and lost and the disease stage at which the latter occurs. These regions might be relevant for carcinogenesis and could serve as diagnostic markers to identify premalignant lesions with high progression risk towards invasive cancer.

METHODS

We performed a systematic literature review and meta-analysis of studies listed in PubMed that analysed chromosomal copy number alterations by comparative genomic hybridisation (CGH) in HPV-positive and -negative cancers or premalignant lesions of the anogenital tract (cervix, anus, vagina, penis and vulva).

FINDINGS

Data were extracted and analysed from 32 studies. The most common alterations in cervical squamous cell carcinoma (SCC) (12 studies, 293 samples) were gains at 3q with a rate of 0.55 (95% confidence interval (CI) 0.43-0.70), losses at 3p (0.36, 95%CI 0.27-0.48) and losses at 11q (0.33, 95%CI 0.26-0.43). Gains at 3q were particularly frequent in HPV16-positive cervical SCC (0.84, 95%CI 0.78-0.90). Also more than one quarter of high grade cervical intraepithelial neoplasia (CIN) harboured gains of 3q (0.27, 95%CI 0.20-0.36), but the rate in low grade CIN was low (0.02, 95%CI 0.00-0.09). For HPV-associated vulvar SCC (four studies, 30 samples) the same common alterations as in cervical SCC were reported. Studies on non-cervical and non-vulvar SCC and premalignant lesions of the lower genital tract are scarce.

INTERPRETATION

3q gains were most frequently found in HPV16-positive cervical SCC. The results suggest the selection of HPV-transformed cell clones harbouring 3q gains in high grade premalignant lesions, while alterations in low grade lesions are rare.