Abstract
We hypothesized that soy isoflavones would attenuate the anticipated increase in androidal fat mass in postmenopausal women during the 36-month treatment, and thereby favorably modify the circulating cardiometabolic risk factors: triacylglycerol, LDL-C, HDL-C, glucose, insulin, uric acid, C-reactive protein, fibrinogen, and homocysteine. We collected data on 224 healthy postmenopausal women at risk for osteoporosis (45.8–65 y, median BMI 24.5) who consumed placebo or soy isoflavones (80 or 120 mg/d) for 36 months and used longitudinal analysis to examine the contribution of isoflavone treatment, androidal fat mass, other biologic factors, and dietary quality to cardiometabolic outcomes. Except for homocysteine, each cardiometabolic outcome model was significant (overall -values from ≤.0001 to .0028). Androidal fat mass was typically the strongest covariate in each model. Isoflavone treatment did not influence any of the outcomes. Thus, androidal fat mass, but not isoflavonetreatment, is likely to alter the cardiometabolic profile in healthy postmenopausal women.
1. Introduction
Menopause is associated with an increase in intra-abdominal fat [1], which is considered a major risk factor of atherosclerotic cardiovascular disease (CVD) [2]. However, whether an increased risk of CVD in postmenopausal women is due to altered body composition, changes in reproductive hormones, or some other physiological process associated with menopause has not been clearly established. It is also uncertain as to what extent androidal fat mass may influence CVD risk factors in healthy nonobese women. The risk of chronic disease is perceived to be considerably higher in obese compared to normal weight adults. However, Gautier et al. [3] recently reported that the effect of waist circumference, reflecting androidal fat mass, in increasing the risk for diabetes was more profound among men and women with a normal body mass index (BMI) of less than 25 kg/m2 than among those with above normal BMI. Thus, it appears that efforts to reduce abdominal fat accumulation would be beneficial regardless of BMI category.
This study was ancillary to the Soy Isoflavones for Reducing Bone Loss (SIRBL) study, a randomized, double-blind, placebo-controlled multicenter (Iowa State University (ISU) and University of California at Davis (UCD)) clinical trial funded by the National Institutes of Health (NIH) [4]. We examined contributors (androidal fat mass, duration of menopause, isoflavone treatment, family history of CVD, and dietary quality) to cardiometabolic risk in primarily normal and overweight healthy postmenopausal women. Cardiometabolic outcomes included circulating low-density lipoprotein cholesterol (LDL-C), high-density lipoprotein cholesterol (HDL-C), triacylglycerol, glucose, homeostatic model assessment (HOMA) of insulin resistance, uric acid, C-reactive protein (CRP), fibrinogen, and homocysteine (Hcy). We hypothesized that consumption of soy isoflavones in tablet form would attenuate the anticipated increase in central fat (androidal) mass during the 36 months of treatment, which in turn may favorably modify circulating concentrations of cardiometabolic risk factors. We also hypothesized that androidal fat mass, other biologic factors, and dietary quality would influence these cardiometabolic risk factors, but that androidal fat mass would predominate.
2. Materials and Methods
2.1. Overall Study Design
The parent study examined the effect of two doses (80 versus 120 mg/day) of soy protein-derived isoflavone versus placebo tablets for 36 months on bone loss in healthy postmenopausal women (45 to 65 years of age) who were at risk for osteoporosis. The parent study included power analysis on the primary outcome (lumbar spine bone mineral density) in the methodology and has been previously reported [4]. This ancillary project focused on the relationship between body composition and cardiometabolic risk factors. Many of these cardiometabolic risk factors were tested throughout the study (LDL-C, HDL-C, TAG, glucose, and uric acid) whereas some were only tested through the 12 month time point (HOMA, CRP, fibrinogen, and Hcy). The respective Institutional Review Boards (IRB) at ISU (ID no. 02-199) and at UCD (ID no. 200210884-2) approved our study protocol, consent form, and all participant-related materials. Approvals for the dual-energy X-ray absorptiometry (DXA) procedures were obtained from each institution’s IRB and State Department of Public Health in Iowa and California. At prebaseline, each participant received a detailed explanation of the study verbally and in writing before signing an informed consent form.
2.2. Participant Recruitment, Screening, and Selection
We recruited (2003 to 2005) women throughout the state of Iowa and in the greater Sacramento and Bay Area regions in northern California primarily through direct mailing lists, stories in local newspapers, and local/regional radio advertisements. Responders () were initially screened via telephone to identify healthy women (without diseases or conditions, not taking hormones or medications) ≤65 years who had undergone natural menopause (cessation of menses 1 to 10 years), were not experiencing excessive vasomotor symptoms, nonsmokers (not currently smoking and who had not smoked in the past 6 months), and had a BMI between 18.5 through 29.9 (except for 9 women from UCD who did not meet the range of inclusion criterion: 8 women had BMI values that ranged from 30 to 32.7 and one woman had a BMI of 17.8, but were enrolled because they were deemed healthy).
The parent SIRBL project established the inclusion/exclusion criteria. We excluded vegans and high alcohol consumers (>7 servings/week), as well as those who were diagnosed with chronic disease, had a first-degree relative with breast cancer, or who chronically used medication (current: cholesterol-lowering and/or antihypertensive; past 3 months: antibiotics; past 6 months: estrogen/progestogen creams, calcitonin; past 12 months: oral hormones/estrogen or selective estrogen receptor modulators; ever: bisphosphonates).
Women who met the initial screening criteria () were invited to the clinic for further eligibility screening, including BMD assessment using DXA. The SIRBL project focused on disease prevention rather than treatment; thus, women with BMD lumbar spine (L1-L4) and/or proximal femur T-scores that were low (>1.5 SD below young adult mean) or high (>1.0 SD above mean) or with evidence of previous or existing spinal fractures were excluded. Once each woman qualified based on BMD, fasting blood was drawn for a clinical chemistry profile. We excluded women with evidence of diabetes mellitus (fasted glucose ≥6.93 mmol/L (126 mg/dl)), abnormal renal, liver (elevated enzymes), and/or thyroid function, or elevated lipids (LDL-cholesterol >4.10 mmol/L (160 mg/dl); triacylglycerol >2.25 mmol/L (200 mg/dl)). Based upon our entry criteria, we randomized 255 women to treatment in the parent trial. We excluded 13 women at UCD from this analysis because they did not meet the entry criteria (11 had thickened endometrium, 1 had breast cancer, 1 could not provide a baseline blood sample). We excluded an additional 19 women who did not have body composition data at either 12, 24, or 36 months because they dropped out of the study, resulting in a sample size of 224 women based upon androidal fat as our primary covariate of interest in this ancillary project.
2.3. Randomization to Treatment and Tablet Formulation
To meet the objectives of the parent project, participants at each location (ISU, UCD) were stratified according to baseline proximal femur BMD (high, medium, and low) [4] based upon NHANES III database population values [5] and randomly assigned to one of three treatment groups: placebo control, 80 mg isoflavones, or 120 mg isoflavones. Tablets were provided by The Archer Daniels Midland Co. (Decatur, IL); tablet composition has been described previously [4]. An independent researcher (Patricia Murphy) at ISU confirmed that the actual isoflavone doses (mean ± SD, mg/d) were similar to those formulated and tested by Archer Daniels Midland, respectively: control = 0 compared with ; 80 mg = compared with ; 120 mg = compared with . Participants in each group were instructed to take three compressed tablets/d. To preserve the double-blind nature of the study, bottles did not indicate treatment assignment.
2.4. Body Size and Composition Measures and Blood Sample Measures
For this study, we used anthropometric and body composition measurements that included weight, height, whole body lean and fat mass, and androidal fat mass, as well as fasting concentrations of lipids/lipoproteins, glucose, and uric acid. These outcomes were assessed at baseline, 12, 24, and 36 months. In addition, circulating insulin, HOMA (calculated as fasting glucose (mg/dL) × fasting insulin (μU/mL)/405), CRP, fibrinogen, Hcy, and red blood cell (RBC) folate concentrations were assessed at baseline and 12 months.
Trained researchers obtained anthropometric measurements according to standard protocols. Standing height was taken twice (average value recorded) with a wall-mounted stadiometer (Model S100; Ayrton Corp., Prior Lake, MN) and weight was measured at ISU using a balance beam scale (ABCO Health-o-Meter; Bridgeview, IL) and at UCD using an electronic scale (Circuits and Systems Inc; E. Rockaway, NY). Women wore hospital scrubs or shorts and a t-shirt, and removed their shoes, belts, watches, and jewelry for the duration of assessment. To ensure standardized data collection, body composition measurements were obtained by certified cross-trained DXA operators using matching DXA instruments (Delphi W Hologic Inc; Bedford, MA) at each site that were calibrated daily. To further ensure quality control, one operator assessed overall composition from the whole body DXA scans for both sites. Regional adiposity analysis was performed by one analyzer (LR) who sectioned each whole body DXA scan into waist, hip, and thigh regions based on bone landmarks [6] using special software (Discovery Version 12.3:7). The waist region included the first lumbar through the fourth lumbar vertebrae. The hip region began below the fourth lumbar vertebrae and extended to the tip of the greater trochanter of the femur. Androidal fat mass (kg) for each participant was the sum of waist and hip fat mass.
Phlebotomists collected fasted (9 h) blood samples between 7:00 and 8:00 am. We separated serum (allowed to clot for 30 min prior to centrifugation) and plasma from whole blood and centrifuged samples for 15 minutes (4°C) at 1000 × g, storing aliquots at −80°C until analyses. Certified clinical laboratories (LabCorp; Kansas City, KS for ISU and UCD Medical Center; Sacramento, CA for UCD) performed a chemistry panel (including serum glucose, lipid profile, and uric acid) on each participant at each time point. We measured the remaining analytes (serum insulin, serum CRP, plasma fibrinogen, plasma Hcy, and RBC folate) for each participant in batch from ISU and UCD samples at baseline and 12 months in duplicate at ISU. We used sufficient in-house sera/plasma as quality-control samples (frozen at −80°C) to run with each kit to calculate interassay coefficient of variation (CV); we used duplicate samples to calculate intraassay CV. The low-to-normal and normal-to-high controls for each kit were well within the acceptable ranges. Serum insulin (μU/mL) concentration was determined with a radioimmunoassay kit (Linco Research, St Charles, MO, USA) using a Cobra II series autogamma counting system (Packard Instrument Company; Meriden, CT, USA). The intra- and interassay CV for insulin were 3.0% and 4.0%, respectively. Serum CRP (mg/L) concentration was determined with a high-sensitivity sandwich enzyme-linked immunosorbent assay kit (ALPCO Diagnostics; Salem, NH) and plasma (heparinized) fibrinogen (mg/mL) concentration was determined with a sandwich enzyme-linked immunosorbent assay kit (AssayPro; St. Charles, MO) using a microtiter plate reader (ELx808; Bio-Tek Instruments, Inc., Winooski, VT). The intra-assay CVs for CRP and fibrinogen were 3.7% and 2.7%, respectively; the interassay CVs for CRP and fibrinogen were 6.0% and 2.3%, respectively. Total Hcy (μmol/L) concentration was determined using a high-performance liquid chromatography (HPLC) method adapted from Araki and Sako [7] and Ubbink et al. [8]. The total Hcy in plasma consists of free Hcy (i.e., reduced plus oxidized Hcy in the nonprotein fraction of plasma) and protein-bound Hcy [9]. N-Acetylcysteine (1 mM) was added as an internal standard to the plasma samples prior to derivatization. The fluorescence intensities were measured with excitation at 385 nm and emission at 515 nm, using a JASCO FP-1520 fluorescence detector. Further assay details have been previously published [9]. The intra- and interassay CV for plasma Hcy were 3.8% and 6.3%, respectively. Intracellular (RBC) folate (ng/mL) was measured using a radioactive immunoassay kit (MP Biomedicals; Irvine, CA). Because four hematocrit values were missing, RBC folate could not be calculated for these four samples. In addition, three baseline samples were incorrectly processed; thus, RBC folate values have been presented for 217 women. The intra- and interassay CV for RBC folate were 3.2% and 11.8%, respectively.
2.5. Interviewer-Administered Questionnaires
During the enrollment phase, trained interviewers administered three questionnaires to participants: a health and medical history [10–12], a reproductive history [13], and a nutrition history [10, 11]. We also gathered data on prescription and over-the-counter medications at each time point, as well as previous and/or current use of herbal therapies or dietary supplements (which they were asked to discontinue prior to baseline testing). We assessed dietary intake at each time point using a semiquantitative food frequency questionnaire from Block Dietary Data Systems (Berkeley, CA). A Healthy Eating Index (HEI) score, with higher scores representing greater adherence to federal dietary guidelines, was calculated for each participant based on this questionnaire and included in the statistical analyses as a covariate.
2.6. Statistical Analysis
We performed statistical analyses using version 2.10.1 of the R software, including version 3.1.96 of the nlme package software and considered results statistically significant (two sided) at . Our 3-year longitudinal analyses of LDL-C, HDL-C, triacylglycerol, and glucose included women with complete data at all time points, baseline through 36 months. Our 1-year longitudinal analyses of fibrinogen, insulin, HOMA, Hcy, and CRP included women with complete data at baseline and 12 months. We reported descriptive statistics for 224 women using median and interquartile range. We constructed longitudinal models to identify significant contributors to each cardiometabolic risk factor (LDL-C, HDL-C, triacylglycerol, uric acid, glucose, insulin, HOMA, CRP, fibrinogen, and Hcy). Each final longitudinal model included these obligatory variables: treatment (control versus combined treatment with 80 mg or 120 mg of isoflavones or, in other words, no treatment versus treatment), time point (baseline versus 12, 24, or 36 months), treatment by time point interaction, and site (ISU versus UCD), as well as potential covariates that included androidal fat mass (kg) adjusted for height, time since last menstrual period (TLMP) (yr) (calculated for each woman: baseline test date—date of her last menstrual period), family history of CVD coded as a categorical variable (none versus positive or none versus unknown), and HEI score. Additionally, the insulin model included glucose as a covariate, and the Hcy model included RBC folate as a covariate. Independent variables in modeling the outcomes of interest included those variables that were biologically plausible. Preliminary models also included dietary fat intake (determined using the Block Food Frequency Questionnaire) and physical activity (determined using the Paffenbarger physical activity recall [14]), but those variables did not emerge as remotely significant contributors to any of the cardiometabolic risk factors and thus were not included in the final models.
Separate height adjustments of fat and fat-free mass have been suggested for children by Wells and Cole [15]; accordingly, we found a significant impact of height on androidal fat mass among the participants in our study based on a log-log regression analysis (parameter estimate 0.8286, ). This suggested the use of height-adjusted androidal fat mass (androidal fat mass (g)/height (cm)) (exponent of 0.8286 for height changed to 1 for ease of interpretation). In other words, we adjusted androidal fat mass for height because taller women typically had greater waist circumferences due to their larger frame size and thus appeared to be at higher risk for CVD compared with shorter women, whereas this is not necessarily the case. We performed log transformation of variables with skewed distributions: triacylglycerol, HOMA, CRP, Hcy, and RBC folate.
Restricted maximum likelihood (REML) estimation was used to obtain estimates of variances and correlations between repeated measures. Model selection was guided by a stepwise backwards selection based on model diagnostics, such as Akaike’s information criterion and Bayes information criterion. An overall model fit was obtained based on a likelihood ratio test of the (maximum likelihood fitted) model at hand and the more parsimonious model of only obligatory covariates. Significance indicates a failure to accept the null hypothesis of covariates contributing to the model in a random fashion. More specifically, significance indicates that the full model with all of the variables included explains a greater proportion of the variability in the outcome than the parsimonious model.
3. Results
3.1. Cardiometabolic Risk Factors Assessed during a 36-Month Period
Baseline values for body composition and cardiometabolic outcomes are summarized in Table 1. Body composition measurements, including androidal fat mass, did not change significantly during the course of the study. The isoflavone treatment did not have an effect on any of the body composition outcomes. Treatment compliance was verified using urinary isoflavone concentrations. Compliance was excellent [4]. The results of longitudinal analysis showed that each cardiometabolic outcome model was highly significant, with overall -values for all but the LDL-C model (). Significant covariates for each analyte are shown in Table 2. Androidal fat mass and site (ISU versus UCD) were consistent predictors of all analytes (36 month data) assessed, including lipids/lipoproteins, glucose, and uric acid. In general, androidal fat mass was typically the strongest covariate (positive) in each model (except the model with HDL-C as the outcome, where it was significantly and negatively associated). Time point also emerged as a significant covariate in the HDL-C and glucose models, indicating that the median concentrations of these analytes increased with time; however, the analyte concentrations at subsequent time points were not significantly different compared to baseline. The median glucose concentration was 85 mg/dL at baseline and 88 mg/dL at 36 months. The median HDL-C concentration was 64 mg/dL at baseline and 65 mg/dL at 36 months. Family history of CVD had a significant association with glucose: women who were unaware of their family history (“do not know,” ), taking into account key factors in the model, had on average a 6.8 mg/dL higher glucose concentration compared with women who indicated that they had a positive family history. Isoflavone treatment, TLMP, or HEI did not influence any of the 36 month analytes.
3.2. Cardiometabolic Risk Factors Assessed during a 12-Month Period
Similar to the 36 month outcomes, each 12 month outcome model was highly significant (overall model -values .0001), except for Hcy (, data not shown). Androidal fat mass was positively associated with each outcome (Table 3) and was the strongest covariate (positive) in each model. The parameter estimate for site was negative, indicating that the women from UCD had lower values than the women at ISU. In addition, TLMP () and TLMP-by-site interaction () contributed significantly, and a positive family history of CVD contributed marginally () to fibrinogen concentration. As TLMP increased, the difference in fibrinogen concentration between the sites decreased. Glucose and time were significant covariates in the insulin model. As expected, a higher glucose concentration was associated with a higher insulin concentration. The unadjusted mean insulin concentration was significantly greater (2.86 μU/mL; ) at 12 months compared with baseline. Isoflavone treatment or HEI did not emerge as significant predictors in any of the 12 month models.
4. Discussion
The main objective of this study was to identify significant contributors to cardiometabolic risk in primarily normal and overweight (BMI < 30) healthy postmenopausal women. In agreement with our hypothesis, androidal fat mass was the strongest and most consistent predictor of all cardiometabolic outcomes (except for Hcy) examined in this study. Our participants, as a group, showed no significant changes in body composition measurements, including androidal fat mass, during the course of the study. It should be noted, however, that the women were instructed to maintain their weight throughout the course of the study by following their usual diet and physical activity patterns, and hence, we did not expect to document overall or androidal fat mass change from baseline to 36 months. Similarly, with the exception of insulin, none of the cardiometabolic risk factors changed significantly between baseline and the time of final assessment (either 12 mo or 36 mo). Nevertheless, strong associations between androidal fat and cardiometabolic risk factors suggested that even a small increase in androidal fat mass may have considerable health consequences. Indeed, Biggs et al. [16] reported that the incidence of type 2 diabetes was 70% higher (a hazard ratio of 1.7) in adults 65 years of age and older who gained at least 10 cm in waist circumference over several years compared with those who remained within 2 cm of their baseline waist circumference. Although the Biggs’s study did not specifically target postmenopausal women, it provided some insight into the magnitude of the impact of central adiposity on the incidence of type 2 diabetes.
Certain factors such as duration of menopause, isoflavone treatment, and diet quality that are recognized in the literature as influential with respect to cardiovascular health did not emerge as significant predictors of the cardiometabolic outcomes examined in this study. Some of our findings are consistent with previous reports. For instance, DeNino et al. [17] found that the relationship between age and blood lipids in nonobese women was abolished after controlling for visceral fat. Similarly, in our study, TLMP did not emerge in any model (except for fibrinogen) as a significant covariate, because time point took precedence over TLMP. Results also suggested that androidal fat mass in combination with other factors predominated in these cardiometabolic risk models. On the other hand, our findings for isoflavone treatment and dietary quality conflict with some of the previously published reports. The lowering effect of soy protein rich in isoflavones on blood lipids/lipoproteins, albeit modest (≤6% reduction), has been documented [18, 19]. The effect is more profound in hypercholesterolemic individuals and in males. Similarly, the effect of diet quality on CVD outcomes is well known: a healthier diet (high in fiber and low in fat and sodium) is associated with a lower risk of CVD [20, 21]. We included both a HEI score as a measure of diet quality and isoflavone treatment in our analysis. None of the models retained either covariate. The lack of association between HEI and cardiometabolic outcomes in our study could be in part explained by a majority of women who had relatively healthy diets (median HEI of 67 out of 100). Very few studies have examined the association between HEI and cardiometabolic outcomes with mixed results. For instance, Kant and Graubard [22] reported that HEI emerged as a negative predictor of serum Hcy, CRP and plasma glucose (), whereas Fung and colleagues [23] did not find an association between HEI and CRP. In a recent study of 125 multiethnic overweight and obese women in early postpartum, the HEI scores were negatively associated with LDL-C and total cholesterol and positively related to HDL-C after adjustment for energy intake, body weight, and lactation status [24].
Based on the results of our study, HEI does not appear to be a predictor of cardiometabolic risk factors. On the other hand, the effects of dietary factors as well as isoflavone treatment on cardiometabolic outcomes may be mediated by androidal fat mass. We previously reported that soy isoflavone treatment for 12 months did not exert an effect on body composition, including androidal fat, in our sample of women [25]. With respect to diet, Fox et al. [26] determined that premenopausal women who participated in a 24 week diet and/or exercise program showed reductions in weight (~7 kg) and total percentage body fat (although it remained greater than 35% for all groups), but no significant improvements in blood lipids, glucose, or insulin concentrations. Further analysis revealed a lack of change in the waist-to-hip ratio, which in turn indicated that body fat distribution was not influenced by the intervention. To summarize, dietary interventions that are not sufficiently potent to reduce androidal fat mass are not likely to produce changes in cardiometabolic outcomes.
In addition to androidal fat mass, the other most common independent predictors of cardiometabolic outcomes in our study included time and site. Although concentrations of HDL-C and glucose did not change significantly during the study, both analytes showed an upward trend from baseline to 36 mo. The effect of time could be related to biological and/or behavioral factors that were not included in our models. In fact, TLMP was no longer significant when time was included, indicating that time was more important than TLMP. The site variable also emerged as a significant covariate in all models: compared with women at ISU, women at UCD had lower concentrations of triacylglycerol, HDL-C, glucose, insulin, HOMA, CRP, and fibrinogen and higher concentrations of LDL-C and uric acid. The site effect may be in part explained by baseline characteristics: ISU women were slightly (albeit not significantly) heavier and younger compared with UCD women. It is also quite possible that the site effect was due to some inherent differences between the two geographic locations. We did not document a relationship between androidal fat mass, as well as other factors, and Hcy likely because the women in our study were well below the clinical cut-off of 15, and Hcy values did not indicate great variability.
Some of the strengths of our study were that we followed a relatively large number of women and monitored longitudinal changes (baseline to 12 or 36 mo) in cardiometabolic risk factors. Limitations were that these women were free-living (both a strength and limitation), lacked ethnic diversity, and underwent infrequent measurements (yearly). Women in this study were mainly nonobese based upon the BMI definition. However, BMI does not adequately reflect body composition. Thus, some of our women may have been identified as obese using other criteria, such as percentage body fat (35–40% = overweight, >40% = obese [27]). In conclusion, although cardiometabolic outcomes examined in this study were not affected by isoflavone treatment, each outcome (except for Hcy) had a strong, significant association with androidal fat mass. Thus, even small changes in androidal fat are likely to alter the cardiometabolic profile in healthy postmenopausal women.
Conflict of Interests/Financial Disclosures
Coauthors have not reported conflicts or disclosures.
Funding
The SIRBL clinical trial was supported mainly by a Grant (RO1 AR046922) from the National Institute of Arthritis and Musculoskeletal and Skin Diseases (NIAMS). The project was also supported by the Nutrition and Wellness Research Center, Iowa State University; USDA/ARS, Western Human Nutrition Research Center, Clinical and Translational Science Center, Clinical Research Center at the University of California (1M01RR19975-01), and National Center for Medical Research (UL1 RR024146). Archer Daniels Midland Company (Decatur, IL) donated soy isoflavone tablets (Novasoy) and GlaxoSmithKline (Moon Township, PA) donated calcium and vitamin D supplements (Os-Cal).
Acknowledgments
Study concept and design; secured funding: Alekel, Matvienko, Van Loan. Acquisition of data: Alekel, Bhupathiraju, Matvienko, Perry, Ritland, Van Loan. Study supervision: Alekel, Van Loan. Administrative, technical support: Alekel, Van Loan. Statistical analysis and support: Hofmann. Analysis and interpretation of data: Alekel, Hofmann, Matvienko. Drafting of paper: Matvienko. Critical revision of paper for important intellectual content: Alekel, Bhupathiraju, Hofmann, Matvienko, Reddy. Final approval of paper: Alekel, Bhupathiraju, Hofmann, Matvienko, Perry, Reddy, Ritland, Van Loan. The SIRBL study team would like to thank all of our participants, since without their dedication, our study could not have been completed. The authors would like to acknowledge our phlebotomists and students (graduate and undergraduate alike) who reported early and steadfastly for testing at our clinic sites. The authors thank the James R. Randall Research Center, Archer Daniels Midland Company (Decatur, IL) that supplied free-of-charge the ingredients, using certified good manufacturing procedures, for the treatment tablets (Novasoyâ), as well as Atrium Biotechnologies Inc. that compressed the ingredients into tablets. The authors thank GlaxoSmithKline (Moon Township, PA) for donating the calcium and vitamin D supplements (Os-Cal). The authors would also like to thank our DSMB (Dennis Black, Ph.D., Chair) and Joan McGowan, Ph.D., Director, Musculoskeletal Diseases Branch at NIAMS, who provided scrutiny, guidance, and valuable feedback throughout the trial.