Identification of susceptibility loci for cervical carcinoma by genome scan of affected sib-pairs

Oral cancers preceded by proliferative verrucous leukoplakia exhibit distinctive molecular features

OBJECTIVE

Proliferative verrucous leukoplakia (PVL) has high rates of malignant transformation into oral squamous cell carcinoma (OSCC), but the clinical and evolutionary pattern of OSCC from PVL (PVL-OSCC) is more favorable than that of OSCC not preceded by PVL (OSCC). Here, we aimed to explore the pathophysiologic differences between PVL-OSCC and OSCC through transcriptomic and DNA methylation analyses.

MATERIALS AND METHODS

In this case-control study, oral biopsies from 8 PVL-OSCC and 10 OSCC patients were obtained for global sequencing using RNAseq and a genome-wide DNA methylation analysis via the Infinium EPIC Platform (graphical abstract).

RESULTS

One hundred and thirty-three differentially expressed genes (DEGs) were detected, 94 of them upregulated in OSCC. Most of these genes were previously described in cancer and associated with prognosis. The integrative analysis revealed 26 DEGs, corresponding to 37 CpGs, whose promoters were regulated by DNA methylation. Twenty-nine of the CpGs were found as hypermethylated in PVL-OSCC. Only 5 of the genes that were aberrantly methylated and differentially expressed were upregulated in PVL-OSCC patients, whereas 21 were underexpressed.

CONCLUSIONS

PVL-OSCC patients presented lower expression of cancer-related genes. Hypermethylation of the promoter region of many genes was also noticed, indicating that DNA methylation could be a regulatory mechanism.

Critical Role of the Transcription Factor AKNA in T-Cell Activation: An Integrative Bioinformatics Approach

The human akna gene encodes an AT-hook transcription factor, the expression of which is involved in various cellular processes. The goal of this study was to identify potential AKNA binding sites in genes that participate in T-cell activation and validate selected genes. Here we analyzed ChIP-seq and microarray assays to determine AKNA-binding motifs and the cellular process altered by AKNA in T-cell lymphocytes. In addition, we performed a validation analysis by RT-qPCR to assess AKNA's role in promoting IL-2 and CD80 expression. We found five AT-rich motifs that are potential candidates as AKNA response elements. We identified these AT-rich motifs in promoter regions of more than a thousand genes in activated T-cells, and demonstrated that AKNA induces the expression of genes involved in helper T-cell activation, such as IL-2. The genomic enrichment and prediction of AT-rich motif analyses demonstrated that AKNA is a transcription factor that can potentially modulate gene expression by recognizing AT-rich motifs in a plethora of genes that are involved in different molecular pathways and processes. Among the cellular processes activated by AT-rich genes, we found inflammatory pathways potentially regulated by AKNA, suggesting AKNA is acting as a master regulator during T-cell activation.

The Human Papillomavirus (HPV) E6 Oncoprotein Regulates CD40 Expression via the AT-Hook Transcription Factor AKNA

Persistent infection with high-risk Human Papillomavirus (HR-HPV) is the main requisite for cervical cancer development. Normally, HPV is limited to the site of infection and regulates a plethora of cellular elements to avoid the immune surveillance by inducing an anti-inflammatory state, allowing the progress through the viral cycle and the carcinogenic process. Recent findings suggest that the AT-hook transcriptional factor AKNA could play a role in the development of cervical cancer. AKNA is strongly related to the expression of co-stimulatory molecules such CD40/CD40L to achieve an anti-tumoral immune response. To date, there is no evidence demonstrating the effect of the HPV E6 oncoprotein on the AT-hook factor AKNA. In this work, minimal expression of AKNA in cervical carcinoma compared to normal tissue was found. We show the ability of E6 from high-risk HPVs 16 and 18 to interact with and down-regulate AKNA as well as its co-stimulatory molecule CD40 in a proteasome dependent manner. We also found that p53 interacts with AKNA and promotes AKNA expression. Our results indicate that the de-regulation of CD40 and AKNA is induced by the HPV E6 oncoprotein, and this event involves the action of p53 suggesting that the axis E6/p53A/AKNA might play an important role in the de-regulation of the immune system during the carcinogenic process induced by HR-HPV.

IL-10 gene promoter and intron polymorphisms as genetic biomarkers of cervical cancer susceptibility among Tunisians

OBJECTIVE

We investigated the association between polymorphisms in the promoter and intron regions of the interleukin-10 (IL-10) gene with the risk of cervical cancer (CC) in Tunisian patients and control women.

METHODS

Study subjects comprised 86 CC cases and 126 control women. Genotyping of IL-10 intron (rs3024491, rs3024490) and promoter (rs1800872, rs1800871, rs1800896) variants was done by real-time PCR, with defined clusters.

RESULTS

The minor allele frequencies of the five tested IL-10 SNPs were not significantly different between cervical cancer cases and control women. However, significantly higher frequencies of homozygous minor allele-carriers in cases was seen for rs3024490 (P=0.023), rs1800872 (P=0.037), and rs1800871 (P=0.028). IL-10 serum levels were significantly reduced in rs3024490 T/T vs. G/G genotype carriers, and in rs1800871 T/T than C/C genotype carriers. While carriage of rs1800872 and rs3024491 minor allele was associated with reduced IL-10 secretion, this was not statistically significant. Haploview analysis demonstrated high linkage disequilibrium (LD) among the IL10 SNPs studied, and only seven haplotypes were common, capturing 98.8% of the total possible haplotypes. Reduced frequency of haplotypes GTCCA (P<0.001) and TGATG (P<0.001) was seen in cervical cancer cases than in control women, thus conferring disease protection nature to these haplotype. This association remained significant for GTCCA (Pc=0.006) and TGATG (P=0.045) after correcting for multiple comparisons.

CONCLUSION

Specific IL-10 variants (rs3024490, rs1800872, and rs1800871) and haplotype (GTCCA and TGATG) may contribute to the development of cervical cancer among Tunisian women.

Expression Analyses Revealed Thymic Stromal Co-Transporter/Slc46A2 Is in Stem Cell Populations and Is a Putative Tumor Suppressor

By combining conventional single cell analysis with flow cytometry and public database searches with bioinformatics tools, we extended the expression profiling of thymic stromal cotransporter (TSCOT), Slc46A2/Ly110, that was shown to be expressed in bipotent precursor and cortical thymic epithelial cells. Genome scale analysis verified TSCOT expression in thymic tissue- and cell type- specific fashion and is also expressed in some other epithelial tissues including skin and lung. Coexpression profiling with genes, Foxn1 and Hoxa3, revealed the role of TSCOT during the organogenesis. TSCOT expression was detected in all thymic epithelial cells (TECs), but not in the CD31(+) endothelial cell lineage in fetal thymus. In addition, ABC transporter-dependent side population and Sca-1(+) fetal TEC populations both contain TSCOT-expressing cells, indicating TEC stem cells express TSCOT. TSCOT expression was identified as early as in differentiating embryonic stem cells. TSCOT expression is not under the control of Foxn1 since TSCOT is present in the thymic rudiment of nude mice. By searching variations in the expression levels, TSCOT is positively associated with Grhl3 and Irf6. Cytokines such as IL1b, IL22 and IL24 are the potential regulators of the TSCOT expression. Surprisingly, we found TSCOT expression in the lung is diminished in lung cancers, suggesting TSCOT may be involved in the suppression of lung tumor development. Based on these results, a model for TEC differentiation from the stem cells was proposed in context of multiple epithelial organ formation.

Cervical cancer-associated promoter polymorphism affects akna expression levels

Cervical cancer (CC) is responsible for >260,000 deaths worldwide each year. Efforts are being focused on identifying genetic susceptibility factors, especially in genes related to the immune response. Akna has been proposed to be one of them, but data regarding its functional role in the disease is scarce. Supporting the notion of akna as a CC susceptibility gene, we found two polymorphisms associated with squamous intraepithelial lesion (SIL) and CC; moreover, we identified an association between high akna expression levels and CC and SIL, but its direction differs in each disease stage. To show the potential existence of a cis-acting polymorphism, we assessed akna allelic expression imbalance for the alleles of the -1372C>A polymorphism. We found that, regardless of the study group, the number of transcripts derived from the A allele was significantly higher than those from the C allele. Our results support the hypothesis that akna is a CC susceptibility genetic factor and suggest that akna transcriptional regulation has a role in the disease. We anticipate our study to be a starting point for in vitro evaluation of akna transcriptional regulation and for the identification of transcription factors and cis-elements regulating AKNA function that are involved in carcinogenesis.

The zinc finger protein ZNF268 is overexpressed in human cervical cancer and contributes to tumorigenesis via enhancing NF-κB signaling

Cervical cancer is one of the most common tumors affecting women's health worldwide. Although human papillomavirus can be detected in nearly all cases, the mechanism of cervical carcinogenesis remains to be further addressed. Here, we demonstrated that ZNF268, a Krüppel-associated box-containing zinc finger protein, might contribute to the development of cervical cancer. We found that ZNF268b2, an isoform of ZNF268, was overexpressed in human squamous cervical cancer specimens. Knockdown of ZNF268 in cervical cancer cells caused cell cycle arrest at the G(0)/G(1) phase, reduced colony formation, and increased sensitivity to TNFα-induced apoptosis. In addition, HeLa cell growth in xenograft nude mice was suppressed by ZNF268 knockdown, with increased apoptosis. Furthermore, ZNF268b2 was shown to increase NF-κB signaling in vitro and in vivo. Reconstitution of NF-κB activity restored proliferation in ZNF268 knockdown HeLa cells. Of note, we observed a high frequency of NF-κB activation in ZNF268-overexpressing cervical cancer tissues, suggesting a pathological coincidence of ZNF268b2 overexpression and NF-κB activation. Taken together, our results reveal a novel role of ZNF268b2 that contributes to cervical carcinogenesis in part through enhancing NF-κB signaling.

Correlation between IL-10 gene expression and HPV infection in cervical cancer: a mechanism for immune response escape

The present study was performed to determine IL-10 expression in cervical tissues in Mexican women according to the severity of the malignity and its association with HPV infection. IL-10 expression showed a clear tendency to increase during the different cervical cancer stages: 37% in LGSIL; 62% in HGSIL; and 84% in cancer. However, all the patients that expressed IL-10 were HPV positives; we found an association with HPV 16. These results suggest a clear relationship between IL-10, HPV and the stage of cervical cancer disease; this event could contribute to the immunosuppressive micro-environment in the tumor site.

SNP frequency estimation using massively parallel sequencing of pooled DNA

Resequencing of genomic regions that have been implicated by linkage and/or association studies to harbor genetic susceptibility loci represents a necessary step to identify causal variants. Massively parallel sequencing (MPS) offers the possibility of SNP discovery and frequency determination among pooled DNA samples. The strategies of pooling DNA samples and pooling PCR amplicons generated from individual DNA samples were evaluated, and both were found to return accurate estimates of SNP frequencies across varying levels of sequence coverage.

Genetic studies on the APOA1-C3-A5 gene cluster in Asian Indians with premature coronary artery disease

BACKGROUND

The APOA1-C3-A5 gene cluster plays an important role in the regulation of lipids. Asian Indians have an increased tendency for abnormal lipid levels and high risk of Coronary Artery Disease (CAD). Therefore, the present study aimed to elucidate the relationship of four single nucleotide polymorphisms (SNPs) in the Apo11q cluster, namely the -75G>A, +83C>T SNPs in the APOA1 gene, the Sac1 SNP in the APOC3 gene and the S19W variant in the APOA5 gene to plasma lipids and CAD in 190 affected sibling pairs (ASPs) belonging to Asian Indian families with a strong CAD history.

METHODS & RESULTS

Genotyping and lipid assays were carried out using standard protocols. Plasma lipids showed a strong heritability (h2 48% - 70%; P < 0.0001). A subset of 77 ASPs with positive sign of Logarithm of Odds (LOD) score showed significant linkage to CAD trait by multi-point analysis (LOD score 7.42, P < 0.001) and to Sac1 (LOD score 4.49) and -75G>A (LOD score 2.77) SNPs by single-point analysis (P < 0.001). There was significant proportion of mean allele sharing (pi) for the Sac1 (pi 0.59), -75G>A (pi 0.56) and +83C>T (pi 0.52) (P < 0.001) SNPs, respectively. QTL analysis showed suggestive evidence of linkage of the Sac1 SNP to Total Cholesterol (TC), High Density Lipoprotein-cholesterol (HDL-C) and Apolipoprotein B (ApoB) with LOD scores of 1.42, 1.72 and 1.19, respectively (P < 0.01). The Sac1 and -75G>A SNPs along with hypertension showed maximized correlations with TC, TG and Apo B by association analysis.

CONCLUSION

The APOC3-Sac1 SNP is an important genetic variant that is associated with CAD through its interaction with plasma lipids and other standard risk factors among Asian Indians.

The power of two-locus affected sib-pair linkage analysis to detect interacting disease loci

It has been shown that two-locus linkage analysis can, for some two-locus disease models, be used to detect effects at disease loci that do not reach significance in a genome scan. However, few examples exist where two-locus linkage has been successfully used to map genes. We study the possible gain in power of affected sib-pair nonparametric two-locus linkage analysis for two-locus models which fulfil the two-locus triangle constraints. Using a new parameterization of the two-locus joint identity-by-descent sharing probabilities we can, for fixed marginal sharing at both of two unlinked disease loci, derive a two-locus distribution such that the power of a two-locus analysis is maximized. In a simulation study we look at two test statistics, the two-locus maximum likelihood score and the correlation between nonparametric linkage scores, and study power as a function of marginal sharing. We show that in a best-case scenario two-locus linkage can have considerable power to detect pairs of interacting loci if there is a moderate increase in allele sharing at one of the two loci, even if there is a very small increase in allele sharing at the other locus. But we also show that the power to detect interacting loci in a two-locus analysis decreases as the marginal sharing at the two loci decreases and for any distribution with a small increase in allele sharing at both loci the power of a two-locus analysis is always low.