Most of our knowledge about the ancient evolutionary history of organisms has been derived from data associated with specific known organisms (i.e., organisms that we can study directly such as plants, metazoans, and culturable microbes). Recently, however, a new source of data for such studies has arrived: DNA sequence data generated directly from environmental samples. Such metagenomic data has enormous potential in a variety of areas including, as we argue here, in studies of very early events in the evolution of gene families and of species.

Dental caries, one of the most globally widespread infectious diseases, is intimately linked to pH dynamics. In supragingival plaque, after the addition of a carbohydrate source, bacterial metabolism decreases the pH which then subsequently recovers. Molecular mechanisms supporting this important homeostasis are poorly characterized in part due to the fact that there are hundreds of active species in dental plaque. Only a few mechanisms (for example, lactate fermentation, the arginine deiminase system) have been identified and studied in detail. Here, we conducted what is to our knowledge, the first full transcriptome and metabolome analysis of a diverse oral plaque community by using a functionally and taxonomically robust in vitro model system greater than 100 species. Differential gene expression analyses from the complete transcriptome of 14 key community members revealed highly varied regulation of both known and previously unassociated pH-neutralizing pathways as a response to the pH drop. Unique expression and metabolite signatures from 400 detected metabolites were found for each stage along the pH curve suggesting it may be possible to define healthy and diseased states of activity. Importantly, for the maintenance of healthy plaque pH, gene transcription activity of known and previously unrecognized pH-neutralizing pathways was associated with the genera Lactobacillus, Veillonella and Streptococcus during the pH recovery phase. Our in vitro study provides a baseline for defining healthy and disease-like states and highlights the power of moving beyond single and dual species applications to capture key players and their orchestrated metabolic activities within a complex human oral microbiome model.

Colonic polyps are common tumors occurring in ~50% of Western populations with ~10% risk of malignant progression. Dietary agents have been considered the primary environmental exposure to promote colorectal cancer (CRC) development. However, the colonic mucosa is permanently in contact with the microbiota and its metabolic products including toxins that also have the potential to trigger oncogenic transformation.

Antibiotics are a mainstay of treatment for bacterial infections worldwide, yet the effects of typical antibiotic prescriptions on human indigenous microbiota have not been thoroughly evaluated. We examined the effects of the two most commonly prescribed antibiotics (amoxicillin and azithromycin) in the USA to discern whether short-term antibiotic courses may have prolonged effects on human microbiota.

Analyses of metagenome data (MG) and metatranscriptome data (MT) are often challenged by a paucity of complete reference genome sequences and the uneven/low sequencing depth of the constituent organisms in the microbial community, which respectively limit the power of reference-based alignment and de novo sequence assembly. These limitations make accurate protein family classification and abundance estimation challenging, which in turn hamper downstream analyses such as abundance profiling of metabolic pathways, identification of differentially encoded/expressed genes, and de novo reconstruction of complete gene and protein sequences from the protein family of interest. The profile hidden Markov model (HMM) framework enables the construction of very useful probabilistic models for protein families that allow for accurate modeling of position specific matches, insertions, and deletions. We present a novel homology detection algorithm that integrates banded Viterbi algorithm for profile HMM parsing with an iterative simultaneous alignment and assembly computational framework. The algorithm searches a given profile HMM of a protein family against a database of fragmentary MG/MT sequencing data and simultaneously assembles complete or near-complete gene and protein sequences of the protein family. The resulting program, HMM-GRASPx, demonstrates superior performance in aligning and assembling homologs when benchmarked on both simulated marine MG and real human saliva MG datasets. On real supragingival plaque and stool MG datasets that were generated from healthy individuals, HMM-GRASPx accurately estimates the abundances of the antimicrobial resistance (AMR) gene families and enables accurate characterization of the resistome profiles of these microbial communities. For real human oral microbiome MT datasets, using the HMM-GRASPx estimated transcript abundances significantly improves detection of differentially expressed (DE) genes. Finally, HMM-GRASPx was used to reconstruct comprehensive sets of complete or near-complete protein and nucleotide sequences for the query protein families. HMM-GRASPx is freely available online from http://sourceforge.net/projects/hmm-graspx.

Analysis of human body microbial diversity is fundamental to understanding community structure, biology and ecology. The National Institutes of Health Human Microbiome Project (HMP) has provided an unprecedented opportunity to examine microbial diversity within and across body habitats and individuals through pyrosequencing-based profiling of 16 S rRNA gene sequences (16 S) from habits of the oral, skin, distal gut, and vaginal body regions from over 200 healthy individuals enabling the application of statistical techniques. In this study, two approaches were applied to elucidate the nature and extent of human microbiome diversity. First, bootstrap and parametric curve fitting techniques were evaluated to estimate the maximum number of unique taxa, S(max), and taxa discovery rate for habitats across individuals. Next, our results demonstrated that the variation of diversity within low abundant taxa across habitats and individuals was not sufficiently quantified with standard ecological diversity indices. This impact from low abundant taxa motivated us to introduce a novel rank-based diversity measure, the Tail statistic, ("τ"), based on the standard deviation of the rank abundance curve if made symmetric by reflection around the most abundant taxon. Due to τ's greater sensitivity to low abundant taxa, its application to diversity estimation of taxonomic units using taxonomic dependent and independent methods revealed a greater range of values recovered between individuals versus body habitats, and different patterns of diversity within habitats. The greatest range of τ values within and across individuals was found in stool, which also exhibited the most undiscovered taxa. Oral and skin habitats revealed variable diversity patterns, while vaginal habitats were consistently the least diverse. Collectively, these results demonstrate the importance, and motivate the introduction, of several visualization and analysis methods tuned specifically for next-generation sequence data, further revealing that low abundant taxa serve as an important reservoir of genetic diversity in the human microbiome.

Protein complexes are typically analyzed by affinity purification and subsequent mass spectrometric analysis. However, in most cases the structure and topology of the complexes remains elusive from such studies. Here we investigate how the yeast two-hybrid system can be used to analyze direct interactions among proteins in a complex. First we tested all pairwise interactions among the seven proteins of Escherichia coli DNA polymerase III as well as an uncharacterized complex that includes MntR and PerR. Four and seven interactions were identified in these two complexes, respectively. In addition, we review Y2H data for three other complexes of known structure which serve as "gold-standards", namely Varicella Zoster Virus (VZV) ribonucleotide reductase (RNR), the yeast proteasome, and bacteriophage lambda. Finally, we review an Y2H analysis of the human spliceosome which may serve as an example for a dynamic mega-complex.

The YbeB (DUF143) family of uncharacterized proteins is encoded by almost all bacterial and eukaryotic genomes but not archaea. While they have been shown to be associated with ribosomes, their molecular function remains unclear. Here we show that YbeB is a ribosomal silencing factor (RsfA) in the stationary growth phase and during the transition from rich to poor media. A knock-out of the rsfA gene shows two strong phenotypes: (i) the viability of the mutant cells are sharply impaired during stationary phase (as shown by viability competition assays), and (ii) during transition from rich to poor media the mutant cells adapt slowly and show a growth block of more than 10 hours (as shown by growth competition assays). RsfA silences translation by binding to the L14 protein of the large ribosomal subunit and, as a consequence, impairs subunit joining (as shown by molecular modeling, reporter gene analysis, in vitro translation assays, and sucrose gradient analysis). This particular interaction is conserved in all species tested, including Escherichia coli, Treponema pallidum, Streptococcus pneumoniae, Synechocystis PCC 6803, as well as human mitochondria and maize chloroplasts (as demonstrated by yeast two-hybrid tests, pull-downs, and mutagenesis). RsfA is unrelated to the eukaryotic ribosomal anti-association/60S-assembly factor eIF6, which also binds to L14, and is the first such factor in bacteria and organelles. RsfA helps cells to adapt to slow-growth/stationary phase conditions by down-regulating protein synthesis, one of the most energy-consuming processes in both bacterial and eukaryotic cells.

Ureaplasma urealyticum (UUR) and Ureaplasma parvum (UPA) are sexually transmitted bacteria among humans implicated in a variety of disease states including but not limited to: nongonococcal urethritis, infertility, adverse pregnancy outcomes, chorioamnionitis, and bronchopulmonary dysplasia in neonates. There are 10 distinct serotypes of UUR and 4 of UPA. Efforts to determine whether difference in pathogenic potential exists at the ureaplasma serovar level have been hampered by limitations of antibody-based typing methods, multiple cross-reactions and poor discriminating capacity in clinical samples containing two or more serovars.

Viruses are the most abundant biological entities on our planet. Interactions between viruses and their hosts impact several important biological processes in the world's oceans such as horizontal gene transfer, microbial diversity and biogeochemical cycling. Interrogation of microbial metagenomic sequence data collected as part of the Sorcerer II Global Ocean Expedition (GOS) revealed a high abundance of viral sequences, representing approximately 3% of the total predicted proteins. Cluster analyses of the viral sequences revealed hundreds to thousands of viral genes encoding various metabolic and cellular functions. Quantitative analyses of viral genes of host origin performed on the viral fraction of aquatic samples confirmed the viral nature of these sequences and suggested that significant portions of aquatic viral communities behave as reservoirs of such genetic material. Distributional and phylogenetic analyses of these host-derived viral sequences also suggested that viral acquisition of environmentally relevant genes of host origin is a more abundant and widespread phenomenon than previously appreciated. The predominant viral sequences identified within microbial fractions originated from tailed bacteriophages and exhibited varying global distributions according to viral family. Recruitment of GOS viral sequence fragments against 27 complete aquatic viral genomes revealed that only one reference bacteriophage genome was highly abundant and was closely related, but not identical, to the cyanomyovirus P-SSM4. The co-distribution across all sampling sites of P-SSM4-like sequences with the dominant ecotype of its host, Prochlorococcus supports the classification of the viral sequences as P-SSM4-like and suggests that this virus may influence the abundance, distribution and diversity of one of the most dominant components of picophytoplankton in oligotrophic oceans. In summary, the abundance and broad geographical distribution of viral sequences within microbial fractions, the prevalence of genes among viral sequences that encode microbial physiological function and their distinct phylogenetic distribution lend strong support to the notion that viral-mediated gene acquisition is a common and ongoing mechanism for generating microbial diversity in the marine environment.

Respiratory syncytial virus (RSV) is an established cause of serious lower respiratory disease in infants, elderly and high-risk populations. The OUTSMART surveillance program aims to characterize patient populations and currently circulating RSV strains, and monitor temporal and geographic evolution of RSV F and G proteins in the U.S.

Necrotizing enterocolitis (NEC) is the most common surgical emergency in preterm infants, and pathogenesis associates with changes in the fecal microbiome. As fecal samples incompletely represent microbial communities in intestinal mucosa, we sought to determine the NEC tissue-specific microbiome and assess its contribution to pathogenesis.

A crucial task in metagenomic analysis is to annotate the function and taxonomy of the sequencing reads generated from a microbiome sample. In general, the reads can either be assembled into contigs and searched against reference databases, or individually searched without assembly. The first approach may suffer from fragmentary and incomplete assembly, while the second is hampered by the reduced functional signal contained in the short reads. To tackle these issues, we have previously developed GRASP (Guided Reference-based Assembly of Short Peptides), which accepts a reference protein sequence as input and aims to assemble its homologs from a database containing fragmentary protein sequences. In addition to a gene-centric assembly tool, GRASP also serves as a homolog search tool when using the assembled protein sequences as templates to recruit reads. GRASP has significantly improved recall rate (60-80% vs. 30-40%) compared to other homolog search tools such as BLAST. However, GRASP is both time- and space-consuming. Subsequently, we developed GRASPx, which is 30X faster than GRASP. Here, we present a completely redesigned algorithm, GRASP2, for this computational problem.

While insulin replacement therapy restores the health and prevents the onset of diabetic complications (DC) for many decades, some T1D patients have elevated hemoglobin A1c values suggesting poor glycemic control, a risk factor of DC. We surveyed the stool microbiome and urinary proteome of a cohort of 220 adolescents and children, half of which had lived with T1D for an average of 7 years and half of which were healthy siblings. Phylogenetic analysis of the 16S rRNA gene did not reveal significant differences in gut microbial alpha-diversity comparing the two cohorts. The urinary proteome of T1D patients revealed increased abundances of several lysosomal proteins that correlated with elevated HbA1c values. In silico protein network analysis linked such proteins to extracellular matrix components and the glycoprotein LRG1. LRG1 is a prominent inflammation and neovascularization biomarker. We hypothesize that these changes implicate aberrant glycation of macromolecules that alter lysosomal function and metabolism in renal tubular epithelial cells, cells that line part of the upper urinary tract.

Effective treatment of Mycobacterium tuberculosis (Mtb) infection requires at least 6 months of daily therapy with multiple orally administered antibiotics. Although this drug regimen is administered annually to millions worldwide, the impact of such intensive antimicrobial treatment on the host microbiome has never been formally investigated. Here, we characterized the longitudinal outcome of conventional isoniazid-rifampin-pyrazinamide (HRZ) TB drug administration on the diversity and composition of the intestinal microbiota in Mtb-infected mice by means of 16S rRNA sequencing. We also investigated the effects of each of the individual antibiotics alone and in different combinations.

Despite being commonly used to collect upper airway epithelial lining fluid, nasal washes are poorly reproducible, not suitable for serial sampling, and limited by a dilution effect. In contrast, nasal filters lack these limitations and are an attractive alternative. To examine whether nasal filters are superior to nasal washes as a sampling method for the characterization of the upper airway microbiome and immune response, we collected paired nasal filters and washes from a group of 40 healthy children and adults. To characterize the upper airway microbiome, we used 16S ribosomal RNA and shotgun metagenomic sequencing. To characterize the immune response, we measured total protein using a BCA assay and 53 immune mediators using multiplex magnetic bead-based assays. We conducted statistical analyses to compare common microbial ecology indices and immune-mediator median fluorescence intensities (MFIs) between sample types. In general, nasal filters were more likely to pass quality control in both children and adults. There were no significant differences in microbiome community richness, α-diversity, or structure between pediatric samples types; however, these were all highly dissimilar between adult sample types. In addition, there were significant differences in the abundance of amplicon sequence variants between sample types in children and adults. In adults, total proteins were significantly higher in nasal filters than nasal washes; consequently, the immune-mediator MFIs were not well detected in nasal washes. Based on better quality control sequencing metrics and higher immunoassay sensitivity, our results suggest that nasal filters are a superior sampling method to characterize the upper airway microbiome and immune response in both children and adults.

Staphylococcus aureus nasal carriage provides the bacterial reservoir for opportunistic infection. In comparing the nasal microbiomes of culture-defined persistent S. aureus carriers versus noncarriers, we detected S. aureus DNA in all noses, including those with an established history of S. aureus negativity based on culture. Colonization with Gammaproteobacteria, including Klebsiella aerogenes, Citrobacter koseri, Moraxella lincolnii, and select Acinetobacter spp., was associated with S. aureus noncarriage. We next developed physiological competition assays for testing anti-S. aureus activity of isolated nasal species, utilizing medium modeling the nutrient-limited fluid of the nasal mucosa, polarized primary nasal epithelia, and nasal secretions. K. aerogenes from the nose of an S. aureus noncarrier demonstrated >99% inhibition of S. aureus recovery in all assays, even when S. aureus was coincubated in 9-fold excess. Secreted S. aureus inhibitory proteins from K. aerogenes and M. lincolnii were heat-stable and <30 kDa, fitting the profile of antimicrobial peptides. C. koseri, Acinetobacter haemolyticus, Acinetobacter junii, and Acinetobacter schindleri inhibited S. aureus recovery on nasal epithelia in a contact-dependent manner, while several other species either had no effect or promoted S. aureus growth. Collectively, this project is one of the first to identify resident nasal microbial species that impede S. aureus survival, and it implies that detectable nasal S. aureus results from shifts in microbial community composition.IMPORTANCE Nasal carriage of Staphylococcus aureus is a risk factor for infection, but it is not yet understood why some individuals carry nasal S. aureus persistently, intermittently, or seemingly not at all when tested via culture methods. This study compared the nasal microbiomes of established S. aureus carriers and noncarriers, identified species associated with noncarriage, and tested them for anti-S. aureus activity using assays developed to model the nutrient-limited nasal mucosa. We determined that all nostril swabs contain S. aureus DNA, even swabs from hosts considered to be long-term noncarriers. Select members of the Gammaproteobacteria class were more prevalent in noncarrier than carrier nostrils and demonstrated potent activity against multiple strains of S. aureus The results described here provide a better understanding of how the nasal microbiome controls S. aureus growth and viability and may be useful in the design of improved S. aureus decolonization strategies.

Clinical severity of Staphylococcus aureus respiratory infection correlates with alpha toxin (AT) expression. AT activates the NLRP3 inflammasome; deletion of Nlrp3, or AT neutralization, protects mice from lethal S. aureus pneumonia. We tested the hypothesis that this protection is not due to a reduction in inflammasome-dependent cytokines (IL-1β/IL-18) but increased bactericidal function of macrophages. In vivo, neutralization of AT or NLRP3 improved bacterial clearance and survival, while blocking IL-1β/IL-18 did not. Primary human monocytes were used in vitro to determine the mechanism through which NLRP3 alters bacterial killing. In cells treated with small interfering RNA (siRNA) targeting NLRP3 or infected with AT-null S. aureus, mitochondria co-localize with bacterial-containing phagosomes. Mitochondrial engagement activates caspase-1, a process dependent on complex II of the electron transport chain, near the phagosome, promoting its acidification. These data demonstrate a mechanism utilized by S. aureus to sequester itself from antimicrobial processes within the cell.

We present two standards developed by the Genomic Standards Consortium (GSC) for reporting bacterial and archaeal genome sequences. Both are extensions of the Minimum Information about Any (x) Sequence (MIxS). The standards are the Minimum Information about a Single Amplified Genome (MISAG) and the Minimum Information about a Metagenome-Assembled Genome (MIMAG), including, but not limited to, assembly quality, and estimates of genome completeness and contamination. These standards can be used in combination with other GSC checklists, including the Minimum Information about a Genome Sequence (MIGS), Minimum Information about a Metagenomic Sequence (MIMS), and Minimum Information about a Marker Gene Sequence (MIMARKS). Community-wide adoption of MISAG and MIMAG will facilitate more robust comparative genomic analyses of bacterial and archaeal diversity.

Despite the medical significance to humans and important ecological roles filled by vipers, few high-quality genomic resources exist for these snakes outside of a few genera of pitvipers. Here we sequence, assemble, and annotate the genome of Fea's Viper (Azemiops feae). This taxon is distributed in East Asia and belongs to a monotypic subfamily, sister to the pitvipers. The newly sequenced genome resulted in a 1.56 Gb assembly, a contig N50 of 1.59 Mb, with 97.6% of the genome assembly in contigs >50 Kb, and a BUSCO completeness of 92.4%. We found that A. feae venom is primarily composed of phospholipase A2 (PLA2) proteins expressed by genes that likely arose from lineage-specific PLA2 gene duplications. Additionally, we show that renin, an enzyme associated with blood pressure regulation in mammals and known from the venoms of two viper species including A. feae, is expressed in the venom gland at comparative levels to known toxins and is present in the venom proteome. The cooption of this gene as a toxin may be more widespread in viperids than currently known. To investigate the historical population demographics of A. feae, we performed coalescent-based analyses and determined that the effective population size has remained stable over the last 100 kyr. This suggests Quaternary glacial cycles likely had minimal influence on the demographic history of A. feae. This newly assembled genome will be an important resource for studying the genomic basis of phenotypic evolution and understanding the diversification of venom toxin gene families.