Phylogenetic analysis of the plant-specific zinc finger-homeobox and mini zinc finger gene families

Zinc finger-homeodomain proteins (ZHD) are present in many plants; however, the evolutionary history of the ZHD gene family remains largely unknown. We show here that ZHD genes are plant-specific, nearly all intronless, and related to MINI ZINC FINGER ( MIF ) genes that possess only the zinc finger. Phylogenetic analyses of ZHD genes from representative land plants suggest that non-seed plant ZHD genes occupy basal positions and angiosperm homologs form seven distinct clades. Several clades contain genes from two or more major angiosperm groups, including eudicots, monocots, magnoliids, and other basal angiosperms, indicating that several duplications occurred before the diversification of flowering plants. In addition, specific lineages have experienced more recent duplications. Unlike the ZHD genes, MIF s are found only from seed plants, possibly derived from ZHD s by loss of the homeodomain before the divergence of seed plants. Moreover, the MIF genes have also undergone relatively recent gene duplications. Finally, genome duplication might have contributed substantially to the expansion of family size in angiosperms and caused a high level of functional redundancy/overlap in these genes.     

Molecular evolution of plant haemoglobin: two haemoglobin genes in Nymphaeaceae Euryale ferox

We isolated and sequenced two haemoglobin genes from the early-branching angiosperm Euryale ferox (Nymphaeaceae). The two genes belong to the two known classes of plant haemoglobin. Their existence in Nymphaeaceae supports the theory that class 1 haemoglobin was ancestrally present in all angiosperms, and is evidence for class 2 haemoglobin being widely distributed. These sequences allowed us to unambiguously root the angiosperm haemoglobin phylogeny, and to corroborate the hypothesis that the class 1/class 2 duplication event occurred before the divergence between monocots and eudicots. We addressed the molecular evolution of plant haemoglobin by comparing the synonymous and nonsynonymous substitution rates in various groups of genes. Class 2 haemoglobin genes of legumes (functionally involved in a symbiosis with nitrogen-fixing bacteria) show a higher nonsynonymous substitution rate than class 1 (nonsymbiotic) haemoglobin genes. This suggests that a change in the selective forces applying to plant haemoglobins has occurred during the evolutionary history of this gene family, potentially in relation with the evolution of symbiosis.     

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.     

Patterns of gene duplication and functional evolution during the diversification of the AGAMOUS subfamily of MADS box genes in angiosperms

Members of the AGAMOUS (AG) subfamily of MIKC-type MADS-box genes appear to control the development of reproductive organs in both gymnosperms and angiosperms. To understand the evolution of this subfamily in the flowering plants, we have identified 26 new AG-like genes from 15 diverse angiosperm species. Phylogenetic analyses of these genes within a large data set of AG-like sequences show that ancient gene duplications were critical in shaping the evolution of the subfamily. Before the radiation of extant angiosperms, one event produced the ovule-specific D lineage and the well-characterized C lineage, whose members typically promote stamen and carpel identity as well as floral meristem determinacy. Subsequent duplications in the C lineage resulted in independent instances of paralog subfunctionalization and maintained functional redundancy. Most notably, the functional homologs AG from Arabidopsis and PLENA (PLE) from Antirrhinum are shown to be representatives of separate paralogous lineages rather than simple genetic orthologs. The multiple subfunctionalization events that have occurred in this subfamily highlight the potential for gene duplication to lead to dissociation among genetic modules, thereby allowing an increase in morphological diversity.     

Portrait of the expansin superfamily in Physcomitrella patens: comparisons with angiosperm expansins

BACKGROUND AND AIMS: Expansins are plant cell wall loosening proteins important in a variety of physiological processes. They comprise a large superfamily of genes consisting of four families (EXPA, EXPB, EXLA and EXLB) whose evolutionary relationships have been well characterized in angiosperms, but not in basal land plants. This work attempts to connect the expansin superfamily in bryophytes with the evolutionary history of this superfamily in angiosperms. METHODS: The expansin superfamily in Physcomitrella patens has been assembled from the Physcomitrella sequencing project data generated by the Joint Genome Institute and compared with angiosperm expansin superfamilies. Phylogenetic, motif, intron and distance analyses have been used for this purpose. KEY RESULTS: A gene superfamily is revealed that contains similar numbers of genes as found in arabidopsis, but lacking EXLA or EXLB genes. This similarity in gene numbers exists even though expansin evolution in Physcomitrella diverged from the angiosperm line approx. 400 million years ago. Phylogenetic analyses suggest that there were a minimum of two EXPA genes and one EXPB gene in the last common ancestor of angiosperms and Physcomitrella. Motif analysis seems to suggest that EXPA protein function is similar in bryophytes and angiosperms, but that EXPB function may be altered. CONCLUSIONS: The EXPA genes of Physcomitrella are likely to have maintained the same biochemical function as angiosperm expansins despite their independent evolutionary history. Changes seen at normally conserved residues in the Physcomitrella EXPB family suggest a possible change in function as one mode of evolution in this family.     

Unraveling the Distribution and Evolution of miR156 targeted SPLs in Plants by Phylogenetic Analysis

Squamosa promoter binding protein like genes (SPLs) are critical during plant development and mostly regulated by miR156. However, little is known about phylogenetic distribution and evolutionary patterns of miR156 targeted SPLs. In this study, 183 SPLs from nine genome sequenced species representing algae, bryophytes, lycophyte, monocots, and eudicots were computationally analyzed. Our results showed that miR156 responsive elements (MREs) on SPLs were present in land plants but absent from unicellular green algae. Phylogenetic analysis revealed that miR156 targeted SPLs only distributed in group II not group I of land plants, suggesting they originated from a common ancestor. In addition, group II were further divided into seven subgroups (IIa IIg) and miR156 targeted SPLs distributed in some specific members of SPLs from six subgroups except subgroup IId. Such distribution pattern was well elucidated by gene structure evolution of miR156 targeted SPLs based on the correlation of phylogenetic classification and gene structure. They could suffer from the exon loss events combined with MREs loss during evolution. Moreover, gene duplication contributed to the abundance of miR156 targeted SPLs, which had significantly increased after angiosperms and lower plants split. With Arabidopsis as the model species, we found segmental and tandem gene duplications predominated during miR156 targeted SPLs expansion. Taken together, these results provide better insights in understanding the function diversity and evolution of miR156 targeted SPLs in plants.     

Angiosperm phylogeny inferred from sequences of four mitochondrial genes

An angiosperm phylogeny was reconstructed in a maximum likelihood analysis of sequences of four mitochondrial genes, atpl, matR, had5, and rps3, from 380 species that represent 376 genera and 296 families of seed plants. It is largely congruent with the phylogeny of angiosperms reconstructed from chloroplast genes atpB, matK, and rbcL, and nuclear 18S rDNA. The basalmost lineage consists of Amborella and Nymphaeales (including Hydatellaceae). Austrobaileyales follow this clade and are sister to the mesangiosperms, which include Chloranthaceae, Ceratophyllum, magnoliids, monocots, and eudicots. With the exception of Chloranthaceae being sister to Ceratophyllum, relationships among these five lineages are not well supported. In eudicots, Ranunculales, Sabiales, Proteales, Trochodendrales, Buxales, Gunnerales, Saxifragales, Vitales, Berberidopsidales, and Dilleniales form a basal grade of lines that diverged before the diversification of rosids and asterids. Within rosids, the COM (Celastrales-Oxalidales-Malpighiales) clade is sister to malvids (or rosid Ⅱ), instead of to the nitrogen-fixing clade as found in all previous large-scale molecular analyses of angiosperms. Santalales and Caryophyllales are members of an expanded asterid clade. This study shows that the mitochondrial genes are informative markers for resolving relationships among genera, families, or higher rank taxa across angiosperms. The low substitution rates and low homoplasy levels of the mitochondrial genes relative to the chloroplast genes, as found in this study, make them particularly useful for reconstructing ancient phylogenetic relationships. A mitochondrial gene-based angiosperm phylogeny provides an independent and essential reference for comparison with hypotheses of angiosperm phylogeny based on chloroplast genes, nuclear genes, and non-molecular data to reconstruct the underlying organismal phylogeny.     

Evolution of Oleosin in Land Plants

Oleosins form a steric barrier surface on lipid droplets in cytoplasm, preventing them from contacting and coalescing with adjacent droplets. Oleosin genes have been detected in numerous plant species. However, the presence of oleosin genes in the most basally diverging lineage of land plants, liverworts, has not been reported previously. Thus we explored whether liverworts have an oleosin gene. In Marchantia polymorpha L., a thalloid liverwort, one predicted sequence was found that could encode oleosin, possessing the hallmark of oleosin, a proline knot (-PX₅SPX₃P-) motif. The phylogeny of the oleosin gene family in land plants was reconstructed based on both nucleotide and amino acid sequences of oleosins, from 31 representative species covering almost all the main lineages of land plants. Based on our phylogenetic trees, oleosin genes were classified into three groups: M-oleosins (defined here as a novel group distinct from the two previously known groups), low molecular weight isoform (L-oleosin), and high molecular weight isoform (H-oleosin), according to their amino-acid organization, phylogenetic relationships, expression tissues, and immunological characteristics. In liverworts, mosses, lycophytes, and gymnosperms, only M-oleosins have been described. In angiosperms, however, while this isoform remains and is highly expressed in the gametophyte pollen tube, two other isoforms also occur, L-oleosins and H-oleosins. Phylogenetic analyses suggest that the M-oleosin isoform is the precursor to the ancestor of L-oleosins and H-oleosins. The later two isoforms evolved by successive gene duplications in ancestral angiosperms. At the genomic level, most oleosins possess no introns. If introns are present, in both the L-isoform and the M-isoform a single intron inserts behind the central region, while in the H-isoform, a single intron is located at the 5′-terminus. This study fills a major gap in understanding functional gene evolution of oleosin in land plants, shedding new light on evolutionary transitions of lipid storage strategies.     

Expression of floral MADS-box genes in basal angiosperms: implications for the evolution of floral regulators

The ABC model of floral organ identity is based on studies of Arabidopsis and Antirrhinum, both of which are highly derived eudicots. Most of the genes required for the ABC functions in Arabidopsis and Antirrhinum are members of the MADS-box gene family, and their orthologs are present in all major angiosperm lineages. Although the eudicots comprise 75% of all angiosperms, most of the diversity in arrangement and number of floral parts is actually found among basal angiosperm lineages, for which little is known about the genes that control floral development. To investigate the conservation and divergence of expression patterns of floral MADS-box genes in basal angiosperms relative to eudicot model systems, we isolated several floral MADS-box genes and examined their expression patterns in representative species, including Amborella (Amborellaceae), Nuphar (Nymphaeaceae) and Illicium (Austrobaileyales), the successive sister groups to all other extant angiosperms, plus Magnolia and Asimina, members of the large magnoliid clade. Our results from multiple methods (relative-quantitative RT-PCR, real-time PCR and RNA in situ hybridization) revealed that expression patterns of floral MADS-box genes in basal angiosperms are broader than those of their counterparts in eudicots and monocots. In particular, (i) AP1 homologs are generally expressed in all floral organs and leaves, (ii) AP3/PI homologs are generally expressed in all floral organs and (iii) AG homologs are expressed in stamens and carpels of most basal angiosperms, in agreement with the expectations of the ABC model; however, an AG homolog is also expressed in the tepals of Illicium. The broader range of strong expression of AP3/PI homologs is inferred to be the ancestral pattern for all angiosperms and is also consistent with the gradual morphological intergradations often observed between adjacent floral organs in basal angiosperms.     

Functional analyses of AGAMOUS family members in Nicotiana benthamiana clarify the evolution of early and late roles of C-function genes in eudicots

The C-function, according to the ABC model of floral organ identity, is required for stamen and carpel development and to provide floral meristem determinacy. Members of the AG lineage of the large MADS box gene family specify the C-function in a broadly conserved manner in angiosperms. In core eudicots, two sub-lineages co-exist, euAG and PLE, which have been extensively characterized in Antirrhinum majus and Arabidopsis thaliana, where strong sub-functionalization has led to highly divergent contributions of the respective paralogs to the C-function. Various scenarios have been proposed to reconstruct the evolutionary history of the euAG and PLE lineages in eudicots, but detailed functional analyses of the roles of these genes in additional representative species to validate evolutionary hypotheses are scarce. Here, we report functional characterization of euAG- and PLE-like genes in Nicotiana benthamiana through expression analyses and phenotypic characterization of the defects caused by their specific down-regulation. We show that both paralogs redundantly contribute to the C-function in this species, providing insights on the likely evolution of these gene lineages following divergence of the major groups within the eudicots (rosids and asterids). Moreover, we have demonstrated a conserved role for the PLE-like genes in controlling fruit dehiscence, which strongly supports the ancestral role of PLE-like genes in late fruit development and suggests a common evolutionary origin of late developmental processes in dry (dehiscent) and fleshy (ripening) fruits.     

Angiosperm divergence times: the effect of genes, codon positions, and time constraints

An understanding of the evolution of modern terrestrial ecosystems requires an understanding of the dynamics associated with angiosperm evolution, including the timing of their origin and diversification into their extraordinary present-day diversity. Molecular estimates of angiosperm age have varied widely, and many substantially predate the Early Cretaceous fossil appearance of the group. In this study, the effect of different genes, codon positions, and chronological constraints on node ages are examined on divergence time estimates across seed plants, with a special focus on angiosperms. Penalized likelihood was used to estimate divergence times on a phylogenetic hypothesis for seed plants derived from Bayesian analysis, with branch lengths estimated with maximum likelihood. The plastid genes atpB, psaA, psbB, and rbcL were used individually and in combination, using first and second, third, and the three codon positions, including and excluding age constraints on 20 nodes derived from a critical examination of the land-plant fossil record. The optimal level of rate smoothing according to each unconstrained and constrained dataset was obtained with penalized likelihood. Tests for a molecular clock revealed significantly unclocklike rates in all datasets. Addition of fossil constraints resulted in even greater departures from constancy. Consistently with significant deviations from a clock, estimated optimal smoothing values were low, but a strict correlation between rate heterogeneity and optimal smoothing value was not found. Age estimates for nodes across the phylogeny varied, sometimes substantially, with gene and codon position. Nevertheless, estimates based on the four concatenated genes are very similar to the mean of the four individual gene estimates. For any given node, unconstrained age estimates are more variable than constrained estimates and are frequently younger than well-substantiated fossil members of the clade. Constrained estimates of ages of clades are older than unconstrained estimates and oldest fossil representatives, sometimes substantially so. Angiosperm age estimates decreased as rate smoothing increased. Whereas the range of unconstrained angiosperm age estimates spans the fossil age of the clade, the range of constrained estimates is narrower (and older) than the earliest angiosperm fossils. Results unambiguously indicate the relevance of constraints in reducing the variability of ages derived from different partitions of the data and diminishing the effect of the smoothing parameter. Constrained optimizations of divergence times and substitution rates across the phylogeny suggest appreciably different evolutionary dynamics for angiosperms and for gymnosperms. Whereas the gymnosperm crown group originated shortly after the origin of seed plants, a long time elapsed before the origin of crown group angiosperms. Although absolute age estimates of angiosperms and angiosperm clades are older than their earliest fossils, the estimated pace of phylogenetic diversification largely agrees with the rapid appearance of angiosperm lineages in stratigraphic sequences.     

The origin and diversification of angiosperms

The angiosperms, one of five groups of extant seed plants, are the largest group of land plants. Despite their relatively recent origin, this clade is extremely diverse morphologically and ecologically. However, angiosperms are clearly united by several synapomorphies. During the past 10 years, higher-level relationships of the angiosperms have been resolved. For example, most analyses are consistent in identifying Amborella, Nymphaeaceae, and Austrobaileyales as the basalmost branches of the angiosperm tree. Other basal lineages include Chloranthaceae, magnoliids, and monocots. Approximately three quarters of all angiosperm species belong to the eudicot clade, which is strongly supported by molecular data but united morphologically by a single synapomorphy-triaperturate pollen. Major clades of eudicots include Ranunculales, which are sister to all other eudicots, and a clade of core eudicots, the largest members of which are Saxifragales, Caryophyllales, rosids, and asterids. Despite rapid progress in resolving angiosperm relationships, several significant problems remain: (1) relationships among the monocots, Chloranthaceae, magnoliids, and eudicots, (2) branching order among basal eudicots, (3) relationships among the major clades of core eudicots, (4) relationships within rosids, (5) relationships of the many lineages of parasitic plants, and (6) integration of fossils with extant taxa into a comprehensive tree of angiosperm phylogeny.     

Functional diversification of B MADS-box homeotic regulators of flower development: Adaptive evolution in protein-protein interaction domains after major gene duplication events

B-class MADS-box genes have been shown to be the key regulators of petal and stamen specification in several eudicot model species such as Arabidopsis thaliana, Antirrhinum majus, and Petunia hybrida. Orthologs of these genes have been found across angiosperms and gymnosperms, and it is thought that the basic regulatory function of B proteins is conserved in seed plant lineages. The evolution of B genes is characterized by numerous duplications that might represent key elements fostering the functional diversification of duplicates with a deep impact on their role in the evolution of the floral developmental program. To evaluate this, we performed a rigorous statistical analysis with B gene sequences. Using maximum likelihood and Bayesian methods, we estimated molecular substitution rates and determined the selective regimes operating at each residue of B proteins. We implemented tests that rely on phylogenetic hypotheses and codon substitution models to detect significant differences in substitution rates (DSRs) and sites under positive adaptive selection (PS) in specific lineages before and after duplication events. With these methods, we identified several protein residues fixed by PS shortly after the origin of PISTILLATA-like and APETALA3-like lineages in angiosperms and shortly after the origin of the euAP3-like lineage in core eudicots, the 2 main B gene duplications. The residues inferred to have been fixed by positive selection lie mostly within the K domain of the protein, which is key to promote heterodimerization. Additionally, we used a likelihood method that accommodates DSRs among lineages to estimate duplication dates for AP3-PI and euAP3-TM6, calibrating with data from the fossil record. The dates obtained are consistent with angiosperm origins and diversification of core eudicots. Our results strongly suggest that novel multimer formation with other MADS proteins could have been crucial for the functional divergence of B MADS-box genes. We thus propose a mechanism of functional diversification and persistence of gene duplicates by the appearance of novel multimerization capabilities after duplications. Multimer formation in different combinations of regulatory proteins can be a mechanistic basis for the origin of novel regulatory functions and a gene regulatory mechanism for the appearance of morphological innovations.     

Geographic patterns and climatic correlates of deep evolutionary legacies for angiosperm assemblages in China

Deep evolutionary histories can play an important role in assembling species into communities, but few studies have explored the effects of deep evolutionary histories on species assembly of angiosperms (flowering plants). Here we explore patterns of family divergence and diversification times (stem and crown ages, respectively) and phylogenetic fuses for angiosperm assemblages in 100 × 100 km grid cells across geographic and ecological gradients in China. We found that both family stem and crown ages of angiosperm assemblages are older in southeastern China with warm and humid climates than in northwestern China with cold and dry climates; these patterns are stronger for family stem ages than for family crown ages; families in colder and drier climates are more closely related across the family-level angiosperm phylogeny; and family phylogenetic fuses are, on average, longer for angiosperm assemblages in warm and humid climates than in cold and dry climate. We conclude that the fact that deep evolutionary histories, which were measured as family stem and crown ages and family phylogenetic fuses in this study, have shown strong geographic and ecological patterns suggests that deep evolutionary histories of angiosperms have profound effects on assembling angiosperm species into ecological communities.