Rare genomic changes as a tool for phylogenetics

High-Throughput Sequencing of Complete Mitochondrial Genomes

Next-generation sequencing has revolutionized mitogenomics, turning a cottage industry into a high throughput process. This chapter outlines methodologies used to sequence, assemble, and annotate mitogenomes of non-model organisms using Illumina sequencing technology, utilizing either long-range PCR amplicons or gDNA as starting template. Instructions are given on how to extract DNA, conduct long-range PCR amplifications, generate short Sanger barcode tag sequences, prepare equimolar sample pools, construct and assess quality library preparations, assemble Illumina reads using either seeded reference mapping or de novo assembly, and annotate mitogenomes in the absence of an automated pipeline. Notes and recommendations, derived from our own experience, are given throughout this chapter.

Determining the Position of Storks on the Phylogenetic Tree of Waterbirds by Retroposon Insertion Analysis

Despite many studies on avian phylogenetics in recent decades that used morphology, mitochondrial genomes, and/or nuclear genes, the phylogenetic positions of several birds (e.g., storks) remain unsettled. In addition to the aforementioned approaches, analysis of retroposon insertions, which are nearly homoplasy-free phylogenetic markers, has also been used in avian phylogenetics. However, the first step in the analysis of retroposon insertions, that is, isolation of retroposons from genomic libraries, is a costly and time-consuming procedure. Therefore, we developed a high-throughput and cost-effective protocol to collect retroposon insertion information based on next-generation sequencing technology, which we call here the STRONG (Screening of Transposons Obtained by Next Generation Sequencing) method, and applied it to 3 waterbird species, for which we identified 35,470 loci containing chicken repeat 1 retroposons (CR1). Our analysis of the presence/absence of 30 CR1 insertions demonstrated the intra- and interordinal phylogenetic relationships in the waterbird assemblage, namely 1) Loons diverged first among the waterbirds, 2) penguins (Sphenisciformes) and petrels (Procellariiformes) diverged next, and 3) among the remaining families of waterbirds traditionally classified in Ciconiiformes/Pelecaniformes, storks (Ciconiidae) diverged first. Furthermore, our genome-scale, in silico retroposon analysis based on published genome data uncovered a complex divergence history among pelican, heron, and ibis lineages, presumably involving ancient interspecies hybridization between the heron and ibis lineages. Thus, our retroposon-based waterbird phylogeny and the established phylogenetic position of storks will help to understand the evolutionary processes of aquatic adaptation and related morphological convergent evolution.

The complete mitochondrial genome of Anoplocephala perfoliata, the first representative for the family Anoplocephalidae

Background

Mitochondrial (mt) genome sequences are widely used to understand phylogenetic relationships among parasites. However, no complete mt genome sequence is available in the family Anoplocephalidae to date. This study sequenced and annotated the complete mt genome of Anoplocephala perfoliata (Anoplocephalidae), and investigated its phylogenetic relationships with other species from the families Hymenolepididae, Dipylidiidae and Taeniidae of the order Cyclophyllidea using the amino acid sequences of the 12 proteins in their mt genomes.

Methods

The complete mt genome of A. perfoliata was amplified by Long-range PCR, sequenced using primer walking and annotated by comparing with those of other cestodes. Its phylogenetic relationship with the species from the families Hymenolepididae, Dipylidiidae and Taeniidae was inferred using the 12 protein sequences based on Maximum likelihood and Bayesian methods.

Results

The complete circular mt genome sequence for A. perfoliata is 14,459 bp in size, and includes 12 protein-coding genes, 2 rRNA genes and 22 tRNA genes. The mt gene arrangement of A. perfoliata is identical to those of previously reported Hymenolepis diminuta (Hymenolepididae) and Dipylidium caninum (Dipylidiidae), but slightly different from those of other taeniids due to an order switch between tRNA(S2) and tRNA(L1). The phylogenetic analyses showed that the Dipylidiidae was more closely related to Anoplocephalidae and Hymenolepididae than to Taeniidae. The relationship among the four families obtained by Maximum likelihood and Bayesian inferences based on predicted amino acid sequences of protein-coding genes is consistent with that based on their mt gene arrangement similarities.

Conclusions

This study determined the first mt genome for the family Anoplocephalidae, providing rich sources for selecting useful molecular markers for ecological and phylogenetic studies. Analyses on mt genome sequences of the four families of cestodes provide novel insights into their phylogenetic relationships. Of couse, more taxon sampling is necessary for future phylogenetic studies of these cestodes using mt genome sequences.

Electronic supplementary material

The online version of this article (doi:10.1186/s13071-015-1172-z) contains supplementary material, which is available to authorized users.

Molecular signatures that are distinctive characteristics of the vertebrates and chordates and supporting a grouping of vertebrates with the tunicates

Members of the phylum Chordata and the subphylum Vertebrata are presently distinguished solely on the basis of morphological characteristics. The relationship of the vertebrates to the two non-vertebrate chordate subphyla is also a subject of debate. Analyses of protein sequences have identified multiple conserved signature indels (CSIs) that are specific for Chordata or for Vertebrata. Five CSIs in 4 important proteins are specific for the Vertebrata, whereas two other CSIs are uniquely found in all sequenced chordate species including Ciona intestinalis and Oikapleura dioica (Tunicates) as well as Branchiostoma floridae (Cephalochordates). The shared presence of these molecular signatures by all vertebrates/chordate species, but in no other animal taxa, strongly indicates that the genetic changes represented by the identified CSIs diagnose monophyletic groups. Two other discovered CSIs are uniquely shared by different vertebrate species and by either one (Ciona intestinalis) or both tunicate (Ciona and Oikapleura) species, but they are not found in Branchiostoma or other animal species. Specific presence of these CSIs in different vertebrates and either one or both tunicate species provides strong independent evidence that the vertebrate species are more closely related to the urochordates (tunicates) than to the cephalochordates.

ddRAD-seq phylogenetics based on nucleotide, indel, and presence-absence polymorphisms: Analyses of two avian genera with contrasting histories

Genotype-by-sequencing (GBS) methods have revolutionized the field of molecular ecology, but their application in molecular phylogenetics remains somewhat limited. In addition, most phylogenetic studies based on large GBS data sets have relied on analyses of concatenated data rather than species tree methods that explicitly account for genealogical stochasticity among loci. We explored the utility of "double-digest" restriction site-associated DNA sequencing (ddRAD-seq) for phylogenetic analyses of the Lagonosticta firefinches (family Estrildidae) and the Vidua brood parasitic finches (family Viduidae). As expected, the number of homologous loci shared among samples was negatively correlated with genetic distance due to the accumulation of restriction site polymorphisms. Nonetheless, for each genus, we obtained data sets of ∼3000 loci shared in common among all samples, including a more distantly related outgroup taxon. For all samples combined, we obtained >1000 homologous loci despite ∼20my divergence between estrildid and parasitic finches. In addition to nucleotide polymorphisms, the ddRAD-seq data yielded large sets of indel and locus presence-absence polymorphisms, all of which had higher consistency indices than mtDNA sequence data in the context of concatenated parsimony analyses. Species tree methods, using individual gene trees or single nucleotide polymorphisms as input, generated results broadly consistent with analyses of concatenated data, particularly for Lagonosticta, which appears to have a well resolved, bifurcating history. Results for Vidua were also generally consistent across methods and data sets, although nodal support and results from different species tree methods were more variable. Lower gene tree congruence in Vidua is likely the result of its unique evolutionary history, which includes rapid speciation by host shift and occasional hybridization and introgression due to incomplete reproductive isolation. We conclude that ddRAD-seq is a cost-effective method for generating robust phylogenetic data sets, particularly for analyses of closely related species and genera.

A phylogenomic and molecular markers based analysis of the phylum Chlamydiae: proposal to divide the class Chlamydiia into two orders, Chlamydiales and Parachlamydiales ord. nov., and emended description of the class Chlamydiia

The phylum Chlamydiae contains nine ecologically and genetically diverse families all placed within a single order. In this work, we have completed a comprehensive comparative analysis of 36 sequenced Chlamydiae genomes in order to identify shared molecular characteristics, namely conserved signature insertions/deletions (CSIs) and conserved signature proteins (CSPs), which can serve as distinguishing characteristics of supra-familial clusters within the phylum Chlamydiae. Our analysis has led to the identification of 32 CSIs which are specific to clusters within the phylum Chlamydiae at various phylogenetic depths. Importantly, 17 CSIs and 98 CSPs were found to be specific for the family Chlamydiaceae while another 3 CSI variants and 15 CSPs were specific for a grouping of the families Criblamydiaceae, Parachlamydiaceae, Simkaniaceae and Waddliaceae. These two clusters were also found to be distinguishable in 16S rRNA based phylogenetic trees, concatenated protein based phylogenetic trees, character compatibility based phylogenetic analyses, and on the basis of 16S rRNA gene sequence identity and average amino acid identity values. On the basis of the identified molecular characteristics, branching in phylogenetic trees, and the genetic distance between the two clusters within the phylum Chlamydiae we propose a division of the class Chlamydiia into two orders: an emended order Chlamydiales, containing the family Chlamydiaceae and the closely related Candidatus family Clavichlamydiaceae, and the novel order Parachlamydiales ord. nov. containing the families Parachlamydiaceae, Simkaniaceae and Waddliaceae and the Candidatus families Criblamydiaceae, Parilichlamydiaceae, Piscichlamydiaceae, and Rhabdochlamydiaceae. We also include a brief discussion of the reunification of the genera Chlamydia and Chlamydophila.

Genome-based taxonomic framework for the class Negativicutes: division of the class Negativicutes into the orders Selenomonadales emend., Acidaminococcales ord. nov. and Veillonellales ord. nov

The class Negativicutes is currently divided into one order and two families on the basis of 16S rRNA gene sequence phylogenies. We report here comprehensive comparative genomic analyses of the sequenced members of the class Negativicutes to demarcate its different evolutionary groups in molecular terms, independently of phylogenetic trees. Our comparative genomic analyses have identified 14 conserved signature indels (CSIs) and 48 conserved signature proteins (CSPs) that either are specific for the entire class or differentiate four main groups within the class. Two CSIs and nine CSPs are shared uniquely by all or most members of the class Negativicutes, distinguishing this class from all other sequenced members of the phylum Firmicutes. Four other CSIs and six CSPs were specific characteristics of the family Acidaminococcaceae, two CSIs and four CSPs were uniquely present in the family Veillonellaceae, six CSIs and eight CSPs were found only in Selenomonas and related genera, and 17 CSPs were identified uniquely in Sporomusa and related genera. Four additional CSPs support a pairing of the groups containing the genera Selenomonas and Sporomusa. We also report detailed phylogenetic analyses for the Negativicutes based on core protein sequences and 16S rRNA gene sequences, which strongly support the four main groups identified by CSIs and by CSPs. Based on the results from different lines of investigation, we propose a division of the class Negativicutes into an emended order Selenomonadales containing the new families Selenomonadaceae fam. nov. and Sporomusaceae fam. nov. and two new orders, Acidaminococcales ord. nov. and Veillonellales ord. nov., respectively containing the families Acidaminococcaceae and Veillonellaceae.

Phylogenetic reconstruction by cross-species chromosome painting and G-banding in four species of Phyllostomini tribe (Chiroptera, Phyllostomidae) in the Brazilian Amazon: an independent evidence for monophyly

The subfamily Phyllostominae comprises taxa with a variety of feeding strategies. From the cytogenetic point of view, Phyllostominae shows different rates of chromosomal evolution between genera, with Phyllostomus hastatus probably retaining the ancestral karyotype for the subfamily. Since chromosomal rearrangements occur rarely in the genome and have great value as phylogenetic markers and in taxonomic characterization, we analyzed three species: Lophostoma silvicola (LSI), Phyllostomus discolor (PDI) and Tonatia saurophila (TSA), representing the tribe Phyllostomini, collected in the Amazon region, by classic and molecular cytogenetic techniques in order to reconstruct the phylogenetic relationships within this tribe. LSA has a karyotype of 2n=34 and FN=60, PDI has 2n=32 and FN=60 and TSA has 2n=16 and FN=20. Comparative analysis using G-banding and chromosome painting show that the karyotypic complement of TSA is highly rearranged relative to LSI and PHA, while LSI, PHA and PDI have similar karyotypes, differing by only three chromosome pairs. Nearly all chromosomes of PDI and PHA were conserved in toto, except for chromosome 15 that was changed by a pericentric inversion. A strongly supported phylogeny (bootstrap=100 and Bremer=10 steps), confirms the monophyly of Phyllostomini. In agreement with molecular topologies, TSA was in the basal position, while PHA and LSI formed sister taxa. A few ancestral syntenies are conserved without rearrangements and most associations are autapomorphic traits for Tonatia or plesiomorphic for the three genera analyzed here. The karyotype of TSA is highly derived in relation to that of other phyllostomid bats, differing from the supposed ancestral karyotype of Phyllostomidae by multiple rearrangements. Phylogenies based on chromosomal data are independent evidence for the monophyly of tribe Phyllostomini as determined by molecular topologies and provide additional support for the paraphyly of the genus Tonatia by the exclusion of the genus Lophostoma.

Phylogenomic and molecular demarcation of the core members of the polyphyletic pasteurellaceae genera actinobacillus, haemophilus, and pasteurella

The genera Actinobacillus, Haemophilus, and Pasteurella exhibit extensive polyphyletic branching in phylogenetic trees and do not represent coherent clusters of species. In this study, we have utilized molecular signatures identified through comparative genomic analyses in conjunction with genome based and multilocus sequence based phylogenetic analyses to clarify the phylogenetic and taxonomic boundary of these genera. We have identified large clusters of Actinobacillus, Haemophilus, and Pasteurella species which represent the “sensu stricto” members of these genera. We have identified 3, 7, and 6 conserved signature indels (CSIs), which are specifically shared by sensu stricto members of Actinobacillus, Haemophilus, and Pasteurella, respectively. We have also identified two different sets of CSIs that are unique characteristics of the pathogen containing genera Aggregatibacter and Mannheimia, respectively. It is now possible to demarcate the genera Actinobacillus sensu stricto, Haemophilus sensu stricto, and Pasteurella sensu stricto on the basis of discrete molecular signatures. The other members of the genera Actinobacillus, Haemophilus, and Pasteurella that do not fall within the “sensu stricto” clades and do not contain these molecular signatures should be reclassified as other genera. The CSIs identified here also provide useful diagnostic targets for the identification of current and novel members of the indicated genera.

The phylogenetic position of ctenophores and the origin(s) of nervous systems

Ctenophores have traditionally been treated as eumetazoans, but some recent whole genome studies have revived the idea that they are, rather, the sister group to all other metazoans. This deep branching position implies either that nervous systems have evolved twice, in Ctenophora and in Eumetazoa, or that an ancestral metazoan nervous system has been lost in sponges and placozoans. We caution, however, that phylogenetic-tree construction artifacts may have placed ctenophores too deep in the metazoan tree. We discuss nervous system origins under these alternative phylogenies and in light of comparative data of ctenophore and eumetazoan nervous systems. We argue that characters like neuropeptide signaling, ciliary photoreceptors, gap junctions and presynaptic molecules are consistent with a shared ancestry of nervous systems. However, if ctenophores are the sister group to all other metazoans, this ancestral nervous system was likely very simple. Further studies are needed to resolve the deep phylogeny of metazoans and to have a better understanding of the early steps of nervous system evolution.

Molecular signatures and phylogenomic analysis of the genus Burkholderia: proposal for division of this genus into the emended genus Burkholderia containing pathogenic organisms and a new genus Paraburkholderia gen. nov. harboring environmental species

The genus Burkholderia contains large number of diverse species which include many clinically important organisms, phytopathogens, as well as environmental species. However, currently, there is a paucity of biochemical or molecular characteristics which can reliably distinguish different groups of Burkholderia species. We report here the results of detailed phylogenetic and comparative genomic analyses of 45 sequenced species of the genus Burkholderia. In phylogenetic trees based upon concatenated sequences for 21 conserved proteins as well as 16S rRNA gene sequence based trees, members of the genus Burkholderia grouped into two major clades. Within these main clades a number of smaller clades including those corresponding to the clinically important Burkholderia cepacia complex (BCC) and the Burkholderia pseudomallei groups were also clearly distinguished. Our comparative analysis of protein sequences from Burkholderia spp. has identified 42 highly specific molecular markers in the form of conserved sequence indels (CSIs) that are uniquely found in a number of well-defined groups of Burkholderia spp. Six of these CSIs are specific for a group of Burkholderia spp. (referred to as Clade I in this work) which contains all clinically relevant members of the genus (viz. the BCC and the B. pseudomallei group) as well as the phytopathogenic Burkholderia spp. The second main clade (Clade II), which is composed of environmental Burkholderia species, is also distinguished by 2 identified CSIs that are specific for this group. Additionally, our work has also identified multiple CSIs that serve to clearly demarcate a number of smaller groups of Burkholderia spp. including 3 CSIs that are specific for the B. cepacia complex, 4 CSIs that are uniquely found in the B. pseudomallei group, 5 CSIs that are specific for the phytopathogenic Burkholderia spp. and 22 other CSI that distinguish two groups within Clade II. The described molecular markers provide highly specific means for the demarcation of different groups of Burkholderia spp. and they also offer novel and useful targets for the development of diagnostic assays for the clinically important members of the BCC or the pseudomallei groups. Based upon the results of phylogenetic analyses, the identified CSIs and the pathogenicity profile of Burkholderia species, we are proposing a division of the genus Burkholderia into two genera. In this new proposal, the emended genus Burkholderia will correspond to the Clade I and it will contain only the clinically relevant and phytopathogenic Burkholderia species. All other Burkholderia spp., which are primarily environmental, will be transferred to a new genus Paraburkholderia gen. nov.

Nanger, Eudorcas, Gazella, and Antilope form a well-supported chromosomal clade within Antilopini (Bovidae, Cetartiodactyla)

The evolutionary clade comprising Nanger, Eudorcas, Gazella, and Antilope, defined by an X;BTA5 translocation, is noteworthy for the many autosomal Robertsonian fusions that have driven the chromosome number variation from 2n = 30 observed in Antilope cervicapra, to the 2n = 58 in present Eudorcas thomsoni and Eudorcas rufifrons. This work reports the phylogenetic relationships within the Antilopini using comprehensive cytogenetic data from A. cervicapra, Gazella leptoceros, Nanger dama ruficollis, and E. thomsoni together with corrected karyotypic data from an additional nine species previously reported in the literature. Fluorescence in situ hybridization using BAC and microdissected cattle painting probes, in conjunction with differential staining techniques, provide the following: (i) a detailed analysis of the E. thomsoni chromosomes, (ii) the identification and fine-scale analysis the BTA3 orthologue in species of Antilopini, and (iii) the location of the pseudoautosomal regions on sex chromosomes of the four species. Our phylogenetic analysis of the chromosomal data supports monophyly of Nanger and Eudorcas and suggests an affiliation between A. cervicapra and some of the Gazella species. This renders Gazella paraphyletic and emphasizes a closer relationship between Antilope and Gazella than what has previously been considered.

Indel reliability in indel-based phylogenetic inference

It is often assumed that it is unlikely that the same insertion or deletion (indel) event occurred at the same position in two independent evolutionary lineages, and thus, indel-based inference of phylogeny should be less subject to homoplasy compared with standard inference which is based on substitution events. Indeed, indels were successfully used to solve debated evolutionary relationships among various taxonomical groups. However, indels are never directly observed but rather inferred from the alignment and thus indel-based inference may be sensitive to alignment errors. It is hypothesized that phylogenetic reconstruction would be more accurate if it relied only on a subset of reliable indels instead of the entire indel data. Here, we developed a method to quantify the reliability of indel characters by measuring how often they appear in a set of alternative multiple sequence alignments. Our approach is based on the assumption that indels that are consistently present in most alternative alignments are more reliable compared with indels that appear only in a small subset of these alignments. Using simulated and empirical data, we studied the impact of filtering and weighting indels by their reliability scores on the accuracy of indel-based phylogenetic reconstruction. The new method is available as a web-server at http://guidance.tau.ac.il/RELINDEL/.

Phylogenetic analysis and molecular signatures defining a monophyletic clade of heterocystous cyanobacteria and identifying its closest relatives

Detailed phylogenetic and comparative genomic analyses are reported on 140 genome sequenced cyanobacteria with the main focus on the heterocyst-differentiating cyanobacteria. In a phylogenetic tree for cyanobacteria based upon concatenated sequences for 32 conserved proteins, the available cyanobacteria formed 8-9 strongly supported clades at the highest level, which may correspond to the higher taxonomic clades of this phylum. One of these clades contained all heterocystous cyanobacteria; within this clade, the members exhibiting either true (Nostocales) or false (Stigonematales) branching of filaments were intermixed indicating that the division of the heterocysts-forming cyanobacteria into these two groups is not supported by phylogenetic considerations. However, in both the protein tree as well as in the 16S rRNA gene tree, the akinete-forming heterocystous cyanobacteria formed a distinct clade. Within this clade, the members which differentiate into hormogonia or those which lack this ability were also separated into distinct groups. A novel molecular signature identified in this work that is uniquely shared by the akinete-forming heterocystous cyanobacteria provides further evidence that the members of this group are specifically related and they shared a common ancestor exclusive of the other cyanobacteria. Detailed comparative analyses on protein sequences from the genomes of heterocystous cyanobacteria reported here have also identified eight conserved signature indels (CSIs) in proteins involved in a broad range of functions, and three conserved signature proteins, that are either uniquely or mainly found in all heterocysts-forming cyanobacteria, but generally not found in other cyanobacteria. These molecular markers provide novel means for the identification of heterocystous cyanobacteria, and they provide evidence of their monophyletic origin. Additionally, this work has also identified seven CSIs in other proteins which in addition to the heterocystous cyanobacteria are uniquely shared by two smaller clades of cyanobacteria, which form the successive outgroups of the clade comprising of the heterocystous cyanobacteria in the protein trees. Based upon their close relationship to the heterocystous cyanobacteria, the members of these clades are indicated to be the closest relatives of the heterocysts-forming cyanobacteria.

A retroposon-based view on the temporal differentiation of sex chromosomes

Retroposon presence/absence patterns in orthologous genomic loci are known to be strong and almost homoplasy-free phylogenetic markers of common ancestry. This is evidenced by the comprehensive reconstruction of various species trees of vertebrate lineages in recent years, as well as the inference of the evolution of genes via retroposon-based gene trees of paralogous genes. Recently, it has been shown that retroposon markers are also suitable for the inference of differentiation events of gametologous genes, i.e., homologous genes on opposite sex chromosomes. This is because sex chromosomes evolved via stepwise cessation of recombination, making the presence or absence of a particular retroposon insertion among the two different gametologs in more or less closely related species a clear-cut indicator of the timing of differentiation events. Here, I examine the advantages and current limitations of this novel perspective for understanding avian sex chromosome evolution, compare the retroposon-based and sequence-based insights into gametolog differentiation and show that retroposons promise to be equally applicable to other sex chromosomal systems, such as the human X and Y chromosomes.

Discriminating between deleterious and neutral non-frameshifting indels based on protein interaction networks and hybrid properties

More than ten thousand coding variants are contained in each human genome; however, our knowledge of the way genetic variants underlie phenotypic differences is far from complete. Small insertions and deletions (indels) are one of the most common types of human genetic variants, and indels play a significant role in human inherited disease. To date, we still lack a comprehensive understanding of how indels cause diseases. Therefore, identification and analysis of such deleterious variants is a key challenge and has been of great interest in the current research in genome biology. Increasing numbers of computational methods have been developed for discriminating between deleterious indels and neutral indels. However, most of the existing methods are based on traditional sequential or structural features, which cannot completely explain the association between indels and the resulting induced inherited disease. In this study, we establish a novel method to predict deleterious non-frameshifting indels based on features extracted from both protein interaction networks and traditional hybrid properties. Each indel was coded by 1,246 features. Using the maximum relevance minimum redundancy method and the incremental feature selection method, we obtained an optimal feature set containing 42 features, of which 21 features were derived from protein interaction networks. Based on the optimal feature set, an 88 % accuracy and a 0.76 MCC value were achieved by a Random Forest as evaluated by the Jackknife cross-validation test. This method outperformed existing methods of predicting deleterious indels, and can be applied in practice for deleterious non-frameshifting indel predictions in genome research. The analysis of the optimal features selected in the model revealed that network interactions play more important roles and could be informative for better illustrating an indel's function and disease associations than traditional sequential or structural features. These results could shed some light on the genetic basis of human genetic variations and human inherited diseases.

Complex patterns of gene fission in the eukaryotic folate biosynthesis pathway

Shared derived genomic characters can be useful for polarizing phylogenetic relationships, for example, gene fusions have been used to identify deep-branching relationships in the eukaryotes. Here, we report the evolutionary analysis of a three-gene fusion of folB, folK, and folP, which encode enzymes that catalyze consecutive steps in de novo folate biosynthesis. The folK-folP fusion was found across the eukaryotes and a sparse collection of prokaryotes. This suggests an ancient derivation with a number of gene losses in the eukaryotes potentially as a consequence of adaptation to heterotrophic lifestyles. In contrast, the folB-folK-folP gene is specific to a mosaic collection of Amorphea taxa (a group encompassing: Amoebozoa, Apusomonadida, Breviatea, and Opisthokonta). Next, we investigated the stability of this character. We identified numerous gene losses and a total of nine gene fission events, either by break up of an open reading frame (four events identified) or loss of a component domain (five events identified). This indicates that this three gene fusion is highly labile. These data are consistent with a growing body of data indicating gene fission events occur at high relative rates. Accounting for these sources of homoplasy, our data suggest that the folB-folK-folP gene fusion was present in the last common ancestor of Amoebozoa and Opisthokonta but absent in the Metazoa including the human genome. Comparative genomic data of these genes provides an important resource for designing therapeutic strategies targeting the de novo folate biosynthesis pathway of a variety of eukaryotic pathogens such as Acanthamoeba castellanii.

A critical appraisal of the use of microRNA data in phylogenetics

Recent progress in resolving the tree of life continues to expose relationships that resist resolution, which drives the search for novel sources of information to solve these difficult phylogenetic problems. A recent example, the presence and absence of microRNA families, has been vigorously promoted as an ideal source of phylogenetic data and has been applied to several perennial phylogenetic problems. The utility of such data for phylogenetic inference hinges critically both on developing stochastic models that provide a reasonable description of the process that give rise to these data, and also on the careful validation of those models in real inference scenarios. Remarkably, however, the statistical behavior and phylogenetic utility of microRNA data have not yet been rigorously characterized. Here we explore the behavior and performance of microRNA presence/absence data under a variety of evolutionary models and reexamine datasets from several previous studies. We find that highly heterogeneous rates of microRNA gain and loss, pervasive secondary loss, and sampling error collectively render microRNA-based inference of phylogeny difficult. Moreover, our reanalyses fundamentally alter the conclusions for four of the five studies that we reexamined. Our results indicate that the capacity of miRNA data to resolve the tree of life has been overstated, and we urge caution in their application and interpretation.

A transcriptome approach to ecdysozoan phylogeny

The monophyly of Ecdysozoa, which comprise molting phyla, has received strong support from several lines of evidence. However, the internal relationships of Ecdysozoa are still contended. We generated expressed sequence tags from a priapulid (penis worm), a kinorhynch (mud dragon), a tardigrade (water bear) and five chelicerate taxa by 454 transcriptome sequencing. A multigene alignment was assembled from 63 taxa, which comprised after matrix optimization 24,249 amino acid positions with high data density (2.6% gaps, 19.1% missing data). Phylogenetic analyses employing various models support the monophyly of Ecdysozoa. A clade combining Priapulida and Kinorhyncha (i.e. Scalidophora) was recovered as the earliest branch among Ecdysozoa. We conclude that Cycloneuralia, a taxon erected to combine Priapulida, Kinorhyncha and Nematoda (and others), are paraphyletic. Rather Arthropoda (including Onychophora) are allied with Nematoda and Tardigrada. Within Arthropoda, we found strong support for most clades, including monophyletic Mandibulata and Pancrustacea. The phylogeny within the Euchelicerata remained largely unresolved. There is conflicting evidence on the position of tardigrades: While Bayesian and maximum likelihood analyses of only slowly evolving genes recovered Tardigrada as a sister group to Arthropoda, analyses of the full data set, and of subsets containing genes evolving at fast and intermediate rates identified a clade of Tardigrada and Nematoda. Notably, the latter topology is also supported by the analyses of indel patterns.

Characterization of the complete mitochondrial genome of Chilo auricilius and comparison with three other rice stem borers

The mitogenome of Chilo auricilius (Lepidoptera: Pyraloidea: Crambidae) was a circular molecule made up of 15,367 bp. Sesamia inferens, Chilo suppressalis, Tryporyza incertulas, and C. auricilius, are closely related, well known rice stem borers that are widely distributed in the main rice-growing regions of China. The gene order and orientation of all four stem borers were similar to that of other insect mitogenomes. Among the four stem borers, all AT contents were below 83%, while all AT contents of tRNA genes were above 80%. The genomes were compact, with only 121-257 bp of non-coding intergenic spacer. There are 56 or 62-bp overlapping nucleotides in Crambidae moths, but were only 25-bp overlapping nucleotides in the noctuid moth S. inferens. There was a conserved motif 'ATACTAAA' between trnS2 (UCN) and nad1 in Crambidae moths, but this same region was 'ATCATA' in the noctuid S. inferens. And there was a 6-bp motif 'ATGATAA' of overlapping nucleotides, which was conserved in Lepidoptera, and a 14-bp motif 'TAAGCTATTTAAAT' conserved in the three Crambidae moths (C. suppressalis, C. auricilius and T. incertulas), but not in the noctuid. Finally, there were no stem-and-loop structures in the two Chilo moths.

Origin and evolution of glutamyl-prolyl tRNA synthetase WHEP domains reveal evolutionary relationships within Holozoa

Repeated domains in proteins that have undergone duplication or loss, and sequence divergence, are especially informative about phylogenetic relationships. We have exploited divergent repeats of the highly structured, 50-amino acid WHEP domains that join the catalytic subunits of bifunctional glutamyl-prolyl tRNA synthetase (EPRS) as a sequence-informed repeat (SIR) to trace the origin and evolution of EPRS in holozoa. EPRS is the only fused tRNA synthetase, with two distinct aminoacylation activities, and a non-canonical translation regulatory function mediated by the WHEP domains in the linker. Investigating the duplications, deletions and divergence of WHEP domains, we traced the bifunctional EPRS to choanozoans and identified the fusion event leading to its origin at the divergence of ichthyosporea and emergence of filozoa nearly a billion years ago. Distribution of WHEP domains from a single species in two or more distinct clades suggested common descent, allowing the identification of linking organisms. The discrete assortment of choanoflagellate WHEP domains with choanozoan domains as well as with those in metazoans supported the phylogenetic position of choanoflagellates as the closest sister group to metazoans. Analysis of clustering and assortment of WHEP domains provided unexpected insights into phylogenetic relationships amongst holozoan taxa. Furthermore, observed gaps in the transition between WHEP domain groupings in distant taxa allowed the prediction of undiscovered or extinct evolutionary intermediates. Analysis based on SIR domains can provide a phylogenetic counterpart to palaentological approaches of discovering “missing links” in the tree of life.

Phylogeny and molecular signatures for the phylum Fusobacteria and its distinct subclades

The members of the phylum Fusobacteria and its two families, Fusobacteriaceae and Leptotrichiaceae, are distinguished at present mainly on the basis of their branching in the 16S rRNA gene trees and analysis of the internal transcribed spacer sequences in the 16S-23S rDNA. However, no biochemical or molecular characteristics are known that are uniquely shared by all of most members of these groups of bacteria. We report here detailed phylogenetic and comparative analyses on 45 sequenced Fusobacteria genomes to examine their evolutionary relationships and to identify molecular markers that are specific for the members of this phylum. In phylogenetic trees based on 16S rRNA gene sequences or concatenated sequences for 17 conserved proteins, members of the families Fusobacteriaceae and Leptotrichiaceae formed strongly supported clades and were clearly distinguished. In these trees, the species from the genus Fusobacterium also formed a number of well-supported clades. In parallel, comparative analyses on Fusobacteria genomes have identified 44 conserved signature indels (CSIs) in proteins involved in a broad range of functions that are either specific for the phylum Fusobacteria or a number of distinct subclades within this phylum. Seven of these CSIs in important proteins are uniquely present in the protein homologs of all sequenced Fusobacteria and they provide potential molecular markers for this phylum. Six and three other CSIs in other protein sequences are specific for members of the families Fusobacteriaceae and Leptotrichiaceae, respectively, and they provide novel molecular means for distinguishing members of these two families. Fourteen additional CSIs in different proteins, which are specific for either members of the genera Fusobacterium or Leptotrichia, or a number of other well-supported clades of Fusobacteria at multiple phylogenetic levels, provide molecular markers for these groups and information regarding the evolutionary relationships among the members of this phylum. Lastly, the present work has also identified 14 CSIs in divergent proteins that are specific for three specific subclades of Fusobacterium species, which are also indicated to be distinct by phylogenetic analyses. The members of these three Fusobacterium subclades also differ significantly from each other in their whole genome average nucleotide identities values, suggesting that they are possible candidates for recognition as different genera. The molecular markers reported here provide novel means for the identification of members of the phylum Fusobacteria and for their classification.

A description of the complete mitochondrial genomes of Amphiporus formidabilis, Prosadenoporus spectaculum and Nipponnemertes punctatula (Nemertea: Hoplonemertea: Monostilifera)

We sequenced the complete mitochondrial genomes (mitogenomes) of three Hoplonemertea species, Amphiporus formidabilis, Prosadenoporus spectaculum and Nipponnemertes punctatula, which are 14,616, 14,655 and 15,354 bp in length, respectively. Each of the three circular mitogenomes consists of 37 typical genes and some non-coding regions. The nucleotide composition of the coding strand is biased toward T, almost a half of total nucleotides in these mitogenomes. There are many poly-T tracts across these mitogenomes, which exhibit T-number variation within different clones of protein-coding genes, mainly resulting from false PCR amplification. The major non-coding regions have tandem repeat motifs and hairpin-like structures that may be associated with the initiation of replication or transcription. Data published to date for nemerteans show that Palaeonemertea species usually bear the largest mitogenomes, while representatives in the more recently derived Distromatonemertea clade bear the smallest ones; and that the gene arrangement of mitogenomes seems to be variable within the phylum Nemertea, but stable within either of Heteronemertea and Hoplonemertea.

Physella acuta: atypical mitochondrial gene order among panpulmonates (Gastropoda)

Mitochondrial (mt) sequences are frequently used for phylogenetic reconstruction and for identification of species of molluscs. This study expands the phylogenetic range of Hygrophila (Panpulmonata) for which such sequence data are available by characterizing the full mt genome of the invasive freshwater snail Physella acuta (Physidae). The mt genome sequences of two P. acuta isolates from Stubblefield Lake, New Mexico, USA, differed in length (14,490 vs 14,314 bp) and showed 11.49% sequence divergence, whereas ITS1 and ITS2 sequences from the nuclear genome differed by 1.75%. The mt gene order of P. acuta (cox1, P, nad6, nad5, nad1, D, F, cox2, Y, W, nad4L, C, Q, atp6, R, E, rrnS, M, T, cox3, I, nad2, K, V, rrnL, L1, A, cytb, G, H, L2, atp8, N, nad2, S1, S2, nad4) differs considerably from the relatively conserved gene order within Panpulmonata. Phylogenetic trees show that the 13 protein-encoding mt gene sequences (equivalent codons) of P. acuta group according to gastropod phylogeny, yet branch lengths and dN/dS ratios for P. acuta indicate elevated amino acid substitutions relative to other gastropods. This study indicates that mt sequences of P. acuta are phylogenetically informative despite a considerable intraspecific divergence and the atypical gene order in its mt genome.

A phylogenomic and molecular marker based proposal for the division of the genus Borrelia into two genera: the emended genus Borrelia containing only the members of the relapsing fever Borrelia, and the genus Borreliella gen. nov. containing the members of the Lyme disease Borrelia (Borrelia burgdorferi sensu lato complex)

The genus Borrelia contains two groups of organisms: the causative agents of Lyme disease and their relatives and the causative agents of relapsing fever and their relatives. These two groups are morphologically indistinguishable and are difficult to distinguish biochemically. In this work, we have carried out detailed comparative genomic analyses on protein sequences from 38 Borrelia genomes to identify molecular markers in the forms of conserved signature inserts/deletions (CSIs) that are specifically found in the Borrelia homologues, and conserved signature proteins (CSPs) which are uniquely present in Borrelia species. Our analyses have identified 31 CSIs and 82 CSPs that are uniquely shared by all sequenced Borrelia species, providing molecular markers for this group of organisms. In addition, our work has identified 7 CSIs and 21 CSPs which are uniquely found in the Lyme disease Borrelia species and eight CSIs and four CSPs that are specific for members of the relapsing fever Borrelia group. Additionally, 38 other CSIs, in proteins which are uniquely found in Borrelia species, also distinguish these two groups of Borrelia. The identified CSIs and CSPs provide novel and highly specific molecular markers for identification and distinguishing between the Lyme disease Borrelia and the relapsing fever Borrelia species. We also report the results of average nucleotide identity (ANI) analysis on Borrelia genomes and phylogenetic analysis for these species based upon 16S rRNA sequences and concatenated sequences for 25 conserved proteins. These analyses also support the distinctness of the two Borrelia clades. On the basis of the identified molecular markers, the results from ANI and phylogenetic studies, and the distinct pathogenicity profiles and arthropod vectors used by different Borrelia spp. for their transmission, we are proposing a division of the genus Borrelia into two separate genera: an emended genus Borrelia, containing the causative agents of relapsing fever and a novel genus, Borreliella gen. nov., containing the causative agents of Lyme disease.

Analyses of charophyte chloroplast genomes help characterize the ancestral chloroplast genome of land plants

Despite the significance of the relationships between embryophytes and their charophyte algal ancestors in deciphering the origin and evolutionary success of land plants, few chloroplast genomes of the charophyte algae have been reconstructed to date. Here, we present new data for three chloroplast genomes of the freshwater charophytes Klebsormidium flaccidum (Klebsormidiophyceae), Mesotaenium endlicherianum (Zygnematophyceae), and Roya anglica (Zygnematophyceae). The chloroplast genome of Klebsormidium has a quadripartite organization with exceptionally large inverted repeat (IR) regions and, uniquely among streptophytes, has lost the rrn5 and rrn4.5 genes from the ribosomal RNA (rRNA) gene cluster operon. The chloroplast genome of Roya differs from other zygnematophycean chloroplasts, including the newly sequenced Mesotaenium, by having a quadripartite structure that is typical of other streptophytes. On the basis of the improbability of the novel gain of IR regions, we infer that the quadripartite structure has likely been lost independently in at least three zygnematophycean lineages, although the absence of the usual rRNA operonic synteny in the IR regions of Roya may indicate their de novo origin. Significantly, all zygnematophycean chloroplast genomes have undergone substantial genomic rearrangement, which may be the result of ancient retroelement activity evidenced by the presence of integrase-like and reverse transcriptase-like elements in the Roya chloroplast genome. Our results corroborate the close phylogenetic relationship between Zygnematophyceae and land plants and identify 89 protein-coding genes and 22 introns present in the chloroplast genome at the time of the evolutionary transition of plants to land, all of which can be found in the chloroplast genomes of extant charophytes.

Evaluating rare amino acid substitutions (RGC_CAMs) in a yeast model clade

When inferring phylogenetic relationships, not all sites in a sequence alignment are equally informative. One recently proposed approach that takes advantage of this inequality relies on sites that contain amino acids whose replacement requires multiple substitutions. Identifying these so-called RGC_CAM substitutions (after Rare Genomic Changes as Conserved Amino acids-Multiple substitutions) requires that, first, at any given site in the amino acid sequence alignment, there must be a minimum of two different amino acids; second, each amino acid must be present in at least two taxa; and third, the amino acids must require a minimum of two nucleotide substitutions to replace each other. Although theory suggests that RGC_CAM substitutions are expected to be rare and less likely to be homoplastic, the informativeness of RGC_CAM substitutions has not been extensively evaluated in biological data sets. We investigated the quality of RGC_CAM substitutions by examining their degree of homoplasy and internode certainty in nearly 2.7 million aligned amino acid sites from 5,261 proteins from five species belonging to the yeast Saccharomyces sensu stricto clade whose phylogeny is well-established. We identified 2,647 sites containing RGC_CAM substitutions, a number that contrasts sharply with the 100,887 sites containing RGC_non-CAM substitutions (i.e., changes between amino acids that require only a single nucleotide substitution). We found that RGC_CAM substitutions had significantly lower homoplasy than RGC_non-CAM ones; specifically RGC_CAM substitutions showed a per-site average homoplasy index of 0.100, whereas RGC_non-CAM substitutions had a homoplasy index of 0.215. Internode certainty values were also higher for sites containing RGC_CAM substitutions than for RGC_non-CAM ones. These results suggest that RGC_CAM substitutions possess a strong phylogenetic signal and are useful markers for phylogenetic inference despite their rarity.

Novel information theory-based measures for quantifying incongruence among phylogenetic trees

Phylogenies inferred from different data matrices often conflict with each other necessitating the development of measures that quantify this incongruence. Here, we introduce novel measures that use information theory to quantify the degree of conflict or incongruence among all nontrivial bipartitions present in a set of trees. The first measure, internode certainty (IC), calculates the degree of certainty for a given internode by considering the frequency of the bipartition defined by the internode (internal branch) in a given set of trees jointly with that of the most prevalent conflicting bipartition in the same tree set. The second measure, IC All (ICA), calculates the degree of certainty for a given internode by considering the frequency of the bipartition defined by the internode in a given set of trees in conjunction with that of all conflicting bipartitions in the same underlying tree set. Finally, the tree certainty (TC) and TC All (TCA) measures are the sum of IC and ICA values across all internodes of a phylogeny, respectively. IC, ICA, TC, and TCA can be calculated from different types of data that contain nontrivial bipartitions, including from bootstrap replicate trees to gene trees or individual characters. Given a set of phylogenetic trees, the IC and ICA values of a given internode reflect its specific degree of incongruence, and the TC and TCA values describe the global degree of incongruence between trees in the set. All four measures are implemented and freely available in version 8.0.0 and subsequent versions of the widely used program RAxML.

Conserved signature indels and signature proteins as novel tools for understanding microbial phylogeny and systematics: identification of molecular signatures that are specific for the phytopathogenic genera Dickeya, Pectobacterium and Brenneria

Genome sequences are enabling applications of different approaches to more clearly understand microbial phylogeny and systematics. Two of these approaches involve identification of conserved signature indels (CSIs) and conserved signature proteins (CSPs) that are specific for different lineages. These molecular markers provide novel and more definitive means for demarcation of prokaryotic taxa and for identification of species from these groups. Genome sequences are also enabling determination of phylogenetic relationships among species based upon sequences for multiple proteins. In this work, we have used all of these approaches for studying the phytopathogenic bacteria belonging to the genera Dickeya , Pectobacterium and Brenneria . Members of these genera, which cause numerous diseases in important food crops and ornamental plants, are presently distinguished mainly on the basis of their branching in phylogenetic trees. No biochemical or molecular characteristic is known that is uniquely shared by species from these genera. Hence, detailed studies using the above approaches were carried out on proteins from the genomes of these bacteria to identify molecular markers that are specific for them. In phylogenetic trees based upon concatenated sequences for 23 conserved proteins, members of the genera Dickeya , Pectobacterium and Brenneria formed a strongly supported clade within the other Enterobacteriales . Comparative analysis of protein sequences from the Dickeya , Pectobacterium and Brenneria genomes has identified 10 CSIs and five CSPs that are either uniquely or largely found in all genome-sequenced species from these genera, but not present in any other bacteria in the database. In addition, our analyses have identified 10 CSIs and 17 CSPs that are specifically present in either all or most sequenced Dickeya species/strains, and six CSIs and 19 CSPs that are uniquely found in the sequenced Pectobacterium genomes. Finally, our analysis also identified three CSIs and one CSP that are specifically shared by members of the genera Pectobacterium and Brenneria , but absent in species of the genus Dickeya , indicating that the former two genera shared a common ancestor exclusive of Dickeya . The identified CSIs and CSPs provide novel tools for identification of members of the genera Dickeya and Pectobacterium and for delimiting these taxa in molecular terms. Descriptions of the genera Dickeya and Pectobacterium have been revised to provide information for these molecular markers. Biochemical studies on these CSIs and CSPs, which are specific for these genera, may lead to discovery of novel properties that are unique to these bacteria and which could be targeted to develop antibacterial agents that are specific for these plant-pathogenic bacteria.

Comparison of RNA isolation methods from insect larvae

Isolating RNA from insects is becoming increasingly important in molecular entomology. Four methods including three commercial kits RNeasy Mini Kit (Qiagen), SV Total RNA isolation system (Promega), TRIzol reagent (Invitrogen), and a cetyl trimethylammonium bromide (CTAB)-based method were compared regarding their ability to isolate RNA from whole-body larvae of Thaumatotibia leucotreta (Meyrick), Thanatophilus micans (F.), Plutella xylostella (L.), and Tenebrio molitor (L.). A difference was observed among the four methods regarding RNA quality but not quantity. However, RNA quality and quantity obtained was not dependent on the insect species. The CTAB-based method produced low-quality RNA and the Trizol reagent produced partially degraded RNA, whereas the RNeasy Mini Kit and SV Total RNA isolation system produced RNA of consistently high quality. However, after reverse transcription to cDNA, RNA produced using all four extraction methods could be used to successfully amplify a 708 bp fragment of the cytochrome oxidase I gene. Of the four methods, the SV Total RNA isolation system showed the least amount of DNA contamination with the highest RNA integrity number and is thus recommended for stringent applications where high-quality RNA is required. This is the first comparison of RNA isolation methods among different insect species and the first to compare RNA isolation methods in insects in the last 20 years.

Tandem duplications in the C-terminal domain of the mesotocin receptor exclusively identified among East Eurasian thrushes

Mesotocin is a neurohypophyseal hormone found in some non-mammalian vertebrates, including birds, reptiles, and amphibians. In this study, we identified and characterized 18-amino acid duplications in the C-terminal domain of the mesotocin receptor (MTR), specifically found in Turdus thrushes (Aves: Passeriforms: Turdidae). These duplicated elements are located in the distal part of the C-terminal tails of MTR and consist of amino acids that are highly conserved among major vertebrates. Intraspecific polymorphisms in a variable number of tandem duplications are commonly found in East Eurasian Turdus, but not in any other genus of Turdidae. Moreover, the genus Turdus can be further classified into 2 groups according to the presence or absence of a 3-amino acid deletion just adjacent to the putative palmitoylation site in the cytoplasmic C-terminal tail. The phylogeny presented here strongly supports the conspecific group of 4 East Eurasian thrushes (Turdus pallidus, T. chrysolaus, T. obscurus, and T. celaenops). Our findings, therefore, provide a new synapomorphy that can be used for phylogenetic assumptions and shed a light on the history of diversification within Eurasian Turdus clades.

Molecular signatures for the phylum (class) Thermotogae and a proposal for its division into three orders (Thermotogales, Kosmotogales ord. nov. and Petrotogales ord. nov.) containing four families (Thermotogaceae, Fervidobacteriaceae fam. nov., Kosmotogaceae fam. nov. and Petrotogaceae fam. nov.) and a new genus Pseudothermotoga gen. nov. with five new combinations

All species from the phylum Thermotogae, class Thermotogae, are currently part of a single family, Thermotogaceae. Using genomic data from 17 Thermotogae species, detailed phylogenetic and comparative genomic analyses were carried out to understand their evolutionary relationships and identify molecular markers that are indicative of species relationships within the phylum. In the 16S rRNA gene tree and phylogenetic trees based upon two different large sets of proteins, members of the phylum Thermotogae formed a number of well-resolved clades. Character compatibility analysis on the protein sequence data also recovered a single largest clique that exhibited similar topology to the protein trees and where all nodes were supported by multiple compatible characters. Comparative genomic analyses have identified 85 molecular markers, in the form of conserved signature indels (CSIs), which are specific for different observed clades of Thermotogae at multiple phylogenetic depths. Eleven of these CSIs were specific for the phylum Thermotogae whereas nine others supported a clade comprising of the genera Thermotoga, Thermosipho and Fervidobacterium. Ten other CSIs provided evidence that the genera Thermosipho and Fervidobacterium shared a common ancestor exclusive of the other Thermotogae and four and eight CSIs in other proteins were specific for the genera Thermosipho and Fervidobacterium, respectively. Two other deep branching clades, one consisting of the genera Kosmotoga and Mesotoga and the other comprising of the genera Petrotoga and Marinitoga, were also supported by multiple CSIs. Based upon the consistent branching of the Thermotogae species using different phylogenetic approaches, and numerous identified CSIs supporting the distinctness of different clades, it is proposed that the class Thermotogae should be divided into three orders (Thermotogales, Kosmotogales ord. nov. and Petrotogales ord. nov.) containing four families (Thermotogaceae, Fervidobacteriaceae fam. nov., Kosmotogaceae fam. nov. and Petrotogaceae fam. nov.). Additionally, the results of our phylogenetic/compatibility studies along with the species distribution patterns of 22 identified CSIs, provide compelling evidence that the current genus Thermotoga is comprised of two evolutionary distinct groups and that it should be divided into two genera. It is proposed that the emended genus Thermotoga should retain only the species Thermotoga maritima, Tt. neapolitana, Tt. petrophila, Tt. naphthophila, Thermotoga sp. EMP, Thermotoga sp. A7A and Thermotoga sp. RQ2 while the other Thermotoga species (viz. Tt. lettingae, Tt. thermarum, Tt. elfii, Tt. subterranean and Tt. hypogea) be transferred to a new genus, Pseudothermotoga gen. nov.

Insect mitochondrial genomics: implications for evolution and phylogeny

The mitochondrial (mt) genome is, to date, the most extensively studied genomic system in insects, outnumbering nuclear genomes tenfold and representing all orders versus very few. Phylogenomic analysis methods have been tested extensively, identifying compositional bias and rate variation, both within and between lineages, as the principal issues confronting accurate analyses. Major studies at both inter- and intraordinal levels have contributed to our understanding of phylogenetic relationships within many groups. Genome rearrangements are an additional data type for defining relationships, with rearrangement synapomorphies identified across multiple orders and at many different taxonomic levels. Hymenoptera and Psocodea have greatly elevated rates of rearrangement offering both opportunities and pitfalls for identifying rearrangement synapomorphies in each group. Finally, insects are model systems for studying aberrant mt genomes, including truncated tRNAs and multichromosomal genomes. Greater integration of nuclear and mt genomic studies is necessary to further our understanding of insect genomic evolution.

Phylogenetic resolution power of microsatellites and various single-nucleotide polymorphism types assessed in 10 divergent chicken populations

There has been some debate over the question of which types of DNA variation are most appropriate to accurately reconstruct evolutionary events. We compared the capacity of microsatellites (STRs) and various types of single-nucleotide polymorphism (SNP) loci in the chicken genome. The SNP types differ in their location: in exons, introns and promoters. Genetic distances between all possible pairs of 10 populations were calculated for each marker type. STR loci, which are much more polymorphic than are SNPs, are considered to have occurred at recent time compared with old evolutionary events of SNPs. Using structure software, STR loci assigned individuals to their population much more correctly than did any other marker types, whereas SNPs at promoter regions gave the poorest ascription. Furthermore, 29 STR markers were even better than all 152 SNPs together. Ancient evolutionary events that produced genetic differences between the most distant populations such as Red Jungle Fowl and domestic chicken were detected better by exons and introns than by STR loci and promoters. The significant interactions found between marker types and populations suggest that marker types had different phylogenetic histories, possibly related to a different timescale.

miRNAs: small genes with big potential in metazoan phylogenetics

microRNAs (miRNAs) are a key component of gene regulatory networks and have been implicated in the regulation of virtually every biological process found in multicellular eukaryotes. What makes them interesting from a phylogenetic perspective is the high conservation of primary sequence between taxa, their accrual in metazoan genomes through evolutionary time, and the rarity of secondary loss in most metazoan taxa. Despite these properties, the use of miRNAs as phylogenetic markers has not yet been discussed within a clear conceptual framework. Here we highlight five properties of miRNAs that underlie their utility in phylogenetics: 1) The processes of miRNA biogenesis enable the identification of novel miRNAs without prior knowledge of sequence; 2) The continuous addition of miRNA families to metazoan genomes through evolutionary time; 3) The low level of secondary gene loss in most metazoan taxa; 4) The low substitution rate in the mature miRNA sequence; and 5) The small probability of convergent evolution of two miRNAs. Phylogenetic analyses using both Bayesian and parsimony methods on a eumetazoan miRNA data set highlight the potential of miRNAs to become an invaluable new tool, especially when used as an additional line of evidence, to resolve previously intractable nodes within the tree of life.

Probing the boundaries of orthology: the unanticipated rapid evolution of Drosophila centrosomin

The rapid evolution of essential developmental genes and their protein products is both intriguing and problematic. The rapid evolution of gene products with simple protein folds and a lack of well-characterized functional domains typically result in a low discovery rate of orthologous genes. Additionally, in the absence of orthologs it is difficult to study the processes and mechanisms underlying rapid evolution. In this study, we have investigated the rapid evolution of centrosomin (cnn), an essential gene encoding centrosomal protein isoforms required during syncytial development in Drosophila melanogaster. Until recently the rapid divergence of cnn made identification of orthologs difficult and questionable because Cnn violates many of the assumptions underlying models for protein evolution. To overcome these limitations, we have identified a group of insect orthologs and present conserved features likely to be required for the functions attributed to cnn in D. melanogaster. We also show that the rapid divergence of Cnn isoforms is apparently due to frequent coding sequence indels and an accelerated rate of intronic additions and eliminations. These changes appear to be buffered by multi-exon and multi-reading frame maximum potential ORFs, simple protein folds, and the splicing machinery. These buffering features also occur in other genes in Drosophila and may help prevent potentially deleterious mutations due to indels in genes with large coding exons and exon-dense regions separated by small introns. This work promises to be useful for future investigations of cnn and potentially other rapidly evolving genes and proteins.

Characterization of bud emergence 46 (BEM46) protein: sequence, structural, phylogenetic and subcellular localization analyses

The bud emergence 46 (BEM46) protein from Neurospora crassa belongs to the α/β-hydrolase superfamily. Recently, we have reported that the BEM46 protein is localized in the perinuclear ER and also forms spots close by the plasma membrane. The protein appears to be required for cell type-specific polarity formation in N. crassa. Furthermore, initial studies suggested that the BEM46 amino acid sequence is conserved in eukaryotes and is considered to be one of the widespread conserved "known unknown" eukaryotic genes. This warrants for a comprehensive phylogenetic analysis of this superfamily to unravel origin and molecular evolution of these genes in different eukaryotes. Herein, we observe that all eukaryotes have at least a single copy of a bem46 ortholog. Upon scanning of these proteins in various genomes, we find that there are expansions leading into several paralogs in vertebrates. Usingcomparative genomic analyses, we identified insertion/deletions (indels) in the conserved domain of BEM46 protein, which allow to differentiate fungal classes such as ascomycetes from basidiomycetes. We also find that exonic indels are able to differentiate BEM46 homologs of different eukaryotic lineage. Furthermore, we unravel that BEM46 protein from N. crassa possess a novel endoplasmic-retention signal (PEKK) using GFP-fusion tagging experiments. We propose that three residues namely a serine 188S, a histidine 292H and an aspartic acid 262D are most critical residues, forming a catalytic triad in BEM46 protein from N. crassa. We carried out a comprehensive study on bem46 genes from a molecular evolution perspective with combination of functional analyses. The evolutionary history of BEM46 proteins is characterized by exonic indels in lineage specific manner.

A phylogenomic and molecular signature based approach for characterization of the phylum Spirochaetes and its major clades: proposal for a taxonomic revision of the phylum

The Spirochaetes species cause many important diseases including syphilis and Lyme disease. Except for their containing a distinctive endoflagella, no other molecular or biochemical characteristics are presently known that are specific for either all Spirochaetes or its different families. We report detailed comparative and phylogenomic analyses of protein sequences from Spirochaetes genomes to understand their evolutionary relationships and to identify molecular signatures for this group. These studies have identified 38 conserved signature indels (CSIs) that are specific for either all members of the phylum Spirochaetes or its different main clades. Of these CSIs, a 3 aa insert in the FlgC protein is uniquely shared by all sequenced Spirochaetes providing a molecular marker for this phylum. Seven, six, and five CSIs in different proteins are specific for members of the families Spirochaetaceae, Brachyspiraceae, and Leptospiraceae, respectively. Of the 19 other identified CSIs, 3 are uniquely shared by members of the genera Sphaerochaeta, Spirochaeta, and Treponema, whereas 16 others are specific for the genus Borrelia. A monophyletic grouping of the genera Sphaerochaeta, Spirochaeta, and Treponema distinct from the genus Borrelia is also strongly supported by phylogenetic trees based upon concatenated sequences of 22 conserved proteins. The molecular markers described here provide novel and more definitive means for identification and demarcation of different main groups of Spirochaetes. To accommodate the extensive genetic diversity of the Spirochaetes as revealed by different CSIs and phylogenetic analyses, it is proposed that the four families of this phylum should be elevated to the order level taxonomic ranks (viz. Spirochaetales, Brevinematales ord. nov., Brachyspiriales ord. nov., and Leptospiriales ord. nov.). It is further proposed that the genera Borrelia and Cristispira be transferred to a new family Borreliaceae fam. nov. within the order Spirochaetales.

Evolution of protein indels in plants, animals and fungi

BACKGROUND

Insertions/deletions (indels) in protein sequences are useful as drug targets, protein structure predictors, species diagnostics and evolutionary markers. However there is limited understanding of indel evolutionary patterns. We sought to characterize indel patterns focusing first on the major groups of multicellular eukaryotes.

RESULTS

Comparisons of complete proteomes from a taxonically broad set of primarily Metazoa, Fungi and Viridiplantae yielded 299 substantial (>250aa) universal, single-copy (in-paralog only) proteins, from which 901 simple (present/absent) and 3,806 complex (multistate) indels were extracted. Simple indels are mostly small (1-7aa) with a most frequent size class of 1aa. However, even these simple looking indels show a surprisingly high level of hidden homoplasy (multiple independent origins). Among the apparently homoplasy-free simple indels, we identify 69 potential clade-defining indels (CDIs) that may warrant closer examination. CDIs show a very uneven taxonomic distribution among Viridiplante (13 CDIs), Fungi (40 CDIs), and Metazoa (0 CDIs). An examination of singleton indels shows an excess of insertions over deletions in nearly all examined taxa. This excess averages 2.31 overall, with a maximum observed value of 7.5 fold.

CONCLUSIONS

We find considerable potential for identifying taxon-marker indels using an automated pipeline. However, it appears that simple indels in universal proteins are too rare and homoplasy-rich to be used for pure indel-based phylogeny. The excess of insertions over deletions seen in nearly every genome and major group examined maybe useful in defining more realistic gap penalties for sequence alignment. This bias also suggests that insertions in highly conserved proteins experience less purifying selection than do deletions.

Molecular signatures for the phylum Aquificae and its different clades: proposal for division of the phylum Aquificae into the emended order Aquificales, containing the families Aquificaceae and Hydrogenothermaceae, and a new order Desulfurobacteriales ord. nov., containing the family Desulfurobacteriaceae

We report here detailed phylogenetic and comparative analyses on 11 sequenced genomes from the phylum Aquificae to identify molecular markers that are specific for the species from this phylum or its different families (viz. Aquificaceae, Hydrogenothermaceae and Desulfurobacteriaceae). In phylogenetic trees based on 16S rRNA gene or concatenated sequences for 32 conserved proteins, species from the three Aquificae families formed distinct clades. These trees also supported a strong relationship between the Aquificaceae and Hydrogenothermaceae families. In parallel, comparative analyses on protein sequences from Aquificae genomes have identified 46 conserved signature indels (CSIs) in broadly distributed proteins that are either exclusively or mainly found in members of the phylum Aquificae or its different families and subclades. Four of these CSIs, which are found in all sequenced Aquificae species, provide potential molecular markers for this phylum. Twelve, six and thirteen other CSIs that respectively are specific for the sequenced Aquificaceae, Hydrogenothermaceae and Desulfurobacteriaceae species provide molecular markers and novel tools for the identification of members of these families and for genetic and biochemical studies on them. Lastly, these studies have identified 11 CSIs in divergent proteins that are uniquely shared by members of the Aquificaceae and Hydrogenothermaceae families providing strong evidence that these two groups of bacteria shared a common ancestor exclusive of all other Aquificae (bacteria). The species from these two families are also very similar in their metabolic and physiological properties and they consist of aerobic or microaerophilic bacteria, which generally obtain energy by oxidation of hydrogen or reduced sulfur compounds by molecular oxygen. Based upon their strong association in phylogenetic trees, unique shared presence of large numbers of CSIs in different proteins, and similarities in their metabolic and physiological properties, it is proposed that the order Aquificales should be emended to include only the members of the families Aquificaceae and Hydrogenothermaceae. The members of the family Desulfurobacteriaceae, which are obligate anaerobes that strictly use hydrogen as electron donor, are now transferred to a new order Desulfurobacteriales ord. nov. The emended descriptions of the phylum Aquificae and its three families incorporating information for different molecular signatures are also provided.

The complete mitochondrial genomes of deep-sea squid (Bathyteuthis abyssicola), bob-tail squid (Semirossia patagonica) and four giant cuttlefish (Sepia apama, S. latimanus, S. lycidas and S. pharaonis), and their application to the phylogenetic analysis of Decapodiformes

We determined the complete mitochondrial (mt) genomes of the deep-sea squid (Bathyteuthis abyssicola; supperfamily Bathyteuthoidea), the bob-tail squid (Semirossia patagonica; order Sepiolida) and four giant cuttlefish (Sepia apama, S. latimanus, S. lycidas and S. pharaonis; order Sepiida). The unique structures of the mt genomes of Bathyteuthis and Semirossia provide new information about the evolution of decapodiform mt genomes. We show that the mt genome of B. abyssicola, like those of other oegopsids studied so far, has two long duplicated regions that include seven genes (COX1-3, ATP6 and ATP8, tRNA(Asn), and either ND2 or ND3) and that one of the duplicated COX3 genes has lost its function. The mt genome of S. patagonica is unlike any other decapodiforms and, like Nautilus, its ATP6 and ATP8 genes are not adjacent to each other. The four giant cuttlefish have identical mt gene order to other cuttlefish determined to date. Molecular phylogenetic analyses using maximum likelihood and Bayesian methods suggest that traditional order Sepioidea (Sepiolida+Sepiida) is paraphyletic and Sepia (cuttlefish) has the sister-relationship with all other decapodiforms. Taking both the phylogenetic analyses and the mt gene order analyses into account, it is likely that the octopus-type mt genome is an ancestral state and that it had maintained from at least the Cephalopoda ancestor to the common ancestor of Oegopsida, Myopsida and Sepiolida.

The complete mitochondrial genome of Solemya velum (Mollusca: Bivalvia) and its relationships with conchifera

Background

Bivalve mitochondrial genomes exhibit a wide array of uncommon features, like extensive gene rearrangements, large sizes, and unusual ways of inheritance. Species pertaining to the order Solemyida (subclass Opponobranchia) show many peculiar evolutionary adaptations, f.i. extensive symbiosis with chemoautotrophic bacteria. Despite Opponobranchia are central in bivalve phylogeny, being considered the sister group of all Autobranchia, a complete mitochondrial genome has not been sequenced yet.

Results

In this paper, we characterized the complete mitochondrial genome of the Atlantic awning clam Solemya velum: A-T content, gene arrangement and other features are more similar to putative ancestral mollusks than to other bivalves. Two supranumerary open reading frames are present in a large, otherwise unassigned, region, while the origin of replication could be located in a region upstream to the cox3 gene.

Conclusions

We show that S. velum mitogenome retains most of the ancestral conchiferan features, which is unusual among bivalve mollusks, and we discuss main peculiarities of this first example of an organellar genome coming from the subclass Opponobranchia. Mitochondrial genomes of Solemya (for bivalves) and Haliotis (for gastropods) seem to retain the original condition of mollusks, as most probably exemplified by Katharina.