|
|
||||||||
EXPERIMENTAL STUDY |
Department of Andrology, Concord Hospital and ANZAC Research Institute, University of Sydney, Sydney NSW 2139, Australia
(Correspondence should be addressed to D J Handelsman; Email: djh{at}anzac.edu.au)
| Abstract |
|---|
|
|
|---|
Methods: We analyzed a large dataset comprising over 4000 consecutive blood samples in which FT as well as TT and SHBG were measured. Dividing the dataset into samples with blood TT above and below 5 nM, using a bootstrap regression modeling approach guided by Akaike Information Criterion for model selection to balance parsimony against reduction of residual error, empirical equations were developed for FT in terms of TT and SHBG.
Results: Comparison between the empirical FT equations with the laboratory FT measurements as well as three widely used calculated FT methods showed the empirical FT formulae had superior fidelity with laboratory measurements while previous FT formulae overestimated and deviated systematically from the laboratory FT values.
Conclusion: We conclude that these simple, assumption-free empirical FT equations can estimate accurately blood FT from TT and SHBG measured in the same samples with the present assay methods and have suitable properties for wider application to evaluate the clinical utility of blood FT measurements.
| Introduction |
|---|
|
|
|---|
|
| Materials and methods |
|---|
|
|
|---|
Data including age and sex from all blood samples submitted between 1999 and 2003 to the Central Sydney Area Health Service Endocrinology Diagnostic Laboratory requesting testosterone measurements were routinely assayed for FT by centrifugal ultrafiltration as well as for TT and SHBG. All laboratory data for this study were provided in de-identified form after discharge of all diagnostic reporting responsibilities so that no ethics approval was required.
Assays
Hormone assays were performed in a single teaching hospital-based routine diagnostic endocrinology laboratory as described previously (2022) using established commercial immunoassays routinely monitored by participation in external quality-control programs. Plasma TT and SHBG were measured by commercial immunoassays (Immulyte, Los Angeles, CA, USA). During the period of this study, coefficients of variation for TT assay (n = 3444 assays) were 8.2% at high (mean, 28.5 nM), 12.4% at mid-range (13.6 nM) and 30.2% at low (2.9 nM) concentrations and for SHBG assay (n = 104 assays) were 5.6% at high (71.4 nM) and 5.6% at low (5.1 nM) concentrations.
FT was measured by an in-house adaptation of the centrifugal ultrafiltration assay (13). Samples and controls (600 µl) were well mixed with 50 µl tritiated testosterone tracer in glass tubes and incubated for 1.5 h in a 37 °C waterbath. Specimens (500 µl) were then added to Centrifree columns (Millipore) and centrifuged at 2000 r.p.m. for 20 min (30 °C) following which 50 µl specimens of filtrate and totals were counted in a liquid scintillation counter. The proportion of unbound testosterone was then calculated and the actual FT calculated from the TT in the same sample. Coefficients of variation for the proportion of unbound testosterone were 9.4% at low (mean, 1.2% unbound) and 12.8% at high (2.2% unbound) quality-control samples.
Data analysis
FT was calculated by published methods described by Sodergard et al. (16) and Vermeulen et al. (2) based on equilibrium binding equations using TT and SHBG results from the same blood sample. Free androgen index (FAI) was calculated as described in (17).
Regression of laboratory or calculated FT on other variables was performed using SPSS version 12 software. Bootstrap resampling (23) was performed with S Plus software. Agreement between estimates and identification of systematic discrepancies between estimates of FT was analysed by deviance plots, modified from the BlandAltman approach (24), with calculation of mean bias and limits of agreement. The deviance plot places the difference between the test and standard method on the y axis whereas the x axis is the laboratory measurement, rather than the mean of test and standard as in the original BlandAltman approach. This enhances interpretation of deviations according to the true result as the laboratory measurement is the recognized gold standard. This adaptation is valid as the correlation between methods is very high (>0.95), thereby avoiding problems of interpretation due to spurious correlation where correlations are lower (25). For evaluation of FAI, a dimensionless ratio, a quantile/quantile plot using standardized normal deviates (Z scores) was used.
Goodness of fit for competing regression models was evaluated by AIC value (26), a maximum-likelihood approach to global model entropy reduction in balancing the number of parameters against the reduction in residual error. This was defined as
![]() |
where
2 =
[(calculated FT laboratory FT)2/laboratory FT] and df (degrees of freedom) = number of cases 1.
| Results |
|---|
|
|
|---|
|
After preliminary analysis indicating the need for two formulae to span the full Range covered by these specimens, data were divided into samples with TT above or below 5 nM for all further analyses (Table 1
). Prediction equations for FT from TT and SHBG considered as main effects as well as including an interaction and quadratic terms for main effects were estimated according to the following general formula:
![]() |
The following hierarchy models was considered:
![]() | (Model 1) |
![]() | (Model 2) |
![]() | (Model 3) |
![]() | (Model 4) |
Full data regression models
Using the full dataset (n = 3975), in the low-TT range model 1 was the best (AIC, 5147) followed by model 2 (6804), model 3 (8308) and model 4 (13 532). In the high-TT range, the sequence was reversed with model 4 the best (AIC, 122 247), followed by model 3 (129 146), model 2 (136 045) and model 1 (142 945). As it was not possible to evaluate the quality of fit of these equations within the same dataset, a modified split-sampling strategy was adopted for subsequent model development and evaluation against laboratory FT and other calculated FT formulae.
Bootstrap resampling regression models
A limitation of using the full dataset for a regression model is that there is no independent mechanism to test the adequacy of its fit in different datasets. For this reason, rather than using a standard split-sampling approach (allowing only a single test), we proceeded to develop a bootstrap resampling methodology using 1000 replications of a random 60% selection (with replacement) from the full dataset.
From the bootstrap models, for low testosterone (TT <5 nM) model 1 was the most favorable followed by model 3, model 2 and model 4 (AIC values of 4874, 4937, 5032 and 6441, respectively). For high testosterone (TT
5 nM), model 2 was preferred followed by model 4, model 1 and model 3 (AIC values of 142 214, 142 359, 177 512 and 178 684, respectively).
The empirical FT (EFT) formulae for further evaluation against other methods had the following coefficients (with TT and SHBG in nM) as specified below:
![]() |
![]() |
Comparison of empirical equations with laboratory FT measurements and other calculated FT values
The empirical FT formulae were then evaluated for agreement and deviation from the laboratory FT measurements and the two other calculated FT estimates (Sodergard et al. (16) and Vermeulen et al. (2)). The distribution of values by centiles is displayed in Table 2
(for low TT) and Table 3
(for high TT). For the empirical FT formulae, the deviance plots (Fig. 2
) show good agreement with the laboratory FT measurements across the full range of specimens. The Sodergard- and Vermeulen-calculated FT formulae show over-estimation and wider limits of agreement compared with the laboratory FT values.
|
|
|
In order to compare FAI, a dimensionless ratio, to laboratory FT measurements (in concentration units), both were transformed into standardized normal deviates (Z scores) and the discrepancy between Z scores (ZFAIZLabFT) displayed in a deviance plot (Fig. 3
). These plots showed severe influence of SHBG on FAI, especially when SHBG was low, as well as poor agreement with laboratory FT through much of the range.
|
| Discussion |
|---|
|
|
|---|
The present study developed and undertook internal validation of a simple empirical, model-free method to calculate FT, requiring only TT and SHBG concentrations in the same sample. Our simplified approach recognizes the non-linearity of the relationship between FT with TT and SHBG as predicted by well-known theoretical equilibrium binding equations (2, 16) while also freeing it from specific assumptions and simplifications inherent in the equilibrium binding equations model. Preliminary modeling showed that no single equation could cover the full range of TT and SHBG concentrations if samples included not just men but also women and children, as well as severely hypogonadal men. However, dividing the samples with an arbitrary threshold of 5 nM allows for development of a pair of excellent, versatile prediction equations. This threshold corresponds to mostly male samples above the threshold as well as females, children or hypogonadal male samples below but was more practical than using sex or age as classification criteria. The inability of a single empirical equation to cover all TT and SHBG concentrations raises questions also about the validity of the theoretical equilibrium binding equation model used in other calculational approaches (2, 16), especially when applied to extreme blood TT and SHBG concentrations. Such invalidity is most likely to stem from the many assumptions required to practically implement equilibrium binding equations rather than any challenge to the Law of Mass Action. The dichotomy between high and low TT levels is reminiscent of the FAI which is invalid for both theoretical reasons as well as empirical evidence when applied to male samples whereas it is valid and useful in women or children who have much lower blood TT (17). Another major contributing factor to the requirement for a distinct model for samples with lower TT concentrations is the recently established invalidity of automated platform testosterone assays for samples in that low range (29, 30), a failing that has been likened to random-number generation (31).
In evaluating our empirical approach, we also reviewed the validity of the three most frequently used calculational estimates of FT against the centrifugal ultrafiltration reference method. While the empirical equation has satisfactory agreement relative to the laboratory FT measurements, the three other widely used calculated FT estimates deviate systematically from the laboratory FT values. The two calculated FT methods based on theoretical binding equations that result in second-degree equations in TT and SHBG both overestimated FT levels throughout the range. The discrepancy between these findings and the previous calculational FT methods may be because the present approach avoids assumptions inherent in them of a nominal, fixed affinity constant for testosterone binding to SHBG, an assumed albumin concentration and affinity and neglecting other blood testosterone-binding proteins. A recent study noted flaws in the suitability of the two calculated FT equations attributable to the notional SHBG-binding affinity and other assumptions (32). It is also notable that the original validation studies for these calculated FT methods were based on very few samples, namely 11 (16) and 28 (2). Another comparative study involved 50 samples (3). On the other hand an important caveat on the present empirical approach is its reliance on the specific TT and SHBG assays used in this study. It cannot be assumed that the specific equations would be portable to other assay combinations, or even the same assays if they were significantly modified.
The third calculated measure, the so-called free androgen (or testosterone) index (FAI), is the ratio of TT/SHBG usually expressed as a percentage. The FAI correlated poorly with laboratory FT measurements and demonstrated extreme influence of low SHBG concentrations. This is not surprising as this ratio over-simplifies and ignores the non-linearity of FT in terms of TT and SHBG as predicted by theoretical binding equations, whereby FT may be approximated by this ratio if, and only if, the TT concentration is negligible in relation to the concentrations of SHBG-binding sites (17). While this is a reasonable approximation for samples from women and children where blood testosterone concentrations are rarely above 10% of SHBG concentration, this is not true in samples from most men in whom the higher blood TT concentrations are comparable with the SHBG concentrations and cannot be neglected. Thus our more extensive evaluation confirms the previous data that FAI does not correspond to actual FT measurements (17).
The free hormone hypothesis (810, 3335) remains unproven (36) and requires further external validation. In asserting that the non-protein-bound free fraction is the most biologically active moiety of a circulating steroid hormone with the protein-bound moiety a reserve, biologically inactive buffer, this concept lacks theoretical validity or empirical proof. In theoretical terms, if non-protein-bound (or lightly bound) circulating steroids are more readily transported to tissues, this applies equally to target tissues (enhancing bioactivity) as well as to hepatic sites of steroid degradation (terminating bioactivity). The net balance between these two countervailing effects is inherently unpredictable depending on many dynamic factors including relative tissue mass and blood flow of target and metabolizing tissues. Hence, whether free hormone measurements represent a more active or less (more rapidly inactivated) biologically active moiety of a circulating steroid cannot be assumed a priori. The demonstrations that SHBG-bound testosterone is biologically active, via binding to cell-surface SHBG receptor (37), and that rodents have no circulating SHBG (38), further question the free hormone hypothesis which predicts that testosterone tightly bound to SHBG would form a biologically inactive buffer reservoir. Yet, despite deficient theoretical rationale, derived testosterone measures might still be useful empirically if they provided improved prediction over standard blood TT measurements to identify rectifiable androgen deficiency. The conspicuous lack of such empirical validation may be at least partly due to the limited availability of valid calculated FT methods. The present findings are not directly informative regarding other derived testosterone measures that involve an additional assay step (e.g. free analog assay or bioavailable testosterone). While both assays may be technically reproducible (2, 3, 39), the free analog assay clearly does not measure free testosterone (18, 19, 40) while the bioavailable (corresponding to free plus loosely albumin-bound fractions) is not intended to correspond to free testosterone measurements. For both assays, empirical validation relies on convincing evidence of providing superior or additional biological information to TT measurements as judged by independent biological effects. This remains contentious and almost entirely unproven.
We conclude that a simple, empirical calculation can provide a robust, reliable and valid estimate of blood FT from TT and SHBG in the same sample. This calculated FT method has improved fidelity with empirical laboratory measurements of FT by a reference method and obviates the need for additional laborious manual assays. This study also highlights the limited validation and systematic deviations from reference methods of the three most widely used calculated FT methods. Further studies using the new formulae developed in this study may be helpful in determining whether FT measurements provide any additional clinical value over the well-established clinical gold standard of blood TT measurements to confirm the clinical diagnosis of treatable androgen deficiency.
| Acknowledgements |
|---|
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
S A Paul Chubb, Z. Hyde, O. P Almeida, L. Flicker, P. E Norman, K. Jamrozik, G. J Hankey, and B. B Yeap Lower sex hormone-binding globulin is more strongly associated with metabolic syndrome than lower total testosterone in older men: the Health in Men Study Eur. J. Endocrinol., June 1, 2008; 158(6): 785 - 792. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. J. Handelsman Update in Andrology J. Clin. Endocrinol. Metab., December 1, 2007; 92(12): 4505 - 4511. [Abstract] [Full Text] [PDF] |
||||
![]() |
F. Giton, S. Urien, C. Born, J. Tichet, J. Guechot, J. Callebert, F. Bronsard, J. P. Raynaud, and J. Fiet Determination of Bioavailable Testosterone [Non Sex Hormone Binding Globulin (SHBG)-Bound Testosterone] in a Population of Healthy French Men: Influence of Androstenediol on Testosterone Binding to SHBG Clin. Chem., December 1, 2007; 53(12): 2160 - 2168. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. Y. Liu, J. Beilin, C. Meier, T. V. Nguyen, J. R. Center, P. J. Leedman, M. J. Seibel, J. A. Eisman, and D. J. Handelsman Age-Related Changes in Serum Testosterone and Sex Hormone Binding Globulin in Australian Men: Longitudinal Analyses of Two Geographically Separate Regional Cohorts J. Clin. Endocrinol. Metab., September 1, 2007; 92(9): 3599 - 3603. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. B. Araujo, V. Kupelian, S. T. Page, D. J. Handelsman, W. J. Bremner, and J. B. McKinlay Sex Steroids and All-Cause and Cause-Specific Mortality in Men Arch Intern Med, June 25, 2007; 167(12): 1252 - 1260. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. B Yeap, O. P Almeida, Z. Hyde, P. E Norman, S A P. Chubb, K. Jamrozik, and L. Flicker In men older than 70 years, total testosterone remains stable while free testosterone declines with age. The Health in Men Study Eur. J. Endocrinol., May 1, 2007; 156(5): 585 - 594. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. Ferk, N. Teran, and K. Gersak The (TAAAA)n microsatellite polymorphism in the SHBG gene influences serum SHBG levels in women with polycystic ovary syndrome Hum. Reprod., April 1, 2007; 22(4): 1031 - 1036. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. Dechaud, A. Denuziere, S. Rinaldi, J. Bocquet, H. Lejeune, and M. Pugeat Age-Associated Discrepancy between Measured and Calculated Bioavailable Testosterone in Men Clin. Chem., April 1, 2007; 53(4): 723 - 728. [Abstract] [Full Text] [PDF] |
||||
![]() |
W. Rosner, R. J. Auchus, R. Azziz, P. M. Sluss, and H. Raff Utility, Limitations, and Pitfalls in Measuring Testosterone: An Endocrine Society Position Statement J. Clin. Endocrinol. Metab., February 1, 2007; 92(2): 405 - 413. [Abstract] [Full Text] [PDF] |
||||
![]() |
W. de Ronde, Y. T. van der Schouw, H. A.P. Pols, L. J.G. Gooren, M. Muller, D. E. Grobbee, and F. H. de Jong Calculation of Bioavailable and Free Testosterone in Men: A Comparison of 5 Published Algorithms Clin. Chem., September 1, 2006; 52(9): 1777 - 1784. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. Jakobsson, L. Ekstrom, N. Inotsume, M. Garle, M. Lorentzon, C. Ohlsson, H.-K. Roh, K. Carlstrom, and A. Rane Large Differences in Testosterone Excretion in Korean and Swedish Men Are Strongly Associated with a UDP-Glucuronosyl Transferase 2B17 Polymorphism J. Clin. Endocrinol. Metab., February 1, 2006; 91(2): 687 - 693. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. N. Zhuang, L. P. Ly, R. G. Cumming, and D. J. Handelsman Growth and Development during Early Manhood as Determinants of Prostate Size in Later Life J. Clin. Endocrinol. Metab., November 1, 2005; 90(11): 6055 - 6063. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |