Iwao Kukimoto

Molecular epidemiology of human papillomavirus 18 infections in Japanese Women

Among the oncogenic genotypes of human papillomavirus (HPV), HPV18 is the second most common type detected in cervical cancer worldwide and is primarily involved in the generation of cervical adenocarcinoma. Although HPV intra-type variants confer different risks of cervical carcinogenesis, there is little information on the genetic diversity of HPV18 compared to the most prevalent type, HPV16. In this study, we investigated the genetic variation of HPV18 in cervical specimens obtained from Japanese women with normal cervices or cervical cancers and precancers. Of the 101 HPV18-positive samples analyzed, viral whole genome amplification followed by next-generation sequencing led to the determination of viral complete genome sequences of 18 samples. Phylogenetic analysis of these HPV18 whole genome sequences identified a distinct variant cluster consisting of only Japanese samples (n = 7) belonging to sublineage A1. Viral genome sequences were also analyzed for the E6/E7 (n = 66) and E2 (n = 27) genes by Sanger sequencing. Phylogenetic analyses of these regions showed that the variant distribution among Japanese women was strongly biased toward sublineage A1 (72 of 87; 82.8%). No significant differences were observed in the prevalence of specific sublineages between cervical cancer/precancer cases and controls, and between squamous cell carcinoma and adenocarcinoma cases. These data contribute to our understanding of the genetic diversity of HPV18 in Japanese women.

Regional differences in human papillomavirus type 52 prevalence among Japanese women with cervical intraepithelial neoplasia

Abstract Although geographical differences in the distribution of human papillomavirus genotypes have been observed worldwide, no studies have reported on national differences in the prevalence of human papillomavirus types in Japan. Here, we report a cross-sectional study to explore regional differences in the prevalence of human papillomavirus types among Japanese women with cervical intraepithelial neoplasia or invasive cervical cancer. Using human papillomavirus genotyping data from the nationwide prospective study on human papillomavirus vaccine effectiveness, we compared the frequency of detection of 15 high-risk and two low-risk human papillomavirus types in each disease category between the women who visited hospitals located in eastern Japan and those who visited hospitals located in western Japan. The risk of cervical intraepithelial neoplasia progression was assessed by calculating a prevalence ratio of each human papillomavirus type for cervical intraepithelial neoplasia grade 2/3 versus grade 1. Among the human papillomavirus types studied, human papillomavirus 52 was detected significantly more frequently in western hospitals than in eastern hospitals in cervical intraepithelial neoplasia grade 1 patients, but was less frequent in cervical intraepithelial neoplasia grade 2/3. The prevalence of particular human papillomavirus types was not significantly different between patients in hospitals in eastern Japan and those in hospitals in western Japan for invasive cervical cancer. In both eastern and western hospitals, a higher risk of cervical intraepithelial neoplasia progression was observed in patients infected with human papillomavirus 16, 31 or 58. In contrast, there was a significantly higher prevalence of human papillomavirus 52 infection in women with cervical intraepithelial neoplasia grade 2/3 than in those with cervical intraepithelial neoplasia grade 1 in eastern hospitals (prevalence ratio, 1.93; 95% confidence interval, 1.48–2.58), but not in western hospitals (prevalence ratio, 1.03; 95% confidence interval, 0.83–1.30). Regional differences of human papillomavirus 52 prevalence in cervical intraepithelial neoplasia lesions may exist and emphasize the importance of continuous monitoring of human papillomavirus type prevalence throughout the country in order to accurately assess the efficacy of human papillomavirus vaccines.

Ancient Evolutionary History of Human Papillomavirus Type 16, 18 and 58 Variants Prevalent Exclusively in Japan

Human papillomavirus (HPV) is a sexually transmitted virus with an approximately 8-kilo base DNA genome, which establishes long-term persistent infection in anogenital tissues. High levels of genetic variations, including viral genotypes and intra-type variants, have been described for HPV genomes, together with geographical differences in the distribution of genotypes and variants. Here, by employing a maximum likelihood method, we performed phylogenetic analyses of the complete genome sequences of HPV16, HPV18 and HPV58 available from GenBank (n = 627, 146 and 157, respectively). We found several characteristic clusters that exclusively contain HPV genomes from Japan: two for HPV16 (sublineages A4 and A5), one for HPV18 (sublineage A1) and two for HPV58 (sublineages A1 and A2). Bayesian phylogenetic analyses of concatenated viral gene sequences showed that divergence of the most recent common ancestor of these Japan-specific clades was estimated to have occurred ~98,000 years before present (YBP) for HPV16 A4, ~39,000 YBP for HPV16 A5, ~38,000 YBP for HPV18 A1, ~26,000 for HPV58 A1 and ~25,000 YBP for HPV58 A2. This estimated timeframe for the divergence of the Japan-specific clades suggests that the introduction of these HPV variants into the Japanese archipelago dates back to at least ~25,000 YBP and provides a scenario of virus co-migration with ancestral Japanese populations from continental Asia during the Upper Paleolithic period.

Intra-Patient Genomic Variations of Human Papillomavirus Type 31 in Cervical Cancer and Precancer

Human papillomavirus type 31 (HPV31) is detected less frequently in cervical cancer than two major causative types, HPV16 and HPV18. Here, we report a comprehensive analysis of HPV31 genome sequences in cervical lesions collected from Japanese women. Of 52 HPV31-positive cervical specimens analyzed by deep sequencing, 43 samples yielded complete genome sequences of around 7900 base pairs and 9 samples yielded partially deleted genome sequences. Phylogenetic analysis showed that HPV31 variant distribution was lineage A in 19 samples (36.5%), lineage B in 28 samples (53.8%), and lineage C in 5 samples (9.6%), indicating that lineage B variants are dominant among HPV31 infections in Japan. Deletions in the viral genome were found in the region from the E1 to L1 genes, but all the deleted genomes retained the E6/E7 genes. Among intra-patient nucleotide variations relative to a consensus genome sequence in each sample, C-to-T substitutions were most frequently detected, followed by T-to-C and C-to-A substitutions. High-frequency, intra-patient mutations (>10%) in cervical cancer samples were found in the E1, E2, and E7 genes, and all of them were nonsynonymous substitutions. The enrichment of high-frequency nonsynonymous substitutions strongly suggests that these intra-patient mutations are positively selected during the development of cervical cancer/precancer.

Downregulation of HLA Class I Expression through HLA-A DNA Methylation Is Associated with Reduced CD8+ T-cell Infiltration in Cervical Cancer

Abstract Human leukocyte antigen class I (HLA-I) is central to tumor immune recognition, but its regulatory mechanisms in cervical cancer remain poorly understood. This study aimed to elucidate the impact of HLA-I regulatory mechanisms on CD8+ T-cell infiltration and identify distinct histotype-specific immune escape strategies across cervical cancer subtypes. Using 98 cervical cancer cases, including squamous cell carcinoma (SCC; n = 53), adenocarcinoma (n = 32), gastric-type adenocarcinoma (GAS; n = 5), small cell carcinoma (Small, n = 4), and mixed histologic types (n = 4), we examined the relationship between CD8+ T-cell infiltration patterns (categorized as infiltrated, excluded, or absent) and HLA-I expression, HLA-A DNA methylation, and HLA-I loss of heterozygosity (LOH). CD8+ T-cell infiltration patterns varied significantly by histologic subtype (P < 0.0001). SCC showed the highest frequency of the infiltrated pattern (73.6%), whereas GAS and Small predominantly displayed an absent pattern. Reduced CD8+ T-cell infiltration correlated with poor survival (P < 0.0001). HLA-I expression mirrored these trends being highest in SCC and lowest in Small and GAS. HLA-A DNA methylation emerged as a key driver of HLA-I downregulation, leading to reduced CD8+ infiltration (P < 0.05). In SCC, both HLA-A methylation and HLA-I LOH contributed to immune evasion; cases lacking these alterations exhibited the highest CD8+ T-cell infiltration levels (P < 0.01). This study identifies distinct HLA-I regulatory mechanisms in cervical cancer, highlighting HLA-A methylation—and particularly HLA-I LOH in SCC—as key drivers of immune evasion. These findings provide a foundation for developing predictive biomarkers and suggest that targeting these specific HLA-I regulatory mechanisms could enhance immunotherapy efficacy.