From algae to angiosperms-inferring the phylogeny of green plants (Viridiplantae) from 360 plastid genomes

Floral biology and ovule and seed ontogeny of Nymphaea thermarum, a water lily at the brink of extinction with potential as a model system for basal angiosperms

BACKGROUND AND AIMS

Nymphaea thermarum is a member of the Nymphaeales, of one of the most ancient lineages of flowering plants. This species was only recently described and then declared extinct in the wild, so little is known about its reproductive biology. In general, the complete ontogeny of ovules and seeds is not well documented among species of Nymphaea and has never been studied in the subgenus Brachyceras, the clade to which N. thermarum belongs.

METHODS

Flowers and fruits were processed for brightfield, epifluorescence and confocal microscopy. Flower morphology, with emphasis on the timing of male and female functions, was correlated with key developmental stages of the ovule and the female gametophyte. Development of the seed tissues and dynamics of polysaccharide reserves in the endosperm, perisperm and embryo were examined.

KEY RESULTS

Pollen release in N. thermarum starts before the flower opens. Cell walls of the micropylar nucellus show layering of callose and cellulose in a manner reminiscent of transfer cell wall patterning. Endosperm development is ab initio cellular, with micropylar and chalazal domains that embark on distinct developmental trajectories. The surrounding maternal perisperm occupies the majority of seed volume and accumulates starch centrifugally. In mature seeds, a minute but fully developed embryo is surrounded by a single, persistent layer of endosperm.

CONCLUSIONS

Early male and female function indicate that N. thermarum is predisposed towards self-pollination, a phenomenon that is likely to have evolved multiple times within Nymphaea. While formation of distinct micropylar and chalazal developmental domains in the endosperm, along with a copious perisperm, characterize the seeds of most members of the Nymphaeales, seed ontogenies vary between and among the constituent families. Floral biology, life history traits and small genome size make N. thermarum uniquely promising as an early-diverging angiosperm model system for genetic and molecular studies.

Variable frequency of plastid RNA editing among ferns and repeated loss of uridine-to-cytidine editing from vascular plants

The distinct distribution and abundance of C-to-U and U-to-C RNA editing among land plants suggest that these two processes originated and evolve independently, but the paucity of information from several key lineages limits our understanding of their evolution. To examine the evolutionary diversity of RNA editing among ferns, we sequenced the plastid transcriptomes from two early diverging species, Ophioglossum californicum and Psilotum nudum. Using a relaxed automated approach to minimize false negatives combined with manual inspection to eliminate false positives, we identified 297 C-to-U and three U-to-C edit sites in the O. californicum plastid transcriptome but only 27 C-to-U and no U-to-C edit sites in the P. nudum plastid transcriptome. A broader comparison of editing content with the leptosporangiate fern Adiantum capillus-veneris and the hornwort Anthoceros formosae uncovered large variance in the abundance of plastid editing, indicating that the frequency and type of RNA editing is highly labile in ferns. Edit sites that increase protein conservation among species are more abundant and more efficiently edited than silent and non-conservative sites, suggesting that selection maintains functionally important editing. The absence of U-to-C editing from P. nudum plastid transcripts and other vascular plants demonstrates that U-to-C editing loss is a recurrent phenomenon in vascular plant evolution.

Direction of illumination controls gametophyte orientation in seedless plants and related algae

The environmental influences that determine dorsiventral or axial gametophyte orientation are unknown for most modern seedless plants. To fill this gap, an experimental laboratory system was employed to evaluate the relative effects of light direction and gravity on body orientation of the dorsiventral green alga Coleochaete orbicularis, and gametophytes of liverworts Blasia pusilla and Marchantia polymorpha, early-diverging moss Sphagnum compactum, and fern Ceratopteris richardii, the latter functioning as experimental control. Replicate clonal cultures were experimentally illuminated only from above, only from below, or from multiple directions, with the same near-saturation PAR level for periods brief enough to minimize nutrient limitation effects, and orientation of new growth was evaluated. For all species tested, direction of illumination exerted stronger control over gametophyte body orientation than gravity. When illuminated only from below: 1) axial Sphagnum gametophores that had initially grown into an overlying air space inverted growth by 180°, burrowing into the substrate; 2) new growth of dorsiventral Blasia, Marchantia, and Ceratopteris gametophytes-whose ventral rhizoids initially penetrated agar substrate and dorsal surfaces initially faced overlying airspace-twisted 180° so that ventral surfaces bearing rhizoids faced overlying air space and rhizoids extended into the air; and 3) Coleochaete lost typical dorsiventral organization and diagnostic dorsal hairs. Direction of illumination also exerted stronger control over orientation of liverwort new growth than surface contact did. These results indicate that early land plants likely inherited light-directed gametophyte body orientation from ancestral streptophyte algae and suggest a mechanism for reorientation of gametophyte-dominant land plants after spatial disturbance.

Ferns: the missing link in shoot evolution and development

Shoot development in land plants is a remarkably complex process that gives rise to an extreme diversity of forms. Our current understanding of shoot developmental mechanisms comes almost entirely from studies of angiosperms (flowering plants), the most recently diverged plant lineage. Shoot development in angiosperms is based around a layered multicellular apical meristem that produces lateral organs and/or secondary meristems from populations of founder cells at its periphery. In contrast, non-seed plant shoots develop from either single apical initials or from a small population of morphologically distinct apical cells. Although developmental and molecular information is becoming available for non-flowering plants, such as the model moss Physcomitrella patens, making valid comparisons between highly divergent lineages is extremely challenging. As sister group to the seed plants, the monilophytes (ferns and relatives) represent an excellent phylogenetic midpoint of comparison for unlocking the evolution of shoot developmental mechanisms, and recent technical advances have finally made transgenic analysis possible in the emerging model fern Ceratopteris richardii. This review compares and contrasts our current understanding of shoot development in different land plant lineages with the aim of highlighting the potential role that the fern C. richardii could play in shedding light on the evolution of underlying genetic regulatory mechanisms.

The Phenotypic and Genetic Underpinnings of Flower Size in Polemoniaceae

Corolla length is a labile flower feature and has strong implications for pollinator success. However, the phenotypic and genetic bases of corolla elongation are not well known, largely due to a lack of good candidate genes for potential genetic exploration and functional work. We investigate both the cellular phenotypic differences in corolla length, as well as the genetic control of this trait, in Saltugilia (Polemoniaceae). Taxa in this clade exhibit a large range of flower sizes and differ dramatically in pollinator guilds. Flowers of each species were collected from multiple individuals during four stages of flower development to ascertain if cell number or cell size is more important in determining flower size. In Saltugilia, increased flower size during development appears to be driven more by cell size than cell number. Differences in flower size between species are governed by both cell size and cell number, with the large-flowered S. splendens subsp. grantii having nearly twice as many cells as the small-flowered species. Fully mature flowers of all taxa contain jigsaw cells similar to cells seen in sepals and leaves; however, these cells are not typically found in the developing flowers of most species. The proportion of this cell type in mature flowers appears to have substantial implications, comprising 17-68% of the overall flower size. To identify candidate genes responsible for differences in cell area and cell type, transcriptomes were generated for two individuals of the species with the smallest (S. australis) and largest (S. splendens subsp. grantii) flowers across the same four developmental stages visualized with confocal microscopy. Analyses identified genes associated with cell wall formation that are up-regulated in the mature flower stage compared to mid-stage flowers (75% of mature size). This developmental change is associated with the origin of jigsaw cells in the corolla tube of mature flowers. Further comparisons between mature flowers in the two species revealed 354 transcripts that are up-regulated in the large-flowered S. splendens subsp. grantii compared to the small-flowered S. australis. These results are likely broadly applicable to Polemoniaceae, a clade of nearly 400 species, with extensive variation in floral form and shape.

Complete chloroplast genome of Macadamia integrifolia confirms the position of the Gondwanan early-diverging eudicot family Proteaceae

BACKGROUND

Sequence data from the chloroplast genome have played a central role in elucidating the evolutionary history of flowering plants, Angiospermae. In the past decade, the number of complete chloroplast genomes has burgeoned, leading to well-supported angiosperm phylogenies. However, some relationships, particulary among early-diverging lineages, remain unresolved. The diverse Southern Hemisphere plant family Proteaceae arose on the ancient supercontinent Gondwana early in angiosperm history and is a model group for adaptive radiation in response to changing climatic conditions. Genomic resources for the family are limited, and until now it is one of the few early-diverging 'basal eudicot' lineages not represented in chloroplast phylogenomic analyses.

RESULTS

The chloroplast genome of the Australian nut crop tree Macadamia integrifolia was assembled de novo from Illumina paired-end sequence reads. Three contigs, corresponding to a collapsed inverted repeat, a large and a small single copy region were identified, and used for genome reconstruction. The complete genome is 159,714 bp in length and was assembled at deep coverage (3.29 million reads; ~2000 x). Phylogenetic analyses based on 83-gene and inverted repeat region alignments, the largest sequence-rich datasets to include the basal eudicot family Proteaceae, provide strong support for a Proteales clade that includes Macadamia, Platanus and Nelumbo. Genome structure and content followed the ancestral angiosperm pattern and were highly conserved in the Proteales, whilst size differences were largely explained by the relative contraction of the single copy regions and expansion of the inverted repeats in Macadamia.

CONCLUSIONS

The Macadamia chloroplast genome presented here is the first in the Proteaceae, and confirms the placement of this family with the morphologically divergent Plantanaceae (plane tree family) and Nelumbonaceae (sacred lotus family) in the basal eudicot order Proteales. It provides a high-quality reference genome for future evolutionary studies and will be of benefit for taxon-rich phylogenomic analyses aimed at resolving relationships among early-diverging angiosperms, and more broadly across the plant tree of life.

A precise chloroplast genome of Nelumbo nucifera (Nelumbonaceae) evaluated with Sanger, Illumina MiSeq, and PacBio RS II sequencing platforms: insight into the plastid evolution of basal eudicots

BACKGROUND

The chloroplast genome is important for plant development and plant evolution. Nelumbo nucifera is one member of relict plants surviving from the late Cretaceous. Recently, a new sequencing platform PacBio RS II, known as 'SMRT (Single Molecule, Real-Time) sequencing', has been developed. Using the SMRT sequencing to investigate the chloroplast genome of N. nucifera will help to elucidate the plastid evolution of basal eudicots.

RESULTS

The sizes of the de novo assembled complete chloroplast genome of N. nucifera were 163,307 bp, 163,747 bp and 163,600 bp with average depths of coverage of 7×, 712× and 105× sequenced by Sanger, Illumina MiSeq and PacBio RS II, respectively. The precise chloroplast genome of N. nucifera was obtained from PacBio RS II data proofread by Illumina MiSeq reads, with a quadripartite structure containing a large single copy region (91,846 bp) and a small single copy region (19,626 bp) separated by two inverted repeat regions (26,064 bp). The genome contains 113 different genes, including four distinct rRNAs, 30 distinct tRNAs and 79 distinct peptide-coding genes. A phylogenetic analysis of 133 taxa from 56 orders indicated that Nelumbo with an age of 177 million years is a sister clade to Platanus, which belongs to the basal eudicots. Basal eudicots began to emerge during the early Jurassic with estimated divergence times at 197 million years using MCMCTree. IR expansions/contractions within the basal eudicots seem to have occurred independently.

CONCLUSIONS

Because of long reads and lack of bias in coverage of AT-rich regions, PacBio RS II showed a great promise for highly accurate 'finished' genomes, especially for a de novo assembly of genomes. N. nucifera is one member of basal eudicots, however, evolutionary analyses of IR structural variations of N. nucifera and other basal eudicots suggested that IR expansions/contractions occurred independently in these basal eudicots or were caused by independent insertions and deletions. The precise chloroplast genome of N. nucifera will present new information for structural variation of chloroplast genomes and provide new insight into the evolution of basal eudicots at the primary sequence and structural level.

Complete chloroplast genome sequence of Omani lime (Citrus aurantiifolia) and comparative analysis within the rosids

The genus Citrus contains many economically important fruits that are grown worldwide for their high nutritional and medicinal value. Due to frequent hybridizations among species and cultivars, the exact number of natural species and the taxonomic relationships within this genus are unclear. To compare the differences between the Citrus chloroplast genomes and to develop useful genetic markers, we used a reference-assisted approach to assemble the complete chloroplast genome of Omani lime (C. aurantiifolia). The complete C. aurantiifolia chloroplast genome is 159,893 bp in length; the organization and gene content are similar to most of the rosids lineages characterized to date. Through comparison with the sweet orange (C. sinensis) chloroplast genome, we identified three intergenic regions and 94 simple sequence repeats (SSRs) that are potentially informative markers with resolution for interspecific relationships. These markers can be utilized to better understand the origin of cultivated Citrus. A comparison among 72 species belonging to 10 families of representative rosids lineages also provides new insights into their chloroplast genome evolution.

Phylotranscriptomic analysis of the origin and early diversification of land plants

Reconstructing the origin and evolution of land plants and their algal relatives is a fundamental problem in plant phylogenetics, and is essential for understanding how critical adaptations arose, including the embryo, vascular tissue, seeds, and flowers. Despite advances in molecular systematics, some hypotheses of relationships remain weakly resolved. Inferring deep phylogenies with bouts of rapid diversification can be problematic; however, genome-scale data should significantly increase the number of informative characters for analyses. Recent phylogenomic reconstructions focused on the major divergences of plants have resulted in promising but inconsistent results. One limitation is sparse taxon sampling, likely resulting from the difficulty and cost of data generation. To address this limitation, transcriptome data for 92 streptophyte taxa were generated and analyzed along with 11 published plant genome sequences. Phylogenetic reconstructions were conducted using up to 852 nuclear genes and 1,701,170 aligned sites. Sixty-nine analyses were performed to test the robustness of phylogenetic inferences to permutations of the data matrix or to phylogenetic method, including supermatrix, supertree, and coalescent-based approaches, maximum-likelihood and Bayesian methods, partitioned and unpartitioned analyses, and amino acid versus DNA alignments. Among other results, we find robust support for a sister-group relationship between land plants and one group of streptophyte green algae, the Zygnematophyceae. Strong and robust support for a clade comprising liverworts and mosses is inconsistent with a widely accepted view of early land plant evolution, and suggests that phylogenetic hypotheses used to understand the evolution of fundamental plant traits should be reevaluated.

Chloroplast DNA sequence utility for the lowest phylogenetic and phylogeographic inferences in angiosperms: the tortoise and the hare IV

PREMISE OF THE STUDY

Noncoding chloroplast DNA (NC-cpDNA) sequences are the staple data source of low-level phylogeographic and phylogenetic studies of angiosperms. We followed up on previous papers (tortoise and hare II and III) that sought to identify the most consistently variable regions of NC-cpDNA. We used an exhaustive literature review and newly available whole plastome data to assess applicability of previous conclusions at low taxonomic levels.

METHODS

We aligned complete plastomes of 25 species pairs from across angiosperms, comparing the number of genetic differences found in 107 NC-cpDNA regions and matK. We surveyed Web of Science for the plant phylogeographic literature between 2007 and 2013 to assess how NC-cpDNA has been used at the intraspecific level.

KEY RESULTS

Several regions are consistently the most variable across angiosperm lineages: ndhF-rpl32, rpl32-trnL((UAG)), ndhC-trnV((UAC)), 5'rps16-trnQ((UUG)), psbE-petL, trnT((GGU))-psbD, petA-psbJ, and rpl16 intron. However, there is no universally best region. The average number of regions applied to low-level studies is ∼2.5, which may be too little to access the full discriminating power of this genome.

CONCLUSIONS

Plastome sequences have been used successfully at lower and lower taxonomic levels. Our findings corroborate earlier works, suggesting that there are regions that are most likely to be the most variable. However, while NC-cpDNA sequences are commonly used in plant phylogeographic studies, few of the most variable regions are applied in that context. Furthermore, it appears that in most studies too few NC-cpDNAs are used to access the discriminating power of the cpDNA genome.

Arabinogalactan proteins mark stigmatic receptivity in the protogynous flowers of Magnolia virginiana (Magnoliaceae)

PREMISE OF THE STUDY

Factors affecting floral receptivity in angiosperms remain opaque, but recent studies suggest that the acquisition of stigmatic receptivity associated with cell-wall-related arabinogalactan proteins (AGPs) may be a widespread feature of flowering plants. Here, the time during which a stigma is receptive is evaluated and related to the secretion of AGPs in Magnolia virginiana, a protogynous member of an early-divergent angiosperm clade (magnoliids) with a clearly discernible female receptive phase.

METHODS

Magnolia virginiana flower phenology was documented, and histochemical changes in the stigma before and after pollination were examined. Stigmatic receptivity was evaluated in relation to the secretion of AGPs detected in whole mounts and immunolocalized in sectioned stigmas.

KEY RESULTS

Protogynous Magnolia flowers had a precise window of stigmatic receptivity, which is concomitant with the secretion of two AGPs labeled for different epitopes. After pollen germination and tube growth, these two AGPs could no longer be detected in the stigmas, suggesting that these AGPs interact with the growing male gametophytes and could be markers of stigmatic receptivity.

CONCLUSIONS

These results show that the period of stigmatic receptivity is finely coordinated with the secretion of two arabinogalactan proteins on stigmas of flowers of M. virginiana. This first report of AGP presence in stigmatic tissues in a member of the magnoliids, together with recently described similar patterns in eudicots, monocots, and members of early-divergent lineages of flowering plants, suggests an ancient and widespread role for AGPs on stigmatic receptivity in angiosperms.

The long goodbye: the rise and fall of flavodoxin during plant evolution

Ferredoxins are electron shuttles harbouring iron-sulfur clusters that connect multiple oxido-reductive pathways in organisms displaying different lifestyles. Some prokaryotes and algae express an isofunctional electron carrier, flavodoxin, which contains flavin mononucleotide as cofactor. Both proteins evolved in the anaerobic environment preceding the appearance of oxygenic photosynthesis. The advent of an oxygen-rich atmosphere proved detrimental to ferredoxin owing to iron limitation and oxidative damage to the iron-sulfur cluster, and many microorganisms induced flavodoxin expression to replace ferredoxin under stress conditions. Paradoxically, ferredoxin was maintained throughout the tree of life, whereas flavodoxin is absent from plants and animals. Of note is that flavodoxin expression in transgenic plants results in increased tolerance to multiple stresses and iron deficit, through mechanisms similar to those operating in microorganisms. Then, the question remains open as to why a trait that still confers plants such obvious adaptive benefits was not retained. We compare herein the properties of ferredoxin and flavodoxin, and their contrasting modes of expression in response to different environmental stimuli. Phylogenetic analyses suggest that the flavodoxin gene was already absent in the algal lineages immediately preceding land plants. Geographical distribution of phototrophs shows a bias against flavodoxin-containing organisms in iron-rich coastal/freshwater habitats. Based on these observations, we propose that plants evolved from freshwater macroalgae that already lacked flavodoxin because they thrived in an iron-rich habitat with no need to back up ferredoxin functions and therefore no selective pressure to keep the flavodoxin gene. Conversely, ferredoxin retention in the plant lineage is probably related to its higher efficiency as an electron carrier, compared with flavodoxin. Several lines of evidence supporting these contentions are presented and discussed.

Chloroplast phylogenomic analysis resolves deep-level relationships within the green algal class Trebouxiophyceae

Background

The green algae represent one of the most successful groups of photosynthetic eukaryotes, but compared to their land plant relatives, surprisingly little is known about their evolutionary history. This is in great part due to the difficulty of recognizing species diversity behind morphologically similar organisms. The Trebouxiophyceae is a species-rich class of the Chlorophyta that includes symbionts (e.g. lichenized algae) as well as free-living green algae. Members of this group display remarkable ecological variation, occurring in aquatic, terrestrial and aeroterrestrial environments. Because a reliable backbone phylogeny is essential to understand the evolutionary history of the Trebouxiophyceae, we sought to identify the relationships among the major trebouxiophycean lineages that have been previously recognized in nuclear-encoded 18S rRNA phylogenies. To this end, we used a chloroplast phylogenomic approach.

Results

We determined the sequences of 29 chlorophyte chloroplast genomes and assembled amino acid and nucleotide data sets derived from 79 chloroplast genes of 61 chlorophytes, including 35 trebouxiophyceans. The amino acid- and nucleotide-based phylogenies inferred using maximum likelihood and Bayesian methods and various models of sequence evolution revealed essentially the same relationships for the trebouxiophyceans. Two major groups were identified: a strongly supported clade of 29 taxa (core trebouxiophyceans) that is sister to the Chlorophyceae + Ulvophyceae and a clade comprising the Chlorellales and Pedinophyceae that represents a basal divergence relative to the former group. The core trebouxiophyceans form a grade of strongly supported clades that include a novel lineage represented by the desert crust alga Pleurastrosarcina brevispinosa. The assemblage composed of the Oocystis and Geminella clades is the deepest divergence of the core trebouxiophyceans. Like most of the chlorellaleans, early-diverging core trebouxiophyceans are predominantly planktonic species, whereas core trebouxiophyceans occupying more derived lineages are mostly terrestrial or aeroterrestrial algae.

Conclusions

Our phylogenomic study provides a solid foundation for addressing fundamental questions related to the biology and ecology of the Trebouxiophyceae. The inferred trees reveal that this class is not monophyletic; they offer new insights not only into the internal structure of the class but also into the lifestyle of its founding members and subsequent adaptations to changing environments.

Phylogeny and divergence times of gymnosperms inferred from single-copy nuclear genes

Phylogenetic reconstruction is fundamental to study evolutionary biology and historical biogeography. However, there was not a molecular phylogeny of gymnosperms represented by extensive sampling at the genus level, and most published phylogenies of this group were constructed based on cytoplasmic DNA markers and/or the multi-copy nuclear ribosomal DNA. In this study, we use LFY and NLY, two single-copy nuclear genes that originated from an ancient gene duplication in the ancestor of seed plants, to reconstruct the phylogeny and estimate divergence times of gymnosperms based on a complete sampling of extant genera. The results indicate that the combined LFY and NLY coding sequences can resolve interfamilial relationships of gymnosperms and intergeneric relationships of most families. Moreover, the addition of intron sequences can improve the resolution in Podocarpaceae but not in cycads, although divergence times of the cycad genera are similar to or longer than those of the Podocarpaceae genera. Our study strongly supports cycads as the basal-most lineage of gymnosperms rather than sister to Ginkgoaceae, and a sister relationship between Podocarpaceae and Araucariaceae and between Cephalotaxaceae-Taxaceae and Cupressaceae. In addition, intergeneric relationships of some families that were controversial, and the relationships between Taxaceae and Cephalotaxaceae and between conifers and Gnetales are discussed based on the nuclear gene evidence. The molecular dating analysis suggests that drastic extinctions occurred in the early evolution of gymnosperms, and extant coniferous genera in the Northern Hemisphere are older than those in the Southern Hemisphere on average. This study provides an evolutionary framework for future studies on gymnosperms.

Overly simplistic substitution models obscure green plant phylogeny

Phylogenetic analysis is an increasingly common and valuable component of plant science. Knowledge of the phylogenetic relationships between plant groups is a prerequisite for understanding the origin and evolution of important plant features, and phylogenetic analysis of individual genes and gene families provides fundamental insights into how those genes and their functions evolved. However, despite an active research community exploring and improving phylogenetic methods, the analytical methods commonly used, and the phylogenetic results they produce, are accorded far more confidence than they warrant. In this opinion article, I emphasise that important parts of the green plant phylogeny are inconsistently resolved and I argue that the lack of consistency arises due to inadequate modelling of changes in the substitution process.

Complete plastome sequences from Glycine syndetika and six additional perennial wild relatives of soybean

Organelle sequences have a long history of utility in phylogenetic analyses. Chloroplast sequences when combined with nuclear data can help resolve relationships among flowering plant genera, and within genera incongruence can point to reticulate evolution. Plastome sequences are becoming plentiful because they are increasingly easier to obtain. Complete plastome sequences allow us to detect rare rearrangements and test the tempo of sequence evolution. Chloroplast sequences are generally considered a nuisance to be kept to a minimum in bacterial artificial chromosome libraries. Here, we sequenced two bacterial artificial chromosomes per species to generate complete plastome sequences from seven species. The plastome sequences from Glycine syndetika and six other perennial Glycine species are similar in arrangement and gene content to the previously published soybean plastome. Repetitive sequences were detected in high frequencies as in soybean, but further analysis showed that repeat sequence numbers are inflated. Previous chloroplast-based phylogenetic trees for perennial Glycine were incongruent with nuclear gene-based phylogenetic trees. We tested whether the hypothesis of introgression was supported by the complete plastomes. Alignment of complete plastome sequences and Bayesian analysis allowed us to date putative hybridization events supporting the hypothesis of introgression and chloroplast "capture."

Did meiosis evolve before sex and the evolution of eukaryotic life cycles?

Biologists have long theorized about the evolution of life cycles, meiosis, and sexual reproduction. We revisit these topics and propose that the fundamental difference between life cycles is where and when multicellularity is expressed. We develop a scenario to explain the evolutionary transition from the life cycle of a unicellular organism to one in which multicellularity is expressed in either the haploid or diploid phase, or both. We propose further that meiosis might have evolved as a mechanism to correct for spontaneous whole-genome duplication (auto-polyploidy) and thus before the evolution of sexual reproduction sensu stricto (i.e. the formation of a diploid zygote via the fusion of haploid gametes) in the major eukaryotic clades. In addition, we propose, as others have, that sexual reproduction, which predominates in all eukaryotic clades, has many different advantages among which is that it produces variability among offspring and thus reduces sibling competition.

Mitochondrial phylogenomics of early land plants: mitigating the effects of saturation, compositional heterogeneity, and codon-usage bias

Phylogenetic analyses using concatenation of genomic-scale data have been seen as the panacea for resolving the incongruences among inferences from few or single genes. However, phylogenomics may also suffer from systematic errors, due to the, perhaps cumulative, effects of saturation, among-taxa compositional (GC content) heterogeneity, or codon-usage bias plaguing the individual nucleotide loci that are concatenated. Here, we provide an example of how these factors affect the inferences of the phylogeny of early land plants based on mitochondrial genomic data. Mitochondrial sequences evolve slowly in plants and hence are thought to be suitable for resolving deep relationships. We newly assembled mitochondrial genomes from 20 bryophytes, complemented these with 40 other streptophytes (land plants plus algal outgroups), compiling a data matrix of 60 taxa and 41 mitochondrial genes. Homogeneous analyses of the concatenated nucleotide data resolve mosses as sister-group to the remaining land plants. However, the corresponding translated amino acid data support the liverwort lineage in this position. Both results receive weak to moderate support in maximum-likelihood analyses, but strong support in Bayesian inferences. Tests of alternative hypotheses using either nucleotide or amino acid data provide implicit support for their respective optimal topologies, and clearly reject the hypotheses that bryophytes are monophyletic, liverworts and mosses share a unique common ancestor, or hornworts are sister to the remaining land plants. We determined that land plant lineages differ in their nucleotide composition, and in their usage of synonymous codon variants. Composition heterogeneous Bayesian analyses employing a nonstationary model that accounts for variation in among-lineage composition, and inferences from degenerated nucleotide data that avoid the effects of synonymous substitutions that underlie codon-usage bias, again recovered liverworts being sister to the remaining land plants but without support. These analyses indicate that the inference of an early-branching moss lineage based on the nucleotide data is caused by convergent compositional biases. Accommodating among-site amino acid compositional heterogeneity (CAT-model) yields no support for the optimal resolution of liverwort as sister to the rest of land plants, suggesting that the robust inference of the liverwort position in homogeneous analyses may be due in part to compositional biases among sites. All analyses support a paraphyletic bryophytes with hornworts composing the sister-group to tracheophytes. We conclude that while genomic data may generate highly supported phylogenetic trees, these inferences may be artifacts. We suggest that phylogenomic analyses should assess the possible impact of potential biases through comparisons of protein-coding gene data and their amino acid translations by evaluating the impact of substitutional saturation, synonymous substitutions, and compositional biases through data deletion strategies and by analyzing the data using heterogeneous composition models. We caution against relying on any one presentation of the data (nucleotide or amino acid) or any one type of analysis even when analyzing large-scale data sets, no matter how well-supported, without fully exploring the effects of substitution models.

Plastid phylogenomics and green plant phylogeny: almost full circle but not quite there

A study in BMC Evolutionary Biology represents the most comprehensive effort to clarify the phylogeny of green plants using sequences from the plastid genome. This study highlights the strengths and limitations of plastome data for resolving the green plant phylogeny, and points toward an exciting future for plant phylogenetics, during which the vast and largely untapped territory of nuclear genomes will be explored.

From algae to angiosperms-inferring the phylogeny of green plants (Viridiplantae) from 360 plastid genomes

BACKGROUND

Next-generation sequencing has provided a wealth of plastid genome sequence data from an increasingly diverse set of green plants (Viridiplantae). Although these data have helped resolve the phylogeny of numerous clades (e.g., green algae, angiosperms, and gymnosperms), their utility for inferring relationships across all green plants is uncertain. Viridiplantae originated 700-1500 million years ago and may comprise as many as 500,000 species. This clade represents a major source of photosynthetic carbon and contains an immense diversity of life forms, including some of the smallest and largest eukaryotes. Here we explore the limits and challenges of inferring a comprehensive green plant phylogeny from available complete or nearly complete plastid genome sequence data.

RESULTS

We assembled protein-coding sequence data for 78 genes from 360 diverse green plant taxa with complete or nearly complete plastid genome sequences available from GenBank. Phylogenetic analyses of the plastid data recovered well-supported backbone relationships and strong support for relationships that were not observed in previous analyses of major subclades within Viridiplantae. However, there also is evidence of systematic error in some analyses. In several instances we obtained strongly supported but conflicting topologies from analyses of nucleotides versus amino acid characters, and the considerable variation in GC content among lineages and within single genomes affected the phylogenetic placement of several taxa.

CONCLUSIONS

Analyses of the plastid sequence data recovered a strongly supported framework of relationships for green plants. This framework includes: i) the placement of Zygnematophyceace as sister to land plants (Embryophyta), ii) a clade of extant gymnosperms (Acrogymnospermae) with cycads + Ginkgo sister to remaining extant gymnosperms and with gnetophytes (Gnetophyta) sister to non-Pinaceae conifers (Gnecup trees), and iii) within the monilophyte clade (Monilophyta), Equisetales + Psilotales are sister to Marattiales + leptosporangiate ferns. Our analyses also highlight the challenges of using plastid genome sequences in deep-level phylogenomic analyses, and we provide suggestions for future analyses that will likely incorporate plastid genome sequence data for thousands of species. We particularly emphasize the importance of exploring the effects of different partitioning and character coding strategies.