Comparison of four different human papillomavirus genotyping methods in cervical samples: Addressing method-specific advantages and limitations

Since human papillomavirus (HPV) is recognized as the causative agent of cervical cancer and associated with anogenital non-cervical and oropharyngeal cancers, the characterization of the HPV types circulating in different geographic regions is an important tool in screening and prevention. In this context, this study compared four methodologies for HPV detection and genotyping: real-time PCR (Cobas® HPV test), nested PCR followed by conventional Sanger sequencing, reverse hybridization (High + Low PapillomaStrip® kit) and next-generation sequencing (NGS) at an Illumina HiSeq2500 platform. Cervical samples from patients followed at the Family Health Strategy from Juiz de Fora, Minas Gerais, Brazil, were collected and subjected to the real-time PCR. Of those, 114 were included in this study according to the results obtained with the real-time PCR, considered herein as the gold standard method. For the 110 samples tested by at least one methodology in addition to real-time PCR, NGS showed the lowest concordance rates of HPV and high-risk HPV identification compared to the other three methods (67-75 %). Real-time PCR and Sanger sequencing showed the highest rates of concordance (97-100 %). All methods differed in their sensitivity and specificity. HPV genotyping contributes to individual risk stratification, therapeutic decisions, epidemiological studies and vaccine development, supporting approaches in prevention, healthcare and management of HPV infection.

Microbiome analysis of Brazilian women cervix reveals specific bacterial abundance correlation to RIG-like receptor gene expression

The relationship among microbiome, immunity and cervical cancer has been targeted by several studies, yet many questions remain unanswered. We characterized herein the virome and bacteriome from cervical samples and correlated these findings with innate immunity gene expression in a Brazilian convenience sample of HPV-infected (HPV+) and uninfected (HPV-) women. For this purpose, innate immune gene expression data were correlated to metagenomic information. Correlation analysis showed that interferon (IFN) is able to differentially modulate pattern recognition receptors (PRRs) expression based on HPV status. Virome analysis indicated that HPV infection correlates to the presence of Anellovirus (AV) and seven complete HPV genomes were assembled. Bacteriome results unveiled that vaginal community state types (CST) distribution was independent of HPV or AV status, although bacterial phyla distribution differed between groups. Furthermore, TLR3 and IFNαR2 levels were higher in the Lactobacillus no iners-dominated mucosa and we detected correlations among RIG-like receptors (RLR) associated genes and abundance of specific anaerobic bacteria. Collectively, our data show an intriguing connection between HPV and AV infections that could foster cervical cancer development. Besides that, TLR3 and IFNαR2 seem to create a protective milieu in healthy cervical mucosa (L. no iners-dominated), and RLRs, known to recognize viral RNA, were correlated to anaerobic bacteria suggesting that they might be related to dysbiosis.

Evidence of recurrent selection of mutations commonly found in SARS-CoV-2 variants of concern in viruses infecting immunocompromised patients

Chronically immunosuppressed patients infected with SARS-CoV-2 often experience prolonged virus shedding, and may pave the way to the emergence of mutations that render viral variants of concern (VOC) able to escape immune responses induced by natural infection or by vaccination. We report herein a SARS-CoV-2⁺ cancer patient from the beginning of the COVID-19 pandemic whose virus quasispecies across multiple timepoints carried several immune escape mutations found in more contemporary VOC, such as alpha, delta and omicron, that appeared to be selected for during infection. We hypothesize that immunosuppressed patients may represent the source of VOC seen throughout the COVID-19 pandemics.

New infections by SARS-CoV-2 variants of concern after natural infections and post-vaccination in Rio de Janeiro, Brazil

After a one-year rollout of the pandemic caused by SARS-CoV-2, the continuous dissemination of the virus has generated a number of variants with increased transmissibility and infectivity, called variants of concern (VOC), which now predominate worldwide. Concerns about the susceptibility of humans that have already been infected before or those already vaccinated to infection by VOC rise among scientists and clinicians. Herein, we assessed the prevalence of different VOC among recent infections at the Brazilian National Cancer Institute (Rio de Janeiro, Brazil). By using a Sanger-based sequencing approach targeting the viral S gene to identify VOC, we have analyzed 72 recent infections. The overall prevalence of VOC was 97%. Among the subjects analyzed, six had been vaccinated with the ChAdOx1-S/nCoV-19 (n = 4; one with two doses and three with one dose) or the CoronaVac (n = 2; both with 2 doses) vaccine, while five subjects represented reinfection cases, being two of them also part of the vaccinated group (each one with one vaccine type). All vaccinated and re-infected subjects carried VOC irrespective of the vaccine type taken, the number of doses taken, IgG titers or being previously infected during the first wave of the Brazilian pandemic. Importantly, all six vaccinees only had mild symptoms. We present here several examples of how natural infections or vaccination may not be fully capable of conferring sterilizing immunity against VOC.

Distinguishing SARS-CoV-2 bonafide re-infection from pre-existing minor variant reactivation

Different groups have recently reported events of SARS-CoV-2 reinfection, where patients had a sequence of positive-negative-positive RT-PCR tests. However, such events could be explained by different scenarios such as intermittent viral shedding, bonafide re-infection or multiple infection with alternating predominance of different viruses. Analysis of minor variants is an important tool to distinguish between these scenarios. Using ARTIC network PCR amplification and next-generation sequencing, we obtained SARS-CoV-2 sequences from two timepoints (with a time span of 102 days) of a patient followed at the Brazilian National Cancer Institute. Within-host variant analysis evidenced three single nucleotide variants (SNVs) at the consensus viral sequence in the second timepoint that were already present in the first timepoint as minor variants. Another five SNVs found in the second timepoint were not detected in the first sample sequenced, suggesting an additional infection by a yet another new virus. Our observation shed light into the existence of different viral populations that are present in dynamic frequencies and fluctuate during the course of SARS-CoV-2 infection. The detection of these variants in distinct disease events of an individual highlights a complex interplay between viral reactivation from a pre-existing minority variant and reinfection by a different virus.

SARS-CoV-2 genomic analyses in cancer patients reveal elevated intrahost genetic diversity

Numerous factors have been identified to influence susceptibility to Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) infection and disease severity. Cancer patients are more prone to clinically evolve to more severe COVID-19 conditions, but the determinants of such a more severe outcome remain largely unknown. We have determined the full-length SARS-CoV-2 genomic sequences of cancer patients and healthcare workers (non-cancer controls) by deep sequencing and investigated the within-host viral population of each infection, quantifying intrahost genetic diversity. Naso- and oropharyngeal SARS-CoV-2+ swabs from 57 cancer patients and 14 healthcare workers from the Brazilian National Cancer Institute were collected in April to May 2020. Complete genome amplification using ARTIC network V3 multiplex primers was performed followed by next-generation sequencing. Assemblies were conducted in Geneious R11, where consensus sequences were extracted and intrahost single nucleotide variants were identified. Maximum likelihood phylogenetic analysis was performed using PhyMLv.3.0 and lineages were classified using Pangolin and CoV-GLUE. Phylogenetic analysis showed that all but one strain belonged to clade B1.1. Four genetically linked mutations known as the globally dominant SARS-CoV-2 haplotype (C241T, C3037T, C14408T and A23403G) were found in the majority of consensus sequences. SNV signatures of previously characterized Brazilian genomes were also observed in most samples. Another 85 SNVs were found at a lower frequency (1.4%-19.7%) among the consensus sequences. Cancer patients displayed a significantly higher intrahost viral genetic diversity compared to healthcare workers. This difference was independent of SARS-CoV-2 Ct values obtained at the diagnostic tests, which did not differ between the two groups. The most common nucleotide changes of intrahost SNVs in both groups were consistent with APOBEC and ADAR activities. Intrahost genetic diversity in cancer patients was not associated with disease severity, use of corticosteroids, or use of antivirals, characteristics that could influence viral diversity. Moreover, the presence of metastasis, either in general or specifically in the lung, was not associated with intrahost diversity among cancer patients. Cancer patients carried significantly higher numbers of minor variants compared to non-cancer counterparts. Further studies on SARS-CoV-2 diversity in especially vulnerable patients will shed light onto the understanding of the basis of COVID-19 different outcomes in humans.

SARS-CoV-2 genomic and quasispecies analyses in cancer patients reveal relaxed intrahost virus evolution

Numerous factors have been identified to influence susceptibility to SARS-CoV-2 infection and disease severity. Cancer patients are more prone to clinically evolve to more severe COVID-19 conditions, but the determinants of such a more severe outcome remain largely unknown. We have determined the full-length SARS-CoV-2 genomic sequences of cancer patients and healthcare workers (HCW; non-cancer controls) by deep sequencing and investigated the within-host viral quasispecies of each infection, quantifying intrahost genetic diversity. Naso- and oropharyngeal SARS-CoV-2 + swabs from 57 cancer patients and 14 healthcare workers (HCW) from the Brazilian Cancer Institute were collected in April-May 2020. Complete genome amplification using ARTIC network V3 multiplex primers was performed followed by next-generation sequencing. Assemblies were conducted in Geneious R11, where consensus sequences were extracted and intrahost single nucleotide variants (iSNVs) were identified. Maximum likelihood phylogenetic analysis was performed using PhyMLv.3.0 and lineages were classified using Pangolin and CoV-GLUE. Phylogenetic analysis showed that all but one strain belonged to clade B1.1. Four genetically linked mutations known as the globally dominant SARS-CoV-2 haplotype (C241T, C3037T, C14408T and A23403G) were found in the majority of consensus sequences. SNV signatures of previously characterized Brazilian genomes were also observed in most samples. Another 85 SNVs were found at a lower frequency (1.4-19.7%). Cancer patients displayed a significantly higher intrahost viral genetic diversity compared to HCW (p = 0.009). Intrahost genetic diversity in cancer patients was independent of SARS-CoV-2 Ct values, and was not associated with disease severity, use of corticosteroids, or use of antivirals, characteristics that could influence viral diversity. Such a feature may explain, at least in part, the more adverse outcomes to which cancer/COVID-19 patients experience.

AUTHOR SUMMARY

Cancer patients are more prone to clinically evolve to more severe COVID-19 conditions, but the determinants of such a more severe outcome remain largely unknown. In this study, phylogenetic and variation analysis of SARS-CoV-2 genomes from cancer patients and non-cancer healthcare workers at the Brazilian National Cancer Institute were characterized by deep sequencing. Viral genomes showed signatures characteristic of Brazilian viruses, consistent with the hypothesis of local, community transmission rather than virus importation from abroad. Despite most genomes in patients and healthcare workers belonging to the same lineage, intrahost variability was higher in cancer patients when compared to non-cancer counterparts. The intrahost genomic diversity analysis presented in our study highlights the relaxed evolution of SARS-CoV-2 in a vulnerable population of cancer patients. The high number of minor variations can result in the selection of immune escape variants, resistance to potential drugs, and/or increased pathogenicity. The impact of this higher intrahost variability over time warrants further investigation.

Composite Analysis of the Virome and Bacteriome of HIV/HPV Co-Infected Women Reveals Proxies for Immunodeficiency

The human cervical microbiome is complex, and its role in health and disease has just begun to be elucidated. In this study, 57 cervical swab samples from 19 HIV/HPV co-infected women were analyzed for both virome and bacteriome composition. Virome analysis focused on circular DNA viruses through rolling circle amplification followed by next-generation sequencing (NGS). Data were assigned to virus families and genera, and HPV types were identified. NGS data of bacterial 16S from a subset of 24 samples were assigned to operational taxonomic units and classified according to vaginal microbiome community state types (CSTs). Four viral families were found: Papillomaviridae, Anelloviridae, Genomoviridae, and Herpesviridae. Papillomavirus reads were more abundant in women with premalignant cervical lesions, which were also strongly associated with multiple (≥3) high-risk HPV infection. Anellovirus read abundance was negatively correlated with host CD4+ T-cell counts. The bacteriome revealed the presence of CST III and CST IV, and women with ≥1% frequency of genomovirus or herpesvirus reads displayed an increased risk of carrying CST IV. By characterizing the composition of the cervical circular DNA viruses and the bacteriome of HIV/HPV co-infected women, we identified putative interactions between these two microorganism communities and their associations with patients’ clinical characteristics, notably immunodeficiency status.

Estimating HIV-1 Genetic Diversity in Brazil Through Next-Generation Sequencing

Approximately 36.7 million people were living with the human immunodeficiency virus (HIV) at the end of 2016 according to UNAIDS, representing a global prevalence rate of 0.8%. In Brazil, an HIV prevalence of 0.24% has been estimated, which represents approximately 830,000 individuals living with the virus. As a touristic and commercial hub in Latin America, Brazil harbors an elevated HIV genetic variability, further contributed by the selective pressure exerted by the host immune system and by antiretroviral treatment. Through the progress of the next-generation sequencing (NGS) techniques, it has been possible to expand the study of HIV genetic diversity, evolutionary, and epidemic processes, allowing the generation of HIV complete or near full-length genomes (NFLG) and improving the characterization of intra- and interhost diversity of viral populations. Greater sensitivity in the detection of viral recombinant forms represents one of the major improvements associated with this development. It is possible to identify unique or circulating recombinant forms using the near full-length viral genomes with increasing accuracy. It also permits the characterization of multiple viral infections within individual hosts. Previous Brazilian studies using NGS to analyze HIV diversity were able to identify several distinct unique and circulating recombinant forms and evidenced dual infections. These data unveiled unprecedented high rates of viral recombination and highlighted that novel recombinants are continually arising in the Brazilian epidemic. In the pooled analysis depicted in this report, HIV subtypes have been determined from HIV-positive patients in five states of Brazil with some of the highest HIV prevalence, three in the Southeast (Rio de Janeiro, São Paulo, and Minas Gerais), one in the Northeast (Pernambuco) and one in the South (Rio Grande do Sul). Combined data analysis showed a significant prevalence of recombinant forms (29%; 101/350), and a similar 26% when only NFLGs were considered. Moreover, the analysis was able to evidence the occurrence of multiple infections in some individuals. Our data highlight the great HIV genetic diversity found in Brazil and unveils a more accurate scenario of the HIV evolutionary dynamics in the region.

Complex virome in feces from Amerindian children in isolated Amazonian villages

The number of viruses circulating in small isolated human populations may be reduced by viral extinctions and rare introductions. Here we used viral metagenomics to characterize the eukaryotic virome in feces from healthy children from a large urban center and from three Amerindian villages with minimal outside contact. Numerous human enteric viruses, mainly from the Picornaviridae and Caliciviridae families, were sequenced from each of the sites. Multiple children from the same villages shed closely related viruses reflecting frequent transmission clusters. Feces of isolated villagers also contained multiple viral genomes of unknown cellular origin from the Picornavirales order and CRESS-DNA group and higher levels of nematode and protozoan DNA. Despite cultural and geographic isolation, the diversity of enteric human viruses was therefore not reduced in these Amazonian villages. Frequent viral introductions and/or increased susceptibility to enteric infections may account for the complex fecal virome of Amerindian children in isolated villages.

Assessment of the gorilla gut virome in association with natural simian immunodeficiency virus infection

Background

Simian immunodeficiency viruses (SIVs) of chimpanzees and gorillas from Central Africa crossed the species barrier at least four times giving rise to human immunodeficiency virus type 1 (HIV-1) groups M, N, O and P. The paradigm of non-pathogenic lentiviral infections has been challenged by observations of naturally infected chimpanzees with SIVcpz associated with a negative impact on their life span and reproduction, CD4⁺ T-lymphocyte loss and lymphoid tissue destruction. With the advent and dissemination of new generation sequencing technologies, novel promising markers of immune deficiency have been explored in human and nonhuman primate species, showing changes in the microbiome (dysbiosis) that might be associated with pathogenic conditions. The aim of the present study was to identify and compare enteric viromes of SIVgor-infected and uninfected gorillas using noninvasive sampling and ultradeep sequencing, and to assess the association of virome composition with potential SIVgor pathogenesis in their natural hosts.

Results

We analyzed both RNA and DNA virus libraries of 23 fecal samples from 11 SIVgor-infected (two samples from one animal) and 11 uninfected western lowland gorillas from Campo-Ma’an National Park (CP), in southwestern Cameroon. Three bacteriophage families (Siphoviridae, Myoviridae and Podoviridae) represented 67.5 and 68% of the total annotated reads in SIVgor-infected and uninfected individuals, respectively. Conversely, mammalian viral families, such as Herpesviridae and Reoviridae, previously associated with gut- and several mammalian diseases were significantly more abundant (p < 0.003) in the SIVgor-infected group. In the present study, we analyzed, for the first time, the enteric virome of gorillas and their association with SIVgor status. This also provided the first evidence of association of specific mammalian viral families and SIVgor in a putative dysbiosis context.

Conclusions

Our results suggested that viromes might be potentially used as markers of lentiviral disease progression in wild gorilla populations. The diverse mammalian viral families, herein described in SIVgor-infected gorillas, may play a pivotal role in a disease progression still unclear in these animals but already well characterized in pathogenic lentiviral infections in other organisms. Larger sample sets should be further explored to reduce intrinsic sampling variation.

Characterization of HIV-1 Near Full-Length Proviral Genome Quasispecies from Patients with Undetectable Viral Load Undergoing First-Line HAART Therapy

Increased access to highly active antiretroviral therapy (HAART) by human immunodeficiency virus postive (HIV⁺) individuals has become a reality worldwide. In Brazil, HAART currently reaches over half of HIV-infected subjects. In the context of a remarkable HIV-1 genetic variability, highly related variants, called quasispecies, are generated. HIV quasispecies generated during infection can influence virus persistence and pathogenicity, representing a challenge to treatment. However, the clinical relevance of minority quasispecies is still uncertain. In this study, we have determined the archived proviral sequences, viral subtype and drug resistance mutations from a cohort of HIV⁺ patients with undetectable viral load undergoing HAART as first-line therapy using next-generation sequencing for near full-length virus genome (NFLG) assembly. HIV-1 consensus sequences representing NFLG were obtained for eleven patients, while for another twelve varying genome coverage rates were obtained. Phylogenetic analysis showed the predominance of subtype B (83%; 19/23). Considering the minority variants, 18 patients carried archived virus harboring at least one mutation conferring antiretroviral resistance; for six patients, the mutations correlated with the current ARVs used. These data highlight the importance of monitoring HIV minority drug resistant variants and their clinical impact, to guide future regimen switches and improve HIV treatment success.

Analysis of the cervical microbiome and potential biomarkers from postpartum HIV-positive women displaying cervical intraepithelial lesions

The cervical microbiota composition and diversity of HIV-positive women in the postpartum period is unknown. Using a high-throughput bacterial 16S rRNA gene sequencing, we identified four community state types (CSTs). CST III (Lactobacillusdominant) and CST IV (IV-A, IV-B.1, IV-B.2; high-diversity) were found in 41% and 59% of samples, respectively. We did not find association of any CST to postpartum period (six or twelve months), HPV infection or cytology (normal or lesion). However, five bacterial genera were associated with cervical lesions (Gardnerella, Aerococcus, Schlegelella, Moryella and Bifidobacterium), with significant odds ratio (OR) of 40 (2.28–706) for the presence of Moryella and 3.5 (1.36–8.9) for Schlegelella. Longitudinal analysis of samples at postpartum that regressed (lesion to normal), progressed (normal to lesion) and maintained the cytology (lesion or normal) evidenced Gardnerella with a significantly higher abundance in regressing lesions. In the current study, we report the first data on the cervical microbiota of HIV-positive women in the postpartum period. Consistent with previous studies of HIV-negative cohorts, HIV-positive women present a stable cervical microbiota of high-diversity in the postpartum period. Our results highlight that specific microbiota species may serve as sensors for changes in the cervical microenvironment associated with cervical lesions.

Characterisation of complete high- and low-risk human papillomavirus genomes isolated from cervical specimens in southern Brazil

The classification of human papillomavirus (HPV) intratypic lineages by complete genome sequencing is a determinant in understanding biological differences in association with this disease. In this work, we have characterised complete HPV genomes from southern Brazil. Fifteen cervicovaginal Pap smear negative samples previously categorised as HPV-positive were sequenced using ultradeep sequencing, and 18 complete genomes from 13 different HPV types were assembled. Phylogenetic and genetic distance analyses were performed to classify the HPV genomes into lineages and sublineages. This is the first report describing the distribution of HPV intratype lineages of high and low oncogenic risk in asymptomatic women from southern Brazil.

Bonafide, type-specific human papillomavirus persistence among HIV-positive pregnant women: predictive value for cytological abnormalities, a longitudinal cohort study

This study investigated the rate of human papillomavirus (HPV) persistence, associated risk factors, and predictors of cytological alteration outcomes in a cohort of human immunodeficiency virus-infected pregnant women over an 18-month period. HPV was typed through L1 gene sequencing in cervical smears collected during gestation and at 12 months after delivery. Outcomes were defined as nonpersistence (clearance of the HPV in the 2nd sample), re-infection (detection of different types of HPV in the 2 samples), and type-specific HPV persistence (the same HPV type found in both samples). An unfavourable cytological outcome was considered when the second exam showed progression to squamous intraepithelial lesion or high squamous intraepithelial lesion. Ninety patients were studied. HPV DNA persistence occurred in 50% of the cases composed of type-specific persistence (30%) or re-infection (20%). A low CD4⁺T-cell count at entry was a risk factor for type-specific, re-infection, or HPV DNA persistence. The odds ratio (OR) was almost three times higher in the type-specific group when compared with the re-infection group (OR = 2.8; 95% confidence interval: 0.43-22.79). Our findings show that bonafide (type-specific) HPV persistence is a stronger predictor for the development of cytological abnormalities, highlighting the need for HPV typing as opposed to HPV DNA testing in the clinical setting.

Identification of novel human papillomavirus lineages and sublineages in HIV/HPV-coinfected pregnant women by next-generation sequencing

Infection by human papillomavirus (HPV) is a necessary condition for development of cervical cancer, and has also been associated with malignancies of other body anatomical sites. Specific HPV types have been associated with premalignant lesions and invasive carcinoma, but mounting evidence suggests that within-type lineages and sublineages also display distinct biological characteristics associated with persistent infections and evolution to cervical cancer. In the present study, we have assessed HPV multiple infection and variation from a cohort of highly susceptible, HIV(+) pregnant women using next-generation sequencing and an in-house pipeline for HPV full-length genome assembly. Seventy-two consensus sequences representing complete or near-complete (>97%) HPV genomes were assembled, spanning 28 different types. Genetic distance and phylogenetic analyses allowed us to propose the classification of novel HPV lineages and sublineages across nine HPV types, including two high-risk types. HPV diversity may be a hallmark of immunosuppressed patients upon HIV infection and AIDS progression.

The effect of human leukocyte antigen G alleles on human papillomavirus infection and persistence in a cohort of HIV-positive pregnant women from Brazil

Patients with compromised immune systems have more severe intraepithelial lesions and more rapid disease progression, in addition to increased risk for cervical cancer. Persistent infection by the human papillomavirus (HPV) is a necessary step in that process. By inducing expression of inhibitory ligands of natural killer cells, like HLA-G, HPV avoids the elimination of infected cells. Recent studies have investigated polymorphisms in HLA-G that may be associated with susceptibility to HPV infection and persistence. One hundred-forty HIV(+) pregnant women from Brazil had a DNA fragment comprising HLA-G exons 2-4 PCR-amplified, cloned, sequenced and analyzed for allele determination. Altogether, 22 alleles comprising 52 different genotypes were found. Four novel HLA-G alleles were characterized. We have not observed association of specific HLA-G alleles with HPV infection, but found a protective effect of the G:01:01:02 allele against the occurrence of intraepithelial lesions. In addition to describing new HLA-G alleles and defining new reference sequences, our data provide a better understanding of the impact of HLA-G alleles on HPV-related disease.

HIV/HPV co-infection during pregnancy in southeastern Brazil: prevalence, HPV types, cytological abnormalities and risk factors

OBJECTIVE

HIV(+) pregnant women are at a higher risk of HPV infection and development of cervical cancer. Our objectives were to assess the prevalence and HPV types in HIV(+) pregnant women and to identify risk factors for HPV infection and cytological abnormalities.

METHODS

Cervicovaginal smears were collected during pregnancy from 140 women. Partial HPV L1 gene and the exon 4 of the human TP53 gene (containing codon 72) were PCR-amplified and sequenced. Amplified products indicating multiple HPV infection were further cloned and sequenced. The association of demographic, obstetric and HIV-related clinical variables with HPV infection and cervical lesions was tested by univariate analyses, and significant factors were subsequently tested by logistic regression multivariate analysis.

RESULTS

HPV DNA tested positive for 118 patients and HPV types were identified in 104 samples. Twenty-eight different types were found, HPV-16 and HPV-58 being the most prevalent. High-risk types were present in 79.8% of samples and multiple infections in 16.3%. Abnormal cervical smears were found in 44 patients (31.4%). Absolute CD4(+) T-cell counts below 350 were associated with HPV infection. Younger age was associated with cervical abnormalities and higher CD4(+) T-cell count was an apparent protective factor.

CONCLUSIONS

We found a high prevalence of HPV infection and high-risk types in this cohort. Our results highlighted the relevance of immune system integrity rather than TP53 variants for protecting this highly vulnerable population to HPV infection and carcinogenesis.