Publications

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

Posted on by Sylva Donaldson in Publications

Lauren LeMay-Nedjelski, Chloe Yonemitsu, Michelle R Asbury, James Butcher, Sylvia H Ley, Anthony J Hanley, Alex Kiss, Sharon Unger, Julia K Copeland, Pauline W Wang, Alain Stintzi, Lars Bode, and Deborah L O’Connor

J Nutr 2021 Nov 2;151(11):3431-3441. doi: 10.1093/jn/nxab270.

PMID: 34510198

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.

The Impact of Migration on the Gut Metagenome of South Asian Canadians

Posted on by Sylva Donaldson in Publications

Julia K. Copeland, Gary Chao, Shelley Vanderhout, Erica Acton, Pauline W. Wang, Eric I. Benchimol, Ahmed El-Sohemy, Ken Croitoru, Jennifer L. Gommerman, David S. Guttman & the GEMINI Research Team

Gut Microbes 2021 Jan-Dec;13(1):1-29. doi: 10.1080/19490976.2021.1902705.

PMID: 33794735

Abstract

South Asian (SA) Canadian immigrants have a higher risk of developing certain immune-mediated inflammatory diseases compared to non-migrant SAs. We sought to investigate the effect of migration on the gut metagenome and to identify microbiological associations between migration and conditions that may influence the development of immune-mediated inflammatory diseases. Metagenomic analysis of 58 first-generation (GEN1) SA immigrants and 38 unrelated Canadian born children-of-immigrants (GEN2) determined that the time lived in Canada was associated with continued changes in gut microbial communities. Migration of GEN1 to Canada early in life results in a gut community with similarities to GEN2 SA Canadians and non-SA North Americans. Conversely, GEN1 immigrants who arrived recently to Canada exhibited pronounced differences from GEN2, while displaying microbial similarities to a non-migrating SA cohort. Multivariate analysis identified that community composition was primarily influenced by high abundance taxa. Prevotella copri dominated in GEN1 and non-migrant SAs. Clostridia and functionally related Bacteroidia spp. replaced P. copri dominance over generations in Canada. Mutually exclusive Dialister species occurred at differing relative abundances over time and generations in Canada. This shift in species composition is accompanied by a change in genes associated with carbohydrate utilization and short-chain fatty acid production. Total energy derived from carbohydrates compared to protein consumption was significantly higher for GEN1 recent immigrants, which may influence the functional requirements of the gut community. This study demonstrates the associations between migration and the gut microbiome, which may be further associated with the altered risk of immune-mediated inflammatory diseases observed for SA Canadians.

Assessment of Inter-Laboratory Variation in the Characterization and Analysis of the Mucosal Microbiota in Crohn’s Disease and Ulcerative Colitis

Posted on by Sylva Donaldson in Publications

Szamosi Jake C., Forbes Jessica D., Copeland Julia K., Knox Natalie C., Shekarriz Shahrokh, Rossi Laura, Graham Morag, Bonner Christine, Guttman David S., Van Domselaar Gary, Surette Michael G., Bernstein Charles N.

Front. Microbiol., 21 August 2020 Sec. Systems Microbiology Volume 11 – 2020. doi: 10.3389/fmicb.2020.02028

PMID: 32973734

Background: In studies evaluating the microbiome, numerous factors can contribute to technical variability. These factors include DNA extraction methodology, sequencing protocols, and data analysis strategies. We sought to evaluate the impact these factors have on the results obtained when the sequence data are independently generated and analyzed by different laboratories.

Methods: To evaluate the effect of technical variability, we used human intestinal biopsy samples resected from individuals diagnosed with an inflammatory bowel disease (IBD), including Crohn’s disease (n = 12) and ulcerative colitis (n = 10), and those without IBD (n = 10). Matched samples from each participant were sent to three laboratories and studied using independent protocols for DNA extraction, library preparation, targeted-amplicon sequencing of a 16S rRNA gene hypervariable region, and processing of sequence data. We looked at two measures of interest – Bray–Curtis PERMANOVA R2 values and log2 fold-change estimates of the 25 most-abundant taxa – to assess variation in the results produced by each laboratory, as well the relative contribution to variation from the different extraction, sequencing, and analysis steps used to generate these measures.

Results: The R2 values and estimated differential abundance associated with diagnosis were consistent across datasets that used different DNA extraction and sequencing protocols, and within datasets that pooled samples from multiple protocols; however, variability in bioinformatic processing of sequence data led to changes in R2 values and inconsistencies in taxonomic assignment and abundance estimates.

Conclusion: Although the contribution of DNA extraction and sequencing methods to variability were observable, we find that results can be robust to the various extraction and sequencing approaches used in our study. Differences in data processing methods have a larger impact on results, making comparison among studies less reliable and the combined analysis of bioinformatically processed samples nearly impossible. Our results highlight the importance of making raw sequence data available to facilitate combined and comparative analyses of published studies using common data processing protocols. Study methodologies should provide detailed data processing methods for validation, interpretability, reproducibility, and comparability.

Mycobiome Sequencing and Analysis Applied to Fungal Community Profiling of the Lower Respiratory Tract During Fungal Pathogenesis

Posted on by Sylva Donaldson in Publications

Lisa R. McTaggart, Julia K. Copeland, Anuradha Surendra, Pauline W. Wang, Shahid Husain, Bryan Coburn, David S. Guttman and Julianne V. Kus

Front Microbiol. 2019 Mar 15;10:512. doi: 10.3389/fmicb.2019.00512. eCollection 2019.

PMID: 30930884

Abstract

Invasive fungal infections are an increasingly important cause of human morbidity and mortality. We generated a next-generation sequencing (NGS)-based method designed to detect a wide range of fungi and applied it to analysis of the fungal microbiome (mycobiome) of the lung during fungal infection. Internal transcribed spacer 1 (ITS1) amplicon sequencing and a custom analysis pipeline detected 96% of species from three mock communities comprised of potential fungal lung pathogens with good recapitulation of the expected species distributions (Pearson correlation coefficients r = 0.63, p = 0.004; r = 0.71, p < 0.001; r = 0.62, p = 0.002). We used this pipeline to analyze mycobiomes of bronchoalveolar lavage (BAL) specimens classified as culture-negative (n = 50) or culture-positive (n = 39) for Blastomyces dermatitidis/gilchristii, the causative agent of North America blastomycosis. Detected in 91.4% of the culture-positive samples, Blastomyces dominated (>50% relative abundance) the mycobiome in 68.6% of these culture-positive samples but was absent in culture-negative samples. To overcome any bias in relative abundance due to between-sample variation in fungal biomass, an abundance-weighting calculation was used to normalize the data by accounting for sample-specific PCR cycle number and PCR product concentration data utilized during sample preparation. After normalization, there was a statistically significant greater overall abundance of ITS1 amplicon in the Blastomyces-culture-positive samples versus culture-negative samples. Moreover, the normalization revealed a greater biomass of yeast and environmental fungi in several Blastomyces-culture-positive samples than in the culture-negative samples. Successful detection of Coccidioides, Scedosporium, Phaeoacremonium, and Aspergillus in 6 additional culture-positive BALs by ITS1 amplicon sequencing demonstrates the ability of this method to detect a broad range of fungi from clinical specimens, suggesting that it may be a potentially useful adjunct to traditional fungal microbiological testing for the diagnosis of respiratory mycoses.

Comparison of Co-housing and Littermate Methods for Microbiota Standardization in Mouse Models

Posted on by biota in Publications

Susan J Robertson, Paul Lemire, Heather Maughan, Ashleigh Goethel, Williams Turpin, Larbi Bedrani, David S Guttman, Kenneth Croitoru, Stephen E Girardin, Dana J Philpott

Cell Rep. 2019 May 7;27(6):1910-1919.e2. doi: 10.1016/j.celrep.2019.04.023

PMID: 31067473

Abstract

The intestinal microbiota is a fundamental factor that broadly influences physiology. Thus, studies using transgenic animals should be designed to limit the confounding effects of microbiota variation between strains. Here, we report the impact on intestinal microbiota of co-housed versus F2-generation littermates, two commonly used techniques to standardize microbiota in animal models. Our results establish that while fecal microbiota is partially normalized by extended co-housing, mucosal communities associated with the proximal colon and terminal ileum remain stable and distinct. In contrast, strain inter-crossing to generate F2 littermates allows robust microbiota standardization in fecal, colon, and ileum sampling locations. Using reciprocal inter-crosses of P1 parents, we identify dissymmetry in F2 community structures caused by maternal transmission, in particular of the Verrucomicrobiaceae. Thus, F2 littermate animals from a unidirectional P1 cross should be used as a standard method to minimize the influence of the microbiota in genotype-phenotype studies.

Entropy and Information within Intrinsically Disordered Protein Regions

Posted on by biota in Publications

Iva Pritišanac, Robert M. Vernon, Alan M. Moses, Julie D. Forman Kay

Entropy 2019 21(7):622. doi: 10.3390/e21070662

PMID: 33267376

Abstract

Bioinformatics and biophysical studies of intrinsically disordered proteins and regions (IDRs) note the high entropy at individual sequence positions and in conformations sampled in solution. This prevents application of the canonical sequence-structure-function paradigm to IDRs and motivates the development of new methods to extract information from IDR sequences. We argue that the information in IDR sequences cannot be fully revealed through positional conservation, which largely measures stable structural contacts and interaction motifs. Instead, considerations of evolutionary conservation of molecular features can reveal the full extent of information in IDRs. Experimental quantification of the large conformational entropy of IDRs is challenging but can be approximated through the extent of conformational sampling measured by a combination of NMR spectroscopy and lower-resolution structural biology techniques, which can be further interpreted with simulations. Conformational entropy and other biophysical features can be modulated by post-translational modifications that provide functional advantages to IDRs by tuning their energy landscapes and enabling a variety of functional interactions and modes of regulation. The diverse mosaic of functional states of IDRs and their conformational features within complexes demands novel metrics of information, which will reflect the complicated sequence-conformational ensemble-function relationship of IDRs.

Limiting oxidative DNA damage reduces microbe-induced colitis-associated colorectal cancer.

Posted on by biota in Publications

Irrazabal T, Thakur BK, Kang M, Malaise Y, Streutker C, Wong EOY, Copeland J, Gryfe R, Guttman DS, Navarre WW, Martin A

Nature Communications 2020 11(1):1802. doi: 10.1038/s41467-020-15549-6

PMID: 32286276

Abstract

Inflammatory bowel disease patients have a greatly increased risk of developing colitis-associated colon cancer (CAC); however, the basis for inflammation-induced genetic damage requisite for neoplasia is unclear. Using three models of CAC, we find that sustained inflammation triggers 8-oxoguanine DNA lesions. Strikingly, antioxidants or iNOS inhibitors reduce 8-oxoguanine and polyps in CAC models. Because the mismatch repair (MMR) system repairs 8-oxoguanine and is frequently defective in colorectal cancer (CRC), we test whether 8-oxoguanine mediates oncogenesis in a Lynch syndrome (MMR-deficient) model. We show that microbiota generates an accumulation of 8-oxoguanine lesions in MMR-deficient colons. Accordingly, we find that 8-oxoguanine is elevated in neoplastic tissue of Lynch syndrome patients compared to matched untransformed tissue or non-Lynch syndrome neoplastic tissue. While antioxidants reduce 8-oxoguanine, they do not reduce CRC in Lynch syndrome models. Hence, microbe-induced oxidative/nitrosative DNA damage play causative roles in inflammatory CRC models, but not in Lynch syndrome models.

Tung Tree (Vernicia fordii) Genome Provides A Resource for Understanding Genome Evolution and Improved Oil Production.

Posted on by biota in Publications

Zhang L, Liu M, Long H, Dong W, Pasha A, Esteban E, Li W, Yang X, Li Z, Song A, Ran D, Zhao G, Zeng Y, Chen H, Zou M, Li J, Liang F, Xie M, Hu J, Wang D, Cao H, Provart NJ, Zhang L, Tan X.

Genomics Proteomics Bioinformatics 2020 S1672-0229(18):30216. doi: 10.1016/j.gpb.2019.03.006

PMID: 32224189

Abstract

Tung tree (Vernicia fordii) is an economically important woody oil plant that produces tung oil rich in eleostearic acid. Here, we report a high-quality chromosome-scale genome sequence of tung tree. The genome sequence was assembled by combining Illumina short reads, Pacific Biosciences single-molecule real-time long reads, and Hi-C sequencing data. The size of tung tree genome is 1.12 Gb, with 28,422 predicted genes and over 73% repeat sequences. The V. fordii underwent an ancient genome triplication event shared by core eudicots but no further whole-genome duplication in the subsequent ca. 34.55 million years of evolutionary history of the tung tree lineage. Insertion time analysis revealed that repeat-driven genome expansion might have arisen as a result of long-standing long terminal repeat retrotransposon bursts and lack of efficient DNA deletion mechanisms. The genome harbors 88 resistance genes encoding nucleotide-binding sites; 17 of these genes may be involved in early-infection stage of Fusarium wilt resistance. Further, 651 oil-related genes were identified, 88 of which are predicted to be directly involved in tung oil biosynthesis. Relatively few phosphoenolpyruvate carboxykinase genes, and synergistic effects between transcription factors and oil biosynthesis-related genes might contribute to the high oil content of tung seed. The tung tree genome constitutes a valuable resource for understanding genome evolution, as well as for molecular breeding and genetic improvements for oil production.

An abscisic acid-responsive protein interaction network for sucrose non-fermenting related kinase1 in abiotic stress response

Posted on by biota in Publications

Carina Steliana Carianopol, Aaron Lorheed Chan, Shaowei Dong, Nicholas J. Provart, Shelley Lumba & Sonia Gazzarrini  

Nature Communications Biology 2020 3(1):145. 10.1038/s42003-020-0866-8 PMID:32218501

Abstract

Yeast Snf1 (Sucrose non-fermenting1), mammalian AMPK (5′ AMP-activated protein kinase) and plant SnRK1 (Snf1-Related Kinase1) are conserved heterotrimeric kinase complexes that re-establish energy homeostasis following stress. The hormone abscisic acid (ABA) plays a crucial role in plant stress response. Activation of SnRK1 or ABA signaling results in overlapping transcriptional changes, suggesting these stress pathways share common targets. To investigate how SnRK1 and ABA interact during stress response in Arabidopsis thaliana, we screened the SnRK1 complex by yeast two-hybrid against a library of proteins encoded by 258 ABA-regulated genes. Here, we identify 125 SnRK1- interacting proteins (SnIPs). Network analysis indicates that a subset of SnIPs form signaling modules in response to abiotic stress. Functional studies show the involvement of SnRK1 and select SnIPs in abiotic stress responses. This targeted study uncovers the largest set of SnRK1 interactors, which can be used to further characterize SnRK1 role in plant survival under stress.

QuiC2 represents a functionally distinct class of dehydroshikimate dehydratases identified in Listeria species including Listeria monocytogenes

Posted on by biota in Publications

Kevin Xue, Stephanie M. Prezioso, Dinesh Christendat

Environmental Microbiology 2020 ():. 10.1111/1462-2920.14987 PMID:32190965

Abstract

Many Listeria species including L. monocytogenes contain the pathway for the biosynthesis of protocatechuate from shikimate and quinate. The qui1 and qui2 operons within these Listeria spp. encode enzymes for this pathway. The diversion of shikimate pathway intermediates in some Listeria species to produce protocatechuate suggests an important biological role for this compound to these organisms. A total of seven ORFs, including quiC2, were identified within qui1 and qui2, however only three proteins encoded by the operons have been functionally annotated. The final step in Listeria’s protocatechuate biosynthesis involves the conversion of dehydroshikimate by a dehydroshikimate dehydratase (DSD). In this study, we demonstrate that QuiC2 functions as a DSD in Listeria spp. through biochemical and structural analyses. Moreover, we show that QuiC2 forms a phylogenetic cluster distinct from other functionally annotated DSDs. The individual phylogenetic clusters of DSD are represented by enzymes that produce protocatechuate for distinct biological processes. Similarly, QuiC2 is expected to produce protocatechuate for a novel biological process. We postulate that protocatechuate produced by DSDs found within the QuiC2 phylogenetic cluster provides an ecological niche for representative organisms.