Cluster-specific associations between the gut microbiota and behavioral outcomes in preschool-aged children

BACKGROUND

The gut microbiota is recognized as a regulator of brain development and behavioral outcomes during childhood. Nonetheless, associations between the gut microbiota and behavior are often inconsistent among studies in humans, perhaps because many host-microbe relationships vary widely between individuals. This study aims to stratify children based on their gut microbiota composition (i.e., clusters) and to identify novel gut microbiome cluster-specific associations between the stool metabolomic pathways and child behavioral outcomes.

METHODS

Stool samples were collected from a community sample of 248 typically developing children (3-5 years). The gut microbiota was analyzed using 16S sequencing while LC-MS/MS was used for untargeted metabolomics. Parent-reported behavioral outcomes (i.e., Adaptive Skills, Internalizing, Externalizing, Behavioral Symptoms, Developmental Social Disorders) were assessed using the Behavior Assessment System for Children (BASC-2). Children were grouped based on their gut microbiota composition using the Dirichlet multinomial method, after which differences in the metabolome and behavioral outcomes were investigated.

RESULTS

Four different gut microbiota clusters were identified, where the cluster enriched in both Bacteroides and Bifidobacterium (Ba2) had the most distinct stool metabolome. The cluster characterized by high Bifidobacterium abundance (Bif), as well as cluster Ba2, were associated with lower Adaptive Skill scores and its subcomponent Social Skills. Cluster Ba2 also had significantly lower stool histidine to urocanate turnover, which in turn was associated with lower Social Skill scores in a cluster-dependent manner. Finally, cluster Ba2 had increased levels of compounds involved in Galactose metabolism (i.e., stachyose, raffinose, alpha-D-glucose), where alpha-D-glucose was associated with the Adaptive Skill subcomponent Daily Living scores (i.e., ability to perform basic everyday tasks) in a cluster-dependent manner.

CONCLUSIONS

These data show novel associations between the gut microbiota, its metabolites, and behavioral outcomes in typically developing preschool-aged children. Our results support the concept that cluster-based groupings could be used to develop more personalized interventions to support child behavioral outcomes. Video Abstract.

Exploring Associations Between the Gut Microbiota and Full-Scale Intelligence in Preschool Children

The relationship between the gut microbiota and neurocognitive outcomes is becoming increasingly recognized; however, findings in humans are inconsistent. In addition, few studies have investigated the gut microbial metabolites that may mediate this relationship. The objective of this study was to investigate associations between full-scale intelligence (FSIQ) and the composition of the gut microbiota and metabolome in preschool children. Stool samples were collected from a community sample of 245 typically developing children (3-5 years) from the Alberta Pregnancy Outcomes and Nutrition (APrON) cohort. The faecal microbiome was assessed using 16S rRNA sequencing and the metabolome using LC-MS/MS. FSIQ and scores on the Verbal Comprehension, Visual Spatial, Working Memory indices of the Wechsler Preschool and Primary Scale of Intelligence-IV were used to assess neurocognition. Associations between the gut microbiota and FSIQ were determined using Pearson and Spearman correlations, which were corrected for multiple testing and relevant covariates. Verbal Comprehension negatively correlated with both Shannon alpha diversity (r=-0.14, p=0.032) and the caffeine-derived metabolite paraxanthine (r=-0.22, p<0.001). No other significant correlations were observed. Overall, the weak to modest correlations between Verbal Comprehension with alpha diversity and paraxanthine provide limited evidence of an association between the gut microbiota and neurocognitive outcomes in typically developing preschool children.

Next-generation sequencing approach to investigate genome variability of Parapoxvirus in Canadian muskoxen (Ovibos moschatus)

In 2016, the first orf virus, a double-stranded DNA (dsDNA) virus of the genus parapoxvirus, from a muskox was isolated on Victoria Island, Nunavut (NU), Canada. We used deep sequencing on DNA extracted from orf virus-positive tissues from wild muskoxen from locations on Victoria Island and the adjacent mainland. Orf virus sequence reads derived from four samples were nearly identical. The consensus sequences generated from pooled reads of MxOV comprises of a large contiguous sequence (contig) of 131,759 bp and a smaller right terminal contig of 3552 bp, containing all coding sequences identified as Parapoxvirus. Individual gene comparisons reveal that MxOV shares genetic characteristics with reference strains from both sheep and goat origin. Recombination analysis using Bootscan, MAXCHI, GENECONV, CHIMAERA, SISCAN, and RDP algorithms within the RDP4 software predicted recombination events in two virulence factors, and a large 3000 bp segment of the MxOV genome. Partial B2L nucleotide sequences from strains around the world and other North American isolates were compared to MxOV using MUSCLE alignments and RAxML phylogenetic trees. MxOV was identical to our previously characterized isolate, and shared similarity with orf virus isolated from sheep and goats. The phylogenetic grouping of partial B2L nucleotide sequences did not follow the sample geographic distribution. More full genomes of orf virus, or at least full B2L gene squences, in wildlife are needed especially in North America to better understand the epidemiology of the disease in muskoxen.

Associations Between the Gut Microbiota and Internalizing Behaviors in Preschool Children

OBJECTIVE

Emerging evidence points toward a connection between mental health and the gut microbiota and its metabolites (e.g., short-chain fatty acids). It is unknown whether the gut microbiota is associated with the development of mental health problems (e.g., internalizing or externalizing behaviors) in preschool children. The objective of this study was to evaluate associations between the gut microbiota and internalizing and externalizing behaviors in preschool-aged children.

METHODS

A community sample of 248 typically developing children (3-5 years of age) provided a stool sample for gut microbiota and SCFA analysis. Parents reported child internalizing and externalizing behaviors using the Child Behavior Checklist. Associations between child behaviors and gut microbiota measures were analyzed using Spearman correlations followed by an adjustment for multiple testing, with subanalysis conducted in children clinically "at risk" for behavioral problems compared with those who were not.

RESULTS

There was a correlation between Shannon alpha diversity with internalizing behaviors (rs = -0.134, p = .035) and its subscale somatic complaints (rs = -0.144, p = .023). In addition, children clinically "at risk" for internalizing problems had decreased alpha diversity (U = 551, p = .017). Internalizing behaviors correlated with valerate and isobutyrate (rs = -0.147, p = .021; rs = -0.140, p = .028, respectively). Furthermore the somatic complaints subscale additionally correlated with acetate and butyrate (rs = -0.219, p = .001; rs = -0.241, p < .001, respectively). These findings were also present in children "at risk" for internalizing problems (U = 569, p = .026; U = 571, p = .028) and somatic complaints (U = 164, p = .004; U = 145, p = .001).

CONCLUSIONS

These analyses reveal novel associations between internalizing behaviors and the gut microbiota in preschool children. Furthermore, a relationship between somatic complaints and acetate and butyrate was identified, indicating that interventions that increase SCFA production warrant future investigation.

Characterization of microglial transcriptomes in the brain and spinal cord of mice in early and late experimental autoimmune encephalomyelitis using a RiboTag strategy

Microglia play an important role in the pathogenesis of multiple sclerosis and the mouse model of MS, experimental autoimmune encephalomyelitis (EAE). To more fully understand the role of microglia in EAE we characterized microglial transcriptomes before the onset of motor symptoms (pre-onset) and during symptomatic EAE. We compared the transcriptome in brain, where behavioral changes are initiated, and spinal cord, where damage is revealed as motor and sensory deficits. We used a RiboTag strategy to characterize ribosome-bound mRNA only in microglia without incurring possible transcriptional changes after cell isolation. Brain and spinal cord samples clustered separately at both stages of EAE, indicating regional heterogeneity. Differences in gene expression were observed in the brain and spinal cord of pre-onset and symptomatic animals with most profound effects in the spinal cord of symptomatic animals. Canonical pathway analysis revealed changes in neuroinflammatory pathways, immune functions and enhanced cell division in both pre-onset and symptomatic brain and spinal cord. We also observed a continuum of many pathways at pre-onset stage that continue into the symptomatic stage of EAE. Our results provide additional evidence of regional and temporal heterogeneity in microglial gene expression patterns that may help in understanding mechanisms underlying various symptomology in MS.

A repeatable and quantitative DNA metabarcoding assay to characterize mixed strongyle infections in horses

Horses are ubiquitously infected by a diversity of gastro-intestinal parasitic helminths. Of particular importance are nematodes of the family Strongylidae, which can significantly impact horse health and performance. However, knowledge about equine strongyles remains limited due to our inability to identify most species non-invasively using traditional morphological techniques. We developed a new internal transcribed spacer 2 (ITS2) DNA metabarcoding 'nemabiome' assay to characterise mixed strongyle infections in horses and assessed its performance by applying it to pools of infective larvae from fecal samples from an experimental herd in Kentucky, USA and two feral horse populations from Sable Island and Alberta, Canada. In addition to reporting the detection of 33 different species with high confidence, we illustrate the assay's repeatability by comparing results generated from aliquots from the same fecal samples and from individual horses sampled repeatedly over multiple days or months. We also validate the quantitative potential of the assay by demonstrating that the proportion of amplicon reads assigned to different species scales linearly with the number of larvae present. This new tool significantly improves equine strongyle diagnostics, presenting opportunities for research on species-specific anthelmintic resistance and the causes and consequences of variation in mixed infections.

Crohn's Disease Pathobiont Adherent-Invasive E coli Disrupts Epithelial Mitochondrial Networks With Implications for Gut Permeability

BACKGROUND & AIMS

Adherent-invasive Escherichia coli are implicated in inflammatory bowel disease, and mitochondrial dysfunction has been observed in biopsy specimens from patients with inflammatory bowel disease. As a novel aspect of adherent-invasive E coli-epithelial interaction, we hypothesized that E coli (strain LF82) would elicit substantial disruption of epithelial mitochondrial form and function.

METHODS

Monolayers of human colon-derived epithelial cell lines were exposed to E coli-LF82 or commensal E coli and RNA sequence analysis, mitochondrial function (adenosine triphosphate synthesis) and dynamics (mitochondrial network imaging, immunoblotting for fission and fusion proteins), and epithelial permeability (transepithelial resistance, flux of fluorescein isothiocyanate-dextran and bacteria) were assessed.

RESULTS

E coli-LF82 significantly affected epithelial expression of ∼8600 genes, many relating to mitochondrial function. E coli-LF82-infected epithelia showed swollen mitochondria, reduced mitochondrial membrane potential and adenosine triphosphate, and fragmentation of the mitochondrial network: events not observed with dead E coli-LF82, medium from bacterial cultures, or control E coli. Treatment with Mitochondrial Division Inhibitor 1 (Mdivi1, inhibits dynamin-related peptide 1, guanosine triphosphatase principally responsible for mitochondrial fission) or P110 (prevents dynamin-related peptide 1 binding to mitochondrial fission 1 protein) partially reduced E coli-LF82-induced mitochondrial fragmentation in the short term. E coli-LF82-infected epithelia showed loss of the long isoform of optic atrophy factor 1, which mediates mitochondrial fusion. Mitochondrial Division Inhibitor 1 reduced the magnitude of E coli-LF82-induced increased transepithelial flux of fluorescein isothiocyanate dextran. By 8 hours after infection, increased cytosolic cytochrome C and DNA fragmentation were apparent without evidence of caspase-3 or apoptosis inducing factor activation.

CONCLUSIONS

Epithelial mitochondrial fragmentation caused by E coli-LF82 could be targeted to maintain cellular homeostasis and mitigate infection-induced loss of epithelial barrier function. Data have been deposited in NCBI's Gene Expression Omnibus and are accessible through GEO series accession numbers GSE154121 and GSE154122 (https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE154121).

Human interleukin-4-treated regulatory macrophages promote epithelial wound healing and reduce colitis in a mouse model

Murine alternatively activated macrophages can exert anti-inflammatory effects. We sought to determine if IL-4-treated human macrophages [i.e., hM(IL4)] would promote epithelial wound repair and can serve as a cell transfer treatment for inflammatory bowel disease (IBD). Blood monocytes from healthy volunteers and patients with active and inactive IBD were converted to hM(IL4)s. IL-4 treatment of blood-derived macrophages from healthy volunteers and patients with inactive IBD resulted in a characteristic CD206⁺CCL18⁺CD14^low/- phenotype (RNA-seq revealed IL-4 affected expression of 996 genes). Conditioned media from freshly generated or cryopreserved hM(IL4)s promoted epithelial wound healing in part by TGF, and reduced cytokine-driven loss of epithelial barrier function in vitro. Systemic delivery of hM(IL4) to dinitrobenzene sulphonic acid (DNBS)-treated Rag1^-/- mice significantly reduced disease. These findings from in vitro and in vivo analyses provide proof-of-concept support for the development of autologous M(IL4) transfer as a cellular immunotherapy for IBD.

Cystic Fibrosis Patients Infected With Epidemic Pseudomonas aeruginosa Strains Have Unique Microbial Communities

Pseudomonas aeruginosa is the archetypal cystic fibrosis (CF) pathogen. However, the clinical course experienced by infected individuals varies markedly. Understanding these differences is imperative if further improvements in outcomes are to be achieved. Multiple studies have found that patients infected with epidemic P. aeruginosa (ePA) strains may have a worse clinical prognosis than those infected with unique, non-clonal strains. Additionally, the traditionally uncultured CF lung bacterial community (i.e., CF microbiome) may further influence the outcome. We sought to identify if these two important variables, not identified through routine culture, associate and together may contribute to disease pathogenesis. Patients were classified as being infected with Prairie Epidemic ePA (PES) or a non-clonal strain, unique PA strains (uPA), through a retrospective assessment of a comprehensive strain biobank using a combination of PFGE and PES-specific PCR. Patients were matched to age, sex, time-period controls and sputum samples from equivalent time periods were identified from a sputum biobank. Bacterial 16S rRNA gene profiling and Pseudomonas qPCR was used to characterize the respiratory microbiome. We identified 31 patients infected with PES and matched them with uPA controls. Patients infected with PES at baseline have lower microbial diversity (P = 0.02) and higher P. aeruginosa relative abundance (P < 0.005). Microbial community structure was found to cluster by PA strain type, although it was not the main determinant of community structure as additional factors were also found to be drivers of CF community structure. Communities from PES infected individuals were enriched with Pseudomonas, Streptococcus and Prevotella OTUs. The disproportionate disease experienced by ePA infected CF patients may be mediated through a combination of pathogen-pathogen factors as opposed to strictly enhanced virulence of infecting P. aeruginosa strains.

Th1-Th2 Cross-Regulation Controls Early Leishmania Infection in the Skin by Modulating the Size of the Permissive Monocytic Host Cell Reservoir

The impact of T helper (Th) 1 versus Th2 immunity on intracellular infections is attributed to classical versus alternative activation of macrophages leading to resistance or susceptibility. However, observations in multiple infectious settings demonstrate deficiencies in mediators of Th1-Th2 immunity, which have paradoxical or no impact. We report that prior to influencing activation, Th1/Th2 immunity first controls the size of the permissive host cell reservoir. During early Leishmania infection of the skin, IFN-γ- or STAT6-mediated changes in phagocyte activation were counteracted by changes in IFN-γ-mediated recruitment of permissive CCR2⁺ monocytes. Monocytes were required for early parasite expansion and acquired an alternatively activated phenotype despite the Th1 dermal environment required for their recruitment. Surprisingly, STAT6 did not enhance intracellular parasite proliferation, but rather modulated the size and permissiveness of the monocytic host cell reservoir via regulation of IFN-γ and IL-10. These observations expand our understanding of the Th1-Th2 paradigm during infection.

Benchmarking hybrid assemblies of Giardia and prediction of widespread intra-isolate structural variation

Background

Currently available short read genome assemblies of the tetraploid protozoan parasite Giardia intestinalis are highly fragmented, highlighting the need for improved genome assemblies at a reasonable cost. Long nanopore reads are well suited to resolve repetitive genomic regions resulting in better quality assemblies of eukaryotic genomes. Subsequent addition of highly accurate short reads to long-read assemblies further improves assembly quality. Using this hybrid approach, we assembled genomes for three Giardia isolates, two with published assemblies and one novel, to evaluate the improvement in genome quality gained from long reads. We then used the long reads to predict structural variants to examine this previously unexplored source of genetic variation in Giardia.

Methods

With MinION reads for each isolate, we assembled genomes using several assemblers specializing in long reads. Assembly metrics, gene finding, and whole genome alignments to the reference genomes enabled direct comparison to evaluate the performance of the nanopore reads. Further improvements from adding Illumina reads to the long-read assemblies were evaluated using gene finding. Structural variants were predicted from alignments of the long reads to the best hybrid genome for each isolate and enrichment of key genes was analyzed using random genome sampling and calculation of percentiles to find thresholds of significance.

Results

Our hybrid assembly method generated reference quality genomes for each isolate. Consistent with previous findings based on SNPs, examination of heterozygosity using the structural variants found that Giardia BGS was considerably more heterozygous than the other isolates that are from Assemblage A. Further, each isolate was shown to contain structural variant regions enriched for variant-specific surface proteins, a key class of virulence factor in Giardia.

Conclusions

The ability to generate reference quality genomes from a single MinION run and a multiplexed MiSeq run enables future large-scale comparative genomic studies within the genus Giardia. Further, prediction of structural variants from long reads allows for more in-depth analyses of major sources of genetic variation within and between Giardia isolates that could have effects on both pathogenicity and host range.

Complex patterns of cell growth in the placenta in normal pregnancy and as adaptations to maternal diet restriction

The major milestones in mouse placental development are well described, but our understanding is limited to how the placenta can adapt to damage or changes in the environment. By using stereology and expression of cell cycle markers, we found that the placenta grows under normal conditions not just by hyperplasia of trophoblast cells but also through extensive polyploidy and cell hypertrophy. In response to feeding a low protein diet to mothers prior to and during pregnancy, to mimic chronic malnutrition, we found that this normal program was altered and that it was influenced by the sex of the conceptus. Male fetuses showed intrauterine growth restriction (IUGR) by embryonic day (E) 18.5, just before term, whereas female fetuses showed IUGR as early as E16.5. This difference was correlated with differences in the size of the labyrinth layer of the placenta, the site of nutrient and gas exchange. Functional changes were implied based on up-regulation of nutrient transporter genes. The junctional zone was also affected, with a reduction in both glycogen trophoblast and spongiotrophoblast cells. These changes were associated with increased expression of Phlda2 and reduced expression of Egfr. Polyploidy, which results from endoreduplication, is a normal feature of trophoblast giant cells (TGC) but also spongiotrophoblast cells. Ploidy was increased in sinusoidal-TGCs and spongiotrophoblast cells, but not parietal-TGCs, in low protein placentas. These results indicate that the placenta undergoes a range of changes in development and function in response to poor maternal diet, many of which we interpret are aimed at mitigating the impacts on fetal and maternal health.

foxc1 is required for embryonic head vascular smooth muscle differentiation in zebrafish

Vascular smooth muscle of the head derives from neural crest, but developmental mechanisms and early transcriptional drivers of the vSMC lineage are not well characterized. We find that in early development, the transcription factor foxc1b is expressed in mesenchymal cells that associate with the vascular endothelium. Using timelapse imaging, we observe that foxc1b expressing mesenchymal cells differentiate into acta2 expressing vascular mural cells. We show that in zebrafish, while foxc1b is co-expressed in acta2 positive smooth muscle cells that associate with large diameter vessels, it is not co-expressed in capillaries where pdgfrβ positive pericytes are located. In addition to being an early marker of the lineage, foxc1 is essential for vSMC differentiation; we find that foxc1 loss of function mutants have defective vSMC differentiation and that early genetic ablation of foxc1b or acta2 expressing populations blocks vSMC differentiation. Furthermore, foxc1 is expressed upstream of acta2 and is required for acta2 expression in vSMCs. Using RNA-Seq we determine an enriched intersectional gene expression profile using dual expression of foxc1b and acta2 to identify novel vSMC markers. Taken together, our data suggests that foxc1 is a marker of vSMCs and plays a critical functional role in promoting their differentiation.

The effects of cycled inhaled aztreonam on the cystic fibrosis (CF) lung microbiome

BACKGROUND

To improve clinical outcomes, cystic fibrosis (CF) patients with chronic Pseudomonas aeruginosa infections are prescribed inhaled anti-pseudomonal antibiotics. Although, a diverse microbial community exists within CF airways, little is known about how the CF microbiota influences patient outcomes. We hypothesized that organisms within the CF microbiota are affected by inhaled-antibiotics and baseline microbiome may be used to predict therapeutic response.

METHODS

Adults with chronic P. aeruginosa infection from four clinics were observed during a single 28-day on/off inhaled-aztreonam cycle. Patients performed serial sputum collection, CF-respiratory infection symptom scores (CRISS), and spirometry. Patients achieving a decrease of ≥2 CRISS by day 28 were categorized as subjective responders (SR). The airway microbiome was defined by Illumina MiSeq analysis of the 16S rRNA gene.

RESULTS

Thirty-seven patients (median 37.4 years and FEV₁ 44% predicted) were enrolled. No significant cohort-wide changes in the microbiome were observed between on/off AZLI cycles in either alpha- or beta-diversity metrics. However, at an individual level shifts were apparent. Twenty-one patients (57%) were SR and fourteen patients did not subjectively respond. While alpha-diversity metrics did not associate with response, patients who did not subjectively respond had a higher abundance of Staphylococcus and Streptococcus, and lower abundance of Haemophilus.

CONCLUSIONS

The CF microbiome is relatively resilient to AZLI perturbations. However, associated changes were observed at the individual patient level. The relative abundance of key "off-target" organisms associated with subjective improvements suggesting that the microbiome may be used as a tool to predict patient response - potentially improving outcomes.

Host genetics and diet composition interact to modulate gut microbiota and predisposition to metabolic syndrome in spontaneously hypertensive stroke-prone rats

Metabolic syndrome encompasses obesity, glucose intolerance, hypertension, and dyslipidemia; however, the interactions between diet and host physiology that predispose to metabolic syndrome are incompletely understood. Here, we explored the effects of a high-fat diet (HFD) on energy balance, gut microbiota, and risk factors of metabolic syndrome in spontaneously hypertensive stroke-prone (SHRSP) and Wistar-Kyoto (WKY) rats. We found that the SHRSP rats were hypertensive, hyperphagic, less sensitive to hypophagic effects of exendin-4, and expended more energy with diminished sensitivity to sympathetic blockade compared to WKY rats. Notably, key thermogenic markers in brown and retroperitoneal adipose tissues and skeletal muscle were up-regulated in SHRSP than WKY rats. Although HFD promoted weight gain, adiposity, glucose intolerance, hypertriglyceridemia, hepatic lipidosis, and hyperleptinemia in both SHRSP and WKY rats, the SHRSP rats weighed less but had comparable percent adiposity to WKY rats, which supports the use of HFD-fed SHRSP rats as a unique model for studying the metabolically obese normal weight (MONW) phenotype in humans. Despite distinct strain differences in gut microbiota composition, diet had a preponderant impact on gut flora with some of the taxa being strongly associated with key metabolic parameters. Together, we provide evidence that interactions between host genetics and diet modulate gut microbiota and predispose SHRSP rats to develop metabolic syndrome.-Singh, A., Zapata, R. C., Pezeshki, A., Workentine, M. L., Chelikani, P. K. Host genetics and diet composition interact to modulate gut microbiota and predisposition to metabolic syndrome in spontaneously hypertensive stroke-prone rats.

Macrophages treated with antigen from the tapeworm Hymenolepis diminuta condition CD25+ T cells to suppress colitis

Macrophages play central roles in immunity as early effectors and modulating adaptive immune reponses; we implicated macrophages in the anticolitic effect of infection with the tapeworm Hymenolepis diminuta. Here, gene arrays revealed that H. diminuta antigen (HdAg) evoked a program in murine macrophages distinct from that elicited by IL-4. Further, HdAg suppressed LPS-evoked release of TNF-α and IL-1β from macrophages via autocrine IL-10 signaling. In assessing the ability of macrophages treated in vitro with an extract of H. diminuta [M^(HdAg)] to affect disease, intravenous, but not peritoneal, injection of M^(HdAg) protected wild-type but not RAG1^-/- mice from dinitrobenzene sulphonic acid (DNBS)-induced colitis. Administration of splenic CD4⁺ T cells from in vitro cocultures with M^(HdAg), but not those cocultured with M^(IL-4) cells, inhibited DNBS-induced colitis; fractionation of the T-cell population indicated that the CD4⁺CD25⁺ T cells from cocultures with M^(HdAg) drove the suppression of DNBS-induced colitis. Use of IL-4^-/- or IL-10^-/- CD4⁺ T cells revealed that neither cytokine alone from the donor cells was essential for the anticolitic effect. These data illustrate that HdAg evokes a unique regulatory program in macrophages, identifies HdAg-evoked IL-10 suppression of macrophage activation, and reveals the ability of HdAg-treated macrophages to educate ( i.e., condition) and mobilize CD4⁺CD25⁺ T cells, which could be deployed to treat colonic inflammation.-Reyes, J. L., Lopes, F., Leung, G., Jayme, T. S., Matisz, C. E., Shute, A., Burkhard, R., Carneiro, M., Workentine, M. L., Wang, A., Petri, B., Beck, P. L., Geuking, M. B., McKay, D. M., Macrophages treated with antigen from the tapeworm Hymenolepis diminuta condition CD25⁺ T cells to suppress colitis.

Sputum microbiota is predictive of long-term clinical outcomes in young adults with cystic fibrosis

BACKGROUND

Complex polymicrobial communities infect cystic fibrosis (CF) lower airways. Generally, communities with low diversity, dominated by classical CF pathogens, associate with worsened patient status at sample collection. However, it is not known if the microbiome can predict future outcomes. We sought to determine if the microbiome could be adapted as a biomarker for patient prognostication.

METHODS

We retrospectively assessed prospectively collected sputum from a cohort of 104 individuals aged 18-22 to determine factors associated with progression to early end-stage lung disease (eESLD; death/transplantation <25 years) and rapid pulmonary function decline (>-3%/year FEV₁ over the ensuing 5 years). Illumina MiSeq paired-end sequencing of the V3-V4 region of the 16S rRNA was used to define the airway microbiome.

RESULTS

Based on the primary outcome analysed, 17 individuals (16%) subsequently progressed to eESLD. They were more likely to have sputum with low alpha diversity, dominated by specific pathogens including Pseudomonas. Communities with abundant Streptococcus were observed to be protective. Microbial communities clustered together by baseline lung disease stage and subsequent progression to eESLD. Multivariable analysis identified baseline lung function and alpha diversity as independent predictors of eESLD. For the secondary outcomes, 58 and 47 patients were classified as rapid progressors based on absolute and relative definitions of lung function decline, respectively. Patients with low alpha diversity were similarly more likely to be classified as experiencing rapid lung function decline over the ensuing 5 years when adjusted for baseline lung function.

CONCLUSIONS

We observed that the diversity of microbial communities in CF airways is predictive of progression to eESLD and disproportionate lung function decline and may therefore represent a novel biomarker.

Effects of transportation to and co-mingling at an auction market on nasopharyngeal and tracheal bacterial communities of recently weaned beef cattle

The objective was to study effects of transportation to and co-mingling at an auction market on nasopharyngeal and tracheal bacterial communities of feedlot cattle. Two groups of 30 Angus-cross heifers were studied from weaning to 28 d after arrival at a feedlot. For each group, half the heifers were either transported directly to a feedlot after weaning (RANC) or transported to and co-mingled at an auction market for 24 h before being placed in a feedlot (AUCT). Deep nasal swabs (DNS) and trans-tracheal aspirates (TTA) were collected at weaning (d0) and at on-arrival processing at the feedlot (d2). At 7 (d9) and 28 d (d30) after arrival, DNS were repeated. The DNA was extracted from DNS and TTA and the V4 region of the 16S rRNA gene sequenced (MiSeq). Alpha diversity analysis did not reveal differences between AUCT and RANC. However, bacterial diversity decreased over time in the nasopharynx, especially at d9. Although beta-diversity was not different between AUCT and RANC, interval after arrival and feedlot where heifers were placed affected composition of the nasopharyngeal bacterial communities. In both groups, a large increase in Mycoplasma was observed after arrival; in one group, Mycoplasma bovis was dominant at d9 and remained dominant until d30. However, in the other group, Mycoplasma dispar dominated at d9 and was replaced by Moraxella at d30. We concluded that transportation to and co-mingling at an auction market for 24 h did not significantly influence diversity or composition of nasopharyngeal or tracheal bacterial communities.

An intact microbiota is required for the gastrointestinal toxicity of the immunosuppressant mycophenolate mofetil

BACKGROUND

Mycophenolate mofetil (MMF) is commonly prescribed after transplantation and has major advantages over other immunosuppressive drugs, but frequent gastrointestinal (GI) side-effects limit its use. The mechanism(s) underlying MMF-related GI toxicity have yet to be elucidated.

METHODS

To investigate MMF-related GI toxicity, experimental mice were fed chow containing MMF (0.563%) and multiple indices of toxicity, including weight loss and colonic inflammation, were measured. Changes in intestinal microbial composition were detected using 16S rRNA Illumina sequencing, and downstream PICRUSt analysis was used to predict metagenomic pathways involved. Germ-free (GF) mice and mice treated with orally administered broad-spectrum antibiotics (ABX) were utilized to interrogate the importance of the microbiota in MMF-induced GI toxicity.

RESULTS

Mice treated with MMF exhibited significant weight loss, related to loss of body fat and muscle, and marked colonic inflammation. MMF exposure was associated with changes in gut microbial composition, as demonstrated by a loss of overall diversity, expansion of Proteobacteria (specifically Escherichia/Shigella), and enrichment of genes involved in lipopolysaccharide (LPS) biosynthesis, which paralleled increased levels of LPS in the feces and serum. MMF-related GI toxicity was dependent on the intestinal microbiota, as MMF did not induce weight loss or colonic inflammation in GF mice. Furthermore, ABX prevented and reversed MMF-induced weight loss and colonic inflammation.

CONCLUSIONS

An intact intestinal microbiota is required to initiate and sustain the GI toxicity of MMF. MMF treatment causes dynamic changes in the composition of the intestinal microbiota that may be a targetable driver of the GI side-effects of MMF.

The effects of inhaled aztreonam on the cystic fibrosis lung microbiome

BACKGROUND

Aztreonam lysine for inhalation (AZLI) is an inhaled antibiotic used to treat chronic Pseudomonas aeruginosa infection in CF. AZLI improves lung function and quality of life, and reduces exacerbations-improvements attributed to its antipseudomonal activity. Given the extremely high aztreonam concentrations achieved in the lower airways by nebulization, we speculate this may extend its spectrum of activity to other organisms. As such, we sought to determine if AZLI affects the CF lung microbiome and whether community constituents can be used to predict treatment responsiveness.

METHODS

Patients were included if they had chronic P. aeruginosa infection and repeated sputum samples collected before and after AZLI. Sputum DNA was extracted, and the V3-hypervariable region of the 16S ribosomal RNA (rRNA) gene amplified and sequenced.

RESULTS

Twenty-four patients naïve to AZLI contributed 162 samples. The cohort had a median age of 37.1 years, and a  median FEV₁ of 44% predicted. Fourteen patients were a priori defined as responders for achieving ≥3% FEV₁ improvement following initiation. No significant changes in alpha diversity were noted following AZLI. Furthermore, beta diversity demonstrated clustering with respect to patients, but had no association with AZLI use. However, we did observe a decline in the relative abundance of several individual operational taxonomic units (OTUs) following AZLI initiation suggesting that specific sub-populations of organisms may be impacted. Patients with higher abundance of Staphylococcus and anaerobic organisms including Prevotella and Fusobacterium were less likely to respond to therapy.

CONCLUSIONS

Results from our study suggest potential alternate/additional mechanisms by which AZLI functions. Moreover, our study suggests that the CF microbiota may be used as a biomarker to predict patient responsiveness to therapy suggesting the microbiome may be harnessed for the personalization of therapies.

Changes in Composition of the Gut Bacterial Microbiome after Fecal Microbiota Transplantation for Recurrent Clostridium difficile Infection in a Pediatric Heart Transplant Patient

The microbiome is increasingly recognized as an important influence on human health and many of the comorbidities that affect patients after solid organ transplantation (SOT) have been shown to involve changes in gut bacterial populations. Thus, microbiome changes in an individual patient may have important health implications after SOT but this area remains understudied. We describe changes in the composition of the fecal microbiome from a pediatric heart transplant recipient before and >2.5 years after he underwent repeated fecal microbiota transplantation (FMT) for recurrent Clostridium difficile infection (CDI). With both documented episodes of CDI, there was marked loss of bacterial diversity with overgrowth of Proteobacteria (>98.9% of phyla identified) associated with symptomatic colitis that was corrected after FMT. We hypothesize that a second CDI occurring after FMT was related to incomplete restoration of normal bowel flora post-FMT with relative deficiencies of the phyla Firmicutes and Bacteroidetes and the families Lachnospiraceae and Ruminococcaceae. Following the second FMT, there was a gradual shift in gut bacterial composition coincident with the recipient developing lymphonodular hyperplasia of the colon and painless hematochezia that resolved with discontinuation of mycophenolate mofetil (MMF). This case documents dynamic changes in the bacterial microbiome after FMT and suggests that MMF may influence the gut microbiome with consequences for the patient.

Giardia duodenalis induces pathogenic dysbiosis of human intestinal microbiota biofilms

Giardia duodenalis is a prevalent cause of acute diarrheal disease worldwide. However, recent outbreaks in Italy and Norway have revealed a link between giardiasis and the subsequent development of chronic post-infectious irritable bowel syndrome. While the mechanisms underlying the causation of post-infectious irritable bowel syndrome remain obscure, recent findings suggest that alterations in gut microbiota communities are linked to the pathophysiology of irritable bowel syndrome. In the present study, we use a laboratory biofilm system to culture and enrich mucosal microbiota from human intestinal biopsies. Subsequently, we show that co-culture with Giardia induces disturbances in biofilm species composition and biofilm structure resulting in microbiota communities that are intrinsically dysbiotic - even after the clearance of Giardia. These microbiota abnormalities were mediated in part by secretory-excretory Giardia cysteine proteases. Using in vitro cell culture and germ-free murine infection models, we show that Giardia-induced disruptions of microbiota promote bacterial invasion, resulting in epithelial apoptosis, tight junctional disruption, and bacterial translocation across an intestinal epithelial barrier. Additionally, these dysbiotic microbiota communities resulted in increased activation of the Toll-like receptor 4 signalling pathway, and overproduction of the pro-inflammatory cytokine IL-1beta in humanized germ-free mice. Previous studies that have sought explanations and risk factors for the development of post-infectious irritable bowel syndrome have focused on features of enteropathogens and attributes of the infected host. We propose that polymicrobial interactions involving Giardia and gut microbiota may cause persistent dysbiosis, offering a new interpretation of the reasons why those afflicted with giardiasis are predisposed to gastrointestinal disorders post-infection.

Genomic analysis of the multi-host pathogen Erysipelothrix rhusiopathiae reveals extensive recombination as well as the existence of three generalist clades with wide geographic distribution

Background

Knowledge about how bacterial populations are structured is an important prerequisite for studying their ecology and evolutionary history and facilitates inquiry into host specificity, pathogenicity, geographic dispersal and molecular epidemiology. Erysipelothrix rhusiopathiae is an opportunistic pathogen that is currently reemerging in both the swine and poultry industries globally. This bacterium sporadically causes mortalities in captive marine mammals, and has recently been implicated in large-scale wildlife die-offs. However, despite its economic relevance and broad geographic and host distribution, including zoonotic potential, the global diversity, recombination rates, and population structure of this bacterium remain poorly characterized. In this study, we conducted a broad-scale genomic comparison of E. rhusiopathiae based on a diverse collection of isolates in order to address these knowledge gaps.

Results

Eighty-three E. rhusiopathiae isolates from a range of host species and geographic origins, isolated between 1958 and 2014, were sequenced and assembled using both reference-based mapping and de novo assembly. We found that a high proportion of the core genome (58 %) had undergone recombination. Therefore, we used three independent methods robust to the presence of recombination to define the population structure of this species: a phylogenetic tree based on a set of conserved protein sequences, in silico chromosome painting, and network analysis. All three methods were broadly concordant and supported the existence of three distinct clades within the species E. rhusiopathiae. Although we found some evidence of host and geographical clustering, each clade included isolates from diverse host species and from multiple continents.

Conclusions

Using whole genome sequence data, we confirm recent suggestions that E. rhusiopathiae is a weakly clonal species that has been shaped extensively by homologous recombination. Despite frequent recombination, we can reliably identify three distinct clades that do not clearly segregate by host species or geographic origin. Our results provide an essential baseline for future molecular epidemiological, ecological and evolutionary studies of E. rhusiopathiae and facilitate comparisons to other recombinogenic, multi-host bacteria.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2643-0) contains supplementary material, which is available to authorized users.

Cyanide Toxicity to Burkholderia cenocepacia Is Modulated by Polymicrobial Communities and Environmental Factors

Microbes within polymicrobial communities can establish positive and negative interactions that have the potential to influence the overall behavior of the community. Pseudomonas aeruginosa and species of the Burkholderia cepacia complex (Bcc) can co-exist in the lower airways, however several studies have shown that P. aeruginosa can effectively kill the Bcc in vitro, for which hydrogen cyanide (HCN) was recently proposed to play a critical role. Here we show that modification of the environment (i.e., culture medium), long-term genetic adaptation of P. aeruginosa to the cystic fibrosis (CF) lung, or the addition of another bacterial species to the community can alter the sensitivity of Burkholderia cenocepacia to P. aeruginosa toxins. We specifically demonstrate that undefined rich media leads to higher susceptibility of B. cenocepacia to P. aeruginosa toxins like cyanide as compared to a synthetic medium (SCFM), that mimics the CF lung nutritional content. Overall, our study shows that the polymicrobial environment can have profound effects on negative interactions mediated by P. aeruginosa against B. cenocepacia. In fact, evolved P. aeruginosa or the presence of other species such as Staphylococcus aureus can directly abolish the direct competition mediated by cyanide and consequently maintaining a higher level of species diversity within the community.

The Challenge and Potential of Metagenomics in the Clinic

The bacteria, fungi, and viruses that live on and in us have a tremendous impact on our day-to-day health and are often linked to many diseases, including autoimmune disorders and infections. Diagnosing and treating these disorders relies on accurate identification and characterization of the microbial community. Current sequencing technologies allow the sequencing of the entire nucleic acid complement of a sample providing an accurate snapshot of the community members present in addition to the full genetic potential of that microbial community. There are a number of clinical applications that stand to benefit from these data sets, such as the rapid identification of pathogens present in a sample. Other applications include the identification of antibiotic-resistance genes, diagnosis and treatment of gastrointestinal disorders, and many other diseases associated with bacterial, viral, and fungal microbiomes. Metagenomics also allows the physician to probe more complex phenotypes such as microbial dysbiosis with intestinal disorders and disruptions of the skin microbiome that may be associated with skin disorders. Many of these disorders are not associated with a single pathogen but emerge as a result of complex ecological interactions within microbiota. Currently, we understand very little about these complex phenotypes, yet clearly they are important and in some cases, as with fecal microbiota transplants in Clostridium difficile infections, treating the microbiome of the patient is effective. Here, we give an overview of metagenomics and discuss a number of areas where metagenomics is applicable in the clinic, and progress being made in these areas. This includes (1) the identification of unknown pathogens, and those pathogens particularly hard to culture, (2) utilizing functional information and gene content to understand complex infections such as Clostridium difficile, and (3) predicting antimicrobial resistance of the community using genetic determinants of resistance identified from the sequencing data. All of these applications rely on sophisticated computational tools, and we also discuss the importance of skilled bioinformatic support for the implementation and use of metagenomics in the clinic.

Spatial distributions of Pseudomonas fluorescens colony variants in mixed-culture biofilms

Background

The emergence of colony morphology variants in structured environments is being recognized as important to both niche specialization and stress tolerance. Pseudomonas fluorescens demonstrates diversity in both its natural environment, the rhizosphere, and in laboratory grown biofilms. Sub-populations of these variants within a biofilm have been suggested as important contributors to antimicrobial stress tolerance given their altered susceptibility to various agents. As such it is of interest to determine how these variants might be distributed in the biofilm environment.

Results

Here we present an analysis of the spatial distribution of Pseudomonas fluorescens colony morphology variants in mixed-culture biofilms with the wildtype phenotype. These findings reveal that two variant colony morphotypes demonstrate a significant growth advantage over the wildtype morphotype in the biofilm environment. The two variant morphotypes out-grew the wildtype across the entire biofilm and this occurred within 24 h and was maintained through to 96 h. This competitive advantage was not observed in homogeneous broth culture.

Conclusions

The significant advantage that the variants demonstrate in biofilm colonization over the wildtype denotes the importance of this phenotype in structured environments.

Phenotypic heterogeneity of Pseudomonas aeruginosa populations in a cystic fibrosis patient

The opportunistic pathogen Pseudomonas aeruginosa chronically infects the lower airways of patients with cystic fibrosis. Throughout the course of infection this organism undergoes adaptations that contribute to its long-term persistence in the airways. While P. aeruginosa diversity has been documented, it is less clear to what extent within-patient diversity contributes to the overall population structure as most studies have been limited to the analysis of only a few isolates per patient per time point. To examine P. aeruginosa population structure in more detail we collected multiple isolates from individual sputum samples of a patient chronically colonized with P. aeruginosa. This strain collection, comprised of 169 clonal isolates and representing three pulmonary exacerbations as well as clinically stable periods, was assayed for a wide selection of phenotypes. These phenotypes included colony morphology, motility, quorum sensing, protease activity, auxotrophy, siderophore levels, antibiotic resistance, and growth profiles. Each phenotype displayed significant variation even within isolates of the same colony morphotype from the same sample. Isolates demonstrated a large degree of individuality across phenotypes, despite being part of a single clonal lineage, suggesting that the P. aeruginosa population in the cystic fibrosis airways is being significantly under-sampled.

Metabolomics and its application to studying metal toxicity

Here we explain the omics approach of metabolomics and how it can be applied to study a physiological response to toxic metal exposure. This review aims to educate the metallomics field to the tool of metabolomics. Metabolomics is becoming an increasingly used tool to compare natural and challenged states of various organisms, from disease states in humans to toxin exposure to environmental systems. This approach is key to understanding and identifying the cellular or biochemical targets of metals and the underlying physiological response. Metabolomics steps are described and overviews of its application to metal toxicity to organisms are given. As this approach is very new there are yet only a small number of total studies and therefore only a brief overview of some metal metabolomics studies is described. A frank critical evaluation of the approach is given to provide newcomers to the method a clear idea of the challenges and the rewards of applying metabolomics to their research.

Differences in metabolism between the biofilm and planktonic response to metal stress

Bacterial biofilms are known to withstand the effects of toxic metals better than planktonic cultures of the same species. This phenomenon has been attributed to many features of the sessile lifestyle not present in free-swimming populations, but the contribution of intracellular metabolism has not been previously examined. Here, we use a combined GC-MS and (1)H NMR metabolomic approach to quantify whole-cell metabolism in biofilm and planktonic cultures of the multimetal resistant bacterium Pseudomonas fluorescens exposed to copper ions. Metabolic changes in response to metal exposure were found to be significantly different in biofilms compared to planktonic cultures. Planktonic metabolism indicated an oxidative stress response that was characterized by changes to the TCA cycle, glycolysis, pyruvate and nicotinate and niacotinamide metabolism. Similar metabolic changes were not observed in biofilms, which were instead dominated by shifts in exopolysaccharide related metabolism suggesting that metal stress in biofilms induces a protective response rather than the reactive changes observed for the planktonic cells. From these results, we conclude that differential metabolic shifts play a role in biofilm-specific multimetal resistance and tolerance. An altered metabolic response to metal toxicity represents a novel addition to a growing list of biofilm-specific mechanisms to resist environmental stress.

Metabolomic investigation of the bacterial response to a metal challenge

Pseudomonas pseudoalcaligenes KF707 is naturally resistant to the toxic metalloid tellurite, but the mechanisms of resistance are not known. In this study we report the isolation of a KF707 mutant (T5) with hyperresistance to tellurite. In order to characterize the bacterial response and the pathways leading to tolerance, we utilized Phenotype MicroArray technology (Biolog) and a metabolomic technique based on nuclear magnetic resonance spectroscopy. The physiological states of KF707 wild-type and T5 cells exposed to tellurite were also compared in terms of viability and reduced thiol content. Our analyses showed an extensive change in metabolism upon the addition of tellurite to KF707 cultures as well as different responses when the wild-type and T5 strains were compared. Even in the absence of tellurite, T5 cells displayed a "poised" physiological status, primed for tellurite exposure and characterized by altered intracellular levels of glutathione, branched-chain amino acids, and betaine, along with increased resistance to other toxic metals and metabolic inhibitors. We conclude that hyperresistance to tellurite in P. pseudoalcaligenes KF707 is correlated with the induction of the oxidative stress response, resistance to membrane perturbation, and reconfiguration of cellular metabolism.

The GacS-GacA two-component regulatory system of Pseudomonas fluorescens: a bacterial two-hybrid analysis

The two-component regulatory system comprised of the sensor kinase, GacS, and its response regulator, GacA, is involved in regulation of secondary metabolism and many other aspects of bacterial physiology. Although it is known that the sensor kinases RetS and LadS feed into the GacS/GacA system, the mechanism through which this occurs is unknown, as are the protein-protein interactions in this system. To characterize and define these interactions, we utilized a bacterial two-hybrid system to study the interactions of GacS and GacA from Pseudomonas fluorescens CHA0. Domains of GacA and GacS, identified through bioinformatics, were subcloned and their ability to interact in vivo was investigated. We found that the entire GacA molecule is required for GacA to interact with itself or GacS. Furthermore, the HisKA/HATPase/REC domains of GacS together are responsible for GacS interacting with GacA, while the HAMP domain of GacS is responsible for GacS interacting with itself. In addition, homologs of Pseudomonas aeruginosa hybrid sensor kinases, RetS and LadS, were identified in P. fluorescens, and shown to interact with GacS, but not GacA.

Pseudomonas fluorescens' view of the periodic table

Growth in a biofilm modulates microbial metal susceptibility, sometimes increasing the ability of microorganisms to withstand toxic metal species by several orders of magnitude. In this study, a high-throughput metal toxicity screen was initiated with the aim of correlating biological toxicity data in planktonic and biofilm cells to the physiochemical properties of metal ions. To this end, Pseudomonas fluorescens ATCC 13525 was grown in the Calgary Biofilm Device (CBD) and biofilms and planktonic cells of this microorganism were exposed to gradient arrays of different metal ions. These arrays included 44 different metals with representative compounds that spanned every group of the periodic table (except for the halogens and noble gases). The minimum inhibitory concentration (MIC), minimum bactericidal concentration (MBC) and minimum biofilm eradication concentration (MBEC) values were obtained after exposing the biofilms to metal ions for 4 h. Using these values, metal ion toxicity was correlated to the following ion-specific physicochemical parameters: standard reduction-oxidation potential, electronegativity, the solubility product of the corresponding metal-sulfide complex, the Pearson softness index, electron density and the covalent index. When the ions were grouped according to outer shell electron structure, we found that heavy metal ions gave the strongest correlations to these parameters and were more toxic on average than the other classes of the ions. Correlations were different for biofilms than for planktonic cells, indicating that chemical mechanisms of metal ion toxicity differ between the two modes of growth. We suggest that biofilms can specifically counter the toxic effects of certain physicochemical parameters, which may contribute to the increased ability of biofilms to withstand metal toxicity.

Physical nature of signal peptide binding to DmsD

Here we describe the biophysical characterization of the interaction of the redox enzyme maturation protein DmsD with the signal peptide of its target protein, DmsA. Isothermal titration calorimetry (ITC), size exclusion chromatography (SEC), and an in vitro Far-Western assay is used to show that DmsD binds the twin-arginine signal peptide from DmsA in the micromolar range and in a 1:1 molar ratio. The SEC also shows that there is no oligomerization upon binding. Urea and guanidium hydrochloride denaturation profiles demonstrate the stability of DmsD and give insights on how electrostatic and hydrophobic interactions are important within this binding process. Furthermore, by use of N- and C-terminal fusions of DmsA signal peptide to GST, we observe that N-terminal display of the peptide is important for binding DmsD. In addition, all the folding forms of DmsD were found to bind the DmsA signal peptide as observed with the Far-Western assay.

Twin-arginine translocase may have a role in the chaperone function of NarJ from Escherichia coli

NarJ is a chaperone involved in folding, maturation, and molybdenum cofactor insertion of nitrate reductase A from Escherichia coli. It has also been shown that NarJ exhibits sequence homology to a family of chaperones involved in maturation and cofactor insertion of E. coli redox enzymes that are mediated by twin-arginine translocase (Tat) dependent translocation. In this study, we show that NarJ binds the N-terminal region of NarG through Far Western studies and isothermal titration calorimetry, and the binding event occurs towards a short peptide sequence that contains a homologous twin-arginine motif. Fractionation experiments also show that the interaction of NarJ to the cytoplasmic membrane exhibits Tat-dependence. Upon further investigation through Far Western blots, the interactome of NarJ also exhibits Tat-dependence. Together the data suggest that the Tat system may play a role in the maturation pathway of nitrate reductase A.