FOXP3 immunoregulatory gene variants are independent predictors of human papillomavirus infection and cervical cancer precursor lesions

Haplotype Structures and Protein Levels of TGFB1 in HPV Infection and Cervical Lesion: A Case-Control Study

This study aimed to verify the role of TGFB1 variants (c.-1638G>A, c.-1347C>T, c.29C>T, and c.74G>C) in HPV infection susceptibility and cervical lesions development, and their impact on TGFB1 cervical and plasma levels. TGFB1 genotypes were assessed with PCR-RFLP and haplotypes were inferred for 190 HPV-uninfected and 161 HPV-infected women. TGFB1 levels were determined with immunofluorimetric assay. Case-control analyses were performed with logistic regression adjusted for possible confounders. Women carrying -1347TT or -1347CT+TT as well as those with 29CT, 29CC, or 29CT+CC were more likely to have HPV than -1347CC and 29TT carriers, respectively. Regarding haplotypes, the most frequent were *4 (GCTG) and *3 (GTCG). Women *4/*4 were less likely to have HPV than those with no *4 copy. Comparing the inheritance of *3 and *4, carriers of *3/*4 or *3/*3 were more susceptible to HPV than *4/*4. The TGFB1 plasma and cervical levels were higher in the infected patients. Plasma levels were also higher in infected women with low-grade lesions. HPV-infected patients carrying *3/Other and *3/Other+*3/*3 presented lower TGFB1 plasma levels than those with no copy of *3. TGFB1 variants could contribute to the comprehension of the TGFB1 role in HPV-caused cervical disease.

Integration of scRNA-Seq and TCGA RNA-Seq to Analyze the Heterogeneity of HPV+ and HPV- Cervical Cancer Immune Cells and Establish Molecular Risk Models

Background

Numerous studies support that Human papillomavirus (HPV) can cause cervical cancer. However, few studies have surveyed the heterogeneity of HPV infected or uninfected (HPV+ and HPV-) cervical cancer (CESC) patients. Integration of scRNA-seq and TCGA data to analyze the heterogeneity of HPV+ and HPV- cervical cancer patients on a single-cell level could improve understanding of the cellular mechanisms during HPV-induced cervical cancer.

Methods

CESC scRNA-seq data obtained from the Gene Expression Omnibus (GEO) database and the Seurat, Monocle3 package were used for scRNA-seq data analysis. The ESTIMATE package was used for single-sample gene immune score, CIBERSORT package was used to identify immune scores of cells, and the “WGCNA” package for the weighted correlation network analysis. Univariate Cox and LASSO regression were performed to establish survival and relapse signatures. KEGG and GO analyses were performed for the signature gene. Gene Expression Profiling Interactive Analysis was used for Pan-cancer analysis.

Results

In the HPV+ CESC group, CD8+ T cells and B cells were down-regulated, whereas T reg cells, CD4+ T cells, and epithelial cells were up-regulated according to scRNA-seq data. Survival analysis of TCGA-CESC revealed that increased expression of naive B cells or CD8+ T cells favors the survival probability of CESC patients. WGCNA, univariate Cox, and LASSO Cox regression established a 9-genes survival signature and a 7-gene relapse model. Pan-cancer analysis identified IKZF3, FOXP3, and JAK3 had a similar distribution and effects in HPV-associated HNSC.

Conclusion

Analysis of scRNA-seq and bulk RNA-seq of HPV+ and HPV- CESC samples revealed heterogeneity from transcriptional state to immune infiltration. Survival and relapse models were adjusted according to the heterogeneity of HPV+ and HPV- CESC immune cells to assess the prognostic risk accurately. Hub genes represent similar protection in HPV- associated HNSC while showing irrelevant to other potential HPV-related cancers.

Genetic association study of intron variants in the forkhead box protein P3 gene in Chinese patients diagnosed with cervical cancer

The aim of this study was to investigate the effects of forkhead box protein P3 (FOXP3) intron single nucleotide variants (SNVs) in high-risk human papilloma virus (HR-HPV) infection and cervical cancer (CC) malignant lesions. We performed FOXP3 genotyping in 350 patients with CC and 350 healthy controls using the ImLDR multiple single nucleotide polymorphism genotyping technology. The heterozygous mutation TC in rs2294021 decreased the risk of HR-HPV infection and CC malignant lesions (TC vs. TT: OR = 0.71, 95% CI = 0.51-0.99); the dominant model TC+CC and allele C in rs2294021 decreased the risk of CC malignant lesions (TC+CC vs. TT: OR = 0.69, 95% CI = 0.50-0.95; C vs. T: OR = 0.78, 95% CI = 0.63-0.97). The heterozygous mutation GA, dominant model GA+AA and allele A in rs3761549 also decreased the risk of HR-HPV infection and CC malignant lesions (GA vs. GG: OR = 0.70, 95% CI = 0.51-0.96; GA+AA vs. GG: OR = 0.69, 95% CI = 0.51-0.94; A vs. G: OR = 0.75, 95% CI = 0.58-0.96). Patients with CC and HR-HPV infection carrying rs2294021 TC and rs3761549 GA had lower expression of FOXP3 protein. Haplotype analysis revealed that T-C-A decreased the risk of HR-HPV infection. Furthermore, we found a significant association between immune cells infiltration and prognosis in patients with CC. Our findings demonstrated that rs2294021 and rs3761549 variants may protect against HR-HPV and CC malignant lesions by downregulating FOXP3 and that FOXP3 was associated with immune cells infiltration, which affected the prognosis of CC.

Sex and FOXP3 gene rs2232365 polymorphism may be associated with the clinical and pathological aspects of chronic viral diseases

Background

The forkhead box protein 3 (FOXP3) transcription factor is one of the main markers of immunological suppression in different pathological profiles, and the presence of polymorphic variants may alter the gene expression of this factor. Despite descriptions of an association between the presence of the rs2232365 polymorphism and chronic diseases, the role of the sex variant in this context has not yet been elucidated, as the FOXP3 gene is located on the human sex chromosome X.

Results

To contribute to this topic, 323 women and 373 men were enrolled in the study, of which 101 were diagnosed with chronic viral liver diseases (39 women and 62 men), 67 with HTLV-1 infection (44 women and 23 men), 230 with coronary artery disease (91 women and 139 men) and 298 healthy and uninfected blood donors (149 women and men). They were genotyped for the rs2232365 polymorphism. The rs2232365 polymorphism was associated with clinical and pathological aspects and biomarkers of viral infections only in men, with functional differences between different infections.

Conclusions

A relationship is suggested between sex and FOXP3 rs2232365 polymorphism, resulting in different biological repercussions.

FOXP3 Genetic Variants Do Not Impact Circulating and Cervical Interleukin-10 Levels in Human Papillomavirus Infection in Women

Regulatory T cell (Treg) lineage plays a central role in inflammation and autoimmunity control. Interleukin-10 (IL-10) has been described as a pleiotropic cytokine that is mainly released by CD4⁺ CD25⁺ FOXP3⁺ Treg cells and has a potent immunosuppressive activity. Forkhead box P3 (FOXP3) transcription factor expression is crucial for Treg to function as a suppressor cell, and FOXP3 gene single nucleotide variants (SNVs) have already been shown to influence on viral pathogenesis. This study was conducted to evaluate the plasmatic and cervical levels of IL-10 in human papillomavirus-infected and uninfected patients and investigate whether the FOXP3 intron -1 SNVs rs3761548 and rs2232365 might alter IL-10 secretion. SNVs were genotyped by the characterization of polymerase chain reaction (PCR) products based on sequence-specific enzymatic cleavage using restriction fragment length polymorphism (RFLP) method. IL-10 levels were determined by quantitative enzyme-linked immunosorbent assay (ELISA). In conclusion, the data indicate that there is no association between FOXP3 SNVs and circulating and cervical IL-10 levels. This finding provides a rationale that IL-10 gene activation is independent of FOXP3 transcription factor activities on Treg cells.

Identification of key genes and pathways of diagnosis and prognosis in cervical cancer by bioinformatics analysis

Background

Cervical cancer as one of the most common malignant tumors lead to bad prognosis among women. Some researches already focus on the carcinogenesis and pathogenesis of cervical cancer, but it is still necessary to identify more key genes and pathways.

Methods

Differentially expressed genes were identified by GEO2R from the gene expression omnibus (GEO) website, then gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyzed by DAVID. Meanwhile, protein–protein interaction network was constructed by STRING, and both key genes and modules were found in visualizing network through Cytoscape. Besides, GEPIA did the differential expression of key genes and survival analysis. Finally, the expression of genes related to prognosis was further explored by UNLCAN, oncomine, and the human protein atlas.

Results

Totally 57 differentially expressed genes were founded, not only enriched in G1/S transition of mitotic cell cycle, mitotic nuclear division, and cell division but also participated in cytokine–cytokine receptor interaction, toll‐like receptor signaling pathway, and amoebiasis. Additionally, 12 hub genes and 3 key modules were screened in the Cytoscape visualization network. Further survival analysis showed that TYMS (OMIM accession number 188350), MCM2 (OMIM accession number 116945), HELLS (OMIM accession number 603946), TOP2A (OMIM accession number 126430), and CXCL8 (OMIM accession number 146930) were associated with the prognosis of cervical cancer.

Conclusion

This study aim to better understand the characteristics of some genes and signaling pathways about cervical cancer by bioinformatics, and could provide further research ideas to find new mechanism, more prognostic factors, and potential therapeutic targets for cervical cancer.