Duplicate genes and the root of angiosperms, with an example using phytochrome sequences

STRIDE: Species Tree Root Inference from Gene Duplication Events

The correct interpretation of any phylogenetic tree is dependent on that tree being correctly rooted. We present STRIDE, a fast, effective, and outgroup-free method for identification of gene duplication events and species tree root inference in large-scale molecular phylogenetic analyses. STRIDE identifies sets of well-supported in-group gene duplication events from a set of unrooted gene trees, and analyses these events to infer a probability distribution over an unrooted species tree for the location of its root. We show that STRIDE correctly identifies the root of the species tree in multiple large-scale molecular phylogenetic data sets spanning a wide range of timescales and taxonomic groups. We demonstrate that the novel probability model implemented in STRIDE can accurately represent the ambiguity in species tree root assignment for data sets where information is limited. Furthermore, application of STRIDE to outgroup-free inference of the origin of the eukaryotic tree resulted in a root probability distribution that provides additional support for leading hypotheses for the origin of the eukaryotes.

A duplicate gene rooting of seed plants and the phylogenetic position of flowering plants

Flowering plants represent the most significant branch in the tree of land plants, with respect to the number of extant species, their impact on the shaping of modern ecosystems and their economic importance. However, unlike so many persistent phylogenetic problems that have yielded to insights from DNA sequence data, the mystery surrounding the origin of angiosperms has deepened with the advent and advance of molecular systematics. Strong statistical support for competing hypotheses and recent novel trees from molecular data suggest that the accuracy of current molecular trees requires further testing. Analyses of phytochrome amino acids using a duplicate gene-rooting approach yield trees that unite cycads and angiosperms in a clade that is sister to a clade in which Gingko and Cupressophyta are successive sister taxa to gnetophytes plus Pinaceae. Application of a cycads + angiosperms backbone constraint in analyses of a morphological dataset yields better resolved trees than do analyses in which extant gymnosperms are forced to be monophyletic. The results have implications both for our assessment of uncertainty in trees from sequence data and for our use of molecular constraints as a way to integrate insights from morphological and molecular evidence.

Phytochromes differentially regulate seed germination responses to light quality and temperature cues during seed maturation

The ratio of red to far-red light (R : FR) experienced by seeds during maturation affects germination, but the genetic regulation of this effect is poorly understood. In Arabidopsis thaliana, responses to R : FR are governed by five phytochrome photoreceptors, PHYA-PHYE. PHYA, PHYB and PHYE mediate germination, but their roles in germination response to the seed maturation environment are largely unknown. Seeds of A. thaliana phytochrome mutants and natural accessions were matured in a factorial combination of cold (16 degrees C) and warm (24 degrees C) temperatures and high (R : FR = 1) and low (R : FR = 0.6) R : FR environments, resembling sunlight and foliar shade, respectively. Germination was observed in resulting seeds. All five phytochromes mediated germination responses to seed maturation temperature and/or R : FR environment. PHYA suppressed germination in seeds matured under cold temperature, and PHYB promoted germination under the same conditions. PHYD and PHYE promoted germination of seeds matured under warm temperature, but this effect diminished when seeds matured under reduced R : FR. The A. thaliana natural accessions exhibited interesting variation in germination responses to the experimental conditions. Our results suggest that the role of individual PHY loci in regulating plant responses to R : FR varies depending on temperature and provide novel insights into the genetic basis of maternal effects.

Phytochrome gene expression and phylogenetic analysis in the short-day plant Pharbitis nil (Convolvulaceae): Differential regulation by light and an endogenous clock

To investigate the role of distinct phytochrome pools in photoperiodic timekeeping, we characterized four phytochrome genes in the short-day plant Pharbitis nil. Each PHY gene had different photosensory properties and sensitivity to night break that inhibits flowering. During extended dark periods, PHYE, PHYB, and PHYC mRNA accumulation exhibited a circadian rhythmicity indicative of control by an endogenous clock. Phylogenetic analysis recovered four clades of angiosperm phytochrome genes, phyA, phyB, phyC, and phyE. All except the phyE clade included sequences from both monocots and eudicots. In addition, phyA is sister to phyC and phyE sister to phyB, with gymnosperm sequences sister to either the phyA-phyC clade or to the phyB-phyE clade. These results suggest that a single duplication occurred in an ancestral seed plant before the divergence of extant gymnosperms from angiosperms and that two subsequent duplications occurred in an ancestral angiosperm before the divergence of monocots from eudicots. Thus in P. nil, a multigene family with different patterns of mRNA abundance in light and darkness contributes to the total phytochrome pool: one pool is light labile (phyA), whereas the other is light stable (phyB and phyE). In addition, PHYC mRNA represents a third phytochrome pool with intermediate photosensory properties.

Complete plastid genome sequences of Drimys, Liriodendron, and Piper: implications for the phylogenetic relationships of magnoliids

Background

The magnoliids with four orders, 19 families, and 8,500 species represent one of the largest clades of early diverging angiosperms. Although several recent angiosperm phylogenetic analyses supported the monophyly of magnoliids and suggested relationships among the orders, the limited number of genes examined resulted in only weak support, and these issues remain controversial. Furthermore, considerable incongruence resulted in phylogenetic reconstructions supporting three different sets of relationships among magnoliids and the two large angiosperm clades, monocots and eudicots. We sequenced the plastid genomes of three magnoliids, Drimys (Canellales), Liriodendron (Magnoliales), and Piper (Piperales), and used these data in combination with 32 other angiosperm plastid genomes to assess phylogenetic relationships among magnoliids and to examine patterns of variation of GC content.

Results

The Drimys, Liriodendron, and Piper plastid genomes are very similar in size at 160,604, 159,886 bp, and 160,624 bp, respectively. Gene content and order are nearly identical to many other unrearranged angiosperm plastid genomes, including Calycanthus, the other published magnoliid genome. Overall GC content ranges from 34–39%, and coding regions have a substantially higher GC content than non-coding regions. Among protein-coding genes, GC content varies by codon position with 1st codon > 2nd codon > 3rd codon, and it varies by functional group with photosynthetic genes having the highest percentage and NADH genes the lowest. Phylogenetic analyses using parsimony and likelihood methods and sequences of 61 protein-coding genes provided strong support for the monophyly of magnoliids and two strongly supported groups were identified, the Canellales/Piperales and the Laurales/Magnoliales. Strong support is reported for monocots and eudicots as sister clades with magnoliids diverging before the monocot-eudicot split. The trees also provided moderate or strong support for the position of Amborella as sister to a clade including all other angiosperms.

Conclusion

Evolutionary comparisons of three new magnoliid plastid genome sequences, combined with other published angiosperm genomes, confirm that GC content is unevenly distributed across the genome by location, codon position, and functional group. Furthermore, phylogenetic analyses provide the strongest support so far for the hypothesis that the magnoliids are sister to a large clade that includes both monocots and eudicots.

Viburnum phylogeny: evidence from the duplicated nuclear gene GBSSI

DNA sequencing studies of the granule-bound starch synthase gene (GBSSI) indicate the presence of two loci in Viburnum. Gene trees from separate and combined phylogenetic analyses of the GBSSI paralogues are generally congruent with each other and with trees from previous analyses, especially those of Donoghue et al. [Syst. Bot. 29 (2004) 188] based on nuclear ribosomal ITS and chloroplast trnK intron DNA sequences. Specifically, our GBSSI trees confirm (i) the monophyly of some and non-monophyly of other traditionally recognized taxonomic sections, (ii) the presence of three major supra-sectional lineages within Viburnum, and (iii) the resolution of many species relationships within the section-level clades. Analyses of GBSSI also provide greater resolution of relationships within the largest supra-sectional lineage. Relationships at the base of the Viburnum phylogeny remain uncertain; in particular, the position of the root, relationships among the supra-sectional clades, and the exact placement of several smaller groups (e.g., Viburnum clemensiae, Viburnum urceolatum, and section Pseudotinus). In two lineages each GBSSI paralogue is represented by two distinct sequences. The presence of additional copies appears to be correlated with polyploidy in these clades. Placement of the homoeologues in our gene trees suggests the possibility of a hybrid origin for these polyploids.

The loss of the hemoglobin H2S-binding function in annelids from sulfide-free habitats reveals molecular adaptation driven by Darwinian positive selection

The hemoglobin of the deep-sea hydrothermal vent vestimentiferan Riftia pachyptila (annelid) is able to bind toxic hydrogen sulfide (H(2)S) to free cysteine residues and to transport it to fuel endosymbiotic sulfide-oxidising bacteria. The cysteine residues are conserved key amino acids in annelid globins living in sulfide-rich environments, but are absent in annelid globins from sulfide-free environments. Synonymous and nonsynonymous substitution analysis from two different sets of orthologous annelid globin genes from sulfide rich and sulfide free environments have been performed to understand how the sulfide-binding function of hemoglobin appeared and has been maintained during the course of evolution. This study reveals that the sites occupied by free-cysteine residues in annelids living in sulfide-rich environments and occupied by other amino acids in annelids from sulfide-free environments, have undergone positive selection in annelids from sulfide-free environments. We assumed that the high reactivity of cysteine residues became a disadvantage when H(2)S disappeared because free cysteines without their natural ligand had the capacity to interact with other blood components, disturb homeostasis, reduce fitness and thus could have been counterselected. To our knowledge, we pointed out for the first time a case of function loss driven by molecular adaptation rather than genetic drift. If constraint relaxation (H(2)S disappearance) led to the loss of the sulfide-binding function in modern annelids from sulfide-free environments, our work suggests that adaptation to sulfide-rich environments is a plesiomorphic feature, and thus that the annelid ancestor could have emerged in a sulfide-rich environment.

Whole-genome comparison of leucine-rich repeat extensins in Arabidopsis and rice. A conserved family of cell wall proteins form a vegetative and a reproductive clade

We have searched the Arabidopsis and rice (Oryza sativa) genomes for homologs of LRX1, an Arabidopsis gene encoding a novel type of cell wall protein containing a leucine-rich repeat (LRR) and an extensin domain. Eleven and eight LRX (LRR/EXTENSIN) genes have been identified in these two plant species, respectively. The LRX gene family encodes proteins characterized by a short N-terminal domain, a domain with 10 LRRs, a cysteine-rich motif, and a variable C-terminal extensin-like domain. Phylogenetic analysis performed on the conserved domains indicates the existence of two major clades of LRX proteins that arose before the eudicot/monocot divergence and then diversified independently in each lineage. In Arabidopsis, gene expression studies by northern hybridization and promoter::uidA fusions showed that the two phylogenetic clades represent a specialization into "reproductive" and "vegetative" LRXs. The four Arabidopsis genes of the "reproductive" clade are specifically expressed in pollen, whereas the seven "vegetative" genes are predominantly expressed in various sporophytic tissues. This separation into two expression classes is also supported by previous studies on maize (Zea mays) and tomato (Lycopersicon esculentum) LRX homologs and by information on available rice ESTs. The strong conservation of the amino acids responsible for the putative recognition specificity of the LRR domain throughout the family suggests that the LRX proteins interact with similar ligands.

Evolution of the sulfide-binding function within the globin multigenic family of the deep-sea hydrothermal vent tubeworm Riftia pachyptila

The giant extracellular hexagonal bilayer hemoglobin (HBL-Hb) of the deep-sea hydrothermal vent tube worm Riftia pachyptila is able to transport simultaneously O(2) and H(2)S in the blood from the gills to a specific organ: the trophosome that harbors sulfide-oxidizing endosymbionts. This vascular HBL-Hb is made of 144 globins from which four globin types (A1, A2, B1, and B2) coevolve. The H(2)S is bound at a specific location (not on the heme site) onto two of these globin types. In order to understand how such a function emerged and evolved in vestimentiferans and other related annelids, six partial cDNAs corresponding to the six globins known to compose the multigenic family of R. pachyptila have been identified and sequenced. These partial sequences (ca. 120 amino acids, i.e., 80% of the entire protein) were used to reconstruct molecular phylogenies in order to trace duplication events that have led to the family organization of these globins and to locate the position of the free cysteine residues known to bind H(2)S. From these sequences, only two free cysteine residues have been found to occur, at positions Cys + 1 (i.e., 1 a.a. from the well-conserved distal histidine) and Cys + 11 (i.e., 11 a.a. from the same histidine) in globins B2 and A2, respectively. These two positions are well conserved in annelids, vestimentiferans, and pogonophorans, which live in sulfidic environments. The structural comparison of the hydrophobic environment that surrounds these cysteine residues (the sulfide-binding domain) using hydrophobic cluster analysis plots, together with the cysteine positions in paralogous strains, suggests that the sulfide-binding function might have emerged before the annelid radiation in order to detoxify this toxic compound. Moreover, globin evolutionary rates are highly different between paralogous strains. This suggests that either the two globin subfamilies involved in the sulfide-binding function (A2 and B2) have evolved under strong directional selective constraints (negative selection) and that the two other globins (A1 and B1) have accumulated more substitutions through positive selection or have evolved neutrally after a relaxation of selection pressures. A likely scenario on the evolution of this multigenic family is proposed and discussed from this data set.

Uninode coding vs gene tree parsimony for phylogenetic reconstruction using duplicate genes

Two different methods of using paralogous genes for phylogenetic inference have been proposed: reconciled trees (or gene tree parsimony) and uninode coding. Gene tree parsimony suffers from 10 serious problems, including differential weighting of nucleotide and gap characters, undersampling which can be misinterpreted as synapomorphy, all of the characters not being allowed to interact, and conflict between gene trees being given equal weight, regardless of branch support. These problems are largely avoided by using uninode coding. The uninode coding method is elaborated to address multiple gene duplications within a single gene tree family and handle problems caused by lack of gene tree resolution. An example of vertebrate phylogeny inferred from nine genes is reanalyzed using uninode coding. We suggest that uninode coding be used instead of gene tree parsimony for phylogenetic inference from paralogous genes.

Was the ANITA rooting of the angiosperm phylogeny affected by long-branch attraction? Amborella, Nymphaeales, Illiciales, Trimeniaceae, and Austrobaileya

Five groups of basal angiosperms, Amborella, Nymphaeales, Illiciales, Trimeniaceae, and Austrobaileya (ANITA), were identified in several recent studies as representing a series of the earliest-diverging lineages of the angiosperm phylogeny. All of these studies except one employed a multigene analysis approach and used gymnosperms as the outgroup to determine the ingroup topology. The high level of divergence between gymnosperms and angiosperms, however, has long been implicated in the difficulty of reconstructing relationships at the base of angiosperm phylogeny using DNA sequences, for fear of long-branch attraction (LBA). In this study, we replaced the gymnosperm sequences from the five-gene matrix (mitochondrial atp1 and matR, plastid atpB and rbcL, and nuclear 18S rDNA) used in our earlier study with four categories of divergent sequences--random sequences with equal base frequencies or equally AT- and GC-rich contents, homopolymers and heteropolymers, misaligned gymnosperm sequences, and aligned lycopod and bryophyte sequences--to evaluate whether the gymnosperms were an appropriate outgroup to angiosperms in our earlier study that identified the ANITA rooting. All 24 analyses performed rooted the angiosperm phylogeny at either Acorus or Alisma (or Alisma-Triglochin-Potamogeton in one case due to use of a slightly different alignment) and placed the monocots as a basal grade, producing genuine LBA results. These analyses demonstrate that the identification of ANITA as the basalmost extant angiosperms was based on historical signals preserved in the gymnosperm sequences and that the gymnosperms were an appropriate outgroup with which to root the angiosperm phylogeny in the multigene sequence analysis. This strategy of evaluating the appropriateness of an outgroup using artificial sequences and a series of outgroups with increments of divergence levels can be applied to investigations of phylogenetic patterns at the bases of other major clades, such as land plants, animals, and eukaryotes.

Rooting a phylogeny with homologous genes on opposite sex chromosomes (gametologs): a case study using avian CHD

We describe a previously unrecognized form of gene homology using the term "gametology," which we define as homology arising through lack of recombination and subsequent differentiation of sex chromosomes. We demonstrate use of gametologous genes to root each other in phylogenetic analyses of sex-specific avian Chromo-helicase-DNA binding gene (CHD) sequences. Phylogenetic analyses of a set of neognath bird sequences yield monophyletic groups for CHD-W and CHD-Z gametologs, as well as congruent relationships between these two clades and between them and current views of avian taxonomy. Phylogenetic analyses including paleognath bird CHD sequences and rooting with crocodilian CHD sequences, suggest an early divergence for paleognath CHD within the avian CHD clade. Based on our CHD analyses calibrated with avian fossil dates, we estimate the divergence between CHD-W and CHD-Z at 123 MYA, suggesting an early differentiation of sex chromosomes that predates most extant avian orders. In agreement with the notion of male-driven evolution, we find a faster rate of change in male-linked CHD-Z sequences.

The granule-bound starch synthase (GBSSI) gene in the Rosaceae: multiple loci and phylogenetic utility

We sampled the 5' end of the granule-bound starch synthase gene (GBSSI or waxy) in Rosaceae, sequencing 108 clones from 18 species in 14 genera representing all four subfamilies (Amygdaloideae, Maloideae, Rosoideae, and Spiraeoideae), as well as four clones from Rhamnus catharticus (Rhamnaceae). This is the first phylogenetic study to use the 5' portion of this nuclear gene. Parsimony and maximum-likelihood analyses of 941 bases from seven complete and two partial exons demonstrate the presence of two loci (GBSSI-1 and GBSSI-2) in the Rosaceae. Southern hybridization analyses with locus-specific probes confirm that all four Rosaceae subfamilies have at least two GBSSI loci, even though only one locus has been reported in all previously studied diploid flowering plants. Phylogenetic analyses also identify four clades representing four loci in the Maloideae. Phylogenetic relationships inferred from GBSSI sequences are largely compatible with those from chloroplast (cpDNA: ndhF, rbcL) and nuclear ribosomal internal transcribed spacer (nrITS) DNA. Large clades are marked by significant intron variation: a long first intron plus no sixth intron in Maloideae GBSSI-1, a long fourth intron in Rosoideae GBSSI-1, and a GT to GC mutation in the 5' splice site of the fourth intron in all GBSSI-2 sequences. Our data do not support the long-held hypothesis that Maloideae originated from an ancient hybridization between amygdaloid and spiraeoid ancestors. Instead, Spiraeoideae genera (Kageneckia and Vauquelinia) are their closest relatives in all four GBSSI clades.

Phytochromes and light signal perception by plants--an emerging synthesis

For plants, the sensing of light in the environment is as important as vision is for animals. Fluctuations in light can be crucial to competition and survival. One way plants sense light is through the phytochromes, a small family of diverse photochromic protein photoreceptors whose origins have been traced to the photosynthetic prokaryotes. During their evolution, the phytochromes have acquired sophisticated mechanisms to monitor light. Recent advances in understanding the molecular mechanisms of phytochromes and their significance to evolutionary biology make possible an interim synthesis of this rapidly advancing branch of photobiology.

Phylogeny reconstruction using duplicate genes

In this paper, we propose a new method (uninode coding) for coding duplicate (paralogous) genes to infer species trees. Uninode coding incorporates data from duplicated and unduplicated gene copies in phylogenetic analyses of taxa. Uninode coding utilizes global parsimony through the inclusion of both duplicated and unduplicated gene copies, allows one to code all data sources from a taxon into a single terminal, and overcomes problems of character dependence among duplicated and unduplicated gene copies. We present an example of uninode coding using the phytochrome A and phytochrome C data from a study by Donoghue and Mathews.

Multigene analyses identify the three earliest lineages of extant flowering plants

Flowering plants (angiosperms) are by far the largest, most diverse, and most important group of land plants, with over 250,000 species and a dominating presence in most terrestrial ecosystems. Understanding the origin and early diversification of angiosperms has posed a long-standing botanical challenge [1]. Numerous morphological and molecular systematic studies have attempted to reconstruct the early history of this group, including identifying the root of the angiosperm tree. There is considerable disagreement among these studies, with various groups of putatively basal angiosperms from the subclass Magnoliidae having been placed at the root of the angiosperm tree (reviewed in [2-4]). We investigated the early evolution of angiosperms by conducting combined phylogenetic analyses of five genes that represent all three plant genomes from a broad sampling of angiosperms. Amborella, a monotypic, vessel-less dioecious shrub from New Caledonia, was clearly identified as the first branch of angiosperm evolution, followed by the Nymphaeales (water lillies), and then a clade of woody vines comprising Schisandraceae and Austrobaileyaceae. These findings are remarkably congruent with those from several concurrent molecular studies [5-7] and have important implications for whether or not the first angiosperms were woody and contained vessels, for interpreting the evolution of other key characteristics of basal angiosperms, and for understanding the timing and pattern of angiosperm origin and diversification.

The root of angiosperm phylogeny inferred from duplicate phytochrome genes

An analysis of duplicate phytochrome genes (PHYA and PHYC) is used to root the angiosperms, thereby avoiding the inclusion of highly diverged outgroup sequences. The results unambiguously place the root near Amborella (one species, New Caledonia) and resolve water lilies (Nymphaeales, approximately 70 species, cosmopolitan), followed by Austrobaileya (one species, Australia), as early branches. These findings bear directly on the interpretation of morphological evolution and diversification within angiosperms.

How good are deep phylogenetic trees?

Great interest is given to species emerging early in phylogenetic reconstruction because they are often assumed to represent an ancestor. Recent studies indicate, however, that species branching deep in molecular trees are often fast-evolving ones, misplaced because of the long-branch artefact. The detection of genuinely deep-branching organisms remains an elusive task.