1
|
LeMay-Nedjelski L, Yonemitsu C, Asbury MR, Butcher J, Ley SH, Hanley AJ, Kiss A, Unger S, Copeland JK, Wang PW, Stintzi A, Bode L, O'Connor DL. Oligosaccharides and Microbiota in Human Milk Are Interrelated at 3 Months Postpartum in a Cohort of Women with a High Prevalence of Gestational Impaired Glucose Tolerance. J Nutr 2021; 151:3431-3441. [PMID: 34510198 PMCID: PMC8562078 DOI: 10.1093/jn/nxab270] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/02/2021] [Revised: 06/07/2021] [Accepted: 07/22/2021] [Indexed: 12/15/2022] Open
Abstract
BACKGROUND Human milk is a rich source of human milk oligosaccharides (HMOs) and bacteria. It is unclear how these components interact within the breast microenvironment. OBJECTIVES The objectives were first, to investigate the association between maternal characteristics and HMOs, and second, to assess the association between HMOs and microbial community composition and predicted function in milk from women with high rates of gestational glucose intolerance. METHODS This was an exploratory analysis of a previously completed prospective cohort study (NCT01405547) where milk samples (n = 107) were collected at 3 mo postpartum. Milk microbiota composition was analyzed by V4-16S ribosomal RNA gene sequencing and HMOs by rapid high-throughput HPLC. Data were stratified and analyzed by maternal secretor status phenotype and associations between HMOs and microbiota were determined using linear regression models (ɑ-diversity), Adonis (B-diversity), Poisson regression models (differential abundance), and general linear models (predicted microbial function). RESULTS Prepregnancy BMI, race, and frequency of direct breastfeeding, but not gestational glucose intolerance, were found to be significantly associated with a number of HMOs among secretors and non-secretors. Fucosyllacto-N-hexaose was negatively associated with microbial richness (Chao1) among secretors [B-estimate (SE): -9.3 × 102 (3.4 × 102); P = 0.0082] and difucosyllacto-N-hexaose was negatively associated with microbiota diversity (Shannon index) [-1.7 (0.78); P = 0.029] among secretors. Lacto-N-neotetraose (LNnT) was associated with both microbial B-diversity (weighted UniFrac R2 = 0.040, P = 0.036) and KEGG ortholog B-diversity (Bray-Curtis R2 = 0.039, P = 0.043) in secretors. Additionally, difucosyllactose in secretors and disialyllacto-N-hexaose and LNnT in non-secretors were associated with enrichment of predicted microbial genes encoding for metabolism- and infection-related pathways (P-false discovery rate < 0.1). CONCLUSIONS HMOs are associated with the microbial composition and predicted microbial functions in human milk at 3 mo postpartum. Further research is needed to investigate the role these relations play in maternal and infant health.
Collapse
Affiliation(s)
- Lauren LeMay-Nedjelski
- Department of Nutritional Sciences, University of Toronto, Toronto, Ontario, Canada
- Translational Medicine Program, The Hospital for Sick Children, Toronto, Ontario, Canada
| | - Chloe Yonemitsu
- Department of Pediatrics and Mother-Milk-Infant Center of Research Excellence (MOMI CORE), University of California, San Diego, La Jolla, CA, USA
| | - Michelle R Asbury
- Department of Nutritional Sciences, University of Toronto, Toronto, Ontario, Canada
- Translational Medicine Program, The Hospital for Sick Children, Toronto, Ontario, Canada
| | - James Butcher
- Ottawa Institute of Systems Biology and Department of Biochemistry, Microbiology, and Immunology, University of Ottawa, Ottawa, Ontario, Canada
| | - Sylvia H Ley
- Department of Epidemiology, Tulane University School of Public Health and Tropical Medicine, New Orleans, LA, USA
| | - Anthony J Hanley
- Department of Nutritional Sciences, University of Toronto, Toronto, Ontario, Canada
| | - Alex Kiss
- Department of Research Design and Biostatistics, Sunnybrook Research Institute, Toronto, Ontario, Canada
| | - Sharon Unger
- Department of Nutritional Sciences, University of Toronto, Toronto, Ontario, Canada
- Translational Medicine Program, The Hospital for Sick Children, Toronto, Ontario, Canada
- Department of Pediatrics, University of Toronto, Toronto, Ontario, Canada
- Department of Pediatrics, Sinai Health, Toronto, Ontario, Canada
- Lunenfeld-Tanenbaum Research Institute, Sinai Health, Toronto, Ontario, Canada
| | - Julia K Copeland
- Centre for the Analysis of Genome Evolution and Function, University of Toronto, Ontario, Canada
| | - Pauline W Wang
- Centre for the Analysis of Genome Evolution and Function, University of Toronto, Ontario, Canada
| | - Alain Stintzi
- Ottawa Institute of Systems Biology and Department of Biochemistry, Microbiology, and Immunology, University of Ottawa, Ottawa, Ontario, Canada
| | - Lars Bode
- Department of Pediatrics and Mother-Milk-Infant Center of Research Excellence (MOMI CORE), University of California, San Diego, La Jolla, CA, USA
| | - Deborah L O'Connor
- Department of Nutritional Sciences, University of Toronto, Toronto, Ontario, Canada
- Translational Medicine Program, The Hospital for Sick Children, Toronto, Ontario, Canada
| |
Collapse
|
3
|
Fumagalli M, Cagliani R, Pozzoli U, Riva S, Comi GP, Menozzi G, Bresolin N, Sironi M. Widespread balancing selection and pathogen-driven selection at blood group antigen genes. Genes Dev 2009; 19:199-212. [PMID: 18997004 PMCID: PMC2652214 DOI: 10.1101/gr.082768.108] [Citation(s) in RCA: 116] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2008] [Accepted: 11/04/2008] [Indexed: 12/31/2022]
Abstract
Historically, allelic variations in blood group antigen (BGA) genes have been regarded as possible susceptibility factors for infectious diseases. Since host-pathogen interactions are major determinants in evolution, BGAs can be thought of as selection targets. In order to verify this hypothesis, we obtained an estimate of pathogen richness for geographic locations corresponding to 52 populations distributed worldwide; after correction for multiple tests and for variables different from selective forces, significant correlations with pathogen richness were obtained for multiple variants at 11 BGA loci out of 26. In line with this finding, we demonstrate that three BGA genes, namely CD55, CD151, and SLC14A1, have been subjected to balancing selection, a process, rare outside MHC genes, which maintains variability at a locus. Moreover, we identified a gene region immediately upstream the transcription start site of FUT2 which has undergone non-neutral evolution independently from the coding region. Finally, in the case of BSG, we describe the presence of a highly divergent haplotype clade and the possible reasons for its maintenance, including frequency-dependent balancing selection, are discussed. These data indicate that BGAs have been playing a central role in the host-pathogen arms race during human evolutionary history and no other gene category shows similar levels of widespread selection, with the only exception of loci involved in antigen recognition.
Collapse
Affiliation(s)
- Matteo Fumagalli
- Scientific Institute IRCCS E. Medea, Bioinformatic Lab, 23842 Bosisio Parini (LC), Italy
- Bioengineering Department, Politecnico di Milano, 20133 Milan, Italy
| | - Rachele Cagliani
- Scientific Institute IRCCS E. Medea, Bioinformatic Lab, 23842 Bosisio Parini (LC), Italy
| | - Uberto Pozzoli
- Scientific Institute IRCCS E. Medea, Bioinformatic Lab, 23842 Bosisio Parini (LC), Italy
| | - Stefania Riva
- Scientific Institute IRCCS E. Medea, Bioinformatic Lab, 23842 Bosisio Parini (LC), Italy
| | - Giacomo P. Comi
- Dino Ferrari Centre, Department of Neurological Sciences, University of Milan, IRCCS Ospedale Maggiore Policlinico, Mangiagalli and Regina Elena Foundation, 20100 Milan, Italy
| | - Giorgia Menozzi
- Scientific Institute IRCCS E. Medea, Bioinformatic Lab, 23842 Bosisio Parini (LC), Italy
| | - Nereo Bresolin
- Scientific Institute IRCCS E. Medea, Bioinformatic Lab, 23842 Bosisio Parini (LC), Italy
- Dino Ferrari Centre, Department of Neurological Sciences, University of Milan, IRCCS Ospedale Maggiore Policlinico, Mangiagalli and Regina Elena Foundation, 20100 Milan, Italy
| | - Manuela Sironi
- Scientific Institute IRCCS E. Medea, Bioinformatic Lab, 23842 Bosisio Parini (LC), Italy
| |
Collapse
|