Genomic Encyclopedia of Type Strains of the Genus Bifidobacterium

Bifidobacteria represent one of the dominant microbial groups that are present in the gut of various animals, being particularly prevalent during the suckling stage of life of humans and other mammals. However, the overall genome structure of this group of microorganisms remains largely unexplored. Here, we sequenced the genomes of 42 representative (sub)species across the Bifidobacterium genus and used this information to explore the overall genetic picture of this bacterial group. Furthermore, the genomic data described here were used to reconstruct the evolutionary development of the Bifidobacterium genus. This reconstruction suggests that its evolution was substantially influenced by genetic adaptations to obtain access to glycans, thereby representing a common and potent evolutionary force in shaping bifidobacterial genomes.     

Conjugation genes are common throughout the genus Rickettsia and are transmitted horizontally

Rickettsia are endosymbionts of arthropods, some of which are vectored to vertebrates where they cause disease. Recently, it has been found that some Rickettsia strains harbour conjugative plasmids and others encode some conjugative machinery within the bacterial genome. We investigated the distribution of these conjugation genes in a phylogenetically diverse collection of Rickettsia isolated from arthropods. We found that these genes are common throughout the genus and, in stark contrast to other genes in the genome, conjugation genes are frequently horizontally transmitted between strains. There is no evidence to suggest that these genes are preferentially transferred between phylogenetically related strains, which is surprising given that closely related strains infect similar host species. In addition to detecting patterns of horizontal transmission between diverse Rickettsia species, these findings have implications for the evolution of pathogenicity, the evolution of Rickettsia genomes and the genetic manipulation of intracellular bacteria.     

Distribution of genes of toxin-antitoxin systems of MazEF and RelBE families in bifidobacteria from human intestinal microbiota

The in silico analysis of 36 sequenced genomes of bacteria of the Bifidobacterium genus determined the presence of 19 genes of toxin-antitoxin (TA) system that belong to the MazEF and RelBE families, including five mazF and two relE genes that encode toxins and 12 relB genes that encode antitoxins. A high level of gene (at the level of nucleotide changes) and genomic (presence or absence of genes in distinct genomes) polymorphism in the investigated genes was revealed. The highest level of polymorphism was observed in strains of the Bifidobacterium longum species, primarily in relB1-relB10 genes. Gene and genomic polymorphism might be used to identify the strain of B. longum species. PCR analysis of genomic DNA of 30 bifidobacteria strains belonging to the three species, B. longum, B. adolescentis, and B. bifidum, isolated from the intestinal microbiota of astronauts demonstrated the presence of mazF and relB genes. The observed polymorphism of TA genes indicates the presence of differences in bifidobacteria strains isolated from the intestinal microbiota of astronauts before and after space flight and the control group.     

Origin and Genetic Diversity of Diploid Parthenogenetic Artemia in Eurasia

There is wide interest in understanding how genetic diversity is generated and maintained in parthenogenetic lineages, as it will help clarify the debate of the evolution and maintenance of sexual reproduction. There are three mechanisms that can be responsible for the generation of genetic diversity of parthenogenetic lineages: contagious parthenogenesis, repeated hybridization and microorganism infections (e.g. Wolbachia). Brine shrimps of the genus Artemia (Crustacea, Branchiopoda, Anostraca) are a good model system to investigate evolutionary transitions between reproductive systems as they include sexual species and lineages of obligate parthenogenetic populations of different ploidy level, which often co-occur. Diploid parthenogenetic lineages produce occasional fully functional rare males, interspecific hybridization is known to occur, but the mechanisms of origin of asexual lineages are not completely understood. Here we sequenced and analysed fragments of one mitochondrial and two nuclear genes from an extensive set of populations of diploid parthenogenetic Artemia and sexual species from Central and East Asia to investigate the evolutionary origin of diploid parthenogenetic Artemia, and geographic origin of the parental taxa. Our results indicate that there are at least two, possibly three independent and recent maternal origins of parthenogenetic lineages, related to A. urmiana and Artemia sp. from Kazakhstan, but that the nuclear genes are very closely related in all the sexual species and parthenogegetic lineages except for A. sinica, who presumable took no part on the origin of diploid parthenogenetic strains. Our data cannot rule out either hybridization between any of the very closely related Asiatic sexual species or rare events of contagious parthenogenesis via rare males as the contributing mechanisms to the generation of genetic diversity in diploid parthenogenetic Artemia lineages.     

Comparative Genomic and Phylogenetic Approaches to Characterize the Role of Genetic Recombination in Mycobacterial Evolution

The genus Mycobacterium encompasses over one hundred named species of environmental and pathogenic organisms, including the causative agents of devastating human diseases such as tuberculosis and leprosy. The success of these human pathogens is due in part to their ability to rapidly adapt to their changing environment and host. Recombination is the fastest way for bacterial genomes to acquire genetic material, but conflicting results about the extent of recombination in the genus Mycobacterium have been reported. We examined a data set comprising 18 distinct strains from 13 named species for evidence of recombination. Genomic regions common to all strains (accounting for 10% to 22% of the full genomes of all examined species) were aligned and concatenated in the chromosomal order of one mycobacterial reference species. The concatenated sequence was screened for evidence of recombination using a variety of statistical methods, with each proposed event evaluated by comparing maximum-likelihood phylogenies of the recombinant section with the non-recombinant portion of the dataset. Incongruent phylogenies were identified by comparing the site-wise log-likelihoods of each tree using multiple tests. We also used a phylogenomic approach to identify genes that may have been acquired through horizontal transfer from non-mycobacterial sources. The most frequent associated lineages (and potential gene transfer partners) in the Mycobacterium lineage-restricted gene trees are other members of suborder Corynebacterinae, but more-distant partners were identified as well. In two examined cases of potentially frequent and habitat-directed transfer (M. abscessus to Segniliparus and M. smegmatis to Streptomyces), observed sequence distances were small and consistent with a hypothesis of transfer, while in a third case (M. vanbaalenii to Streptomyces) distances were larger. The analyses described here indicate that whereas evidence of recombination in core regions within the genus is relatively sparse, the acquisition of genes from non-mycobacterial lineages is a significant feature of mycobacterial evolution.     

Comparative genomics of Bifidobacterium species isolated from marmosets and humans

The genus Bifidobacterium is purported to have beneficial consequences for human health and is a major component of many gastrointestinal probiotics. Although species of Bifidobacterium are generally at low relative frequency in the adult human gastrointestinal tract, they can constitute high proportions of the gastrointestinal communities of adult marmosets. To identify genes that might be important for the maintenance of Bifidobacterium in adult marmosets, ten strains of Bifidobacterium were isolated from the feces of seven adult marmosets, and their genomes were sequenced. There were six B. reuteri strains, two B. callitrichos strains, one B. myosotis sp. nov. and one B. tissieri sp. nov. among our isolates. Phylogenetic analysis showed that three of the four species we isolated were most closely related to B. bifidum, B. breve and B. longum, which are species found in high abundance in human infants. There were 1357 genes that were shared by at least one strain of B. reuteri, B. callitrichos, B. breve, and B. longum, and 987 genes that were found in all strains of the four species. There were 106 genes found in B. reuteri and B. callitrichos but not in human bifidobacteria, and several of these genes were involved in nutrient uptake. These pathways for nutrient uptake appeared to be specific to Bifidobacterium from New World monkeys. Additionally, the distribution of Bifidobacterium in fecal samples from captive adult marmosets constituted as much as 80% of the gut microbiome, although this was variable between individuals and colonies. We suggest that nutrient transporters may be important for the maintenance of Bifidobacterium during adulthood in marmosets.     

Deoxyribonucleic acid homology among the species of the genus Bifidobacterium isolated from animals

Summary Representatives of the species of the genus Bifidobacterium isolated from bovine rumen, intestinal content of honey bees and pig feces were studied for their deoxyribonucleic acid similarities with DNA-DNA filter hybridization-competition experiments. The species were: B. asteroides, B. indicum and B. coryneforme from honey bees; B. ruminale and B. globosum from rumen; B. suis from pig feces; Mitsuoka's species B. thermophilum and B. pseudolongum. Strains of B. bifidum, B. infantis, B. longum and B. breve from human sources were included in some comparative experiments. The group of bacteria investigated is largely heterogeneous: virtually no homology exists among the species from rumen and pig feces; the species or types from bees, although unrelated with rumen and pig types, share with these a significant portion of the genome; the competitor DNA's from human strains did not at all react with any of the homologous systems. The annealing reaction was carried out under stringent conditions to ensure against mismatching of nucleotide sequences (70°C in Denhardt's reactive mixture). The separation within the genus Bifidobacterium of several specific entities is fully substantiated.This investigation was supported by a research grant of C. N. R. (Consiglio Nazionale delle Ricerche, Roma).     

Pan-genomics of Ochrobactrum species from clinical and environmental origins reveals distinct populations and possible links

Ochrobactrum genus is comprised of soil-dwelling Gram-negative bacteria mainly reported for bioremediation of toxic compounds. Since last few years, mainly two species of this genus, O. intermedium and O. anthropi were documented for causing infections mostly in the immunocompromised patients. Despite such ubiquitous presence, study of adaptation in various niches is still lacking. Thus, to gain insights into the niche adaptation strategies, pan-genome analysis was carried out by comparing 67 genome sequences belonging to Ochrobactrum species. Pan-genome analysis revealed it is an open pan-genome indicative of the continuously evolving nature of the genus. The presence/absence of gene clusters also illustrated the unique presence of antibiotic efflux transporter genes and type IV secretion system genes in the clinical strains while the genes of solvent resistance and exporter pumps in the environmental strains. A phylogenomic investigation based on 75 core genes depicted better and robust phylogenetic resolution and topology than the 16S rRNA gene. To support the pan-genome analysis, individual genomes were also investigated for the mobile genetic elements (MGE), antibiotic resistance genes (ARG), metal resistance genes (MRG) and virulence factors (VF). The analysis revealed the presence of MGE, ARG, and MRG in all the strains which play an important role in the species evolution which is in agreement with the pan-genome analysis. The average nucleotide identity (ANI) based on the genetic relatedness between the Ochrobactrum species indicated a distinction between individual species. Interestingly, the ANI tool was able to classify the Ochrobactrum genomes to the species level which were assigned till the genus level on the NCBI database.     

Specific Identification and Targeted Characterization of Bifidobacterium lactis from Different Environmental Isolates by a Combined Multiplex-PCR Approach

The species Bifidobacterium lactis, with its main representative strain Bb12 (DSM 10140), is a yoghurt isolate used as a probiotic strain and is commercially applied in different types of yoghurts and infant formulas. In order to ensure the genetic identity and safety of this bacterial isolate, species- and strain-specific molecular tools for genetic fingerprinting must be available to identify isolated bifidobacteria or lactic acid bacteria from, e.g., various clinical environments of relevance in medical microbiology. Two opposing rRNA gene-targeted primers have been developed for specific detection of this microorganism by PCR. The specificity of this approach was evaluated and verified with DNA samples isolated from single and mixed cultures of bifidobacteria and lactobacilli (48 isolates, including the type strains of 29 Bifidobacterium and 9 Lactobacillus species). Furthermore, we performed a Multiplex-PCR using oligonucleotide primers targeting a specific region of the 16S rRNA gene for the genus Bifidobacterium and a conserved eubacterial 16S rDNA sequence. The specificity and sensitivity of this detection with a pure culture of B. lactis were, respectively, 100 bacteria/ml after 25 cycles of PCR and 1 to 10 bacteria/ml after a 50-cycle nested-PCR approach.     

Intraspecies Genomic Diversity and Long-Term Persistence of Bifidobacterium longum

Members of genus Bifidobacterium are Gram-positive bacteria, representing a large part of the human infant microbiota and moderately common in adults. However, our knowledge about their diversity, intraspecific phylogeny and long-term persistence in humans is still limited. Bifidobacterium longum is generally considered to be the most common and prevalent species in the intestinal microbiota. In this work we studied whole genome sequences of 28 strains of B. longum, including 8 sequences described in this paper. Part of these strains were isolated from healthy children during a long observation period (up to 10 years between isolation from the same patient). The three known subspecies (longum, infantis and suis) could be clearly divided using sequence-based phylogenetic methods, gene content and the average nucleotide identity. The profiles of glycoside hydrolase genes reflected the different ecological specializations of these three subspecies. The high impact of horizontal gene transfer on genomic diversity was observed, which is possibly due to a large number of prophages and rapidly spreading plasmids. The pan-genome characteristics of the subspecies longum corresponded to the open pan-genome model. While the major part of the strain-specific genetic loci represented transposons and phage-derived regions, a large number of cell envelope synthesis genes were also observed within this category, representing high variability of cell surface molecules. We observed the cases of isolation of high genetically similar strains of B. longum from the same patients after long periods of time, however, we didn’t succeed in the isolation of genetically identical bacteria: a fact, reflecting the high plasticity of microbiota in children.     

Complete mitochondrial genomes of the human follicle mites Demodex brevis and D. folliculorum: novel gene arrangement,truncated tRNA genes,and ancient divergence between species

Background

Follicle mites of the genus Demodex are found on a wide diversity of mammals, including humans; surprisingly little is known, however, about the evolution of this association. Additional sequence information promises to facilitate studies of Demodex variation within and between host species. Here we report the complete mitochondrial genome sequences of two species of Demodex known to live on humans—Demodex brevis and D. folliculorum—which are the first such genomes available for any member of the genus. We analyzed these sequences to gain insight into the evolution of mitochondrial genomes within the Acariformes. We also used relaxed molecular clock analyses, based on alignments of mitochondrial proteins, to estimate the time of divergence between these two species.

Results

Both Demodex genomes shared a novel gene order that differs substantially from the ancestral chelicerate pattern, with transfer RNA (tRNA) genes apparently having moved much more often than other genes. Mitochondrial tRNA genes of both species were unusually short, with most of them unable to encode tRNAs that could fold into the canonical cloverleaf structure; indeed, several examples lacked both D- and T-arms. Finally, the high level of sequence divergence observed between these species suggests that these two lineages last shared a common ancestor no more recently than about 87 mya.

Conclusions

Among Acariformes, rearrangements involving tRNA genes tend to occur much more often than those involving other genes. The truncated tRNA genes observed in both Demodex species would seem to require the evolution of extensive tRNA editing capabilities and/or coevolved interacting factors. The molecular machinery necessary for these unusual tRNAs to function might provide an avenue for developing treatments of skin disorders caused by Demodex. The deep divergence time estimated between these two species sets a lower bound on the time that Demodex have been coevolving with their mammalian hosts, and supports the hypothesis that there was an early split within the genus Demodex into species that dwell in different skin microhabitats.

Electronic supplementary material

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

Coexistence and Within-Host Evolution of Diversified Lineages of Hypermutable Pseudomonas aeruginosa in Long-term Cystic Fibrosis Infections

The advent of high-throughput sequencing techniques has made it possible to follow the genomic evolution of pathogenic bacteria by comparing longitudinally collected bacteria sampled from human hosts. Such studies in the context of chronic airway infections by Pseudomonas aeruginosa in cystic fibrosis (CF) patients have indicated high bacterial population diversity. Such diversity may be driven by hypermutability resulting from DNA mismatch repair system (MRS) deficiency, a common trait evolved by P. aeruginosa strains in CF infections. No studies to date have utilized whole-genome sequencing to investigate within-host population diversity or long-term evolution of mutators in CF airways. We sequenced the genomes of 13 and 14 isolates of P. aeruginosa mutator populations from an Argentinian and a Danish CF patient, respectively. Our collection of isolates spanned 6 and 20 years of patient infection history, respectively. We sequenced 11 isolates from a single sample from each patient to allow in-depth analysis of population diversity. Each patient was infected by clonal populations of bacteria that were dominated by mutators. The in vivo mutation rate of the populations was ∼100 SNPs/year–∼40-fold higher than rates in normo-mutable populations. Comparison of the genomes of 11 isolates from the same sample showed extensive within-patient genomic diversification; the populations were composed of different sub-lineages that had coexisted for many years since the initial colonization of the patient. Analysis of the mutations identified genes that underwent convergent evolution across lineages and sub-lineages, suggesting that the genes were targeted by mutation to optimize pathogenic fitness. Parallel evolution was observed in reduction of overall catabolic capacity of the populations. These findings are useful for understanding the evolution of pathogen populations and identifying new targets for control of chronic infections.     

Genomic evolution of the class Acidithiobacillia: deep-branching Proteobacteria living in extreme acidic conditions

Members of the genus Acidithiobacillus, now ranked within the class Acidithiobacillia, are model bacteria for the study of chemolithotrophic energy conversion under extreme conditions. Knowledge of the genomic and taxonomic diversity of Acidithiobacillia is still limited. Here, we present a systematic analysis of nearly 100 genomes from the class sampled from a wide range of habitats. Some of these genomes are new and others have been reclassified on the basis of advanced genomic analysis, thus defining 19 Acidithiobacillia lineages ranking at different taxonomic levels. This work provides the most comprehensive classification and pangenomic analysis of this deep-branching class of Proteobacteria to date. The phylogenomic framework obtained illuminates not only the evolutionary past of this lineage, but also the molecular evolution of relevant aerobic respiratory proteins, namely the cytochrome bo₃ ubiquinol oxidases.Subject terms: Bacterial evolution, Bacterial genomics, Phylogenetics     

Diversity of effector genes in plant pathogenic bacteria of genus Xanthomonas

Gram-negative plant pathogenic bacteria are secreting into plant cell a special type of pathogeni city-related proteins called effectors. They are capable of suppressing plant innate immunity or stimulating synthesis and export of metabolites desired by the pathogen. We identified a number of effector-coding genes typical of xanthomonads analyzing 8 completely sequenced genomes of genus Xanthomonas. Using representative collection provided by Russian Research Institute of Phytopathology we identified genetic diversity of effector gene loci in population of Xanthomonas bacteria. Patterns of effector genes were identified for individual strains and statistic linkage between particular genes and race of the pathogen was established. For the first time several untypical effector genes were found in strains of Xanthomonas campestris pv. campestris.     

Intracellular parasitism,the driving force of evolution of Legionella pneumophila and the genus Legionella

Legionella pneumophila is an intracellular pathogen that causes a severe pneumonia called Legionnaires' disease that is often fatal when not promptly diagnosed and treated. Legionella parasitize aquatic protozoa with which it co-evolved over an evolutionary long time. The close relationship between hosts and pathogens, their co-evolution, led to molecular interactions such as the exchange of genetic material through horizontal gene transfer (HGT). Genome sequencing of L. pneumophila and of the entire genus Legionella that comprises over 60 species revealed that Legionellae have co-opted genes and thus cellular functions from their eukaryotic hosts to a surprisingly high extent. Acquisition and loss of these eukaryotic-like genes and domains is an on-going process underlining the highly dynamic nature of the Legionella genomes. Although the large amount and diversity of HGT in Legionella seems to be unique in the prokaryotic world the analyses of more and more genomes from environmental organisms and symbionts of amoeba revealed that such genetic exchanges occur among all amoeba associated bacteria and also among the different microorganisms that infect amoeba. This dynamic reshuffling and gene-acquisition has led to the emergence of Legionella as human pathogen and may lead to the emergence of new human pathogens from the environment.     

Complete Chloroplast Genome Sequence of Omani Lime (Citrus aurantiifolia) and Comparative Analysis within the Rosids

The genus Citrus contains many economically important fruits that are grown worldwide for their high nutritional and medicinal value. Due to frequent hybridizations among species and cultivars, the exact number of natural species and the taxonomic relationships within this genus are unclear. To compare the differences between the Citrus chloroplast genomes and to develop useful genetic markers, we used a reference-assisted approach to assemble the complete chloroplast genome of Omani lime (C. aurantiifolia). The complete C. aurantiifolia chloroplast genome is 159,893 bp in length; the organization and gene content are similar to most of the rosids lineages characterized to date. Through comparison with the sweet orange (C. sinensis) chloroplast genome, we identified three intergenic regions and 94 simple sequence repeats (SSRs) that are potentially informative markers with resolution for interspecific relationships. These markers can be utilized to better understand the origin of cultivated Citrus. A comparison among 72 species belonging to 10 families of representative rosids lineages also provides new insights into their chloroplast genome evolution.     

Recombination in Thermotoga: implications for species concepts and biogeography

Here we characterize regions of the genomes of eight members of the hyperthermophilic genus Thermotoga. These bacteria differ from each other physiologically and by 3-20% in gene content and occupy physically distinct environments in widely disparate regions of the globe. Among the four different lineages (represented by nine different strains) that we compare, no two are closer than 96% in the average sequences of their genes. By most accepted recent definitions these are different "ecotypes" and different "species." And yet we find compelling evidence for recombination between them. We suggest that no single prokaryotic species concept can accommodate such uncoupling of ecotypic and genetic aspects of cohesion and diversity, and that without a single concept, the question of whether or not prokaryotic species might in general be cosmopolitan cannot be sensibly addressed. We can, however, recast biogeographical questions in terms of the distribution of genes and their alleles.     

Comparative Mitochondrial Genomics within and among Yeast Species of the Lachancea Genus

Yeasts are leading model organisms for mitochondrial genome studies. The explosion of complete sequence of yeast mitochondrial (mt) genomes revealed a wide diversity of organization and structure between species. Recently, genome-wide polymorphism survey on the mt genome of isolates of a single species, Lachancea kluyveri, was also performed. To compare the mitochondrial genome evolution at two hierarchical levels: within and among closely related species, we focused on five species of the Lachancea genus, which are close relatives of L. kluyveri. Hence, we sequenced the complete mt genome of L. dasiensis, L. nothofagi, L. mirantina, L. fantastica and L. meyersii. The phylogeny of the Lachancea genus was explored using these data. Analysis of intra- and interspecific variability across the whole Lachancea genus led to the same conclusions regarding the mitochondrial genome evolution. These genomes exhibit a similar architecture and are completely syntenic. Nevertheless, genome sizes vary considerably because of the variations of the intergenic regions and the intron content, contributing to mitochondrial genome plasticity. The high variability of the intergenic regions stands in contrast to the high level of similarity of protein sequences. Quantification of the selective constraints clearly revealed that most of the mitochondrial genes are under purifying selection in the whole genus.     

Genus-Wide Comparative Genomics of Malassezia Delineates Its Phylogeny,Physiology, and Niche Adaptation on Human Skin

Malassezia is a unique lipophilic genus in class Malasseziomycetes in Ustilaginomycotina, (Basidiomycota, fungi) that otherwise consists almost exclusively of plant pathogens. Malassezia are typically isolated from warm-blooded animals, are dominant members of the human skin mycobiome and are associated with common skin disorders. To characterize the genetic basis of the unique phenotypes of Malassezia spp., we sequenced the genomes of all 14 accepted species and used comparative genomics against a broad panel of fungal genomes to comprehensively identify distinct features that define the Malassezia gene repertoire: gene gain and loss; selection signatures; and lineage-specific gene family expansions. Our analysis revealed key gene gain events (64) with a single gene conserved across all Malassezia but absent in all other sequenced Basidiomycota. These likely horizontally transferred genes provide intriguing gain-of-function events and prime candidates to explain the emergence of Malassezia. A larger set of genes (741) were lost, with enrichment for glycosyl hydrolases and carbohydrate metabolism, concordant with adaptation to skin’s carbohydrate-deficient environment. Gene family analysis revealed extensive turnover and underlined the importance of secretory lipases, phospholipases, aspartyl proteases, and other peptidases. Combining genomic analysis with a re-evaluation of culture characteristics, we establish the likely lipid-dependence of all Malassezia. Our phylogenetic analysis sheds new light on the relationship between Malassezia and other members of Ustilaginomycotina, as well as phylogenetic lineages within the genus. Overall, our study provides a unique genomic resource for understanding Malassezia niche-specificity and potential virulence, as well as their abundance and distribution in the environment and on human skin.     

Investigation of the Evolutionary Development of the Genus Bifidobacterium by Comparative Genomics

The Bifidobacterium genus currently encompasses 48 recognized taxa, which have been isolated from different ecosystems. However, the current phylogeny of bifidobacteria is hampered by the relative paucity of genotypic data. Here, we reassessed the taxonomy of this bacterial genus using genome-based approaches, which demonstrated that the previous taxonomic view of bifidobacteria contained several inconsistencies. In particular, high levels of genetic relatedness were shown to exist between particular Bifidobacterium taxa which would not justify their status as separate species. The results presented are here based on average nucleotide identity analysis involving the genome sequences for each type strain of the 48 bifidobacterial taxa, as well as phylogenetic comparative analysis of the predicted core genome of the Bifidobacterium genus. The results of this study demonstrate that the availability of complete genome sequences allows the reconstruction of a more robust bifidobacterial phylogeny than that obtained from a single gene-based sequence comparison, thus discouraging the assignment of a new or separate bifidobacterial taxon without such a genome-based validation.