|
|
||||||||





,||
,
,1
* Cancer Prevention Program, Fred Hutchinson Cancer Research Center, Seattle, Washington 98109;
Department of Epidemiology, University of Washington, Seattle, Washington 98195;
Interdisciplinary Program in Nutritional Sciences, University of Washington, Seattle, Washington 98195;
Department of Chemistry, University of Helsinki, Helsinki, Finland 00014; and || Epidemiology Program, Fred Hutchinson Cancer Research Center, Seattle, Washington 98109
1To whom requests for reprints should be addressed at Fred Hutchinson Cancer Research Center, 1100 Fairview Avenue North, M4-B402, P.O. Box 19024, Seattle, WA 981091024. E-mail: jlampe{at}fhcrc.org
| Abstract |
|---|
|
|
|---|
Key Words: daidzein equol O-desmethylangolensin intestinal bacteria soy family
| Introduction |
|---|
|
|
|---|
The variability in bacterial daidzein metabolism may be a determinant of human disease risk, particularly sex hormonerelated conditions. The capacity to produce equol has been associated with circulating reproductive hormones in premenopausal women (12), inversely associated with mammographic density in postmenopausal women (13) and prostate cancer in Japanese men (10, 14), and positively associated with bone mineral density in response to isoflavone intervention (11). Among postmenopausal women, the capacity to produce O-DMA has been positively associated with urinary 2-hydroxyestrone excretion (15) and with mammographic density (13), an intermediate marker of breast cancer risk (16).
Both genetic predisposition and shared environment have been associated with overall intestinal bacteria profile (1723) and with other markers of the presence of particular intestinal bacteria, such as the breath methane-excreter phenotype (24, 25). However, little is known about the specific factors that influence the capacity to harbor particular intestinal bacteria, including daidzein-metabolizing bacteria. We therefore conducted a family study in order to evaluate genetic and nongenetic factors associated with the equol-producer and O-DMAproducer phenotypes.
| Materials and Methods |
|---|
|
|
|---|
Study Procedures.
When signed consent forms were received at FHCRC from at least three family members, each individual was mailed a urine collection kit (Fisher Scientific, Pittsburgh, PA), a questionnaire, and soy food items. In order to ascertain daidzein-metabolizing phenotypes, each participant supplemented his/her usual diet with a soy food item once per day for three consecutive days (soy bar, Revival Soy, ~83 mg daidzein per day or soy nuts, GeniSoy, ~10 mg daidzein per day). Because phenotype determination was based on the presence or absence of equol/O-DMA, the difference in daidzein dose between these two foods did not bias phenotype determination. On the morning of the fourth day, each participant collected a first-void urine sample (5080 ml). Participants mailed urine samples and completed questionnaires to FHCRC. Urine was aliquoted and stored at 20°C until analysis.
Questionnaire Data.
Questions included date of birth, race, education, general health, height, weight, nutritional supplement use, general dietary restrictions, smoking, alcohol consumption, and in relation to the soy challenge, the date of urine collection and amount of soy bar or nuts consumed each of the 3 days. For additional dietary data, participants were asked to indicate the number of servings (per day, week, or month) of particular categories of food eaten in the past month. Two separate, but similar, questionnaires were used for adults and children (individuals 1018 years of age). The questionnaire for children excluded questions about smoking and alcohol consumption, and, specific to females, pregnancy history and hormonal contraceptive use.
Isoflavonoid Analysis.
Urine samples were ether-extracted (26) and analyzed for isoflavonoids (equol, O-DMA, and daidzein) by gas chromatographymass spectrometry (GC-MS) as described elsewhere (15). The mean intraassay coefficients of variation (CVs) for isoflavonoids in the quality control urine samples, measured in duplicate in each batch, were <7%. The mean interassay CVs in the quality control urine samples were <26%. Given the sensitivity of the assay, urine concentrations less than 182 nmol/l (44 ng/ml) of equol and 170 nmol/l (44 ng/ml) of O-DMA were considered below detectable limit. Equol and O-DMA producers were defined as individuals with any detectable concentration of equol and O-DMA, respectively. Urine samples that took more than 14 days to arrive at the laboratory after the reported date of collection were excluded from analysis (n = 6).
Isoflavone Stability.
In preparation for measuring the isoflavonoid content of urine specimens shipped at ambient temperature, we conducted a study to determine the stability of equol and O-DMA in urine kept at room temperature. Four known equol producers (not from the family study sample) followed the same soy consumption and urine collection protocol as described in Study Procedures section, with the exception that urine samples were brought to FHCRC on the day of collection. One aliquot was frozen on Day 0 and additional aliquots were kept at room temperature and then frozen every 2 days for 14 days. Aliquots were subsequently analyzed for isoflavonoids.
Creatinine Analysis.
To assure sufficient concentration of the urine samples, urinary creatinine concentrations were measured on aliquots from all overnight urine samples based on a kinetic modification of the Jaffe reaction using the Roche Reagent for Creatinine (Roche Diagnostic Systems, Nutley, NJ) on a Roche Cobas Mira Plus chemistry analyzer, after first diluting the samples with distilled water 1:50. The assay was linear to 20 mg/dl and the intraassay and interassay CVs were 1%2%.
Statistical Analysis.
In order to test the stability of isoflavones left at room temperature, between-day CVs were calculated for the eight time points for each of the four urine samples.
Associations between nongenetic factors and daidzein-metabolizing phenotypes were analyzed in 410 individuals, and familial aggregation and segregation was evaluated in 366 individuals. Nine participants who reported consumption of less than half of the soy bar on any day were excluded from analysis. Forty-four individuals who completed the daidzein-metabolizing phenotyping were not included in aggregation and segregation analyses because of insufficient number of informative family members completing the study, that is, not having at least two familial relationships that were being considered in this analysis.
To evaluate familial aggregation of the daidzein-metabolizing phenotypes, uniform weight intraclass and interclass correlations (ICC) and standard errors were calculated for familial relationships between individuals using the FCOR program of S.A.G.E. (version 4.4, Case Western Reserve University, Cleveland, OH). EE software (Quantitative Genetic Epidemiology, FHCRC, Seattle, WA) was used to calculate age-adjusted ICC.
To evaluate inheritance of the equol-producer and O-DMAproducer phenotypes, segregation analysis was performed using the SEGREG program of S.A.G.E. under the Class A regressive logistic model (27). Mendelian inheritance was assumed to be through an autosomal locus with two alleles A and B, where the A allele was associated with the relevant phenotype. The likelihood for family data (28) was calculated as a function of the genotype-specific baseline susceptibility parameters (ßAA, ßAB, ßBB); an age parameter (
); the population allele frequency (qA) assuming Hardy-Weinberg equilibrium; and the probability that a parent with each genotype will transmit the A allele (
AA
AB,
BB). Residual associations (
) among relative pairs were allowed. Seven transmission models were fit to the data: a Mendelian dominant model, a Mendelian recessive model, a Mendelian codominant, a "no major type"model that assumes random susceptibility, an environmental model with two susceptibility groups, an environmental model with three susceptibility groups, and a general (full) model in which all parameters in the likelihood function were estimated without restriction. Parameters in each model were estimated using maximum likelihood methods, and the likelihood of each nested model was compared with that under a general model using a likelihood ratio test. Models with likelihoods that were not statistically significantly worse than the likelihood of the general model may be considered to fit the data better than models that were statistically significantly different than the general model. Akaikes information criterion (AIC), calculated as 2(ln L) + 2nparameters, identifies the model with the most parsimo-nious fit to the data.
Fractional polynomial regression was also used to model the association between age and each daidzein-metabolizing phenotype (29); these models and associated plots provide a good visual representation of potentially nonlinear trends. Eight power terms (2, 1, 0.5, 0, 0.5, 1, 2, 3) for age were evaluated, and the best fit model was selected based on deviance value. A first-degree model was chosen because no significant gains were achieved with additional terms in the models. Associations between daidzein-metabolizing phenotypes and other nongenetic factors were analyzed using generalized estimating equations (GEE). Because it was not possible to assign separate correlations for each of the types of relative pairings, an exchangeable correlation across familial relationships was used to control for familial correlation. Because the hypothesis to be tested was that particular nongenetic factors give rise to daidzein-metabolizing phenotypes, equol-producer and O-DMAproducer phenotypes were considered as dependent variables and nongenetic factors were considered as independent variables. Effect modification by sex of several factors potentially related to sex hormone concentrations, such as smoking, anthropometry, age, and diet, was modeled for each phenotype by including an interaction term between sex and each factor using GEE as described for main effects analyses. SAS (version 8.0, SAS Institute Inc., Cary, NC) and Stata (version 8.0, Stata Corp, College Station, TX) were used.
All confidence intervals (CI) are expressed as 95% CI. Because each phenotype is a marker of the capacity to harbor particular bacteria that is distinct from the other phenotype, the factors associated with each phenotype are not necessarily expected to be the same; thus, the results and discussion are presented in relation to each phenotype rather than in relation to the exposures.
| Results |
|---|
|
|
|---|
|
Daidzein-Metabolizing Phenotypes.
All participants had urinary daidzein concentrations above the assay detection limit (all were
109 ng/ml) and had acceptable urinary creatinine concentrations (all were >8 mg/dl).
One hundred seventy-two (42%) individuals were equol producers, and 337 (82%) individuals were O-DMA producers. In our study, 202 (85%) equol nonproducers were O-DMA producers and 135 (78%) equol producers were O-DMA producers. The correlation between the equol-producer and O-DMAproducer phenotypes within an individual was weak (ICC = 0.09 [standard error = 0.11]), with no difference in relation to sex (ICCmale = 0.12 [SE = 0.09], ICCfemale = 0.07 [SE = 0.07]).
Equol-Producer Phenotype.
Weak correlations for main familial relationships were observed in relation to the equol-producer phenotype (Table 1
). Correlations were stronger for father or mother with son than father or mother with daughter and positive correlations between sisters and between sister and brother and an inverse correlation between brothers were observed. Adjustment for age did not alter the overall interpretation of these results.
|
|
|
|
|
O-DMAProducer Phenotype.
An inverse correlation in O-DMAproducer phenotype between parents was observed, and the correlation between siblings was stronger than the correlation between parents and children (Table 1
). There was no correlation between father and son, but positive correlations for other parent:offspring subtypes were observed. Positive correlations between sisters and between sister and brother and an inverse correlation between brothers were observed. Adjustment for age did not alter the overall interpretation of these results.
In segregation analyses, the fit of the Mendelian dominant and recessive models, no major gene, and the environmental models were all statistically significantly worse than the general model (Table 5
). A better fit was achieved when spousal residual associations were also estimated for these models; thus, a likelihood ratio comparison for the codominant model could not be made because the number of estimated parameters was the same as for the general model. Within the Mendelian models, the dominant model was not significantly different than the codominant model (P = 0.86), nor was the recessive model significantly different than the codominant model (p = 0.74).
|
Smoking, having companion animals or livestock, and most dietary factors were not significantly associated with the O-DMAproducer phenotype (Table 4
). However, lower consumption of caffeinated beverages and greater consumption of fried foods were associated with being O-DMA producer.
There were no interactions between sex and age, age-adjusted height or BMI, or, in adults, years of education or alcohol consumption in relation to O-DMAproducer phenotype (P interaction > 0.2). Associations were stronger among females than among males for a 10-cm increase in age-adjusted height (ORfemale = 0.12, CI = 0.03, 0.51; ORmale = 0.67, CI = 0.25, 1.81; P interaction = 0.08), a 10-kg increase in age-adjusted weight (ORfemale = 0.77, CI = 0.68, 0.87; ORmale = 0.94, CI = 0.84, 1.04; P interaction = 0.03), and weekly servings of caffeinated beverages (ORfemale = 0.95, CI = 0.92, 0.99; ORmale = 0.99, CI = 0.96, 1.02; P interaction = 0.11). Compared with never-smokers, former-smoking females were less likely to be O-DMA producers and former-smoking men were more likely to be O-DMA producers (ORfemale = 0.51, CI = 0.25, 1.06; ORmale = 2.10, CI = 0.64, 6.82; P interaction = 0.07), and current smoking was not associated with O-DMAproducer phenotype in relation to sex (ORfemale = 0.97, CI = 0.96, 0.99; ORmale = 0.99, CI = 0.98, 1.00; P interaction = 0.89).
| Discussion |
|---|
|
|
|---|
We observed prevalences of equol producers and O-DMA producers similar to those observed in other studies (510). Our results suggest that these phenotypes are not strongly correlated within an individual. Similarly, we observed in another study, a weak, nonsignificant within-individual association; among postmenopausal women, 95% of equol producers were O-DMA producers, and 79% of equol nonproducers were O-DMA producers (15).
The results of our segregation analyses suggest that the inheritance of equol-producer phenotype may be consistent with autosomal dominant transmission, despite the fact that the familial correlations we observed were weak. These observations are not necessarily inconsistent; observing a high degree of familial correlation does not prove the existence of a genetic mechanism, and, likewise, a low degree of familial correlation does not exclude a genetic mechanism (30). Segregation analysis may be considered strong statistical evidence regarding, although not proof of, genetic control (30).
There was some suggestion that older adults were less likely to be equol producers than children or young to middle-aged adults. These results may be explained by either a cohort effect or age-related changes in intestinal bacteria populations. A cohort effect has been hypothesized for Helicobacter pylori (31), which colonizes in the gastrointestinal tract; lower prevalence has been observed in younger individuals and may be associated with improved hygiene over time. We would expect to observe a higher prevalence of equol producers in older individuals if changes in hygiene were associated with eradication of the equol-producing bacteria; thus, the lower prevalence we observed in older individuals is less likely to be a cohort effect. In two cross-sectional studies, bacterial diversity in healthy young adults was greater than the diversity in healthy elderly adults (32, 33), suggesting that there are age-related changes in intestinal microflora, but this has not yet been evaluated prospectively.
We observed a positive association of greater education with being an equol producer. There is no apparent biologic mechanism by which education may be associated with equol-producer phenotype. However, education may be serving as a marker for other exposures, such as geographic location, medication use, or exposure to environmental contaminants and infectious agents.
Because in other classes of mammals, all appear to be equol producers (11), we evaluated an association with having companion animals. Equol-producer phenotype was not associated with having companion animals, suggesting that shared environment with or frequent exposure to other mammals may not influence equol-producer phenotype.
Several observational studies have reported differences in dietary intake, particularly fat and fiber, in relation to the equol-producer phenotype (5, 7). However, in intervention studies of fiber and soy, no alterations in equol-producer phenotype have been observed in response to interventions (8, 3436). We observed small or no differences in dietary factors and equol-producer phenotype. However, our evaluation of diet was not very sensitive and would have likely detected only large differences between producers and nonproducers.
The familial correlations observed for the O-DMAproducer phenotype are consistent with some degree of genetic predisposition. However, segregation analyses did not provide evidence for genetic or for environmental influence on the O-DMAproducer phenotype. We observed sex-related patterns in the correlations in relationship subtypes, suggesting a possible role for sex-linked heritability or differences in shared environmental factors or differences in response to environmental factors between male and female parents with their children and between male and female siblings.
Age was inversely associated with the O-DMAproducer phenotype. Specifically, children were more likely to be O-DMA producers than were adults. A similar observation has been seen for intestinal colonization by Oxalobacter formigenes; O. formigenes colonization has been detected in nearly all children (37), but in only about 75%80% of adults (38). These results suggest that factors associated with late childhood may be related to change in the intestinal microflora composition. There are several substantial changes that can occur during this time of life. Perhaps the most substantial change that occurs during this time of life is puberty. Although no studies have evaluated an association between the onset of puberty and changes in intestinal bacteria, studies of oral microflora in females in relation to the onset of puberty and during pregnancy suggest that changes in sex hormones during these times may result in alterations of the oral microbial environment (3941).
We observed that Asian individuals were less likely to be O-DMA producers. Racial differences may reflect differences in genetic predisposition or in environment that influence the likelihood of being an O-DMA producer. However, we had few non-Caucasian participants in our study, and the difference in producer prevalence between Caucasians and Asians may have been observed by chance.
Height and BMI were significantly and independently inversely associated with O-DMAproducer phenotype, and these associations were stronger in females. There are no studies known to have evaluated intestinal bacteria in relation to anthropometric characteristics. In a recent hypothesis paper, Beard and Blaser (42) propose that indigenous microflora are among the many factors that influence adult height, but there is currently little evidence to support this hypothesis. Among women, because of their relationship with hormone concentrations, height and BMI may be serving as surrogates for hormone exposure (43, 44).
Smoking has been associated with oral microflora composition (45), but little is known about smoking and intestinal microflora. We observed no association between O-DMAproducer phenotype and smoking. However, when stratified by sex, women who were former smokers were half as likely to be O-DMA producers as women who were never-smokers, and men who were former smokers were twice as likely to be O-DMA producers as men who were never-smokers. In contrast, in another study, a significant positive association was observed between the O-DMAproducer phenotype and former-smoking status in post-menopausal women (13); former smokers were approximately six times as likely to be O-DMA producers than O-DMA nonproducers. Although the potential mechanism for these associations is unclear, these observations suggest that starting and stopping smoking may be associated with intestinal bacteria in a manner that is different from nonsmoking or continual smoking, and that these associations may be different in relation to sex and, among women, to menopausal status. However, no studies have yet evaluated intestinal bacteria changes in relation to smoking cessation. The existence and magnitude of an association between smoking and the capacity to produce O-DMA remains unclear.
There are a few considerations in the interpretation of our study. First, no correction for ascertainment bias was necessary because participants were not ascertained by phenotype. Second, a relatively simple assessment of dietary factors was used to minimize participant burden, and this may have contributed to the lack of association we observed between phenotypes and dietary factors. Third, in order to better maintain individual privacy, we did not ascertain cohabitation of family members; thus, we were unable to evaluate the degree to which cohabitation or length of cohabitation influences these phenotypes. Fourth, we did not to test for statistical significance of the ICC for equol-producer and O-DMAproducer phenotypes. The pairs included in the calculation of these correlations are not independent; for example, if there are three siblings, each sibling would be included twice to represent all the possible pairings of these siblings. As such, standard methods of significance testing cannot be applied. Because this is the first large-scale family study of these phenotypes and our analyses are largely hypothesis generating, we did not believe it was appropriate to apply more complex methods of significance testing. Nonetheless, the standard errors presented provide information about the precision of our estimates.
Equol-producer and O-DMAproducer phenotypes are markers of the presence of particular intestinal bacteria capable of metabolizing daidzein to equol and O-DMA, respectively. Although results from association studies have suggested potential candidate bacteria responsible for the equol-producer phenotype (46, 47) and the O-DMAproducer phenotype (48), the specific bacteria that are responsible for this conversion in humans have not yet been definitively identified. It is likely that there are interindividual differences in the bacteria responsible for these phenotypes (1).
An understanding of whether and how daidzein-metabolizing phenotypes can be modified may provide an avenue for prevention research in targeted populations. Identifying genetic and nongenetic factors associated with and the bacteria responsible for these daidzein-metabolizing phenotypes is a step toward such an understanding. The results of our study provide information regarding the genetic and nongenetic influences on these markers of intestinal bacteria profile. If the genetic association is confirmed by other studies, these results would provide a foundation for identifying genes associated with capacity to harbor bacteria capable of producing equol. Further evaluation is also warranted to evaluate the associations, and underlying mechanisms, we observed for several nongenetic factors with both daidzein-metabolizing phenotypes. Overall, these results, from the largest study of these phenotypes to date, provide preliminary directions for further study of genetic and nongenetic influence on daidzein-metabolizing phenotypes.
| Acknowledgments |
|---|
| Footnotes |
|---|
Received for publication May 11, 2004. Accepted for publication June 11, 2004.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
K. Jaceldo-Siegl, G. E Fraser, J. Chan, A. Franke, and J. Sabate Validation of soy protein estimates from a food-frequency questionnaire with repeated 24-h recalls and isoflavonoid excretion in overnight urine in a Western population with a wide range of soy intakes Am. J. Clinical Nutrition, May 1, 2008; 87(5): 1422 - 1427. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. Atkinson, K. M Newton, E. J A. Bowles, M. Yong, and J. W Lampe Demographic, anthropometric, and lifestyle factors and dietary intakes in relation to daidzein-metabolizing phenotypes among premenopausal women in the United States Am. J. Clinical Nutrition, March 1, 2008; 87(3): 679 - 687. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Bolca, S. Possemiers, A. Herregat, I. Huybrechts, A. Heyerick, S. De Vriese, M. Verbruggen, H. Depypere, D. De Keukeleire, M. Bracke, et al. Microbial and Dietary Factors Are Associated with the Equol Producer Phenotype in Healthy Postmenopausal Women J. Nutr., October 1, 2007; 137(10): 2242 - 2246. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. D. R. Setchell and S. J. Cole Method of Defining Equol-Producer Status and Its Frequency among Vegetarians J. Nutr., August 1, 2006; 136(8): 2188 - 2193. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. Gu, S. E. House, R. L. Prior, N. Fang, M. J. J. Ronis, T. B. Clarkson, M. E. Wilson, and T. M. Badger Metabolic Phenotype of Isoflavones Differ among Female Rats, Pigs, Monkeys, and Women J. Nutr., May 1, 2006; 136(5): 1215 - 1221. [Abstract] [Full Text] [PDF] |
||||
![]() |
K. B. Song, C. Atkinson, C. L. Frankenfeld, T. Jokela, K. Wahala, W. K. Thomas, and J. W. Lampe Prevalence of Daidzein-Metabolizing Phenotypes Differs between Caucasian and Korean American Women and Girls J. Nutr., May 1, 2006; 136(5): 1347 - 1351. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. M. Harris, E. Besselink, S. M. Henning, V. L. W. Go, and D. Heber Phytoestrogens Induce Differential Estrogen Receptor Alpha- or Beta-Mediated Responses in Transfected Breast Cancer Cells Experimental Biology and Medicine, September 1, 2005; 230(8): 558 - 568. [Abstract] [Full Text] [PDF] |
||||
![]() |
Y. Ma, D. Chiriboga, B. C. Olendzki, R. Nicolosi, P. A. Merriam, and I. S. Ockene Effect of Soy Protein Containing Isoflavones on Blood Lipids in Moderately Hypercholesterolemic Adults: A Randomized Controlled Trial J. Am. Coll. Nutr., August 1, 2005; 24(4): 275 - 285. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. Atkinson, C. L. Frankenfeld, and J. W. Lampe Gut Bacterial Metabolism of the Soy Isoflavone Daidzein: Exploring the Relevance to Human Health Experimental Biology and Medicine, March 1, 2005; 230(3): 155 - 170. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |