Lara S. Yoon1
Jonathan P. Jacobs1,2,3
Jessica Hoehner4
Ana Pereira5Juan Cristóbal Gana6Camila Corvalán5Karin B. Michels1,7*- 1Department of Epidemiology, Fielding School of Public Health, University of California, Los Angeles, CA, United States
- 2Vatche and Tamar Manoukian Division of Digestive Diseases, Department of Medicine, David Geffen School of Medicine, University of California, Los Angeles, CA, United States
- 3Division of Gastroenterology, Hepatology and Parenteral Nutrition, Veterans Administration Greater Los Angeles Healthcare System, Los Angeles, CA, United States
- 4Leidos, Inc., Atlanta, GA, United States
- 5Institute of Nutrition and Food Technology, University of Chile, Santiago, Chile
- 6Department of Pediatric Gastroenterology and Nutrition, Division of Pediatrics, School of Medicine, Pontificia Universidad Católica de Chile, Santiago, Chile
- 7Institute for Prevention and Cancer Epidemiology, Faculty of Medicine and Medical Center, University of Freiburg, Freiburg, Germany
The gut microbiome has been linked to breast cancer via immune, inflammatory, and hormonal mechanisms. We examined the relation between adolescent breast density and gut microbial composition and function in a cohort of Chilean girls. This cross-sectional study included 218 female participants in the Growth and Obesity Cohort Study who were 2 years post-menarche. We measured absolute breast fibroglandular volume (aFGV) and derived percent FGV (%FGV) using dual energy X-ray absorptiometry. All participants provided a fecal sample. The gut microbiome was characterized using 16S ribosomal RNA sequencing of the V3-V4 hypervariable region. We examined alpha diversity and beta diversity across terciles of %FGV and aFGV. We used MaAsLin2 for multivariable general linear modeling to assess differential taxa and predicted metabolic pathway abundance (MetaCyc) between %FGV and aFGV terciles. All models were adjusted for potential confounding variables and corrected for multiple comparisons. The mean %FGV and aFGV was 49.5% and 217.0 cm3, respectively, among study participants. Similar median alpha diversity levels were found across %FGV and aFGV terciles when measured by the Shannon diversity index (%FGV T1: 4.0, T2: 3.9, T3: 4.1; aFGV T1: 4.0, T2: 4.0, T3: 4.1). %FGV was associated with differences in beta diversity (R2 =0.012, p=0.02). No genera were differentially abundant when comparing %FGV nor aFGV terciles after adjusting for potential confounders (q > 0.56 for all genera). We found no associations between predicted MetaCyc pathway abundance and %FGV and aFGV. Overall, breast density measured at 2 years post-menarche was not associated with composition and predicted function of the gut microbiome among adolescent Chilean girls.
Introduction
Breast cancer is the most common cancer among women worldwide, however, one third of affected women have no known or suspected risk factors. Greater breast density is strongly associated with increased risk of breast cancer among adult women (Boyd et al., 1998). Breast density, measured as the relative proportion of fibroglandular tissue to fatty tissue in the breast, is inversely related to age, parity, and later menopause and also associated with childhood and adolescent body fatness (Martin and Boyd, 2008; Harris et al., 2011). Initial peak breast density is hypothesized to be established during adolescence (Boyd et al., 2009; Ghadge et al., 2020), and breast tissue may be particularly vulnerable to exposures during puberty, a period of rapid breast development (Ruder et al., 2008). Adolescent breast density has been associated with pubertal maturation and body fatness (Novotny et al., 2011). However, few studies have characterized other important associates of breast development and composition during this vulnerable period.
A novel mechanistic pathway that may contribute to development of the breast is the gut microbiome. Collectively the trillions of microbes living in the human intestinal tract, the gut microbiome plays important roles in numerous biological processes including immune regulation, dietary metabolism, epithelial barrier function, and hormone regulation (Ley et al., 2006; Turnbaugh et al., 2006; Belkaid and Hand, 2014; Martin et al., 2019). Experimental and epidemiologic evidence suggests that the gut microbiome may be associated with different diseases such as metabolic disorders and cancers, including breast cancer, through direct and indirect mechanisms (Selber-Hnatiw et al., 2017; Fernández et al., 2018; Chen et al., 2019; Eslami-S et al., 2020; Laborda-Illanes et al., 2020; Sepich-Poore et al., 2021). For instance, perturbations in microbiota composition, or dysbiosis, may lead to systemic inflammation resulting in an increased vulnerability to pathogens (Zitvogel et al., 2015). Carcinogenesis may follow from dysbiosis-induced permeability of the intestinal epithelium and resulting extracellular vesicle circulation throughout the body in biofluids (Ricci et al., 2020). Changes in the composition of gut microbiome, and specifically in microbes involved in estrogen metabolism, may influence breast cancer development through an increase in circulating estrogen levels (Goedert et al., 2015). Evidence for an association between circulating estrogens and adult mammographic breast density is mixed (Boyd et al., 2002; Fuhrman et al., 2012; Gierach et al., 2015; Jones et al., 2019). However, we have shown that levels of prepubertal estrogen, measured with an ultra-sensitive method, are associated with earlier thelarche, which is in turn related to breast density at the end of puberty (Pereira et al., 2015). Estrogen levels are directly associated with breast development; therefore, the gut microbiome may play a mechanistic role in the development of breast density during puberty and adolescence (Pereira et al., 2015; Baker et al., 2017).
It is unknown whether the gut microbiome is associated with breast density in adolescence. The current study examined the association between breast density at two years post-menarche and the gut microbiome in a cohort of adolescent Chilean girls. We hypothesized that the microbial composition and function of the gut would differ across densities of the breast.
Materials and Methods
Study Design and Population
We conducted a cross-sectional analysis of a subset of female participants in the Growth and Obesity Cohort Study (GOCS). The original GOCS, which began in 2006, has been described previously (Corvalán et al., 2009). In brief, 1,196 children aged 2.5 to 4 years from low- and middle- income families in the South East region of the Santiago metropolitan area in Chile and enrolled in preschool at the National Board of Preschool Council Program (Junta Nacional de Jardines Infantiles) were enrolled in the study; approximately half (601) were girls. Participants in the study visit the Institute of Nutrition and Food Technology (INTA) Health Clinic at the Universidad de Chile in Santiago, Chile at least once per year for anthropometric assessments, pubertal (Tanner) evaluation, collection of biospecimens, and to complete 24-hour dietary recall interviews. A limited set of behavioral and demographic information was also collected via questionnaire. Breast composition was measured when the participants were two years post-menarche (2PM). The current study included a subset of 218 girls randomly selected to provide a fecal sample and who had a breast composition measurement at 2PM. The 2PM breast assessment and fecal sample collection occurred between 2018 and 2019. The study protocol and written consent forms were approved by the University of Chile Ethics Committee at INTA.
Breast Composition Measurement
Breast composition was measured at the clinic visit corresponding to a timepoint of 2PM for each girl using the dual energy X-ray absorptiometry (DXA) breast scanning protocol developed by Shepherd et al. at the University of California, San Francisco (Shepherd et al., 2006). The Prodigy DXA system software (version 13.6, series 200674; GE Healthcare) was used to scan each breast for quantification of adipose fat and fibroglandular (FG) tissue. Stable calibration of the system was continuously performed using a quality control breast phantom. Absolute fibroglandular volume (FGV; cm3) and total breast volume for each breast were derived from a two-compartment model of adipose fat and fibroglandular tissue. The percentage of FGV tissue in the breast was derived by dividing the absolute FGV by total breast volume and multiplying by 100. Percent FGV and absolute FGV for the left and right breast were averaged to obtain two single measures of breast density: percent FGV (%FGV) and absolute FGV (aFGV). DXA is frequently used in studies of bone density in children; exposure to ionizing radiation from the DXA protocol is low (Crabtree et al., 2014). The DXA approach for breast composition assessment has high validity and precision among adolescent girls (Shepherd et al., 2008). Both breast density outcomes were categorized into terciles (‘T1’, ‘T2’, ‘T3’) based on the distribution of the sample for statistical analyses.
Fecal Collection
Fecal samples were collected after annual visits to the INTA health clinic occurring when the girls were between 13 and 15 years old. During the visit, girls were provided with materials (sealed plastic bag, stool catcher, a plastic, sterile container with a spoon to manipulate the sample and procedure gloves) and instructions for at-home fecal collection. Briefly, after the stool deposit, girls were asked to collect a part of the stool approximately the size of a walnut and place it in the plastic sterile container and label with the date and time of sample collection. Within 15 minutes of collection, the samples were sealed and then stored temporarily in the participant’s freezer, at which time study personnel were contacted to schedule the sample pick up. Samples were retrieved from the participant’s home no later than three days after initial collection, transported on ice to INTA, labeled with a de-identified key, and stored in a -80°C freezer prior to shipment. All stool samples were shipped on dry ice to the National Exposure Assessment Laboratory at Emory University.
Fecal Processing and 16s rRNA Sequencing
Fecal processing of 279 samples (266 unique samples and 13 duplicates) was performed at the National Exposure Assessment Laboratory at Emory University, a Children’s Health Exposure Analysis Resource Laboratory. DNA was extracted from fecal samples using the Qiagen DNeasy PowerSoil Kit (Qiagen; 12888) with the TissueLyser at maximum speed (30 beats per second) for 10 minutes. Composition of the gut microbiome was determined through sequencing and amplification of the V3-V4 hypervariable region of the 16S rRNA gene (PCR primer pair: 341F 5’-GTGCCAGCMGCCGCGGTAA-3’; 805R 5’-GACTACHVGGGTWTCTAAT-3’). Libraries were made from 12.5 ng of DNA following a standard 16S Metagenomic Library Preparation Workflow from Illumina, Inc. Libraries were pooled in equal amounts based on fluorescence quantification, resulting in a 630 bp amplicon. Final library pools were quantitated via qPCR (Kapa Biosystems; KK4824). The pooled library was sequenced on an Illumina MiSeq using MiSeq v3 600 cycle chemistry (Illumina MS-102-3003) at a loading density of 6-8 pM with 20% PhiX, generating roughly 20M, 300 bp paired-end reads. In addition to the 279 experimental samples, traditional negative, no template control (NTC) negative, positive, and ZymoBIOMICS mock microbial community controls were included in the assays (Pollock et al., 2018). The Emory Integrated Genomics Core performed the assays.
After sequencing, demultiplexed raw amplicon sequences were processed using QIIME2 (Quantitative Insights Into Microbial Ecology) version 2019.4 (Bolyen et al., 2019). Denoising and dereplication, including chimera removal and trimming of reads based on Phred quality scores, were performed using the Divisive Amplicon Denoising Algorithm 2 (DADA2) module (Callahan et al., 2016a). The majority of Phred scores were greater than 25. To improve sensitivity of the algorithm, we trimmed the first 30 base pairs and truncated the reads at position 290, where the mean read quality distribution dropped for the overall dataset. Amplicon sequence variants (ASVs) were inferred using DADA2 to increase resolution and allow for intrinsic biologic meaning (Callahan et al., 2017). Taxonomy was assigned using a naive Bayes classifier on the SILVA database (SILVA 132 release) (Quast et al., 2013). We excluded samples from girls that did not have breast assessments at 2PM (n=48). The end product yielded a total of 6,270 unique ASVs from 218 samples. Standard preprocessing (filtering, subsetting, agglomeration) with the R package ‘phyloseq’ was used to exclude undefined (NA) or ambiguous (uncharacterized) phyla, phyla below a prevalence of 1, and ASVs below a prevalence threshold of 2% of total samples (McMurdie and Holmes, 2013; Callahan et al., 2016b). The final ASV feature table comprised 18,628,903 total reads (mean per sample = 87,872, range 28,118 to 278,065) and a total of 1,600 unique ASVs across 218 samples.
Covariates
Demographic, anthropometric, and nutritional data were collected by trained dietitians during the annual study visit at the health clinic. Age- and sex- adjusted body mass index (BMI; kg/m2) Z-scores were calculated using the World Health Organization growth reference data and categorized into Norma/Underweight (Z-score ≤ 1), Overweight (1 < Z-score ≤ 2), and Obese (2 < Z-score). Body fat percentage was estimated using Tanita-BC-418 MA bioelectrical impedance measurements (Tanita-Corporation, Tokyo, Japan) and categorized into Healthy, Overweight, and Obese based on Tanita children’s age- and sex-specific body fat reference curves (McCarthy et al., 2006; Cediel et al., 2016). Age at menarche was determined via phone interviews completed by study dietitians every three months during puberty and dichotomized (≤ 12 years, > 12 years). Mothers of the participants were present at the clinic visits and were asked to complete short questionnaires with information on their highest level of education (secondary education or less, post-secondary education), birth mode of the participant (cesarean, vaginal), and months of exclusive breast feeding (<3 months, 3-6 months, >6 months). Dietary 24-hour recalls were collected longitudinally beginning in April 2014 by trained dietitians using the USDA multiple-pass method. Food and nutrient information were obtained using a harmonization process that mapped Chilean foods to the USDA Food and Nutrient Database for Dietary Studies. Daily intake of five major food groups (vegetables [g], fruit [g], red or processed meat [g], yogurt [g], and whole grains [g]) and total energy intake (kCal) were averaged over the data collection period and up to the date of the 2PM breast composition assessment to reduce random-measurement error and to obtain a more accurate assessment of long-term diet (Willett, 2012). We also collected data on ethnicity (Mapuche [native Chilean indigenous] and non-Mapuche, according to last name), average daily hours of television during the week (≤1 hour, 1-3 hours, >3 hours) as a proxy of physical activity, and antibiotic use in the 6 months prior to fecal sample (yes, no). Missing covariate data were imputed using last observation carried forward for body fat percentage and dietary variables. We assumed the remaining covariates with missing values, including maternal education, ethnicity, breast feeding, birth mode, antibiotic use, and television hours were missing completely at random. We used mean imputation (breast feeding) or median imputation (maternal education, ethnicity, birth mode, antibiotic use, television hours) in the R package ‘mice’ (van Buuren and Groothuis-Oudshoorn, 2011). The variable with the highest proportion of missing values was television hours (n=22, 10.1%); all other covariates were missing at a proportion less than 3%.
Statistical Analysis
Relative abundance at the phylum level was plotted for all samples and across %FGV and aFGV terciles. We estimated alpha diversity as the observed richness (the number of species per sample) and the Shannon index (a measure of richness and evenness) on unfiltered data (Shannon, 1948; Chao, 1987). Rarefaction without replacement was used to standardize library sizes and account for uneven sampling depth prior to estimating alpha and beta diversity metrics and has been suggested to be suitable for microbiome data (McKnight et al., 2019; Estaki et al., 2020). Overall differences in observed richness and the Shannon index by %FGV and aFGV terciles were tested using the Kruskal-Wallis (KW) followed by post-hoc pairwise testing with the Wilcoxon Rank Sum in cases where data were compatible with evidence to reject the KW null hypothesis. To evaluate whether our results were sensitive to body composition, we performed analysis of covariance (ANCOVA) and adjusted for body fat percentage. Beta diversity was visualized using principal coordinate analysis (PCoA) plots with Bray-Curtis dissimilarity. Single and multivariate permutational analysis of variance (PERMANOVA) were used to test for overall differences in microbial composition between %FGV and aFGV terciles and produced marginal effects using the adonis2 function in the ‘vegan’ package (McArdle and Anderson, 2001; Oksanen et al., 2007). Homogeneity of variance was used to test whether differences in community structure were due to dissimilar dispersions. Supplemental PCoA visualization and PERMANOVA tests were done using weighted UniFrac, which differs from Bray-Curtis in that it incorporates phylogenetic relationships (Lozupone et al., 2007). Dietary intake of whole grains, fruit, vegetables, and yogurt are important predictors of gut microbial composition and may influence breast density. Therefore, we performed supplementary PERMANOVA analyses with added terms for four food-specific intake groups, including average whole grain intake (g/day), average fruit intake (g/day), average vegetable intake (g/day), and average yogurt intake (g/day). Microbial diversity metrics were estimated using the R packages ‘phyloseq’ and ‘vegan’ (Oksanen et al., 2007; McMurdie and Holmes, 2013).
ASVs from the filtered final feature table were agglomerated to the genus level prior to associating microbial community features with breast density outcomes. A total of 223 genera remained after agglomeration of the 1,600 ASVs. We used the R package ‘MaAsLin2’ (Microbiome Multivariable Associations with Linear Models), which relies on a modified generalized linear model for compositional data, to identify differentially abundant microbe genera across %FGV and aFGV terciles (Mallick et al., 2021). All models specified minimum abundance of 0.01% in 10% of samples, total sum scaling normalization, and arcsine square root-transformation (Mallick et al., 2021). The models were adjusted for a set of potential confounders including age, body fat percentage, antibiotic use, maternal education, total calories, hours of TV watching, ethnicity, birth mode, and breast feeding. Where appropriate, the Benjamini and Hochberg (BH) correction method was used to control the false discovery rate (FDR) and produced q-values (Benjamini and Hochberg, 1995). Supplementary MaAsLin2 analyses were also performed with adjustment for the four food-specific dietary intake groups.
Functional metabolic pathways of microbial communities were predicted from 16S rRNA marker sequencing data using PICRUSt2 implemented as a QIIME2 plugin (Douglas et al., 2020, 2). Default parameters were used in the pipeline, including SATe-enabled phylogenetic placement, hidden-state prediction, and a distance cut-off of 2. PCoA ordination with Bray-Curtis dissimilarity was used to visualize predicted abundance of metabolic pathways present in the MetaCyc Metabolic Pathway Database. Differences in MetaCyc pathway abundance between %FGV and aFGV terciles were testing using multivariable PERMANOVA. Multivariable MaAsLin2 was used to evaluate differential abundance of 458 predicted MetaCyc pathways with default parameters and BH correction for multiple comparisons. As in prior analyses, we additionally adjusted for the four food-specific dietary intake groups in supplementary analyses.
Results
Study Sample Characteristics
The association between breast density and gut microbial composition and function was assessed among 218 GOCS participants. The mean age of the participants was 14 years at the time of breast assessment and 12 years at reported menarche (Table 1). The participants were primarily non-indigenous (Non-Mapuche, 83.0%) and 78.4% of the mothers had a secondary school educational attainment or less. The mean %FGV among the sample was 49.5% (SD=14.5); for tercile 1 (T1), mean %FGV was 34.0%; tercile 2 (T2), 48.4%; and tercile 3 (T3), 65.3%. The mean aFGV among the sample was 217.0 cm3; for T1, mean aFGV was 140.2 cm3; T2, 210.8 cm3; and T3, 292.7 cm3. More than half the cohort was considered overweight or obese. Overall, participants had a mean BMI Z-score of 1.0 and a mean body fat percentage of 32.6%. Additional study population characteristics are reported in Table 1.
Table 1 Population characteristics of 218 girls participating in the Growth and Obesity Cohort Study.
Phylum-Level Microbial Composition of the Sample
The relative composition of the gut microbiome for each %FGV and aFGV tercile is displayed in Figure 1. Overall, Firmicutes was the most highly represented bacterial phylum and comprised 66% of the sample abundance. Microbes from the Bacteroidetes and Actinobacteria phyla represented 18.1% and 10.4% of the abundance, respectively. We found minor differences in relative abundance across %FGV terciles for Bacteroidetes (T1: 20.4%, T2: 15.4%, T3: 18.5%), Actinobacteria (T1: 8.9%, T2: 10.9%, T3: 11.5%) and Euryarchaeota (T1: 1.4%, T2: 0.9%, T3: 0.9%) (Supplementary Figure 1). There were no differences in relative abundance at the phylum level across aFGV terciles (Supplementary Figure 2).
Figure 1 Microbial mean relative abundance at the phylum level stratified by (A) %FGV tercile and (B) aFGV tercile in fecal microbiota samples from 218 GOCS participants.
Global Microbial Diversity Analyses
We calculated alpha diversity for each sample using observed richness and the Shannon diversity index. Overall, the mean number of observed species present in the fecal samples was 214 and the mean Shannon index was 4.0 (Figure 2). We observed differences in median observed species richness across %FGV terciles (T1: 226; T2: 200, T3: 223). There were no differences in the observed number of species across aFGV terciles nor in the Shannon index for both %FGV and aFGV terciles. Results from ANCOVA did not suggest a difference in alpha diversity across %FGV nor aFGV terciles when adjusting for body fat percentage (data not shown).
Figure 2 Box plots of alpha diversity metrics for observed richness and Shannon index across terciles of %FGV (A, B) and aFGV (C, D). P-values are presented for overall differences in alpha diversity metrics (Kruskal-Wallis) and post-hoc pairwise differences. Boxes represent the lower, median, and upper quartile of the data and whiskers are 1.5*interquartile range.
Beta diversity analysis for microbial ASV abundance was performed using Bray-Curtis dissimilarity and visualized by principal coordinate analysis (PCoA) (Figure 3). PCoA axes 1 and 2 represented 9.2% and 5.8% of the total variance, respectively. Results from multivariable PERMANOVA analyses suggested that a small but statistically significant portion of the variability in microbial composition was associated with %FGV (R2 = 0.012, p=0.02) (Table 2). No significant contribution to variability was associated with aFGV (R2 = 0.09, p=0.53). We also did not note any visual shifts in global predicted microbial pathway abundances (Figure 4) nor any significant predictors of Bray-Curtis dissimilarity in PERMANOVA analyses (Table 3). In supplementary analyses of beta diversity using weighted UniFrac, rather than Bray-Curtis, we found no significant association with %FGV and aFGV (Supplementary Table 1). Supplementary PERMANOVA analyses with additional adjustment for dietary predictors did not substantially modify the variation in microbial nor pathway diversity (Supplementary Tables 2, 3).
Figure 3 Principal coordinate analysis (PCoA) plot of microbial composition derived from Bray-Curtis dissimilarity among fecal samples provided by 218 girls in GOCS colored by %FGV (A) and aFGV (B) terciles. Ellipses are 95% confidence regions for each tercile.
Table 2 Multivariable PERMANOVA analyses to identify variation (R2) in microbial beta diversity (Bray-Curtis dissimilarity) explained by study characteristics.
Figure 4 Principal coordinate analysis (PCoA) plot of predicted MetaCyc pathway abundance derived from Bray-Curtis dissimilarity among fecal samples provided by 218 girls in GOCS colored by %FGV (A) and aFGV (B) terciles. Ellipses are 95% confidence regions for each tercile.
Table 3 Multivariable PERMANOVA analyses to identify variation (R2) in MetaCyc predicted pathway beta diversity (Bray-Curtis dissimilarity) explained by study characteristics.
Differential Abundance of Microbial Taxa and Predicted Functional Pathways
In multivariable linear modeling with MaAsLin2, no genera were associated with %FGV nor with aFGV after FDR correction (Supplementary Table 4). These null associations remained after further adjustment for additional dietary predictors, including whole grain, fruit, vegetable, and yogurt intake (Supplementary Table 5). Additionally, there were no major differences in abundance of predicted MetaCyc pathways by %FGV nor aFGV when adjusting for potential confounders (Supplementary Table 6) or important dietary predictors (Supplementary Table 7).
Discussion
In this cross-sectional study of Chilean adolescents, we found primarily null associations of microbial composition and predicted function with breast density. Our results suggest minimal differences in alpha diversity: girls in the lowest and highest breast %FGV terciles had slightly higher observed species richness compared to those in tercile 2. However, this pattern was not observed for aFGV and when examining different alpha diversity metrics (Shannon index). We also noted a small but significant contribution of breast %FGV to beta diversity.
Several recent clinical studies have examined the association between the gut microbiome and breast cancer. In a case-control study of 96 post-menopausal women, pre-treatment breast cancer patients had altered microbial composition (beta diversity) and increased abundance of Clostridiaceae, Faecalibacterium, and Ruminococcaceae compared to controls (Goedert et al., 2015). A cross-sectional study of incident pre- and post-menopausal breast cancer patients reported a less diverse microbiome and differential abundance of Firmicutes in women with human epidermal growth factor receptor 2 (HER2) positive breast cancer compared to HER2 negative breast cancer (Wu et al., 2020). Other studies also support associations of specific microbial taxa (e.g., Bacteroidetes, Blautia spp.) with breast cancer staging and clinical characteristics, including body size (Bard et al., 2015; Luu et al., 2017). Much less information is available from epidemiologic studies on the possible association between the gut microbiome and breast density. A study of healthy menopausal women in the United States found no association between mammographic breast density and gut microbial beta diversity and Firmicutes to Bacteroidetes (F/B) ratio, and suggestive but non-significant differences in alpha and beta diversity (Yaghjyan et al., 2021). These results are comparable to those from a study of cancer-free postmenopausal women, which found that alpha diversity and relative abundance did not differ in women with high versus low mammographic density (Jones et al., 2019). To our knowledge, our study is the first that examines the association between the gut microbiome and breast density in adolescents. Our study differs in several ways from the prior studies of the breast density and gut microbiome relation, notably in study population and breast density assessment method. In contrast to the studies from Yyghjyan and Jones, we found small but significant differences in beta diversity across terciles of breast density. However, we also found similar null associations for alpha diversity and relative abundance with respect to breast density alone. Null findings specifically in studies of the gut microbiome-breast density association might reflect the complex interaction between body composition (e.g., body fatness), the gut microbiome, and breast density. Obesity is strongly inversely associated with breast density and with diversity of the gut microbiome (Ley et al., 2006; Turnbaugh et al., 2006; Dorgan et al., 2012). It is plausible that an association is mediated by body composition, such that any effect is cancelled out when controlling for body fat percentage (Boyd et al., 2011). However, we did not note any breast density and microbial composition associations in unadjusted analyses. It is also possible that studies of the microbiome-breast cancer association are reflective of alterations to the composition of the gut following the disease state, rather than the hypothesized effects of the gut microbiome in contributing to breast cancer pathogenesis through early alteration of estrogen metabolism.
Our study has several limitations. The cross-sectional design does not allow for a temporal specification of the gut microbiome – breast density association. However, the human gut microbiome is thought to be relatively stable after infancy and early childhood (Lozupone et al., 2012). We cannot fully preclude that other perturbations of the microbiome, such as major dietary changes, occurred in the time just prior to sample collection. However, we collected information on antibiotic usage and clinical diagnoses of disease in the six months prior to stool collection, and less than 2% of the sample reported any probiotic consumption. We did not have comprehensive information on nutritional supplements; however, use of vitamins or nutritional supplements is uncommon in Chile, particularly among girls. We had limited information on physical activity, which is associated with gut microbial composition and many anthropometric characteristics (e.g., body fatness) associated with breast density. We also lack direct information on functional microbiome data which may show relations with breast density that were not apparent when examining composition alone. However, we were able to approximate functional potential of the community using PICRUSt2.
Our study also has several strengths, including a large sample of geographically, socioeconomically, and behaviorally similar girls. We were also able to comprehensively assess the gut microbial community using fecal 16s rRNA gene sequencing. The GOCS study has collected longitudinal data on lifestyle, socioeconomic, and anthropometric factors, allowing for specific control of potential confounders. MaAsLin2, a novel general linear modeling approach, allows for preserved statistical power with multiple covariates and control for false-discovery rate.
In conclusion, we found no important association between the gut microbiome and breast density at 2 years post-menarche in female adolescents.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ncbi.nlm.nih.gov/ (SRA, BioProject PRJNA744878), https://www.doi.org/10.36043/1977_490.
Ethics Statement
The studies involving human participants were reviewed and approved by Universidad de Chile Instituto de Nutrición y Tecnología de los Alimentos and the University of California, Los Angeles. Written informed consent to participate in this study was provided by the participants’ legal guardian/next of kin.
Author Contributions
LSY designed the study, analyzed the data, interpreted the results, and wrote the manuscript. KBM conceived the study, contributed resources, and interpreted the study results. JPJ and JH analyzed the data and interpreted the study results. AP, JC, and CC collected biospecimen samples and other study data and contributed resources. All authors contributed to the article and approved the submitted version.
Funding
This research was supported by research grant U01ES026130 from the National Institute of Environmental Health and the National Cancer Institute, National Institutes of Health (to KBM). This research was also supported by the National Institute of Environmental Health Sciences of the National Institutes of Health under Award Number U2CES026560 to Emory University and the Emory Integrated Computational Core (EICC), which is subsidized by the Emory University School of Medicine and is one of the Emory Integrated Core Facilities. Additional support was provided by the National Center for Georgia Clinical & Translational Science Alliance of the National Institutes of Health under Award Number UL1TR002378. LSY was supported by T32 training grant 5T32CA009142 from the National Cancer Institute, National Institutes of Health and the Karen Toffler Charitable Trust. JPJ was supported by VA CDA2 IK2CX001717. The content is solely the responsibility of the authors and does not necessarily reflect the official views of the National Institutes of Health. The article processing charge was funded by the Baden-Wuerttemberg Ministry of Science, Research and Art and the University of Freiburg in the funding programme Open Access Publishing.
Author Disclaimer
The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.
Conflict of Interest
Author JH is employed by Leidos, Inc.
The remaining 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
We would like to thank the participants of the Growth and Obesity Cohort Study and all staff who helped collect, store, and process the samples.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcimb.2021.794610/full#supplementary-material
References
Baker, J. M., Al-Nakkash, L., Herbst-Kralovetz, M. M. (2017). Estrogen–gut Microbiome Axis: Physiological and Clinical Implications. Maturitas 103, 45–53. doi: 10.1016/j.maturitas.2017.06.025
Bard, J.-M., Luu, H. T., Dravet, F., Michel, C., Moyon, T., Pagniez, A., et al. (2015). Relationship Between Intestinal Microbiota and Clinical Characteristics of Patients With Early Stage Breast Cancer. FASEB J. 29, 914.2. doi: 10.1096/fasebj.29.1_supplement.914.2
Belkaid, Y., Hand, T. W. (2014). Role of the Microbiota in Immunity and Inflammation. Cell 157, 121–141. doi: 10.1016/j.cell.2014.03.011
Benjamini, Y., Hochberg, Y. (1995). Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. J. R. Stat. Society Ser. B (Methodol) 57, 289–300. doi: 10.1111/j.2517-6161.1995.tb02031.x
Bolyen, E., Rideout, J. R., Dillon, M. R., Bokulich, N. A., Abnet, C. C., Al-Ghalith, G. A., et al. (2019). Reproducible, Interactive, Scalable and Extensible Microbiome Data Science Using QIIME 2. Nat. Biotechnol. 37, 852–857. doi: 10.1038/s41587-019-0209-9
Boyd, N. F., Lockwood, G. A., Byng, J. W., Tritchler, D. L., Yaffe, M. J. (1998). Mammographic Densities and Breast Cancer Risk. Cancer Epidemiol. Biomarkers Prev. 7, 1133–1144. doi: 10.3233/BD-1998-103-412
Boyd, N., Martin, L., Chavez, S., Gunasekara, A., Salleh, A., Melnichouk, O., et al. (2009). Breast-Tissue Composition and Other Risk Factors for Breast Cancer in Young Women: A Cross-Sectional Study. Lancet Oncol. 10, 569–580. doi: 10.1016/S1470-2045(09)70078-6
Boyd, N. F., Martin, L. J., Yaffe, M. J., Minkin, S. (2011). Mammographic Density and Breast Cancer Risk: Current Understanding and Future Prospects. Breast Cancer Res. 13, 223. doi: 10.1186/bcr2942
Boyd, N. F., Stone, J., Martin, L. J., Jong, R., Fishell, E., Yaffe, M., et al. (2002). The Association of Breast Mitogens With Mammographic Densities. Br. J. Cancer 87, 876–882. doi: 10.1038/sj.bjc.6600537
Callahan, B. J., McMurdie, P. J., Holmes, S. P. (2017). Exact Sequence Variants Should Replace Operational Taxonomic Units in Marker-Gene Data Analysis. ISME J. 11, 2639–2643. doi: 10.1038/ismej.2017.119
Callahan, B. J., McMurdie, P. J., Rosen, M. J., Han, A. W., Johnson, A. J. A., Holmes, S. P. (2016a). DADA2: High-Resolution Sample Inference From Illumina Amplicon Data. Nat. Methods 13, 581–583. doi: 10.1038/nmeth.3869
Callahan, B. J., Sankaran, K., Fukuyama, J. A., McMurdie, P. J., Holmes, S. P. (2016b). Bioconductor Workflow for Microbiome Data Analysis: From Raw Reads to Community Analyses. F1000Res 5, 1492. doi: 10.12688/f1000research.8986.2
Cediel, G., Corvalán, C., Aguirre, C., de Romaña, D. L., Uauy, R. (2016). Serum 25-Hydroxyvitamin D Associated With Indicators of Body Fat and Insulin Resistance in Prepubertal Chilean Children. Int. J. Obes. 40, 147–152. doi: 10.1038/ijo.2015.148
Chao, A. (1987). Estimating the Population Size for Capture-Recapture Data With Unequal Catchability. Biometrics 43, 783–791. doi: 10.2307/2531532
Chen, J., Douglass, J., Prasath, V., Neace, M., Atrchian, S., Manjili, M. H., et al. (2019). The Microbiome and Breast Cancer: A Review. Breast Cancer Res. Treat 178, 493–496. doi: 10.1007/s10549-019-05407-5
Corvalán, C., Uauy, R., Stein, A. D., Kain, J., Martorell, R. (2009). Effect of Growth on Cardiometabolic Status at 4 Y of Age. Am. J. Clin. Nutr. 90, 547–555. doi: 10.3945/ajcn.2008.27318
Crabtree, N. J., Arabi, A., Bachrach, L. K., Fewtrell, M., El-Hajj Fuleihan, G., Kecskemethy, H. H., et al. (2014). Dual-Energy X-Ray Absorptiometry Interpretation and Reporting in Children and Adolescents: The Revised 2013 ISCD Pediatric Official Positions. J. Clin. Densitometry 17, 225–242. doi: 10.1016/j.jocd.2014.01.003
Dorgan, J. F., Klifa, C., Shepherd, J. A., Egleston, B. L., Kwiterovich, P. O., Himes, J. H., et al. (2012). Height, Adiposity and Body Fat Distribution and Breast Density in Young Women. Breast Cancer Res. 14, R107. doi: 10.1186/bcr3228
Douglas, G. M., Maffei, V. J., Zaneveld, J., Yurgel, S. N., Brown, J. R., Taylor, C. M., et al. (2020). PICRUSt2: An Improved and Customizable Approach for Metagenome Inference. bioRxiv 672295. doi: 10.1101/672295
Eslami-S, Z., Majidzadeh-A, K., Halvaei, S., Babapirali, F., Esmaeili, R. (2020). Microbiome and Breast Cancer: New Role for an Ancient Population. Front. Oncol. 10, 120. doi: 10.3389/fonc.2020.00120
Estaki, M., Jiang, L., Bokulich, N. A., McDonald, D., González, A., Kosciolek, T., et al. (2020). QIIME 2 Enables Comprehensive End-To-End Analysis of Diverse Microbiome Data and Comparative Studies With Publicly Available Data. Curr. Protoc. Bioinf. 70, e100. doi: 10.1002/cpbi.100
Fernández, M. F., Reina-Pérez, I., Astorga, J. M., Rodríguez-Carrillo, A., Plaza-Díaz, J., Fontana, L. (2018). Breast Cancer and Its Relationship With the Microbiota. Int. J. Environ. Res. Public Health 15 (8), 1747. doi: 10.3390/ijerph15081747
Fuhrman, B. J., Brinton, L. A., Pfeiffer, R. M., Xu, X., Veenstra, T. D., Teter, B. E., et al. (2012). Estrogen Metabolism and Mammographic Density in Postmenopausal Women: A Cross-Sectional Study. Cancer Epidemiol. Biomarkers Prev. 21, 1582–1591. doi: 10.1158/1055-9965.EPI-12-0247
Ghadge, A. G., Dasari, P., Stone, J., Thompson, E. W., Robker, R. L., Ingman, W. V. (2020). Pubertal Mammary Gland Development is a Key Determinant of Adult Mammographic Density. Semin. Cell Dev. Biol. 114, 143–158. doi: 10.1016/j.semcdb.2020.11.011
Gierach, G. L., Patel, D. A., Falk, R. T., Pfeiffer, R. M., Geller, B. M., Vacek, P. M., et al. (2015). Relationship of Serum Estrogens and Metabolites With Area and Volume Mammographic Densities. Horm. Cancer 6, 107–119. doi: 10.1007/s12672-015-0216-3
Goedert, J. J., Jones, G., Hua, X., Xu, X., Yu, G., Flores, R., et al. (2015). Investigation of the Association Between the Fecal Microbiota and Breast Cancer in Postmenopausal Women: A Population-Based Case-Control Pilot Study. J. Natl. Cancer Inst. 107 (8). doi: 10.1093/jnci/djv147
Harris, H. R., Tamimi, R. M., Willett, W. C., Hankinson, S. E., Michels, K. B. (2011). Body Size Across the Life Course, Mammographic Density, and Risk of Breast Cancer. Am. J. Epidemiol. 174, 909–918. doi: 10.1093/aje/kwr225
Jones, G. S., Feigelson, H. S., Falk, R. T., Hua, X., Ravel, J., Yu, G., et al. (2019). Mammographic Breast Density and its Association With Urinary Estrogens and the Fecal Microbiota in Postmenopausal Women. PloS One 14, e0216114. doi: 10.1371/journal.pone.0216114
Laborda-Illanes, A., Sanchez-Alcoholado, L., Dominguez-Recio, M. E., Jimenez-Rodriguez, B., Lavado, R., Comino-Méndez, I., et al. (2020). Breast and Gut Microbiota Action Mechanisms in Breast Cancer Pathogenesis and Treatment. Cancers 12, 2465. doi: 10.3390/cancers12092465
Ley, R. E., Turnbaugh, P. J., Klein, S., Gordon, J. I. (2006). Microbial Ecology: Human Gut Microbes Associated With Obesity. Nature 444, 1022–1023. doi: 10.1038/4441022a
Lozupone, C. A., Hamady, M., Kelley, S. T., Knight, R. (2007). Quantitative and Qualitative β Diversity Measures Lead to Different Insights Into Factors That Structure Microbial Communities. Appl. Environ. Microbiol. 73, 1576–1585. doi: 10.1128/AEM.01996-06
Lozupone, C. A., Stombaugh, J. I., Gordon, J. I., Jansson, J. K., Knight, R. (2012). Diversity, Stability and Resilience of the Human Gut Microbiota. Nature 489, 220–230. doi: 10.1038/nature11550
Luu, T. H., Michel, C., Bard, J.-M., Dravet, F., Nazih, H., Bobin-Dubigeon, C. (2017). Intestinal Proportion of Blautia Sp. Is Associated With Clinical Stage and Histoprognostic Grade in Patients With Early-Stage Breast Cancer. Nutr. Cancer 69, 267–275. doi: 10.1080/01635581.2017.1263750
Mallick, H., Rahnavard, A., McIver, L. J., Ma, S., Zhang, Y., Nguyen, L. H., et al. (2021). Multivariable Association Discovery in Population-Scale Meta-Omics Studies. bioRxiv 01.20.427420. doi: 10.1101/2021.01.20.427420
Martin, L. J., Boyd, N. F. (2008). Mammographic Density. Potential Mechanisms of Breast Cancer Risk Associated With Mammographic Density: Hypotheses Based on Epidemiological Evidence. Breast Cancer Res. 10, 201. doi: 10.1186/bcr1831
Martin, A. M., Sun, E. W., Rogers, G. B., Keating, D. J. (2019). The Influence of the Gut Microbiome on Host Metabolism Through the Regulation of Gut Hormone Release. Front. Physiol. 10, 428. doi: 10.3389/fphys.2019.00428
McArdle, B. H., Anderson, M. J. (2001). Fitting Multivariate Models to Community Data: A Comment on Distance-Based Redundancy Analysis. Ecology 82, 290–297. doi: 10.1890/0012-9658(2001)082[0290:FMMTCD]2.0.CO;2
McCarthy, H. D., Cole, T. J., Fry, T., Jebb, S. A., Prentice, A. M. (2006). Body Fat Reference Curves for Children. Int. J. Obes. 30, 598–602. doi: 10.1038/sj.ijo.0803232
McKnight, D. T., Huerlimann, R., Bower, D. S., Schwarzkopf, L., Alford, R. A., Zenger, K. R. (2019). Methods for Normalizing Microbiome Data: An Ecological Perspective. Methods Ecol. Evol. 10, 389–400. doi: 10.1111/2041-210X.13115
McMurdie, P. J., Holmes, S. (2013). Phyloseq: An R Package for Reproducible Interactive Analysis and Graphics of Microbiome Census Data. PloS One 8, e61217. doi: 10.1371/journal.pone.0061217
Novotny, R., Daida, Y., Morimoto, Y., Shepherd, J., Maskarinec, G. (2011). Puberty, Body Fat, and Breast Density in Girls of Several Ethnic Groups. Am. J. Hum. Biol. 23, 359–365. doi: 10.1002/ajhb.21145
Oksanen, J., Kindt, R., Legendre, P., O’Hara, B., Stevens, M. H. (2007). The Vegan Package. Community Ecol. Package 10, 719.
Pereira, A., Corvalán, C., Uauy, R., Klein, K. O., Mericq, V. (2015). Ultrasensitive Estrogen Levels at 7 Years of Age Predict Earlier Thelarche: Evidence From Girls of the Growth and Obesity Chilean Cohort. Eur. J. Endocrinol. 173, 835–842. doi: 10.1530/EJE-15-0327
Pollock, J., Glendinning, L., Wisedchanwet, T., Watson, M. (2018). The Madness of Microbiome: Attempting To Find Consensus “Best Practice” for 16S Microbiome Studies. Appl. Environ. Microbiol. 84 (7), e02627–17. doi: 10.1128/AEM.02627-17
Quast, C., Pruesse, E., Yilmaz, P., Gerken, J., Schweer, T., Yarza, P., et al. (2013). The SILVA Ribosomal RNA Gene Database Project: Improved Data Processing and Web-Based Tools. Nucleic Acids Res. 41, D590–D596. doi: 10.1093/nar/gks1219
Ricci, V., Carcione, D., Messina, S., Colombo, G. I., D’Alessandra, Y. (2020). Circulating 16s RNA in Biofluids: Extracellular Vesicles as Mirrors of Human Microbiome? Int. J. Mol. Sci. 21, 8959. doi: 10.3390/ijms21238959
Ruder, E. H., Dorgan, J. F., Kranz, S., Kris-Etherton, P. M., Hartman, T. J. (2008). Examining Breast Cancer Growth and Lifestyle Risk Factors: Early Life, Childhood, and Adolescence. Clin. Breast Cancer 8, 334–342. doi: 10.3816/CBC.2008.n.038
Selber-Hnatiw, S., Rukundo, B., Ahmadi, M., Akoubi, H., Al-Bizri, H., Aliu, A. F., et al. (2017). Human Gut Microbiota: Toward an Ecology of Disease. Front. Microbiol. 8, 1265. doi: 10.3389/fmicb.2017.01265
Sepich-Poore, G. D., Zitvogel, L., Straussman, R., Hasty, J., Wargo, J. A., Knight, R. (2021). The Microbiome and Human Cancer. Science 371 (6536), eabc4552. doi: 10.1126/science.abc4552
Shannon, C. E. (1948). A Mathematical Theory of Communication. Bell system Tech. J. 27, 379–423. doi: 10.1002/j.1538-7305.1948.tb01338.x
Shepherd, J. A., Herve, L., Landau, J., Fan, B., Kerlikowske, K., Cummings, S. R. (2006). Clinical Comparison of a Novel Breast DXA Technique to Mammographic Density. Med. Phys. 33, 1490–1498. doi: 10.1118/1.2193691
Shepherd, J. A., Malkov, S., Fan, B., Laidevant, A., Novotny, R., Maskarinec, G. (2008). Breast Density Assessment in Adolescent Girls Using DXA: A Feasibility Study. Cancer Epidemiol. Biomarkers Prev. 17, 1709–1713. doi: 10.1158/1055-9965.EPI-08-0006
Turnbaugh, P. J., Ley, R. E., Mahowald, M. A., Magrini, V., Mardis, E. R., Gordon, J. I. (2006). An Obesity-Associated Gut Microbiome With Increased Capacity for Energy Harvest. Nature 444, 1027–1031. doi: 10.1038/nature05414
van Buuren, S., Groothuis-Oudshoorn, K. (2011). Mice: Multivariate Imputation by Chained Equations in R. J. Stat. Software 45, 1–67. doi: 10.18637/jss.v045.i03
Wu, A. H., Tseng, C., Vigen, C., Yu, Y., Cozen, W., Garcia, A. A., et al. (2020). Gut Microbiome Associations With Breast Cancer Risk Factors and Tumor Characteristics: A Pilot Study. Breast Cancer Res. Treat 182, 451–463. doi: 10.1007/s10549-020-05702-6
Yaghjyan, L., Mai, V., Wang, X., Ukhanova, M., Tagliamonte, M., Martinez, Y. C., et al. (2021). Gut Microbiome, Body Weight, and Mammographic Breast Density in Healthy Postmenopausal Women. Cancer Causes Control 32 (7), 681–692. doi: 10.1007/s10552-021-01420-6
Keywords: adolescence, human, epidemiology, breast density, breast cancer, gut microbiota, 16s rRNA
Citation: Yoon LS, Jacobs JP, Hoehner J, Pereira A, Gana JC, Corvalán C and Michels KB (2021) The Association Between Breast Density and Gut Microbiota Composition at 2 Years Post-Menarche: A Cross-Sectional Study of Adolescents in Santiago, Chile. Front. Cell. Infect. Microbiol. 11:794610. doi: 10.3389/fcimb.2021.794610
Received: 13 October 2021; Accepted: 22 November 2021;
Published: 17 December 2021.
Edited by:
Wayne Young, AgResearch Ltd, New ZealandReviewed by:
Robert Kraaij, Erasmus Medical Center, NetherlandsBrandilyn Peters, Albert Einstein College of Medicine, United States
Copyright © 2021 Yoon, Jacobs, Hoehner, Pereira, Gana, Corvalán and Michels. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Karin B. Michels, k.michels@ucla.edu