Multiple sexual partners and vaginal microecological disorder are associated with HPV infection and cervical carcinoma development.

There is an indirect link between multiple sexual partners (MSP) and cervical intraepithelial neoplasia (CIN) or even cervical cancer (CC). MSP may also lead to bacterial vaginosis (BV). The relationship among MSP, BV, human papillomavirus (HPV) infection and CIN/CC development in Chinese women remains unclear. The present study was designed to clarify their association. The study retrospectively analyzed 549 female patients who had visited a physical examination center. The MSP information was acquired, and vaginal microecology, HPV and cervical conization pathology (CCP) tests were performed when necessary. MSP status was distinct among patients with different levels of BV severity. In addition, as the severity of BV progressed, the HPV-positive ratio increased. Meanwhile, MSP was significantly associated with a positive HPV outcome, including HPV 16, HPV 18 and other high-risk HPV infections. The MSP group had a significantly higher percentage of positive CCP outcomes (particularly cases with CIN-II and CIN-III). Similarly, higher BV severity meant more severe CIN/CC progression. A logistic regression model based on age, MSP status and the Nugent score level was used in order to predict the CCP outcome. Furthermore, a receiver operating characteristic curve analysis resulted in an area under the curve of 0.834. In conclusion, the combination of MSP information and BV examination may provide a rapid, economic and accurate prediction of CIN/CC. Health education on sexual behavior and timely detection/treatment of BV should be conducted to reduce the risk of CC. Copyright: © Huang et al.

Introduction

Cervical cancer (CC) is a serious public health problem in China, with ~102,000 new cases estimated according to the National Central Cancer Registry in 2013 (1). In addition, CC is the second most common malignant tumor in women worldwide (2,3). It is widely known that human papillomavirus (HPV) infection (4), family history of CC (5) and HIV infection (6) contribute to the development of CC. HPV infection is one of the main causes of CC. Nearly all cases of CC are causally associated with HPV infection (7,8), particularly high-risk types of HPV, including HPV 16 and 18 (9–11). Simultaneously, sexual contact is the dominant risk factor for HPV infection, particularly in those with multiple sexual partners (MSP) (12). Thus, there is an indirect link between MSP and cervical intraepithelial neoplasia (CIN)/CC.

CIN is a precancerous lesion characterized by abnormal cell growth in the lining of the cervix. CIN can be divided into three stages, according to the degree of the lesion: CIN-I, CIN-II and CIN-III. CIN has the potential to progress to CC if left untreated; this usually takes >10 years (13). Moreover, having MSP has been reported as an independent risk factor of CC after adjusting for HPV infection (14,15). Furthermore, the link between MSP and CIN/CC is considered to be partially associated with geographical region, as a meta-analysis of epidemiological studies revealed (16). Another consequence of MSP is an imbalance in the vaginal microecology, which theoretically directly leads to bacterial vaginosis (BV). A previous study reported a positive association between BV and cervical HPV infection (17). Dysregulated vaginal microbiota can enhance the risk of transient and persistent HPV infections (18), thus BV and HPV may jointly increase the risk of CIN. It was reported that BV is a cofactor of HPV and the development of CIN (19), and a significantly increased risk of CIN has been detected when BV is present (20,21). However, controversial results also exist, for example, a study in 1995 reported that, in women with dyskaryotic cervical smears, the prevalence of bacterial vaginosis did not seem to be increased, and bacterial vaginosis did not influence the histologic changes. In addition, in women with low grade squamous intraepithelial lesion, aerobic vaginitis is very common but is not an indicator of HPV infection (22,23). Since previous studies have only included a small number of samples (17–21), and the relationship among MSP, BV, HPV and CIN/CC development in Chinese women remains unclear, the present study was designed to clarify their association. Notably, the current study indicated that when information regarding MSP has been obtained, it is important to assess the severity of BV to better predict HPV and CIN/CC development in Chinese women.

Materials and methods

Patients

The present study retrospectively analyzed 549 female patients who had visited the physical examination center, outpatient and inpatient departments of Fujian Provincial Hospital (Fuzhou, China) from January 2017 to December 2019. The inclusion criteria were as follows: i) Aged between 20 and 55 years old, with regular menstruation and sexual activity history; ii) suspicion of high risk of HPV infection and BV infection; iii) 1 week before diagnosis, there was no history of sexual activity, the use of vaginal plugs, vaginal irrigation or other treatments. The exclusion criteria were: i) Menopause; ii) a history of total hysterectomy due to other diseases; iii) being pregnant; or iv) presenting with severe malignant tumors, AIDS and other immunodeficiencies.

The MSP information was acquired using a questionnaire; accordingly, recruited patients were divided into two groups: No MSP and MSP. The MSP group consisted of patients who had more than one male sexual partner in the last 6 months. The no MSP group comprised: i) Patients with a fixed male sexual partner, ii) patients who lived alone and had not had sex in the last 6 months and iii) those who lived alone but occasionally had sex with a fixed male.

Vaginal microecology test

The severity of BV was determined using a gram-stained slide and was detected with the Cro vaginal microbiological evaluation system (ShanDong Shida Si Biological Industry Co., Ltd.). The results were hierarchically represented by the Nugent Gram stain scoring system (24–27). A Nugent score of 7 to 10 was defined as severe BV, a score of 4 to 6 was defined as mild BV, and a score of 1 to 3 was defined as negative BV.

Cobas high-risk HPV test

The fully automated Cobas ×480 nucleic acid extraction instrument and Cobas z480 PCR analyzer instrument [both from Roche Diagnostics (Shanghai) Co., Ltd.] were utilized according to the manufacturer's protocol to detect HPV types 16, 18 and 12 other high-risk subtypes (types 31, 33, 35, 39, 45, 51, 52, 56, 58, 59, 66 and 68).

Cervical conization pathology (CCP)

Colposcopy was ordered if a patient had abnormal thin-layer cytology test results or was highly possibly with CIN/CC development. The colposcope magnified the visual field, and acetic acid was used to swab the vaginal surface to discolor abnormal tissues and easily identify lesions. When abnormal areas were identified, colposcopy was followed by cervical conization to determine the final diagnosis (24,25). Briefly, a cone-shaped tissue sample around and below the zone of the cervical squamocolumnar junction was resected. The resection included normal tissue ranging from 3 to 5 mm outside the abnormal tissue, and the cone depth ranged from 15 to 25 mm. Cone cutting depth was determined by the site of the lesion. The tissue samples were then immediately sent for pathological examination. According to the CCP results, these patients were subdivided into CIN-I, CIN-II, CIN-III and CC groups. The remaining patients who had normal outcomes and did not receive CCP were divided into the No-Need-to-Test group.

Statistical analysis

All categorical data were expressed in frequencies and percentages; continuous data were expressed as the mean ± standard deviation. SPSS software (version 25.0; IBM Corp.) was used for statistical analysis and data visualization. The χ2 test was used to compare frequencies in different groups, and one-way ANOVA was used for comparison of continuous data. The logistic regression model was established to predict the outcomes of different levels of CCP outcomes (No-need-to-test, CIN-I, CIN-II, CIN-III and CC). In addition, receiver operating characteristic (ROC) curves were generated for the prediction of CCP-positive results, and the area under the curve (AUC), as well as recommended threshold values in distinguishing different CCP outcomes, were calculated. P<0.05 was considered to indicate a statistically significant difference.

Results

Age and MSP distribution

The demographic characteristics, including age and MSP status, are presented in Table I. A total of 549 women were included in the present study, with an average age of 36.01±8.07 years. Of these, 309 women tested negative for BV, and among women that tested positive for BV, 141 were in the BV mild group and 99 were in the BV severe group, according to the Nugent score. The three groups displayed a similar age distribution with no significant differences across BV groups, whereas MSP status was distinct among patients with different levels of BV severity. For women who did not have MSP, the number in the BV-negative group was 298, whereas the number in the BV mild group was 116, and in the BV severe group was 75. For those with MSP, 11 women tested negative for BV, 25 women had mild BV and 24 women had severe BV. The percentage of patients with mild and severe BV were significantly increased in the MSP group (P<0.001).

Table I.

Age and MSP status of enrolled patients.

Index	Total (n=549)	BV negative (n=309)	BV mild (n=141)	BV severe (n=99)	F or χ2	P-value
Age	36.01±8.07	36.20±7.92	35.44±7.92	36.24±8.75	0.476	0.622
MSP, n (%)					41.964	<0.001
No	489 (89.1%)	298 (96.4%)	116 (82.3%)	75 (75.8%)
Yes	60 (10.9)	11 (3.6%)	25 (17.7%)	24 (24.2%)

BV, bacterial vaginosis; MSP, multiple sexual partners.

Association between MSP, BV severity and HPV infection outcomes

Next, the present study investigated the relationship between BV severity and HPV infection. As BV severity progressed, the HPV-positive ratio increased (Table II). This association was consistent across all HPV groups, including total HPV, HPV 16, HPV 18 and other high-risk HPV infections. Furthermore, MSP was significantly associated with a HPV-positive outcome; this association was observed for HPV 16, HPV 18 and other high-risk HPV infections (Table III).

Table II.

Association between BV severity and HPV infection.

HPV infection	BV negative (%)	BV mild (%)	BV severe (%)	χ2	P-value
Positive				60.833	<0.001
No	174 (56.3)	45 (31.9)	15 (15.2)
Yes	135 (43.7)	96 (68.1)	84 (84.8)
HPV 16				39.000	<0.001
No	289 (93.5)	111 (78.7)	70 (70.7)
Yes	20 (6.5)	30 (21.3)	29 (29.3)
HPV 18				20.946	<0.001
No	298 (96.4)	128 (90.8)	82 (82.8)
Yes	11 (3.6)	13 (9.2)	17 (17.2)
Other HPV types				41.225	<0.001
No	204 (66.0)	70 (49.6)	30 (30.3)
Yes	105 (34.0)	71 (50.4)	69 (69.7)

BV, bacterial vaginosis; HPV, human papillomavirus.

Table III.

Association between MSP and HPV infection.

HPV infection	No MSP (n=489) (%)	MSP (n=60) (%)	χ2	P-value
Positive			23.630	<0.001
No	226 (46.2)	8 (13.3)
Yes	263 (53.8)	52 (86.7)
HPV 16			31.348	<0.001
No	433 (88.5)	37 (61.7)
Yes	56 (11.5)	23 (38.3)
HPV 18			24.536	<0.001
No	462 (94.5)	46 (76.7)
Yes	27 (5.5)	14 (23.3)
Other HPV types			15.320	<0.001
No	285 (58.3)	19 (31.7)
Yes	204 (41.7)	41 (68.3)

HPV, human papillomavirus; MSP, multiple sexual partners.

Association between MSP, BV severity and CIN/CC development

Analysis of the association between MSP and CCP outcome demonstrated that the MSP group had a significantly higher percentage of positive CCP outcomes (Table IV). Notably, the MSP group had a higher proportion of CIN-II and CIN-III cases. Similarly, Nugent score classification was significantly associated with the CCP outcome (Table V); higher BV severity was associated with a more severe CIN/CC progression. To understand the association between MSP, BV and CCP outcome, six groups were established according to BV severity and MSP status: No-MSP-BV-negative, No-MSP-BV-mild, No-MSP-BV-severe, MSP-BV-negative, MSP-BV-mild and MSP-BV-severe. Positive CCP outcomes were significantly associated with MSP and BV; the percentage of positive CCP results was successively increased from No-MSP-BV-negative, to No-MSP-BV-mild, MSP-BV-negative, No-MSP-BV-severe, MSP-BV-mild and MSP-BV-severe (Table VI). These univariate results suggested that the combination of MSP and the Nugent score level may be helpful to predict the final CCP results.

Table IV.

Association between MSP and CCP outcome.

CCP outcome	No MSP (n=489) (%)	MSP (n=60) (%)	χ2	P-value
Result			38.683	<0.001
Negative	445 (91.0)	38 (63.3)
Positive	44 (9.0)	22 (36.7)
Level			39.991	<0.001
No-need-to-test	443 (90.6)	38 (63.3)
CIN-I	2 (0.4)	0 (0.0)
CIN-II	24 (4.9)	12 (20.0)
CIN-III	19 (3.9)	10 (16.7)
CC	1 (0.2)	0 (0.0)

CCP, cervical conization pathology; CC, cervical cancer; CIN, cervical intraepithelial neoplasia; MSP, multiple sexual partners.

Table V.

Association between Nugent score level and CCP outcome.

CCP outcome	BV negative (%)	BV mild (%)	BV severe (%)	χ2	P-value
Result				74.886	<0.001
Positive	297 (96.1)	123 (87.2)	63 (63.6)
Negative	12 (3.9)	18 (12.8)	36 (36.4)
Level				89.104	<0.001
No-need-to-test	297 (96.1)	123 (87.2)	61 (61.6)
CIN-I	0 (0.0)	0 (0.0)	2 (2.0)
CIN-II	8 (2.6)	11 (7.8)	17 (17.2)
CIN-III	4 (1.3)	7 (5.0)	18 (18.2)
CC	0 (0.0)	0 (0.0)	1 (1.0)

CCP, cervical conization pathology; BV, bacterial vaginosis; CC, cervical cancer; CIN, cervical intraepithelial neoplasia.

Table VI.

Frequencies of positive cervical conization pathology outcome in six subgroups.

Group	Negative (%)	Positive (%)	χ2	P-value
No-MSP-BV-negative	287 (96.3	11 (3.7)	95.301	<0.001
No-MSP-BV-mild	106 (91.4)	10 (8.6)
No-MSP-BV-severe	52 (69.3)	23 (30.7)
MSP-BV-negative	10 (90.9)	1 (9.1)
MSP-BV-mild	17 (68.0)	8 (32.0)
MSP-BV-severe	11 (45.8)	13 (54.2)

BV, bacterial vaginosis; MSP, multiple sexual partners.

Next, we used for the logistic two-ends regression model (Table VII) based on age (continuous variable), MSP status (0, no; 1, yes) and the Nugent score level (1, negative; 2, mild; 3, severe) in order to predict the CCP outcome (negative or positive). The calculated value of the logistic regression model using three variables (age, MSP and Nugent score, namely AMN) revealed an AUC of 0.834 in the ROC curve. This predictor had a sensitivity level of 0.803 and specificity level of 0.702 at the cutoff value of 0.0818. When only two indexes (MSP and Nugent score, namely MN) were applied, an AUC of 0.789 was revealed. Meanwhile, the AUC of the single index Nugent score was 0.768, which highlights the informative role of BV examination (Fig. 1). Together, the severity of BV infection may be used to predict CIN/CC development, with satisfactory sensitivity and accuracy, and vaginal microecology may strongly impact infection and cervical carcinogenesis.

Table VII.

Variables in the logistic regression model of cervical conization pathology outcomes.

Variable	B	Wald	P-value	Odd ratio	95% CI of odd ratio
Age	0.095	22.282	<0.0001	1.099	1.057–1.143
MSP or not	1.708	20.551	<0.0001	5.518	2.637–11.548
BV severity levels		44.181	<0.0001
Mild vs. negative	1.070	6.779	0.0090	2.915	1.303–6.523
Severe vs. negative	2.466	41.582	<0.0001	11.780	5.566–24.930

Figure 1.

Receiver operating characteristic curve for efficient analysis of risk factors in CIN/CC prediction. Two variables, AMN and MN, were used to predict the outcomes of cervical conization pathology. AMN, age, MSP and Nugent score; AUC, area under the curve; MN, MSP and Nugent score; MSP: Multiple sexual partners.

Discussion

Sexually transmitted infections, including HPV, have been reported to be associated with the etiology of cervical carcinogenesis. Although HPV infection can be eliminated within 2 years without any complications, 1% of HPV-positive women develop CC (in the United States of America) (4). Thus, other potential risk factors need to be investigated to better predict CC. It is known that the normal vaginal microecology plays an indispensable role in preventing HPV infection and accelerating HPV clearance (26,27). Abnormal vaginal microecology may be a synergistic factor of persistent HPV infection (28), which ultimately increases the risk of CIN/CC. A study with a twin cohort reported that increased vaginal microbiota diversity was closely related to HPV infection (29). HPV is a sexually transmitted infection, and increased multiple-partner sexual contact may enhance HPV transmission and CC risk (30). A collaborative study that collected data from >15,000 women with CIN-III/CC confirmed the relationship between major indicators of sexual behavior and the risk of CC (31). Given that regional disparity can influence the relationship between MSP, BV severity, HPV and CIN/CC (16), it is worth exploring the association between these variables in Western and Eastern countries. A Spanish study reported that MSP increased the risk of HPV infection (12). In addition, a recent study revealed that marital status, age at first sexual intercourse, number of lifetime sexual partners and condom use were related to the presence of cervical pathology and the presence of oral HPV (32). Still, there have not been sufficient studies in other countries, particularly Eastern countries, investigating these associations. In the current study, a cohort of the Chinese population was used to confirm that MSP and BV were associated with CIN/CC in Chinese women.

MSP has been associated with an increased risk of CIN/CC. One of the key causes of CIN/CC development is long-term infection with HPV (33). In addition, vaginal microbiota impact the acquisition, persistence and clearance of HPV, and hence affect the risk of developing CC (34). MSP status contributes to the abnormality of vaginal microecology, and this abnormal intravaginal microbial environment can be the basis for CIN/CC development. In addition, MSP may result in the introduction of other sexually transmitted pathogens, such as HIV, which is known to increase the risk of CIN/CC (35).

In the present study, patients potentially infected with HPV and BV were enrolled to clarify the relationships between MSP, BV and CIN/CC development. MSP and different levels of BV severity were found to be associated with HPV infection and CIN/CC. Through logistic regression analysis of CCP outcomes, MSP status and Nugent score were revealed to be risk factors for CIN/CC prediction. Using two-parameter (MN) or three-parameter (AMN) prediction models, results obtained showed a satisfactory accuracy and sensitivity for the prediction of CIN/CC.

Overall, the findings of the current study are consistent with previously published articles. According to a study of 211 cervical samples, BV and inflammatory response were independently associated with the severity of cervical neoplasia in HPV-positive women (36). A meta-analysis also provided the same association between BV and uterine cervical HPV infection (17). The mechanisms underlying this association may include the production of carcinogenic nitrosamines and an altered inflammatory cytokine profile based on higher vaginal pH in BV-infected patients with persistent HPV (37,38). However, more mechanisms need to be explored. It has been suggested that HPV examination provides 60–70% higher protection against invasive CC and it is usually regarded as a more direct way of early warning of CC than cytological methods (39). The present work highlights the role of vaginal microecological testing in the monitoring of CIN/CC. In addition, it was confirmed that BV might be involved in the pathological process of CIN/CC. In the BV mild stage, the percentage of positive CCP cases sharply increased in comparison with BV negative. Therefore, even for patients with mild BV, timely detection and treatment should be conducted to halt the development of CIN/CC.

There are some limitations in the present study. Firstly, this cross-sectional study lacks longitudinal data and the causality of the MSP-BV-HPV-CIN axis (as well as other potential causal chains) remains unclear. Secondly, the single-index MSP used in this study is not informative enough, which only concerns the total number of sexual partners. Other sexual behaviors, such as age at first intercourse and the definite number of sexual partners, were not included in this study. Moreover, the MSP index only focused on sexual activity with men, but not women. Finally, HPV 16 and HPV 18 were primarily used as the first-choice risk factors, but other high-risk HPV strains, such as HPV 31, HPV 33 and HPV 35 that are implicated in 99% of CC cases have not been analyzed (9,40).

In conclusion, MSP and BV severity may be the main risk factors for HPV and CC development among Chinese women. A combination of MSP information and BV examination could provide a rapid, economic and accurate prediction of CIN/CC. Therefore, appropriate sexual health education, and prompt detection and treatment of BV could reduce the risk of CC.