Mariam El-Zein1, Agnihotram V Ramanakumar1, Paulo Naud2, Cecilia M Roteli-Martins3, Newton S de Carvalho4, Paola Colares de Borba5, Julio C Teixeira6, Anna-Barbara Moscicki7, Diane M Harper8, Stephen K Tyring9, Brian Ramjattan10, Gary Dubin11, Eduardo L Franco1. 1. From the Division of Cancer Epidemiology, McGill University, Montreal, Quebec, Canada. 2. Federal University of Rio Grande do Sul-UFRGS/HCPA, Hospital de Clínicas de Porto Alegre, Porto Alegre. 3. Leonor Mendes de Barros Hospital-Secretaria da Saude de São Paulo, São Paulo. 4. Department of Gynecology and Obstetrics, Gynecology and Obstetrics Infectious Diseases Sector, Federal University of Parana, Curitiba, Parana. 5. Instituto de Prevenção do Câncer, Fortaleza. 6. Department of Gynecology and Obstetrics, Oncology Division-CAISM, State University of Campinas, Campinas, Brazil. 7. University of California, Division of Adolescent Medicine, San Francisco, CA. 8. Department of Obstetrics and Gynecology and Community and Family Medicine, Geisel Medical School at Dartmouth, Norris Cotton Cancer Center, Hanover, NH. 9. Departments of Microbiology/Molecular Genetics, Dermatology & Internal. 10. First Line Medical Services Ltd., St. John's, Newfoundland and Labrador, Canada. 11. Takeda Pharmaceuticals, Zurich, Switzerland.
Abstract
BACKGROUND: Global variation in human papillomavirus (HPV) prevalence and persistence may be explained by differences in risk factors, such as sexual activity, oral contraceptive use, and behavioral factors. We evaluated determinants of acquisition and clearance of HPV infection among young women previously unexposed to HPV. METHODS:Five hundred thirty-four women aged 15 to 25 years who were cytology and HPV DNA negative, and seronegative for anti-HPV-16/18 antibodies, were recruited (July 2000-September 2001) from study centers in Brazil, the United States, and Canada (NCT00689741/NCT00120848). They were followed up for 76 months. Cervical samples were HPV genotyped via polymerase chain reaction. We used multivariable (forward stepwise, P = 0.15) Cox proportional hazards regression to estimate rate ratios (RR) and 95% confidence intervals (CI), separately according to length of follow-up time. RESULTS: On short-term follow-up (0-27 months), 257 (48%; 8535.80 person-months; incidence rate = 30.11; 95% CI, 26.64-34.02) incident HPV infections were detected. Marital status, lifetime number of sex partners, history of any sexually transmitted disease, and occasional use of oral contraceptives were strongly associated with acquisition of any HPV. Having 2 or more lifetime sex partners (RR, 2.03; 95% CI, 1.37-3.02) and a history of any sexually transmitted disease (RR, 1.98; 95% CI, 1.19-3.29) were the most important determinants of high-risk HPV (hrHPV) incidence. During the entire follow-up (0-76 months), an increased hrHPV clearance was found among women in North America (RR, 1.38; 95% CI, 1.08-1.78) and black women (RR, 1.64; 95% CI, 1.04-2.60). Greater number of lifetime partners was associated with reduced clearance rates for any HPV (RR, 0.65; 95% CI, 0.43-0.98). CONCLUSIONS: We identified variation in risk of HPV acquisition and clearance among women unexposed to HPV at baseline.
RCT Entities:
BACKGROUND: Global variation in human papillomavirus (HPV) prevalence and persistence may be explained by differences in risk factors, such as sexual activity, oral contraceptive use, and behavioral factors. We evaluated determinants of acquisition and clearance of HPV infection among young women previously unexposed to HPV. METHODS: Five hundred thirty-four women aged 15 to 25 years who were cytology and HPV DNA negative, and seronegative for anti-HPV-16/18 antibodies, were recruited (July 2000-September 2001) from study centers in Brazil, the United States, and Canada (NCT00689741/NCT00120848). They were followed up for 76 months. Cervical samples were HPV genotyped via polymerase chain reaction. We used multivariable (forward stepwise, P = 0.15) Cox proportional hazards regression to estimate rate ratios (RR) and 95% confidence intervals (CI), separately according to length of follow-up time. RESULTS: On short-term follow-up (0-27 months), 257 (48%; 8535.80 person-months; incidence rate = 30.11; 95% CI, 26.64-34.02) incident HPV infections were detected. Marital status, lifetime number of sex partners, history of any sexually transmitted disease, and occasional use of oral contraceptives were strongly associated with acquisition of any HPV. Having 2 or more lifetime sex partners (RR, 2.03; 95% CI, 1.37-3.02) and a history of any sexually transmitted disease (RR, 1.98; 95% CI, 1.19-3.29) were the most important determinants of high-risk HPV (hrHPV) incidence. During the entire follow-up (0-76 months), an increased hrHPV clearance was found among women in North America (RR, 1.38; 95% CI, 1.08-1.78) and black women (RR, 1.64; 95% CI, 1.04-2.60). Greater number of lifetime partners was associated with reduced clearance rates for any HPV (RR, 0.65; 95% CI, 0.43-0.98). CONCLUSIONS: We identified variation in risk of HPV acquisition and clearance among women unexposed to HPV at baseline.
Human papillomavirus (HPV) is one of the most common sexually transmitted infections worldwide. Although most infections resolve spontaneously within 3 years,[1] persistent infection with high-risk types can lead to oral and anogenital cancers via consecutive premalignant stages.[2] Low-risk genotypes (eg, HPVs 6 and 11), also known as non-oncogenic HPV, are associated with benign conditions, such as genital warts and respiratory papillomatosis.[3] Knowledge of determinants of HPV acquisition and persistence is essential in understanding its epidemiology and natural history to develop future effective interventions aimed at reducing HPV-related disease burden.[4]Epidemiologic studies have identified risk factors associated with HPV infection. There is clear evidence that some patterns of sexual behavior, such as age at sexual debut (ie, commencement of sexual activity at an early age), multiple sexual partners, and partners' prior sexual history are directly associated with an increased risk of genital HPV acquisition.[4,5] Age, smoking, co-infection with other sexually transmitted agents, long-term oral contraceptive use and parity are also key cofactors to cervical carcinogenesis in high-risk HPV (hrHPV) positive women.[6]Most research on HPV-induced carcinogenesis included women over the age of 30 years, focusing on those attending cervical cancer screening. HPV clearance of prevalent (ie, existing cases) or mixed prevalent and incident infections was mainly studied. In an earlier article by our group that assessed a population with no evidence of prior exposure to HPV infection, the overall detection of incident HPV infections was 20.61% (95% confidence interval [CI], 18.47–22.99) per 1000 person-months during a 6-year follow-up among the placebo group of young women participating in a vaccine clinical trial.[7] In the current report, we present results on determinants of acquisition and clearance of HPV infections in the same cohort of women. Because participants had no serological or virological evidence of an HPV infection at study entry, we also stratified the analyses by length of follow-up to identify key determinants according to follow-up time and evaluate the robustness of findings.
MATERIALS AND METHODS
Study Design and Population
A detailed description of the clinical trial methodology including eligibility, sample collection methods, study sites, and recruitment and follow-up dates, has been published previously.[8,9] Briefly, a cohort of 1,113 women (age, 15–25 years) from 32 study centers in Brazil, the United States, and Canada were enrolled in a phase II double-blind, randomized controlled trial assessing the HPV-16/18 AS04-adjuvanted vaccine. Recruitment took place between July 2000 and September 2001. Figure 1 provides a schematic of study visits and follow-up phases. Participants followed up until 18 months were asked to participate in an extended follow-up period until 27 months after receiving the first vaccine dose [HPV-001 (NCT00689741, Clinicaltrials.gov)]. Upon completion of short-term follow-up (0–27 months), vaccine and placebo recipients were invited to participate in a long-term follow-up study of ∼3-year duration (up to 76 months) [HPV-007 (NCT00120848, Clinicaltrials.gov)]. Anonymized individual participant data and study documents can be requested for further research from www.clinicalstudydatarequest.com.
Figure 1
Schematic timeline and specimen collection. Cervicovaginal testing was done at every visit between months 0 and 27 (short-term follow-up), but was optional at month 21 and month 27 (Visits 10 and 12, respectively). Cervical testing was done at every visit between months 40 and 76 (long-term follow-up), but was also done at the screening visit and at month 6, month 12, and month 18 (visits 1, 4, 7, and 9, respectively). Cervicovaginal testing was optional at month 46, month 58, and month 70 (visits 14, 16, and 18, respectively).
Schematic timeline and specimen collection. Cervicovaginal testing was done at every visit between months 0 and 27 (short-term follow-up), but was optional at month 21 and month 27 (Visits 10 and 12, respectively). Cervical testing was done at every visit between months 40 and 76 (long-term follow-up), but was also done at the screening visit and at month 6, month 12, and month 18 (visits 1, 4, 7, and 9, respectively). Cervicovaginal testing was optional at month 46, month 58, and month 70 (visits 14, 16, and 18, respectively).The current ad-hoc analysis was restricted to 534 women randomized to the placebo arm (received only aluminum hydroxide adjuvant), who, at the screening visit, were seronegative for HPV16 and HPV18 antibodies and HPV DNA-negative, and had no more than 6 sexual partners and no prior abnormal Pap smear or treatment for cervical disease. Women who participated in the initial phase (0–18 months) and received 3 doses of placebo were eligible to participate in the extended and long-term follow-up phases. This study design provided a unique opportunity to specifically study determinants of incident infections as early events at the beginning of follow-up.The guidelines of the Declaration of Helsinki, the International Conference on Harmonization-Good Clinical Practices were followed. The study protocol was approved by independent ethics committees or institutional review boards of the respective study centers. Written informed consent was obtained at the screening visit from all participants or from a legally acceptable representative for those younger than the legal age of consent.
Specimen Collection and HPV DNA Testing
As illustrated in Figure 1, cervical samples were collected for cytology and HPV DNA testing at the screening visit; at months 6, 12, and 18 in the initial phase and at 6-month intervals in the long-term follow-up period. Self-collected cervicovaginal samples were obtained at months 0, 6, and then every 3 months until the end of the extended follow-up phase.HPV DNA testing for 25 genotypes was based on the broad-spectrum polymerase chain reaction (PCR) SPF10 LiPA25 system.[10] The genotypes tested included hrHPV (16, 18, 31, 33, 35, 39, 45, 51, 52, 56, 58, 59, 66, and 68) and low-risk HPV (lrHPV; 6, 11, 34, 40, 42, 43, 44, 53, 54, 70, and 74) types. Type-specific PCR was also done for HPVs 16 and 18.[11]
Risk Determinants
Participants completed a self-administered baseline questionnaire at the enrollment visit. Information on sociodemographic characteristics as well as behavioral, sexual, contraceptive, and other gynecological medical history was collected, covering known risk factors for HPV infection. The questionnaire was translated to the appropriate local language.
Statistical Analysis
Descriptive statistics were used to examine characteristics of the population by study site. Differences in distribution were assessed using the χ2 test. Univariate (adjusted for age and study site) and multivariable Cox proportional hazards models were used to estimate rate ratios (RRs) and 95% CI of HPV acquisition, separately for the short-term (0–27 months) and long-term (40–76 months) follow-up phases as well as for the entire follow-up (0–76 months). A parsimonious multivariable model was constructed using forward selection (stepwise, P < 0.15).The association with HPV clearance over the entire follow-up was modeled through Cox proportional hazards regression with adjustment for age and region. Clearance was defined in 2 ways: as a single HPV-negative visit after 1 or more HPV-positive visits (liberal clearance) or conditioning on at least 2 consecutively negative visits definition (conservative clearance).The following covariates were considered a priori as potential correlates of HPV acquisition and clearance: age (<20 years, ≥20 years), country/region (Brazil, North America), marital status (single/divorced, married/partner), ethnicity (White, black, others), education (sexually transmitted disease (no, yes), condom use (never, regularly, sometimes), and oral contraceptive use (never, regularly, sometimes). Analyses were performed for any HPV infection and grouped HPV infections (lrHPV, hrHPV) using combined clinician-collected cervical and self-collected cervicovaginal samples as we have previously shown comparable HPV detection between the 2 in this population.[7] Months since infection acquisition and clearance were used as the time metric. Data were analyzed with Stata (StataCorp. 2013. Stata Statistical Software: Release 13. College Station, TX: StataCorp LP.).
RESULTS
Table 1 presents baseline characteristics of the population stratified by study site. Mean age of participants was 20.79 years ±2.70 SD (standard deviation); women in Brazil were slightly younger (SD, 20.21 ± 2.93 years) compared with those in North America (21.26 years ±2.40 SD). Overall, the majority were single or divorced (77.3%), 69.1% were white, and 55.4% had beyond high school education. Regarding behavioral characteristics, 15.7% of women were current smokers, 10.7% had 2 or more lifetime sexual partners, 60.5% were 15 years and older at first sexual intercourse, 5.2% had a history of sexually transmitted diseases, 27.7% regularly used condoms, and 53.2% regularly used oral contraceptive. Comparison by study site indicated significant differences in certain characteristics (ie, age group, marital status, ethnicity, education, lifetime number of sex partners and age at first sexual intercourse), but not for smoking, history of any sexually transmitted disease, condom use, and oral contraceptive use.
TABLE 1
Baseline Characteristics of the Study Population (N = 534) by Study Region, Control Arm of a Phase II HPV-16/18 Vaccine Trial
Baseline Characteristics of the Study Population (N = 534) by Study Region, Control Arm of a Phase II HPV-16/18 Vaccine TrialTable 2 presents results from age- and region-adjusted Cox regression models of risk factors associated with new HPV infections within the initial, short-term follow-up. During this period, the cumulative incidence of any HPV infection was 48.1% (8535.80 person-months; incidence rate, 30.11; 95% CI, 26.64–34.02); 19.1% for lrHPV (10783.57 person-months; incidence rate, 9.46; 95% CI, 7.79–11.48) and 45.3% for hrHPV (8926.5 person-months; incidence rate, 27.11, 95% CI: 23.90–30.75) types. Married women had a significantly lower incidence of HPV (any, hrHPV, lrHPV) compared with single/divorced women. Women who reported 2 or more lifetime sex partners had more than twice the risk of incident infections with any type and with hrHPV types relative to those with 1 partner or none. We did not find any association between the number of lifetime sexual partners and lrHPV infection. Having a history of any sexually transmitted disease, regular condom use, and occasional oral contraceptive use were also significantly associated with any HPVinfection acquisition. We found no significant associations with age, study site, ethnicity, education, smoking, or age at first sexual intercourse.
TABLE 2
Age- and Region-Adjusted Associations Between Baseline Sociodemographic and Behavioral Factors and HPV Acquisition, Short-Term Follow-Up (0–27 Months)
Age- and Region-Adjusted Associations Between Baseline Sociodemographic and Behavioral Factors and HPV Acquisition, Short-Term Follow-Up (0–27 Months)Table 3 shows multivariable analysis results. Marital status, lifetime number of sex partners, history of any sexually transmitted disease, and occasional use of oral contraceptives retained their independent associations with acquisition of any HPV infection, in addition to ethnicity. Having 2 or more lifetime sex partners (adjusted RR, 2.03; 95% CI, 1.37–3.02) and a history of any sexually transmitted disease (adjusted RR, 1.98; 95% CI, 1.19–3.29) were also positively associated with risk of incident hrHPV infections. A protective association was still observed with marital status, whereas age, study site, education, and smoking were generally unrelated to HPV acquisition.
TABLE 3
Multivariable* Association Between Baseline Sociodemographic and Behavioral Factors and HPV Acquisition, Short-Term Follow-Up (0–27 Months)
Multivariable* Association Between Baseline Sociodemographic and Behavioral Factors and HPV Acquisition, Short-Term Follow-Up (0–27 Months)Supplemental Tables 1 http://links.lww.com/OLQ/A399 and 2 http://links.lww.com/OLQ/A400 present, respectively, age- and region- adjusted as well as multivariable analyses considering the entire follow-up (0–76 months), whereas corresponding findings restricted to long-term follow-up (40–76 months, ignoring outcomes prior to 40 months) are presented in Supplemental Tables 3 http://links.lww.com/OLQ/A401 and 4 http://links.lww.com/OLQ/A402. Overall, similar patterns of determinants were found irrespective of length of follow-up.Table 4 presents age- and region-adjusted associations of baseline sociodemographic and behavioral factors with clearance of incident HPV infections. A total of 226 (71.1%) of 318 incident infections during follow-up cleared, including 146 (88.5%) of 165 lrHPV and 265 (86.9%) of 305 hrHPV infections. No significant associations were observed between each of age, marital status, education, smoking, age at first sexual intercourse, history of any sexually transmitted disease, condom use, and oral contraceptive use with HPV clearance (any, lrHPV, and hrHPV types). There was a tendency toward increased hrHPV clearance among women in North America (adjusted RR, 1.38; 95% CI, 1.08–1.78) and among black women (adjusted RR, 1.64; 95% CI, 1.04–2.60). A greater number of lifetime sex partners was significantly associated with reduced clearance rates only for any HPV infection. Only the negative effect of multiple sexual partners remained when defining clearance conservatively as 2 or more consecutive negative visits with the index HPV outcome (Supplemental Table 5 http://links.lww.com/OLQ/A403).
TABLE 4
Age- and Region-Adjusted Associations Between Baseline Sociodemographic and Behavioral Factors and HPV Clearance (Liberal), Entire Follow-Up (0–76 Months)
Age- and Region-Adjusted Associations Between Baseline Sociodemographic and Behavioral Factors and HPV Clearance (Liberal), Entire Follow-Up (0–76 Months)
DISCUSSION
Understanding the role of risk factors for the acquisition and clearance of HPV infection will support the development of effective prevention and screening strategies. We undertook the present analysis to better understand factors associated with HPV acquisition and clearance in a prospectively followed cohort of young women with no evidence of previous exposure to HPV (based on well-characterized serological and cervical samples results) and the restriction to having 6 or fewer lifetime sexual partners. We also attempted to identify a time frame when these factors were most influential. We found lifetime number of sex partners, marital status, and history of sexually transmitted disease at baseline to be the most important determinants of HPV acquisition. There were no remarkable differences in patterns of risk determinants when stratifying the analyses by length of follow-up.Sexual behavior and age are well-recognized risk determinants for HPV infections. The number of lifetime sex partners emerged as a determinant of incident HPV infection, in line with previous studies. A high risk of acquiring a new HPV infection was reported in women aged 18–24 years soon after sexual debut with a first male sex partner, and the 1-year cumulative incidence was 28.5% which increased to 39.2% by 2 years and to 49.1% by 3 years.[12] In a longitudinal study on 97 virgin women and 105 monogamous women (ie, only had 1 sex partner) at enrollment, HPV acquisition was strongly determined by the accumulation of sexual partners; women with ≥3 partners had a 9-fold increased risk for acquiring HPV compared with women with 1 partner during follow-up.[13] Unfortunately, we did not collect information on new sexual partners during follow-up.In a meta-analysis of 194 studies published between 1995 and 2009, comprising more than one million women with normal cytological findings, the age-specific HPV distribution had a first peak in prevalence at younger ages (<25 years) just after sexual debut, a lower prevalence plateau at middle ages, followed by a varying rebound at older ages (≥45 years).[14] Similarly, HPV acquisition predominated at younger ages in a population-based cohort of 7,237 women in Costa Rica.[15] Although age in the current analysis was not a strong predictor of HPV acquisition, it is possible that the restricted age range of participants in this cohort (15–25 years) may have been too narrow to allow age-specific comparisons in risk. Various other studies found younger age and higher number of lifetime sex partners to be significant risk factors for new HPV infections,[16-19] but the patterns are not always consistent across populations.[20]Two early prospective studies on the epidemiology of HPV acquisition reported differences in risk-factor profiles according to HPV oncogenicity.[21,22] In a subset of 1425 participants aged 18 to 60 years in the Ludwig-McGill cohort longitudinal investigation of the natural history of HPV infection in Brazil, an association was found between age, age at first intercourse and oral contraceptive use with risk of hrHPV, but not with lrHPV, whereas number of sexual partners was associated with both hrHPV and lrHPV types.[21] In a cohort of 1610 women aged 15 to 85 years who were HPV negative and with normal cytological results at enrollment in Colombia, age and having new sex partners during follow-up were consistent predictors of oncogenic, but not non-oncogenic, HPV in multivariable analyses.[22] We found a strong protective association between married women (ie, more likely to be in monogamous relationship than divorced/separated women) and incident hrHPV and lrHPV, even after adjustment whereas for having a history of any sexually transmitted disease the observed increased risk remained statistically significant in the multivariable model only for hrHPV types.We did not find an association between smoking and HPV acquisition, consistent with some[18,22,23] but not with all studies.[21,24] Interestingly, condom use was associated with an increased risk of HPV (any, lrHPV, hrHPV) in multivariable analyses; statistical significance was attained only for lrHPV among regular condom users. It is also likely that women with an increased number of sex partners were more likely to use condoms.As for determinants of clearance of grouped HPV infections, these have been less frequently studied. Apart from finding an inverse association between an increased number of lifetime sex partners and clearance of any HPV, and a marginal increased risk of clearance of hrHPV types among black women and those in North America, no associations were found between other sociodemographic and behavioral factors and HPV persistence. Our finding that black women cleared their hrHPV infection more rapidly than whites is in keeping with the results from another study.[25] A cohort of 416 Canadian Inuit women also reported no association between sociodemographic, lifestyle or sexual factors with clearance of any incident HPV infections.[26] Other studies found that condom and oral contraceptive use conferred varying levels of protection in terms of HPV clearance based on oncogenic potential.[25,27]A limitation of our study is the lack of information on other epidemiological risk factors such as parity, genital hygiene, nutritional status, behaviors of sexual partners, and male circumcision.[28] The main strength of our study is that our sample had a low risk profile; participants were young healthy women who had a maximum of 6 lifetime sex partners. The fact that participants were young and ascertained to be HPV negative upon entry supports the contention that what we detected as incident infections may be presumed as true new infections. However, given the uncertainty that exists concerning latency, instances of new detection may represent the unmasking of a previously undetectable infection because of better sampling via exfoliation or fluctuations in productivity of infection foci.[29]It is plausible to assume that universal vaccination against some of the most clinically relevant HPV types will transform the epidemiology of these infections, a change that is already noticeably strong even in settings with suboptimal vaccination coverage.[30] Establishing an epidemiologic baseline, such as our study does, will help in assisting the design of surveillance mechanisms to monitor infection outcomes and behavioral risk factors in vaccinated populations.
Authors: Diane M Harper; Eduardo L Franco; Cosette M Wheeler; Anna-Barbara Moscicki; Barbara Romanowski; Cecilia M Roteli-Martins; David Jenkins; Anne Schuind; Sue Ann Costa Clemens; Gary Dubin Journal: Lancet Date: 2006-04-15 Impact factor: 79.321
Authors: M F Baay; W G Quint; J Koudstaal; H Hollema; J M Duk; M P Burger; E Stolz; P Herbrink Journal: J Clin Microbiol Date: 1996-03 Impact factor: 5.948
Authors: Rachel Bennett; Helen Cerigo; François Coutlée; Michel Roger; Eduardo L Franco; Paul Brassard Journal: Sex Transm Dis Date: 2015-05 Impact factor: 2.830
Authors: S K Kjaer; B Chackerian; A J van den Brule; E I Svare; G Paull; J M Walbomers; J T Schiller; J E Bock; M E Sherman; D R Lowy; C L Meijer Journal: Cancer Epidemiol Biomarkers Prev Date: 2001-02 Impact factor: 4.254
Authors: Harriet Richardson; Michal Abrahamowicz; Pierre-Paul Tellier; Gail Kelsall; Roxane du Berger; Alex Ferenczy; François Coutlée; Eduardo L Franco Journal: Cancer Epidemiol Biomarkers Prev Date: 2005-05 Impact factor: 4.254
Authors: Mahboobeh Safaeian; Mohammad Kiddugavu; Patti E Gravitt; Stephen J Gange; Joseph Ssekasanvu; Dan Murokora; Marc Sklar; David Serwadda; Maria J Wawer; Keerti V Shah; Ron Gray Journal: Cancer Epidemiol Biomarkers Prev Date: 2008-06 Impact factor: 4.254
Authors: Jaimie Z Shing; Marie R Griffin; Rachel S Chang; Alicia Beeghly-Fadiel; Staci L Sudenga; James C Slaughter; Manideepthi Pemmaraju; Edward F Mitchel; Pamela C Hull Journal: Am J Prev Med Date: 2021-10-29 Impact factor: 5.043
Authors: Marisa C Eisenberg; Rafael Meza; Andrew F Brouwer; Lora P Campredon; Heather M Walline; Brittany M Marinelli; Christine M Goudsmit; Trey B Thomas; Rachel L Delinger; Yan Kwan Lau; Emily C Andrus; Thankam Nair; Thomas E Carey Journal: BMJ Open Date: 2022-01-03 Impact factor: 2.692