A case-control study based on the National Health and Nutrition Examination Survey to evaluate the effects of human papilloma virus on bone health in women
===========================================================================================================================================================

* Xiang Li
* Guangjun Jiao
* Yunzhen Chen

## Abstract

**Background** Globally, both human papillomavirus (HPV) infection and osteoporosis (OP) are more prevalent in women than in men. It remains unclear whether HPV has an impact on bone health.

**Methods** This case-control study was based on data from the National Health and Nutrition Examination Survey (NHANES). Comparable datasets were created via the nearest neighbor propensity score matching (PSM) method (1:2). The Welch two-sample t test was used to analyze the association between HPV infection and bone mineral density (BMD). Restricted cubic spline (RCS) and Kendall’s tau-b tests were used to explore the effect of HPV infection type on BMD.

**Results** BMDs in the legs and lumbar spine were lower in subjects infected with HPV than in uninfected subjects. RCS analysis showed that the larger the number of cooccurring HPV types in a woman, the lower the BMD was. In addition, four HPV types were negatively associated with leg BMD, and 14 HPV types were negatively associated with lumbar spine BMD. HPV types 53, 59, and 89 had effects on both leg and lumbar spine BMDs.

**Conclusions** HPV infection was associated with a decrease in BMD. Appropriately designed trials can help determine whether interventions to prevent HPV infection will have a protective effect on BMD.

**Funding** This research was not supported by any specific grant from any funding agency in the public, commercial or not-for-profit sector.

Keywords
*   Women
*   Human papillomavirus
*   Bone mineral density
*   Correlation

## 1. Introduction

Osteoporosis (OP) is one of the most common chronic diseases of the skeleton. It is characterized by a decrease in bone mineral density (BMD) and the destruction of bone structure[1, 2]. OP is of great clinical importance because OP is a major risk factor for fracture and greatly affects human health[3, 4]. It has been reported that there are more than 200 million OP patients worldwide[5, 6]. In some regions, the prevalence of OP is approximately 40% in the elderly population[7-9]. OP and associated fractures impose substantial economic and social burdens. OP, a complex disease, is associated with a variety of factors. A better understanding of emerging risk factors contributing to the disease may enable the development of prevention and treatment strategies targeting high-risk individuals.

Human papillomavirus (HPV) is one of the most recognized viruses in humans, with a fairly high infection rate (2-44%)[10-12]. Most individuals are exposed to HPV infection at some point in their lifetime. HPV resulting in lesions on the genitalia has become the most common sexually transmitted infection in the world. Moreover, HPV is also recognized as a causative factor for cervical cancer[13]. It is also a risk factor for several tumors[14-19]. The rate of HPV infection is significantly higher among women than among men[20, 21]. This is similar to OP, for which female sex is a risk factor[22-24]. However, the role of HPV in the bone health of women is unclear, and no studies have confirmed the association between HPV and BMD.

Here, we conducted a large case-control study based on National Health and Nutrition Examination Survey (NHANES) to explore the association between HPV and BMD. Elucidating the strength and pattern of this association will provide evidence to guide appropriate targeted public health interventions.

## 2. Subjects and methods

### 2.1 General information about the NHANES

The National Center for Health Statistics of the Centers for Disease Control and Prevention in the United States conducted a multistage, stratified, large-scale, nationally representative study named the NHANES[25]. This cross-sectional survey collected demographic, dietary, examination, and questionnaire data. The NHANES protocol is described in detail elsewhere[26]. In this study, the data were obtained from the 2009–2014 NHANES database.

### 2.2 Definition of variables

#### 2.2.1 HPV infection

All eligible female participants underwent vaginal swab sample collected. The obtained samples were tested with the Roche linear array HPV genotyping test to determine the presence of HPV infection and the type of HPV in cases of infection. This method can detect 37 HPV types: 6, 11, 16, 18, 26, 31, 33, 35, 39, 40, 42, 45, 51, 52, 53, 54, 55, 56, 58, 59, 61, 62, 64, 66, 67, 68, 69, 70, 71, 72, 73, 81, 82, 83, 84, 89, and IS39.

We included subjects who were not infected with HPV in group A and those who were infected (all types) into group B. Subjects with missing data were excluded.

#### 2.2.2 Measurement of BMD

Dual energy X-ray absorptiometry was performed to detect the BMD of eligible female subjects. The scanning area included the whole body. However, in this study, we focused on only the BMD of the subjects’ legs and lumbar spines.

For the data from 2009 to 2010, the BMD of the leg was the average of bilateral total femur BMD, and the BMD of the lumbar spine was the average of L1–4 BMD. For the data from 2011 to 2014, the BMD value of the leg was the average BMD in both legs. In addition, if the BMD value of one leg was missing, the BMD in that leg was considered the value of the other leg. Subjects for whom the BMD in the legs or lumbar spine could not be calculated were excluded.

#### 2.2.3 Covariates

We included the following covariates of interest: age, race, education level, marital status, poverty-to-income ratio, smoking status, alcohol consumption status, diabetes, depression score and sleeping difficulty.

Smoking was defined as smoking ≥ 100 cigarettes in a lifetime[27]. Alcohol consumption was defined as having 4/5 or more drinks every day[28, 29].

Subjects were considered to have diabetes if they met one of the following criteria: told by doctors or other health professionals they had diabetes, insulin use, or diabetes medication use to lower their blood sugar[30].

Subjects with a Patient Health Questionnaire (PHQ)-9 total score of ≥ 10 were considered to have clinically relevant depression[8, 31]. Subjects who answered “yes” to the question “Have you ever been told by doctor that you have a sleep disorder?” were considered to have sleeping difficulty.

### 2.3 Ethical statement

The National Center for Health Statistics Research Ethics Review Board had approved the NHAHES. Therefore, NHANES data could be accessed without ethical or administrative permission.

### 2.4 Statistical analyses

We used normal Q-Q plots to test the distribution of the data. Welch two-sample t tests or nonparametric Mann–Whitney U tests were used to analyze continuous variables. Chi-square tests or Fisher’s exact tests were used to compare categorical variables.

Propensity score matching (PSM) has been shown to be effective in reducing selection bias in retrospective studies[7, 32, 33]. We first transformed BMD into a categorical variable based on whether BMD was ≥1. Subsequently, PSM, using the 1:2 nearest neighbor matching algorithm, was used to match subjects.

In the case of nonlinear a priori associations, restricted cubic spline (RCS) models were used to describe dose-response relationships between the continuous exposure and outcomes[34, 35]. In this study, RCS with five knots was used to further explore the association between HPV infection and BMD.

Finally, Kendall’s tau-b test was used to analyze the association between HPV type and BMD. In this analysis, we calculated only the individual effects of each HPV type and did not consider interactions between different types. Since there were no subjects infected with only HPV type 26, 64 or 82, we could not analyze their effects. P < 0.05 was considered statistically significant, and all tests were one-sided. Data extraction and merging, statistical analyses and figures were completed with R (version 4.1.1).

## 3. Results

### 3.1 Participant characteristics before matching

After screening, a total of 4673 female subjects were included in the study. The mean age of group A (n = 2614) was 39 ± 11.97 years, and the mean age of group B (n = 2059) was 36.67 ± 12.29 years. Supplemental table 1 shows the comparisons of the baseline characteristics between the two groups. The results indicated statistically significant differences for the vast majority of characteristics (P < 0.001). This suggested significant selection bias between the two groups of subjects. Therefore, performed PSM to eliminate this phenomenon.

View this table:
[Table 1.](http://medrxiv.org/content/early/2022/09/14/2022.09.13.22279899/T1)

Table 1. 
Comparison of BMD between the two groups after matching. BMD, bone mineral density.

### 3.2 Construction of matched datasets

After 1:2 matching, we obtained datasets for BMDs in the legs and lumbar spine. In total, 1066 subjects in group A and 533 subjects in group B were included. Analysis of baseline characteristics showed that the selection bias had been eliminated (P > 0.05). More detailed data are presented in Supplemental table 2 and Supplemental table 3. The BMD values of the legs and lumbar spine in the different age groups are shown in Fig. 1A. Normal Q-Q plots showed a significant normal distribution of BMD (Fig. 1B - D).

View this table:
[Table 2.](http://medrxiv.org/content/early/2022/09/14/2022.09.13.22279899/T2)

Table 2. 
The correlation between HPV types and BMD of leg.

View this table:
[Table 3.](http://medrxiv.org/content/early/2022/09/14/2022.09.13.22279899/T3)

Table 3. 
The correlation between HPV types and BMD of lumbar spine.

![Figure 1.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2022/09/14/2022.09.13.22279899/F1.medium.gif)

[Figure 1.](http://medrxiv.org/content/early/2022/09/14/2022.09.13.22279899/F1)

Figure 1. 
After matching, the bone mineral density (BMD) of leg and lumbar spine in different age groups (A). The normal Q-Q plots of BMD in group A (B, C) and group B (D, E).

### 3.3 Lower BMD in subjects with HPV infection

Welch two-sample t tests showed a lower BMD in subjects infected with HPV in both the legs and lumbar spine (P < 0.001). The decrease in BMD was greater in the lumbar spine than in the legs (Table 1). We performed RCS analysis to more clearly elucidate the association between HPV infection and BMD. In Fig. 2, we found a nonlinear association between the number of cooccurring HPV types and BMD (P for nonlinear association < 0.001). The larger the number of cooccurring HPV types, the lower the BMD was.

![Figure 2.](http://medrxiv.org/http://medrxiv.stage.highwire.org/content/medrxiv/early/2022/09/14/2022.09.13.22279899/F2.medium.gif)

[Figure 2.](http://medrxiv.org/content/early/2022/09/14/2022.09.13.22279899/F2)

Figure 2. 
Restricted cubic spline models for the relationship between the number of infected HPV types and bone mineral density (BMD) after matching.

Finally, we analyzed each HPV type in an attempt to identify those with the greatest impact on BMD. The data showed that four HPV types (06, 53, 59 and 89) were negatively correlated with leg BMD (P < 0.05). The strongest correlation was found for type 89 (r = -0.060). Fourteen HPV types were negatively correlated with lumbar spine BMD. Among them, type 53 had the strongest correlation (r = -0.085). There were three HPV types that had an effect on both leg and lumbar spine BMDs, namely, 53, 59 and\ 89. More detailed data are shown in Tables 2 and 3.

## 4. Discussion

To date, the literature on the link between viral infection and OP is scarce. Viruses that have been shown to be associated with OP include human immunodeficiency virus (HIV), hepatitis B virus, hepatitis C virus, and severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)[36-40]. In the present study, we demonstrated that HPV infection led to a decrease in BMD. Combining the results of the PSM and RCS analyses, we concluded that the association between the two showed a nonlinear negative correlation. The larger the number of cooccurring HPV types, the lower the BMD was.

As mentioned earlier, we confirmed the association between HPV and BMD. However, it is important to determine the causal relationship between them. It often takes a long time for BMD to decline and eventually progress to OP. Compared to this process, HPV infection can be considered an immediate event. Therefore, we believe it is more likely that HPV infection caused the decline in BMD rather than the reverse. In addition, the results after statistical adjustment support our view. PSM eliminated selection bias, prevented our findings from being confounded by other factors, and improved the level of evidence.

However, the biological mechanisms underlying the decline in BMD due to HPV remain unclear. It is known that persistent HPV infection triggers a range of responses, including immunosuppression and abnormal cell destruction[41]. Scott et al. conducted a prospective study and found increased IL-12 concentrations in cervicovaginal lavage fluid specimens in women with persistent HPV infection[42]. This may reflect IL-12-induced macrophage activation and associated subsequent chronic inflammation[43]. The tissue damage caused by inflammation is self-explanatory. In addition, cervical HPV infection also causes infiltration of macrophages and dendritic cells. The activation of these cells also leads to the release of TNF, IL-8 and IL-12[44, 45]. In addition to effects on the immune system, effects of HPV on lipid metabolism have also generated concerns. Lipid metabolism (digestion, uptake, synthesis and peroxidation) is crucial in high-risk HPV infection. The lipid component of the cell membrane makes HPV infection more efficient[46]. Oxidative stress and lipid peroxidation induced by HPV infection also promote cellular damage that may lead to cell death[47]. Although these effects are more pronounced in HPV-infected local tissues, prolonged effects likely lead to systemic damage, including bone damage.

We also noted a larger decline in BMD in the lumbar spine than in the legs in HPV-infected individuals. Fourteen HPV types were negatively associated with lumbar spine BMD, and this number was much higher than that associated with the legs (four types). This finding suggests that lumbar spine BMD is more sensitive than leg BMD. This finding is consistent with the findings of some previous studies to some extent. Snyder et al. found that testosterone increased BMD in older men and that the increase in spinal BMD was larger than that in the hip[48].

Finally, we would like to discuss HPV types 53, 59 and 89. All three types were negatively associated with BMDs in the lumbar spine and legs. Public HPV preventive measures and the detection of HPV infection should additionally focus on these three HPV types. These three HPV types could also help reveal the biological mechanisms underlying the association of HPV with BMD.

Although our statistical analysis attenuated selection bias, some limitations remain in the present study. First, prospective evidence is still needed to truly determine the causal association between HPV infection and BMD decline. Second, elucidating the biological mechanisms of the association between HPV and BMD may require additional molecular experiments. Third, only 37 HPV types were included in this study. However, there are many more types of HPV in existence. How HPV types not detected in the NHANES affect bones health in women is still unknown.

In conclusion, we used data from a nationally representative database of US adults and found a negative association between HPV infection and BMD in female subjects. The larger the number of cooccurring HPV types, the lower the BMD was.

## Contributors

Yunzhen Chen and Guangjun Jiao planned the study and interpreted the analyses. Xiang Li collected data, planned the study, performed the analyses, produced the figures and tables, and drafted the manuscript. All authors contributed to the interpretation of the results and revision of the paper and approved the final manuscript.

## Data Availability

The data generated during and/or analyzed during the current study are available from the corresponding author upon reasonable request.

## Declaration of interests

The authors state that they have no conflicts of interest.

## Data sharing statement

The data generated during and/or analyzed during the current study are available from the corresponding author upon reasonable request.

View this table:
[Supplemental table 1.](http://medrxiv.org/content/early/2022/09/14/2022.09.13.22279899/T4)

Supplemental table 1. Baseline characteristics for the total subjects before matching.

View this table:
[Supplemental table 2.](http://medrxiv.org/content/early/2022/09/14/2022.09.13.22279899/T5)

Supplemental table 2. Baseline characteristics for subjects of leg BMD after matching.

View this table:
[Supplemental table 3.](http://medrxiv.org/content/early/2022/09/14/2022.09.13.22279899/T6)

Supplemental table 3. Baseline characteristics for subjects of lumbar

## Acknowledgments

Not applicable.

*   Received September 13, 2022.
*   Revision received September 13, 2022.
*   Accepted September 14, 2022.


*   © 2022, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution 4.0 International), CC BY 4.0, as described at [http://creativecommons.org/licenses/by/4.0/](http://creativecommons.org/licenses/by/4.0/)

## References

1.  1.Zhu, X.,  W. Bai, and  H. Zheng, Twelve years of GWAS discoveries for osteoporosis and related traits: advances, challenges and applications. Bone Res, 2021. 9(1): p. 23.
    
    

2.  2.Bellavia, D., et al., Non-flavonoid polyphenols in osteoporosis: preclinical evidence. Trends Endocrinol Metab, 2021. 32(7): p. 515–529.
    
    

3.  3.Yong, E.L. and  S. Logan, Menopausal osteoporosis: screening, prevention and treatment. Singapore Med J, 2021. 62(4): p. 159–166.
    
    

4.  4.Yu, X.H., et al., Birth weight is positively associated with adult osteoporosis risk: observational and Mendelian randomization studies. J Bone Miner Res, 2021. 36(8): p. 1469–1480.
    
    

5.  5.Khatib, J.,  K. Stote, and  A.R. Gosmanov, Utility of DXA screening for diagnosis of osteoporosis in US veterans aged 70 years and older. J Investig Med, 2018. 66(2): p. 298–303.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamltIjtzOjU6InJlc2lkIjtzOjg6IjY2LzIvMjk4IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDkvMTQvMjAyMi4wOS4xMy4yMjI3OTg5OS5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

6.  6.Yang, T.S., et al., Osteoporosis: prevalence in Taiwanese women. Osteoporos Int, 2004. 15(4): p. 345–7.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=14872301&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F09%2F14%2F2022.09.13.22279899.atom) 

7.  7.Rosenbaum, P.R. and  D.B. Rubin, The central role of the propensity score in observational studies for causal effects. Biometrika, 1983. 70(1): p. 41–55.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/biomet/70.1.41&link_type=DOI) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1983QH66900005&link_type=ISI) 

8.  8.Yu, B., et al., Trends in depression among Adults in the United States, NHANES 2005-2016. J Affect Disord, 2020. 263: p. 609–620.
    
    

9.  9.Zhang, Q., et al., Prevalence and contributing factors of osteoporosis in the elderly over 70 years old: an epidemiological study of several community health centers in Shanghai. Ann Palliat Med, 2020. 9(2): p. 231–238.
    
    

10. 10.Baseman, J.G. and  L.A. Koutsky, The epidemiology of human papillomavirus infections. J Clin Virol, 2005. 32 Suppl 1: p. S16–24.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jcv.2004.12.008&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15753008&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F09%2F14%2F2022.09.13.22279899.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000228272000003&link_type=ISI) 

11. 11.Gjøoen, K., et al., Prevalence of human papillomavirus in cervical scrapes, as analyzed by PCR, in a population-based sample of women with and without cervical dysplasia. Apmis, 1996. 104(1): p. 68–74.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=8645461&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F09%2F14%2F2022.09.13.22279899.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1996UA58700011&link_type=ISI) 

12. 12.Bosch, F.X. and  S. de Sanjosé, Chapter 1: Human papillomavirus and cervical cancer--burden and assessment of causality. J Natl Cancer Inst Monogr, 2003(31): p. 3–13.
    
    

13. 13.Peltecu, G., et al., Human papilloma virus and cervical preinvasive disease. J Med Life, 2009. 2(4): p. 373–7.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20108750&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F09%2F14%2F2022.09.13.22279899.atom) 

14. 14.Rezazadeh, A., et al., The role of human papilloma virus in lung cancer: a review of the evidence. Am J Med Sci, 2009. 338(1): p. 64–7.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19593105&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F09%2F14%2F2022.09.13.22279899.atom) 

15. 15.Torrente, M.C. and  J.M. Ojeda, Exploring the relation between human papilloma virus and larynx cancer. Acta Otolaryngol, 2007. 127(9): p. 900–6.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1080/00016480601110238&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17712666&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F09%2F14%2F2022.09.13.22279899.atom) 

16. 16.Jarzyński, A., et al., Occurrence of BK Virus and Human Papilloma Virus in colorectal cancer. Ann Agric Environ Med, 2017. 24(3): p. 440–445.
    
    

17. 17.Lightner, A.L., Human Papilloma Virus and Anal Squamous Cell Cancer in IBD: Is It Time to Update Our Practice Parameters? Dis Colon Rectum, 2017. 60(12): p. 1231–1232.
    
    

18. 18.Stratton, K.L. and  D.J. Culkin, A Contemporary Review of HPV and Penile Cancer. Oncology (Williston Park), 2016. 30(3): p. 245–9.
    
    

19. 19.de Martel, C., et al., Worldwide burden of cancer attributable to HPV by site, country and HPV type. Int J Cancer, 2017. 141(4): p. 664–670.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/ijc.30716&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28369882&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F09%2F14%2F2022.09.13.22279899.atom) 

20. 20.Schiffman, M.H., Recent progress in defining the epidemiology of human papillomavirus infection and cervical neoplasia. J Natl Cancer Inst, 1992. 84(6): p. 394–8.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/jnci/84.6.394&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=1311392&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F09%2F14%2F2022.09.13.22279899.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1992HU83300011&link_type=ISI) 

21. 21.Kjaer, S.K., et al., High-risk human papillomavirus is sexually transmitted: evidence from a follow-up study of virgins starting sexual activity (intercourse). Cancer Epidemiol Biomarkers Prev, 2001. 10(2): p. 101–6.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiY2VicCI7czo1OiJyZXNpZCI7czo4OiIxMC8yLzEwMSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIyLzA5LzE0LzIwMjIuMDkuMTMuMjIyNzk4OTkuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

22. 22.Wang, Y., et al., Associated Factors for Osteoporosis and Fracture in Chinese Elderly. Med Sci Monit, 2019. 25: p. 5580–5588.
    
    

23. 23.Gregson, C.L., et al., UK clinical guideline for the prevention and treatment of osteoporosis. Arch Osteoporos, 2022. 17(1): p. 58.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F09%2F14%2F2022.09.13.22279899.atom) 

24. 24.Porter, J.L. and  M. Varacallo, Osteoporosis, in StatPearls. 2022, StatPearls Publishing Copyright © 2022, StatPearls Publishing LLC.: Treasure Island (FL).
    
    

25. 25.Juul, F., et al., Ultra-processed food consumption among US adults from 2001 to 2018. Am J Clin Nutr, 2022. 115(1): p. 211–221.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ajcn/nqab305&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=34647997&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F09%2F14%2F2022.09.13.22279899.atom) 

26. 26.Mei, Z., et al., Physiologically based serum ferritin thresholds for iron deficiency in children and non-pregnant women: a US National Health and Nutrition Examination Surveys (NHANES) serial cross-sectional study. Lancet Haematol, 2021. 8(8): p. e572–e582.
    
    

27. 27.Semenov, Y.R., et al., Psoriasis and mortality in the United States: Data from the National Health and Nutrition Examination Survey. J Am Acad Dermatol, 2021. 85(2): p. 396–403.
    
    

28. 28.Gay, I.C.,  D.T. Tran, and  D.W. Paquette, Alcohol intake and periodontitis in adults aged ≥30 years: NHANES 2009-2012. J Periodontol, 2018. 89(6): p. 625–634.
    
    

29. 29.Courtney, J.B.,  M.A. Russell, and  D.E. Conroy, Tobacco and cannabis use as moderators of the association between physical activity and alcohol use across the adult lifespan in the United States: NHANES, 2005–2016. Preventive Medicine, 2022. 155: p. 106931.
    
    

30. 30.Al-Ibrahim, A.A. and  R.T. Jackson, Healthy eating index versus alternate healthy index in relation to diabetes status and health markers in U.S. adults: NHANES 2007-2010. Nutr J, 2019. 18(1): p. 26.
    
    

31. 31.Kroenke, K.,  R.L. Spitzer, and  J.B. Williams, The PHQ-9: validity of a brief depression severity measure. J Gen Intern Med, 2001. 16(9): p. 606–13.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1046/j.1525-1497.2001.016009606.x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=11556941&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F09%2F14%2F2022.09.13.22279899.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000171184700005&link_type=ISI) 

32. 32.Cho, S.M., et al., Cerebrovascular Events in Patients With Centrifugal-Flow Left Ventricular Assist Devices: Propensity Score-Matched Analysis From the Intermacs Registry. Circulation, 2021. 144(10): p. 763–772.
    
    

33. 33.Kaster, T.S., et al., Risk of serious medical events in patients with depression treated with electroconvulsive therapy: a propensity score-matched, retrospective cohort study. Lancet Psychiatry, 2021. 8(8): p. 686–695.
    
    

34. 34.Desquilbet, L. and  F. Mariotti, Dose-response analyses using restricted cubic spline functions in public health research. Stat Med, 2010. 29(9): p. 1037–57.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/sim.3841&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20087875&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F09%2F14%2F2022.09.13.22279899.atom) 

35. 35.Li, S.,  W. Sun, and  D. Zhang, Association of Zinc, Iron, Copper, and Selenium Intakes with Low Cognitive Performance in Older Adults: A Cross-Sectional Study from National Health and Nutrition Examination Survey (NHANES). J Alzheimers Dis, 2019. 72(4): p. 1145–1157.
    
    

36. 36.Wijarnpreecha, K., et al., Hepatitis C virus infection and risk of osteoporosis: A meta-analysis. Saudi J Gastroenterol, 2017. 23(4): p. 216–221.
    
    

37. 37.Disser, N.P., et al., Musculoskeletal Consequences of COVID-19. J Bone Joint Surg Am, 2020. 102(14): p. 1197–1204.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.2106/jbjs.20.00847&link_type=DOI) 

38. 38.Hampson, G., et al., Diagnosis and Management of Osteoporosis During COVID-19: Systematic Review and Practical Guidance. Calcif Tissue Int, 2021. 109(4): p. 351–362.
    
    

39. 39.Dong, H.V., et al., Osteoporosis and fractures in HIV/hepatitis C virus coinfection: a systematic review and meta-analysis. Aids, 2014. 28(14): p. 2119–31.
    
    

40. 40.Xie, X., et al., Association between hepatitis B virus infection and risk of osteoporosis: a systematic review and meta-analysis: A protocol for systematic review. Medicine (Baltimore), 2020. 99(16): p. e19719.
    
    

41. 41.Moerman-Herzog, A. and  M. Nakagawa, Early Defensive Mechanisms against Human Papillomavirus Infection. Clin Vaccine Immunol, 2015. 22(8): p. 850–7.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiY2RsaSI7czo1OiJyZXNpZCI7czo4OiIyMi84Lzg1MCI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIyLzA5LzE0LzIwMjIuMDkuMTMuMjIyNzk4OTkuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

42. 42.Scott, M.E., et al., Cervical cytokines and clearance of incident human papillomavirus infection: Hawaii HPV cohort study. Int J Cancer, 2013. 133(5): p. 1187–96.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/ijc.28119&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23436563&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F09%2F14%2F2022.09.13.22279899.atom) 

43. 43.Gately, M.K., et al., The interleukin-12/interleukin-12-receptor system: role in normal and pathologic immune responses. Annu Rev Immunol, 1998. 16: p. 495–521.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1146/annurev.immunol.16.1.495&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9597139&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F09%2F14%2F2022.09.13.22279899.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000073129400018&link_type=ISI) 

44. 44.Brennan, K. and  A.G. Bowie, Activation of host pattern recognition receptors by viruses. Curr Opin Microbiol, 2010. 13(4): p. 503–7.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.mib.2010.05.007&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20538506&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F09%2F14%2F2022.09.13.22279899.atom) 

45. 45.Coleman, N., et al., Immunological events in regressing genital warts. Am J Clin Pathol, 1994. 102(6): p. 768–74.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ajcp/102.6.768&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=7801889&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F09%2F14%2F2022.09.13.22279899.atom) 

46. 46.Siddiqa, A.,  J. Broniarczyk, and  L. Banks, Papillomaviruses and Endocytic Trafficking. Int J Mol Sci, 2018. 19(9).
    
    

47. 47.Cruz-Gregorio, A., et al., Lipid metabolism and oxidative stress in HPV-related cancers. Free Radical Biology and Medicine, 2021. 172: p. 226–236.
    
    

48. 48.Snyder, P.J., et al., Effect of Testosterone Treatment on Volumetric Bone Density and Strength in Older Men With Low Testosterone: A Controlled Clinical Trial. JAMA Intern Med, 2017. 177(4): p. 471–479.