Assessment of methods for amino acid matrix selection and their use on empirical data shows that ad hoc assumptions for choice of matrix are not justified

Comparative chloroplast genomics and phylogenetics of Fagopyrum esculentum ssp. ancestrale -a wild ancestor of cultivated buckwheat

BACKGROUND

Chloroplast genome sequences are extremely informative about species-interrelationships owing to its non-meiotic and often uniparental inheritance over generations. The subject of our study, Fagopyrum esculentum, is a member of the family Polygonaceae belonging to the order Caryophyllales. An uncertainty remains regarding the affinity of Caryophyllales and the asterids that could be due to undersampling of the taxa. With that background, having access to the complete chloroplast genome sequence for Fagopyrum becomes quite pertinent.

RESULTS

We report the complete chloroplast genome sequence of a wild ancestor of cultivated buckwheat, Fagopyrum esculentum ssp. ancestrale. The sequence was rapidly determined using a previously described approach that utilized a PCR-based method and employed universal primers, designed on the scaffold of multiple sequence alignment of chloroplast genomes. The gene content and order in buckwheat chloroplast genome is similar to Spinacia oleracea. However, some unique structural differences exist: the presence of an intron in the rpl2 gene, a frameshift mutation in the rpl23 gene and extension of the inverted repeat region to include the ycf1 gene. Phylogenetic analysis of 61 protein-coding gene sequences from 44 complete plastid genomes provided strong support for the sister relationships of Caryophyllales (including Polygonaceae) to asterids. Further, our analysis also provided support for Amborella as sister to all other angiosperms, but interestingly, in the bayesian phylogeny inference based on first two codon positions Amborella united with Nymphaeales.

CONCLUSION

Comparative genomics analyses revealed that the Fagopyrum chloroplast genome harbors the characteristic gene content and organization as has been described for several other chloroplast genomes. However, it has some unique structural features distinct from previously reported complete chloroplast genome sequences. Phylogenetic analysis of the dataset, including this new sequence from non-core Caryophyllales supports the sister relationship between Caryophyllales and asterids.

Unraveling the evolutionary history of the phosphoryl-transfer chain of the phosphoenolpyruvate:phosphotransferase system through phylogenetic analyses and genome context

Background

The phosphoenolpyruvate phosphotransferase system (PTS) plays a major role in sugar transport and in the regulation of essential physiological processes in many bacteria. The PTS couples solute transport to its phosphorylation at the expense of phosphoenolpyruvate (PEP) and it consists of general cytoplasmic phosphoryl transfer proteins and specific enzyme II complexes which catalyze the uptake and phosphorylation of solutes. Previous studies have suggested that the evolution of the constituents of the enzyme II complexes has been driven largely by horizontal gene transfer whereas vertical inheritance has been prevalent in the general phosphoryl transfer proteins in some bacterial groups. The aim of this work is to test this hypothesis by studying the evolution of the phosphoryl transfer proteins of the PTS.

Results

We have analyzed the evolutionary history of the PTS phosphoryl transfer chain (PTS-ptc) components in 222 complete genomes by combining phylogenetic methods and analysis of genomic context. Phylogenetic analyses alone were not conclusive for the deepest nodes but when complemented with analyses of genomic context and functional information, the main evolutionary trends of this system could be depicted.

Conclusion

The PTS-ptc evolved in bacteria after the divergence of early lineages such as Aquificales, Thermotogales and Thermus/Deinococcus. The subsequent evolutionary history of the PTS-ptc varied in different bacterial lineages: vertical inheritance and lineage-specific gene losses mainly explain the current situation in Actinobacteria and Firmicutes whereas horizontal gene transfer (HGT) also played a major role in Proteobacteria. Most remarkably, we have identified a HGT event from Firmicutes or Fusobacteria to the last common ancestor of the Enterobacteriaceae, Pasteurellaceae, Shewanellaceae and Vibrionaceae. This transfer led to extensive changes in the metabolic and regulatory networks of these bacteria including the development of a novel carbon catabolite repression system. Hence, this example illustrates that HGT can drive major physiological modifications in bacteria.

Spectral diversity of fluorescent proteins from the anthozoan Corynactis californica

Color morphs of the temperate, nonsymbiotic corallimorpharian Corynactis californica show variation in pigment pattern and coloring. We collected seven distinct color morphs of C. californica from subtidal locations in Monterey Bay, California, and found that tissue- and color-morph-specific expression of at least six different genes is responsible for this variation. Each morph contains at least three to four distinct genetic loci that code for these colors, and one morph contains at least five loci. These genes encode a subfamily of new GFP-like proteins, which fluoresce across the visible spectrum from green to red, while sharing between 75% to 89% pairwise amino-acid identity. Biophysical characterization reveals interesting spectral properties, including a bright yellow protein, an orange protein, and a red protein exhibiting a "fluorescent timer" phenotype. Phylogenetic analysis indicates that the FP genes from this species evolved together but that diversification of anthozoan fluorescent proteins has taken place outside of phylogenetic constraints, especially within the Corallimorpharia. The discovery of more examples of fluorescent proteins in a non-bioluminescent, nonsymbiotic anthozoan highlights possibilities of adaptive ecological significance unrelated to light regulation for algal symbionts. The patterns and colors of fluorescent proteins in C. californica and similar species may hold meaning for organisms that possess the visual pigments to distinguish them.

Dominance of Giardia assemblage B in León, Nicaragua

Giardiasis is a major problem in León, Nicaragua, yet despite this no data are available regarding the prevalence of different Giardia genotypes in this area. To address this question, a molecular analysis of Giardia isolates from humans and dogs living in the same area in León, Nicaragua was performed. Giardia isolates from 119 Nicaraguan patients and 8 dogs were successfully genotyped using single and/or nested beta-giardin PCR with subsequent restriction length fragment polymorphism (RFLP) analysis. The analyses of human samples yielded 94 (79%) assemblage B isolates and 25 (21%) assemblage A isolates. Only the non-human-associated assemblages C and D were found in the dog samples. Sixteen isolates with assemblage A pattern, 26 isolates with assemblage B pattern and all dog isolates were further characterized by sequencing the nested beta-giardin PCR product and by molecular analyses of the glutamate dehydrogenase (gdh) gene. Within the study area the assemblage A isolates were highly genetically homogenous, showing only sub-genotypes A2 (n=3) or A3 (n=13) at the beta-giardin locus and AII only at the gdh locus while assemblage B showed a high genetic polymorphism at both loci. Seven different sub-genotypes were identified within 13 of the sequenced assemblage B beta-giardin isolates. The remaining 13 sequenced assemblage B-isolates appeared to contain several different variants of the beta-giardin gene since the chromatograms displayed one to seven double peaks. The gdh sequences showed an even higher polymorphism since only 2 of 26 assemblage B isolates were without double peaks. Two mixed infections between assemblage A and B were found when the gdh gene was analyzed. Polymorphisms were also observed in the dog-associated assemblages C and D, but to a lesser extent than in assemblage B.

An improved general amino acid replacement matrix

Amino acid replacement matrices are an essential basis of protein phylogenetics. They are used to compute substitution probabilities along phylogeny branches and thus the likelihood of the data. They are also essential in protein alignment. A number of replacement matrices and methods to estimate these matrices from protein alignments have been proposed since the seminal work of Dayhoff et al. (1972). An important advance was achieved by Whelan and Goldman (2001) and their WAG matrix, thanks to an efficient maximum likelihood estimation approach that accounts for the phylogenies of sequences within each training alignment. We further refine this method by incorporating the variability of evolutionary rates across sites in the matrix estimation and using a much larger and diverse database than BRKALN, which was used to estimate WAG. To estimate our new matrix (called LG after the authors), we use an adaptation of the XRATE software and 3,912 alignments from Pfam, comprising approximately 50,000 sequences and approximately 6.5 million residues overall. To evaluate the LG performance, we use an independent sample consisting of 59 alignments from TreeBase and randomly divide Pfam alignments into 3,412 training and 500 test alignments. The comparison with WAG and JTT shows a clear likelihood improvement. With TreeBase, we find that 1) the average Akaike information criterion gain per site is 0.25 and 0.42, when compared with WAG and JTT, respectively; 2) LG is significantly better than WAG for 38 alignments (among 59), and significantly worse with 2 alignments only; and 3) tree topologies inferred with LG, WAG, and JTT frequently differ, indicating that using LG impacts not only the likelihood value but also the output tree. Results with the test alignments from Pfam are analogous. LG and a PHYML implementation can be downloaded from http://atgc.lirmm.fr/LG.

Revealing pancrustacean relationships: phylogenetic analysis of ribosomal protein genes places Collembola (springtails) in a monophyletic Hexapoda and reinforces the discrepancy between mitochondrial and nuclear DNA markers

BACKGROUND

In recent years, several new hypotheses on phylogenetic relations among arthropods have been proposed on the basis of DNA sequences. One of the challenged hypotheses is the monophyly of hexapods. This discussion originated from analyses based on mitochondrial DNA datasets that, due to an unusual positioning of Collembola, suggested that the hexapod body plan evolved at least twice. Here, we re-evaluate the position of Collembola using ribosomal protein gene sequences.

RESULTS

In total 48 ribosomal proteins were obtained for the collembolan Folsomia candida. These 48 sequences were aligned with sequence data on 35 other ecdysozoans. Each ribosomal protein gene was available for 25% to 86% of the taxa. However, the total sequence information was unequally distributed over the taxa and ranged between 4% and 100%. A concatenated dataset was constructed (5034 inferred amino acids in length), of which ~66% of the positions were filled. Phylogenetic tree reconstructions, using Maximum Likelihood, Maximum Parsimony, and Bayesian methods, resulted in a topology that supports monophyly of Hexapoda.

CONCLUSION

Although ribosomal proteins in general may not evolve independently, they once more appear highly valuable for phylogenetic reconstruction. Our analyses clearly suggest that Hexapoda is monophyletic. This underpins the inconsistency between nuclear and mitochondrial datasets when analyzing pancrustacean relationships. Caution is needed when applying mitochondrial markers in deep phylogeny.

Analysis of Acropora muricata calmodulin (CaM) indicates that scleractinian corals possess the ancestral exon/intron organization of the eumetazoan CaM gene

Calmodulin (CaM), belonging to the tropinin C (TnC) superfamily, is one of the calcium-binding proteins that are highly conserved in their protein and gene structure. Based on the structure comparison among published vertebrate and invertebrate CaM, it is proposed that the ancestral form of eumetazoan CaM genes should have five exons and four introns (four-intron hypothesis). In this study, we determined the gene structure of CaM in the coral Acropora muricata, an anthozoan cnidarian representing the basal position in animal evolution. A CaM clone was isolated from a cDNA library constructed from the spawned eggs of A. muricata. This clone was composed of 908 nucleotides, including 162 base pairs (bp) of 5'-untranslated region (UTR), 296 bp of 3'-UTR, and an open reading frame 450 bp in length. The deduced amino acid indicated that the Acropora CaM protein is identical to that of the actiniarian, Metridinium senile, and has four putative calcium-binding domains highly similar to those of other vertebrate or invertebrate CaMs. Southern blot analysis revealed that Acropora CaM is a putative single-copy gene in the nuclear genome. Genomic sequencing showed that Acropora CaM was composed of five exons and four introns, with intron II not corresponding to any region in the actiniarian CaM gene, which possesses only four exons and three introns. Our results highlight that the coral CaM gene isolated from A. muricata has four introns at the predicted positions of the early metazoan CaM gene organization, providing the first evidence from the basal eumetazoan phylum to support the four-intron hypothesis.

The human phylome

The human phylome, which includes evolutionary relationships of all human proteins and their homologs among thirty-nine fully sequenced eukaryotes, is reconstructed.

Background:

Phylogenomics analyses serve to establish evolutionary relationships among organisms and their genes. A phylome, the complete collection of all gene phylogenies in a genome, constitutes a valuable source of information, but its use in large genomes still constitutes a technical challenge. The use of phylomes also requires the development of new methods that help us to interpret them.

Results:

We reconstruct here the human phylome, which includes the evolutionary relationships of all human proteins and their homologs among 39 fully sequenced eukaryotes. Phylogenetic techniques used include alignment trimming, branch length optimization, evolutionary model testing and maximum likelihood and Bayesian methods. Although differences with alternative topologies are minor, most of the trees support the Coelomata and Unikont hypotheses as well as the grouping of primates with laurasatheria to the exclusion of rodents. We assess the extent of gene duplication events and their relationship with the functional roles of the protein families involved. We find support for at least one, and probably two, rounds of whole genome duplications before vertebrate radiation. Using a novel algorithm that is independent from a species phylogeny, we derive orthology and paralogy relationships of human proteins among eukaryotic genomes.

Conclusion:

Topological variations among phylogenies for different genes are to be expected, highlighting the danger of gene-sampling effects in phylogenomic analyses. Several links can be established between the functions of gene families duplicated at certain phylogenetic splits and major evolutionary transitions in those lineages. The pipeline implemented here can be easily adapted for use in other organisms.

Climate-driven diversification in two widespread Galerida larks

BACKGROUND

The major impact of Plio-Pleistocene climatic oscillations on the current genetic structure of many species is widely recognised but their importance in driving speciation remains a matter of controversies. In addition, since most studies focused on Europe and North America, the influence of many other biogeographic barriers such as the Sahara remains poorly understood. In this paper, climate-driven diversification was investigated by using a comparative phylogeographic approach in combination with phenotypic data in two avian species groups distributed on both sides of the deserts belt of Africa and Asia. In particular, we tested whether: 1) vicariance diversification events are concomitant with past climatic events; and 2) current ecological factors (using climate and competition as proxies) contribute to phenotypic divergence between allopatric populations.

RESULTS

Mitochondrial and nuclear sequence data indicated that the crested and Thekla lark species groups diverged in the early Pliocene and that subsequent speciation events were congruent with major late Pliocene and Pleistocene climatic events. In particular, steep increase in aridity in Africa near 2.8 and 1.7 million years ago were coincident with two north-south vicariance speciation events mediated by the Sahara. Subsequent glacial cycles of the last million years seem to have shaped patterns of genetic variation within the two widespread species (G. cristata and G. theklae). The Sahara appears to have allowed dispersal from the tropical areas during climatic optima but to have isolated populations north and south of it during more arid phases. Phenotypic variation did not correlate with the history of populations, but was strongly influenced by current ecological conditions. In particular, our results suggested that (i) desert-adapted plumage evolved at least three times and (ii) variation in body size was mainly driven by interspecific competition, but the response to competition was stronger in more arid areas.

CONCLUSION

Climatic fluctuations of the Plio-Pleistocene strongly impacted diversification patterns in the Galerida larks. Firstly, we found that cladogenesis coincides with major climatic changes, and the Sahara appears to have played a key role in driving speciation events. Secondly, we found that morphology and plumage were strongly determined by ecological factors (interspecific competition, climate) following vicariance.

Bacterial communities of Bartonella-positive fleas: diversity and community assembly patterns

We investigated the bacterial communities of nine Bartonella-positive fleas (n = 6 Oropsylla hirsuta fleas and n = 3 Oropsylla montana fleas), using universal primers, clone libraries, and DNA sequencing. DNA sequences were used to classify bacteria detected in a phylogenetic context, to explore community assembly patterns within individual fleas, and to survey diversity patterns in dominant lineages.

Yeast forms dominate fungal diversity in the deep oceans

Fungi are the principal degraders of biomass in most terrestrial ecosystems. In contrast to surface environments, deep-sea environmental gene libraries have suggested that fungi are rare and non-diverse in high-pressure marine environments. Here, we report the diversity of fungi from 11 deep-sea samples from around the world representing depths from 1,500 to 4,000 m (146-388 atm) and two shallower water column samples (250 and 500m). We sequenced 239 clones from 10 fungal-specific 18S rRNA gene libraries constructed from these samples, from which we detected only 18 fungal 18S-types in deep-sea samples. Our phylogenetic analyses show that a total of only 32 fungal 18S-types have so far been recovered from deep-sea habitats, and our results suggest that fungi, in general, are relatively rare in the deep-sea habitats we sampled. The fungal diversity detected suggests that deep-sea environments host an evolutionarily diverse array of fungi dominated by groups of distantly related yeasts, although four putative filamentous fungal 18S-types were detected. The majority of our new sequences branch close to known fungi found in surface environments. This pattern contradicts the proposal that deep-sea and hydrothermal vent habitats represent ancient ecosystems, and demonstrates a history of frequent dispersal between terrestrial and deep-sea habitats.

Pathogen evolution and disease emergence in carnivores

Emerging infectious diseases constitute some of the most pressing problems for both human and domestic animal health, and biodiversity conservation. Currently it is not clear whether the removal of past constraints on geographical distribution and transmission possibilities for pathogens alone are sufficient to give rise to novel host-pathogen combinations, or whether pathogen evolution is also generally required for establishment in novel hosts. Canine distemper virus (CDV) is a morbillivirus that is prevalent in the world dog population and poses an important conservation threat to a diverse range of carnivores. We performed an extensive phylogenetic and molecular evolution analysis on complete sequences of all CDV genes to assess the role of selection and recombination in shaping viral genetic diversity and driving the emergence of CDV in non-dog hosts. We tested the specific hypothesis that molecular adaptation at known receptor-binding sites of the haemagglutinin gene is associated with independent instances of the spread of CDV to novel non-dog hosts in the wild. This hypothesis was upheld, providing compelling evidence that repeated evolution at known functional sites (in this case residues 530 and 549 of the haemagglutinin molecule) is associated with multiple independent occurrences of disease emergence in a range of novel host species.

Molecular phylogenetics of central Canadian Physidae (Pulmonata: Basommatophora)

The phylogenetic relationships of four nominal south-central Canadian freshwater physids ( Physa (sensu lato) skinneri Taylor, 1954, Physa integra (Haldeman, 1841), Physa gyrina (Say, 1821), and the endemic Physa winnipegensis Pip, 2004) were studied by analyses of combined partial sequences coding for mitochondrial 16S and cytochrome c oxidase (COI) genes. Maximum parsimony and neighbour joining analyses, as well as comparisons with published sequences, supported four major clades of physids, of which three are represented in central Canada. Physa skinneri and P. integra were placed within the Physa fontinalis (L., 1758) and Physa acuta Draparnaud, 1805 clades, respectively. Physa winnipegensis formed a distinct branch within the P. acuta group. An additional, previously unreported and unclassified, morph within the P. acuta group was identified from Lake Winnipeg.

Phylogenetic analysis of the NS5 gene of dengue viruses isolated in Ecuador

Dengue virus (DENV) is a member of the genus Flavivirus of the family Flaviviridae. DENV causes a wide range of diseases in humans, from the acute febrile illness dengue fever (DF) to life-threatening dengue hemorrhagic fever/dengue shock syndrome (DHF/DSS). There is not knowledge of the genetic relations among DENV circulating in Ecuador. Given the emerging behaviour of DENV, a single tube RT-PCR assay using a pair of consensus primers to target the NS5 coding region has been recently validated for rapid detection of flaviviruses. In order to gain insight into the degree of genetic variation of DENV strains isolated in Ecuador, DENV NS5 sequences from 23 patients were obtained by direct sequencing of PCR fragments using the mentioned one step RT-PCR assay. Phylogenetic analysis carried out using the 23 Ecuadorian DENV NS5 sequences, as well as 56 comparable sequences from DENV strains isolated elsewhere, revealed a close genetic relation among Ecuadorian strains and DENV isolates of Caribbean origin. The use of partial NS5 gene sequences may represent a useful alternative for a rapid phylogenetic analysis of DENV outbreaks.

Global isolation by distance despite strong regional phylogeography in a small metazoan

BACKGROUND

Small vagile eukaryotic organisms, which comprise a large proportion of the Earth's biodiversity, have traditionally been thought to lack the extent of population structuring and geographic speciation observed in larger taxa. Here we investigate the patterns of genetic diversity, amongst populations of the salt lake microscopic metazoan Brachionus plicatilis s. s. (sensu stricto) (Rotifera: Monogononta) on a global scale. We examine the phylogenetic relationships of geographic isolates from four continents using a 603 bp fragment of the mitochondrial COI gene to investigate patterns of phylogeographic subdivision in this species. In addition we investigate the relationship between genetic and geographic distances on a global scale to try and reconcile the paradox between the high vagility of this species and the previously reported patterns of restricted gene flow, even over local spatial scales.

RESULTS

Analysis of global sequence diversity of B. plicatilis s. s. reveals the presence of four allopatric genetic lineages: North American-Far East Asian, Western Mediterranean, Australian, and an Eastern Mediterranean lineage represented by a single isolate. Geographically orientated substructure is also apparent within the three best sampled lineages. Surprisingly, given this strong phylogeographic structure, B. plicatilis s. s. shows a significant correlation between geographic and genetic distance on a global scale ('isolation by distance' - IBD).

CONCLUSION

Despite its cosmopolitan distribution and potential for high gene flow, B. plicatilis s. s. is strongly structured at a global scale. IBD patterns have traditionally been interpreted to indicate migration-drift equilibrium, although in this system equilibrium conditions are incompatible with the observed genetic structure. Instead, we suggest the pattern may have arisen through persistent founder effects, acting in a similar fashion to geographic barriers for larger organisms. Our data indicates that geographic speciation, contrary to historical views, is likely to be very important in microorganisms. By presenting compelling evidence for geographic speciation in a small eukaryote we add to the growing body of evidence that is forcing us to rethink our views of global biodiversity.

ChsVb, a class VII chitin synthase involved in septation, is critical for pathogenicity in Fusarium oxysporum

A new myosin motor-like chitin synthase gene, chsVb, has been identified in the vascular wilt fungus Fusarium oxysporum f. sp. lycopersici. Phylogenetic analysis of the deduced amino acid sequence of the chsVb chitin synthase 2 domain (CS2) revealed that ChsVb belongs to class VII chitin synthases. The ChsVb myosin motor-like domain (MMD) is shorter than the MMD of class V chitin synthases and does not contain typical ATP-binding motifs. Targeted disrupted single (DeltachsVb) and double (DeltachsV DeltachsVb) mutants were unable to infect and colonize tomato plants or grow invasively on tomato fruit tissue. These strains were hypersensitive to compounds that interfere with fungal cell wall assembly, produced lemon-like shaped conidia, and showed swollen balloon-like structures in hyphal subapical regions, thickened walls, aberrant septa, and intrahyphal hyphae. Our results suggest that the chsVb gene is likely to function in polarized growth and confirm the critical importance of cell wall integrity in the complex infection process of this fungus.

Bayesian coalescent inference of hepatitis A virus populations: evolutionary rates and patterns

Hepatitis A virus (HAV) is a hepatotropic member of the family Picornaviridae. Previous studies suggested that HAV may evolve more slowly than other members of the family. To estimate HAV substitution rates precisely, we used a Bayesian Markov chain Monte Carlo (MCMC) approach on temporally sampled HAV VP1 full-length sequences from strains isolated in France. A mean rate of evolutionary change of 9.76×10−4 nucleotide substitution per site per year was found. The results also revealed that the synonymous rate found for HAV is lower than that of other members of the family. Bayesian skyline plots revealed a sharp decline in the effective number of infections in 1996, coinciding with the introduction of HAV vaccine.

Evolution of the cutinase gene family: evidence for lateral gene transfer of a candidate Phytophthora virulence factor

Lateral gene transfer (LGT) can facilitate the acquisition of new functions in recipient lineages, which may enable them to colonize new environments. Several recent publications have shown that gene transfer between prokaryotes and eukaryotes occurs with appreciable frequency. Here we present a study of interdomain gene transfer of cutinases -- well documented virulence factors in fungi -- between eukaryotic plant pathogens Phytophthora species and prokaryotic bacterial lineages. Two putative cutinase genes were cloned from Phytophthora brassicae and Northern blotting experiments showed that these genes are expressed early during the infection of the host Arabidopsis thaliana and induced during cyst germination of the pathogen. Analysis of the gene organisation of this gene family in Phytophthora ramorum and P. sojae showed three and ten copies in tight succession within a region of 5 and 25 kb, respectively, probably indicating a recent expansion in Phytophthora lineages by gene duplications. Bioinformatic analyses identified orthologues only in three genera of Actinobacteria, and in two distantly related eukaryotic groups: oomycetes and fungi. Together with phylogenetic analyses this limited distribution of the gene in the tree of life strongly support a scenario where cutinase genes originated after the origin of land plants in a microbial lineage living in proximity of plants and subsequently were transferred between distantly related plant-degrading microbes. More precisely, a cutinase gene was likely acquired by an ancestor of P. brassicae, P. sojae, P. infestans and P. ramorum, possibly from an actinobacterial source, suggesting that gene transfer might be an important mechanism in the evolution of their virulence. These findings could indeed provide an interesting model system to study acquisition of virulence factors in these important plant pathogens.

Evolution of receptors for growth hormone and somatolactin in fish and land vertebrates: lessons from the lungfish and sturgeon orthologues

Two cognate hormones, growth hormone (GH) and somatolactin (SL), control several important physiological processes in vertebrates. Knowledge about GH and its receptor (GHR) has accumulated over the last decades. However, much less is known about SL and its receptor (SLR). SL is found only in fish (including lungfish), suggesting that it was present in the common ancestor of vertebrates, but was lost secondarily in the lineage leading to land vertebrates after the lungfish branched off. SLR was suggested to be a duplicated copy of GHR acquired only in teleosts via the fish-specific genome duplication (FSGD). This scenario (i.e., the existence of SL but not SLR in the vertebrate ancestors) is intriguing but contested. In this study, we first evaluated the plausibility of this scenario through synteny analyses and found that the loci for GHR and SLR are located in syntenic genomic positions, whereas the loci for GH and SL are not. Next, we cloned GHRs of lungfish and sturgeon, which possess SL but did not undergo the FSGD (i.e., they should not possess SLR). Their phylogenetic positions in the GHR/SLR gene tree further support the fish-specific scenario for the GHR-SLR duplication. Interestingly, their sequences share greater similarity with teleost SLRs and reptilian/amphibian GHRs than with the GHRs of mammals, birds, and teleosts. On the basis of these results, we discuss the validity of the nomenclature of the teleost-specific copy of GHR as SLR and an ancestral receptor(s) for SL before the evolution of SLR during the FSGD.

Avian UCP: the killjoy in the evolution of the mitochondrial uncoupling proteins

The understanding of mitochondrial functioning is of prime importance since it combines the production of energy as adenosine triphosphate (ATP) with an efficient chain of redox reactions, but also with the unavoidable production of reactive oxygen species (ROS) involved in aging. Mitochondrial respiration may be uncoupled from ATP synthesis by a proton leak induced by the thermogenic uncoupling protein 1 (UCP1). Mild uncoupling activity, as proposed for UCP2, UCP3, and avian UCP could theoretically control ROS production, but the nature of their transport activities is far from being definitively understood. The recent discovery of a UCP1 gene in fish has balanced the evolutionary view of uncoupling protein history. The thermogenic proton transport of mammalian UCP1 seems now to be a late evolutionary characteristic and the hypothesis that ancestral UCPs may carry other substrates is tempting. Using in silico genome analyses among taxa and a biochemical approach, we present a detailed phylogenetic analysis of UCPs and investigate whether avian UCP is a good candidate for pleiotropic mitochondrial activities, knowing that only one UCP has been characterized in the avian genome, unlike all other vertebrates. We show, here, that the avian class seems to be the only vertebrate lineage lacking two of the UCP1/2/3 homologues present in fish and mammals. We suggest, based on phylogenetic evidence and synteny of the UCP genes, that birds have lost UCP1 and UCP2. The phylogeny also supports the history of two rounds of duplication during vertebrate evolution. The avian uncoupling protein then represents a unique opportunity to explore how UCPs' activities are controlled, but also to understand why birds exhibit such a particular relationship between high metabolism and slow rate of aging.

Molecular phylogeny of Malagasy reed frogs, Heterixalus, and the relative performance of bioacoustics and color-patterns for resolving their systematics

The members of the genus Heterixalus constitute one of the endemic frog radiations in Madagascar. Here we present a complete species-level phylogeny based on DNA sequences (4876 base pairs) of three nuclear and four mitochondrial markers to clarify the phylogenetic relationships among and within all known species of this genus, as well as the phylogenetic position of the monospecific Seychellean Tachycnemis seychellensis. Although the performance to resolve supported clades of Heterixalus species differed among the investigated gene fragments when analyzed separately, we could identify five well-supported species groups within Heterixalus in the combined analysis of all gene fragments. Our data strongly support a Heterixalus-Tachycnemis clade, and indicate the probable monophyly of Heterixalus placed sister to Tachycnemis. However, the diversification of these lineages may have happened in a short interval of time, leading to an unstable placement of Tachycnemis in the single-gene fragment phylogenies. Referring to the hitherto existing classification of Heterixalus, which is predominantly based on chromatic and bioacoustic characters, we examined the relative performance of these data sets relative to our molecular phylogeny. A Bayesian tree reconstructed with a bioacoustic data set yielded a higher resemblance to the molecular phylogeny than a tree constructed using a chromatic data set, which supports the importance of bioacoustic characters for systematic analyses of these anurans.

Comparative genomics of ParaHox clusters of teleost fishes: gene cluster breakup and the retention of gene sets following whole genome duplications

Background

The evolutionary lineage leading to the teleost fish underwent a whole genome duplication termed FSGD or 3R in addition to two prior genome duplications that took place earlier during vertebrate evolution (termed 1R and 2R). Resulting from the FSGD, additional copies of genes are present in fish, compared to tetrapods whose lineage did not experience the 3R genome duplication. Interestingly, we find that ParaHox genes do not differ in number in extant teleost fishes despite their additional genome duplication from the genomic situation in mammals, but they are distributed over twice as many paralogous regions in fish genomes.

Results

We determined the DNA sequence of the entire ParaHox C1 paralogon in the East African cichlid fish Astatotilapia burtoni, and compared it to orthologous regions in other vertebrate genomes as well as to the paralogous vertebrate ParaHox D paralogons. Evolutionary relationships among genes from these four chromosomal regions were studied with several phylogenetic algorithms. We provide evidence that the genes of the ParaHox C paralogous cluster are duplicated in teleosts, just as it had been shown previously for the D paralogon genes. Overall, however, synteny and cluster integrity seems to be less conserved in ParaHox gene clusters than in Hox gene clusters. Comparative analyses of non-coding sequences uncovered conserved, possibly co-regulatory elements, which are likely to contain promoter motives of the genes belonging to the ParaHox paralogons.

Conclusion

There seems to be strong stabilizing selection for gene order as well as gene orientation in the ParaHox C paralogon, since with a few exceptions, only the lengths of the introns and intergenic regions differ between the distantly related species examined. The high degree of evolutionary conservation of this gene cluster's architecture in particular – but possibly clusters of genes more generally – might be linked to the presence of promoter, enhancer or inhibitor motifs that serve to regulate more than just one gene. Therefore, deletions, inversions or relocations of individual genes could destroy the regulation of the clustered genes in this region. The existence of such a regulation network might explain the evolutionary conservation of gene order and orientation over the course of hundreds of millions of years of vertebrate evolution. Another possible explanation for the highly conserved gene order might be the existence of a regulator not located immediately next to its corresponding gene but further away since a relocation or inversion would possibly interrupt this interaction. Different ParaHox clusters were found to have experienced differential gene loss in teleosts. Yet the complete set of these homeobox genes was maintained, albeit distributed over almost twice the number of chromosomes. Selection due to dosage effects and/or stoichiometric disturbance might act more strongly to maintain a modal number of homeobox genes (and possibly transcription factors more generally) per genome, yet permit the accumulation of other (non regulatory) genes associated with these homeobox gene clusters.

Phylogenetic analysis of mitochondrial protein coding genes confirms the reciprocal paraphyly of Hexapoda and Crustacea

BACKGROUND

The phylogeny of Arthropoda is still a matter of harsh debate among systematists, and significant disagreement exists between morphological and molecular studies. In particular, while the taxon joining hexapods and crustaceans (the Pancrustacea) is now widely accepted among zoologists, the relationships among its basal lineages, and particularly the supposed reciprocal paraphyly of Crustacea and Hexapoda, continues to represent a challenge. Several genes, as well as different molecular markers, have been used to tackle this problem in molecular phylogenetic studies, with the mitochondrial DNA being one of the molecules of choice. In this study, we have assembled the largest data set available so far for Pancrustacea, consisting of 100 complete (or almost complete) sequences of mitochondrial genomes. After removal of unalignable sequence regions and highly rearranged genomes, we used nucleotide and inferred amino acid sequences of the 13 protein coding genes to reconstruct the phylogenetic relationships among major lineages of Pancrustacea. The analysis was performed with Bayesian inference, and for the amino acid sequences a new, Pancrustacea-specific, matrix of amino acid replacement was developed and used in this study.

RESULTS

Two largely congruent trees were obtained from the analysis of nucleotide and amino acid datasets. In particular, the best tree obtained based on the new matrix of amino acid replacement (MtPan) was preferred over those obtained using previously available matrices (MtArt and MtRev) because of its higher likelihood score. The most remarkable result is the reciprocal paraphyly of Hexapoda and Crustacea, with some lineages of crustaceans (namely the Malacostraca, Cephalocarida and, possibly, the Branchiopoda) being more closely related to the Insecta s.s. (Ectognatha) than two orders of basal hexapods, Collembola and Diplura. Our results confirm that the mitochondrial genome, unlike analyses based on morphological data or nuclear genes, consistently supports the non monophyly of Hexapoda.

CONCLUSION

The finding of the reciprocal paraphyly of Hexapoda and Crustacea suggests an evolutionary scenario in which the acquisition of the hexapod condition may have occurred several times independently in lineages descending from different crustacean-like ancestors, possibly as a consequence of the process of terrestrialization. If this hypothesis was confirmed, we should therefore re-think our interpretation of the evolution of the Arthropoda, where terrestrialization may have led to the acquisition of similar anatomical features by convergence. At the same time, the disagreement between reconstructions based on morphological, nuclear and mitochondrial data sets seems to remain, despite the use of larger data sets and more powerful analytical methods.

"Candidatus Bacilloplasma," a novel lineage of Mollicutes associated with the hindgut wall of the terrestrial isopod Porcellio scaber (Crustacea: Isopoda)

Pointed, rod-shaped bacteria colonizing the cuticular surface of the hindgut of the terrestrial isopod crustacean Porcellio scaber (Crustacea: Isopoda) were investigated by comparative 16S rRNA gene sequence analysis and electron microscopy. The results of phylogenetic analysis, and the absence of a cell wall, affiliated these bacteria with the class Mollicutes, within which they represent a novel and deeply branched lineage, sharing less than 82.6% sequence similarity to known Mollicutes. The lineage has been positioned as a sister group to the clade comprising the Spiroplasma group, the Mycoplasma pneumoniae group, and the Mycoplasma hominis group. The specific signature sequence was identified and used as a probe in in situ hybridization, which confirmed that the retrieved sequences originate from the attached rod-shaped bacteria from the hindgut of P. scaber and made it possible to detect these bacteria in their natural environment. Scanning and transmission electron microscopy revealed a spherically shaped structure at the tapered end of the rod-shaped bacteria, enabling their specific and exclusive attachment to the tip of the cuticular spines on the inner surface of the gut. Specific adaptation to the gut environment, as well as phylogenetic positioning, indicate the long-term association and probable coevolution of the bacteria and the host. Taking into account their pointed, rod-shaped morphology and their phylogenetic position, the name "Candidatus Bacilloplasma" has been proposed for this new lineage of bacteria specifically associated with the gut surface of P. scaber.

Persistent genetic signatures of colonization inBrachionus manjavacasrotifers in the Iberian Peninsula

Recent phylogeographical assessments have consistently shown that continental zooplankton display high levels of population subdivision, despite the high dispersal capacity of their diapausing propagules. As such, there is an apparent paradox between observed cosmopolitanism in the zooplankton that is associated with long-distance dispersal, and strong phylogeographical structures at a regional scale. Such population dynamics, far from migration-drift equilibrium, have been shown in the rotifer species complex Brachionus plicatilis, a group of over a dozen species inhabiting salt lakes and coastal lagoons worldwide. Here we present the mitochondrial DNA phylogeography of one of these species, Brachionus manjavacas, in the Iberian Peninsula, where it often co-occurs with the morphologically similar species B. plicatilis sensu stricto. We obtained sequences from 233 individuals from diapausing eggs and clonal cultures from 16 lakes in the Iberian Peninsula, and a Tunisian lake. Two strongly supported deep mitochondrial DNA clades were found (A and B). Phylogenetic and nested clade analysis showed that clade A has a strong phylogeographical structure, with a strong similarity of phylogeographical patterns between B. manjavacas clade A and B. plicatilis s.s. These include (i) signatures of allopatric fragmentation between central and southern populations, and (ii) range expansions in the Iberian Peninsula, both likely to have occurred during the Pleistocene. We find evidence for a glacial refugium in the Guadiana basin. Clades A and B co-occurred in several of these lakes because of range expansion and secondary contact between both clades. The co-occurrence between B. plicatilis s.s. and B. manjavacas is not recent, and both species might have experienced similar environmental challenges during the Pleistocene. The strong correlation of genetic and geographical distance found suggests that historical events can lead to such correlation, mirroring the effects of 'isolation by distance' in equilibrium populations.

Dynamics of the evolution of Batesian mimicry: molecular phylogenetic analysis of ant-mimicking Myrmarachne (Araneae: Salticidae) species and their ant models

Batesian mimicry is seen as an example of evolution by natural selection, with predation as the main driving force. The mimic is under selective pressure to resemble its model, whereas it is disadvantageous for the model to be associated with the palatable mimic. In consequence one might expect there to be an evolutionary arms race, similar to the one involving host-parasite coevolution. In this study, the evolutionary dynamics of a Batesian mimicry system of model ants and ant-mimicking salticids is investigated by comparing the phylogenies of the two groups. Although Batesian mimics are expected to coevolve with their models, we found the phylogenetic patterns of the models and the mimics to be indicative of adaptive radiation by the mimic rather than co-speciation between the mimic and the model. This shows that there is strong selection pressure on Myrmarachne, leading to a high degree of polymorphism. There is also evidence of sympatric speciation in Myrmarachne, the reproductive isolation possibly driven by female mate choice in polymorphic species.

Reconstructing the evolution of the mitochondrial ribosomal proteome

For production of proteins that are encoded by the mitochondrial genome, mitochondria rely on their own mitochondrial translation system, with the mitoribosome as its central component. Using extensive homology searches, we have reconstructed the evolutionary history of the mitoribosomal proteome that is encoded by a diverse subset of eukaryotic genomes, revealing an ancestral ribosome of alpha-proteobacterial descent that more than doubled its protein content in most eukaryotic lineages. We observe large variations in the protein content of mitoribosomes between different eukaryotes, with mammalian mitoribosomes sharing only 74 and 43% of its proteins with yeast and Leishmania mitoribosomes, respectively. We detected many previously unidentified mitochondrial ribosomal proteins (MRPs) and found that several have increased in size compared to their bacterial ancestral counterparts by addition of functional domains. Several new MRPs have originated via duplication of existing MRPs as well as by recruitment from outside of the mitoribosomal proteome. Using sensitive profile-profile homology searches, we found hitherto undetected homology between bacterial and eukaryotic ribosomal proteins, as well as between fungal and mammalian ribosomal proteins, detecting two novel human MRPs. These newly detected MRPs constitute, along with evolutionary conserved MRPs, excellent new screening targets for human patients with unresolved mitochondrial oxidative phosphorylation disorders.

Horizontal transfer of a eukaryotic plastid-targeted protein gene to cyanobacteria

Background

Horizontal or lateral transfer of genetic material between distantly related prokaryotes has been shown to play a major role in the evolution of bacterial and archaeal genomes, but exchange of genes between prokaryotes and eukaryotes is not as well understood. In particular, gene flow from eukaryotes to prokaryotes is rarely documented with strong support, which is unusual since prokaryotic genomes appear to readily accept foreign genes.

Results

Here, we show that abundant marine cyanobacteria in the related genera Synechococcus and Prochlorococcus acquired a key Calvin cycle/glycolytic enzyme from a eukaryote. Two non-homologous forms of fructose bisphosphate aldolase (FBA) are characteristic of eukaryotes and prokaryotes respectively. However, a eukaryotic gene has been inserted immediately upstream of the ancestral prokaryotic gene in several strains (ecotypes) of Synechococcus and Prochlorococcus. In one lineage this new gene has replaced the ancestral gene altogether. The eukaryotic gene is most closely related to the plastid-targeted FBA from red algae. This eukaryotic-type FBA once replaced the plastid/cyanobacterial type in photosynthetic eukaryotes, hinting at a possible functional advantage in Calvin cycle reactions. The strains that now possess this eukaryotic FBA are scattered across the tree of Synechococcus and Prochlorococcus, perhaps because the gene has been transferred multiple times among cyanobacteria, or more likely because it has been selectively retained only in certain lineages.

Conclusion

A gene for plastid-targeted FBA has been transferred from red algae to cyanobacteria, where it has inserted itself beside its non-homologous, functional analogue. Its current distribution in Prochlorococcus and Synechococcus is punctate, suggesting a complex history since its introduction to this group.

MultiPhyl: a high-throughput phylogenomics webserver using distributed computing

With the number of fully sequenced genomes increasing steadily, there is greater interest in performing large-scale phylogenomic analyses from large numbers of individual gene families. Maximum likelihood (ML) has been shown repeatedly to be one of the most accurate methods for phylogenetic construction. Recently, there have been a number of algorithmic improvements in maximum-likelihood-based tree search methods. However, it can still take a long time to analyse the evolutionary history of many gene families using a single computer. Distributed computing refers to a method of combining the computing power of multiple computers in order to perform some larger overall calculation. In this article, we present the first high-throughput implementation of a distributed phylogenetics platform, MultiPhyl, capable of using the idle computational resources of many heterogeneous non-dedicated machines to form a phylogenetics supercomputer. MultiPhyl allows a user to upload hundreds or thousands of amino acid or nucleotide alignments simultaneously and perform computationally intensive tasks such as model selection, tree searching and bootstrapping of each of the alignments using many desktop machines. The program implements a set of 88 amino acid models and 56 nucleotide maximum likelihood models and a variety of statistical methods for choosing between alternative models. A MultiPhyl webserver is available for public use at: http://www.cs.nuim.ie/distributed/multiphyl.php.

Supertrees Disentangle the Chimerical Origin of Eukaryotic Genomes

Eukaryotes are traditionally considered to be one of the three natural divisions of the tree of life and the sister group of the Archaebacteria. However, eukaryotic genomes are replete with genes of eubacterial ancestry, and more than 20 mutually incompatible hypotheses have been proposed to account for eukaryote origins. Here we test the predictions of these hypotheses using a novel supertree-based phylogenetic signal-stripping method, and recover supertrees of life based on phylogenies for up to 5,741 single gene families distributed across 185 genomes. Using our signal-stripping method, we show that there are three distinct phylogenetic signals in eukaryotic genomes. In order of strength, these link eukaryotes with the Cyanobacteria, the Proteobacteria, and the Thermoplasmatales, an archaebacterial (euryarchaeotes) group. These signals correspond to distinct symbiotic partners involved in eukaryote evolution: plastids, mitochondria, and the elusive host lineage. According to our whole-genome data, eukaryotes are hardly the sister group of the Archaebacteria, because up to 83% of eukaryotic genes with a prokaryotic homolog have eubacterial, not archaebacterial, origins. The results reject all but two of the current hypotheses for the origin of eukaryotes: those assuming a sulfur-dependent or hydrogen-dependent syntrophy for the origin of mitochondria.

Mitochondrial DNA phylogeography and mating compatibility reveal marked genetic structuring and speciation in the NE Atlantic bryozoan Celleporella hyalina

The marine bryozoan Celleporella hyalina is a species complex composed of many highly divergent and mostly allopatric genetic lineages that are reproductively isolated but share a remarkably similar morphology. One such lineage commonly encrusts macroalgae throughout the NE Atlantic coast. To explore the processes leading to geographical diversification, reproductive isolation and speciation in this taxon, we (i) investigated NE Atlantic C. hyalina mitochondrial DNA phylogeography, and (ii) used breeding trials between geographical isolates to ascertain reproductive isolation. We find that haplotype diversity is geographically variable and there is a strong population structure, with significant isolation by distance. NE Atlantic C. hyalina is structured into two main parapatric lineages that appear to have had independent Pleistocene histories. Range expansions have resulted in two contact zones in Spain and W Ireland. Lineage 1 is found from Ireland to Spain and has low haplotype diversity, with closely related haplotypes, suggesting a recent population expansion into the Irish Sea, S Ireland, S England and Spain. Lineage 2 is found from Iceland to Spain and has high haplotype diversity. Complete reproductive isolation was found between some geographical isolates representing both lineages, whereas it was incomplete or asymmetric between others, suggesting these latter phylogeographical groups probably represent incipient species. The phylogeographical distribution of NE Atlantic C. hyalina does not fall easily into a pattern of southern refugia, and we discuss likely differences between terrestrial and marine system responses to Pleistocene glacial cycles.

Testing the Neutral Fixation of Hetero-Oligomerism in the Archaeal Chaperonin CCT

The evolutionary transition from homo-oligomerism to hetero-oligomerism in multimeric proteins and its contribution to function innovation and organism complexity remain to be investigated. Here, we undertake the challenge of contributing to this theoretical ground by investigating the hetero-oligomerism in the molecular chaperonin cytosolic chaperonin containing tailless complex polypeptide 1 (CCT) from archaea. CCT is amenable to this study because, in contrast to eukaryotic CCTs where sub-functionalization after gene duplication has been taken to completion, archaeal CCTs present no evidence for subunit functional specialization. Our analyses yield additional information to previous reports on archaeal CCT paralogy by identifying new duplication events. Analyses of selective constraints show that amino acid sites from 1 subunit have fixed slightly deleterious mutations at inter-subunit interfaces after gene duplication. These mutations have been followed by compensatory mutations in nearby regions of the same subunit and in the interface contact regions of its paralogous subunit. The strong selective constraints in these regions after speciation support the evolutionary entrapment of CCTs as hetero-oligomers. In addition, our results unveil different evolutionary dynamics depending on the degree of CCT hetero-oligomerism. Archaeal CCT protein complexes comprising 3 distinct classes of subunits present 2 evolutionary processes. First, slightly deleterious and compensatory mutations were fixed neutrally at inter-subunit regions. Second, sub-functionalization may have occurred at substrate-binding and adenosine triphosphate-binding regions after the 2nd gene duplication event took place. CCTs with 2 distinct types of subunits did not present evidence of sub-functionalization. Our results provide the 1st in silico evidence for the neutral fixation of hetero-oligomerism in archaeal CCTs and provide information on the evolution of hetero-oligomerism toward sub-functionalization in archaeal CCTs.

A canonical FtsZ protein in Verrucomicrobium spinosum, a member of the Bacterial phylum Verrucomicrobia that also includes tubulin-producing Prosthecobacter species

Background

The origin and evolution of the homologous GTP-binding cytoskeletal proteins FtsZ typical of Bacteria and tubulin characteristic of eukaryotes is a major question in molecular evolutionary biology. Both FtsZ and tubulin are central to key cell biology processes – bacterial septation and cell division in the case of FtsZ and in the case of tubulins the function of microtubules necessary for mitosis and other key cytoskeleton-dependent processes in eukaryotes. The origin of tubulin in particular is of significance to models for eukaryote origins. Most members of domain Bacteria possess FtsZ, but bacteria in genus Prosthecobacter of the phylum Verrucomicrobia form a key exception, possessing tubulin homologs BtubA and BtubB. It is therefore of interest to know whether other members of phylum Verrucomicrobia possess FtsZ or tubulin as their FtsZ-tubulin gene family representative.

Results

Verrucomicrobium spinosum, a member of Phylum Verrucomicrobia of domain Bacteria, has been found to possess a gene for a protein homologous to the cytoskeletal protein FtsZ. The deduced amino acid sequence has sequence signatures and predicted secondary structure characteristic for FtsZ rather than tubulin, but phylogenetic trees and sequence analysis indicate that it is divergent from all other known FtsZ sequences in members of domain Bacteria. The FtsZ gene of V. spinosum is located within a dcw gene cluster exhibiting gene order conservation known to contribute to the divisome in other Bacteria and comparable to these clusters in other Bacteria, suggesting a similar functional role.

Conclusion

Verrucomicrobium spinosum has been found to possess a gene for a protein homologous to the cytoskeletal protein FtsZ. The results suggest the functional as well as structural homology of the V. spinosum FtsZ to the FtsZs of other Bacteria implying its involvement in cell septum formation during division. Thus, both bacteria-like FtsZ and eukaryote-like tubulin cytoskeletal homologs occur in different species of the phylum Verrucomicrobia of domain Bacteria, a result with potential major implications for understanding evolution of tubulin-like cytoskeletal proteins and the origin of eukaryote tubulins.

Origin and diversification of the basic helix-loop-helix gene family in metazoans: insights from comparative genomics

BACKGROUND

Molecular and genetic analyses conducted in model organisms such as Drosophila and vertebrates, have provided a wealth of information about how networks of transcription factors control the proper development of these species. Much less is known, however, about the evolutionary origin of these elaborated networks and their large-scale evolution. Here we report the first evolutionary analysis of a whole superfamily of transcription factors, the basic helix-loop-helix (bHLH) proteins, at the scale of the whole metazoan kingdom.

RESULTS

We identified in silico the putative full complement of bHLH genes in the sequenced genomes of 12 different species representative of the main metazoan lineages, including three non-bilaterian metazoans, the cnidarians Nematostella vectensis and Hydra magnipapillata and the demosponge Amphimedon queenslandica. We have performed extensive phylogenetic analyses of the 695 identified bHLHs, which has allowed us to allocate most of these bHLHs to defined evolutionary conserved groups of orthology.

CONCLUSION

Three main features in the history of the bHLH gene superfamily can be inferred from these analyses: (i) an initial diversification of the bHLHs has occurred in the pre-Cambrian, prior to metazoan cladogenesis; (ii) a second expansion of the bHLH superfamily occurred early in metazoan evolution before bilaterians and cnidarians diverged; and (iii) the bHLH complement during the evolution of the bilaterians has been remarkably stable. We suggest that these features may be extended to other developmental gene families and reflect a general trend in the evolution of the developmental gene repertoires of metazoans.

Origin and evolution of the mitochondrial aminoacyl-tRNA synthetases

Many theories favor a fusion of 2 prokaryotic genomes for the origin of the Eukaryotes, but there are disagreements on the origin, timing, and cellular structures of the cells involved. Equally controversial is the source of the nuclear genes for mitochondrial proteins, although the alpha-proteobacterial contribution to the mitochondrial genome is well established. Phylogenetic inferences show that the nuclearly encoded mitochondrial aminoacyl-tRNA synthetases (aaRSs) occupy a position in the tree that is not close to any of the currently sequenced alpha-proteobacterial genomes, despite cohesive and remarkably well-resolved alpha-proteobacterial clades in 12 of the 20 trees. Two or more alpha-proteobacterial clusters were observed in 8 cases, indicative of differential loss of paralogous genes or horizontal gene transfer. Replacement and retargeting events within the nuclear genomes of the Eukaryotes was indicated in 10 trees, 4 of which also show split alpha-proteobacterial groups. A majority of the mitochondrial aaRSs originate from within the bacterial domain, but none specifically from the alpha-Proteobacteria. For some aaRS, the endosymbiotic origin may have been erased by ongoing gene replacements on the bacterial as well as the eukaryotic side. For others that accurately resolve the alpha-proteobacterial divergence patterns, the lack of affiliation with mitochondria is more surprising. We hypothesize that the ancestral eukaryotic gene pool hosted primordial "bacterial-like" genes, to which a limited set of alpha-proteobacterial genes, mostly coding for components of the respiratory chain complexes, were added and selectively maintained.

A fungal phylogeny based on 42 complete genomes derived from supertree and combined gene analysis

BACKGROUND

To date, most fungal phylogenies have been derived from single gene comparisons, or from concatenated alignments of a small number of genes. The increase in fungal genome sequencing presents an opportunity to reconstruct evolutionary events using entire genomes. As a tool for future comparative, phylogenomic and phylogenetic studies, we used both supertrees and concatenated alignments to infer relationships between 42 species of fungi for which complete genome sequences are available.

RESULTS

A dataset of 345,829 genes was extracted from 42 publicly available fungal genomes. Supertree methods were employed to derive phylogenies from 4,805 single gene families. We found that the average consensus supertree method may suffer from long-branch attraction artifacts, while matrix representation with parsimony (MRP) appears to be immune from these. A genome phylogeny was also reconstructed from a concatenated alignment of 153 universally distributed orthologs. Our MRP supertree and concatenated phylogeny are highly congruent. Within the Ascomycota, the sub-phyla Pezizomycotina and Saccharomycotina were resolved. Both phylogenies infer that the Leotiomycetes are the closest sister group to the Sordariomycetes. There is some ambiguity regarding the placement of Stagonospora nodurum, the sole member of the class Dothideomycetes present in the dataset. Within the Saccharomycotina, a monophyletic clade containing organisms that translate CTG as serine instead of leucine is evident. There is also strong support for two groups within the CTG clade, one containing the fully sexual species Candida lusitaniae, Candida guilliermondii and Debaryomyces hansenii, and the second group containing Candida albicans, Candida dubliniensis, Candida tropicalis, Candida parapsilosis and Lodderomyces elongisporus. The second major clade within the Saccharomycotina contains species whose genomes have undergone a whole genome duplication (WGD), and their close relatives. We could not confidently resolve whether Candida glabrata or Saccharomyces castellii lies at the base of the WGD clade.

CONCLUSION

We have constructed robust phylogenies for fungi based on whole genome analysis. Overall, our phylogenies provide strong support for the classification of phyla, sub-phyla, classes and orders. We have resolved the relationship of the classes Leotiomyctes and Sordariomycetes, and have identified two classes within the CTG clade of the Saccharomycotina that may correlate with sexual status.

Divergence of the Dof gene families in poplar, Arabidopsis, and rice suggests multiple modes of gene evolution after duplication

It is widely accepted that gene duplication is a primary source of genetic novelty. However, the evolutionary fate of duplicated genes remains largely unresolved. The classical Ohno's Duplication-Retention-Non/Neofunctionalization theory, and the recently proposed alternatives such as subfunctionalization or duplication-degeneration-complementation, and subneofunctionalization, each can explain one or more aspects of gene fate after duplication. Duplicated genes are also affected by epigenetic changes. We constructed a phylogenetic tree using Dof (DNA binding with one finger) protein sequences from poplar (Populus trichocarpa) Torr. & Gray ex Brayshaw, Arabidopsis (Arabidopsis thaliana), and rice (Oryza sativa). From the phylogenetic tree, we identified 27 pairs of paralogous Dof genes in the terminal nodes. Analysis of protein motif structure of the Dof paralogs and their ancestors revealed six different gene fates after gene duplication. Differential protein methylation was revealed between a pair of duplicated poplar Dof genes, which have identical motif structure and similar expression pattern, indicating that epigenetics is involved in evolution. Analysis of reverse transcription-PCR, massively parallel signature sequencing, and microarray data revealed that the paralogs differ in expression pattern. Furthermore, analysis of nonsynonymous and synonymous substitution rates indicated that divergence of the duplicated genes was driven by positive selection. About one-half of the motifs in Dof proteins were shared by non-Dof proteins in the three plants species, indicating that motif co-option may be one of the forces driving gene diversification. We provided evidence that the Ohno's Duplication-Retention-Non/Neofunctionalization, subfunctionalization/duplication-degeneration-complementation, and subneofunctionalization hypotheses are complementary with, not alternative to, each other.

Evolution of filamentous plant pathogens: gene exchange across eukaryotic kingdoms

Filamentous fungi and oomycetes are eukaryotic microorganisms that grow by producing networks of thread-like hyphae, which secrete enzymes to break down complex nutrients, such as wood and plant material, and recover the resulting simple sugars and amino acids by osmotrophy. These organisms are extremely similar in both appearance and lifestyle and include some of the most economically important plant pathogens . However, the morphological similarity of fungi and oomycetes is misleading because they represent some of the most distantly related eukaryote evolutionary groupings, and their shared osmotrophic growth habit is interpreted as being the result of convergent evolution . The fungi branch with the animals, whereas the oomycetes branch with photosynthetic algae as part of the Chromalveolata . In this report, we provide strong phylogenetic evidence that multiple horizontal gene transfers (HGT) have occurred from filamentous ascomycete fungi to the distantly related oomycetes. We also present evidence that a subset of the associated gene families was initially the product of prokaryote-to-fungi HGT. The predicted functions of the gene products associated with fungi-to-oomycete HGT suggest that this process has played a significant role in the evolution of the osmotrophic, filamentous lifestyle on two separate branches of the eukaryote tree.

Identification of rDNA-specific non-LTR retrotransposons in Cnidaria

Ribosomal RNA genes are abundant repetitive sequences in most eukaryotes. Ribosomal DNA (rDNA) contains many insertions derived from mobile elements including non-long terminal repeat (non-LTR) retrotransposons. R2 is the well-characterized 28S rDNA-specific non-LTR retrotransposon family that is distributed over at least 4 bilaterian phyla. R2 is a large family sharing the same insertion specificity and classified into 4 clades (R2-A, -B, -C, and -D) based on the N-terminal domain structure and the phylogeny. There is no observation of horizontal transfer of R2; therefore, the origin of R2 dates back to before the split between protostomes and deuterostomes. Here, we in silico identified 1 R2 element from the sea anemone Nematostella vectensis and 2 R2-like retrotransposons from the hydrozoan Hydra magnipapillata. R2 from N. vectensis was inserted into the 28S rDNA like other R2, but the R2-like elements from H. magnipapillata were inserted into the specific sequence in the highly conserved region of the 18S rDNA. We designated the Hydra R2-like elements R8. R8 is inserted at 37 bp upstream from R7, another 18S rDNA-specific retrotransposon family. There is no obvious sequence similarity between targets of R2 and R8, probably because they recognize long DNA sequences. Domain structure and phylogeny indicate that R2 from N. vectensis is the member of the R2-D clade, and R8 from H. magnipapillata belongs to the R2-A clade despite its different sequence specificity. These results suggest that R2 had been generated before the split between cnidarians and bilaterians and that R8 is a retrotransposon family that changed its target from the 28S rDNA to the 18S rDNA.

Evolution of four gene families with patchy phylogenetic distributions: influx of genes into protist genomes

BACKGROUND

Lateral gene transfer (LGT) in eukaryotes from non-organellar sources is a controversial subject in need of further study. Here we present gene distribution and phylogenetic analyses of the genes encoding the hybrid-cluster protein, A-type flavoprotein, glucosamine-6-phosphate isomerase, and alcohol dehydrogenase E. These four genes have a limited distribution among sequenced prokaryotic and eukaryotic genomes and were previously implicated in gene transfer events affecting eukaryotes. If our previous contention that these genes were introduced by LGT independently into the diplomonad and Entamoeba lineages were true, we expect that the number of putative transfers and the phylogenetic signal supporting LGT should be stable or increase, rather than decrease, when novel eukaryotic and prokaryotic homologs are added to the analyses.

RESULTS

The addition of homologs from phagotrophic protists, including several Entamoeba species, the pelobiont Mastigamoeba balamuthi, and the parabasalid Trichomonas vaginalis, and a large quantity of sequences from genome projects resulted in an apparent increase in the number of putative transfer events affecting all three domains of life. Some of the eukaryotic transfers affect a wide range of protists, such as three divergent lineages of Amoebozoa, represented by Entamoeba, Mastigamoeba, and Dictyostelium, while other transfers only affect a limited diversity, for example only the Entamoeba lineage. These observations are consistent with a model where these genes have been introduced into protist genomes independently from various sources over a long evolutionary time.

CONCLUSION

Phylogenetic analyses of the updated datasets using more sophisticated phylogenetic methods, in combination with the gene distribution analyses, strengthened, rather than weakened, the support for LGT as an important mechanism affecting the evolution of these gene families. Thus, gene transfer seems to be an on-going evolutionary mechanism by which genes are spread between unrelated lineages of all three domains of life, further indicating the importance of LGT from non-organellar sources into eukaryotic genomes.