{"paper_id":"6c36b650-d227-417b-aaf0-afb5c0f31a38","body_text":"Cervical cancer is a common malignant tumor in women worldwide  1 . According to the American Cancer Society, cervical cancer is the second most common cancer-related death among females aged 20 to 39  2 . In China, age-standardized mortality has increased in the past decade  3 . Chronic infection with high-risk human papillomavirus (HPV) promotes the progression from pre-cancerous lesions to benign cervical dysplasia and cervical intraepithelial neoplasia (CIN) and finally to cervical cancer  4 .\nApolipoprotein B mRNA-editing enzyme catalytic (APOBEC) is a family of cytidine deaminases that converts cytosines in DNA and RNA to uracil, resulting in DNA/RNA mutations or damage  5 . Deep sequencing demonstrated that APOBEC mutations are critical in HPV-induced carcinogenesis and are prevalent in cervical cancer tissues and cancer cell lines  6 - 9 . APOBEC3B expression in cervical cancer is associated with HPV-18 infection  10  and is often associated with cell proliferation, migration, and cell cycle progression in cervical cancers, except for gastric-type cervical adenocarcinoma  11 - 13 . APOBEC3A overexpression leads to deamination and degradation of transfected exogenous plasmid DNA, suggesting that APOBEC3A may mediate the clearance of HPV DNA in persistently infected cells  14 .\nAPOBEC3A is the enzyme with the highest catalytic activity among APOBEC3 proteins and is a significant source of APOBEC-signature mutations  15 . APOBEC3A-induced cytidine deamination-associated somatic mutations are enriched in cervical cancer genomes, where the C-to-T mutational signature predominates  16 ,  17 . C-to-T or G-to-A mutations also exist in the mitochondrial DNA of cervical cancer cells  18 . Uracil-DNA glycosylase (UNG) is a critical molecule in the nucleic acid base excision repair pathways, which can recognize and excise the C-to-U mutation induced by APOBECs, triggering the hydrolysis of the nucleic acid chain for mismatch repair  19 . We suggested that the imbalance between APOBEC3B and UNG caused by chronic non-resolving inflammation may promote hepatocellular carcinoma  20 . APOBEC3A has potent deamination activity in breast cancer cell lines and is strongly associated with mutational burden in breast tumors, suggesting that APOBEC3A is a significant driver of genomic mutations in cancer of distinct histotypes  21 .\nIn contrast, APOBEC3A is a protective marker and a prognostic biomarker for predicting survival and immunotherapy response in ovarian cancer  22 . However, the roles of APOBEC3A in cervical cancer are unclear. Therefore, we explored the APOBEC3A expression pattern and its effects on biological behaviors such as proliferation, apoptosis, invasion, and migration in cervical cancer.\n\nThe RNA-Seq data of cervical squamous cell carcinoma and endocervical adenocarcinoma (CESC), adjacent normal samples (TCGA-CESC), and normal cervical samples (GTEx database) were downloaded from UCSC Xena ( https://xenabrowser.net/ ). The TCGA-CESC and GTEx cohorts were combined to form a larger cohort (TCGA and GTEx) consisting of 13 normal cervical tissues and 294 CESC samples. The  GSE7803  dataset was obtained from the GEO database ( https://www.ncbi.nlm.nih.gov/geo/ ). All downloaded patients had no other malignant tumors or diseases. No patient had undergone adjuvant treatment before surgery. The clinical phenotypes were grouped and analyzed. Gene set enrichment analysis (GSEA) plots were generated using the R packages  GSEABase ,  DESeq2 ,  clusterProfiler , and  enrichplot . Survival plots were drawn using the R packages  survivor  and  survminer , and box plots were drawn using the R package  ggplot2 . All analyses used the R Foundation for Statistical Computing (2021) version 4.1.1.\nWe used 77 cases from the First Affiliated Hospital of Second Military Medical University between January 2014 and January 2017. These included 38 patients with cervical cancer, 16 with CIN who underwent total hysterectomy or cervical conization, and 23 histologically confirmed normal patients who underwent total hysterectomy for symptomatic uterine fibroids or adenomyosis. All patients were hospitalized for surgery; no other malignant tumors were combined. No patient had undergone adjuvant treatment before surgery. Cervical formalin-fixed paraffin-embedded tissues were collected from each participant. The Second Military Medical University Ethics Committee approved the study, and all patients provided written informed consent.\nCervical cancer, CIN, and normal cervical tissues were processed using standard techniques. Anti-APOBEC3A (1:50 dilution; abcam, ab262853, San Diego, CA, USA) was used for IHC staining. Three investigators (Tan XJ, Yu YW, and Lu SL) evaluated the IHC staining independently. APOBEC3A-stained slides were scanned using Aperio ScanScope XT (Leica Biosystems) software, and the images were magnified 40 times for histological scoring, as reported previously  23 . In selected areas of the specimen, the Aperio Nuclear Algorithm ranked the immunostaining as negative (-), slightly positive (+), moderately positive (++), or strongly positive (+++). The APOBEC3A H-score was calculated using the linear formula: H-score = 1 × (% slightly positive cells) + 2 × (% moderately positive cells) + 3 × (% strongly positive cells). Two of the three investigators usually reached a close agreement (90%) on H-score. Otherwise, a consensus was reached by discussion.\nHeLa, SiHa, and HEK293T cells were purchased from the Cell Bank of the Chinese Academy of Sciences (Shanghai, China). HeLa and HEK293T cells were cultured in Dulbecco's modified Eagle's medium (DMEM) (BasalMedia, Shanghai, China), and SiHa cells were cultured in Minimum Essential Medium (GIBCO, Grand Island, NY, USA). The media were supplemented with 10% fetal bovine serum (Gibco, New York, NY, USA) and 1% penicillin/streptomycin (Bio-light, Shanghai, China). Cell cultures were performed in a humidified 5% CO 2  environment at 37 °C. Mycoplasma contamination test of cell cultures was performed every 3 months. All experiments were performed with mycoplasma-free cells. The APOBEC3A-expressing vector containing human APOBEC3A coding sequence in the pLenti-CMV-mCherry-PGK-blasticidin expression vector was purchased from Obio Technology (Shanghai, China). The structure diagrams of empty control and APOBEC3A expressing plasmids are shown in  Figure S1 . The APOBEC3A expressing and empty control plasmids were transfected using a Lipofectamine 3000 kit (Invitrogen, Carlsbad, CA, USA).\nWestern blotting was conducted as described  24 . The primary antibodies were anti-APOBEC3A (dilution 1:500, PA5-78800, Invitrogen) and anti-GAPDH (dilution 1:1000, sc-25778, Santa Cruz, Dallas, TX, USA) rabbit polyclonal antibodies. GAPDH was used as the loading and normalization control. ImageJ software (version 1.51,  https://imagej.nih.gov/ij/ ) was used to quantify the signal strength of the bands. We performed at least three independent experiments.\nA CCK8 kit (Dojindo, Osaka, Japan) was used to assess cell proliferation. Cells were seeded into 96-well plates (Corning Incorporated Coster, Kennebunk, ME) at 2000 cells per well. Every 24 h, 10 µL of CCK8 was added to wells and incubated for 2 h before measuring optical density at 450 nm. For cell cycle assay, cells were digested and fixed at -20 °C in 75% ethanol overnight. Propidium iodide/RNase staining buffer (BD, San Diego, CA, USA) measured the cell cycle distribution using flow cytometry (Merck Millipore, Rockville, MD, USA). ModFit LT software (Verity Software House, Topsham, ME) analyzed the raw data. For apoptosis assay, cells were harvested in 0.25% trypsin without EDTA (Gibco, Grand Island, NY, USA) and resuspended in binding buffer (BD Pharmingen, San Diego, CA, USA), stained with Annexin V-FITC and 7-AAD (BD Pharmingen, San Diego, CA, USA) according to the manufacturer's instructions. A Beckman Coulter CytoFlEX S flow cytometer (Beckman Coulter, Miami, FL) was used to measure apoptosis.\nThe  in vitro  wound-healing assay was performed to study directional cell migration. Briefly, HeLa and SiHa cells were seeded in six-well plates at 5 × 10 5  cells/well. After incubation for 24 h, APOBEC3A overexpression and empty control plasmids were transfected using Lipofectamine 3000 kit (Invitrogen). We created a scrape with a 100-µL pipette tip when the monolayers reached confluence. Cells were washed with phosphate-buffered saline, and images of migrated cells were captured at zero and 48 h using a fluorescence microscope (Leica, Germany). The wound-healing rate was analyzed using Image J software. The formula was % Cure Rate = (New Surface - Previous Surface) / New Surface × 100.\nCell invasion and migration were measured using Transwell inserts (Corning, NY, USA). We added 100 µL of Matrigel (Corning) vertically to the bottom center of the upper chamber, incubated at 37 °C for 5 h, and discarded the residual liquid for the invasion assay. The migration assay did not require Matrigel (Corning).\nWe seeded 4 × 10 5  cells into the upper chamber with 100 µL serum-free DMEM containing 0.1% bovine serum albumin (Bio-light, Shanghai, China). The lower chamber was filled with 500 µL DMEM containing 10% fetal bovine serum. After incubation for 24 h, the upper chamber was transferred to 500 μL of 4% paraformaldehyde (Bio-light) and fixed at room temperature for 15 min. Then, the upper chamber was transferred to crystalline violet staining solution (Bio-light, Shanghai, China), stained for 20 min, dried with a cotton swab, and photographed under the microscope. Five fields were randomly selected and photographed with a microscope at 10 × magnification. ImageJ software ( https://imagej.net/software/imagej/ ) was used to calculate the cell numbers. Each assay was performed in triplicate.\nThe APOBEC3A-expressing and empty control plasmids were transfected into HEK293T cells with the MD2G packaging plasmid and PAX2 envelope plasmid using a Lipofectamine 3000 kit (Invitrogen). After incubation for 48 h, lentivirus-containing supernatants were harvested, and lentiviral particles were concentrated using Lenti-X (Clontech, CA, USA). SiHa cells were infected with lentivirus for 24 h, and two stable cell strains (SiHa-APOBEC3A and SiHa-Vector) were screened using 10 μg/mL Blasticidin S (Yeasen, Shanghai, China).\nWe purchased 20 female mice (BALB/c-nude) from Cavens Laboratory Animal (Changzhou, China) to examine the tumorigenicity. Twenty mice were randomly assigned to two groups. SiHa-APOBEC3A cells (5 × 10 6 ) or SiHa-Vector cells (5 × 10 6 ) were mixed with Matrigel and injected into mice. The tumor volumes were measured twice weekly and calculated as volume = 0.5 × length × width 2 . After 38 days, mice were sacrificed, and xenograft tumors and mice were weighed. Tumors were harvested for hematoxylin-eosin staining. The Ethics Committee of the Second Military Medical University approved the protocol.\nHeLa and SiHa cells were seeded into 96-well plates (Corning) at 2000 cells per well. The APOBEC3A-expressing and empty control plasmids were transfected using a Lipofectamine 3000 kit (Invitrogen). After incubated for 48 h, the phosphorylated gamma-H2AX staining were performed using the DNA Damage Assay Kit by γ-H2AX Immunofluorescence (Beyotime Biotechnology, China) following the manufacturer's instructions. The images were taken on a fluorescence microscope (Leica, Germany) with an exposure time of 2s, and analyzed with ImageJ software ( https://imagej.net/software/imagej/ ).\nTotal RNA was isolated from two cell lines (SiHa-APOBEC3A and SiHa-Vector) using TRIzol (Invitrogen) and preserved at -80 °C. A NanoDrop spectrophotometer (Thermo Fisher Scientific, Waltham, MA, USA) determined the concentration, quality, and integrity. Three micrograms of RNA were used as input material for the RNA sample preparations. Sequencing libraries were generated using a TruSeq RNA Sample Preparation Kit (Illumina, San Diego, CA, USA). RNA-Seq and analysis were performed at Shanghai Personal Biotechnology Co., Ltd. (Shanghai, China).\nAnalyses were performed using GraphPad Prism version 8.3.0 (GraphPad Software, San Diego, CA, USA). The Student's  t -test was used to determine differences in continuous variables, while the chi-square test was used to evaluate differences in categorical variables. Transcriptome data and survival analysis were performed using R platform version 4.1.1 (Core Development Team,  http://www.r-project.org/ ). Survival plots were generated using the R packages  survivor  and  survminer . The box plots were created using the R package  ggplot2 . A two-sided  P  < 0.05 was considered statistically significant for all analyses.\n\nAPOBEC3A and UNG expression levels were significantly higher in cervical cancer than in adjacent normal tissues, according to TCGA and GTEx ( P  < 0.001, Figures  1 A and B). The expression of interleukin-6 (IL-6) and nuclear factor kappa-B (NF-κB) was significantly higher in tumor tissues than normal tissues ( P  = 0.021 and  P  = 0.00072, respectively), suggesting more inflammation in tumor tissues than in normal tissues. (Figures  1 C and D). According to  GSE7803 , APOBEC3A expression was significantly higher in cervical cancer than in intraepithelial lesion tissues ( P  = 0.001, Figure  1 E). The APOBEC3A to UNG expression level ratio was significantly higher in cervical cancer than in intraepithelial lesion tissues ( P  = 0.002, Figure  1 F).\nRepresentative IHC images were shown in Figure  2 A. IHC revealed that APOBEC3A expression was significantly higher in cervical cancer than in normal cervical tissues ( P  < 0.001), and APOBEC3A expression level in CIN tissues was higher than in normal cervical tissues ( P  < 0.001) (Figure  2 B). High-grade squamous intraepithelial lesions (HSILs) had higher APOBEC3A expression than low-grade squamous intraepithelial lesions (LSILs) ( P  < 0.01) (Figure  2 C). Poorly differentiated cervical cancer tissues had higher APOBEC3A expression than moderately well-differentiated cervical cancer ( P  < 0.001) (Figure  2 D).\nWe classified lesions above or equal to the median H-score as having high APOBEC3A expression (H-APOBEC3A) and those below the median as having low expression (L-APOBEC3A). High APOBEC3A expression was associated with advanced International Federation of Gynecology and Obstetrics (FIGO) staging (χ² = 12.931,  P  = 0.002), lymph node metastasis (χ² = 4.156,  P  = 0.041), and pathological histological hypofractionation (χ² = 5.911,  P  = 0.015) in cervical cancer patients (Supplementary  Table S1 ). After removing patients containing missing values, we also categorized the TCGA downloaded patients into high APOBEC3A expression (H-APOBEC3A) and low expression (L-APOBEC3A) groups (Supplementary  Table S2 ). The downloaded data did not contain the information of histological grade and tumor size. The results showed that APOBEC3A expression was not associated with FIGO stage (χ² =6.5951, P =0.08599) and lymph node metastasis (χ² = 3.6467, P =0.05618), which were inconsistent with the results of our own cohort. A total of 183 patients containing missing values in the downloaded data were excluded, which creates a serious bias. By comparing the common clinical parameters between our own cohort and the TCGA cohort, we found no significant statistical difference between the two cohorts (Supplementary  Table S3 ).\nHigher expression of APOBEC3A and UNG in cervical cancer was associated with better outcomes ( P  = 0.039 and  P  = 0.019, respectively) (Figures  2 E and F), according to TCGA. A higher ratio of APOBEC3A to UNG was associated with better outcomes ( P  = 0.04) (Figure  2 G).\nEctopic expression of APOBEC3A in SiHa and HeLa cells was demonstrated using western blots (Figure  3 A). The CCK8 assay indicated that the ectopic expression of APOBEC3A significantly inhibited proliferation in HeLa and SiHa cells (Figure  3 B and Figure  3 C). Flow cytometry results demonstrated that upregulating APOBEC3A led to a significant upregulation of cells at the S-phase and significant downregulation at the G2/M-phase (Figure  3 D- 3 G).\nEctopic expression of APOBEC3A inhibited the wound-healing rate at 24 and 48 h in HeLa and SiHa cells (Figure  4 ). The Transwell assay demonstrated that compared with HeLa-Vector and SiHa-Vector cells, HeLa-APOBEC3A and SiHa-APOBEC3A cells exhibited weaker invasion and migration ability, respectively (Figure  5 A- 5 D). Ectopic expression of APOBEC3A significantly increased cell apoptosis in SiHa and HeLa cells ( P  < 0.01) (Figure  5 E- 5 G). Staining for phosphorylated gamma histone 2AX (γ-H2AX) revealed that the ectopic expression of APOBEC3A significantly increased numbers of γ-H2AX-positive cells in in HeLa and SiHa cells ( P  < 0.0001) (Figure  6 ), suggesting that overexpression of APOBEC3A is cytotoxic.\nThe effects of APOBEC3A overexpression were also investigated  in vivo . Tumor volumes and weights were smaller in the APOBEC3A overexpression group than in the control group (Figure  7 A). At day 38, the average tumor size in the SiHa-APOBEC3A group was 0.83 ± 0.11 cm 3  and 1.49 ± 0.18 cm 3  in the SiHa-Vector group (Figure  7 C). The average tumor weight in the SiHa-APOBEC3A group was 0.80 ± 0.35 g, and 1.57 ± 0.59 g in the SiHa-Vector group (Figure  7 D). The representative hematoxylin-eosin staining pattern is shown in Figure  7 B.\nWe used RNA-Seq analysis to explore the pathways for APOBEC3A action in cervical cancer cells. There were 910 differentially expressed genes (fold change >2,  P <0.05) (Figure  8 A). Of these, 93.52% (851) were significantly downregulated in the SiHa-APOBEC3A cells compared with the SiHa-Vector cells. According to the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis, differentially expressed genes were involved in “Pathways in cancer” (KEGG-hsa05200) and were enriched in the phosphatidylinositide 3-kinases-protein kinase B (PI3K-Akt) signaling pathway (KEGG-hsa04151) (Figure  9 A). The significantly enriched gene ontology (GO) terms were related to biological processes, including anatomical structure morphogenesis (GO:0009653), system development (GO:0048731), and multicellular organism development (GO:0007275) (Figure  9 B). The PI3K-Akt signaling pathway, the mitogen-activated protein kinase (MAPK) signaling pathway, and pathways in cancer were significantly downregulated by APOBEC3A (Figure  8 B- 8 D). The top ten GSEA pathways are shown in Supplementary  Figure S2 .\n\nAPOBEC3A and pro-inflammatory factors IL-6 and NF-κB were more highly expressed in cervical cancer tissues than in adjacent normal tissues (Figure  1 ). APOBEC3A is a single-domain enzyme (i.e., only one zinc-dependent cytidine deaminase domain) and an influential restriction factor for HPV; APOBEC3 upregulation occurs throughout disease progression to suppress the infection  25 ,  26 .\nNF-κB activation upregulates the expression of APOBEC proteins in cervical cancer  27 . Under the persistent stimulation of inflammatory factors, APOBEC3 is significantly upregulated in hepatocellular carcinoma  28 . Our data support that HPV infection and higher inflammation levels in cervical cancer promote APOBEC3A expression.\nWe also observed an elevated level of UNG in cervical cancer tissues (Figure  1 B). APOBEC cytosine deamidation is managed by the base excision repair (BER) system. UNG-initiated BER is the primary mechanism to counteract ingAPOBEC3-induced genomic mutations  29 . Hence, the high level of APOBEC3A activates the BER system and upregulates UNG expression. However, the ratio of APOBEC3A to UNG is also increased in tumor tissues (Figure  1 F), suggesting that the DNA damage-repair balance is disturbed in cervical cancer.\nSomatic mutations are the basis of carcinogenesis in various cancers, and APOBEC3 is a significant source of somatic mutations. APOBEC3A-signature mutations are enriched in cervical cancer genome  16 . However, our results indicated that although APOBEC3A is significantly upregulated in cancer tissues, it exerts an anti-cancer effect and correlates with better outcomes in cervical cancer. These findings are in consistent with previous studies  30 ,  31 . But these previous related studies didn't explain this contradiction. Here, we presented and validated a hypothesis that may clarify this contradiction: Factors including HPV infection and high inflammation levels promote the expression of APOBEC3A in cervical cancer. However, the upregulation of APOBEC3A disturbs the DNA damage-repair balance because of its powerful mutagenic efficiency.\nThe APOBEC3A-mediated DNA damage exceeds the cellular tolerance, which promoting cancer cell apoptosis. Therefore, APOBEC3A conversely exerts an anti-cancer effect in cervical cancer. This hypothesis is also reported in renal cell carcinoma and biliary tract cancer  32 ,  33 . The cancer evolution-development hypothesis holds that in virus-associated cancers, the virus and host cells undergo a mutation-selection-adaptation evolutionary process; the balance between the strength of mutagenesis and the tolerance of cells to genotoxicity is critical for cancer evolution and development  15 . Steady-state APOBEC3A levels are very low in most cells because of tight regulation  34 . Ectopic APOBEC3A expression may hinder cancer evolution at the mutagenesis stage because of its powerful mutagenic role.\nAPOBEC3A interacts with genomic DNA, resulting in genomic instability, DNA mutations, breaks, and damage  35 . APOBEC3A functional polymorphisms increase APOBEC3A expression levels and significantly decrease the risk of renal cell carcinoma and biliary tract cancers. APOBEC3A overexpression exerts an anti-cancer effect by causing apoptosis in these cells  32 ,  33 . Our study also showed that ectopically expressed APOBEC3A inhibits migration and invasion, induces DNA double strand break, and promotes apoptosis. RNA-seq indicated that most genes were downregulated, and several cancer-related signaling pathways were inhibited when APOBEC3A was upregulated. These findings suggest that APOBEC3A inhibits majority of genes transcription by promoting DNA/RNA damage. APOBEC3A exerts anti-tumor effects in cervical cancer is partly because the APOBEC3A-mediated DNA damage exceeds the cellular tolerance, thus promoting cancer cell apoptosis.\nAPOBEC3A overexpression causes cell-cycle arrest in a deaminase-dependent manner  36 . The G2-phase checkpoint can prevent the initiation of mitosis and repair G2-phase and previously unrepaired DNA damage, primarily mediated by the ataxia-telangiectasia-mutated-and-Rad3-related kinase-checkpoint kinase 1 (ATR-Chk1) pathway. High APOBEC3A-expressing leukemia cells require the ATR-Chk1 checkpoint for survival, while blocking this pathway causes cell death  37 . Our study also showed that ectopically expressed APOBEC3A caused cell-cycle arrest at the S-phase. Our study demonstrated that ectopically expressed APOBEC3A caused cell-cycle arrest in the S-phase and downregulated the G2 stage, suggesting that APOBEC3A inhibits the ATR-Chk1 pathway in cervical cancer cells, leading to cell death.\nOur study has limitations. We focused on the effects of APOBEC3A on the biological behavior of cervical cancer cells and the relationship between APOBEC3A and outcomes. However, the somatic mutagenicity and antiviral effects of APOBEC3A are also essential for the occurrence and development of cervical cancer. First, HPV infection is a significant risk factor for cervical cancer. The SiHa cell line integrates a segment of the HPV16 genome (one to two copies per cell), while the HeLa cell line contains the HPV18 sequence (10 to 50 copies/cell)  38 . APOBEC3A plays an important role in the innate immune response to viral infection by editing the viral genome. The relationship between APOBEC3A and HPV in cervical cancer development remains to be further investigated. Second, APOBEC3A has potent deamination activity in cervical cancer cell lines and is strongly associated with mutational burden in cervical cancers. The complex roles of APOBEC3A-induced somatic mutations in tumorigenesis and development should be investigated. Third, genetic polymorphisms may also influence individual susceptibility to cervical cancer and its outcome. However, several studies reported that genetic polymorphisms of APOBEC3A is not significantly associated with cervical cancer risk, cervical carcinogenesis and prognosis  31 ,  39 - 40 . The association between the APOBEC3A polymorphism and cervical cancer needs to be further investigated.\n\nAPOBEC3A is significantly overexpressed in cervical cancer, and higher expression of APOBEC3A is associated with better outcomes. APOBEC3A inhibits cell migration and invasion, arrests cells in S-phase, promotes cervical cancer cell apoptosis by inducing DNA damage that exceeds cellular tolerance, and inhibits several cancer-associated signaling pathways.\n\nSupplementary figures and tables.\nClick here for additional data file.","source_license":"CC-BY-4.0","license_restricted":false}