A Longitudinal Metagenomic Comparative Analysis of Oral Microbiome Shifts in Patients Receiving Proton Radiation versus Photon Radiation for Head and Neck Cancer

Introduction

Due to the radiation-sparing effects on salivary gland acini, changes in the composition of the oral microbiome may be a driver for improved outcomes in patients receiving proton radiation, with potentially worse outcomes in patients exposed to photon radiation therapy. To date, a head-to-head comparison of oral microbiome changes at a metagenomic level with longitudinal sampling has yet to be performed in these patient cohorts.

Methods and Materials

To comparatively analyze oral microbiome shifts during head and neck radiation therapy, a prospective pilot cohort study was performed at the Maryland Proton Treatment Center and the University of Maryland Marlene and Stewart Greenebaum Comprehensive Cancer Center. A longitudinal metagenomic comparative analysis of oral microbiome shifts was performed at three time points (pre-radiation, during radiation, and immediately post-radiation). Head and neck cancer patients receiving proton radiation (n = 4) were compared to photon radiation (n = 4). Additional control groups included healthy age- and sex-matched controls (n = 5), head and neck cancer patients who never received radiation therapy (n = 8), and patients with oral inflammatory disease (n = 3).

Results

Photon therapy patients presented with lower microbial alpha diversity at all timepoints, and there was a trend towards reduced species richness as compared with proton therapy. Healthy controls and proton patients exhibited overall higher and similar diversity. A more dysbiotic state was observed in patients receiving photon therapy as compared to proton therapy, in which oral microbial homeostasis was maintained. Mucositis was observed in 3/4 photon patients and was not observed in any proton patients during radiation therapy. The bacterial de novo pyrimidine biosynthesis pathway and the nitrate reduction V pathway were comparatively higher following photon exposure. These functional changes in bacterial metabolism may suggest that photon exposure produces a more permissive environment for the proliferation of pathogenic bacteria.

Conclusion

Oral microbiome dysbiosis in patients receiving photon radiation may be associated with increased mucositis occurrence. Proton radiation therapy for head and neck cancer demonstrates a safer side effect profile in terms of oral complications, oral microbiome dysbiosis, and functional metabolic status.

Genomic Confirmation of Borrelia garinii, United States

Lyme disease is a multisystem disorder primarily caused by Borrelia burgdorferi sensu lato. However, B. garinii, which has been identified on islands off the coast of Newfoundland and Labrador, Canada, is a cause of Lyme disease in Eurasia. We report isolation and whole-genome nucleotide sequencing of a B. garinii isolate from a cotton mouse (Peromyscus gossypinus) in South Carolina, USA. We identified a second B. garinii isolate from the same repository. Phylogenetic analysis does not associate these isolates with the previously described isolates of B. garinii from Canada.

Genomic hybrid capture assay to detect Borrelia burgdorferi: an application to diagnose neuroborreliosis in horses

Antemortem diagnosis of neuroborreliosis in horses has been hindered by both the low sensitivity of PCR testing for Borrelia burgdorferi in CSF and the low specificity of serum:CSF ELISA ratios used to determine intrathecal antibody production against the bacterium. PCR testing of the CSF of an adult horse with acute neurologic disease for the B. burgdorferi flagellin gene was negative. However, we enriched B. burgdorferi DNA through nucleic acid hybrid capture, followed by next-generation sequencing, and identified B. burgdorferi in the CSF of the horse, confirming a diagnosis of neuroborreliosis.

Dynamics of the infant gut microbiota in the first 18 months of life: the impact of maternal HIV infection and breastfeeding

BACKGROUND

Access to antiretroviral therapy (ART) during pregnancy and breastfeeding for mothers with HIV has resulted in fewer children acquiring HIV peri- and postnatally, resulting in an increase in the number of children who are exposed to the virus but are not infected (HEU). HEU infants have an increased likelihood of childhood infections and adverse growth outcomes, as well as increased mortality compared to their HIV-unexposed (HUU) peers. We explored potential differences in the gut microbiota in a cohort of 272 Nigerian infants born to HIV-positive and negative mothers in this study during the first 18 months of life.

RESULTS

The taxonomic composition of the maternal vaginal and gut microbiota showed no significant differences based on HIV status, and the composition of the infant gut microbiota at birth was similar between HUU and HEU. Longitudinal taxonomic composition of the infant gut microbiota and weight-for-age z-scores (WAZ) differed depending on access to breast milk. HEU infants displayed overall lower WAZ than HUU infants at all time points. We observed a significantly lower relative abundance of Bifidobacterium in HEU infants at 6 months postpartum. Breast milk composition also differed by time point and HIV infection status. The antiretroviral therapy drugs, lamivudine and nevirapine, as well as kynurenine, were significantly more abundant in the breast milk of mothers with HIV. Levels of tiglyl carnitine (C5) were significantly lower in the breast milk of mothers without HIV. ART drugs in the breast milk of mothers with HIV were associated with a lower relative abundance of Bifidobacterium longum.

CONCLUSIONS

Maternal HIV infection was associated with adverse growth outcomes of HEU infants in this study, and these differences persist from birth through at least 18 months, which is a critical window for the development of the immune and central nervous systems. We observed that the relative abundance of Bifidobacterium spp. was significantly lower in the gut microbiota of all HEU infants over the first 6 months postpartum, even if HEU infants were receiving breast milk. Breastfeeding was of benefit in our HEU infant cohort in the first weeks postpartum; however, ART drug metabolites in breast milk were associated with a lower abundance of Bifidobacterium. Video abstract.

Gut microbiome and metabolome in a non-human primate model of chronic excessive alcohol drinking

A relationship between the gut microbiome and alcohol use disorder has been suggested. Excessive alcohol use produces changes in the fecal microbiome and metabolome in both rodents and humans. Yet, these changes can be observed only in a subgroup of the studied populations, and reversal does not always occur after abstinence. We aimed to analyze fecal microbial composition and function in a translationally relevant baboon model of chronic heavy drinking that also meets binge criteria (drinking too much, too fast, and too often), i.e., alcohol ~1 g/kg and blood alcohol levels (BALs) ≥ 0.08 g/dL in a 2-hour period, daily, for years. We compared three groups of male baboons (Papio anubis): L = Long-term alcohol drinking group (12.1 years); S = Short-term alcohol drinking group (2.7 years); and C = Control group, drinking a non-alcoholic reinforcer (Tang®) (8.2 years). Fecal collection took place during 3 days of Drinking (D), followed by a short period (3 days) of Abstinence (A). Fecal microbial alpha- and beta-diversity were significantly lower in L vs. S and C (p's < 0.05). Members of the commensal families Lachnospiraceae and Prevotellaceae showed a relative decrease, whereas the opportunistic pathogen Streptococcus genus showed a relative increase in L vs. S and C (p's < 0.05). Microbiota-related metabolites of aromatic amino acids, tricarboxylic acid cycle, and pentose increased in L vs. S and C (FDR-corrected p < 0.01), with the latter two suggesting high energy metabolism and enhanced glycolysis in the gut lumen in response to alcohol. Consistent with the long-term alcohol exposure, mucosal damage and oxidative stress markers (N-acetylated amino acids, 2-hydroxybutyrate, and metabolites of the methionine cycle) increased in L vs. S and C (FDR-corrected p < 0.01). Overall, S showed few differences vs. C, possibly due to the long-term, chronic alcohol exposure needed to alter the normal gut microbiota. In the three groups, the fecal microbiome barely differed between conditions D and A, whereas the metabolome shifted in the transition from condition D to A. In conclusion, changes in the fecal microbiome and metabolome occur after significant long-term excessive drinking and are only partially affected by acute forced abstinence from alcohol. These results provide novel information on the relationship between the fecal microbiome and metabolome in a controlled experimental setting and using a unique non-human primate model of chronic excessive alcohol drinking.

Bifidobacterium animalis subsp. lactis BB-12 Protects against Antibiotic-Induced Functional and Compositional Changes in Human Fecal Microbiome

The administration of broad-spectrum antibiotics is often associated with antibiotic-associated diarrhea (AAD), and impacts gastrointestinal tract homeostasis, as evidenced by the following: (a) an overall reduction in both the numbers and diversity of the gut microbiota, and (b) decreased short-chain fatty acid (SCFA) production. Evidence in humans that probiotics may enhance the recovery of microbiota populations after antibiotic treatment is equivocal, and few studies have addressed if probiotics improve the recovery of microbial metabolic function. Our aim was to determine if Bifidobacterium animalis subsp. lactis BB-12 (BB-12)-containing yogurt could protect against antibiotic-induced fecal SCFA and microbiota composition disruptions. We conducted a randomized, allocation-concealed, controlled trial of amoxicillin/clavulanate administration (days 1-7), in conjunction with either BB-12-containing or control yogurt (days 1-14). We measured the fecal levels of SCFAs and bacterial composition at baseline and days 7, 14, 21, and 30. Forty-two participants were randomly assigned to the BB-12 group, and 20 participants to the control group. Antibiotic treatment suppressed the fecal acetate levels in both the control and probiotic groups. Following the cessation of antibiotics, the fecal acetate levels in the probiotic group increased over the remainder of the study and returned to the baseline levels on day 30 (-1.6% baseline), whereas, in the control group, the acetate levels remained suppressed. Further, antibiotic treatment reduced the Shannon diversity of the gut microbiota, for all the study participants at day 7. The magnitude of this change was larger and more sustained in the control group compared to the probiotic group, which is consistent with the hypothesis that BB-12 enhanced microbiota recovery. There were no significant baseline clinical differences between the two groups. Concurrent administration of amoxicillin/clavulanate and BB-12 yogurt, to healthy subjects, was associated with a significantly smaller decrease in the fecal SCFA levels and a more stable taxonomic profile of the microbiota over time than the control group.

Differences in the Platelet mRNA Landscape Portend Racial Disparities in Platelet Function and Suggest Novel Therapeutic Targets

The African American (AA) population displays a 1.6 to 3-fold higher incidence of thrombosis and stroke mortality compared with European Americans (EAs). Current antiplatelet therapies target the ADP-mediated signaling pathway, which displays significant pharmacogenetic variation for platelet reactivity. The focus of this study was to define underlying population differences in platelet function in an effort to identify novel molecular targets for future antiplatelet therapy. We performed deep coverage RNA-Seq to compare gene expression levels in platelets derived from a cohort of healthy volunteers defined by ancestry determination. We identified > 13,000 expressed platelet genes of which 480 were significantly differentially expressed genes (DEGs) between AAs and EAs. DEGs encoding proteins known or predicted to modulate platelet aggregation, morphology, or platelet count were upregulated in AA platelets. Numerous G-protein coupled receptors, ion channels, and pro-inflammatory cytokines not previously associated with platelet function were likewise differentially expressed. Many of the signaling proteins represent potential pharmacologic targets of intervention. Notably, we confirmed the differential expression of cytokines IL32 and PROK2 in an independent cohort by quantitative real-time polymerase chain reaction, and provide functional validation of the opposing actions of these two cytokines on collagen-induced AA platelet aggregation. Using Genotype-Tissue Expression whole blood data, we identified 516 expression quantitative trait locuses with Fst values > 0.25, suggesting that population-differentiated alleles may contribute to differences in gene expression. This study identifies gene expression differences at the population level that may affect platelet function and serve as potential biomarkers to identify cardiovascular disease risk. Additionally, our analysis uncovers candidate novel druggable targets for future antiplatelet therapies.

Transcriptional dysregulation in developing trigeminal sensory neurons in the LgDel mouse model of DiGeorge 22q11.2 deletion syndrome

LgDel mice, which model the heterozygous deletion of genes at human chromosome 22q11.2 associated with DiGeorge/22q11.2 deletion syndrome (22q11DS), have cranial nerve and craniofacial dysfunction as well as disrupted suckling, feeding and swallowing, similar to key 22q11DS phenotypes. Divergent trigeminal nerve (CN V) differentiation and altered trigeminal ganglion (CNgV) cellular composition prefigure these disruptions in LgDel embryos. We therefore asked whether a distinct transcriptional state in a specific population of early differentiating LgDel cranial sensory neurons, those in CNgV, a major source of innervation for appropriate oropharyngeal function, underlies this departure from typical development. LgDel versus wild-type (WT) CNgV transcriptomes differ significantly at E10.5 just after the ganglion has coalesced. Some changes parallel altered proportions of cranial placode versus cranial neural crest-derived CNgV cells. Others are consistent with a shift in anterior-posterior patterning associated with divergent LgDel cranial nerve differentiation. The most robust quantitative distinction, however, is statistically verifiable increased variability of expression levels for most of the over 17 000 genes expressed in common in LgDel versus WT CNgV. Thus, quantitative expression changes of functionally relevant genes and increased stochastic variation across the entire CNgV transcriptome at the onset of CN V differentiation prefigure subsequent disruption of cranial nerve differentiation and oropharyngeal function in LgDel mice.

Changes in the Gut Microbiota Following Bariatric Surgery Are Associated with Increased Alcohol Intake in a Female Rat Model

AIMS

We aimed to investigate if differences in gut microbiota diversity and composition are associated with post-operative alcohol intake following bariatric surgery in a rat model.

METHODS

Twenty-four female rats were randomized to three treatment groups: sham surgery, vertical sleeve gastrectomy (VSG) or Roux-en-Y gastric bypass (RYGB). Stool was collected pre- and post-operatively and 16S rRNA gene amplification and sequencing was performed. Analysis focused on correlating microbial diversity, type of surgery and alcohol (EtOH) intake.

RESULTS

Pre-operative stools samples on regular diet showed similar taxonomic composition and Shannon diversity among the three treatment groups. There was a significant decrease in Shannon diversity and a change in taxonomic composition of the gut microbiota after rats was fed high fat diet. Post-operatively, the RYGB group showed significantly lower taxonomic diversity than the VSG and sham groups, while the VSG and sham groups diversity were not significantly different. Taxonomic composition and function prediction based on PICRUSt analysis showed the RYGB group to be distinct from the VSG and sham groups. Shannon diversity was found to be negatively associated with EtOH intake.

CONCLUSIONS

Changes in the taxonomic profile of the gut microbiota following bariatric surgery, particularly RYGB, are associated with increased EtOH intake and may contribute to increased alcohol use disorder risk through the gut-brain-microbiome axis.

Nasal Microbiota and Infectious Complications After Elective Surgical Procedures

IMPORTANCE

The association of the nasal microbiome with outcomes in surgical patients is poorly understood.

OBJECTIVE

To characterize the composition of nasal microbiota in patients undergoing clean elective surgical procedures and to examine the association between characteristics of preoperative nasal microbiota and occurrence of postoperative infection.

DESIGN, SETTING, AND PARTICIPANTS

Using a nested matched case-control design, 53 individuals who developed postoperative infection were matched (approximately 3:1 by age, sex, and surgical procedure) with 144 individuals who were not infected (ie, the control group). The 2 groups were selected from a prospective cohort of patients undergoing surgical procedures at 2 tertiary care university hospitals in Baltimore, Maryland, who were at high risk for postoperative infectious complications. Included individuals were aged 40 years or older; had no history of autoimmune disease, immunocompromised state, immune-modulating medication, or active infection; and were scheduled to undergo elective cardiac, vascular, spinal, or intracranial surgical procedure. Data were analyzed from October 2015 through September 2020.

EXPOSURES

Nasal microbiome cluster class served as the main exposure. An unsupervised clustering method (ie, grades of membership modeling) was used to classify nasal microbial samples into 2 groups based on features derived from 16S ribosomal RNA gene sequencing. The microbiome cluster groups were derived independently and agnostic of baseline clinical characteristics and infection status.

MAIN OUTCOMES AND MEASURES

Composite of surgical site infection, bacteremia, and pneumonia occurring within 6 months after surgical procedure.

RESULTS

Among 197 participants (mean [SD] age, 64.1 [10.6] years; 63 [37.7%] women), 553 bacterial taxa were identified from preoperative nasal swab samples. A 2-cluster model (with 167 patients in cluster 1 and 30 patients in cluster 2) accounted for the largest proportion of variance in microbial profiles using grades of membership modeling and was most parsimonious. After adjusting for potential confounders, the probability of assignment to cluster 2 was associated with 6-fold higher odds of infection after surgical procedure (odds ratio [OR], 6.18; 95% CI, 3.33-11.7; P < .001) independent of baseline clinical characteristics, including nasal carriage of Staphylococcus aureus. Intrasample (ie, α) diversity was inversely associated with infectious outcome in both clusters (OR, 0.57; 95% CI, 0.42-0.75; P < .001); however, probability of assignment to cluster 2 was associated with higher odds of infection independent of α diversity (OR, 4.61; 95% CI, 2.78-7.86; P < .001).

CONCLUSIONS AND RELEVANCE

These findings suggest that the nasal microbiome was an independent risk factor associated with infectious outcomes among individuals who underwent elective surgical procedures and may serve as a biomarker associated with infection susceptibility in this population.

Cost effective, experimentally robust differential-expression analysis for human/mammalian, pathogen and dual-species transcriptomics

As sequencing read length has increased, researchers have quickly adopted longer reads for their experiments. Here, we examine 14 pathogen or host–pathogen differential gene expression data sets to assess whether using longer reads is warranted. A variety of data sets was used to assess what genomic attributes might affect the outcome of differential gene expression analysis including: gene density, operons, gene length, number of introns/exons and intron length. No genome attribute was found to influence the data in principal components analysis, hierarchical clustering with bootstrap support, or regression analyses of pairwise comparisons that were undertaken on the same reads, looking at all combinations of paired and unpaired reads trimmed to 36, 54, 72 and 101 bp. Read pairing had the greatest effect when there was little variation in the samples from different conditions or in their replicates (e.g. little differential gene expression). But overall, 54 and 72 bp reads were typically most similar. Given differences in costs and mapping percentages, we recommend 54 bp reads for organisms with no or few introns and 72 bp reads for all others. In a third of the data sets, read pairing had absolutely no effect, despite paired reads having twice as much data. Therefore, single-end reads seem robust for differential-expression analyses, but in eukaryotes paired-end reads are likely desired to analyse splice variants and should be preferred for data sets that are acquired with the intent to be community resources that might be used in secondary data analyses.

Molecular mechanisms of probiotic prevention of antibiotic-associated diarrhea

Antibiotic-associated diarrhea (AAD) is a common and unintended adverse effect of antibiotic treatment. It is characterized by the disruption of the gut microbiota, decreased intestinal short chain fatty acid (SCFA) concentrations, accumulation of luminal carbohydrates and colonic bile acids, altered water absorption, and ultimately diarrhea. Probiotics were shown to prevent AAD in numerous clinical trials. This review examines what is currently known about how probiotics reduce the risk for AAD via modulating the gut microbiota, altering nutrient and bile acid metabolism, inducing epithelial solute transporter activity, supporting intestinal barrier function, and influencing the immune system. Although probiotics are frequently prescribed with antibiotic use, mechanistic evidence verifying how they confer protection against AAD is extremely limited. This information is urgently needed for improving recommendations for sustaining probiotic development and for implementing probiotics in clinical settings.

Strains used in whole organism Plasmodium falciparum vaccine trials differ in genome structure, sequence, and immunogenic potential

BACKGROUND

Plasmodium falciparum (Pf) whole-organism sporozoite vaccines have been shown to provide significant protection against controlled human malaria infection (CHMI) in clinical trials. Initial CHMI studies showed significantly higher durable protection against homologous than heterologous strains, suggesting the presence of strain-specific vaccine-induced protection. However, interpretation of these results and understanding of their relevance to vaccine efficacy have been hampered by the lack of knowledge on genetic differences between vaccine and CHMI strains, and how these strains are related to parasites in malaria endemic regions.

METHODS

Whole genome sequencing using long-read (Pacific Biosciences) and short-read (Illumina) sequencing platforms was conducted to generate de novo genome assemblies for the vaccine strain, NF54, and for strains used in heterologous CHMI (7G8 from Brazil, NF166.C8 from Guinea, and NF135.C10 from Cambodia). The assemblies were used to characterize sequences in each strain relative to the reference 3D7 (a clone of NF54) genome. Strains were compared to each other and to a collection of clinical isolates (sequenced as part of this study or from public repositories) from South America, sub-Saharan Africa, and Southeast Asia.

RESULTS

While few variants were detected between 3D7 and NF54, we identified tens of thousands of variants between NF54 and the three heterologous strains. These variants include SNPs, indels, and small structural variants that fall in regulatory and immunologically important regions, including transcription factors (such as PfAP2-L and PfAP2-G) and pre-erythrocytic antigens that may be key for sporozoite vaccine-induced protection. Additionally, these variants directly contributed to diversity in immunologically important regions of the genomes as detected through in silico CD8+ T cell epitope predictions. Of all heterologous strains, NF135.C10 had the highest number of unique predicted epitope sequences when compared to NF54. Comparison to global clinical isolates revealed that these four strains are representative of their geographic origin despite long-term culture adaptation; of note, NF135.C10 is from an admixed population, and not part of recently formed subpopulations resistant to artemisinin-based therapies present in the Greater Mekong Sub-region.

CONCLUSIONS

These results will assist in the interpretation of vaccine efficacy of whole-organism vaccines against homologous and heterologous CHMI.

Establishing What Constitutes a Healthy Human Gut Microbiome: State of the Science, Regulatory Considerations, and Future Directions

On December 17, 2018, the North American branch of the International Life Sciences Institute (ILSI North America) convened a workshop "Can We Begin to Define a Healthy Gut Microbiome Through Quantifiable Characteristics?" with >40 invited academic, government, and industry experts in Washington, DC. The workshop objectives were to 1) develop a collective expert assessment of the state of the evidence on the human gut microbiome and associated human health benefits, 2) see if there was sufficient evidence to establish measurable gut microbiome characteristics that could serve as indicators of "health," 3) identify short- and long-term research needs to fully characterize healthy gut microbiome-host relationships, and 4) publish the findings. Conclusions were as follows: 1) mechanistic links of specific changes in gut microbiome structure with function or markers of human health are not yet established; 2) it is not established if dysbiosis is a cause, consequence, or both of changes in human gut epithelial function and disease; 3) microbiome communities are highly individualized, show a high degree of interindividual variation to perturbation, and tend to be stable over years; 4) the complexity of microbiome-host interactions requires a comprehensive, multidisciplinary research agenda to elucidate relationships between gut microbiome and host health; 5) biomarkers and/or surrogate indicators of host function and pathogenic processes based on the microbiome need to be determined and validated, along with normal ranges, using approaches similar to those used to establish biomarkers and/or surrogate indicators based on host metabolic phenotypes; 6) future studies measuring responses to an exposure or intervention need to combine validated microbiome-related biomarkers and/or surrogate indicators with multiomics characterization of the microbiome; and 7) because static genetic sampling misses important short- and long-term microbiome-related dynamic changes to host health, future studies must be powered to account for inter- and intraindividual variation and should use repeated measures within individuals.

Ghrelin receptor deletion reduces binge-like alcohol drinking in rats

Ghrelin is a gastric hormone that has been implicated in the neurobiology of alcohol drinking. We have recently developed a ghrelin receptor (growth hormone secretagogue receptor; GHSR) knockout (KO) rat model, which exhibits reduced food consumption and body weight. In addition, recent preliminary work suggests that the gut-microbiome, which appears to interact with the ghrelin system, may modulate alcohol drinking. In the present study, we investigated the effects of GHSR deletion on alcohol consumption utilising GHSR KO and wild-type (WT) rats in three separate alcohol consumption paradigms: (i) operant self-administration (30-minute sessions); (ii) drinking in the dark (DID) (4-hour sessions); and (iii) intermittent access (24-hour sessions). These paradigms model varying degrees of alcohol consumption. Furthermore, we aimed to investigate the gut-microbiome composition of GHSR KO and WT rats before and after alcohol exposure. We found that the GHSR KO rats self-administered significantly less alcohol compared to WT rats in the operant paradigm, and consumed less alcohol than WT in the initial stages of the DID paradigm. No genotype differences were found in the intermittent access test. In addition, we found a significant decrease in gut-microbial diversity after alcohol exposure in both genotypes. Thus, the present results indicate that the ghrelin system may be involved in drinking patterns that result in presumably increased alcohol exposure levels. Furthermore, GHSR may constitute a potential pharmacological target for the reduction of binge-alcohol consumption. The potential functional role of the gut-microbiome in alcohol drinking, as well as interaction with the ghrelin system, is an interesting topic for further investigation.

A Distinct Gut Microbiota Exists Within Crohn's Disease-Related Perianal Fistulae

BACKGROUND

Gut bacteria are strongly suspected to play a key role in the pathogenesis of Crohn's disease (CD). Studies have demonstrated alterations in the gut microbiota in this patient population. The purpose of this study was to characterize the gut microbiota of fistulizing perianal CD.

MATERIALS AND METHODS

Stool and fistula samples were obtained from patients undergoing surgery for CD-related anorectal fistulae. Microbial compositions of matched stool and fistula samples were characterized using 16S rRNA gene profiling. The effect of sample type, patient gender, disease classification (Montreal A/B), disease activity (Harvey Bradshaw Index), antibiotic use, and presence of active proctitis on microbial composition was assessed.

RESULTS

Samples were obtained from 18 patients. Bacteroides was the most abundant genera across all samples collected, followed by Streptococcus and Bifidobacterium. Bifidobacterium was present at significantly higher levels in fecal samples than fistula samples, whereas Achromobacter and Corynebacterium were present at significantly higher levels in fistula samples. Antibiotic, but not thiopurine or antitumor necrosis factor medication, exposure affected the gut microbial composition. Patient gender, disease classification, disease activity, and presence of active proctitis did not alter stool or fistula microbiota.

CONCLUSIONS

Our data show that the gut microbiota within CD-related anorectal fistulae is distinct from that in stool samples obtained from the same patients. We also observe a dysbiosis in patients treated with antibiotics compared with those not treated with antibiotics.

Intratumor genetic heterogeneity in squamous cell carcinoma of the oral cavity

BACKGROUND

We sought to evaluate intratumor heterogeneity in squamous cell carcinoma of the oral cavity (OCC) and specifically determine the effect of physical separation and histologic differentiation within the same tumor.

METHODS

We performed whole exome sequencing on five biopsy sites-two from well-differentiated, two from poorly differentiated regions, and one from normal parenchyma-from five primary OCC specimens.

RESULTS

We found high levels of intratumor heterogeneity and, in four primary tumors, identified only 0 to 2 identical mutations in all subsites. We found that the heterogeneity inversely correlated with physical separation and that pairs of well-differentiated samples were more similar to each other than analogous poorly differentiated specimens. Only TP53 mutations, but not other purported "driver mutations" in head and neck squamous cell carcinoma, were found in multiple biopsy sites.

CONCLUSION

These data highlight the challenges to characterization of the mutational landscape of OCC with single site biopsy and have implications for personalized medicine.

Metaproteomics reveals persistent and phylum-redundant metabolic functional stability in adult human gut microbiomes of Crohn's remission patients despite temporal variations in microbial taxa, genomes, and proteomes

BACKGROUND

The gut microbiome plays a fundamental role in the human host's overall health by contributing key biological functions such as expanded metabolism and pathogen defense/immune control. In a healthy individual, the gut microbiome co-exists within the human host in a symbiotic, non-inflammatory relationship that enables mutual benefits, such as microbial degradation of indigestible food products into small molecules that the host can utilize, and enhanced pathogen defense. In abnormal conditions, such as Crohn's disease, this favorable metabolic relationship breaks down and a variety of undesirable activities result, including chronic inflammation and other health-related issues. It has been difficult, however, to elucidate the overall functional characteristics of this relationship because the microbiota can vary substantially in composition for healthy humans and possibly even more in individuals with gut disease conditions such as Crohn's disease. Overall, this suggests that microbial membership composition may not be the best way to characterize a phenotype. Alternatively, it seems to be more informative to examine and characterize the functional composition of a gut microbiome. Towards that end, this study examines 25 metaproteomes measured in several Crohn's disease patients' post-resection surgery across the course of 1 year, in order to examine persistence of microbial taxa, genes, proteins, and metabolic functional distributions across time in individuals whose microbiome might be more variable due to the gut disease condition.

RESULTS

The measured metaproteomes were highly personalized, with all the temporally-related metaproteomes clustering most closely by individual. In general, the metaproteomes were remarkably distinct between individuals and to a lesser extent within individuals. This prompted a need to characterize the metaproteome at a higher functional level, which was achieved by annotating identified protein groups with KEGG orthologous groups to infer metabolic modules. At this level, similar and redundant metabolic functions across multiple phyla were observed across time and between individuals. Tracking through these various metabolic modules revealed a clear path from carbohydrate, lipid, and amino acid degradation to central metabolism and finally the production of fermentation products.

CONCLUSIONS

The human gut metaproteome can vary quite substantially across time and individuals. However, despite substantial intra-individual variation in the metaproteomes, there is a clear persistence of conserved metabolic functions across time and individuals. Additionally, the persistence of these core functions is redundant across multiple phyla but is not always observable in the same sample. Finally, the gut microbiome's metabolism is not driven by a set of discrete linear pathways but a web of interconnected reactions facilitated by a network of enzymes that connect multiple molecules across multiple pathways.

Toward unrestricted use of public genomic data

Publication interests should not limit access to public data

Compositional and Functional Differences in the Human Gut Microbiome Correlate with Clinical Outcome following Infection with Wild-Type Salmonella enterica Serovar Typhi

Insights into disease susceptibility as well as the efficacy of vaccines against typhoid and other enteric pathogens may be informed by better understanding the relationship between the effector immune response and the gut microbiota. In the present study, we characterized the composition (16S rRNA gene profiling) and function (RNA sequencing [RNA-seq]) of the gut microbiota following immunization and subsequent exposure to wild-type Salmonella enterica serovar Typhi in a human challenge model to further investigate the central hypothesis that clinical outcomes may be linked to the gut microbiota. Metatranscriptome analysis of longitudinal stool samples collected from study subjects revealed two stable patterns of gene expression for the human gut microbiota, dominated by transcripts from either Methanobrevibacter or a diverse representation of genera in the Firmicutes phylum. Immunization with one of two live oral attenuated vaccines against S. Typhi had minimal effects on the composition or function of the gut microbiota. It was observed that subjects harboring the methanogen-dominated transcriptome community at baseline displayed a lower risk of developing symptoms of typhoid following challenge with wild-type S. Typhi. Furthermore, genes encoding antioxidant proteins, metal homeostasis and transport proteins, and heat shock proteins were expressed at a higher level at baseline or after challenge with S. Typhi in subjects who did not develop symptoms of typhoid. These data suggest that functional differences relating to redox potential and ion homeostasis in the gut microbiota may impact clinical outcomes following exposure to wild-type S. Typhi.IMPORTANCES. Typhi is a significant cause of systemic febrile morbidity in settings with poor sanitation and limited access to clean water. It has been demonstrated that the human gut microbiota can influence mucosal immune responses, but there is little information available on the impact of the human gut microbiota on clinical outcomes following exposure to enteric pathogens. Here, we describe differences in the composition and function of the gut microbiota in healthy adult volunteers enrolled in a typhoid vaccine trial and report that these differences are associated with host susceptibility to or protection from typhoid after challenge with wild-type S Typhi. Our observations have important implications in interpreting the efficacy of oral attenuated vaccines against enteric pathogens in diverse populations.

Salmonella Typhi Colonization Provokes Extensive Transcriptional Changes Aimed at Evading Host Mucosal Immune Defense During Early Infection of Human Intestinal Tissue

Commensal microorganisms influence a variety of host functions in the gut, including immune response, glucose homeostasis, metabolic pathways and oxidative stress, among others. This study describes how Salmonella Typhi, the pathogen responsible for typhoid fever, uses similar strategies to escape immune defense responses and survive within its human host. To elucidate the early mechanisms of typhoid fever, we performed studies using healthy human intestinal tissue samples and "mini-guts," organoids grown from intestinal tissue taken from biopsy specimens. We analyzed gene expression changes in human intestinal specimens and bacterial cells both separately and after colonization. Our results showed mechanistic strategies that S. Typhi uses to rearrange the cellular machinery of the host cytoskeleton to successfully invade the intestinal epithelium, promote polarized cytokine release and evade immune system activation by downregulating genes involved in antigen sampling and presentation during infection. This work adds novel information regarding S. Typhi infection pathogenesis in humans, by replicating work shown in traditional cell models, and providing new data that can be applied to future vaccine development strategies.

Primordial origin and diversification of plasmids in Lyme disease agent bacteria

Background

With approximately one-third of their genomes consisting of linear and circular plasmids, the Lyme disease agent cluster of species has the most complex genomes among known bacteria. We report here a comparative analysis of plasmids in eleven Borreliella (also known as Borrelia burgdorferi sensu lato) species.

Results

We sequenced the complete genomes of two B. afzelii, two B. garinii, and individual B. spielmanii, B. bissettiae, B. valaisiana and B. finlandensis isolates. These individual isolates carry between seven and sixteen plasmids, and together harbor 99 plasmids. We report here a comparative analysis of these plasmids, along with 70 additional Borreliella plasmids available in the public sequence databases. We identify only one new putative plasmid compatibility type (the 30th) among these 169 plasmid sequences, suggesting that all or nearly all such types have now been discovered. We find that the linear plasmids in the non-B. burgdorferi species have undergone the same kinds of apparently random, chaotic rearrangements mediated by non-homologous recombination that we previously discovered in B. burgdorferi. These rearrangements occurred independently in the different species lineages, and they, along with an expanded chromosomal phylogeny reported here, allow the identification of several whole plasmid transfer events among these species. Phylogenetic analyses of the plasmid partition genes show that a majority of the plasmid compatibility types arose early, most likely before separation of the Lyme agent Borreliella and relapsing fever Borrelia clades, and this, with occasional cross species plasmid transfers, has resulted in few if any species-specific or geographic region-specific Borreliella plasmid types.

Conclusions

The primordial origin and persistent maintenance of the Borreliella plasmid types support their functional indispensability as well as evolutionary roles in facilitating genome diversity. The improved resolution of Borreliella plasmid phylogeny based on conserved partition-gene clusters will lead to better determination of gene orthology which is essential for prediction of biological function, and it will provide a basis for inferring detailed evolutionary mechanisms of Borreliella genomic variability including homologous gene and plasmid exchanges as well as non-homologous rearrangements.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4597-x) contains supplementary material, which is available to authorized users.

Systemic and Terminal Ileum Mucosal Immunity Elicited by Oral Immunization With the Ty21a Typhoid Vaccine in Humans

BACKGROUND & AIMS

Systemic cellular immunity elicited by the Ty21a oral typhoid vaccine has been extensively characterized. However, very limited data are available in humans regarding mucosal immunity at the site of infection (terminal ileum [TI]). Here we investigated the host immunity elicited by Ty21a immunization on terminal ileum-lamina propria mononuclear cells (LPMC) and peripheral blood in volunteers undergoing routine colonoscopy.

METHODS

We characterized LPMC-T memory (T_M) subsets and assessed Salmonella enterica serovar Typhi (S Typhi)-specific responses by multichromatic flow cytometry.

RESULTS

No differences were observed in cell yields and phenotypes in LPMC CD8⁺-T_M subsets following Ty21a immunization. However, Ty21a immunization elicited LPMC CD8⁺ T cells exhibiting significant S Typhi-specific responses (interferon-γ, tumor necrosis factor-α, interleukin-17A, and/or CD107a) in all major T_M subsets (T-effector/memory [T_EM], T-central/memory, and T_EM-CD45RA⁺), although each T_M subset exhibited unique characteristics. We also investigated whether Ty21a immunization elicited S Typhi-specific multifunctional effectors in LPMC CD8⁺ T_EM. We observed that LPMC CD8⁺ T_EM responses were mostly multifunctional, except for those cells exhibiting the characteristics associated with cytotoxic responses. Finally, we compared mucosal with systemic responses and made the important observation that LPMC CD8⁺S Typhi-specific responses were unique and distinct from their systemic counterparts.

CONCLUSIONS

This study provides the first demonstration of S Typhi-specific responses in the human terminal ileum mucosa and provides novel insights into the generation of mucosal immune responses following oral Ty21a immunization.

Aligner optimization increases accuracy and decreases compute times in multi-species sequence data

As sequencing technologies have evolved, the tools to analyze these sequences have made similar advances. However, for multi-species samples, we observed important and adverse differences in alignment specificity and computation time for bwa- mem (Burrows-Wheeler aligner-maximum exact matches) relative to bwa-aln. Therefore, we sought to optimize bwa-mem for alignment of data from multi-species samples in order to reduce alignment time and increase the specificity of alignments. In the multi-species cases examined, there was one majority member (i.e. Plasmodium falciparum or Brugia malayi) and one minority member (i.e. human or the Wolbachia endosymbiont wBm) of the sequence data. Increasing bwa-mem seed length from the default value reduced the number of read pairs from the majority sequence member that incorrectly aligned to the reference genome of the minority sequence member. Combining both source genomes into a single reference genome increased the specificity of mapping, while also reducing the central processing unit (CPU) time. In Plasmodium, at a seed length of 18 nt, 24.1 % of reads mapped to the human genome using 1.7±0.1 CPU hours, while 83.6 % of reads mapped to the Plasmodium genome using 0.2±0.0 CPU hours (total: 107.7 % reads mapping; in 1.9±0.1 CPU hours). In contrast, 97.1 % of the reads mapped to a combined Plasmodium-human reference in only 0.7±0.0 CPU hours. Overall, the results suggest that combining all references into a single reference database and using a 23 nt seed length reduces the computational time, while maximizing specificity. Similar results were found for simulated sequence reads from a mock metagenomic data set. We found similar improvements to computation time in a publicly available human-only data set.

CloVR-Comparative: automated, cloud-enabled comparative microbial genome sequence analysis pipeline

Background

The benefit of increasing genomic sequence data to the scientific community depends on easy-to-use, scalable bioinformatics support. CloVR-Comparative combines commonly used bioinformatics tools into an intuitive, automated, and cloud-enabled analysis pipeline for comparative microbial genomics.

Results

CloVR-Comparative runs on annotated complete or draft genome sequences that are uploaded by the user or selected via a taxonomic tree-based user interface and downloaded from NCBI. CloVR-Comparative runs reference-free multiple whole-genome alignments to determine unique, shared and core coding sequences (CDSs) and single nucleotide polymorphisms (SNPs). Output includes short summary reports and detailed text-based results files, graphical visualizations (phylogenetic trees, circular figures), and a database file linked to the Sybil comparative genome browser. Data up- and download, pipeline configuration and monitoring, and access to Sybil are managed through CloVR-Comparative web interface. CloVR-Comparative and Sybil are distributed as part of the CloVR virtual appliance, which runs on local computers or the Amazon EC2 cloud. Representative datasets (e.g. 40 draft and complete Escherichia coli genomes) are processed in <36 h on a local desktop or at a cost of <$20 on EC2.

Conclusions

CloVR-Comparative allows anybody with Internet access to run comparative genomics projects, while eliminating the need for on-site computational resources and expertise.

Electronic supplementary material

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

New var reconstruction algorithm exposes high var sequence diversity in a single geographic location in Mali

Background

Encoded by the var gene family, highly variable Plasmodium falciparum erythrocyte membrane protein-1 (PfEMP1) proteins mediate tissue-specific cytoadherence of infected erythrocytes, resulting in immune evasion and severe malaria disease. Sequencing and assembling the 40–60 var gene complement for individual infections has been notoriously difficult, impeding molecular epidemiological studies and the assessment of particular var elements as subunit vaccine candidates.

Methods

We developed and validated a novel algorithm, Exon-Targeted Hybrid Assembly (ETHA), to perform targeted assembly of var gene sequences, based on a combination of Pacific Biosciences and Illumina data.

Results

Using ETHA, we characterized the repertoire of var genes in 12 samples from uncomplicated malaria infections in children from a single Malian village and showed them to be as genetically diverse as vars from isolates from around the globe. The gene var2csa, a member of the var family associated with placental malaria pathogenesis, was present in each genome, as were vars previously associated with severe malaria.

Conclusion

ETHA, a tool to discover novel var sequences from clinical samples, will aid the understanding of malaria pathogenesis and inform the design of malaria vaccines based on PfEMP1.

ETHA is available at: https://sourceforge.net/projects/etha/.

Electronic supplementary material

The online version of this article (doi:10.1186/s13073-017-0422-4) contains supplementary material, which is available to authorized users.

Plasmid diversity and phylogenetic consistency in the Lyme disease agent Borrelia burgdorferi

Background

Bacteria from the genus Borrelia are known to harbor numerous linear and circular plasmids. We report here a comparative analysis of the nucleotide sequences of 236 plasmids present in fourteen independent isolates of the Lyme disease agent B. burgdorferi.

Results

We have sequenced the genomes of 14 B. burgdorferi sensu stricto isolates that carry a total of 236 plasmids. These individual isolates carry between seven and 23 plasmids. Their chromosomes, the cp26 and cp32 circular plasmids, as well as the lp54 linear plasmid, are quite evolutionarily stable; however, the remaining plasmids have undergone numerous non-homologous and often duplicative recombination events. We identify 32 different putative plasmid compatibility types among the 236 plasmids, of which 15 are (usually) circular and 17 are linear. Because of past rearrangements, any given gene, even though it might be universally present in these isolates, is often found on different linear plasmid compatibility types in different isolates. For example, the arp gene and the vls cassette region are present on plasmids of four and five different compatibility types, respectively, in different isolates. A majority of the plasmid types have more than one organizationally different subtype, and the number of such variants ranges from one to eight among the 18 linear plasmid types. In spite of this substantial organizational diversity, the plasmids are not so variable that every isolate has a novel version of every plasmid (i.e., there appears to be a limited number of extant plasmid subtypes).

Conclusions

Although there have been many past recombination events, both homologous and nonhomologous, among the plasmids, particular organizational variants of these plasmids correlate with particular chromosomal genotypes, suggesting that there has not been rapid horizontal transfer of whole linear plasmids among B. burgdorferi lineages. We argue that plasmid rearrangements are essentially non-revertable and are present at a frequency of only about 0.65% that of single nucleotide changes, making rearrangement-derived novel junctions (mosaic boundaries) ideal phylogenetic markers in the study of B. burgdorferi population structure and plasmid evolution and exchange.

Electronic supplementary material

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

Genome-wide diversity and gene expression profiling of Babesia microti isolates identify polymorphic genes that mediate host-pathogen interactions

Babesia microti, a tick-transmitted, intraerythrocytic protozoan parasite circulating mainly among small mammals, is the primary cause of human babesiosis. While most cases are transmitted by Ixodes ticks, the disease may also be transmitted through blood transfusion and perinatally. A comprehensive analysis of genome composition, genetic diversity, and gene expression profiling of seven B. microti isolates revealed that genetic variation in isolates from the Northeast United States is almost exclusively associated with genes encoding the surface proteome and secretome of the parasite. Furthermore, we found that polymorphism is restricted to a small number of genes, which are highly expressed during infection. In order to identify pathogen-encoded factors involved in host-parasite interactions, we screened a proteome array comprised of 174 B. microti proteins, including several predicted members of the parasite secretome. Using this immuno-proteomic approach we identified several novel antigens that trigger strong host immune responses during the onset of infection. The genomic and immunological data presented herein provide the first insights into the determinants of B. microti interaction with its mammalian hosts and their relevance for understanding the selective pressures acting on parasite evolution.

Efficient Enrichment of Bacterial mRNA from Host-Bacteria Total RNA Samples

Despite numerous advances in genomics and bioinformatics, technological hurdles remain to examine host-microbe transcriptomics. Sometimes the transcriptome of either or both can be ascertained merely by generating more sequencing reads. However, many cases exist where bacterial mRNA needs to be enriched further to enable cost-effective sequencing of the pathogen or endosymbiont. While a suitable method is commercially available for mammalian samples of this type, development of such methods has languished for invertebrate samples. Furthermore, a common method across multiple taxa would facilitate comparisons between bacteria in invertebrate vectors and their vertebrate hosts. Here, a method is described to concurrently remove polyadenylated transcripts, prokaryotic rRNA, and eukaryotic rRNA, including those with low amounts of starting material (e.g. 100 ng). In a Wolbachia-Drosophila system, this bacterial mRNA enrichment yielded a 3-fold increase in Wolbachia mRNA abundance and a concomitant 3.3-fold increase in the percentage of transcripts detected. More specifically, 70% of the genome could be recovered by transcriptome sequencing compared to 21% in the total RNA. Sequencing of similar bacterial mRNA-enriched samples generated from Ehrlichia-infected canine cells covers 93% of the Ehrlichia genome, suggesting ubiquitous transcription across the entire Ehrlichia chaffeensis genome. This technique can potentially be used to enrich bacterial mRNA in many studies of host-microbe interactions.

Annotated draft genome sequences of three species of Cryptosporidium: Cryptosporidium meleagridis isolate UKMEL1, C. baileyi isolate TAMU-09Q1 and C. hominis isolates TU502_2012 and UKH1

Human cryptosporidiosis is caused primarily by Cryptosporidium hominis, C. parvum and C. meleagridis. To accelerate research on parasites in the genus Cryptosporidium, we generated annotated, draft genome sequences of human C. hominis isolates TU502_2012 and UKH1, C. meleagridis UKMEL1, also isolated from a human patient, and the avian parasite C. baileyi TAMU-09Q1. The annotation of the genome sequences relied in part on RNAseq data generated from the oocyst stage of both C. hominis and C. baileyi. The genome assembly of C. hominis is significantly more complete and less fragmented than that available previously, which enabled the generation of a much-improved gene set for this species, with an increase in average gene length of 500 bp relative to the protein-encoding genes in the 2004 C. hominis annotation. Our results reveal that the genomes of C. hominis and C. parvum are very similar in both gene density and average gene length. These data should prove a valuable resource for the Cryptosporidium research community.

The release of the draft genome sequence, and corresponding annotation, of Cryptosporidium baileyi, C. hominis isolates TU502_2012 and UKH1, and C. meleagridis, will accelerate research on Cryptosporidium parasites.

An integrated genomic and transcriptomic survey of mucormycosis-causing fungi

Mucormycosis is a life-threatening infection caused by Mucorales fungi. Here we sequence 30 fungal genomes, and perform transcriptomics with three representative Rhizopus and Mucor strains and with human airway epithelial cells during fungal invasion, to reveal key host and fungal determinants contributing to pathogenesis. Analysis of the host transcriptional response to Mucorales reveals platelet-derived growth factor receptor B (PDGFRB) signaling as part of a core response to divergent pathogenic fungi; inhibition of PDGFRB reduces Mucorales-induced damage to host cells. The unique presence of CotH invasins in all invasive Mucorales, and the correlation between CotH gene copy number and clinical prevalence, are consistent with an important role for these proteins in mucormycosis pathogenesis. Our work provides insight into the evolution of this medically and economically important group of fungi, and identifies several molecular pathways that might be exploited as potential therapeutic targets.

Genomic insights into the Ixodes scapularis tick vector of Lyme disease

Ticks transmit more pathogens to humans and animals than any other arthropod. We describe the 2.1 Gbp nuclear genome of the tick, Ixodes scapularis (Say), which vectors pathogens that cause Lyme disease, human granulocytic anaplasmosis, babesiosis and other diseases. The large genome reflects accumulation of repetitive DNA, new lineages of retro-transposons, and gene architecture patterns resembling ancient metazoans rather than pancrustaceans. Annotation of scaffolds representing ∼57% of the genome, reveals 20,486 protein-coding genes and expansions of gene families associated with tick-host interactions. We report insights from genome analyses into parasitic processes unique to ticks, including host 'questing', prolonged feeding, cuticle synthesis, blood meal concentration, novel methods of haemoglobin digestion, haem detoxification, vitellogenesis and prolonged off-host survival. We identify proteins associated with the agent of human granulocytic anaplasmosis, an emerging disease, and the encephalitis-causing Langat virus, and a population structure correlated to life-history traits and transmission of the Lyme disease agent.

Pulmonary Nontuberculous Mycobacterial Infection. A Multisystem, Multigenic Disease

RATIONALE

The clinical features of patients infected with pulmonary nontuberculous mycobacteria (PNTM) are well described, but the genetic components of infection susceptibility are not.

OBJECTIVES

To examine genetic variants in patients with PNTM, their unaffected family members, and a control group.

METHODS

Whole-exome sequencing was done on 69 white patients with PNTM and 18 of their white unaffected family members. We performed a candidate gene analysis using immune, cystic fibrosis transmembrance conductance regulator (CFTR), cilia, and connective tissue gene sets. The numbers of patients, family members, and control subjects with variants in each category were compared, as was the average number of variants per person.

MEASUREMENTS AND MAIN RESULTS

A significantly higher number of patients with PNTM than the other subjects had low-frequency, protein-affecting variants in immune, CFTR, cilia, and connective tissue categories (35, 26, 90, and 90%, respectively). Patients with PNTM also had significantly more cilia and connective tissue variants per person than did control subjects (2.47 and 2.55 compared with 1.38 and 1.40, respectively; P = 1.4 × 10(-6) and P = 2.7 × 10(-8), respectively). Patients with PNTM had an average of 5.26 variants across all categories (1.98 in control subjects; P = 2.8 × 10(-17)), and they were more likely than control subjects to have variants in multiple categories. We observed similar results for family members without PNTM infection, with the exception of the immune category.

CONCLUSIONS

Patients with PNTM have more low-frequency, protein-affecting variants in immune, CFTR, cilia, and connective tissue genes than their unaffected family members and control subjects. We propose that PNTM infection is a multigenic disease in which combinations of variants across gene categories, plus environmental exposures, increase susceptibility to the infection.

The Oral Bacterial Communities of Children with Well-Controlled HIV Infection and without HIV Infection

The oral microbial community (microbiota) plays a critical role in human health and disease. Alterations in the oral microbiota may be associated with disorders such as gingivitis, periodontitis, childhood caries, alveolar osteitis, oral candidiasis and endodontic infections. In the immunosuppressed population, the spectrum of potential oral disease is even broader, encompassing candidiasis, necrotizing gingivitis, parotid gland enlargement, Kaposi’s sarcoma, oral warts and other diseases. Here, we used 454 pyrosequencing of bacterial 16S rRNA genes to examine the oral microbiome of saliva, mucosal and tooth samples from HIV-positive and negative children. Patient demographics and clinical characteristics were collected from a cross-section of patients undergoing routine dental care. Multiple specimens from different sampling sites in the mouth were collected for each patient. The goal of the study was to observe the potential diversity of the oral microbiota among individual patients, sample locations, HIV status and various dental characteristics. We found that there were significant differences in the microbiome among the enrolled patients, and between sampling locations. The analysis was complicated by uneven enrollment in the patient cohorts, with only five HIV-negative patients enrolled in the study and by the rapid improvement in the health of HIV-infected children between the time the study was conceived and completed. The generally good oral health of the HIV-negative patients limited the number of dental plaque samples that could be collected. We did not identify significant differences between well-controlled HIV-positive patients and HIV-negative controls, suggesting that well-controlled HIV-positive patients essentially harbor similar oral flora compared to patients without HIV. Nor were significant differences in the oral microbiota identified between different teeth or with different dental characteristics. Additional studies are needed to better characterize the oral microbiome in children and those with poorly-controlled HIV infections.

Potential Influence of Staphylococcus aureus Clonal Complex 30 Genotype and Transcriptome on Hematogenous Infections

Background.  The contemporary Staphylococcus aureus clonal complex (CC) 30 lineage is associated with complicated infections, including endocarditis and osteomyelitis. This lineage diverged from the phage-type 80/81 S aureus clone responsible for a major bacterial epidemic of the 20th century. The genome and transcriptome features that contribute to complicated infections of the CC30 lineage are unknown. Methods.  Twenty-nine clinical methicillin-resistant S aureus (MRSA) strains (8 from CC30 and 21 from other major CCs were evaluated for virulence using murine and Galleria mellonella sepsis models. Genomic features of CC30 were identified by comparative genome sequencing and RNA-Seq transcriptome analysis of the 29 strains and 31 previously sequenced S aureus genomes. Results.  The CC30 isolates displayed lower virulence in the sepsis models compared with other CCs [P < .0001]. Comparisons of orthologous proteins and transcriptome analysis identified genes (eg, nitric oxide reductase) and changes in metabolic pathways (eg, pyrimidine metabolism) that contribute to the distinct CC30 phenotype. Previously reported nonsynonymous single-nucleotide polymorphisms (SNPs) were found in accessory gene regulator C (agrC) and α-hemolysin (hla), molecules important for virulence. Additional nonsynonymous SNPs conserved across clinical CC30 isolates when compared with the first sequenced contemporary CC30 clone, MRSA-16, were identified in multiple genes, suggesting continuing evolutionary divergence in this lineage. Conclusions.  Genomic and transcriptional analyses suggest that the CC30 lineage has acquired metabolic features that contribute to persistent and complicated infections. Absence of sepsis-induced mortality in animal models may be due in part to its unique genomic profile and suggests that specific genotypes of S aureus elicit distinct types of infection types.

Transcriptional attenuation controls macrolide inducible efflux and resistance in Streptococcus pneumoniae and in other Gram-positive bacteria containing mef/mel(msr(D)) elements

Macrolide resistance, emerging in Streptococcus pneumoniae and other Gram-positive bacteria, is increasingly due to efflux pumps encoded by mef/mel(msr) operons found on discrete mobile genetic elements. The regulation of mef/mel(msr) in these elements is not well understood. We identified the mef(E)/mel transcriptional start, localized the mef(E)/mel promoter, and demonstrated attenuation of transcription as a mechanism of regulation of macrolide-inducible mef-mediated macrolide resistance in S. pneumoniae. The mef(E)/mel transcriptional start site was a guanine 327 bp upstream of mef(E). Consensus pneumococcal promoter -10 (5′-TATACT-3′) and -35 (5′-TTGAAC-3′) boxes separated by 17 bp were identified 7 bp upstream of the start site. Analysis of the predicted secondary structure of the 327 5’ region identified four pairs of inverted repeats R1-R8 predicted to fold into stem-loops, a small leader peptide [MTASMRLR, (Mef(E)L)] required for macrolide induction and a Rho-independent transcription terminator. RNA-seq analyses provided confirmation of transcriptional attenuation. In addition, expression of mef(E)L was also influenced by mef(E)L-dependent mRNA stability. The regulatory region 5’ of mef(E) was highly conserved in other mef/mel(msr)-containing elements including Tn1207.1 and the 5612IQ complex in pneumococci and Tn1207.3 in Group A streptococci, indicating a regulatory mechanism common to a wide variety of Gram-positive bacteria containing mef/mel(msr) elements.

High-level relatedness among Mycobacterium abscessus subsp. massiliense strains from widely separated outbreaks

Three recently sequenced strains isolated from patients during an outbreak of Mycobacterium abscessus subsp. massiliense infections at a cystic fibrosis center in the United States were compared with 6 strains from an outbreak at a cystic fibrosis center in the United Kingdom and worldwide strains. Strains from the 2 cystic fibrosis outbreaks showed high-level relatedness with each other and major-level relatedness with strains that caused soft tissue infections during an epidemic in Brazil. We identified unique single-nucleotide polymorphisms in cystic fibrosis and soft tissue outbreak strains, separate single-nucleotide polymorphisms only in cystic fibrosis outbreak strains, and unique genomic traits for each subset of isolates. Our findings highlight the necessity of identifying M. abscessus to the subspecies level and screening all cystic fibrosis isolates for relatedness to these outbreak strains. We propose 2 diagnostic strategies that use partial sequencing of rpoB and secA1 genes and a multilocus sequence typing protocol.

Single molecule sequencing and genome assembly of a clinical specimen of Loa loa, the causative agent of loiasis

Background

More than 20% of the world’s population is at risk for infection by filarial nematodes and >180 million people worldwide are already infected. Along with infection comes significant morbidity that has a socioeconomic impact. The eight filarial nematodes that infect humans are Wuchereria bancrofti, Brugia malayi, Brugia timori, Onchocerca volvulus, Loa loa, Mansonella perstans, Mansonella streptocerca, and Mansonella ozzardi, of which three have published draft genome sequences. Since all have humans as the definitive host, standard avenues of research that rely on culturing and genetics have often not been possible. Therefore, genome sequencing provides an important window into understanding the biology of these parasites. The need for large amounts of high quality genomic DNA from homozygous, inbred lines; the availability of only short sequence reads from next-generation sequencing platforms at a reasonable expense; and the lack of random large insert libraries has limited our ability to generate high quality genome sequences for these parasites. However, the Pacific Biosciences single molecule, real-time sequencing platform holds great promise in reducing input amounts and generating sufficiently long sequences that bypass the need for large insert paired libraries.

Results

Here, we report on efforts to generate a more complete genome assembly for L. loa using genetically heterogeneous DNA isolated from a single clinical sample and sequenced on the Pacific Biosciences platform. To obtain the best assembly, numerous assemblers and sequencing datasets were analyzed, combined, and compared. Quiver-informed trimming of an assembly of only Pacific Biosciences reads by HGAP2 was selected as the final assembly of 96.4 Mbp in 2,250 contigs. This results in ~9% more of the genome in ~85% fewer contigs from ~80% less starting material at a fraction of the cost of previous Roche 454-based sequencing efforts.

Conclusions

The result is the most complete filarial nematode assembly produced thus far and demonstrates the utility of single molecule sequencing on the Pacific Biosciences platform for genetically heterogeneous metazoan genomes.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-788) contains supplementary material, which is available to authorized users.

BorreliaBase: a phylogeny-centered browser of Borrelia genomes

Background

The bacterial genus Borrelia (phylum Spirochaetes) consists of two groups of pathogens represented respectively by B. burgdorferi, the agent of Lyme borreliosis, and B. hermsii, the agent of tick-borne relapsing fever. The number of publicly available Borrelia genomic sequences is growing rapidly with the discovery and sequencing of Borrelia strains worldwide. There is however a lack of dedicated online databases to facilitate comparative analyses of Borrelia genomes.

Description

We have developed BorreliaBase, an online database for comparative browsing of Borrelia genomes. The database is currently populated with sequences from 35 genomes of eight Lyme-borreliosis (LB) group Borrelia species and 7 Relapsing-fever (RF) group Borrelia species. Distinct from genome repositories and aggregator databases, BorreliaBase serves manually curated comparative-genomic data including genome-based phylogeny, genome synteny, and sequence alignments of orthologous genes and intergenic spacers.

Conclusions

With a genome phylogeny at its center, BorreliaBase allows online identification of hypervariable lipoprotein genes, potential regulatory elements, and recombination footprints by providing evolution-based expectations of sequence variability at each genomic locus. The phylo-centric design of BorreliaBase (http://borreliabase.org) is a novel model for interactive browsing and comparative analysis of bacterial genomes online.

Culture-independent evaluation of the appendix and rectum microbiomes in children with and without appendicitis

PURPOSE

The function of the appendix is largely unknown, but its microbiota likely contributes to function. Alterations in microbiota may contribute to appendicitis, but conventional culture studies have not yielded conclusive information. We conducted a pilot, culture-independent 16S rRNA-based microbiota study of paired appendix and rectal samples.

METHODS

We collected appendix and rectal swabs from 21 children undergoing appendectomy, six with normal appendices and fifteen with appendicitis (nine perforated). After DNA extraction, we amplified and sequenced 16S rRNA genes and analyzed sequences using CLoVR. We identified organisms differing in relative abundance using ANOVA (p<0.05) by location (appendix vs. rectum), disease (appendicitis vs. normal), and disease severity (perforated vs. non-perforated).

RESULTS

We identified 290 taxa in the study's samples. Three taxa were significantly increased in normal appendices vs. normal rectal samples: Fusibacter (p = 0.009), Selenomonas (p = 0.026), and Peptostreptococcus (p = 0.049). Five taxa were increased in abundance in normal vs. diseased appendices: Paenibacillaceae (p = 0.005), Acidobacteriaceae GP4 (p = 0.019), Pseudonocardinae (p = 0.019), Bergeyella (p = 0.019) and Rhizobium (p = 0.045). Twelve taxa were increased in the appendices of appendicitis patients vs. normal appendix: Peptostreptococcus (p = 0.0003), Bilophila (p = 0.0004), Bulleidia (p = 0.012), Fusobacterium (p = 0.018), Parvimonas (p = 0.003), Mogibacterium (p = 0.012), Aminobacterium (p = 0.019), Proteus (p = 0.028), Actinomycineae (p = 0.028), Anaerovorax (p = 0.041), Anaerofilum (p = 0.045), Porphyromonas (p = 0.010). Five taxa were increased in appendices in patients with perforated vs. nonperforated appendicitis: Bulleidia (p = 0.004), Fusibacter (p = 0.005), Prevotella (p = 0.021), Porphyromonas (p = 0.030), Dialister (p = 0.035). Three taxa were increased in rectum samples of patients with appendicitis compared to the normal patients: Bulleidia (p = 0.034), Dialister (p = 0.003), and Porphyromonas (p = 0.026).

CONCLUSION

Specific taxa are more abundant in normal appendices compared to the rectum, suggesting that a distinctive appendix microbiota exists. Taxa with altered abundance in diseased and severely diseased (perforated) samples may contribute to appendicitis pathogenesis, and may provide microbial signatures in the rectum useful for guiding both treatment and diagnosis of appendicitis.

Probiotics: achieving a better regulatory fit

The development and marketing of new probiotic products, substances containing live microorganisms that have a beneficial effect on the human body, have dramatically increased over the last few years. This article examines how the Food and Drug Administration and Federal Trade Commission currently regulate probiotics and makes recommendations as to changes that might be made to ensure that probiotic products are made available to the general public in a way that is both safe and effective.

Simultaneous transcriptional profiling of bacteria and their host cells

We developed an RNA-Seq-based method to simultaneously capture prokaryotic and eukaryotic expression profiles of cells infected with intracellular bacteria. As proof of principle, this method was applied to Chlamydia trachomatis-infected epithelial cell monolayers in vitro, successfully obtaining transcriptomes of both C. trachomatis and the host cells at 1 and 24 hours post-infection. Chlamydiae are obligate intracellular bacterial pathogens that cause a range of mammalian diseases. In humans chlamydiae are responsible for the most common sexually transmitted bacterial infections and trachoma (infectious blindness). Disease arises by adverse host inflammatory reactions that induce tissue damage & scarring. However, little is known about the mechanisms underlying these outcomes. Chlamydia are genetically intractable as replication outside of the host cell is not yet possible and there are no practical tools for routine genetic manipulation, making genome-scale approaches critical. The early timeframe of infection is poorly understood and the host transcriptional response to chlamydial infection is not well defined. Our simultaneous RNA-Seq method was applied to a simplified in vitro model of chlamydial infection. We discovered a possible chlamydial strategy for early iron acquisition, putative immune dampening effects of chlamydial infection on the host cell, and present a hypothesis for Chlamydia-induced fibrotic scarring through runaway positive feedback loops. In general, simultaneous RNA-Seq helps to reveal the complex interplay between invading bacterial pathogens and their host mammalian cells and is immediately applicable to any bacteria/host cell interaction.

Inter- and intra-specific pan-genomes of Borrelia burgdorferi sensu lato: genome stability and adaptive radiation

BACKGROUND

Lyme disease is caused by spirochete bacteria from the Borrelia burgdorferi sensu lato (B. burgdorferi s.l.) species complex. To reconstruct the evolution of B. burgdorferi s.l. and identify the genomic basis of its human virulence, we compared the genomes of 23 B. burgdorferi s.l. isolates from Europe and the United States, including B. burgdorferi sensu stricto (B. burgdorferi s.s., 14 isolates), B. afzelii (2), B. garinii (2), B. "bavariensis" (1), B. spielmanii (1), B. valaisiana (1), B. bissettii (1), and B. "finlandensis" (1).

RESULTS

Robust B. burgdorferi s.s. and B. burgdorferi s.l. phylogenies were obtained using genome-wide single-nucleotide polymorphisms, despite recombination. Phylogeny-based pan-genome analysis showed that the rate of gene acquisition was higher between species than within species, suggesting adaptive speciation. Strong positive natural selection drives the sequence evolution of lipoproteins, including chromosomally-encoded genes 0102 and 0404, cp26-encoded ospC and b08, and lp54-encoded dbpA, a07, a22, a33, a53, a65. Computer simulations predicted rapid adaptive radiation of genomic groups as population size increases.

CONCLUSIONS

Intra- and inter-specific pan-genome sizes of B. burgdorferi s.l. expand linearly with phylogenetic diversity. Yet gene-acquisition rates in B. burgdorferi s.l. are among the lowest in bacterial pathogens, resulting in high genome stability and few lineage-specific genes. Genome adaptation of B. burgdorferi s.l. is driven predominantly by copy-number and sequence variations of lipoprotein genes. New genomic groups are likely to emerge if the current trend of B. burgdorferi s.l. population expansion continues.

Genome sequences of 65 Helicobacter pylori strains isolated from asymptomatic individuals and patients with gastric cancer, peptic ulcer disease, or gastritis

Helicobacter pylori, inhabitant of the gastric mucosa of over half of the world population, with decreasing prevalence in the U.S., has been associated with a variety of gastric pathologies. However, the majority of H. pylori-infected individuals remain asymptomatic, and negative correlations between H. pylori and allergic diseases have been reported. Comprehensive genome characterization of H. pylori populations from different human host backgrounds including healthy individuals provides the exciting potential to generate new insights into the open question whether human health outcome is associated with specific H. pylori genotypes or dependent on other environmental factors. We report the genome sequences of 65 H. pylori isolates from individuals with gastric cancer, preneoplastic lesions, peptic ulcer disease, gastritis, and from asymptomatic adults. Isolates were collected from multiple locations in North America (USA and Canada) as well as from Columbia and Japan. The availability of these H. pylori genome sequences from individuals with distinct clinical presentations provides the research community with a resource for detailed investigations into genetic elements that correlate either positively or negatively with the epidemiology, human host adaptation, and gastric pathogenesis and will aid in the characterization of strains that may favor the development of specific pathology, including gastric cancer.

Differential response of the cynomolgus macaque gut microbiota to Shigella infection

Little is known about the role of gut microbiota in response to live oral vaccines against enteric pathogens. We examined the effect of immunization with an oral live-attenuated Shigella dysenteriae 1 vaccine and challenge with wild-type S. dysenteriae 1 on the fecal microbiota of cynomolgus macaques using 16 S rRNA analysis of fecal samples. Multi-dimensional cluster analysis identified different bacterial community types within macaques from geographically distinct locations. The fecal microbiota of Mauritian macaques, observed to be genetically distinct, harbored a high-diversity community and responded differently to Shigella immunization, as well as challenge compared to the microbiota in non-Mauritian macaques. While both macaque populations exhibited anti-Shigella antibody responses, clinical shigellosis was observed only among non-Mauritian macaques. These studies highlight the importance of further investigation into the possible protective role of the microbiota against enteric pathogens and consideration of host genetic backgrounds in conducting vaccine studies.

Impact of oral typhoid vaccination on the human gut microbiota and correlations with s. Typhi-specific immunological responses

The resident microbial consortia of the human gastrointestinal tract play an integral role in modulating immune responses both locally and systemically. However, detailed information regarding the effector immune responses after vaccine administration in relation to the gastrointestinal microbiota is absent. In this study, the licensed oral live-attenuated typhoid vaccine Ty21a was administered in a clinical study to investigate whether oral immunization resulted in alterations of the microbiota and to identify whether a given microbiota composition, or subsets of the community, are associated with defined S. Typhi-specific immunological responses. The fecal microbiota composition and temporal dynamics were characterized using bacterial 16S rRNA pyrosequencing from individuals who were either immunized with the Ty21a typhoid vaccine (n = 13) or served as unvaccinated controls (n = 4). The analysis revealed considerable inter- and intra-individual variability, yet no discernible perturbations of the bacterial assemblage related to vaccine administration were observed. S. Typhi-specific cell mediated immune (CMI) responses were evaluated by measurement of intracellular cytokine production using multiparametric flow cytometry, and humoral responses were evaluated by measurement of serum anti-LPS IgA and IgG titers. Volunteers were categorized according to the kinetics and magnitude of their responses. While differences in microbial composition, diversity, or temporal stability were not observed among individuals able to mount a positive humoral response, individuals displaying multiphasic CMI responses harbored more diverse, complex communities. In line with this preliminary observation, over two hundred operational taxonomic units (OTUs) were found to differentiate multiphasic and late CMI responders, the vast majority of which classified within the order Clostridiales. These results provide an unprecedented view into the dramatic temporal heterogeneity of both the gut microbiota and host immune responses.

Integrated metagenomics/metaproteomics reveals human host-microbiota signatures of Crohn's disease

Crohn's disease (CD) is an inflammatory bowel disease of complex etiology, although dysbiosis of the gut microbiota has been implicated in chronic immune-mediated inflammation associated with CD. Here we combined shotgun metagenomic and metaproteomic approaches to identify potential functional signatures of CD in stool samples from six twin pairs that were either healthy, or that had CD in the ileum (ICD) or colon (CCD). Integration of these omics approaches revealed several genes, proteins, and pathways that primarily differentiated ICD from healthy subjects, including depletion of many proteins in ICD. In addition, the ICD phenotype was associated with alterations in bacterial carbohydrate metabolism, bacterial-host interactions, as well as human host-secreted enzymes. This eco-systems biology approach underscores the link between the gut microbiota and functional alterations in the pathophysiology of Crohn's disease and aids in identification of novel diagnostic targets and disease specific biomarkers.