Sex-Specific Associations of Testosterone and Genetic Factors With Health Span

Background: Previous studies have suggested associations between testosterone, genetic factors, and a series of complex diseases, but the associations with the lifespan phenotype, such as health span, remain unclear.

Methods: In this prospective cohort study, we analyzed 145,481 men and 147,733 women aged 38–73 years old from UK Biobank (UKB) to investigate the sex-specific associations of total testosterone (TT), free testosterone (FT), or polygenic risk score (PRS) with health span termination (HST) risk. At baseline, serum testosterone levels were measured. HST was defined by eight events strongly associated with longevity. PRS, an efficient tool combining the effect of common genetic variants to discriminate genetic risk of complex phenotypes, was constructed by 12 single-nucleotide polymorphisms related to health span from UKB (P ≤ 5.0 × 10⁻⁸). We used multivariable Cox regression models to estimate hazard ratios (HRs) and 95% confidence intervals (CIs).

Results: With a median follow-up time of 7.70 years, 26,748 (18.39%) men and 18,963 (12.84%) women had HST. TT was negatively associated with HST in men [HR per standard deviation (SD) increment of log-TT: 0.92, 95% CI: 0.88–0.97]. Inversely, both TT (HR per SD increment of log-TT: 1.05, 95% CI: 1.02–1.08) and FT (HR per SD increment of log-FT: 1.08, 95% CI: 1.05–1.11) presented an increased risk of HST in women. PRS was positively associated with HST risk (quintile 5 versus quintile 1, men, HR: 1.19, 95% CI: 1.15–1.24; women, HR: 1.21, 95% CI: 1.16–1.27). Moreover, men with high TT and low genetic risk showed the lowest HST risk (HR: 0.80, 95% CI: 0.73–0.88), whereas HST risk for women with both high TT and genetic risk increased obviously (HR: 1.32, 95% CI: 1.19–1.46). Similar joint effects were observed for FT in both genders.

Conclusions: We observed sex-specific associations that testosterone was negatively associated with HST risk in men and positively associated with HST risk in women. Genetic factors increased the HST risk, suggesting that participants with both high genetic risk and abnormal testosterone levels (high level in women or low level in men) should be the target for early intervention. Although our findings highlight the associations between testosterone and health span, further mechanistic studies and prospective trials are warranted to explore the causation behind.

Introduction

Aging is becoming one of the most substantial social transformations of developed countries in the twenty-first century. Although the average life expectancy in most developed countries has nearly doubled over the past two centuries (1), the increase of health span is more significant and meaningful than simply living longer. Briefly, “health span” is a period of good health during the lifetime and is also termed as healthy and morbidity-free lifespan (2). In many supercentenarians, health span approximates lifespan due to a consistent delay in the occurrence of age-related diseases (3), which has stimulated enormous interest in mysteries behind. Hence, health span becomes a promising phenotype for longevity research.

Interestingly, men usually have a shorter life expectancy than women, suggesting men seem to age slightly faster than women (4). Considering the important effect of sex hormones on maintaining individual muscle mass, bone mass, and physical function, it is rational to speculate that sex hormones have a potential relationship with individual aging. Testosterone is a sex hormone mainly produced by testes in men or synthesized by ovaries and adrenal glands of women in small amounts. Majority of testosterone circulating in the blood is tightly bound to sex hormone-binding globulin (SHBG) or weakly bound to albumin. Approximately only 1%–3% of them is unbound and readily available to the tissues, referred to as free testosterone (FT) (5). It is worth noting that the concentrations of serum testosterone and free testosterone decrease with age in both genders (6), and deficiencies in multiple anabolic hormones have been shown to associate with individual longevity (7). Previous studies have revealed the associations between testosterone and several complex diseases, such as cardiovascular diseases (8), metabolic disease (9), and cancer (10, 11), as well as mortality (12, 13). However, direct evaluation on the effect of testosterone for health span is limited. Additionally, given sex differences in serum testosterone concentrations and biological mechanisms, sex-specific analyses for health span are also necessary.

Recently, a genome-wide association study (GWAS) systematically selected eight health events closely connected with aging to construct a reasonable definition for health span, and further identified 12 single-nucleotide polymorphisms (SNPs) related to health span (14). This provides a great chance to establish polygenic risk score (PRS) for health span. To our knowledge, PRS is an efficient tool estimating the cumulative effect of multiple risk-associated variants to discriminate populations at high risk of complex diseases (15), such as lung cancer (16) and breast cancer (17). Considering that an interplay between serum testosterone and genetic factors may exist, investigation for their joint association with health span is essential as well.

Therefore, based on complete clinical and genetic data from UK Biobank (UKB), we explored the sex-specific associations between endogenous testosterone, genetic factors, and health span by a large-scale sample including 293,214 participants.

Materials and Methods

Study Participants

The UK Biobank is an ongoing national prospective cohort study with a detailed protocol publicly available (18). In brief, nearly 500,000 volunteers aged 40 to 70 years were enrolled from 22 assessment centers across England, Wales, and Scotland between 2006 and 2010. At baseline, characteristics of participants are assessed by a self-administered questionnaire, physical assessments, brief interview, and blood collection. The UK Biobank study has approval by the North West Multi-Centre Research Ethics Committee (http://www.ukbiobank.ac.uk/ethics/).

In this study, we successively excluded 72,477 and 29,027 participants whose health span had terminated at baseline according to in-patient hospital admissions data (UKB category 2000) and self-reported diagnoses obtained via verbal interview (UKB category 100074). Additionally, 23,721 non-white participants and 84,068 participants with missing data on testosterone, SHBG, or albumin were both excluded. Finally, a total of 145,481 men and 147,733 women were included for the following analyses (Supplementary Figure 1).

Assessment of Biomarkers

Details about serum biomarker measurements have been described elsewhere (19, 20). In short, serum samples were prepared and stored at −80°C until assayed for testosterone, albumin, and SHBG in the UK Biobank central laboratory. Then, serum concentrations of TT and SHBG were measured on the platform of Beckman Coulter DXI 800 by chemiluminescent immunoassays, with reportable ranges of 0.35–55.52 and 0.33–242 nmol/L, respectively. Albumin was assayed on Beckman Coulter AU5800 by a colorimetric method, with a reportable range of 15–60 g/L. The coefficients of variation of internal quality control samples with known high, medium, and low concentrations for TT, SHBG, and albumin ranged from 3.66% to 8.34%, 5.22% to 5.67%, and 2.09% to 2.20%, respectively. Moreover, the assay of serum biomarkers was registered with an external quality assurance scheme to verify the accuracy of measurement, and the high proportions of participated distributions for TT (99%), SHBG (95%), and albumin (97%) showed good or acceptable assay performance. FT concentrations were calculated from assayed values of TT, SHBG, albumin, and the association constants for the binding of testosterone to SHBG and albumin by using the Vermeulen equation (21).

Ascertainment of Outcome

Following the law of Gompertz dynamics (22), Zenin et al. identified a cluster of the top 8 morbidities with an independent clinical manifestation to define health span. These events included congestive heart failure (CHF), myocardial infarction (MI), chronic obstructive pulmonary disease (COPD), stroke, dementia, diabetes, dancer, and death (14). In this study, participants that had the first experience of any of these events mentioned above during the follow-up were considered to have “health span termination (HST)”, which marked a shift into a biological or clinical phase with a hazard of frailty, multimorbidity, and premature death accumulated progressively.

We extracted data from the National Cancer Registry (UKB data category 100092) and National Death Registry (UKB data category 100093) to define cancer and death events and followed a list of hospital data codes and self-reported data codes to define the other six events (Supplementary Table 1). To ensure consistency of updated time between these three databases, we set “December 14, 2016” as the end of follow-up. We calculated follow-up time from the date of attending an assessment center until the date of HST or “December 14, 2016”, whichever occurred first.

Polygenic Risk Score Calculation

A detailed procedure for genetic data in UKB has been reported elsewhere (23). We included 12 SNPs that were genome-wide significantly (P ≤ 5.0 × 10⁻⁸) associated with health span in Caucasians (Supplementary Table 2) (14). PRS was created using a weighted strategy under an additive model as follows:

$$PRS={\beta }_1\times SN{P}_{i,1}+{\beta }_2\times SN{P}_{i,2}+\cdots +{\beta }_j\times SN{P}_{i,j}$$

Briefly, the dosage of risk allele for each individual (SNP_ij) was summed after weighting with its respective effect size of the site-specific phenotype (ß_j).

Statistical Analyses

Cox proportional hazards regression models were used to estimate hazard ratios (HRs) and 95% confidence intervals (CIs) of the associations between testosterone and HST risk, with adjustment for age, menopause status (women only), Townsend deprivation index, education, body mass index (BMI), smoking status, alcohol intake frequency, physical activity, healthy diet (24), family history of cancer and cardiac-cerebral vascular disease (CCVD), aspirin/ibuprofen use, and hormone replacement therapy use (HRT, women only). To clarify the independent effect of TT, sex hormone-binding globulin (SHBG) was additionally adjusted in the fully adjusted model of TT (see Supplementary Methods for details) (12, 13). We also explored the linear association of testosterone after natural log transformation. The proportional hazard assumption was assessed based on Schoenfeld residuals.

Then, we visualized the concrete associations between testosterone and risk of HST by multivariate cubic regression splines with 4 knots for further exploration. Likelihood ratio was adopted to test potential non-linearity by comparing the model with only a linear term against the model with both linear and cubic spline terms. We also separately assessed the association of testosterone with cause-specific incidence [death, cancer (colorectal, lung, prostate, breast, ovary and endometrial cancer), congestive heart failure, chronic obstructive pulmonary disease, myocardial infarction, dementia, diabetes, and stroke] for better understanding the internal association with health span. Moreover, joint effects of TT and SHBG were also evaluated because of potential biological interplay, and stratification analyses were conducted to explore the heterogeneity of associations across the subgroups.

Models with and without the interaction terms between health span PRS and testosterone were compared using likelihood ratio tests to evaluate the statistical significance of potential effect modification. Participants were categorized into low (the bottom quintile, Q1), moderate (the second to fourth quintiles, Q2–Q4), and high (the top quintile, Q5) genetic strata. A Cox regression model, additionally adjusted for the top 10 genetic principal components (PC1–10) and genotyping chip, was executed to examine the risk of HST among those with a high or intermediate genetic risk, compared with participants with low genetic risk. Furthermore, combined analyses of the risk of HST were performed across different genetic risk and TT/FT categories. Finally, sensitivity analyses were performed to test the robustness of the observed associations by excluding participants who had HST in the first 2 years or with poor self-reported health status at baseline, as well as those with outliers of testosterone (top 1% and bottom 1%). We also additionally adjusted the fasting time and menstrual cycle proxy variables (“Time since last menstrual period”, “Length of menstrual cycle”, and “Menstruating today”) to control the potential bias introduced by these variables.

All analyses were performed in R (version 3.6.3). Cox regression and multivariate cubic regression splines were performed by “survival” and “rms” R packages. P <0.05 was statistically significant.

Results

Totally, the median follow-up time was 7.70 years (interquartile range: 6.80–8.45 years). Nearly half of the HST events were caused by cancer (46.3%), followed by MI (17.8%) and death (11.4%) (Supplementary Table 3).

The baseline characteristics of the participants by sex and health span are described in Table 1. There were 26,748 (18.39%) men and 18,963 (12.84%) women with HST during the follow-up. Men with HST were older and had a lower TT and FT level. They were more likely to have higher BMI and Townsend deprivation index. Moreover, they also tended to be current smokers and have a family history of CCVD or cancer. Similar characteristics were observed in women with HST. The majority of these women had reached menopause, were hormone replacement therapy users, and had higher levels of FT than the healthy women. Additionally, characteristics of participants by quintiles of total testosterone were also performed (Supplementary Table 4). The non-normality distribution of TT and FT was exhibited clearly, and the levels of TT and FT were totally different in both genders (Supplementary Figure 2).

TABLE 1

Table 1 Baseline characteristics of participants with or without terminated health span.

The sex-specific associations between TT, FT, and risk of HST are exhibited in Table 2. For men, those in the second to fourth quintiles (Q2–Q4) of TT all had a lower HST risk in the fully adjusted models with the bottom quintile (Q1) of TT as a reference (Q2, HR: 0.94, 95% CI: 0.91–0.98; Q3, HR: 0.95, 95% CI: 0.91–0.99; Q4, HR: 0.94, 95% CI: 0.90–0.98), and the hazard decreased by 8% per SD (standard deviation) increment of log-TT (HR: 0.92, 95% CI: 0.88–0.97). However, the association of FT in men diminished in the fully adjusted model. Compared with women in the bottom quintile (Q1), only those in the top quintile (Q5) of TT (HR: 1.08, 95% CI: 1.04–1.13) and FT (HR: 1.12, 95% CI: 1.07–1.18) had a significantly increased HST risk.

TABLE 2

Table 2 Associations of total testosterone and free testosterone with risk of health span termination in both genders.

Non-linear relationships of testosterone and the risk of HST were observed for both genders. In men, an approximate L-shape curvilinear relationship was observed, both in TT and FT, using multivariate cubic models (Figures 1A, B). After natural log transformation for TT and FT concentrations, negative associations were observed below the corresponding medians (TT, HR: 0.88, 95% CI: 0.82–0.95; FT, HR: 0.87, 95% CI: 0.81–0.93). Interestingly, TT and FT had an inverse association with HST risk in women. The relationship was linear for TT, while it was more like a J-shape for FT (Figures 1C, D). Above the medians, HST risks were positively correlated with TT (HR: 1.13, 95% CI: 1.05–1.21) and FT (HR: 1.16, 95% CI: 1.09–1.22).

FIGURE 1

Figure 1 Non-linear associations of total testosterone/free testosterone (TT/FT) with risk of health span termination by multivariate cubic regression splines. (A) TT in men, (B) FT in men, (C) TT in women, and (D) FT in women. Solid lines were point estimation of hazard ratio, and shaded regions represented the 95% CIs. The dotted lines were the median values. Fully adjusted model: age, menopause status (for women), college or university degree, deprivation index, body mass index, smoking status, alcohol drinking, IPAQ group, healthy diet, family history of cardiac-cerebral vascular disease (CCVD) or cancer, use of aspirin/ibuprofen, and hormone replacement therapy (for women). SHBG was additionally adjusted in total testosterone. HR, hazard ratio; SD, standard deviation.

Due to the complexity of health span, we further explored the relationship between testosterone and cause-specific incidence. In men, testosterone was observed to have decreased the death, COPD, dementia, diabetes, and lung cancer risk, as well as the overall risk of non-cancer events, but increased the prostate cancer and myocardial infarction risk (Supplementary Table 5). In women, testosterone was positively associated with death, breast cancer, ovary cancer, endometrial cancer, and overall cancer, as well as dementia (total testosterone only) and diabetes (free testosterone only), but negatively related to COPD (Supplementary Table 6). In the stratified analyses, the associations of HST with testosterone were largely consistent across subgroups, except alcohol intake and Townsend deprivation index in men (Supplementary Table 7) and women with a family history of CCVD (Supplementary Table 8). In the joint analyses between TT and SHBG, the HST risk increased from 10% to 25% as the increase of TT level in women with the lowest strata of SHBG (Supplementary Figure 3).

To explore the associations between genetic factors and health span, a polygenetic risk score was constructed by 12 SNPs of health span to evaluate individual genetic risk. As shown by the distribution of polygenetic risk scores, both men and women with HST showed a higher genetic risk (Supplementary Figure 4). Furthermore, we found a positive relationship between PRS and HST risk. Participants from both genders with high genetic risk tended to terminate health span prematurely, when compared with those with low genetic risk (men: HR: 1.19, 95% CI: 1.15–1.24; women: HR: 1.21, 95% CI: 1.16–1.27, Figure 2). However, no interaction between PRS and testosterone was observed (P > 0.05), and the associations of testosterone with HST were unchanged after additionally adjusting for PRS (Supplementary Tables 9, 10).

FIGURE 2

Figure 2 Associations of polygenetic risk score (PRS) with risk of health span termination. (A) PRS in men and (B) PRS in women. Participants were divided into five groups according to the quintiles of PRS, and the HRs of each quintile were compared with those in the bottom quintiles. The level of PRS was further divided into three groups: low (quintile 1), intermediate (quintiles 2–4), and high (quintile 5). The blue squares are point estimation of HRs and error bars are 95% CIs. Fully adjusted model: age, menopause status (for women), college or university degree, deprivation index, body mass index, smoking status, alcohol drinking, IPAQ group, healthy diet, family history of CCVD or cancer, use of aspirin/ibuprofen, use of hormone replacement therapy (for women), the top 10 principal components of ancestry, and genotyping chip. HR, hazard ratio; CI, confidential interval.

Then, we found a joint effect of genetic risk and TT/FT. Specifically, compared with men with high genetic risk and low level of TT/FT, a 20% lower risk was observed in those with both low genetic risk and high level of TT/FT (TT, HR: 0.80, 95% CI: 0.73–0.88; FT, HR: 0.80, 95% CI: 0.73–0.87, Figures 3A, B). Besides, the HST risk increased with both genetic risk and testosterone concentration in women. Women with high genetic risk and testosterone level showed the highest risk in subgroups (TT, HR: 1.32, 95% CI: 1.19–1.46; FT, HR: 1.38, 95% CI: 1.25–1.53, Figures 3C, D).

FIGURE 3

Figure 3 Joint analyses of testosterone and genetic categories for risk of health span termination. The joint effects for health span termination were estimated according to TT of men (A), FT of men (B), TT of women (C), FT of women (D), and PRS categories, respectively. The blue squares are point estimation of HRs and error bars are 95% CIs. The level of testosterone was divided into three groups: low (quintile 1), moderate (quintiles 2–4), and high (quintile 5), and the categories of PRS were similar. Each group was compared with the participants with low level of testosterone and high genetic risk in men, while the combination of low level of testosterone and low genetic risk was the reference in women. Fully adjusted model: age, menopause status (for women), college or university degree, deprivation index, body mass index, smoking status, alcohol drinking, IPAQ group, healthy diet, family history of CCVD or cancer, use of aspirin/ibuprofen, use of hormone replacement therapy (for women), the top 10 principal components of ancestry, and genotyping chip. SHBG was additionally adjusted in total testosterone. HR, hazard ratio; CI, confidential interval; Ref., reference; No., number; TT, total testosterone; FT, free testosterone.

Finally, sensitivity analyses showed consistent results after excluding participants with HST during the first 2 years of follow-up (Supplementary Table 11). When we excluded the participants with poor self-reported health status at baseline, the associations were also consistent (Supplementary Table 12). Moreover, after excluding participants with outliers of testosterone level (top 1% and bottom 1%), we still observed similar associations (Supplementary Table 13). After additional adjustment for fasting time or menstrual cycle proxy factors, similar associations were also presented (Supplementary Tables 14, 15).

Discussion

In this large-scale cohort study, comprised of 145,481 men and 147,733 women from the UK Biobank, we observed sex-specific associations between serum testosterone and HST risk. The findings suggested that total testosterone might decrease the risk of health span termination in men. Inversely, women with a high level of TT or FT were at a higher risk of HST. Furthermore, genetic factors increased the risk of premature end of health span to some extent. Although no interplay between testosterone and genetic factors was observed, participants with both high genetic risk and abnormal levels of testosterone (low level in men and high level in women) had the highest HST risk.

Despite the scarcity of direct evidence in the recent two decades relating testosterone and health span, a growing number of epidemiological studies indicated that high levels of testosterone might be protective for all-cause mortality in men (12, 25), but inversely in women (13), which is consistent with our cause-specific findings. Similar to a prospective study from Denmark (10), we observed that the risk of overall cancer for men was not modified by testosterone levels. However, some researchers indicated that high testosterone might be a risk factor for prostate cancer (11). For women, the effects of endogenous testosterone might increase the risk of common gynecological cancer, such as breast (26) and endometrial cancer (27), supporting our cause-specific results. Metabolic syndrome was considered as a critical risk factor for common cancers (28). A high level of testosterone had been proven beneficial for metabolic syndrome risk in men but harmful in women (9, 29, 30), which supported our findings. Additionally, recent researches pointed out that testosterone deficiency in men was associated with increased risk of non-cancer chronic disease, such as cardiovascular disease (8, 31), type 2 diabetes (29, 32) and dementia (33, 34), as well as exacerbated the condition of COPD patients (35). In women, early studies revealed an inverse association of testosterone with increased CVD risk (36, 37), while it was not observed in the recent Rotterdam study (38). Although testosterone may be safe for myocardial function at physiological levels, a toxic effect also exists at supraphysiologic dose (39). As mentioned above, higher levels of testosterone are unfavorable to the metabolism of women and increase their risk of type 2 diabetes (30, 40). Thus, these pieces of evidence support that a sex-specific association indeed exists between testosterone and HST risk.

Except for endocrine factors, genetic variations also play an important role in the determination of individual health span (14). PRS is a simple approach for summing risks across multiple susceptibility loci, aiming at classifying populations with distinct risk levels to drive clinical or personal decision-making (15). Although similar longevity genome-wide association studies (41, 42) exist for a long time, studies that have assessed the predictive performance of PRS on longevity or health span remain scarce. Besides, the small sample size of a former study could affect the strength of the findings (43). In our study, PRS based on 12 SNPs from nearly 300,000 participants stratifies the participants into distinct HST risk strata clearly, indicating GWAS-derived PRS is a promising tool. The effectiveness of PRS still requires independent validation in other populations, as conducted in a large-scale prospective cohort study that applied the PRS to the CKB cohort for re-evaluation (16).

By summarizing the existing evidence related to death, cancer, and other non-cancer events, we indirectly discuss the sex-specific roles of testosterone in maintaining the period of health span. In view of a potential hazard of high concentration in prostate cancer (11) or other diseases in men, the negative association of testosterone with HST we observed in this study might be offset to some degree, accounting for non-significance in quintile 5 of TT. Thus, we conservatively suggest that this association is more stable at an intermediate dose (the second to fourth quintiles, Q2–Q4) in men, which is consistent with the guidelines by the European Endocrine Society that recommended achieving testosterone concentration in the mid-normal range during treatment (44). Previous studies indicated that a loss of testosterone due to aging contributed to the onset of physical frailty, such as triggering muscle loss, muscle weakness, and decreased functional performance, which eventually limited the lifespan in both genders (7, 45). Hence, an appropriate dose of testosterone is not necessarily harmful to women; our results also showed that HST risk only obviously increased at a high concentration of TT or FT (the top quintile, Q5). Besides, individuals with a very high or very low level of testosterone are strongly associated with health status (46, 47). Some diseases, such as pituitary disorders (48), HIV (49), and hyperthyroidism (50), as well as obesity (51), may interfere the testosterone production, increasing the possibility of extreme level of occurrence. It is worth noting that these diseases have a possible influence on personal health span as well. Therefore, we respectively excluded the participants with self-reported poor health status or outliers in the sensitivity analyses to control the potential bias introduced by the extreme values.

Given that both endocrine and genetic factors may contribute to health span collectively, we assess the joint effects of genetic risk and testosterone. Despite the non-significant interaction effect of testosterone and genetic propensity on health span, we still find the combined associations greater than individually. This suggests that those participants both with abnormal testosterone levels (especially high level in women and low level in men) and high genetic risk are in a more dangerous situation and deserve more attention. Therefore, taking appropriate intervention measures early, such as testosterone replacement therapy (TRT) for men with testosterone deficiency or anti-testosterone therapy for women with a supraphysiologic dose of testosterone, is considerable to extend personal health span.

Over the past two decades, TRT has been widely used for low age-related levels of testosterone in men (52) or even for postmenopausal women with low sexual desire (53), so there is a chance for overuse in the population. Therefore, current mainstream guidelines are cautious about such therapy and are critical about the usage and concentration (54, 55). For a better understanding of the roles played by exogenous testosterone in health span and to better assist clinical or personal decision-making before testosterone utilization, it is a wise choice to set endogenous testosterone as an entry point of research. Hopefully soon, we will be able to screen populations suitable for testosterone therapy by utilizing PRS in risk prediction for precise prevention.

Several limitations deserve to be noticed. Firstly, a single measurement for testosterone at baseline may not reflect long-term exposure very well. However, some researchers consider it reliable in categorizing average levels, as it calculates the intraclass correlation coefficients (ICCs) between two measurements of TT that were 4 years apart in a subcohort of UKB (13). Secondly, although the Vermeulen formula is the most widely used among different methods, calculating FT is still possible to be over- or underestimated relative to laboratory measurement (56). Third, given the aromatization of free testosterone into estradiol (57), further adjustment for estradiol is necessary to evaluate the independent association of FT in women. However, it is infeasible to perform at the current phase because most postmenopausal women in the UK Biobank have estradiol levels below the threshold for detection and premenopausal women also have a high missing rate of estradiol. Finally, there is evidence of a “healthy volunteer” selection bias (58), which suggests that it is warranted to replicate our findings in other cohorts to determine the applicability of our findings to other populations.

In conclusion, this prospective cohort study revealed the sex-specific associations that testosterone was negatively associated with HST risk in men and positively associated with HST risk in women. Genetic factors increased the HST risk, which indicated that participants with abnormal testosterone levels (high level in women or low level in men) and high genetic risk should be the target for early prevention. Although our findings highlight the association between testosterone and health span, further mechanistic studies and prospective trials are warranted to explore the causation behind.

Data Availability Statement

Publicly available datasets were analyzed in this study. Datasets are available from the UK Biobank via application to the UK Biobank.

Ethics Statement

The studies involving human participants were reviewed and approved by the North West Multi-Centre Research Ethics Committee (http://www.ukbiobank.ac.uk/ethics/). The patients/participants provided their written informed consent to participate in this study.

Author Contributions

JD supervised the entire project and designed the work. XZ conducted the statistical analysis and wrote the first draft. SL, NW, and TH helped apply for permission to use data and contributed to data management. MS and JF offered statistical support during the study. MZ, CW, DH, and YJ critically revised the manuscript for important intellectual content. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by the National Natural Science Foundation of China (Grant ID: 81941020) as part of a project entitled “Identification of longevity-related genetic variations and their biological mechanisms based on a natural population cohort”.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

This research was conducted using the UK Biobank Resource (Application Number: 64689). We thank the investigators and participants in UK Biobank for their contributions to this study.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fendo.2021.773464/full#supplementary-material

References

1. Oeppen J, Vaupel JW. Demography. Broken Limits to Life Expectancy. Science (2002) 296(5570):1029–31. doi: 10.1126/science.1069675

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Partridge L, Deelen J, Slagboom PE. Facing Up to the Global Challenges of Ageing. Nature (2018) 561(7721):45–56. doi: 10.1038/s41586-018-0457-8

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Andersen SL, Sebastiani P, Dworkis DA, Feldman L, Perls TT. Health Span Approximates Life Span Among Many Supercentenarians: Compression of Morbidity at the Approximate Limit of Life Span. J Gerontol A Biol Sci Med Sci (2012) 67(4):395–405. doi: 10.1093/gerona/glr223

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Abeliansky AL, Erel D, Strulik H. Aging in the USA: Similarities and Disparities Across Time and Space. Sci Rep (2020) 10(1):14309. doi: 10.1038/s41598-020-71269-3

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Raverot V, Lopez J, Grenot C, Pugeat M, Déchaud H. New Approach for Measurement of Non-SHBG-Bound Testosterone in Human Plasma. Anal Chim Acta (2010) 658(1):87–90. doi: 10.1016/j.aca.2009.10.057

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Feldman HA, Longcope C, Derby CA, Johannes CB, Araujo AB, Coviello AD, et al. Age Trends in the Level of Serum Testosterone and Other Hormones in Middle-Aged Men: Longitudinal Results From the Massachusetts Male Aging Study. J Clin Endocrinol Metab (2002) 87(2):589–98. doi: 10.1210/jcem.87.2.8201

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Horstman AM, Dillon EL, Urban RJ, Sheffield-Moore M. The Role of Androgens and Estrogens on Healthy Aging and Longevity. J Gerontol A Biol Sci Med Sci (2012) 67(11):1140–52. doi: 10.1093/gerona/gls068

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Elagizi A, Köhler TS, Lavie CJ. Testosterone and Cardiovascular Health. Mayo Clin Proc (2018) 93(1):83–100. doi: 10.1016/j.mayocp.2017.11.006

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Rao PM, Kelly DM, Jones TH. Testosterone and Insulin Resistance in the Metabolic Syndrome and T2DM in Men. Nat Rev Endocrinol (2013) 9(8):479–93. doi: 10.1038/nrendo.2013.122

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Ørsted DD, Nordestgaard BG, Bojesen SE. Plasma Testosterone in the General Population, Cancer Prognosis and Cancer Risk: A Prospective Cohort Study. Ann Oncol (2014) 25(3):712–8. doi: 10.1093/annonc/mdt590

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Watts EL, Appleby PN, Perez-Cornago A, Bueno-de-Mesquita HB, Chan JM, Chen C, et al. Low Free Testosterone and Prostate Cancer Risk: A Collaborative Analysis of 20 Prospective Studies. Eur Urol (2018) 74(5):585–94. doi: 10.1016/j.eururo.2018.07.024

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Yeap BB, Marriott RJ, Antonio L, Chan YX, Raj S, Dwivedi G, et al. Serum Testosterone is Inversely and Sex Hormone-Binding Globulin Is Directly Associated With All-Cause Mortality in Men. J Clin Endocrinol Metab (2021) 106(2):e625–37. doi: 10.1210/clinem/dgaa743

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Wang J, Fan X, Yang M, Song M, Wang K, Giovannucci E, et al. Sex-Specific Associations of Circulating Testosterone Levels With All-Cause and Cause-Specific Mortality. Eur J Endocrinol (2021) 184(5):723–32. doi: 10.1530/EJE-20-1253

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Zenin A, Tsepilov Y, Sharapov S, Getmantsev E, Menshikov LI, Fedichev PO, et al. Identification of 12 Genetic Loci Associated With Human Healthspan. Commun Biol (2019) 2:41. doi: 10.1038/s42003-019-0290-0

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Torkamani A, Wineinger NE, Topol EJ. The Personal and Clinical Utility of Polygenic Risk Scores. Nat Rev Genet (2018) 19(9):581–90. doi: 10.1038/s41576-018-0018-x

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Dai J, Lv J, Zhu M, Wang Y, Qin N, Ma H, et al. Identification of Risk Loci and a Polygenic Risk Score for Lung Cancer: A Large-Scale Prospective Cohort Study in Chinese Populations. Lancet Respir Med (2019) 7(10):881–91. doi: 10.1016/S2213-2600(19)30144-4

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Maas P, Barrdahl M, Joshi AD, Auer PL, Gaudet MM, Milne RL, et al. Breast Cancer Risk From Modifiable and Nonmodifiable Risk Factors Among White Women in the United States. JAMA Oncol (2016) 2(10):1295–302. doi: 10.1001/jamaoncol.2016.1025

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Palmer LJ. UK Biobank: Bank on It. Lancet (2007) 369(9578):1980–2. doi: 10.1016/S0140-6736(07)60924-6

PubMed Abstract | CrossRef Full Text | Google Scholar

19. UK Biobank. Biomarker Assay Quality Procedures: Approaches Used to Minimise Systematic and Random Errors (and the Wider Epidemiological Implications). Version 1.2 (2019). Available at: http://www.ukbiobank.ac.uk/ (Accessed January 31, 2020).

Google Scholar

20. Fry D AR, Moffat S, Gordon M, Singh P. UK Biobank Biomarker Project. Companion Document to Accompany Serum Biomarker Data. Version 1.0 (2019). Available at: http://www.ukbiobank.ac.uk/ (Accessed January 31, 2020).

Google Scholar

21. Vermeulen A, Verdonck L, Kaufman JM. A Critical Evaluation of Simple Methods for the Estimation of Free Testosterone in Serum. J Clin Endocrinol Metab (1999) 84(10):3666–72. doi: 10.1210/jcem.84.10.6079

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Kirkwood TB. Deciphering Death: A Commentary on Gompertz (1825) 'On the Nature of the Function Expressive of the Law of Human Mortality, and on a New Mode of Determining the Value of Life Contingencies'. Philos Trans R Soc London Ser B Biol Sci (2015) 370(1666):20140379. doi: 10.1098/rstb.2014.0379

CrossRef Full Text | Google Scholar

23. Bycroft C, Freeman C, Petkova D, Band G, Elliott LT, Sharp K, et al. The UK Biobank Resource With Deep Phenotyping and Genomic Data. Nature (2018) 562(7726):203–9. doi: 10.1038/s41586-018-0579-z

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Rutten-Jacobs LC, Larsson SC, Malik R, Rannikmäe K, Sudlow CL, Dichgans M, et al. Genetic Risk, Incident Stroke, and the Benefits of Adhering to a Healthy Lifestyle: Cohort Study of 306 473 UK Biobank Participants. BMJ (2018) 363:k4168. doi: 10.1136/bmj.k4168

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Pye SR, Huhtaniemi IT, Finn JD, Lee DM, O'Neill TW, Tajar A, et al. Late-Onset Hypogonadism and Mortality in Aging Men. J Clin Endocrinol Metab (2014) 99(4):1357–66. doi: 10.1210/jc.2013-2052

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Key TJ, Appleby PN, Reeves GK, Travis RC, Alberg AJ, Barricarte A, et al. Sex Hormones and Risk of Breast Cancer in Premenopausal Women: A Collaborative Reanalysis of Individual Participant Data From Seven Prospective Studies. Lancet Oncol (2013) 14(10):1009–19. doi: 10.1016/S1470-2045(13)70301-2

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Allen NE, Key TJ, Dossus L, Rinaldi S, Cust A, Lukanova A, et al. Endogenous Sex Hormones and Endometrial Cancer Risk in Women in the European Prospective Investigation Into Cancer and Nutrition (EPIC). Endocr Relat Cancer (2008) 15(2):485–97. doi: 10.1677/ERC-07-0064

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Esposito K, Chiodini P, Colao A, Lenzi A, Giugliano D. Metabolic Syndrome and Risk of Cancer: A Systematic Review and Meta-Analysis. Diabetes Care (2012) 35(11):2402–11. doi: 10.2337/dc12-0336

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Bianchi VE, Locatelli V. Testosterone a Key Factor in Gender Related Metabolic Syndrome. Obes Rev (2018) 19(4):557–75. doi: 10.1111/obr.12633

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Ruth KS, Day FR, Tyrrell J, Thompson DJ, Wood AR, Mahajan A, et al. Using Human Genetics to Understand the Disease Impacts of Testosterone in Men and Women. Nat Med (2020) 26(2):252–8. doi: 10.1038/s41591-020-0751-5

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Kloner RA, Carson C 3rd, Dobs A, Kopecky S, Mohler ER 3rd. Testosterone and Cardiovascular Disease. J Am Coll Cardiol (2016) 67(5):545–57. doi: 10.1016/j.jacc.2015.12.005

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Kelly DM, Jones TH. Testosterone: A Metabolic Hormone in Health and Disease. J Endocrinol (2013) 217(3):R25–45. doi: 10.1530/JOE-12-0455

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Lv W, Du N, Liu Y, Fan X, Wang Y, Jia X, et al. Low Testosterone Level and Risk of Alzheimer's Disease in the Elderly Men: A Systematic Review and Meta-Analysis. Mol Neurobiol (2016) 53(4):2679–84. doi: 10.1007/s12035-015-9315-y

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Verdile G, Laws SM, Henley D, Ames D, Bush AI, Ellis KA, et al. Associations Between Gonadotropins, Testosterone and β Amyloid in Men at Risk of Alzheimer's Disease. Mol Psychiatry (2014) 19(1):69–75. doi: 10.1038/mp.2012.147

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Baillargeon J, Urban RJ, Zhang W, Zaiden MF, Javed Z, Sheffield-Moore M, et al. Testosterone Replacement Therapy and Hospitalization Rates in Men With COPD. Chron Respir Dis (2019) 16:1479972318793004. doi: 10.1177/1479972318793004

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Davis SR. Cardiovascular and Cancer Safety of Testosterone in Women. Curr Opin Endocrinol Diabetes Obes (2011) 18(3):198–203. doi: 10.1097/MED.0b013e328344f449

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Sievers C, Klotsche J, Pieper L, Schneider HJ, März W, Wittchen HU, et al. Low Testosterone Levels Predict All-Cause Mortality and Cardiovascular Events in Women: A Prospective Cohort Study in German Primary Care Patients. Eur J Endocrinol (2010) 163(4):699–708. doi: 10.1530/EJE-10-0307

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Meun C, Franco OH, Dhana K, Jaspers L, Muka T, Louwers Y, et al. High Androgens in Postmenopausal Women and the Risk for Atherosclerosis and Cardiovascular Disease: The Rotterdam Study. J Clin Endocrinol Metab (2018) 103(4):1622–30. doi: 10.1210/jc.2017-02421

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Bianchi VE. Testosterone, Myocardial Function, and Mortality. Heart Fail Rev (2018) 23(5):773–88. doi: 10.1007/s10741-018-9721-0

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Karakas M, Schäfer S, Appelbaum S, Ojeda F, Kuulasmaa K, Brückmann B, et al. Testosterone Levels and Type 2 Diabetes-No Correlation With Age, Differential Predictive Value in Men and Women. Biomolecules (2018) 8(3):76. doi: 10.3390/biom8030076

CrossRef Full Text | Google Scholar

41. Walter S, Atzmon G, Demerath EW, Garcia ME, Kaplan RC, Kumari M, et al. A Genome-Wide Association Study of Aging. Neurobiol Aging (2011) 32(11):2109.e2115–2128. doi: 10.1016/j.neurobiolaging.2011.05.026

CrossRef Full Text | Google Scholar

42. Lunetta KL, D'Agostino RBSR, Karasik D, Benjamin EJ, Guo CY, Govindaraju R, et al. Genetic Correlates of Longevity and Selected Age-Related Phenotypes: A Genome-Wide Association Study in the Framingham Study. BMC Med Genet (2007) 8 Suppl 1(Suppl 1):S13. doi: 10.1186/1471-2350-8-S1-S13

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Zeng Y, Nie C, Min J, Chen H, Liu X, Ye R, et al. Sex Differences in Genetic Associations With Longevity. JAMA Netw Open (2018) 1(4):e181670. doi: 10.1001/jamanetworkopen.2018.1670

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Bhasin S, Brito JP, Cunningham GR, Hayes FJ, Hodis HN, Matsumoto AM, et al. Testosterone Therapy in Men With Hypogonadism: An Endocrine Society Clinical Practice Guideline. J Clin Endocrinol Metab (2018) 103(5):1715–44. doi: 10.1210/jc.2018-00229

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Saad F, Röhrig G, von Haehling S, Traish A. Testosterone Deficiency and Testosterone Treatment in Older Men. Gerontology (2017) 63(2):144–56. doi: 10.1159/000452499

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Mailliez A, Guilbaud A, Puisieux F, Dauchet L, Boulanger É. Circulating Biomarkers Characterizing Physical Frailty: CRP, Hemoglobin, Albumin, 25OHD and Free Testosterone as Best Biomarkers. Results of a Meta-Analysis. Exp Gerontol (2020) 139:111014. doi: 10.1016/j.exger.2020.111014

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Basaria S. Male Hypogonadism. Lancet (London England) (2014) 383(9924):1250–63. doi: 10.1016/S0140-6736(13)61126-5

CrossRef Full Text | Google Scholar

48. Tajar A, Forti G, O'Neill TW, Lee DM, Silman AJ, Finn JD, et al. Characteristics of Secondary, Primary, and Compensated Hypogonadism in Aging Men: Evidence From the European Male Ageing Study. J Clin Endocrinol Metab (2010) 95(4):1810–8. doi: 10.1210/jc.2009-1796

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Pezzaioli LC, Quiros-Roldan E, Paghera S, Porcelli T, Maffezzoni F, Delbarba A, et al. The Importance of SHBG and Calculated Free Testosterone for the Diagnosis of Symptomatic Hypogonadism in HIV-Infected Men: A Single-Centre Real-Life Experience. Infection (2021) 49(2):295–303. doi: 10.1007/s15010-020-01558-6

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Crawford M, Kennedy L. Testosterone Replacement Therapy: Role of Pituitary and Thyroid in Diagnosis and Treatment. Trans Andrology Urol (2016) 5(6):850–8. doi: 10.21037/tau.2016.09.01

CrossRef Full Text | Google Scholar

51. Kelly DM, Jones TH. Testosterone and Obesity. Obes Rev (2015) 16(7):581–606. doi: 10.1111/obr.12282

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Baillargeon J, Kuo YF, Westra JR, Urban RJ, Goodwin JS. Testosterone Prescribing in the United States, 2002-2016. JAMA (2018) 320(2):200–2. doi: 10.1001/jama.2018.7999

PubMed Abstract | CrossRef Full Text | Google Scholar

53. Islam RM, Bell RJ, Green S, Page MJ, Davis SR. Safety and Efficacy of Testosterone for Women: A Systematic Review and Meta-Analysis of Randomised Controlled Trial Data. Lancet Diabetes Endocrinol (2019) 7(10):754–66. doi: 10.1016/S2213-8587(19)30189-5

PubMed Abstract | CrossRef Full Text | Google Scholar

54. Salter CA, Mulhall JP. Guideline of Guidelines: Testosterone Therapy for Testosterone Deficiency. BJU Int (2019) 124(5):722–9. doi: 10.1111/bju.14899

PubMed Abstract | CrossRef Full Text | Google Scholar

55. Qaseem A, Horwitch CA, Vijan S, Etxeandia-Ikobaltzeta I, Kansagara D. Testosterone Treatment in Adult Men With Age-Related Low Testosterone: A Clinical Guideline From the American College of Physicians. Ann Intern Med (2020) 172(2):126–33. doi: 10.7326/M19-0882

PubMed Abstract | CrossRef Full Text | Google Scholar

56. Ly LP, Sartorius G, Hull L, Leung A, Swerdloff RS, Wang C, et al. Accuracy of Calculated Free Testosterone Formulae in Men. Clin Endocrinol (Oxf) (2010) 73(3):382–8. doi: 10.1111/j.1365-2265.2010.03804.x

PubMed Abstract | CrossRef Full Text | Google Scholar

57. Simpson ER. Aromatization of Androgens in Women: Current Concepts and Findings. Fertil Steril (2002) 77(Suppl 4):S6–10. doi: 10.1016/S0015-0282(02)02984-9

PubMed Abstract | CrossRef Full Text | Google Scholar

58. Fry A, Littlejohns TJ, Sudlow C, Doherty N, Adamska L, Sprosen T, et al. Comparison of Sociodemographic and Health-Related Characteristics of UK Biobank Participants With Those of the General Population. Am J Epidemiol (2017) 186(9):1026–34. doi: 10.1093/aje/kwx246

PubMed Abstract | CrossRef Full Text | Google Scholar