Microbial transformation from normal oral microbiota to acute endodontic infections.

BACKGROUND: Endodontic infections are a leading cause of oro-facial pain and tooth loss in western countries, and may lead to severe life-threatening infections. These infections are polymicrobial with high bacterial diversity. Understanding the spatial transition of microbiota from normal oral cavities through the infected root canal to the acute periapical abscess can improve our knowledge of the pathogenesis of endodontic infections and lead to more effective treatment. We obtained samples from the oral cavity, infected root canal and periapical abscess of 8 patients (5 with localized and 3 with systemic infections). Microbial populations in these samples were analyzed using next-generation sequencing of 16S rRNA amplicons. Bioinformatics tools and statistical tests with rigorous criteria were used to elucidate the spatial transition of the microbiota from normal to diseased sites. RESULTS: On average, 10,000 partial 16S rRNA gene sequences were obtained from each sample. All sequences fell into 11 different bacterial phyla. The microbial diversity in root canal and abscess samples was significantly lower than in the oral samples. Streptococcus was the most abundant genus in oral cavities while Prevotella and Fusobacterium were most abundant in diseased samples. The microbiota community structures of root canal and abscess samples were, however, more similar to each other than to the oral cavity microbiota. Using rigorous criteria and novel bioinformatics tools, we found that Granulicatella adiacens, Eubacterium yurii, Prevotella melaninogenica, Prevotella salivae, Streptococcus mitis, and Atopobium rimae were over-represented in diseased samples. CONCLUSIONS: We used a novel approach and high-throughput methodologies to characterize the microbiota associated normal and diseased oral sites in the same individuals.

Genomic diversity of 2010 Haitian cholera outbreak strains.

The millions of deaths from cholera during the past 200 y, coupled with the morbidity and mortality of cholera in Haiti since October 2010, are grim reminders that Vibrio cholerae, the etiologic agent of cholera, remains a scourge. We report the isolation of both V. cholerae O1 and non-O1/O139 early in the Haiti cholera epidemic from samples collected from victims in 18 towns across eight Arrondissements of Haiti. The results showed two distinct populations of V. cholerae coexisted in Haiti early in the epidemic. As non-O1/O139 V. cholerae was the sole pathogen isolated from 21% of the clinical specimens, its role in this epidemic, either alone or in concert with V. cholerae O1, cannot be dismissed. A genomic approach was used to examine similarities and differences among the Haitian V. cholerae O1 and V. cholerae non-O1/O139 strains. A total of 47 V. cholerae O1 and 29 V. cholerae non-O1/O139 isolates from patients and the environment were sequenced. Comparative genome analyses of the 76 genomes and eight reference strains of V. cholerae isolated in concurrent epidemics outside Haiti and 27 V. cholerae genomes available in the public database demonstrated substantial diversity of V. cholerae and ongoing flux within its genome.

Draft genome sequences of the diarrheagenic Escherichia coli collection.

We report the draft genome sequences of the collection referred to as the Escherichia coli DECA collection, which was assembled to contain representative isolates of the 15 most common diarrheagenic clones in humans (http://shigatox.net/new/). These genomes represent a valuable resource to the community of researchers who examine these enteric pathogens.

Genetic deletion of trkB.T1 increases neuromuscular function.

Neurotrophin-dependent activation of the tyrosine kinase receptor trkB.FL modulates neuromuscular synapse maintenance and function; however, it is unclear what role the alternative splice variant, truncated trkB (trkB.T1), may have in the peripheral neuromuscular axis. We examined this question in trkB.T1 null mice and demonstrate that in vivo neuromuscular performance and nerve-evoked muscle tension are significantly increased. In vitro assays indicated that the gain-in-function in trkB.T1(-/-) animals resulted specifically from an increased muscle contractility, and increased electrically evoked calcium release. In the trkB.T1 null muscle, we identified an increase in Akt activation in resting muscle as well as a significant increase in trkB.FL and Akt activation in response to contractile activity. On the basis of these findings, we conclude that the trkB signaling pathway might represent a novel target for intervention across diseases characterized by deficits in neuromuscular function.

Whole-genome sequences of Borrelia bissettii, Borrelia valaisiana, and Borrelia spielmanii.

It has been known for decades that human Lyme disease is caused by the three spirochete species Borrelia burgdorferi, Borrelia afzelii, and Borrelia garinii. Recently, Borrelia valaisiana, Borrelia spielmanii, and Borrelia bissettii have been associated with Lyme disease. We report the complete genome sequences of B. valaisiana VS116, B. spielmanii A14S, and B. bissettii DN127.

Strategies for metagenomic-guided whole-community proteomics of complex microbial environments.

Accurate protein identification in large-scale proteomics experiments relies upon a detailed, accurate protein catalogue, which is derived from predictions of open reading frames based on genome sequence data. Integration of mass spectrometry-based proteomics data with computational proteome predictions from environmental metagenomic sequences has been challenging because of the variable overlap between proteomic datasets and corresponding short-read nucleotide sequence data. In this study, we have benchmarked several strategies for increasing microbial peptide spectral matching in metaproteomic datasets using protein predictions generated from matched metagenomic sequences from the same human fecal samples. Additionally, we investigated the impact of mass spectrometry-based filters (high mass accuracy, delta correlation), and de novo peptide sequencing on the number and robustness of peptide-spectrum assignments in these complex datasets. In summary, we find that high mass accuracy peptide measurements searched against non-assembled reads from DNA sequencing of the same samples significantly increased identifiable proteins without sacrificing accuracy.

Crosstalk between B lymphocytes, microbiota and the intestinal epithelium governs immunity versus metabolism in the gut.

Using a systems biology approach, we discovered and dissected a three-way interaction between the immune system, the intestinal epithelium and the microbiota. We found that, in the absence of B cells, or of IgA, and in the presence of the microbiota, the intestinal epithelium launches its own protective mechanisms, upregulating interferon-inducible immune response pathways and simultaneously repressing Gata4-related metabolic functions. This shift in intestinal function leads to lipid malabsorption and decreased deposition of body fat. Network analysis revealed the presence of two interconnected epithelial-cell gene networks, one governing lipid metabolism and another regulating immunity, that were inversely expressed. Gene expression patterns in gut biopsies from individuals with common variable immunodeficiency or with HIV infection and intestinal malabsorption were very similar to those of the B cell-deficient mice, providing a possible explanation for a longstanding enigmatic association between immunodeficiency and defective lipid absorption in humans.

Whole-genome sequences of two Borrelia afzelii and two Borrelia garinii Lyme disease agent isolates.

Human Lyme disease is commonly caused by several species of spirochetes in the Borrelia genus. In Eurasia these species are largely Borrelia afzelii, B. garinii, B. burgdorferi, and B. bavariensis sp. nov. Whole-genome sequencing is an excellent tool for investigating and understanding the influence of bacterial diversity on the pathogenesis and etiology of Lyme disease. We report here the whole-genome sequences of four isolates from two of the Borrelia species that cause human Lyme disease, B. afzelii isolates ACA-1 and PKo and B. garinii isolates PBr and Far04.

Sparse distance-based learning for simultaneous multiclass classification and feature selection of metagenomic data.

MOTIVATION: Direct sequencing of microbes in human ecosystems (the human microbiome) has complemented single genome cultivation and sequencing to understand and explore the impact of commensal microbes on human health. As sequencing technologies improve and costs decline, the sophistication of data has outgrown available computational methods. While several existing machine learning methods have been adapted for analyzing microbiome data recently, there is not yet an efficient and dedicated algorithm available for multiclass classification of human microbiota. RESULTS: By combining instance-based and model-based learning, we propose a novel sparse distance-based learning method for simultaneous class prediction and feature (variable or taxa, which is used interchangeably) selection from multiple treatment populations on the basis of 16S rRNA sequence count data. Our proposed method simultaneously minimizes the intraclass distance and maximizes the interclass distance with many fewer estimated parameters than other methods. It is very efficient for problems with small sample sizes and unbalanced classes, which are common in metagenomic studies. We implemented this method in a MATLAB toolbox called MetaDistance. We also propose several approaches for data normalization and variance stabilization transformation in MetaDistance. We validate this method on several real and simulated 16S rRNA datasets to show that it outperforms existing methods for classifying metagenomic data. This article is the first to address simultaneous multifeature selection and class prediction with metagenomic count data. AVAILABILITY: The MATLAB toolbox is freely available online at http://metadistance.igs.umaryland.edu/. CONTACT: zliu@umm.edu SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Pervasive recombination and sympatric genome diversification driven by frequency-dependent selection in Borrelia burgdorferi, the Lyme disease bacterium.

How genomic diversity within bacterial populations originates and is maintained in the presence of frequent recombination is a central problem in understanding bacterial evolution. Natural populations of Borrelia burgdorferi, the bacterial agent of Lyme disease, consist of diverse genomic groups co-infecting single individual vertebrate hosts and tick vectors. To understand mechanisms of sympatric genome differentiation in B. burgdorferi, we sequenced and compared 23 genomes representing major genomic groups in North America and Europe. Linkage analysis of >13,500 single-nucleotide polymorphisms revealed pervasive horizontal DNA exchanges. Although three times more frequent than point mutation, recombination is localized and weakly affects genome-wide linkage disequilibrium. We show by computer simulations that, while enhancing population fitness, recombination constrains neutral and adaptive divergence among sympatric genomes through periodic selective sweeps. In contrast, simulations of frequency-dependent selection with recombination produced the observed pattern of a large number of sympatric genomic groups associated with major sequence variations at the selected locus. We conclude that negative frequency-dependent selection targeting a small number of surface-antigen loci (ospC in particular) sufficiently explains the maintenance of sympatric genome diversity in B. burgdorferi without adaptive divergence. We suggest that pervasive recombination makes it less likely for local B. burgdorferi genomic groups to achieve host specialization. B. burgdorferi genomic groups in the northeastern United States are thus best viewed as constituting a single bacterial species, whose generalist nature is a key to its rapid spread and human virulence.

Bacillus anthracis comparative genome analysis in support of the Amerithrax investigation.

Before the anthrax letter attacks of 2001, the developing field of microbial forensics relied on microbial genotyping schemes based on a small portion of a genome sequence. Amerithrax, the investigation into the anthrax letter attacks, applied high-resolution whole-genome sequencing and comparative genomics to identify key genetic features of the letters' Bacillus anthracis Ames strain. During systematic microbiological analysis of the spore material from the letters, we identified a number of morphological variants based on phenotypic characteristics and the ability to sporulate. The genomes of these morphological variants were sequenced and compared with that of the B. anthracis Ames ancestor, the progenitor of all B. anthracis Ames strains. Through comparative genomics, we identified four distinct loci with verifiable genetic mutations. Three of the four mutations could be directly linked to sporulation pathways in B. anthracis and more specifically to the regulation of the phosphorylation state of Spo0F, a key regulatory protein in the initiation of the sporulation cascade, thus linking phenotype to genotype. None of these variant genotypes were identified in single-colony environmental B. anthracis Ames isolates associated with the investigation. These genotypes were identified only in B. anthracis morphotypes isolated from the letters, indicating that the variants were not prevalent in the environment, not even the environments associated with the investigation. This study demonstrates the forensic value of systematic microbiological analysis combined with whole-genome sequencing and comparative genomics.

Whole genome sequence of an unusual Borrelia burgdorferi sensu lato isolate.

Human Lyme disease is caused by a number of related Borrelia burgdorferi sensu lato species. We report here the complete genome sequence of Borrelia sp. isolate SV1 from Finland. This isolate is to date the closest known relative of B. burgdorferi sensu stricto, but it is sufficiently genetically distinct from that species that it and its close relatives warrant its candidacy for new-species status. We suggest that this isolate should be named "Borrelia finlandensis."

Whole-genome sequences of thirteen isolates of Borrelia burgdorferi.

Borrelia burgdorferi is a causative agent of Lyme disease in North America and Eurasia. The first complete genome sequence of B. burgdorferi strain 31, available for more than a decade, has assisted research on the pathogenesis of Lyme disease. Because a single genome sequence is not sufficient to understand the relationship between genotypic and geographic variation and disease phenotype, we determined the whole-genome sequences of 13 additional B. burgdorferi isolates that span the range of natural variation. These sequences should allow improved understanding of pathogenesis and provide a foundation for novel detection, diagnosis, and prevention strategies.

Sequencing of Culex quinquefasciatus establishes a platform for mosquito comparative genomics.

Culex quinquefasciatus (the southern house mosquito) is an important mosquito vector of viruses such as West Nile virus and St. Louis encephalitis virus, as well as of nematodes that cause lymphatic filariasis. C. quinquefasciatus is one species within the Culex pipiens species complex and can be found throughout tropical and temperate climates of the world. The ability of C. quinquefasciatus to take blood meals from birds, livestock, and humans contributes to its ability to vector pathogens between species. Here, we describe the genomic sequence of C. quinquefasciatus: Its repertoire of 18,883 protein-coding genes is 22% larger than that of Aedes aegypti and 52% larger than that of Anopheles gambiae with multiple gene-family expansions, including olfactory and gustatory receptors, salivary gland genes, and genes associated with xenobiotic detoxification.

Draft genome sequence of the oilseed species Ricinus communis.

Castor bean (Ricinus communis) is an oilseed crop that belongs to the spurge (Euphorbiaceae) family, which comprises approximately 6,300 species that include cassava (Manihot esculenta), rubber tree (Hevea brasiliensis) and physic nut (Jatropha curcas). It is primarily of economic interest as a source of castor oil, used for the production of high-quality lubricants because of its high proportion of the unusual fatty acid ricinoleic acid. However, castor bean genomics is also relevant to biosecurity as the seeds contain high levels of ricin, a highly toxic, ribosome-inactivating protein. Here we report the draft genome sequence of castor bean (4.6-fold coverage), the first for a member of the Euphorbiaceae. Whereas most of the key genes involved in oil synthesis and turnover are single copy, the number of members of the ricin gene family is larger than previously thought. Comparative genomics analysis suggests the presence of an ancient hexaploidization event that is conserved across the dicotyledonous lineage.

Bacterial diversity in the oral cavity of 10 healthy individuals.

The composition of the oral microbiota from 10 individuals with healthy oral tissues was determined using culture-independent techniques. From each individual, 26 specimens, each from different oral sites at a single point in time, were collected and pooled. An 11th pool was constructed using portions of the subgingival specimens from all 10 individuals. The 16S ribosomal RNA gene was amplified using broad-range bacterial primers, and clone libraries from the individual and subgingival pools were constructed. From a total of 11,368 high-quality, nonchimeric, near full-length sequences, 247 species-level phylotypes (using a 99% sequence identity threshold) and 9 bacterial phyla were identified. At least 15 bacterial genera were conserved among all 10 individuals, with significant interindividual differences at the species and strain level. Comparisons of these oral bacterial sequences with near full-length sequences found previously in the large intestines and feces of other healthy individuals suggest that the mouth and intestinal tract harbor distinct sets of bacteria. Co-occurrence analysis showed significant segregation of taxa when community membership was examined at the level of genus, but not at the level of species, suggesting that ecologically significant, competitive interactions are more apparent at a broader taxonomic level than species. This study is one of the more comprehensive, high-resolution analyses of bacterial diversity within the healthy human mouth to date, and highlights the value of tools from macroecology for enhancing our understanding of bacterial ecology in human health.

Fast, adaptive evolution at a bacterial host-resistance locus: the PFam54 gene array in Borrelia burgdorferi.

Microbial pathogens have evolved sophisticated mechanisms for evasion of host innate and adaptive immunities. PFam54 is the largest paralogous gene family in the genomes of Borrelia burgdorferi, the Lyme disease bacterium. One member of PFam54, the complement-regulator acquiring surface proteins 1 (BbCrasp-1), is able to abort the alternative pathway of complement activation via binding human complement-regulator factor H (FH). The gene coding for BbCRASP-1 exists in a tandem array of PFam54 genes in the B. burgdorferi genome, a result apparently of repeated gene duplications. To help elucidate the functions of the large number of PFam54 genes, we performed phylogenomic and structural analyses of the PFam54 gene array from ten B. burgdorferi genomes. Analyses based on gene tree, genome synteny, and structural models revealed rapid adaptive evolution of this array through gene duplication, gene loss, and functional diversification. Individual PFam54 genes, however, do not show high intra-population sequence polymorphisms as genes providing evasion from adaptive immunity generally do. PFam54 members able to bind human FH are not monophyletic, suggesting that human FH affinity, however strong, is an incidental rather than main function of these PFam54 proteins. The large number of PFam54 genes existing in any single B. burgdorferi genome may target different innate-immunity proteins of a single host species or the same immune protein of a variety of host species. Genetic variability of the PFam54 gene array suggests that universally present PFam54 lineages such as BBA64, BBA65, BBA66, and BBA73 may be better candidates for the development of broad-spectrum vaccines or drugs than strain-restricted lineages such as BbCRASP-1.