Phylogenetic and phyletic studies of informational genes in genomes highlight existence of a 4 domain of life including giant viruses

The discovery of Mimivirus, with its very large genome content, made it possible to identify genes common to the three domains of life (Eukarya, Bacteria and Archaea) and to generate controversial phylogenomic trees congruent with that of ribosomal genes, branching Mimivirus at its root. Here we used sequences from metagenomic databases, Marseillevirus and three new viruses extending the Mimiviridae family to generate the phylogenetic trees of eight proteins involved in different steps of DNA processing. Compared to the three ribosomal defined domains, we report a single common origin for Nucleocytoplasmic Large DNA Viruses (NCLDV), DNA processing genes rooted between Archaea and Eukarya, with a topology congruent with that of the ribosomal tree. As for translation, we found in our new viruses, together with Mimivirus, five proteins rooted deeply in the eukaryotic clade. In addition, comparison of informational genes repertoire based on phyletic pattern analysis supports existence of a clade containing NCLDVs clearly distinct from that of Eukarya, Bacteria and Archaea. We hypothesize that the core genome of NCLDV is as ancient as the three currently accepted domains of life.     

Sorting through the chaff, nDNA gene trees for phylogenetic inference and hybrid identification of annual sunflowers (Helianthus sect. Helianthus)

The annual sunflowers (Helianthus sect. Helianthus) present a formidable challenge for phylogenetic inference because of ancient hybrid speciation, recent introgression, and suspected issues with deep coalescence. Here we analyze sequence data from 11 nuclear DNA (nDNA) genes for multiple genotypes of species within the section to (1) reconstruct the phylogeny of this group, (2) explore the utility of nDNA gene trees for detecting hybrid speciation and introgression; and (3) test an empirical method of hybrid identification based on the phylogenetic congruence of nDNA gene trees from tightly linked genes. We uncovered considerable topological heterogeneity among gene trees with or without three previously identified hybrid species included in the analyses, as well as a general lack of reciprocal monophyly of species. Nonetheless, partitioned Bayesian analyses provided strong support for the reciprocal monophyly of all species except H. annuus (0.89 PP), the most widespread and abundant annual sunflower. Previous hypotheses of relationships among taxa were generally strongly supported (1.0 PP), except among taxa typically associated with H. annuus, apparently due to the paraphyly of the latter in all gene trees. While the individual nDNA gene trees provided a useful means for detecting recent hybridization, identification of ancient hybridization was problematic for all ancient hybrid species, even when linkage was considered. We discuss biological factors that affect the efficacy of phylogenetic methods for hybrid identification.     

Weighted genome trees: refinements and applications

There are many ways to group completed genome sequences in hierarchical patterns (trees) reflecting relationships between their genes. Such groupings help us organize biological information and bear crucially on underlying processes of genome and organismal evolution. Genome trees make use of all comparable genes but can variously weight the contributions of these genes according to similarity, congruent patterns of similarity, or prevalence among genomes. Here we explore such possible weighting strategies, in an analysis of 142 prokaryotic and 5 eukaryotic genomes. We demonstrate that alternate weighting strategies have different advantages, and we propose that each may have its specific uses in systematic or evolutionary biology. Comparisons of results obtained with different methods can provide further clues to major events and processes in genome evolution.     

Genomic basis for cell-wall diversity in plants. A comparative approach to gene families in rice and Arabidopsis

Monocotyledons and dicotyledons are distinct, not only in their body plans and developmental patterns, but also in the structural features of their cell walls. The recent completion of the rice (Oryza sativa) genomic sequence and publication of the sequence data, together with the completed database of the Arabidopsis thaliana genome, provide the first opportunity to compare the full complement of cell-wall-related genes from the two distinct classes of flowering plants. We made this comparison by exploiting the fact that Arabidopsis and rice have type I and type II walls, respectively, and therefore represent the two extremes in terms of the structural features of plant cell walls. In this review article, we classify all cell-wall-related genes into 32 gene families, and generate their phylogenetic trees. Using these data, we can phylogenetically compare individual genes of particular interest between Arabidopsis and rice. This comparative genome approach shows that the differences in wall architecture in the two plant groups actually mirror the diversity of the individual gene families involved in the cell-wall dynamics of the respective plant species. This study also identifies putative rice orthologs of genes with well-defined functions in Arabidopsis and other plant species.     

Alpha‐Tubulin and Small Subunit rRNA Phylogenies of Peritrichs Are Congruent and Do Not Support the Clustering of Mobilids and Sessilids (Ciliophora,Oligohymenophorea)

ABSTRACT. Peritrich ciliates have been traditionally subdivided into two orders, Sessilida and Mobilida within the subclass Peritrichia. However, all the existing small subunit (SSU) rRNA phylogenetic trees showed that the sessilids and mobilids did not branch together. To shed some light on this disagreement, we tested whether or not the classic Peritrichia is a monophyletic group by assessing the reliability of the SSU rRNA phylogeny in terms of congruency with α‐tubulin phylogeny. For this purpose, we obtained 10 partial α‐tubulin sequences from peritrichs and built phylogenetic trees based on α‐tubulin nucleotide and amino acid data. A phylogenetic tree from the α‐tubulin and SSU rRNA genes in combination was also constructed and compared with that from the SSU rRNA gene using a similar species sampling. Our results show that the mobilids and sessilids are consistently separated in all trees, which reinforces the idea that the peritrichs do not constitute a monophyletic group. However, in all α‐tubulin gene trees, the urceolariids and trichodiniids do not group together, suggested mobilids may not be a monophyletic group.     

Phylo-MCOA: a fast and efficient method to detect outlier genes and species in phylogenomics using multiple co-inertia analysis

Full genome data sets are currently being explored on a regular basis to infer phylogenetic trees, but there are often discordances among the trees produced by different genes. An important goal in phylogenomics is to identify which individual gene and species produce the same phylogenetic tree and are thus likely to share the same evolutionary history. On the other hand, it is also essential to identify which genes and species produce discordant topologies and therefore evolve in a different way or represent noise in the data. The latter are outlier genes or species and they can provide a wealth of information on potentially interesting biological processes, such as incomplete lineage sorting, hybridization, and horizontal gene transfers. Here, we propose a new method to explore the genomic tree space and detect outlier genes and species based on multiple co-inertia analysis (MCOA), which efficiently captures and compares the similarities in the phylogenetic topologies produced by individual genes. Our method allows the rapid identification of outlier genes and species by extracting the similarities and discrepancies, in terms of the pairwise distances, between all the species in all the trees, simultaneously. This is achieved by using MCOA, which finds successive decomposition axes from individual ordinations (i.e., derived from distance matrices) that maximize a covariance function. The method is freely available as a set of R functions. The source code and tutorial can be found online at http://phylomcoa.cgenomics.org.     

Phylogenomic analysis resolves the interordinal relationships and rapid diversification of the laurasiatherian mammals

Although great progress has been made in resolving the relationships of placental mammals, the position of several clades in Laurasiatheria remain controversial. In this study, we performed a phylogenetic analysis of 97 orthologs (46,152 bp) for 15 taxa, representing all laurasiatherian orders. Additionally, phylogenetic trees of laurasiatherian mammals with draft genome sequences were reconstructed based on 1608 exons (2,175,102 bp). Our reconstructions resolve the interordinal relationships within Laurasiatheria and corroborate the clades Scrotifera, Fereuungulata, and Cetartiodactyla. Furthermore, we tested alternative topologies within Laurasiatheria, and among alternatives for the phylogenetic position of Perissodactyla, a sister-group relationship with Cetartiodactyla receives the highest support. Thus, Pegasoferae (Perissodactyla + Carnivora + Pholidota + Chiroptera) does not appear to be a natural group. Divergence time estimates from these genes were compared with published estimates for splits within Laurasiatheria. Our estimates were similar to those of several studies and suggest that the divergences among these orders occurred within just a few million years.     

Lateral gene transfer and phylogenetic assignment of environmental fosmid clones

Metagenomic data, especially sequence data from large insert clones, are most useful when reasonable inferences about phylogenetic origins of inserts can be made. Often, clones that bear phylotypic markers (usually ribosomal RNA genes) are sought, but sometimes phylogenetic assignments have been based on the preponderance of blast hits obtained with predicted protein coding sequences (CDSs). Here we use a cloning method which greatly enriches for ribosomal RNA-bearing fosmid clones to ask two questions: (i) how reliably can we judge the phylogenetic origin of a clone (that is, its RNA phylotype) from the sequences of its CDSs? and (ii) how much lateral gene transfer (LGT) do we see, as assessed by CDSs of different phylogenetic origins on the same fosmid? We sequenced 12 rRNA containing fosmid clones, obtained from libraries constructed using DNA isolated from Baltimore harbour sediments. Three of the clones are from bacterial candidate divisions for which no cultured representatives are available, and thus represent the first protein coding sequences from these major bacterial lineages. The amount of LGT was assessed by making phylogenetic trees of all the CDSs in the fosmid clones and comparing the phylogenetic position of the CDS to the rRNA phylotype. We find that the majority of CDSs in each fosmid, 57-96%, agree with their respective rRNA genes. However, we also find that a significant fraction of the CDSs in each fosmid, 7-44%, has been acquired by LGT. In several cases, we can infer co-transfer of functionally related genes, and generate hypotheses about mechanism and ecological significance of transfer.     

Geographic structure in metabolome and herbivore community co‐occurs with genetic structure in plant defence genes

Plant–herbivore interactions vary across the landscape and have been hypothesised to promote local adaption in plants to the prevailing herbivore regime. Herbivores that feed on European aspen (Populus tremula) change across regional scales and selection on host defence genes may thus change at comparable scales. We have previously observed strong population differentiation in a set of inducible defence genes in Swedish P. tremula. Here, we study the geographic patterns of abundance and diversity of herbivorous insects, the untargeted metabolome of the foliage and genetic variation in a set of wound‐induced genes and show that the geographic structure co‐occurs in all three data sets. In response to this structure, we observe local maladaptation of herbivores, with fewer herbivores on local trees than on trees originated from more distant localities. Finally, we also identify 28 significant associations between single nucleotide polymorphisms SNPs from defence genes and a number of the herbivore traits and metabolic profiles.     

Dealing with saturation at the amino acid level: a case study based on anciently duplicated zebrafish genes

The ray-finned fishes (Actinopterygii) seem to have two copies of many tetrapod (Sarcopterygii) genes. The origin of these duplicate fish genes is the subject of some controversy. One explanation for the existence of these extra fish genes could be an increase in the rate of independent gene duplications in fishes. Alternatively, gene duplicates in fish may have been formed in the ancestor of all or most Actinopterygii during a complete genome duplication event. A third possibility is that tetrapods have lost more genes than fish after gene or genome duplication events in the common ancestor of both lineages. These three hypotheses can be tested by phylogenetic reconstruction. Previously, we found that a large number of anciently duplicated genes of zebrafish are sister sequences in evolutionary trees suggesting that they were produced in Actinopterygii after the divergence of Sarcopterygii [Phil. Trans. R. Soc. Lond. B 356 (2001) 119]. On the other hand, several well-supported trees showed one of the two fish genes as the sister sequence to a monophyletic clade that included the second fish gene and genes from frog, chicken, mouse and human. These so-called outgroup topologies suggest that the origin of many fish duplicates predates the divergence of the Sarcopterygii and Actinopterygii and support the hypothesis that tetrapods have lost duplicates that have been retained in fish. Here we show that many of these 'outgroup' tree topologies are erroneous and can be corrected when mutational saturation is taken into account. To this end, a Java-based application has been developed to visualize the amount of saturation in amino acid sequences. The program graphically displays the number of observed frequent and rare amino acid replacements between pairs of sequences against their overall evolutionary distance. Discrimination between frequent and rare amino acid replacements is based on substitution probability matrices (e.g. PAM and BLOSUM). Evolutionary distances between sequences can be computed from the fraction of unsaturated sites only and evolutionary trees inferred by pairwise distance methods. When trees are computed by omitting the saturated fraction of sites, most fish duplicates are sister sequences.     

Nuclear and mitochondrial genes contain similar phylogenetic signal for pigeons and doves (Aves: Columbiformes)

Molecular systematic studies generally assume that gene trees are reasonable estimates of species trees. We tested the validity of this assumption in the pigeons and doves (Aves: Columbiformes) by comparing phylogenies derived from nuclear (beta-fibrinogen intron 7) and mitochondrial (cytochrome b) genes. Trees derived from the two genes when analyzed separately contained many nodes in common. A partition homogeneity test revealed no significant incongruence between trees derived from the two genes; so, we combined nuclear and mitochondrial data in subsequent phylogenetic analyses. The resulting tree, which was highly resolved and generally well supported, contained a strong biogeographic component. The rate of nucleotide substitution for the nuclear intron was approximately six times slower than that of cytochrome b. This resulted in a much higher consistency index for trees derived from the intron because of the low level of multiple substitution. However, the degree of resolution and support for trees reconstructed from the two genes was similar. We also examined the transition and transversion substitution rates for the genes. Third position transversions for cytochrome b accumulated linearly with intron divergence, suggesting low levels of multiple substitution for third position transversions.     

Estimation of phylogenetic inconsistencies in the three domains of life

Discrepancies in phylogenetic trees of bacteria and archaea are often explained as lateral gene transfer events. However, such discrepancies may also be due to phylogenetic artifacts or orthology assignment problems. A first step that may help to resolve this dilemma is to estimate the extent of phylogenetic inconsistencies in trees of prokaryotes in comparison with those of higher eukaryotes, where no lateral gene transfer is expected. To test this, we used 21 proteomes each of eukaryotes (mainly opisthokonts), proteobacteria, and archaea that spanned equivalent levels of genetic divergence. In each domain of life, we defined a set of putative orthologous sequences using a phylogenetic-based orthology protocol and, as a reference topology, we used a tree constructed with concatenated genes of each domain. Our results show, for most of the tests performed, that the magnitude of topological inconsistencies with respect to the reference tree was very similar in the trees of proteobacteria and eukaryotes. When clade support was taken into account, prokaryotes showed some more inconsistencies, but then all values were very low. Discrepancies were only consistently higher in archaea but, as shown by simulation analysis, this is likely due to the particular tree of the archaeal species used here being more difficult to reconstruct, whereas the trees of proteobacteria and eukaryotes were of similar difficulty. Although these results are based on a relatively small number of genes, it seems that phylogenetic reconstruction problems, including orthology assignment problems, have a similar overall effect over prokaryotic and eukaryotic trees based on single genes. Consequently, lateral gene transfer between distant prokaryotic species may have been more rare than previously thought, which opens the way to obtain the tree of life of bacterial and archaeal species using genomic data and the concatenation of adequate genes, in the same way as it is usually done in eukaryotes.     

Defining the Core of Nontransferable Prokaryotic Genes: The Euryarchaeal Core

If lateral gene transfer (LGT) has affected all genes over the course of prokaryotic evolution, reconstruction of organismal phylogeny is compromised. However, if a core of genes is immune to transfer, then the evolutionary history of that core might be our most reliable guide to the evolution of organisms. Such a core should be preferentially included in the subset of genes shared by all organisms, but where universally conserved genes have been analyzed, there is too little phylogenetic signal to allow determination of whether or not they indeed have the same history (Hansmann and Martin 2000; Teichmann and Mitchison 1999). Here we look at a more restricted set, 521 homologous genes (COGs) simultaneously present in four sequenced euryarchaeal genomes. Although there is overall little robust phylogenetic signal in this data set, there is, among well-supported trees, strong representation of all three possible four-taxon topologies. ``Informational' genes seem no less subject to LGT than are ``operational genes,' within the euryarchaeotes. We conclude that (i) even in this collection of conserved genes there has been extensive LGT (orthologous gene replacement) and (ii) the notion that there is a core of nontransferable genes (the ``core hypothesis') has not been proven and may be unprovable. Received: 7 November 2000 / Accepted: 20 February 2001     

Linking dendroecology and association genetics in natural populations: Stress responses archived in tree rings associate with SNP genotypes in silver fir (Abies alba Mill.)

Genetic association studies in forest trees would greatly benefit from information on the response of trees to environmental stressors over time, which can be provided by dendroecological analysis. Here, we jointly analysed dendroecological and genetic data of surviving silver fir trees to explore the genetic basis of their response to the iconic stress episode of the 1970s and 1980s that led to large‐scale forest dieback in Central Europe and has been attributed to air pollution. Specifically, we derived dendrophenotypic measures from 190 trees in the Bavarian Forest that characterize the resistance, resilience and recovery during this growth depression, and in the drought year in 1976. By focusing on relative growth changes of trees and by standardizing the dendrophenotypes within stands, we accounted for variation introduced by micro‐ and macroscale environmental differences. We associated the dendrophenotypes with single nucleotide polymorphisms (SNPs) in candidate genes using general linear models (GLMs) and the machine learning algorithm random forest with subsequent feature selection. Most trees at our study sites experienced a severe growth decline from 1974 until the mid‐1980s with minimum values during the drought year. Fifteen genes were associated with the dendrophenotypes, including genes linked to photosynthesis and drought stress. With our study, we show that dendrophenotypes can be a powerful resource for genetic association studies that permit to account for micro‐ and macroenvironmental variation when data are derived from natural populations. We call for a wider collaboration of dendroecologists and forest geneticists to integrate individual tree‐level dendrophenotypes in genetic association studies.     

Restoration of threatened species: a noble cause for transgenic trees

Some of the first applications of transgenic trees in North America may be for the conservation or restoration of threatened forest trees that have been devastated by fungal pathogens or insect pests. In some cases, where resistance has yet to be found in the natural population of a tree species, incorporating genes from other organisms may offer the only hope for restoration. In others, transgenics may play a role as part of an integrated approach, along with conventional breeding or biocontrol agents. American chestnut (Castanea dentata) was wiped out as a canopy species by a fungal disease accidentally introduced into the United States around 1900. Similarly, American elm (Ulmus americana) virtually disappeared as a favored street tree from Northeastern U.S. cities after the introduction of the Dutch elm disease fungus in the 1940s. In both cases, progress has been made toward restoration via conventional techniques such as selection and propagation of tolerant cultivars (American elm) or breeding with a related resistant species (American chestnut). Recently, progress has also been made with development of systems for engineering antifungal candidate genes into these “heritage trees.” An Agrobacterium-leaf disk system has been used to produce transgenic American elm trees engineered with an antimicrobial peptide gene that may enhance resistance to Dutch elm disease. Two gene transfer systems have been developed for American chestnut using Agrobacterium-mediated transformation of embryogenic cultures, setting the stage for the first tests of potential antifungal genes for their ability to confer resistance to the chestnut blight fungus. Despite the promise of transgenic approaches for restoration of these heritage trees, a number of technical, environmental, economic, and ethical questions remain to be addressed before such trees can be deployed, and the debate around these questions may be quite different from that associated with transgenic trees developed for other purposes.     

Base-compositional biases and the bat problem. II. DNA-hybridization trees based on AT- and GC-enriched tracers.

We conducted a series of parallel DNA-hybridization experiments on a small group of bats (species of Pteropus, Rhinolophus, Noctilio and Pteronotus) and outgroups (Lemur, Cynocephalus, Didelphis), using whole-genome labels and tracers made from extracts enriched with AT and two levels of GC content. FITCH (additive phylogenetic trees) topologies were constructed from the four sets of comparisons, indexed as both delta Tmode and delta NPHs (normalized percentage of hybridization). Based on our previous work showing that the shared AT bias of pteropodids and some microchiropterans may affect the rank-ordering of taxa based on either AT- or GC-rich labels, our expectation was that the resulting trees would show differing topologies when generated from tracers made with the variously enriched DNA extracts. Whereas there was some variation among the trees, most of them grouped the bats together, and almost all paired the representative megachiropteran and rhinolophoid microchiropteran as sister-taxa in contrast to the other microchiropterans. As the pteropodid-rhinolophoid relationship is an unexpected and unlikely one, we attribute this association to an AT bias that was not obviated even by our most GC-rich labels, and suggest that such a bias may compromise the truth of some molecular trees. Accordingly, we believe the broader issue of bat monophyly remains unresolved by DNA-hybridization and probably also by gene-sequencing studies.     

Gene expression and metabolite profiling of Populus euphratica growing in the Negev desert

Background

Plants growing in their natural habitat represent a valuable resource for elucidating mechanisms of acclimation to environmental constraints. Populus euphratica is a salt-tolerant tree species growing in saline semi-arid areas. To identify genes involved in abiotic stress responses under natural conditions we constructed several normalized and subtracted cDNA libraries from control, stress-exposed and desert-grown P. euphratica trees. In addition, we identified several metabolites in desert-grown P. euphratica trees.

Results

About 14,000 expressed sequence tag (EST) sequences were obtained with a good representation of genes putatively involved in resistance and tolerance to salt and other abiotic stresses. A P. euphratica DNA microarray with a uni-gene set of ESTs representing approximately 6,340 different genes was constructed. The microarray was used to study gene expression in adult P. euphratica trees growing in the desert canyon of Ein Avdat in Israel. In parallel, 22 selected metabolites were profiled in the same trees.

Conclusion

Of the obtained ESTs, 98% were found in the sequenced P. trichocarpa genome and 74% in other Populus EST collections. This implies that the P. euphratica genome does not contain different genes per se, but that regulation of gene expression might be different and that P. euphratica expresses a different set of genes that contribute to adaptation to saline growth conditions. Also, all of the five measured amino acids show increased levels in trees growing in the more saline soil.     

Molecular phylogenetic relationships based on mitochondrial and nuclear gene sequences for the Todies (Todus, Todidae) of the Caribbean

We used mitochondrial/nuclear gene sequence analyses to determine the historical relationships of the endemic species of Todus (Aves: Todidae) from the Caribbean. We collected 1920-bp of nucleotide sequence data from the mitochondrial genes cytochrome b, ATPase 6, ATPase 8, and 591-bp of the single-copy nuclear gene c-mos for all Todus species and representatives of their outgroup taxa (Hylomanes, Barypthengus, Chloroceryle, Ceryle, and Galbula) to reconstruct the evolutionary history (via parsimony and maximum likelihood) of the five Todus species. The substitution rates among the mitochondrial genes were found to be much higher than the substitution rate for the c-mos gene, consequently resulting in higher substitutional saturation for the mitochondrial genes. When we applied weighting schemes to account for the variance in substitutional heterogeneity among the genes then parsimony and likelihood analyses both demonstrate that the genus Todus is monophyletic and closer to the Hylomanes and Barypthengus genera than the Chloroceryle and Ceryle genera. The mitochondrial-gene trees and nuclear-gene trees both show similar results, thus providing support for the relationships among the taxa from loci within two independently evolving genomes. The nuclear gene c-mos was found, therefore, to be a viable nuclear gene candidate for resolving intermediate and deep divergences.     

Girdling affects carbohydrate-related gene expression in leaves,bark and roots of alternate-bearing citrus trees

Effects of girdling on carbohydrate status and carbohydrate-related gene expression in citrus trees were investigated. Alternate-bearing 'Murcott' (a Citrus reticulata hybrid of unknown origin) trees were girdled during autumn (25 Sep. 2001) and examined 10 weeks later. Girdling brought about carbohydrate (soluble sugar and starch) accumulation in leaves and shoot bark above the girdle, in trees during their fruitless, 'off' year. Trees during their heavy fruit load, 'on' year did not accumulate carbohydrates above the girdle due to the high demand for carbohydrates by the developing fruit. Girdling caused a strong decline in soluble sugar and starch concentrations in organs below the girdle (roots), in both 'on' and 'off' trees. Expression of STPH-L and STPH-H (two isoforms of starch phosphorylase), Agps (ADP-glucose pyrophosphorylase, small subunit), AATP (plastidic ADP/ATP transporter), PGM-C (phosphoglucomutase) and CitSuS1 (sucrose synthase), all of which are associated with starch accumulation, was studied. It was found that gene expression is related to starch accumulation in all 'off' tree organs. RNA levels of all the genes examined were high in leaves and bark that accumulated high concentrations of starch, and low in roots with declining starch concentrations. It may be hypothesized that changes in specific sugars signal the up- and down-regulation of genes involved in starch synthesis.