Efficient computation of contributional diversity metrics from microbiome data with FuncDiv.

MOTIVATION: Microbiome datasets with taxa linked to the functions (e.g. genes) they encode are becoming more common as metagenomics sequencing approaches improve. However, these data are challenging to analyze due to their complexity. Summary metrics, such as the alpha and beta diversity of taxa contributing to each function (i.e. contributional diversity), represent one approach to investigate these data, but currently there are no straightforward methods for doing so. RESULTS: We addressed this gap by developing FuncDiv, which efficiently performs these computations. Contributional diversity metrics can provide novel insights that would be impossible to identify without jointly considering taxa and functions. AVAILABILITY AND IMPLEMENTATION: FuncDiv is distributed under a GNU Affero General Public License v3.0 and is available at https://github.com/gavinmdouglas/FuncDiv. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Cross-study analyses of microbial abundance using generalized common factor methods.

BACKGROUND: By creating networks of biochemical pathways, communities of micro-organisms are able to modulate the properties of their environment and even the metabolic processes within their hosts. Next-generation high-throughput sequencing has led to a new frontier in microbial ecology, promising the ability to leverage the microbiome to make crucial advancements in the environmental and biomedical sciences. However, this is challenging, as genomic data are high-dimensional, sparse, and noisy. Much of this noise reflects the exact conditions under which sequencing took place, and is so significant that it limits consensus-based validation of study results. RESULTS: We propose an ensemble approach for cross-study exploratory analyses of microbial abundance data in which we first estimate the variance-covariance matrix of the underlying abundances from each dataset on the log scale assuming Poisson sampling, and subsequently model these covariances jointly so as to find a shared low-dimensional subspace of the feature space. CONCLUSIONS: By viewing the projection of the latent true abundances onto this common structure, the variation is pared down to that which is shared among all datasets, and is likely to reflect more generalizable biological signal than can be inferred from individual datasets. We investigate several ways of achieving this, demonstrate that they work well on simulated and real metagenomic data in terms of signal retention and interpretability, and recommend a particular implementation.

Integrating phylogenetic and functional data in microbiome studies.

MOTIVATION: Microbiome functional data are frequently analyzed to identify associations between microbial functions (e.g. genes) and sample groups of interest. However, it is challenging to distinguish between different possible explanations for variation in community-wide functional profiles by considering functions alone. To help address this problem, we have developed POMS, a package that implements multiple phylogeny-aware frameworks to more robustly identify enriched functions. RESULTS: The key contribution is an extended balance-tree workflow that incorporates functional and taxonomic information to identify functions that are consistently enriched in sample groups across independent taxonomic lineages. Our package also includes a workflow for running phylogenetic regression. Based on simulated data we demonstrate that these approaches more accurately identify gene families that confer a selective advantage compared with commonly used tools. We also show that POMS in particular can identify enriched functions in real-world metagenomics datasets that are potential targets of strong selection on multiple members of the microbiome. AVAILABILITY AND IMPLEMENTATION: These workflows are freely available in the POMS R package at https://github.com/gavinmdouglas/POMS. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Exploring the microbiome of oral epithelial dysplasia as a predictor of malignant progression.

A growing body of research associates the oral microbiome and oral cancer. Well-characterized clinical samples with outcome data are required to establish relevant associations between the microbiota and disease. The objective of this study was to characterize the community variations and the functional implications of the microbiome in low-grade oral epithelial dysplasia (OED) using 16S rRNA gene sequencing from annotated archival swabs in progressing (P) and non-progressing (NP) OED. We characterised the microbial community in 90 OED samples - 30 swabs from low-grade OED that progressed to cancer (cases) and 60 swabs from low-grade OED that did not progress after a minimum of 5 years of follow up (matched control subjects). There were small but significant differences between P and NP samples in terms of alpha diversity as well as beta diversity in conjunction with other clinical factors such as age and smoking status for both taxa and functional predictions. Across all samples, the most abundant genus was Streptococcus, followed by Haemophilus, Rothia, and Neisseria. Taxa and predicted functions were identified that were significantly differentially abundant with progression status (all Ps and NPs), when samples were grouped broadly by the number of years between sampling and progression or in specific time to progression for Ps only. However, these differentially abundant features were typically present only at low abundances. For example, Campylobacter was present in slightly higher abundance in Ps (1.72%) than NPs (1.41%) and this difference was significant when Ps were grouped by time to progression. Furthermore, several of the significantly differentially abundant functions were linked to the Campylobacteraceae family in Ps and may justify further investigation. Larger cohort studies to further explore the microbiome as a potential biomarker of risk in OED are warranted.

The Phosphate Binder Ferric Citrate Alters the Gut Microbiome in Rats with Chronic Kidney Disease.

In chronic kidney disease (CKD), the gut microbiome is altered and bacterial-derived uremic toxins promote systemic inflammation and cardiovascular disease. Ferric citrate complex is a dietary phosphate binder prescribed for patients with end-stage kidney disease to treat hyperphosphatemia and secondary hyperparathyroidism. Iron is an essential nutrient in both microbes and mammals. This study was undertaken to test the hypothesis that the large iron load administered with ferric citrate in CKD may significantly change the gut microbiome. Male Sprague-Dawley rats underwent 5/6 nephrectomy to induce CKD. Normal control and CKD rats were randomized to regular chow or a 4% ferric citrate diet for 6 weeks. Fecal and cecal microbial DNA was analyzed via 16S ribosomal RNA gene sequencing on the Illumina MiSeq system. CKD rats had lower abundances of Firmicutes and Lactobacillus compared with normal rats and had lower overall gut microbial diversity. CKD rats treated with ferric citrate had improved hemoglobin and creatinine clearance and amelioration of hyperphosphatemia and hypertension. Ferric citrate treatment increased bacterial diversity in CKD rats almost to levels observed in control rats. The tryptophanase-possessing families Verrucomicrobia, Clostridiaceae, and Enterobacteriaceae were increased by ferric citrate treatment. The uremic toxins indoxyl sulfate and p-cresyl sulfate were not increased with ferric citrate treatment. Verrucomicrobia was largely represented by Akkermansia muciniphila, which has important roles in mucin degradation and gut barrier integrity. In summary, ferric citrate therapy in CKD rats was associated with significant changes in the gut microbiome and beneficial kidney and blood pressure parameters.

IslandViewer 3: more flexible, interactive genomic island discovery, visualization and analysis.

IslandViewer (http://pathogenomics.sfu.ca/islandviewer) is a widely used web-based resource for the prediction and analysis of genomic islands (GIs) in bacterial and archaeal genomes. GIs are clusters of genes of probable horizontal origin, and are of high interest since they disproportionately encode genes involved in medically and environmentally important adaptations, including antimicrobial resistance and virulence. We now report a major new release of IslandViewer, since the last release in 2013. IslandViewer 3 incorporates a completely new genome visualization tool, IslandPlot, enabling for the first time interactive genome analysis and gene search capabilities using synchronized circular, horizontal and vertical genome views. In addition, more curated virulence factors and antimicrobial resistance genes have been incorporated, and homologs of these genes identified in closely related genomes using strict filters. Pathogen-associated genes have been re-calculated for all pre-computed complete genomes. For user-uploaded genomes to be analysed, IslandViewer 3 can also now handle incomplete genomes, with an improved queuing system on compute nodes to handle user demand. Overall, IslandViewer 3 represents a significant new version of this GI analysis software, with features that may make it more broadly useful for general microbial genome analysis and visualization.

GenomeD3Plot: a library for rich, interactive visualizations of genomic data in web applications.

MOTIVATION: A simple static image of genomes and associated metadata is very limiting, as researchers expect rich, interactive tools similar to the web applications found in the post-Web 2.0 world. GenomeD3Plot is a light weight visualization library written in javascript using the D3 library. GenomeD3Plot provides a rich API to allow the rapid visualization of complex genomic data using a convenient standards based JSON configuration file. When integrated into existing web services GenomeD3Plot allows researchers to interact with data, dynamically alter the view, or even resize or reposition the visualization in their browser window. In addition GenomeD3Plot has built in functionality to export any resulting genome visualization in PNG or SVG format for easy inclusion in manuscripts or presentations. RESULTS: GenomeD3Plot is being utilized in the recently released Islandviewer 3 (www.pathogenomics.sfu.ca/islandviewer/) to visualize predicted genomic islands with other genome annotation data. However, its features enable it to be more widely applicable for dynamic visualization of genomic data in general. AVAILABILITY AND IMPLEMENTATION: GenomeD3Plot is licensed under the GNU-GPL v3 at https://github.com/brinkmanlab/GenomeD3Plot/. CONTACT: brinkman@sfu.ca.

Genomic Analysis of a Serotype 5 Streptococcus pneumoniae Outbreak in British Columbia, Canada, 2005-2009.

Background. Streptococcus pneumoniae can cause a wide spectrum of disease, including invasive pneumococcal disease (IPD). From 2005 to 2009 an outbreak of IPD occurred in Western Canada, caused by a S. pneumoniae strain with multilocus sequence type (MLST) 289 and serotype 5. We sought to investigate the incidence of IPD due to this S. pneumoniae strain and to characterize the outbreak in British Columbia using whole-genome sequencing. Methods. IPD was defined according to Public Health Agency of Canada guidelines. Two isolates representing the beginning and end of the outbreak were whole-genome sequenced. The sequences were analyzed for single nucleotide variants (SNVs) and putative genomic islands. Results. The peak of the outbreak in British Columbia was in 2006, when 57% of invasive S. pneumoniae isolates were serotype 5. Comparison of two whole-genome sequenced strains showed only 10 SNVs between them. A 15.5 kb genomic island was identified in outbreak strains, allowing the design of a PCR assay to track the spread of the outbreak strain. Discussion. We show that the serotype 5 MLST 289 strain contains a distinguishing genomic island, which remained genetically consistent over time. Whole-genome sequencing holds great promise for real-time characterization of outbreaks in the future and may allow responses tailored to characteristics identified in the genome.

IslandViewer update: Improved genomic island discovery and visualization.

IslandViewer (http://pathogenomics.sfu.ca/islandviewer) is a web-accessible application for the computational prediction and analysis of genomic islands (GIs) in bacterial and archaeal genomes. GIs are clusters of genes of probable horizontal origin and are of high interest because they disproportionately encode virulence factors and other adaptations of medical, environmental and industrial interest. Many computational tools exist for the prediction of GIs, but three of the most accurate methods are available in integrated form via IslandViewer: IslandPath-DIMOB, SIGI-HMM and IslandPick. IslandViewer GI predictions are precomputed for all complete microbial genomes from National Center for Biotechnology Information, with an option to upload other genomes and/or perform customized analyses using different settings. Here, we report recent changes to the IslandViewer framework that have vastly improved its efficiency in handling an increasing number of users, plus better facilitate custom genome analyses. Users may also now overlay additional annotations such as virulence factors, antibiotic resistance genes and pathogen-associated genes on top of current GI predictions. Comparisons of GIs between user-selected genomes are now facilitated through a highly requested side-by-side viewer. IslandViewer improvements aim to provide a more flexible interface, coupled with additional highly relevant annotation information, to aid analysis of GIs in diverse microbial species.

Oral microbial signatures associated with age and frailty in Canadian adults.

This study aimed to assess the association between the oral microbiome, age, and frailty. Data and saliva samples were obtained from male and female participants aged 35-70 years (n = 1357). Saliva samples were analysed by 16S rRNA gene sequencing and differences in microbial diversity and community compositions were examined in relation to chronological age and the frailty index (FI). Most alpha diversity measures (Richness, Shannon Diversity, Faith's Phylogenetic Diversity) showed an inverse association with frailty, whereas a positive association was observed with age and Shannon Diversity and Evenness. A further sex-stratified analysis revealed differences in measures of microbial diversity and composition. Multiple genera were detected as significantly differentially abundant with increasing frailty and age by at least two methods. With age, the relative abundance of Veillonella was reduced in both males and females, whereas increases in Corynebacterium appeared specific to males and Aggregatibacter, Fusobacterium, Neisseria, Stomatobaculum, and Porphyromonas specific to females. Beta diversity was significantly associated with multiple mental health components of the FI. This study shows age and frailty are differentially associated with measures of microbial diversity and composition, suggesting the oral microbiome may be a useful indicator of increased risk of frailty or a potential target for improving health in ageing adults.

Nitrogen fixation in the widely distributed marine γ-proteobacterial diazotroph Candidatus Thalassolituus haligoni.

The high diversity and global distribution of heterotrophic bacterial diazotrophs (HBDs) in the ocean has recently become apparent. However, understanding the role these largely uncultured microorganisms play in marine N2 fixation poses a challenge due to their undefined growth requirements and the complex regulation of the nitrogenase enzyme. We isolated and characterized Candidatus Thalassolituus haligoni, a member of a widely distributed clade of HBD belonging to the Oceanospirillales. Analysis of its nifH gene via amplicon sequencing revealed the extensive distribution of Cand. T. haligoni across the Pacific, Atlantic, and Arctic Oceans. Pangenome analysis indicates that the isolate shares >99% identity with an uncultured metagenome-assembled genome called Arc-Gamma-03, recently recovered from the Arctic Ocean. Through combined genomic, proteomic, and physiological approaches, we confirmed that the isolate fixes N2 gas. However, the mechanisms governing nitrogenase regulation in Cand. T. haligoni remain unclear. We propose Cand. T. haligoni as a globally distributed, cultured HBD model species within this understudied clade of Oceanospirillales.

MicrobeDB: a locally maintainable database of microbial genomic sequences.

UNLABELLED: Analysis of microbial genomes often requires the general organization and comparison of tens to thousands of genomes both from public repositories and unpublished sources. MicrobeDB provides a foundation for such projects by the automation of downloading published, completed bacterial and archaeal genomes from key sources, parsing annotations of all genomes (both public and private) into a local database, and allowing interaction with the database through an easy to use programming interface. MicrobeDB creates a simple to use, easy to maintain, centralized local resource for various large-scale comparative genomic analyses and a back-end for future microbial application design. AVAILABILITY: MicrobeDB is freely available under the GNU-GPL at: http://github.com/mlangill/microbedb/

Investigating the oral microbiome in retrospective and prospective cases of prostate, colon, and breast cancer.

The human microbiome has been proposed as a potentially useful biomarker for several cancers. To examine this, we made use of salivary samples from the Atlantic Partnership for Tomorrow's Health (PATH) project and Alberta's Tomorrow Project (ATP). Sample selection was divided into both a retrospective and prospective case control design examining prostate, breast, and colon cancer. In total 89 retrospective and 260 prospective cancer cases were matched to non-cancer controls and saliva samples were sequenced using 16S rRNA gene sequencing. We found no significant differences in alpha diversity. All beta diversity measures were insignificant except for unweighted UniFrac profiles in retrospective breast cancer cases and weighted UniFrac, Bray-Curtis and Robust Atchinson's distances in colon cancer after testing with age and sex adjusted MiRKAT models. Differential abundance (DA) analysis showed several taxa that were associated with previous cancer in all three groupings. Only one genus (Clostridia UCG-014) in breast cancer and one ASV (Fusobacterium periodonticum) in colon cancer was identified by more than one DA tool. In prospective cases three ASVs were associated with colon cancer, one ASV with breast cancer, and one ASV with prostate cancer. Random Forest classification showed low levels of signal in both study designs in breast and prostate cancer. Contrastingly, colon cancer did show signal in our retrospective analysis (AUC: 0.737) and in one of two prospective cohorts (AUC: 0.717). Our results indicate that it is unlikely that reliable microbial oral biomarkers for breast and prostate cancer exist.. However, further research into the oral microbiome and colon cancer could be fruitful.

Changes to Soil Microbiome Resulting from Synergetic Effects of Fungistatic Compounds Pyrimethanil and Fluopyram in Lowbush Blueberry Agriculture, with Nine Fungicide Products Tested.

Lowbush blueberries (Vaccinium spp.) are a crop of economic significance to Atlantic Canada, Quebec, and Maine. The fruit is produced by the management of naturally occurring plant populations. The plants have an intimate relationship with the soil microbiome and depend on it for their health and productivity. Fungicides are an important tool in combatting disease pressure but pose a potential risk to soil health. In this study, amplicon sequencing was used to determine the effects of six fungistatic compounds both alone and in combination via nine commercially available fungicide products on the bacterial and fungal microbiomes associated with lowbush blueberries and to study whether these effects are reflected in crop outcomes and plant phenotypes. One fungicide, Luna Tranquility, a combination of fluopyram and pyrimethanil, was found to impart significant effects to fungal and bacterial community structure, fungal taxonomic abundances, and bacterial functions relative to control. The two fungicides which contained fluopyram and pyrimethanil as single ingredients (Velum Prime and Scala, respectively) did not induce significant changes in any of these regards. These results suggest the possibility that these microbiome changes are the result of the synergistic effect of fluopyram and pyrimethanil on soil microbiomes. While these results suggest a significant disruption to the soil microbiome, no corresponding changes to crop development and outcomes were noted. Ultimately, the majority of the fungicides analysed in this trial did not produce significant changes to the soil microbiome relative to the untreated group (UTG). However, one of the fungicide treatments, Luna Tranquility, did produce significant changes to the soil ecosystem that could have longer-term effects on soil health and its future use may merit additional investigation onto its ecotoxicological properties.

From defaults to databases: parameter and database choice dramatically impact the performance of metagenomic taxonomic classification tools.

In metagenomic analyses of microbiomes, one of the first steps is usually the taxonomic classification of reads by comparison to a database of previously taxonomically classified genomes. While different studies comparing metagenomic taxonomic classification methods have determined that different tools are 'best', there are two tools that have been used the most to-date: Kraken (k-mer-based classification against a user-constructed database) and MetaPhlAn (classification by alignment to clade-specific marker genes), the latest versions of which are Kraken2 and MetaPhlAn 3, respectively. We found large discrepancies in both the proportion of reads that were classified as well as the number of species that were identified when we used both Kraken2 and MetaPhlAn 3 to classify reads within metagenomes from human-associated or environmental datasets. We then investigated which of these tools would give classifications closest to the real composition of metagenomic samples using a range of simulated and mock samples and examined the combined impact of tool-parameter-database choice on the taxonomic classifications given. This revealed that there may not be a one-size-fits-all 'best' choice. While Kraken2 can achieve better overall performance, with higher precision, recall and F1 scores, as well as alpha- and beta-diversity measures closer to the known composition than MetaPhlAn 3, the computational resources required for this may be prohibitive for many researchers, and the default database and parameters should not be used. We therefore conclude that the best tool-parameter-database choice for a particular application depends on the scientific question of interest, which performance metric is most important for this question and the limit of available computational resources.

Hepatotoxicity of polymeric proanthocyanidins is caused by translocation of bacterial lipopolysaccharides through impaired gut epithelium.

Polymeric proanthocyanidins (P-PAC) induced hepatotoxicity in C57BL/6 mice. Mice were supplemented with P-PAC alone or with a mixture of probiotic bacteria (PB), Lactobacillus, Bifidobacterium, and Akkermansia muciniphila for 14 consecutive days. The liver tissues of sacrificed mice were analyzed by mass spectrometry to identify and quantify the P-PAC metabolites. Potential P-PAC metabolites, 2-hydroxyphenylacetic acid and pyrocatechol were detected in higher concentrations and 4-hydroxybenzoic acid was detected exclusively in the mice supplemented with P-PAC and PB. Supplementation with P-PAC alone or with PB caused no shift in the α-diversity of mice gut microbiota. P-PAC induced nonalcoholic steatohepatitis in mice through increasing liver exposure to intestinal bacterial lipopolysaccharides by reducing expression of gut epithelial tight junction proteins, claudin-3 and occludin. Lipopolysaccharide concentrations in the livers of mice supplemented with P-PAC were significantly high compared to the control mice. Furthermore, P-PAC downregulated the expressions of claudin-3 and claudin-4 tight junction proteins in cultured Caco-2 cell monolayers. PB biotransformed P-PAC into bioavailable metabolites and potentially reduced the toxicity of P-PAC. The toxicity of P-PAC and their synbiotics need to be critically evaluated for the safety of human consumption.

The gastrointestinal antibiotic resistome in pediatric leukemia and lymphoma patients.

Introduction: Most children with leukemia and lymphoma experience febrile neutropenia. These are treated with empiric antibiotics that include ß-lactams and/or vancomycin. These are often administered for extended periods, and the effect on the resistome is unknown. Methods: We examined the impact of repeated courses and duration of antibiotic use on the resistome of 39 pediatric leukemia and lymphoma patients. Shotgun metagenome sequences from 127 stool samples of pediatric oncology patients were examined for abundance of antibiotic resistance genes (ARGs) in each sample. Abundances were grouped by repeated courses (no antibiotics, 1-2 courses, 3+ courses) and duration (no use, short duration, long and/or mixed durationg) of ß-lactams, vancomycin and "any antibiotic" use. We assessed changes in both taxonomic composition and prevalence of ARGs among these groups. Results: We found that Bacteroidetes taxa and ß-lactam resistance genes decreased, while opportunistic Firmicutes and Proteobacteria taxa, along with multidrug resistance genes, increased with repeated courses and/or duration of antibiotics. Efflux pump related genes predominated (92%) among the increased multidrug genes. While we found ß-lactam ARGs present in the resistome, the taxa that appear to contain them were kept in check by antibiotic treatment. Multidrug ARGs, mostly efflux pumps or regulators of efflux pump genes, were associated with opportunistic pathogens, and both increased in the resistome with repeated antibiotic use and/or increased duration. Conclusions: Given the strong association between opportunistic pathogens and multidrug-related efflux pumps, we suggest that drug efflux capacity might allow the opportunistic pathogens to persist or increase despite repeated courses and/or duration of antibiotics. While drug efflux is the most direct explanation, other mechanisms that enhance the ability of opportunistic pathogens to handle environmental stress, or other aspects of the treatment environment, could also contribute to their ability to flourish within the gut during treatment. Persistence of opportunistic pathogens in an already dysbiotic and weakened gastrointestinal tract could increase the likelihood of life-threatening blood borne infections. Of the 39 patients, 59% experienced at least one gastrointestinal or blood infection and 60% of bacteremia's were bacteria found in stool samples. Antimicrobial stewardship and appropriate use and duration of antibiotics could help reduce morbidity and mortality in this vulnerable population.

Microbial hitchhikers harbouring antimicrobial-resistance genes in the riverine plastisphere.

BACKGROUND: The widespread nature of plastic pollution has given rise to wide scientific and social concern regarding the capacity of these materials to serve as vectors for pathogenic bacteria and reservoirs for Antimicrobial Resistance Genes (ARG). In- and ex-situ incubations were used to characterise the riverine plastisphere taxonomically and functionally in order to determine whether antibiotics within the water influenced the ARG profiles in these microbiomes and how these compared to those on natural surfaces such as wood and their planktonic counterparts. RESULTS: We show that plastics support a taxonomically distinct microbiome containing potential pathogens and ARGs. While the plastisphere was similar to those biofilms that grew on wood, they were distinct from the surrounding water microbiome. Hence, whilst potential opportunistic pathogens (i.e. Pseudomonas aeruginosa, Acinetobacter and Aeromonas) and ARG subtypes (i.e. those that confer resistance to macrolides/lincosamides, rifamycin, sulfonamides, disinfecting agents and glycopeptides) were predominant in all surface-related microbiomes, especially on weathered plastics, a completely different set of potential pathogens (i.e. Escherichia, Salmonella, Klebsiella and Streptococcus) and ARGs (i.e. aminoglycosides, tetracycline, aminocoumarin, fluoroquinolones, nitroimidazole, oxazolidinone and fosfomycin) dominated in the planktonic compartment. Our genome-centric analysis allowed the assembly of 215 Metagenome Assembled Genomes (MAGs), linking ARGs and other virulence-related genes to their host. Interestingly, a MAG belonging to Escherichia -that clearly predominated in water- harboured more ARGs and virulence factors than any other MAG, emphasising the potential virulent nature of these pathogenic-related groups. Finally, ex-situ incubations using environmentally-relevant concentrations of antibiotics increased the prevalence of their corresponding ARGs, but different riverine compartments -including plastispheres- were affected differently by each antibiotic. CONCLUSIONS: Our results provide insights into the capacity of the riverine plastisphere to harbour a distinct set of potentially pathogenic bacteria and function as a reservoir of ARGs. The environmental impact that plastics pose if they act as a reservoir for either pathogenic bacteria or ARGs is aggravated by the persistence of plastics in the environment due to their recalcitrance and buoyancy. Nevertheless, the high similarities with microbiomes growing on natural co-occurring materials and even more worrisome microbiome observed in the surrounding water highlights the urgent need to integrate the analysis of all environmental compartments when assessing risks and exposure to pathogens and ARGs in anthropogenically-impacted ecosystems. Video Abstract.

Sifting through genomes with iterative-sequence clustering produces a large, phylogenetically diverse protein-family resource.

BACKGROUND: New computational resources are needed to manage the increasing volume of biological data from genome sequencing projects. One fundamental challenge is the ability to maintain a complete and current catalog of protein diversity. We developed a new approach for the identification of protein families that focuses on the rapid discovery of homologous protein sequences. RESULTS: We implemented fully automated and high-throughput procedures to de novo cluster proteins into families based upon global alignment similarity. Our approach employs an iterative clustering strategy in which homologs of known families are sifted out of the search for new families. The resulting reduction in computational complexity enables us to rapidly identify novel protein families found in new genomes and to perform efficient, automated updates that keep pace with genome sequencing. We refer to protein families identified through this approach as "Sifting Families," or SFams. Our analysis of ~10.5 million protein sequences from 2,928 genomes identified 436,360 SFams, many of which are not represented in other protein family databases. We validated the quality of SFam clustering through statistical as well as network topology-based analyses. CONCLUSIONS: We describe the rapid identification of SFams and demonstrate how they can be used to annotate genomes and metagenomes. The SFam database catalogs protein-family quality metrics, multiple sequence alignments, hidden Markov models, and phylogenetic trees. Our source code and database are publicly available and will be subject to frequent updates (http://edhar.genomecenter.ucdavis.edu/sifting_families/).

Antibiotic and antifungal use in pediatric leukemia and lymphoma patients are associated with increasing opportunistic pathogens and decreasing bacteria responsible for activities that enhance colonic defense.

Due to decreased immunity, both antibiotics and antifungals are regularly used in pediatric hematologic-cancer patients as a means to prevent severe infections and febrile neutropenia. The general effect of antibiotics on the human gut microbiome is profound, yielding decreased diversity and changes in community structure. However, the specific effect on pediatric oncology patients is not well-studied. The effect of antifungal use is even less understood, having been studied only in mouse models. Because the composition of the gut microbiome is associated with regulation of hematopoiesis, immune function and gastrointestinal integrity, changes within the patient gut can have implications for the clinical management of hematologic malignancies. The pediatric population is particularly challenging because the composition of the microbiome is age dependent, with some of the most pronounced changes occurring in the first three years of life. We investigated how antibiotic and antifungal use shapes the taxonomic composition of the stool microbiome in pediatric patients with leukemia and lymphoma, as inferred from both 16S rRNA and metagenome data. Associations with age, antibiotic use and antifungal use were investigated using multiple analysis methods. In addition, multivariable differential abundance was used to identify and assess specific taxa that were associated with multiple variables. Both antibiotics and antifungals were linked to a general decline in diversity in stool samples, which included a decrease in relative abundance in butyrate producers that play a critical role in host gut physiology (e.g., Faecalibacterium, Anaerostipes, Dorea, Blautia),. Furthermore, antifungal use was associated with a significant increase in relative abundance of opportunistic pathogens. Collectively, these findings have important implications for the treatment of leukemia and lymphoma patients. Butyrate is important for gastrointestinal integrity; it inhibits inflammation, reinforces colonic defense, mucosal immunity. and decreases oxidative stress. The routine use of broad-spectrum anti-infectives in pediatric oncology patients could simultaneously contribute to a decline in gastrointestinal integrity and colonic defense while promoting increases in opportunistic pathogens within the patient gut. Because the gut microbiome has been linked to both short-term clinical outcomes, and longer-lasting health effects, systematic characterization of the gut microbiome in pediatric patients during, and beyond, treatment is warranted.