From Microbes to Megafauna: The Ecology of Gut Microbiomes in Wild Herbivores

All mammals host communities of commensal microbes in and on their bodies. Recent technological advances, combined with experimental studies in laboratory animals, are beginning to reveal the ubiquitous links between the gut microbiome and host disease, metabolism, immunity, and numerous other host functions. A new challenge of microbiome research is to understand how microbiomes function in natural ecosystems and within free-living hosts, in order to expand our understanding of the ecological principles underlying microbiome-host dynamics, and to facilitate applications that are relevant for population health. Here I describe microbiome diversity and dynamics across multiple ecological scales in three wild herbivore populations. I elucidate mechanisms by which microbiome diversity is maintained in a population throughout temporal resource fluctuations, and identify relationships between diet, disease, and the microbiome in a population of African buffalo. At a larger spatial scale, I describe microbiome variation across a metapopulation of bighorn sheep, and identify correlations with geography, environment, and genetic structure. Finally, I demonstrate the potential utility of microbiome surveillance for wildlife conservation and management by elucidating the impacts of supplemental feeding on the gut microbiomes of elk in the Greater Yellowstone Ecosystem. Taken together, this body of work contributes toward clarifying the ecological processes that shape the microbiomes of free-living hosts across multiple ecological scales. The patterns and processes identified in this dissertation lay the groundwork for understanding how host-associated microbiome communities may respond to changing environmental conditions, disease, and habitat alterations

Serum biochemistry panels in African buffalo: Defining reference intervals and assessing variability across season, age and sex

Serum biochemical parameters can be utilized to evaluate the physiological status of an animal, and relate it to the animal’s health. In order to accurately interpret individual animal biochemical results, species-specific reference intervals (RI) must be established. Reference intervals for biochemical parameters differ between species, and physiological differences including reproductive status, nutritional resource availability, disease status, and age affect parameters within the same species. The objectives of this study were to (1) establish RI for biochemical parameters in managed African buffalo (Syncerus caffer), (2) assess the effects of age, sex, pregnancy, and season on serum biochemistry values, and (3) compare serum biochemistry values from a managed herd to a free-ranging buffalo herd and to values previously published for captive (zoo) buffalo. Season profoundly affected all biochemistry parameters, possibly due to changes in nutrition and disease exposure. Age also affected all biochemical parameters except gamma glutamyl transferase and magnesium, consistent with patterns seen in cattle. Sex and reproductive status had no detectable effects on the parameters that were measured. The biochemical profiles of managed buffalo were distinct from those observed in the free-ranging herd and captive buffalo. Biochemical differences between buffalo from captive, managed, and free-ranging populations may be related to nutritional restriction or lack of predation in the context of management or captivity. The reference intervals provided in this study, in addition to the seasonal and age-related patterns observed, provide a foundation for health investigations that may inform management strategies in this ecologically and economically important species

Fecal Methylmercury Correlates With Gut Microbiota Taxa in Pacific Walruses (Odobenus rosmarus divergens)

ObjectivesMethylmercury metabolism was investigated in Pacific walruses (Odobenus rosmarus divergens) from St. Lawrence Island, Alaska, United States.MethodsTotal mercury and methylmercury concentrations were measured in fecal samples and paired colon samples (n = 16 walruses). Gut microbiota composition and diversity were determined using 16S rRNA gene sequencing. Associations between fecal and colon mercury and the 24 most prevalent gut microbiota taxa were investigated using linear models.ResultsIn fecal samples, the median values for total mercury, methylmercury, and %methylmercury (of total mercury) were 200 ng/g, 4.7 ng/g, and 2.5%, respectively, while in colon samples, the median values for the same parameters were 28 ng/g, 7.8 ng/g, and 26%, respectively. In fecal samples, methylmercury was negatively correlated with one Bacteroides genus, while members of the Oscillospirales order were positively correlated with both methylmercury and %methylmercury (of total mercury). In colon samples, %methylmercury (of total mercury) was negatively correlated with members of two genera, Romboutsia and Paeniclostridium.ConclusionsMedian %methylmercury (of total mercury) was 10 times higher in the colon compared to the fecal samples, suggesting that methylmercury was able to pass through the colon into systemic circulation. Fecal total mercury and/or methylmercury concentrations in walruses were comparable to some human studies despite differences in seafood consumption rates, suggesting that walruses excreted less mercury. There are no members (at this time) of the Oscillospirales order which are known to contain the genes to methylate mercury, suggesting the source of methylmercury in the gut was from diet and not in vivo methylation

Genome-wide interaction analysis of menopausal hormone therapy use and breast cancer risk among 62,370 women

Use of menopausal hormone therapy (MHT) is associated with increased risk for breast cancer. However, the relevant mechanisms and its interaction with genetic variants are not fully understood. We conducted a genome-wide interaction analysis between MHT use and genetic variants for breast cancer risk in 27,585 cases and 34,785 controls from 26 observational studies. All women were post-menopausal and of European ancestry. Multivariable logistic regression models were used to test for multiplicative interactions between genetic variants and current MHT use. We considered interaction p-values < 5 × 10–8 as genome-wide significant, and p-values < 1 × 10–5 as suggestive. Linkage disequilibrium (LD)-based clumping was performed to identify independent candidate variants. None of the 9.7 million genetic variants tested for interactions with MHT use reached genome-wide significance. Only 213 variants, representing 18 independent loci, had p-values < 1 × 105. The strongest evidence was found for rs4674019 (p-value = 2.27 × 10–7), which showed genome-wide significant interaction (p-value = 3.8 × 10–8) with current MHT use when analysis was restricted to population-based studies only. Limiting the analyses to combined estrogen–progesterone MHT use only or to estrogen receptor (ER) positive cases did not identify any genome-wide significant evidence of interactions. In this large genome-wide SNP-MHT interaction study of breast cancer, we found no strong support for common genetic variants modifying the effect of MHT on breast cancer risk. These results suggest that common genetic variation has limited impact on the observed MHT–breast cancer risk association

Five endometrial cancer risk loci identified through genome-wide association analysis.

We conducted a meta-analysis of three endometrial cancer genome-wide association studies (GWAS) and two follow-up phases totaling 7,737 endometrial cancer cases and 37,144 controls of European ancestry. Genome-wide imputation and meta-analysis identified five new risk loci of genome-wide significance at likely regulatory regions on chromosomes 13q22.1 (rs11841589, near KLF5), 6q22.31 (rs13328298, in LOC643623 and near HEY2 and NCOA7), 8q24.21 (rs4733613, telomeric to MYC), 15q15.1 (rs937213, in EIF2AK4, near BMF) and 14q32.33 (rs2498796, in AKT1, near SIVA1). We also found a second independent 8q24.21 signal (rs17232730). Functional studies of the 13q22.1 locus showed that rs9600103 (pairwise r(2) = 0.98 with rs11841589) is located in a region of active chromatin that interacts with the KLF5 promoter region. The rs9600103[T] allele that is protective in endometrial cancer suppressed gene expression in vitro, suggesting that regulation of the expression of KLF5, a gene linked to uterine development, is implicated in tumorigenesis. These findings provide enhanced insight into the genetic and biological basis of endometrial cancer.I.T. is supported by Cancer Research UK and the Oxford Comprehensive Biomedical Research Centre. T.H.T.C. is supported by the Rhodes Trust and the Nuffield Department of Medicine. Funding for iCOGS infrastructure came from the European Community's Seventh Framework Programme under grant agreement 223175 (HEALTH-F2-2009-223175) (COGS), Cancer Research UK (C1287/A10118, C1287/A10710, C12292/A11174, C1281/A12014, C5047/A8384, C5047/A15007, C5047/A10692 and C8197/A16565), the US National Institutes of Health (R01 CA128978, U19 CA148537, U19 CA148065 and U19 CA148112), the US Department of Defense (W81XWH-10-1-0341), the Canadian Institutes of Health Research (CIHR) for the CIHR Team in Familial Risks of Breast Cancer, the Susan G. Komen Foundation for the Cure, the Breast Cancer Research Foundation and the Ovarian Cancer Research Fund. SEARCH recruitment was funded by a programme grant from Cancer Research UK (C490/A10124). Stage 1 and stage 2 case genotyping was supported by the NHMRC (552402 and 1031333). Control data were generated by the WTCCC, and a full list of the investigators who contributed to the generation of the data is available from the WTCCC website. We acknowledge use of DNA from the British 1958 Birth Cohort collection, funded by UK Medical Research Council grant G0000934 and Wellcome Trust grant 068545/Z/02; funding for this project was provided by the Wellcome Trust under award 085475. NSECG was supported by the European Union's Framework Programme 7 CHIBCHA grant and Wellcome Trust Centre for Human Genetics Core Grant 090532/Z/09Z, and CORGI was funded by Cancer Research UK. BCAC is funded by Cancer Research UK (C1287/A10118 and C1287/A12014). OCAC is supported by a grant from the Ovarian Cancer Research Fund thanks to donations by the family and friends of Kathryn Sladek Smith (PPD/RPCI.07) and the UK National Institute for Health Research Biomedical Research Centres at the University of Cambridge.This is the author accepted manuscript. The final version is available from Nature Publishing Group via http://dx.doi.org/10.1038/ng.356

Polygenic Risk Scores for Prediction of Breast Cancer and Breast Cancer Subtypes

Stratification of women according to their risk of breast cancer based on polygenic risk scores (PRSs) could improve screening and prevention strategies. Our aim was to develop PRSs, optimized for prediction of estrogen receptor (ER)-specific disease, from the largest available genome-wide association dataset and to empirically validate the PRSs in prospective studies. The development dataset comprised 94,075 case subjects and 75,017 control subjects of European ancestry from 69 studies, divided into training and validation sets. Samples were genotyped using genome-wide arrays, and single-nucleotide polymorphisms (SNPs) were selected by stepwise regression or lasso penalized regression. The best performing PRSs were validated in an independent test set comprising 11,428 case subjects and 18,323 control subjects from 10 prospective studies and 190,040 women from UK Biobank (3,215 incident breast cancers). For the best PRSs (313 SNPs), the odds ratio for overall disease per 1 standard deviation in ten prospective studies was 1.61 (95%CI: 1.57-1.65) with area under receiver-operator curve (AUC) = 0.630 (95%CI: 0.628-0.651). The lifetime risk of overall breast cancer in the top centile of the PRSs was 32.6%. Compared with women in the middle quintile, those in the highest 1% of risk had 4.37- and 2.78-fold risks, and those in the lowest 1% of risk had 0.16- and 0.27-fold risks, of developing ER-positive and ER-negative disease, respectively. Goodness-of-fit tests indicated that this PRS was well calibrated and predicts disease risk accurately in the tails of the distribution. This PRS is a powerful and reliable predictor of breast cancer risk that may improve breast cancer prevention programs.NovartisEli Lilly and CompanyAstraZenecaAbbViePfizer UKCelgeneEisaiGenentechMerck Sharp and DohmeRocheCancer Research UKGovernment of CanadaArray BioPharmaGenome CanadaNational Institutes of HealthEuropean CommissionMinistère de l'Économie, de l’Innovation et des Exportations du QuébecSeventh Framework ProgrammeCanadian Institutes of Health Researc

Genome-wide association analysis of more than 120,000 individuals identifies 15 new susceptibility loci for breast cancer.

Genome-wide association studies (GWAS) and large-scale replication studies have identified common variants in 79 loci associated with breast cancer, explaining ∼14% of the familial risk of the disease. To identify new susceptibility loci, we performed a meta-analysis of 11 GWAS, comprising 15,748 breast cancer cases and 18,084 controls together with 46,785 cases and 42,892 controls from 41 studies genotyped on a 211,155-marker custom array (iCOGS). Analyses were restricted to women of European ancestry. We generated genotypes for more than 11 million SNPs by imputation using the 1000 Genomes Project reference panel, and we identified 15 new loci associated with breast cancer at P < 5 × 10(-8). Combining association analysis with ChIP-seq chromatin binding data in mammary cell lines and ChIA-PET chromatin interaction data from ENCODE, we identified likely target genes in two regions: SETBP1 at 18q12.3 and RNF115 and PDZK1 at 1q21.1. One association appears to be driven by an amino acid substitution encoded in EXO1.BCAC is funded by Cancer Research UK (C1287/A10118, C1287/A12014) and by the European Community's Seventh Framework Programme under grant agreement 223175 (HEALTH-F2-2009-223175) (COGS). Meetings of the BCAC have been funded by the European Union COST programme (BM0606). Genotyping on the iCOGS array was funded by the European Union (HEALTH-F2-2009-223175), Cancer Research UK (C1287/A10710, C8197/A16565), the Canadian Institutes of Health Research (CIHR) for the CIHR Team in Familial Risks of Breast Cancer program and the Ministry of Economic Development, Innovation and Export Trade of Quebec, grant PSR-SIIRI-701. Combination of the GWAS data was supported in part by the US National Institutes of Health (NIH) Cancer Post-Cancer GWAS initiative, grant 1 U19 CA148065-01 (DRIVE, part of the GAME-ON initiative). For a full description of funding and acknowledgments, see the Supplementary Note.This is the author accepted manuscript. The final version is available from NPG via http://dx.doi.org/10.1038/ng.324

Genome-wide association study identifies 32 novel breast cancer susceptibility loci from overall and subtype-specific analyses.

Breast cancer susceptibility variants frequently show heterogeneity in associations by tumor subtype1-3. To identify novel loci, we performed a genome-wide association study including 133,384 breast cancer cases and 113,789 controls, plus 18,908 BRCA1 mutation carriers (9,414 with breast cancer) of European ancestry, using both standard and novel methodologies that account for underlying tumor heterogeneity by estrogen receptor, progesterone receptor and human epidermal growth factor receptor 2 status and tumor grade. We identified 32 novel susceptibility loci (P < 5.0 × 10-8), 15 of which showed evidence for associations with at least one tumor feature (false discovery rate < 0.05). Five loci showed associations (P < 0.05) in opposite directions between luminal and non-luminal subtypes. In silico analyses showed that these five loci contained cell-specific enhancers that differed between normal luminal and basal mammary cells. The genetic correlations between five intrinsic-like subtypes ranged from 0.35 to 0.80. The proportion of genome-wide chip heritability explained by all known susceptibility loci was 54.2% for luminal A-like disease and 37.6% for triple-negative disease. The odds ratios of polygenic risk scores, which included 330 variants, for the highest 1% of quantiles compared with middle quantiles were 5.63 and 3.02 for luminal A-like and triple-negative disease, respectively. These findings provide an improved understanding of genetic predisposition to breast cancer subtypes and will inform the development of subtype-specific polygenic risk scores

An original phylogenetic approach identified mitochondrial haplogroup T1a1 as inversely associated with breast cancer risk in BRCA2 mutation carriers

Introduction: Individuals carrying pathogenic mutations in the BRCA1 and BRCA2 genes have a high lifetime risk of breast cancer. BRCA1 and BRCA2 are involved in DNA double-strand break repair, DNA alterations that can be caused by exposure to reactive oxygen species, a main source of which are mitochondria. Mitochondrial genome variations affect electron transport chain efficiency and reactive oxygen species production. Individuals with different mitochondrial haplogroups differ in their metabolism and sensitivity to oxidative stress. Variability in mitochondrial genetic background can alter reactive oxygen species production, leading to cancer risk. In the present study, we tested the hypothesis that mitochondrial haplogroups modify breast cancer risk in BRCA1/2 mutation carriers. Methods: We genotyped 22,214 (11,421 affected, 10,793 unaffected) mutation carriers belonging to the Consortium of Investigators of Modifiers of BRCA1/2 for 129 mitochondrial polymorphisms using the iCOGS array. Haplogroup inference and association detection were performed using a phylogenetic approach. ALTree was applied to explore the reference mitochondrial evolutionary tree and detect subclades enriched in affected or unaffected individuals. Results: We discovered that subclade T1a1 was depleted in affected BRCA2 mutation carriers compared with the rest of clade T (hazard ratio (HR) = 0.55; 95% confidence interval (CI), 0.34 to 0.88; P = 0.01). Compared with the most frequent haplogroup in the general population (that is, H and T clades), the T1a1 haplogroup has a HR of 0.62 (95% CI, 0.40 to 0.95; P = 0.03). We also identified three potential susceptibility loci, including G13708A/rs28359178, which has demonstrated an inverse association with familial breast cancer risk. Conclusions: This study illustrates how original approaches such as the phylogeny-based method we used can empower classical molecular epidemiological studies aimed at identifying association or risk modification effects.Peer reviewe

Fine-Scale Mapping of the 4q24 Locus Identifies Two Independent Loci Associated with Breast Cancer Risk

Background: A recent association study identified a common variant (rs9790517) at 4q24 to be associated with breast cancer risk. Independent association signals and potential functional variants in this locus have not been explored. Methods: We conducted a fine-mapping analysis in 55,540 breast cancer cases and 51,168 controls from the Breast Cancer Association Consortium. Results: Conditional analyses identified two independent association signals among women of European ancestry, represented by rs9790517 [conditional P = 2.51 × 10−4; OR, 1.04; 95% confidence interval (CI), 1.02–1.07] and rs77928427 (P = 1.86 × 10−4; OR, 1.04; 95% CI, 1.02–1.07). Functional annotation using data from the Encyclopedia of DNA Elements (ENCODE) project revealed two putative functional variants, rs62331150 and rs73838678 in linkage disequilibrium (LD) with rs9790517 (r2 ≥ 0.90) residing in the active promoter or enhancer, respectively, of the nearest gene, TET2. Both variants are located in DNase I hypersensitivity and transcription factor–binding sites. Using data from both The Cancer Genome Atlas (TCGA) and Molecular Taxonomy of Breast Cancer International Consortium (METABRIC), we showed that rs62331150 was associated with level of expression of TET2 in breast normal and tumor tissue. Conclusion: Our study identified two independent association signals at 4q24 in relation to breast cancer risk and suggested that observed association in this locus may be mediated through the regulation of TET2. Impact: Fine-mapping study with large sample size warranted for identification of independent loci for breast cancer risk