A phylogenetic analysis of the phylum Fibrobacteres

Members of the phylum Fibrobacteres are highly efficient cellulolytic bacteria, best known for their role in rumen function and as potential sources of novel enzymes for bioenergy applications. Despite being key members of ruminants and other digestive microbial communities, our knowledge of this phylum remains incomplete, as much of our understanding is focused on two recognized species, Fibrobacter succinogenes and F. intestinalis. As a result, we lack insights regarding the environmental niche, host range, and phylogenetic organization of this phylum. Here, we analyzed over 1000 16S rRNA Fibrobacteres sequences available from public databases to establish a phylogenetic framework for this phylum. We identify both species- and genus-level clades that are suggestive of previously unknown taxonomic relationships between Fibrobacteres in addition to their putative lifestyles as host-associated or free-living. Our results shed light on this poorly understood phylum and will be useful for elucidating the function, distribution, and diversity of these bacteria in their niches.     

Molecular signatures for the phylum Synergistetes and some of its subclades

Species belonging to the phylum Synergistetes are poorly characterized. Though the known species display Gram-negative characteristics and the ability to ferment amino acids, no single characteristic is known which can define this group. For eight Synergistetes species, complete genome sequences or draft genomes have become available. We have used these genomes to construct detailed phylogenetic trees for the Synergistetes species and carried out comprehensive analysis to identify molecular markers consisting of conserved signature indels (CSIs) in protein sequences that are specific for either all Synergistetes or some of their sub-groups. We report here identification of 32 CSIs in widely distributed proteins such as RpoB, RpoC, UvrD, GyrA, PolA, PolC, MraW, NadD, PyrE, RpsA, RpsH, FtsA, RadA, etc., including a large >300 aa insert within the RpoC protein, that are present in various Synergistetes species, but except for isolated bacteria, these CSIs are not found in the protein homologues from any other organisms. These CSIs provide novel molecular markers that distinguish the species of the phylum Synergistetes from all other bacteria. The large numbers of other CSIs discovered in this work provide valuable information that supports and consolidates evolutionary relationships amongst the sequenced Synergistetes species. Of these CSIs, seven are specifically present in Jonquetella, Pyramidobacter and Dethiosulfovibrio species indicating a cladal relationship among them, which is also strongly supported by phylogenetic trees. A further 15 CSIs that are only present in Jonquetella and Pyramidobacter indicate a close association between these two species. Additionally, a previously described phylogenetic relationship between the Aminomonas and Thermanaerovibrio species was also supported by 9 CSIs. The strong relationships indicated by the indel analysis provide incentives for the grouping of species from these clades into higher taxonomic groups such as families or orders. The identified molecular markers, due to their specificity for Synergistetes and presence in highly conserved regions of important proteins suggest novel targets for evolutionary, genetic and biochemical studies on these bacteria as well as for the identification of additional species belonging to this phylum in different environments.     

Evolutionary Analysis of Mitogenomes from Parasitic and Free-Living Flatworms

Mitochondrial genomes (mitogenomes) are useful and relatively accessible sources of molecular data to explore and understand the evolutionary history and relationships of eukaryotic organisms across diverse taxonomic levels. The availability of complete mitogenomes from Platyhelminthes is limited; of the 40 or so published most are from parasitic flatworms (Neodermata). Here, we present the mitogenomes of two free-living flatworms (Tricladida): the complete genome of the freshwater species Crenobia alpina (Planariidae) and a nearly complete genome of the land planarian Obama sp. (Geoplanidae). Moreover, we have reanotated the published mitogenome of the species Dugesia japonica (Dugesiidae). This contribution almost doubles the total number of mtDNAs published for Tricladida, a species-rich group including model organisms and economically important invasive species. We took the opportunity to conduct comparative mitogenomic analyses between available free-living and selected parasitic flatworms in order to gain insights into the putative effect of life cycle on nucleotide composition through mutation and natural selection. Unexpectedly, we did not find any molecular hallmark of a selective relaxation in mitogenomes of parasitic flatworms; on the contrary, three out of the four studied free-living triclad mitogenomes exhibit higher A+T content and selective relaxation levels. Additionally, we provide new and valuable molecular data to develop markers for future phylogenetic studies on planariids and geoplanids.     

Phylogenetic perspective on the relationships and evolutionary history of the Acipenseriformes

The Acipenseriformes, as one of the earliest extant vertebrates, plays an important role in the evolution of fishes and even the whole vertebrates. Here we collected and analyzed all complete mitochondrial genomes of Acipenseriformes species. Phylogenetic analyses demonstrated that the polytomous branch included Acipenseridae and Polyodontidae formed five clades. The Polyodontidae clade and the Scaphirhynchus clade both were monophyletic group, whereas the Acipenser species and the Huso species both were polyphyletic group. The Bayesian divergence times showed that the origin time for Acipenseriformes was at 318.0 Mya, which was similar to the some previous results of 312.1 Mya, 346.9 Mya and 389.7 Mya. The result was in good consistent with the paleontological data available and the split time of the Pacific and Atlantic Oceans from the Jurassic to the Cretaceous (Laurasia splits in North America and Eurasia). The dN/dS ratios showed the evolutionary rates gradually slow down in five major Acipenseriformes clades from the Clade A (the Pacific sturgeons species) to Clade C (the genus Scaphirhynchus), which was related to the process of geographical formation.     

Phylogenetic Framework and Molecular Signatures for the Main Clades of the Phylum Actinobacteria

Summary: The phylum Actinobacteria harbors many important human pathogens and also provides one of the richest sources of natural products, including numerous antibiotics and other compounds of biotechnological interest. Thus, a reliable phylogeny of this large phylum and the means to accurately identify its different constituent groups are of much interest. Detailed phylogenetic and comparative analyses of >150 actinobacterial genomes reported here form the basis for achieving these objectives. In phylogenetic trees based upon 35 conserved proteins, most of the main groups of Actinobacteria as well as a number of their superageneric clades are resolved. We also describe large numbers of molecular markers consisting of conserved signature indels in protein sequences and whole proteins that are specific for either all Actinobacteria or their different clades (viz., orders, families, genera, and subgenera) at various taxonomic levels. These signatures independently support the existence of different phylogenetic clades, and based upon them, it is now possible to delimit the phylum Actinobacteria (excluding Coriobacteriia) and most of its major groups in clear molecular terms. The species distribution patterns of these markers also provide important information regarding the interrelationships among different main orders of Actinobacteria. The identified molecular markers, in addition to enabling the development of a stable and reliable phylogenetic framework for this phylum, also provide novel and powerful means for the identification of different groups of Actinobacteria in diverse environments. Genetic and biochemical studies on these Actinobacteria-specific markers should lead to the discovery of novel biochemical and/or other properties that are unique to different groups of Actinobacteria.     

Massive comparative genomic analysis reveals convergent evolution of specialized bacteria

Background  

Genome size and gene content in bacteria are associated with their lifestyles. Obligate intracellular bacteria (i.e., mutualists and parasites) have small genomes that derived from larger free-living bacterial ancestors; however, the different steps of bacterial specialization from free-living to intracellular lifestyle have not been studied comprehensively. The growing number of available sequenced genomes makes it possible to perform a statistical comparative analysis of 317 genomes from bacteria with different lifestyles.     

Chromosome-level genomes provide insights into genome evolution,organization and size in Epichloe fungi

Epichloe fungi are endophytes of cool season grasses, both wild species and commercial cultivars, where they may exhibit mutualistic or pathogenic lifestyles. The Epichloe-grass symbiosis is of great interest to agricultural research for the fungal bioprotective properties conferred to host grasses but also serves as an ideal system to study the evolution of fungal plant-pathogens in natural environments. Here, we assembled and annotated gapless chromosome-level genomes of two pathogenic Epichloe sibling species. Both genomes have a bipartite genome organization, with blocks of highly syntenic gene-rich regions separated by blocks of AT-rich DNA. The AT-rich regions show an extensive signature of RIP (repeat-induced point mutation) and the expansion of this compartment accounts for the large difference in genome size between the two species. This study reveals how the rapid evolution of repeat structure can drive divergence between closely related taxa and highlights the evolutionary role of dynamic compartments in fungal genomes.     

Reproductive mode evolution in nematodes: insights from molecular phylogenies and recently discovered species

The Phylum Nematoda has long been known to contain a great diversity of species that vary in reproductive mode, though our understanding of the evolutionary origins, causes and consequences of nematode reproductive mode change have only recently started to mature. Here we bring together and analyze recent progress on reproductive mode evolution throughout the phylum, resulting from the application of molecular phylogenetic approaches and newly discovered nematode species. Reproductive mode variation is reviewed in multiple free-living, animal-parasitic and plant-parasitic nematode groups. Discussion ranges from the model nematode Caenorhabditis elegans and its close relatives, to the plant-parasitic nematodes of the Meloidogyne genus where there is extreme variation in reproductive mode between and even within species, to the vertebrate-parasitic genus Strongyloides and related genera where reproductive mode varies across generations (heterogony). Multiple evolutionary transitions from dioecous (obligately outcrossing) to hermaphroditism and parthenogenesis in the phylum are discussed, along with one case of an evolutionary transition from hermaphroditism to doioecy in the Oscheius genus. We consider the roles of underlying genetic mechanisms in promoting reproductive plasticity in this phylum, as well as the potential evolutionary forces promoting transitions in reproductive mode.     

Genome-Wide Comparative Analysis Reveals Similar Types of NBS Genes in Hybrid Citrus sinensis Genome and Original Citrus clementine Genome and Provides New Insights into Non-TIR NBS Genes

In this study, we identified and compared nucleotide-binding site (NBS) domain-containing genes from three Citrus genomes (C. clementina, C. sinensis from USA and C. sinensis from China). Phylogenetic analysis of all Citrus NBS genes across these three genomes revealed that there are three approximately evenly numbered groups: one group contains the Toll-Interleukin receptor (TIR) domain and two different Non-TIR groups in which most of proteins contain the Coiled Coil (CC) domain. Motif analysis confirmed that the two groups of CC-containing NBS genes are from different evolutionary origins. We partitioned NBS genes into clades using NBS domain sequence distances and found most clades include NBS genes from all three Citrus genomes. This suggests that three Citrus genomes have similar numbers and types of NBS genes. We also mapped the re-sequenced reads of three pomelo and three mandarin genomes onto the C. sinensis genome. We found that most NBS genes of the hybrid C. sinensis genome have corresponding homologous genes in both pomelo and mandarin genomes. The homologous NBS genes in pomelo and mandarin suggest that the parental species of C. sinensis may contain similar types of NBS genes. This explains why the hybrid C. sinensis and original C. clementina have similar types of NBS genes in this study. Furthermore, we found that sequence variation amongst Citrus NBS genes were shaped by multiple independent and shared accelerated mutation accumulation events among different groups of NBS genes and in different Citrus genomes. Our comparative analyses yield valuable insight into the structure, organization and evolution of NBS genes in Citrus genomes. Furthermore, our comprehensive analysis showed that the non-TIR NBS genes can be divided into two groups that come from different evolutionary origins. This provides new insights into non-TIR genes, which have not received much attention.     

Phylogeny and molecular signatures for the phylum Thermotogae and its subgroups

Thermotogae species are currently identified mainly on the basis of their unique toga and distinct branching in the rRNA and other phylogenetic trees. No biochemical or molecular markers are known that clearly distinguish the species from this phylum from all other bacteria. The taxonomic/evolutionary relationships within this phylum, which consists of a single family, are also unclear. We report detailed phylogenetic analyses on Thermotogae species based on concatenated sequences for many ribosomal as well as other conserved proteins that identify a number of distinct clades within this phylum. Additionally, comprehensive analyses of protein sequences from Thermotogae genomes have identified >60 Conserved Signature Indels (CSI) that are specific for the Thermotogae phylum or its different subgroups. Eighteen CSIs in important proteins such as PolI, RecA, TrpRS and ribosomal proteins L4, L7/L12, S8, S9, etc. are uniquely present in various Thermotogae species and provide molecular markers for the phylum. Many CSIs were specific for a number of Thermotogae subgroups. Twelve of these CSIs were specific for a clade consisting of various Thermotoga species except Tt. lettingae, which was separated from other Thermotoga species by a long branch in phylogenetic trees; Fourteen CSIs were specific for a clade consisting of the Fervidobacterium and Thermosipho genera and eight additional CSIs were specific for the genus Thermosipho. In addition, the existence of a clade consisting of the deep branching species Petrotoga mobilis, Kosmotoga olearia and Thermotogales bacterium mesG1 was supported by seven CSIs. The deep branching of this clade was also supported by a number of CSIs that were present in various Thermotogae species, but absent in this clade and all other bacteria. Most of these clades were strongly supported by phylogenetic analyses based on two datasets of protein sequences and they identify potential higher taxonomic grouping (viz. families) within this phylum. We also report 16 CSIs that are shared by either some or all Thermotogae species and some species from other taxa such as Archaea, Aquificae, Firmicutes, Proteobacteria, Deinococcus, Fusobacteria, Dictyoglomus, Chloroflexi and eukaryotes. The shared presence of some of these CSIs could be due to lateral gene transfers between these groups. However, no clear preference for any particular group was observed in this regard. The molecular probes based on different genes/proteins, which contain these Thermotogae-specific CSIs, provide novel and highly specific means for identification of both known as well as previously unknown Thermotogae species in different environments. Additionally, these CSIs also provide valuable tools for genetic and biochemical studies that could lead to discovery of novel properties that are unique to these bacteria.     

Phylum-Level Conservation of Regulatory Information in Nematodes despite Extensive Non-coding Sequence Divergence

Gene regulatory information guides development and shapes the course of evolution. To test conservation of gene regulation within the phylum Nematoda, we compared the functions of putative cis-regulatory sequences of four sets of orthologs (unc-47, unc-25, mec-3 and elt-2) from distantly-related nematode species. These species, Caenorhabditis elegans, its congeneric C. briggsae, and three parasitic species Meloidogyne hapla, Brugia malayi, and Trichinella spiralis, represent four of the five major clades in the phylum Nematoda. Despite the great phylogenetic distances sampled and the extensive sequence divergence of nematode genomes, all but one of the regulatory elements we tested are able to drive at least a subset of the expected gene expression patterns. We show that functionally conserved cis-regulatory elements have no more extended sequence similarity to their C. elegans orthologs than would be expected by chance, but they do harbor motifs that are important for proper expression of the C. elegans genes. These motifs are too short to be distinguished from the background level of sequence similarity, and while identical in sequence they are not conserved in orientation or position. Functional tests reveal that some of these motifs contribute to proper expression. Our results suggest that conserved regulatory circuitry can persist despite considerable turnover within cis elements.     

Genomics and phylogeny of the proposed phylum ‘Candidatus Poribacteria’ associated with the excavating sponge Thoosa mismalolli

Members of the proposed phylum ‘Candidatus Poribacteria’ are among the most abundant microorganisms in the highly diverse microbiome of the sponge mesohyl. Genomic and phylogenetic characteristics of this proposed phylum are barely known. In this study, we analyzed metagenome-assembled genomes (MAGs) obtained from the coral reef excavating sponge Thoosa mismalolli from the Mexican Pacific Ocean. Two MAGs were extracted and analyzed together with 32 MAGs and single-amplified genomes (SAGs) obtained from NCBI. The phylogenetic tree based on the sequences of 139 single-copy genes (SCG) showed two clades. Clade A (23 genomes) represented 67.7% of the total of the genomes, while clade B (11 genomes) comprised 32.3% of the genomes. The Average Nucleotide Identity (ANI) showed values between 66 and 99% for the genomes of the proposed phylum, and the pangenome of genomes revealed a total of 37,234 genes that included 1722 core gene. The number of genes used in the phylogenetic analysis increased from 28 (previous studies) to 139 (this study), which allowed a better resolution of the phylogeny of the proposed phylum. The results supported the two previously described classes, ‘Candidatus Entoporibacteria’ and ‘Candidatus Pelagiporibacteria’, and the genomes SB0101 and SB0202 obtained in this study belong to two new species of the class ‘Candidatus Entoporibacteria’. This is the first comparative study that includes MAGs from a non-sponge host (Porites lutea) to elucidate the taxonomy of the poorly known Candidatus phylum in a polyphasic approach. Finally, our study also contributes to the sponge microbiome project by reporting the first MAGs of the proposed phylum ‘Candidatus Poribacteria’ isolated from the excavating sponge T. mismalolli.

     

Comparative Genomics and the Gene Complement of a Minimal Cell

The concept of a minimal cell is discussed from the viewpoint of comparative genomics. Analysis of published DNA content values determined for 641 different archaeal and bacterial species by pulsed field gel electrophoresis has lead to a more precise definition of the genome size ranges of free-living and host-associated organisms. DNA content is not an indicator of phylogenetic position. However, the smallest genomes in our sample do not have a random distribution in rRNA-based evolutionary trees, and are found mostly in (a) the basal branches of the tree where thermophiles are located; and (b) in late clades, such as those of Gram positive bacteria. While the smallest-known genome size for an endosymbiont is only 450 kb, no free-living prokaryote has been described to have genomes < 1450 kb. Estimates of the size of minimal gene complement can provide important insights in the primary biological functions required for a sustainable, reproducing cell nowadays and throughout evolutionary times, but definitions of the minimum cell is dependent on specific environments.     

The Smallest Known Genomes of Multicellular and Toxic Cyanobacteria: Comparison,Minimal Gene Sets for Linked Traits and the Evolutionary Implications

Cyanobacterial morphology is diverse, ranging from unicellular spheres or rods to multicellular structures such as colonies and filaments. Multicellular species represent an evolutionary strategy to differentiate and compartmentalize certain metabolic functions for reproduction and nitrogen (N₂) fixation into specialized cell types (e.g. akinetes, heterocysts and diazocytes). Only a few filamentous, differentiated cyanobacterial species, with genome sizes over 5 Mb, have been sequenced. We sequenced the genomes of two strains of closely related filamentous cyanobacterial species to yield further insights into the molecular basis of the traits of N₂ fixation, filament formation and cell differentiation. Cylindrospermopsis raciborskii CS-505 is a cylindrospermopsin-producing strain from Australia, whereas Raphidiopsis brookii D9 from Brazil synthesizes neurotoxins associated with paralytic shellfish poisoning (PSP). Despite their different morphology, toxin composition and disjunct geographical distribution, these strains form a monophyletic group. With genome sizes of approximately 3.9 (CS-505) and 3.2 (D9) Mb, these are the smallest genomes described for free-living filamentous cyanobacteria. We observed remarkable gene order conservation (synteny) between these genomes despite the difference in repetitive element content, which accounts for most of the genome size difference between them. We show here that the strains share a specific set of 2539 genes with >90% average nucleotide identity. The fact that the CS-505 and D9 genomes are small and streamlined compared to those of other filamentous cyanobacterial species and the lack of the ability for heterocyst formation in strain D9 allowed us to define a core set of genes responsible for each trait in filamentous species. We presume that in strain D9 the ability to form proper heterocysts was secondarily lost together with N₂ fixation capacity. Further comparisons to all available cyanobacterial genomes covering almost the entire evolutionary branch revealed a common minimal gene set for each of these cyanobacterial traits.     

The Search for Therapeutic Bacteriophages Uncovers One New Subfamily and Two New Genera of Pseudomonas-Infecting Myoviridae

In a previous study, six virulent bacteriophages PAK_P1, PAK_P2, PAK_P3, PAK_P4, PAK_P5 and CHA_P1 were evaluated for their in vivo efficacy in treating Pseudomonas aeruginosa infections using a mouse model of lung infection. Here, we show that their genomes are closely related to five other Pseudomonas phages and allow a subdivision into two clades, PAK_P1-like and KPP10-like viruses, based on differences in genome size, %GC and genomic contents, as well as number of tRNAs. These two clades are well delineated, with a mean of 86% and 92% of proteins considered homologous within individual clades, and 25% proteins considered homologous between the two clades. By ESI-MS/MS analysis we determined that their virions are composed of at least 25 different proteins and electron microscopy revealed a morphology identical to the hallmark Salmonella phage Felix O1. A search for additional bacteriophage homologs, using profiles of protein families defined from the analysis of the 11 genomes, identified 10 additional candidates infecting hosts from different species. By carrying out a phylogenetic analysis using these 21 genomes we were able to define a new subfamily of viruses, the Felixounavirinae within the Myoviridae family. The new Felixounavirinae subfamily includes three genera: Felixounalikevirus, PAK_P1likevirus and KPP10likevirus. Sequencing genomes of bacteriophages with therapeutic potential increases the quantity of genomic data on closely related bacteriophages, leading to establishment of new taxonomic clades and the development of strategies for analyzing viral genomes as presented in this article.     

Old Lineages in a New Ecosystem: Diversification of Arcellinid Amoebae (Amoebozoa) and Peatland Mosses

Arcellinid testate amoebae (Amoebozoa) form a group of free-living microbial eukaryotes with one of the oldest fossil records known, yet several aspects of their evolutionary history remain poorly understood. Arcellinids occur in a range of terrestrial, freshwater and even brackish habitats; however, many arcellinid morphospecies such as Hyalosphenia papilio are particularly abundant in Sphagnum-dominated peatlands, a relatively new ecosystem that appeared during the diversification of Sphagnum species in the Miocene (5–20 Myr ago). Here, we reconstruct divergence times in arcellinid testate amoebae after selecting several fossils for clock calibrations and then infer whether or not arcellinids followed a pattern of diversification that parallels the pattern described for Sphagnum. We found that the diversification of core arcellinids occurred during the Phanerozoic, which is congruent with most arcellinid fossils but not with the oldest known amoebozoan fossil (i.e. at ca. 662 or ca. 750 Myr). Overall, Sphagnum and the Hyalospheniidae exhibit different patterns of diversification. However, an extensive molecular phylogenetic analysis of distinct clades within H. papilio species complex demonstrated a correlation between the recent diversification of H. papilio, the recent diversification of Sphagnum mosses, and the establishment of peatlands.     

Candidatus Prosiliicoccus vernus,a spring phytoplankton bloom associated member of the Flavobacteriaceae

Microbial degradation of algal biomass following spring phytoplankton blooms has been characterised as a concerted effort among multiple clades of heterotrophic bacteria. Despite their significance to overall carbon turnover, many of these clades have resisted cultivation. One clade known from 16S rRNA gene sequencing surveys at Helgoland in the North Sea, was formerly identified as belonging to the genus Ulvibacter. This clade rapidly responds to algal blooms, transiently making up as much as 20% of the free-living bacterioplankton. Sequence similarity below 95% between the 16S rRNA genes of described Ulvibacter species and those from Helgoland suggest this is a novel genus. Analysis of 40 metagenome assembled genomes (MAGs) derived from samples collected during spring blooms at Helgoland support this conclusion. These MAGs represent three species, only one of which appears to bloom in response to phytoplankton. MAGs with estimated completeness greater than 90% could only be recovered for this abundant species. Additional, less complete, MAGs belonging to all three species were recovered from a mini-metagenome of cells sorted via flow cytometry using the genus specific ULV995 fluorescent rRNA probe. Metabolic reconstruction indicates this highly abundant species most likely degrades proteins and the polysaccharide laminarin. Fluorescence in situ hybridisation showed coccoid cells, with a mean diameter of 0.78 mm, with standard deviation of 0.12 μm. Based on the phylogenetic and genomic characteristics of this clade, we propose the novel candidate genus Candidatus Prosiliicoccus, and for the most abundant and well characterised of the three species the name Candidatus Prosiliicoccus vernus.     

Phylogenetic analysis of the Argonaute protein family in platyhelminths

Argonaute proteins (AGOs) are mediators of gene silencing via recruitment of small regulatory RNAs to induce translational regression or degradation of targeted molecules. Platyhelminths have been reported to express microRNAs but the diversity of AGOs in the phylum has not been explored. Phylogenetic relationships of members of this protein family were studied using data from six platyhelminth genomes. Phylogenetic analysis showed that all cestode and trematode AGOs, along with some triclad planarian AGOs, were grouped into the Ago subfamily and its novel sister clade, here referred to as Cluster 1. These were very distant from Piwi and Class 3 subfamilies. By contrast, a number of planarian Piwi-like AGOs formed a novel sister clade to the Piwi subfamily. Extensive sequence searching revealed the presence of an additional locus for AGO2 in the cestode Echinococcus granulosus and exon expansion in this species and E. multilocularis. The current study suggests the absence of the Piwi subfamily and Class 3 AGOs in cestodes and trematodes and the Piwi-like AGO expansion in a free-living triclad planarian and the occurrence of exon expansion prior to or during the evolution of the most-recent common ancestor of the Echinococcus species studied.