BACKGROUND
viridiplantae, or green plants, are a clade of perhaps  <dig>  species  <cit>  that exhibit an astounding diversity of life forms, including some of the smallest and largest eukaryotes  <cit> .
fossil evidence suggests the clade is at least  <dig> million years old  <cit> , while divergence time estimates from molecular data suggest it may be more than one billion years old  <cit> .
reconstructing the phylogenetic relationships across green plants is challenging because of the age of the clade, the extinction of major lineages  <cit> , and extreme molecular rate and compositional heterogeneity  <cit> .
most phylogenetic analyses of viridiplantae have recovered two well-supported subclades, chlorophyta and streptophyta <cit> .
chlorophyta contain most of the traditionally recognized “green algae,” and streptophyta contain the land plants , as well as several other lineages also considered “green algae”.
land plants include the seed plants , which consist of ~ <dig>  to ~ <dig>  species  <cit> .
while many of the major green plant clades are well defined, questions remain regarding the relationships among them.
for example, the closest relatives of land plants have varied among analyses  <cit> , as have the relationships among the three bryophyte lineages   <cit> .
the relationships among extant gymnosperms also remain contentious, particularly with respect to the placement of gnetophyta <cit> .
most broad analyses of green plant relationships based on nuclear gene sequence data have relied largely on 18s/26s rdna sequences  <cit> , although recent analyses have employed numerous nuclear genes  <cit> .
some studies have used mitochondrial gene sequence data, often in combination with other data  <cit> .
however, investigations of green plant phylogeny typically have either largely or exclusively employed chloroplast genes .
sequence data from the plastid genome have transformed plant systematics and contributed greatly to the current view of plant relationships.
with the plastid genome present in high copy numbers in each cell in most plants, and with relatively little variation in gene content and order  <cit> , as well as few reported instances of gene duplication or horizontal gene transfer  <cit> , the plastid genome provides a wealth of phylogenetically informative data that are relatively easy to obtain and use  <cit> .
although early phylogenetic studies using one or a few chloroplast loci provided fundamental insights into relationships within and among green plant clades, these analyses failed to resolve some backbone relationships  <cit> .
these remaining enigmatic portions of the green plant tree of life ultimately motivated the use of entire, or nearly entire, plastid genome sequences for phylogenetic inference.
complete sequencing of the relatively small  plastid genome has been technically feasible since the mid-1980s  <cit> , although few plastid genomes were sequenced prior to  <dig> .
next-generation sequencing  technologies, such as  <dig>  <cit>  and illumina  <cit> , greatly reduced the cost and difficulty of sequencing plastid genomes, and consequently, the number of plastid genomes available on genbank increased nearly six-fold from  <dig> to  <dig>  <cit> .
phylogenetic analyses based on complete plastid genome sequences have provided valuable insights into relationships among and within subclades across the green plant tree of life .
still, studies employing complete plastid genomes generally have either focused on subclades of green plants or have had relatively low taxon sampling.
thus, they have not addressed the major relationships across all green plants simultaneously.
we assembled available plastid genome sequences to build a phylogenetic framework for viridiplantae that reflects the wealth of new plastid genome sequence data.
furthermore, we highlight analytical challenges for resolving the green plant tree of life with this type of data.
we performed phylogenetic analyses of protein-coding data on  <dig> genes from  <dig> taxa, exploring the effects of different partitioning and character-coding protocols for the entire data set as well as subsets of the data.
while our analyses recover many well-supported relationships and reveal strong support for some contentious relationships, several factors, including base composition biases, can affect the results.
we also highlight the challenges of using plastid genome data in deep-level phylogenomic analyses and provide suggestions for future analyses that will incorporate plastid genome data for thousands of species.
RESULTS
data set
we assembled plastid protein-coding sequences from  <dig> species  for which complete or nearly complete plastid genome sequences were available on genbank.
of the  <dig> species, there were  <dig> angiosperms ,  <dig> gymnosperms , seven monilophytes , four lycophytes , three liverworts , one hornwort , two mosses , six taxa from the paraphyletic streptophytic algae, and  <dig> chlorophytic algae .
the phylogenetic character matrices contained sequences from  <dig> genes and the following number of alignment positions:  <dig>  bp for the matrix containing all nucleotide positions  and the ry-coded  version of the ntall matrix;  <dig>  bp in the matrix containing only the first and second codon positions , and  <dig>  amino acids .
the number of genes present per taxon varied from  <dig> to  <dig> , while the number of taxa present per gene ranged from  <dig> to  <dig> .
taxa with few genes present, such as helicosporidium  and rhizanthella , represent highly modified complete plastid genomes of non-photosynthetic species  <cit> .
the percentage of missing data  was ~ <dig> % for each of the four data sets.
the pattern of data across each of the four matrices is decisive, meaning that it can uniquely define a single tree for all taxa  <cit> .
the data contain 100% of all possible triplets of taxa, and are decisive for 100% of all possible trees.
all alignments have been deposited in the dryad data repository  <cit> .
gc bias
gc content varied considerably both among lineages and also within single genomes, and chi-square tests rejected the null hypothesis of homogeneous base frequencies .
the average gc content in the ntall matrix was  <dig> %, and it ranged from  <dig> % in selaginella uncinata to  <dig> % in helicosporidium sp.
.
also, the average gc content varied among first, second, and third codon positions, with by far the most variation among lineages at the third codon position .
although there was extensive heterogeneity in gc content across all species, there was relatively little variation among the seed plant taxa .
there also was significant correlation between nucleotide composition and amino acid composition.
plastid genomes that are gc-rich had a significantly higher percentage  of amino acids that are encoded by gc-rich codons .
similarly, gc-rich plastid genomes had a significantly lower percentage  of amino acids that are coded by at-rich codons .
phylogenetic analyses
in the phylogenetic analyses of all data sets and partitioning schemes, the partitioning strategy with the most partitions consistently fit the data best based on the aicc .
these best-fit models partitioned the aa matrix by gene  and the nucleotide  and ry matrices by codon position and gene .
all a posteriori bootstopping analyses indicated that convergence of support values had been reached after  <dig> replicates, and thus our choice of  <dig> replicates was more than sufficient to obtain reliable bootstrap values.
partitioning strategies judged to be the best by the aicc are in bold.
we will focus on reporting the relationships of major clades of viridiplantae shown in the 50% maximum likelihood  majority-rule bootstrap consensus summary trees for each data set: ntall , ntno3rd , ry , and aa .
these summary trees collapse some clades for ease of viewing the major relationships within viridiplantae.
a summary of important results and conflicts among these four data sets is given in table  <dig>  we provide full majority-rule bootstrap consensus trees for the ntall , ntno3rd , ry , and aa  data sets.
ml trees with branch lengths and bs values are also provided: ntall , ntno3rd , ry , and aa .
average support values among all internal nodes in the ml trees were slightly higher in the ntall phylogeny  compared to the other data sets .
the ntall phylogeny also had the most clades resolved with ≥ 70% bs  while the ntno3rd, ry, and aa data sets had 87%, 87%, and 86% of the possible bipartitions resolved at ≥ 70% bs, respectively.
all resulting trees have been deposited in the dryad data repository  <cit> .
bootstrap support  values >50% are shown as percentages.
when sister groups for the taxon of interest are listed, bootstrap support  values on the left are for the clade including the taxon of interest and its sister group within viridiplantae, while bs values on the right are for the more inclusive clade excluding the taxon of interest.
if only one bs value is given for a sister relationship, only two terminals are involved .
the monophyly of chlorophyta receives 100% bs in all analyses.
prasinophyceae are consistently not monophyletic.
instead, the prasinophyte nephroselmis is sister to all other chlorophyta , while remaining prasinophyceae form a clade that is variously supported  and is sister to a clade of the remaining chlorophyta.
chlorophyceae are monophyletic , but trebouxiophyceae and ulvophyceae are not monophyletic, and the relationship of chlorophyceae to these lineages is unresolved.
we consistently recovered a single set of relationships among the streptophytic algae subtending the land plant clade.
zygnematophyceae are sister to land plants, coleochaetophyceae are sister to zygnematophyceae + embryophyta, charophyceae are sister to coleochaetophyceae + , and a clade of mesostigmatophyceae + chlorokybophyceae is sister to all other streptophyta.
each of these relationships has ≥86% bs support .
the branching order of the non-vascular land plant lineages differs among analyses.
in analyses of the ntall and ry data sets, marchantiophyta , followed by bryophyta , and then anthocerotophyta  are the earliest-branching land plant lineages, with anthocerotophyta the immediate sister to the vascular plants .
in the ntall and ry analyses, these relationships had ≥89% bs support except for the bryophyta +  relationship in the ntall analysis, which received only 69% bs .
in contrast, in the ntno3rd and aa analyses, bryophyta and marchantiophyta formed a clade , followed by anthocerophyta as sister to tracheophyta .
within tracheophyta, the ntno3rd, ry, and aa data sets all place lycopodiophyta sister to a euphyllophyta clade .
however, the analysis of the ntall data set places monilophyta sister to a clade of lycopodiophyta + spermatophyta .
our analyses of monilophyta generally reveal strong support for a clade of equisetales + psilotales as sister to marattiales + leptosporangiate ferns .
the lowest support obtained was for equisetales + psilotales in the ntno3rd analysis  and ntall ; all other nodes in all analyses received > 90% bs, with marattiales + leptosporangiate ferns receiving ≥ 99% bs.
within spermatophyta, all analyses place the extant gymnosperms  sister to angiospermae with 100% bs.
within extant gymnosperms, cycadales and ginkgoales form a clade  that is sister to a clade in which gnetophyta  are nested within the paraphyletic conifers.
there is generally high support  placing gnetophyta as sister to a clade of araucariales + cupressales.
this “gnecup” clade  is then sister to pinales, which has 100% bs in all analyses.
in all analyses, angiospermae receive 100% bs, and amborella  is sister to all other angiosperms, followed by nymphaeales, and then austrobaileyales.
these relationships are mostly supported by 100% bs.
however, nymphaeales +  receives 81% bs  in the ntno3rd analyses and 70% bs  in the aa analyses.
the remaining angiosperms  receive 100% bs in all analyses.
within mesangiospermae, the relationships among monocotyledoneae, magnoliidae, eudicotyledoneae, and ceratophyllum  are not well supported and vary depending on the analysis.
the strongest support for the placement of ceratophyllales is 75% bs as sister to eudicotyledoneae in the ry analysis .
chloranthales receive 61-69% bs as sister to the well-supported  magnoliidae.
however, magnoliidae are not monophyletic in the aa analyses, where piperales are sister to ceratophyllales .
within the monocot clade , acorales, followed by alismatales, have 100% bs in all analyses as subsequent sisters to the remaining monocots.
in three of our analyses , a variously supported clade  of liliales +  is sister to a clade  of the remaining monocots .
however, in the ry-coded analysis, pandanales + dioscoreales  is sister to a clade of liliales + , which receives 69% bs .
here asparagales + commelinidae is supported by 80% bs.
within the eudicots , which receive 100% bs in all analyses, ranunculales are sister to the remaining taxa.
in the ntall, ntno3rd, ry, and aa analyses, the clade of these remaining taxa receives 100%, 85%, 100%, and 62% bs, respectively.
relationships vary among sabiaceae, proteales, and a clade of the remaining taxa, depending on the analysis.
in the ntall and ntno3rd analyses, proteales + sabiaceae are supported as a clade, although with only 63% and 60% bs, respectively.
however, in the ry analysis, proteales are sister to a clade containing sabiaceae plus the remaining taxa, which has 79% bs.
in the aa analysis, relationships among these three clades are unresolved.
among the remaining eudicots, we consistently recovered trochodendrales as sister to buxales + pentapetalae and gunnerales as sister to the remaining lineages of pentapetalae: dilleniaceae, superrosidae, and superasteridae.
the placement of dilleniaceae remains uncertain.
the family is sister to superrosidae in the ntall , ntno3rd , and ry  analyses, but appears as sister to superasteridae  in the aa analysis.
within superrosidae, a clade of vitales + saxifragales is supported in the ntall , ntno3rd , and aa  analyses.
in the ry analysis, the relationship among saxifragales, vitales, and remaining rosidae  is unresolved.
fabidae and malvidae are both recovered with ≥ 99% bs in the ntall and ry analyses.
however, each clade receives only 70% bs in the ntno3rd analysis.
in the aa analysis neither clade is monophyletic; zygophyllales are embedded  within a clade of malvidae taxa.
the com clade  is sister to a clade of fagales, cucurbitales, rosales, and fabales in fabidae in the aa , ry , and ntall  trees and forms a trichotomy with zygophyllales and the clade of fagales, cucurbitales, rosales, and fabales in the ntno3rd tree .
zygophyllales are sister to geraniales  in the aa tree and sister to all other fabidae in the ntall and ry trees .
superasteridae  are recovered in all analyses.
this clade receives 100% bs in the ntall and ry analyses, 95% bs in the ntno3rd analysis, and 66% bs in the aa analysis.
santalales and berberidopsidales are strongly supported as subsequent sisters to caryophyllales + asteridae.
within asteridae, cornales, followed by ericales, are subsequent sisters to a strongly supported clade that comprises strongly supported campanulidae and lamiidae clades.
within lamiidae, the placement of boraginaceae is weak among the various analyses.
boraginaceae are sister to gentianales  in the aa tree, part of a trichotomy  with lamiales and solanales + gentianales in the ntall tree, and sister to a weakly supported clade including gentianales, lamiales, and solanales in the ntno3rd  and ry  trees.
analysis of only the third codon positions  resulted in several very strong conflicts along the backbone of viridiplantae when compared to the topology from the ntno3rd analyses.
these conflicts include the backbone relationships within chlorophyta, the placements of cycadales and lycopodiophyta, the relationships of the three major bryophyte lineages, and backbone relationships within poales.
removal of four taxa  with elevated rates of molecular evolution and few genes present in the data sets did not significantly affect the resulting topologies.
discussion
while the enormous phylogenetic data sets that result from new genome or transcriptome sequencing efforts can ameliorate the effects of random or stochastic error, they also may exacerbate the effects of systematic error, or error resulting from problems in the analysis, such as model inaccuracy.
the high amount of agreement among our various analyses and strong support for results generally consistent with previous studies  suggest that plastid genome sequence data hold much promise for resolving relationships throughout the green plants.
however, several areas of conflict between analyses using different character-coding strategies demonstrate that plastid genome phylogenetics is also susceptible to systematic error.
here we evaluate the phylogenetic results, emphasizing areas of agreement and concern, and then address some of the methodological issues raised by our results.
evaluation of phylogenetic relationships
historically, chlorophyta have been divided into prasinophyceae, trebouxiophyceae, chlorophyceae, and ulvophyceae based on the ultrastructure of the flagellar apparatus and features related to cytokinesis  <cit> .
the current status of green algae phylogenetics  has been reviewed recently  <cit> .
the most comparable study to ours in terms of data and taxon sampling is by lang and nedelcu  <cit> , who constructed a phylogeny of green algae with plastid genome sequence data.
however, they analyzed only an amino acid data set using bayesian inference and the cat model  <cit> .
we found a paraphyletic prasinophyceae , which agrees with previous molecular analyses  <cit> .
however, lang and nedelcu  <cit>  recovered a monophyletic prasinophyceae, albeit with little support.
chlorophyceae are monophyletic , which agrees with the results of lang and nedelcu  <cit> .
we also find that trebouxiophyceae and ulvophyceae are not monophyletic, and that the relationship of chlorophyceae to these lineages is unresolved.
the branching order of the various trebouxiophyceae, ulvophyceae, and chlorophyceae lineages within chlorophyta, unresolved in our analyses, was also uncertain in earlier analyses .
similarly in lang and nedelcu  <cit> , trebouxiophyceae and ulvophyceae were not supported as monophyletic, although unlike our results, almost all nodes in their phylogeny were maximally supported.
our analyses provide consistent, strong support for the relationships of streptophytic algae to land plants, and all analyses support zygnematophyceae as the sister to land plants .
relationships among these lineages and the closest relatives of land plants have varied in previous studies depending on taxon sampling and gene choice.
some studies agree with our results placing zygnematophyceae as sister to land plants  <cit> , while other phylogenetic analyses indicate that charophyceae <cit>  or coleochaetophyceae <cit>  occupy this position.
depending on the analysis, zhong et al.
<cit>  found either zygnematophyceae alone or a clade of zygnematophyceae + coleochaetophyceae as sister to land plants.
in particular, the results of lang and nedelcu  <cit>  conflict with our results regarding the sister group to embryophyta.
while we find a clade of coleochaetophyceae + , their results strongly support zygnematophyceae + .
phylogenetic relationships among bryophytes  are also contentious, and nearly every possible relationship among these lineages has been reported, often with strong support.
most studies have shown the bryophytes as paraphyletic with respect to tracheophyta rather than as a clade  <cit> .
as recovered in our ntall and ry analyses , liverworts  often are placed sister to all other land plants, followed by mosses , and with hornworts  sister to tracheophyta <cit> .
a sister relationship between mosses and liverworts, found in our ntno3rd and aa analyses , was proposed previously based on morphological  <cit>  and molecular data  <cit>  and has been recovered with numerous nuclear genes .
this relationship was also recovered in analyses of complete plastid genome data by karol et al.
<cit>  when divergent taxa  were excluded from phylogenetic analyses and also by wolf and karol  <cit>  when third positions were excluded.
our results placing lycopodiophyta sister to euphyllophyta in all but the ntall analysis agree with most molecular phylogenetic analyses  <cit> .
this split is also supported by analyses of morphological characters in fossil  <cit>  and extant taxa  <cit> .
monilophyta and spermatophyta also possess a 30-kb inversion in the large single-copy region of the plastid genome not found in lycopodiophyta and the three bryophyte clades  <cit> .
in the ntall analysis, euphyllophyta are not monophyletic ; lycopodiophyta, rather than monilophyta, are sister to spermatophyta.
this relationship has been reported previously  <cit> ; however, it likely is a phylogenetic artifact, perhaps related to base composition bias .
the plastid genome of the lycophyte selaginella has an especially high gc content  <cit> , with selaginella unicata having the highest gc content in our ntall data set .
in some previous studies, relationships among lineages of monilophyta have not been well resolved or supported .
as a result, the relationships among equisetales, psilotales, marattiales, and leptosporangiate ferns are often represented as a polytomy .
in contrast, most of our analyses recovered strong support for a clade of equisetales + psilotales as sister to marattiales + leptosporangiate ferns .
these relationships agree with recent studies of monilophyte relationships based on plastid genome sequence data  <cit> , although support is stronger here.
unfortunately, ophioglossales, which often appear as sister to psilotales, lacked a sequenced plastome at the time of our analyses.
however, plastid genome data for ophioglossales have subsequently been published and analyzed in a phylogenetic context  <cit> , with strong support for ophioglossales as sister to psilotales and weak support for this clade as sister to equisetales.
results from that study with regard to marattiales and leptosporangiate ferns agree with the relationships presented here.
relationships among the lineages of extant seed plants, and especially the placement of gnetophyta, have long been debated  <cit> .
gnecup trees, found in all of our analyses, were initially recovered by nickrent et al.
<cit> , and then more recently by zhong et al.
<cit> .
however, zhong et al.
<cit>  suggested that the support for gnecup may be the result of long-branch attraction; by removing highly variable proteins, support for gnecup decreased.
furthermore, by removing what they considered parallel substitutions between lineages leading to gnetophyta and to cryptomeria , a gnepine topology was recovered.
although several different placements for gnetophyta have been recovered and strongly supported, many studies involving multiple genes have placed gnetophyta sister to pinales .
using both coalescent and concatenation analyses, xi et al.
<cit>  found that the phylogenetic placement of gnetophyta differs between the nuclear and plastid genomes.
in their analyses using nuclear data, the gnepine hypothesis is supported, while their analyses of plastid data support the gnecup hypothesis.
in contrast, lee et al.
<cit>  found strong support for gnetophyta sister to the remaining gymnosperms  in an ml analysis of  <dig>  sets of nuclear gene orthologs from  <dig> land plant genera.
the backbone relationships among angiosperm  lineages generally agree with results from recent analyses, including a 17-gene analysis of  <dig> angiosperms  <cit>  and previous analyses of plastid genome data sets  <cit> .
the position of ceratophyllum , and thus the relationships among monocotyledoneae, eudicotyledoneae, and magnoliidae, varies among our analyses, although without strong support.
this contrasts with several other large, multi-gene analyses in which monocotyledoneae are sister to ceratophyllales + eudicotyledoneae <cit> .
interestingly, the strongest support for the placement of ceratophyllales sister to eudicotyledoneae is in the ry analysis .
however, in that analysis, the relationships among ceratophyllales + eudicotyledoneae, monocotyledoneae, and magnoliidae are unresolved.
within the angiosperms, some relationships that have been uncertain, particularly at deep levels , receive moderate to strong support in at least some of our analyses.
for example, the placement of myrtales and geraniales in the malvidae is supported with 70% bs  in the ntno3rd tree and ≥ 99% bs in the ry  and ntall  trees.
myrtales and geraniales are also placed in a clade with the malvidae taxa in the aa analysis ; however, zygophyllales are also included within this clade, making malvidae non-monophyletic.
likewise, chloranthales are sister to magnoliidae in all trees, but with weaker support .
in two cases, all analyses but ry resolve relationships , with ry producing a polytomy that does not conflict with the resolutions found in the other analyses.
these two cases are as follows:  vitales + saxifragales supported by ≥ 70% bs in all analyses but ry, with saxifragales, vitales, and remaining rosidae forming a polytomy in the ry tree ;  dasypogonaceae + arecales in all but ry  and a trichotomy of dasypogonaceae, arecales, and poales +  in the ry tree .
in two additional cases when ry is compared to the other three analyses, the ry analysis produced either stronger support for the placement of a taxon or a different placement altogether.
first, in the ntall, ntno3rd, and aa analyses, the position of sabiaceace among the early-diverging lineages of eudicotyledoneae is weakly supported.
however, in the ry analysis, sabiaceae receive moderate support  as sister to a strongly supported  clade of trochodendrales + ).
this contrasts with previous studies that often place sabiaceae as sister to proteales <cit> .
an example of a different placement of a taxon in the ry analysis when compared to the other analyses involves liliales.
the ntall, ntno3rd, and the aa analyses support liliales as sister to a clade of dioscoreales + pandanales with 72%, 69%, and 80% bs, respectively.
this placement of liliales was also recovered in barrett et al.
<cit> .
in contrast, in the ry analysis, liliales are placed in a clade with asparagales + commelinidae with moderate support .
this latter placement of liliales was strongly supported in an analysis with much better taxon sampling  <cit> .
some taxa that have been problematic in previous studies  continue to defy definitive placement.
their positions vary among our analyses, although they are generally not well supported in some, or all, of the trees.
despite its general placement of the com clade in fabidae in these and other plastid analyses, this clade is more closely related to malvidae in some analyses, particularly those using mitochondrial gene sequences .
recent analyses of plastid, mitochondrial, and nuclear data suggest that the com clade may represent ancient reticulation involving fabidae and malvidae during the rapid radiation of rosidae .
methodological issues of plastid phylogenomic analyses
to address potential systematic error in large-scale phylogenetic analyses, scientists often either try to improve the fit of models to the data or change or remove problematic data.
with increasing sequence length and number of genes, it is more likely that a sequence alignment will contain regions with heterogeneous processes of molecular evolution.
we see evidence of this high heterogeneity with our model-fitting experiments, which always favor the most parameter-rich models .
thus, defining partitioning schemes and models that can accurately reflect the true processes of molecular evolution while not over-parameterizing the analysis remains critically important for phylogenetic analyses of large plastid data sets.
although we assessed models that account for heterogeneity in patterns of molecular evolution among genes and in some cases codon positions, our model selection tests only evaluated a small selection of possible models and partitioning schemes.
it is possible that other partitioning schemes could enable simpler models.
most conventional phylogenetic models, like those used in our analyses, also assume homogeneous processes of evolution throughout the tree.
yet when the branches of the phylogeny encompass over one billion years of evolutionary history, as likely do those in the green plants, the patterns of evolution almost certainly differ among lineages and through time.
this is apparent from the often good fit of covarion models  to plastid genes  <cit>  and the presence of nucleotide compositional heterogeneity, which can confound conventional phylogenetic analyses .
also, our models do not account for shifts in selective pressure or instances of positive selection that will affect nucleotide and amino acid substitution patterns .
nucleotide compositional heterogeneity remains a concern for green plant plastid genome analyses.
this variation is most evident in non-seed plant taxa , and thus it has not been a focus of many previous phylogenetic analyses of plastid genome sequences.
a gc bias in itself is not necessarily problematic for phylogenetic analyses, but nearly all commonly used models for likelihood-based phylogenetic analyses assume single equilibrium nucleotide frequencies.
given that gc content appears to vary by codon position in plants   <cit> , a partitioning scheme that estimates separate nucleotide frequencies for each codon position may account for some of the spatial heterogeneity in gc content in the plastid genome, but it does not address the differences in gc frequency among lineages.
a commonly used strategy to reduce the effects of gc heterogeneity across lineages is ry-coding, in which the purines  are coded as rs and the pyrimidines  are coded as ys  <cit> .
ry-coding can reduce the compositional variability among lineages, improve the fit of models, and increase the signal for internal branches  <cit> .
an obvious disadvantage to ry-coding is that by coding the sequences with two character states instead of four, it reduces the amount of information in the sequences.
in general, we see little overall reduction, and even some gains, in bootstrap support when using ry-coding compared to the use of all nucleotide data , suggesting that the benefits of ry-coding make up for any potential costs of information loss.
perhaps the biggest topological difference in the ry phylogeny  compared to ntall  is the placement of monilophyta rather than lycopodiophyta as sister to seed plants.
the unexpected placement of lycopodiophyta as the sister to seed plants in the ntall analysis  is almost certainly an artifact of systematic error; several other lines of evidence support monilophyta as the sister group of seed plants .
approaches to reducing systematic errors by excluding problematic data, which often include fast-evolving or saturated sites, also have been suggested for plastid genome analyses  <cit> .
with the proper model of molecular evolution and adequate taxon sampling, fast sites are not necessarily problematic; they are only problematic insofar as they are difficult to model.
yet with heterogeneous processes of molecular evolution throughout the tree, the fast-evolving or saturated sites can produce a significant non-phylogenetic signal .
indeed, the third codon positions appear to have especially high levels of compositional heterogeneity, potentially causing systematic error , and an analysis of just the third codon positions  conflicts with the analyses of other data sets in several critical parts of the tree .
however, third codon positions also represent a large proportion of the variable sites in the alignment, and removing them may exclude much of the phylogenetic information in some parts of the tree.
with regard to backbone relationships in our phylogeny, excluding the third position sites  produces several interesting changes in contrast to ntall: 1) it supports the sister relationship of mosses and liverworts, 2) monilophytes, not lycophytes, are placed sister to seed plants as expected, and 3) support for some of the backbone angiosperm relationships is reduced.
thus, the effects of removing the third codon position sites appear to vary in different parts of the tree.
another strategy for overcoming potential error associated with fast-evolving sites is to code the sequences as amino acids rather than nucleotides.
this does not necessarily eliminate problems of compositional heterogeneity, as the gc bias also may bias amino acid composition   <cit> .
regarding backbone green plant relationships, the aa analysis provided similar results to analyses of only first and second codon positions.
aa analysis also produced some weakly supported, questionable relationships among angiosperm lineages .
in previous deep-level plant analyses, analyses of amino acid data have resulted in arguably more problematic or questionable relationships than analyses of nucleotide data  <cit> .
however, these results are likely due to inappropriate models of amino acid evolution  <cit> , and with better models, optimized for plastid evolution, amino acid data may be a valuable source of phylogenetic information.
taxon sampling is also important for plastid phylogenomic studies, especially when the model of evolution is inadequate  <cit> , and genome-scale analyses often have limited taxon sampling.
new methods for rapid and inexpensive plastid genome sequencing  may ameliorate the effects of insufficient sampling of extant taxa; however, many major lineages of green plants are now extinct, precluding their inclusion in analyses of molecular data .
in addition, ancient, rapid radiations abound within portions of the green plant tree of life, creating extremely difficult phylogenetic problems no matter the taxon sampling  <cit> .
furthermore, even in the absence of systematic error, it is possible that a tree built from plastid genome data will not reflect species relationships.
the plastid genome represents a single locus of linked genes .
for phylogenetic analyses, this can be beneficial because combining genes with different evolutionary histories into a single character matrix can lead to phylogenetic error  <cit> .
however, incomplete lineage sorting or ancient reticulation could lead to conflict between the plastid gene tree and the species phylogeny  <cit> .
for this reason, it will be interesting to compare phylogenetic hypotheses from the plastid genome with independent phylogenetic estimates from numerous nuclear and mitochondrial loci.
finally, while full plastid genome sequence data provide much power for resolving difficult phylogenetic relationships, it is not clear that they can resolve all plant relationships.
theoretical work suggests that extremely large data sets may be necessary to resolve some relationships when the internal nodes are separated by very short branches  <cit> , and recent analyses indicate that full plastid genomes are not sufficient to reject alternative topologies among monocots  <cit> .
indeed, the unresolved or conflicting parts of the green plant phylogeny in our analyses are generally associated with short internal branch lengths .
thus, even if the model of evolution accurately reflects the true process of molecular evolution, and there is no systematic error, plastid genome data alone may not be sufficient to resolve all parts of the green plant tree of life.
that is, the topology may not be identifiable with the plastid data alone.
a recent analysis using a new diagnostic test for phylogenetic identifiability based on data cloning suggested that a backbone topology of angiosperms was identifiable from plastid sequence data using the gtr + Γ model  <cit> , but the tree in this paper is much larger and the models more complex.
in any case, it will be necessary to include perspectives from the nuclear genome and phenotypic data before we are confident about all deep-level relationships among green plants.
CONCLUSIONS
our diverse analyses provide a first approach to addressing some of the difficult issues associated with plastid phylogenetic analyses at this evolutionary depth and level of taxon sampling.
the results of the analyses using different models, character-coding strategies, and character subsets suggest that much of the tree is robust to many different phylogenetic approaches, and they highlight regions of the tree that need more scrutiny .
more sophisticated modelling approaches may more accurately characterize the heterogeneous processes of molecular evolution, but it is also crucial that the parameters of these complex models can be estimated by the data at hand  <cit> .
while it may be impossible for any model to reflect perfectly the complexities of molecular evolution, as we better characterize these processes it will be possible to examine through simulations their possible effects on phylogenetic analyses and to recognize phylogenetic error caused by model misspecification.

 viridiplantae, or green plants, are a clade of perhaps species that exhibit an astounding diversity of life forms, including some of the smallest and largest eukaryotes. fossil evidence suggests the clade is at least million years old, while divergence time estimates from molecular data suggest it may be more than one billion years old. reconstructing the phylogenetic relationships across green plants is challenging because of the age of the clade, the extinction of major lineages, and extreme molecular rate and compositional heterogeneity. most phylogenetic analyses of viridiplantae have recovered two well-supported subclades, chlorophyta and streptophyta. chlorophyta contain most of the traditionally recognized “green algae,” and streptophyta contain the land plants, as well as several other lineages also considered “green algae”. land plants include the seed plants, which consist of ~ to ~ species. while many of the major green plant clades are well defined, questions remain regarding the relationships among them. for example, the closest relatives of land plants have varied among analyses, as have the relationships among the three bryophyte lineages. the relationships among extant gymnosperms also remain contentious, particularly with respect to the placement of gnetophyta. most broad analyses of green plant relationships based on nuclear gene sequence data have relied largely on 18s/26s rdna sequences, although recent analyses have employed numerous nuclear genes. some studies have used mitochondrial gene sequence data, often in combination with other data. however, investigations of green plant phylogeny typically have either largely or exclusively employed chloroplast genes. sequence data from the plastid genome have transformed plant systematics and contributed greatly to the current view of plant relationships. with the plastid genome present in high copy numbers in each cell in most plants, and with relatively little variation in gene content and order, as well as few reported instances of gene duplication or horizontal gene transfer, the plastid genome provides a wealth of phylogenetically informative data that are relatively easy to obtain and use. although early phylogenetic studies using one or a few chloroplast loci provided fundamental insights into relationships within and among green plant clades, these analyses failed to resolve some backbone relationships. these remaining enigmatic portions of the green plant tree of life ultimately motivated the use of entire, or nearly entire, plastid genome sequences for phylogenetic inference. complete sequencing of the relatively small plastid genome has been technically feasible since the mid-1980s, although few plastid genomes were sequenced prior to. next-generation sequencing technologies, such as and illumina, greatly reduced the cost and difficulty of sequencing plastid genomes, and consequently, the number of plastid genomes available on genbank increased nearly six-fold from to. phylogenetic analyses based on complete plastid genome sequences have provided valuable insights into relationships among and within subclades across the green plant tree of life. still, studies employing complete plastid genomes generally have either focused on subclades of green plants or have had relatively low taxon sampling. thus, they have not addressed the major relationships across all green plants simultaneously. we assembled available plastid genome sequences to build a phylogenetic framework for viridiplantae that reflects the wealth of new plastid genome sequence data. furthermore, we highlight analytical challenges for resolving the green plant tree of life with this type of data. we performed phylogenetic analyses of protein-coding data on genes from taxa, exploring the effects of different partitioning and character-coding protocols for the entire data set as well as subsets of the data. while our analyses recover many well-supported relationships and reveal strong support for some contentious relationships, several factors, including base composition biases, can affect the results. we also highlight the challenges of using plastid genome data in deep-level phylogenomic analyses and provide suggestions for future analyses that will incorporate plastid genome data for thousands of species.  data set we assembled plastid protein-coding sequences from species for which complete or nearly complete plastid genome sequences were available on genbank. of the species, there were angiosperms, gymnosperms, seven monilophytes, four lycophytes, three liverworts, one hornwort, two mosses, six taxa from the paraphyletic streptophytic algae, and chlorophytic algae. the phylogenetic character matrices contained sequences from genes and the following number of alignment positions:  bp for the matrix containing all nucleotide positions and the ry-coded version of the ntall matrix;  bp in the matrix containing only the first and second codon positions, and amino acids. the number of genes present per taxon varied from to, while the number of taxa present per gene ranged from to. taxa with few genes present, such as helicosporidium and rhizanthella, represent highly modified complete plastid genomes of non-photosynthetic species. the percentage of missing data was ~ % for each of the four data sets. the pattern of data across each of the four matrices is decisive, meaning that it can uniquely define a single tree for all taxa. the data contain 100% of all possible triplets of taxa, and are decisive for 100% of all possible trees. all alignments have been deposited in the dryad data repository. gc bias gc content varied considerably both among lineages and also within single genomes, and chi-square tests rejected the null hypothesis of homogeneous base frequencies. the average gc content in the ntall matrix was %, and it ranged from % in selaginella uncinata to % in helicosporidium sp.. also, the average gc content varied among first, second, and third codon positions, with by far the most variation among lineages at the third codon position. although there was extensive heterogeneity in gc content across all species, there was relatively little variation among the seed plant taxa. there also was significant correlation between nucleotide composition and amino acid composition. plastid genomes that are gc-rich had a significantly higher percentage of amino acids that are encoded by gc-rich codons. similarly, gc-rich plastid genomes had a significantly lower percentage of amino acids that are coded by at-rich codons. phylogenetic analyses in the phylogenetic analyses of all data sets and partitioning schemes, the partitioning strategy with the most partitions consistently fit the data best based on the aicc. these best-fit models partitioned the aa matrix by gene and the nucleotide and ry matrices by codon position and gene. all a posteriori bootstopping analyses indicated that convergence of support values had been reached after replicates, and thus our choice of replicates was more than sufficient to obtain reliable bootstrap values. partitioning strategies judged to be the best by the aicc are in bold. we will focus on reporting the relationships of major clades of viridiplantae shown in the 50% maximum likelihood majority-rule bootstrap consensus summary trees for each data set: ntall, ntno3rd, ry, and aa. these summary trees collapse some clades for ease of viewing the major relationships within viridiplantae. a summary of important results and conflicts among these four data sets is given in table  we provide full majority-rule bootstrap consensus trees for the ntall, ntno3rd, ry, and aa data sets. ml trees with branch lengths and bs values are also provided: ntall, ntno3rd, ry, and aa. average support values among all internal nodes in the ml trees were slightly higher in the ntall phylogeny compared to the other data sets. the ntall phylogeny also had the most clades resolved with ≥ 70% bs while the ntno3rd, ry, and aa data sets had 87%, 87%, and 86% of the possible bipartitions resolved at ≥ 70% bs, respectively. all resulting trees have been deposited in the dryad data repository. bootstrap support values >50% are shown as percentages. when sister groups for the taxon of interest are listed, bootstrap support values on the left are for the clade including the taxon of interest and its sister group within viridiplantae, while bs values on the right are for the more inclusive clade excluding the taxon of interest. if only one bs value is given for a sister relationship, only two terminals are involved. the monophyly of chlorophyta receives 100% bs in all analyses. prasinophyceae are consistently not monophyletic. instead, the prasinophyte nephroselmis is sister to all other chlorophyta, while remaining prasinophyceae form a clade that is variously supported and is sister to a clade of the remaining chlorophyta. chlorophyceae are monophyletic, but trebouxiophyceae and ulvophyceae are not monophyletic, and the relationship of chlorophyceae to these lineages is unresolved. we consistently recovered a single set of relationships among the streptophytic algae subtending the land plant clade. zygnematophyceae are sister to land plants, coleochaetophyceae are sister to zygnematophyceae + embryophyta, charophyceae are sister to coleochaetophyceae +, and a clade of mesostigmatophyceae + chlorokybophyceae is sister to all other streptophyta. each of these relationships has ≥86% bs support. the branching order of the non-vascular land plant lineages differs among analyses. in analyses of the ntall and ry data sets, marchantiophyta, followed by bryophyta, and then anthocerotophyta are the earliest-branching land plant lineages, with anthocerotophyta the immediate sister to the vascular plants. in the ntall and ry analyses, these relationships had ≥89% bs support except for the bryophyta +  relationship in the ntall analysis, which received only 69% bs. in contrast, in the ntno3rd and aa analyses, bryophyta and marchantiophyta formed a clade, followed by anthocerophyta as sister to tracheophyta. within tracheophyta, the ntno3rd, ry, and aa data sets all place lycopodiophyta sister to a euphyllophyta clade. however, the analysis of the ntall data set places monilophyta sister to a clade of lycopodiophyta + spermatophyta. our analyses of monilophyta generally reveal strong support for a clade of equisetales + psilotales as sister to marattiales + leptosporangiate ferns. the lowest support obtained was for equisetales + psilotales in the ntno3rd analysis and ntall; all other nodes in all analyses received > 90% bs, with marattiales + leptosporangiate ferns receiving ≥ 99% bs. within spermatophyta, all analyses place the extant gymnosperms sister to angiospermae with 100% bs. within extant gymnosperms, cycadales and ginkgoales form a clade that is sister to a clade in which gnetophyta are nested within the paraphyletic conifers. there is generally high support placing gnetophyta as sister to a clade of araucariales + cupressales. this “gnecup” clade is then sister to pinales, which has 100% bs in all analyses. in all analyses, angiospermae receive 100% bs, and amborella is sister to all other angiosperms, followed by nymphaeales, and then austrobaileyales. these relationships are mostly supported by 100% bs. however, nymphaeales +  receives 81% bs in the ntno3rd analyses and 70% bs in the aa analyses. the remaining angiosperms receive 100% bs in all analyses. within mesangiospermae, the relationships among monocotyledoneae, magnoliidae, eudicotyledoneae, and ceratophyllum are not well supported and vary depending on the analysis. the strongest support for the placement of ceratophyllales is 75% bs as sister to eudicotyledoneae in the ry analysis. chloranthales receive 61-69% bs as sister to the well-supported magnoliidae. however, magnoliidae are not monophyletic in the aa analyses, where piperales are sister to ceratophyllales. within the monocot clade, acorales, followed by alismatales, have 100% bs in all analyses as subsequent sisters to the remaining monocots. in three of our analyses, a variously supported clade of liliales +  is sister to a clade of the remaining monocots. however, in the ry-coded analysis, pandanales + dioscoreales is sister to a clade of liliales +, which receives 69% bs. here asparagales + commelinidae is supported by 80% bs. within the eudicots, which receive 100% bs in all analyses, ranunculales are sister to the remaining taxa. in the ntall, ntno3rd, ry, and aa analyses, the clade of these remaining taxa receives 100%, 85%, 100%, and 62% bs, respectively. relationships vary among sabiaceae, proteales, and a clade of the remaining taxa, depending on the analysis. in the ntall and ntno3rd analyses, proteales + sabiaceae are supported as a clade, although with only 63% and 60% bs, respectively. however, in the ry analysis, proteales are sister to a clade containing sabiaceae plus the remaining taxa, which has 79% bs. in the aa analysis, relationships among these three clades are unresolved. among the remaining eudicots, we consistently recovered trochodendrales as sister to buxales + pentapetalae and gunnerales as sister to the remaining lineages of pentapetalae: dilleniaceae, superrosidae, and superasteridae. the placement of dilleniaceae remains uncertain. the family is sister to superrosidae in the ntall, ntno3rd, and ry analyses, but appears as sister to superasteridae in the aa analysis. within superrosidae, a clade of vitales + saxifragales is supported in the ntall, ntno3rd, and aa analyses. in the ry analysis, the relationship among saxifragales, vitales, and remaining rosidae is unresolved. fabidae and malvidae are both recovered with ≥ 99% bs in the ntall and ry analyses. however, each clade receives only 70% bs in the ntno3rd analysis. in the aa analysis neither clade is monophyletic; zygophyllales are embedded within a clade of malvidae taxa. the com clade is sister to a clade of fagales, cucurbitales, rosales, and fabales in fabidae in the aa, ry, and ntall trees and forms a trichotomy with zygophyllales and the clade of fagales, cucurbitales, rosales, and fabales in the ntno3rd tree. zygophyllales are sister to geraniales in the aa tree and sister to all other fabidae in the ntall and ry trees. superasteridae are recovered in all analyses. this clade receives 100% bs in the ntall and ry analyses, 95% bs in the ntno3rd analysis, and 66% bs in the aa analysis. santalales and berberidopsidales are strongly supported as subsequent sisters to caryophyllales + asteridae. within asteridae, cornales, followed by ericales, are subsequent sisters to a strongly supported clade that comprises strongly supported campanulidae and lamiidae clades. within lamiidae, the placement of boraginaceae is weak among the various analyses. boraginaceae are sister to gentianales in the aa tree, part of a trichotomy with lamiales and solanales + gentianales in the ntall tree, and sister to a weakly supported clade including gentianales, lamiales, and solanales in the ntno3rd and ry trees. analysis of only the third codon positions resulted in several very strong conflicts along the backbone of viridiplantae when compared to the topology from the ntno3rd analyses. these conflicts include the backbone relationships within chlorophyta, the placements of cycadales and lycopodiophyta, the relationships of the three major bryophyte lineages, and backbone relationships within poales. removal of four taxa with elevated rates of molecular evolution and few genes present in the data sets did not significantly affect the resulting topologies. discussion while the enormous phylogenetic data sets that result from new genome or transcriptome sequencing efforts can ameliorate the effects of random or stochastic error, they also may exacerbate the effects of systematic error, or error resulting from problems in the analysis, such as model inaccuracy. the high amount of agreement among our various analyses and strong support for results generally consistent with previous studies suggest that plastid genome sequence data hold much promise for resolving relationships throughout the green plants. however, several areas of conflict between analyses using different character-coding strategies demonstrate that plastid genome phylogenetics is also susceptible to systematic error. here we evaluate the phylogenetic results, emphasizing areas of agreement and concern, and then address some of the methodological issues raised by our results. evaluation of phylogenetic relationships historically, chlorophyta have been divided into prasinophyceae, trebouxiophyceae, chlorophyceae, and ulvophyceae based on the ultrastructure of the flagellar apparatus and features related to cytokinesis. the current status of green algae phylogenetics has been reviewed recently. the most comparable study to ours in terms of data and taxon sampling is by lang and nedelcu, who constructed a phylogeny of green algae with plastid genome sequence data. however, they analyzed only an amino acid data set using bayesian inference and the cat model. we found a paraphyletic prasinophyceae, which agrees with previous molecular analyses. however, lang and nedelcu recovered a monophyletic prasinophyceae, albeit with little support. chlorophyceae are monophyletic, which agrees with the results of lang and nedelcu. we also find that trebouxiophyceae and ulvophyceae are not monophyletic, and that the relationship of chlorophyceae to these lineages is unresolved. the branching order of the various trebouxiophyceae, ulvophyceae, and chlorophyceae lineages within chlorophyta, unresolved in our analyses, was also uncertain in earlier analyses. similarly in lang and nedelcu, trebouxiophyceae and ulvophyceae were not supported as monophyletic, although unlike our results, almost all nodes in their phylogeny were maximally supported. our analyses provide consistent, strong support for the relationships of streptophytic algae to land plants, and all analyses support zygnematophyceae as the sister to land plants. relationships among these lineages and the closest relatives of land plants have varied in previous studies depending on taxon sampling and gene choice. some studies agree with our results placing zygnematophyceae as sister to land plants, while other phylogenetic analyses indicate that charophyceae or coleochaetophyceae occupy this position. depending on the analysis, zhong et al. found either zygnematophyceae alone or a clade of zygnematophyceae + coleochaetophyceae as sister to land plants. in particular, the results of lang and nedelcu conflict with our results regarding the sister group to embryophyta. while we find a clade of coleochaetophyceae +, their results strongly support zygnematophyceae +. phylogenetic relationships among bryophytes are also contentious, and nearly every possible relationship among these lineages has been reported, often with strong support. most studies have shown the bryophytes as paraphyletic with respect to tracheophyta rather than as a clade. as recovered in our ntall and ry analyses, liverworts often are placed sister to all other land plants, followed by mosses, and with hornworts sister to tracheophyta. a sister relationship between mosses and liverworts, found in our ntno3rd and aa analyses, was proposed previously based on morphological and molecular data and has been recovered with numerous nuclear genes. this relationship was also recovered in analyses of complete plastid genome data by karol et al. when divergent taxa were excluded from phylogenetic analyses and also by wolf and karol when third positions were excluded. our results placing lycopodiophyta sister to euphyllophyta in all but the ntall analysis agree with most molecular phylogenetic analyses. this split is also supported by analyses of morphological characters in fossil and extant taxa. monilophyta and spermatophyta also possess a 30-kb inversion in the large single-copy region of the plastid genome not found in lycopodiophyta and the three bryophyte clades. in the ntall analysis, euphyllophyta are not monophyletic; lycopodiophyta, rather than monilophyta, are sister to spermatophyta. this relationship has been reported previously; however, it likely is a phylogenetic artifact, perhaps related to base composition bias. the plastid genome of the lycophyte selaginella has an especially high gc content, with selaginella unicata having the highest gc content in our ntall data set. in some previous studies, relationships among lineages of monilophyta have not been well resolved or supported. as a result, the relationships among equisetales, psilotales, marattiales, and leptosporangiate ferns are often represented as a polytomy. in contrast, most of our analyses recovered strong support for a clade of equisetales + psilotales as sister to marattiales + leptosporangiate ferns. these relationships agree with recent studies of monilophyte relationships based on plastid genome sequence data, although support is stronger here. unfortunately, ophioglossales, which often appear as sister to psilotales, lacked a sequenced plastome at the time of our analyses. however, plastid genome data for ophioglossales have subsequently been published and analyzed in a phylogenetic context, with strong support for ophioglossales as sister to psilotales and weak support for this clade as sister to equisetales. results from that study with regard to marattiales and leptosporangiate ferns agree with the relationships presented here. relationships among the lineages of extant seed plants, and especially the placement of gnetophyta, have long been debated. gnecup trees, found in all of our analyses, were initially recovered by nickrent et al., and then more recently by zhong et al.. however, zhong et al. suggested that the support for gnecup may be the result of long-branch attraction; by removing highly variable proteins, support for gnecup decreased. furthermore, by removing what they considered parallel substitutions between lineages leading to gnetophyta and to cryptomeria, a gnepine topology was recovered. although several different placements for gnetophyta have been recovered and strongly supported, many studies involving multiple genes have placed gnetophyta sister to pinales. using both coalescent and concatenation analyses, xi et al. found that the phylogenetic placement of gnetophyta differs between the nuclear and plastid genomes. in their analyses using nuclear data, the gnepine hypothesis is supported, while their analyses of plastid data support the gnecup hypothesis. in contrast, lee et al. found strong support for gnetophyta sister to the remaining gymnosperms in an ml analysis of sets of nuclear gene orthologs from land plant genera. the backbone relationships among angiosperm lineages generally agree with results from recent analyses, including a 17-gene analysis of angiosperms and previous analyses of plastid genome data sets. the position of ceratophyllum, and thus the relationships among monocotyledoneae, eudicotyledoneae, and magnoliidae, varies among our analyses, although without strong support. this contrasts with several other large, multi-gene analyses in which monocotyledoneae are sister to ceratophyllales + eudicotyledoneae. interestingly, the strongest support for the placement of ceratophyllales sister to eudicotyledoneae is in the ry analysis. however, in that analysis, the relationships among ceratophyllales + eudicotyledoneae, monocotyledoneae, and magnoliidae are unresolved. within the angiosperms, some relationships that have been uncertain, particularly at deep levels, receive moderate to strong support in at least some of our analyses. for example, the placement of myrtales and geraniales in the malvidae is supported with 70% bs in the ntno3rd tree and ≥ 99% bs in the ry and ntall trees. myrtales and geraniales are also placed in a clade with the malvidae taxa in the aa analysis; however, zygophyllales are also included within this clade, making malvidae non-monophyletic. likewise, chloranthales are sister to magnoliidae in all trees, but with weaker support. in two cases, all analyses but ry resolve relationships, with ry producing a polytomy that does not conflict with the resolutions found in the other analyses. these two cases are as follows: vitales + saxifragales supported by ≥ 70% bs in all analyses but ry, with saxifragales, vitales, and remaining rosidae forming a polytomy in the ry tree; dasypogonaceae + arecales in all but ry and a trichotomy of dasypogonaceae, arecales, and poales +  in the ry tree. in two additional cases when ry is compared to the other three analyses, the ry analysis produced either stronger support for the placement of a taxon or a different placement altogether. first, in the ntall, ntno3rd, and aa analyses, the position of sabiaceace among the early-diverging lineages of eudicotyledoneae is weakly supported. however, in the ry analysis, sabiaceae receive moderate support as sister to a strongly supported clade of trochodendrales + ). this contrasts with previous studies that often place sabiaceae as sister to proteales. an example of a different placement of a taxon in the ry analysis when compared to the other analyses involves liliales. the ntall, ntno3rd, and the aa analyses support liliales as sister to a clade of dioscoreales + pandanales with 72%, 69%, and 80% bs, respectively. this placement of liliales was also recovered in barrett et al.. in contrast, in the ry analysis, liliales are placed in a clade with asparagales + commelinidae with moderate support. this latter placement of liliales was strongly supported in an analysis with much better taxon sampling. some taxa that have been problematic in previous studies continue to defy definitive placement. their positions vary among our analyses, although they are generally not well supported in some, or all, of the trees. despite its general placement of the com clade in fabidae in these and other plastid analyses, this clade is more closely related to malvidae in some analyses, particularly those using mitochondrial gene sequences. recent analyses of plastid, mitochondrial, and nuclear data suggest that the com clade may represent ancient reticulation involving fabidae and malvidae during the rapid radiation of rosidae. methodological issues of plastid phylogenomic analyses to address potential systematic error in large-scale phylogenetic analyses, scientists often either try to improve the fit of models to the data or change or remove problematic data. with increasing sequence length and number of genes, it is more likely that a sequence alignment will contain regions with heterogeneous processes of molecular evolution. we see evidence of this high heterogeneity with our model-fitting experiments, which always favor the most parameter-rich models. thus, defining partitioning schemes and models that can accurately reflect the true processes of molecular evolution while not over-parameterizing the analysis remains critically important for phylogenetic analyses of large plastid data sets. although we assessed models that account for heterogeneity in patterns of molecular evolution among genes and in some cases codon positions, our model selection tests only evaluated a small selection of possible models and partitioning schemes. it is possible that other partitioning schemes could enable simpler models. most conventional phylogenetic models, like those used in our analyses, also assume homogeneous processes of evolution throughout the tree. yet when the branches of the phylogeny encompass over one billion years of evolutionary history, as likely do those in the green plants, the patterns of evolution almost certainly differ among lineages and through time. this is apparent from the often good fit of covarion models to plastid genes and the presence of nucleotide compositional heterogeneity, which can confound conventional phylogenetic analyses. also, our models do not account for shifts in selective pressure or instances of positive selection that will affect nucleotide and amino acid substitution patterns. nucleotide compositional heterogeneity remains a concern for green plant plastid genome analyses. this variation is most evident in non-seed plant taxa, and thus it has not been a focus of many previous phylogenetic analyses of plastid genome sequences. a gc bias in itself is not necessarily problematic for phylogenetic analyses, but nearly all commonly used models for likelihood-based phylogenetic analyses assume single equilibrium nucleotide frequencies. given that gc content appears to vary by codon position in plants, a partitioning scheme that estimates separate nucleotide frequencies for each codon position may account for some of the spatial heterogeneity in gc content in the plastid genome, but it does not address the differences in gc frequency among lineages. a commonly used strategy to reduce the effects of gc heterogeneity across lineages is ry-coding, in which the purines are coded as rs and the pyrimidines are coded as ys. ry-coding can reduce the compositional variability among lineages, improve the fit of models, and increase the signal for internal branches. an obvious disadvantage to ry-coding is that by coding the sequences with two character states instead of four, it reduces the amount of information in the sequences. in general, we see little overall reduction, and even some gains, in bootstrap support when using ry-coding compared to the use of all nucleotide data, suggesting that the benefits of ry-coding make up for any potential costs of information loss. perhaps the biggest topological difference in the ry phylogeny compared to ntall is the placement of monilophyta rather than lycopodiophyta as sister to seed plants. the unexpected placement of lycopodiophyta as the sister to seed plants in the ntall analysis is almost certainly an artifact of systematic error; several other lines of evidence support monilophyta as the sister group of seed plants. approaches to reducing systematic errors by excluding problematic data, which often include fast-evolving or saturated sites, also have been suggested for plastid genome analyses. with the proper model of molecular evolution and adequate taxon sampling, fast sites are not necessarily problematic; they are only problematic insofar as they are difficult to model. yet with heterogeneous processes of molecular evolution throughout the tree, the fast-evolving or saturated sites can produce a significant non-phylogenetic signal. indeed, the third codon positions appear to have especially high levels of compositional heterogeneity, potentially causing systematic error, and an analysis of just the third codon positions conflicts with the analyses of other data sets in several critical parts of the tree. however, third codon positions also represent a large proportion of the variable sites in the alignment, and removing them may exclude much of the phylogenetic information in some parts of the tree. with regard to backbone relationships in our phylogeny, excluding the third position sites produces several interesting changes in contrast to ntall: 1) it supports the sister relationship of mosses and liverworts, 2) monilophytes, not lycophytes, are placed sister to seed plants as expected, and 3) support for some of the backbone angiosperm relationships is reduced. thus, the effects of removing the third codon position sites appear to vary in different parts of the tree. another strategy for overcoming potential error associated with fast-evolving sites is to code the sequences as amino acids rather than nucleotides. this does not necessarily eliminate problems of compositional heterogeneity, as the gc bias also may bias amino acid composition. regarding backbone green plant relationships, the aa analysis provided similar results to analyses of only first and second codon positions. aa analysis also produced some weakly supported, questionable relationships among angiosperm lineages. in previous deep-level plant analyses, analyses of amino acid data have resulted in arguably more problematic or questionable relationships than analyses of nucleotide data. however, these results are likely due to inappropriate models of amino acid evolution, and with better models, optimized for plastid evolution, amino acid data may be a valuable source of phylogenetic information. taxon sampling is also important for plastid phylogenomic studies, especially when the model of evolution is inadequate, and genome-scale analyses often have limited taxon sampling. new methods for rapid and inexpensive plastid genome sequencing may ameliorate the effects of insufficient sampling of extant taxa; however, many major lineages of green plants are now extinct, precluding their inclusion in analyses of molecular data. in addition, ancient, rapid radiations abound within portions of the green plant tree of life, creating extremely difficult phylogenetic problems no matter the taxon sampling. furthermore, even in the absence of systematic error, it is possible that a tree built from plastid genome data will not reflect species relationships. the plastid genome represents a single locus of linked genes. for phylogenetic analyses, this can be beneficial because combining genes with different evolutionary histories into a single character matrix can lead to phylogenetic error. however, incomplete lineage sorting or ancient reticulation could lead to conflict between the plastid gene tree and the species phylogeny. for this reason, it will be interesting to compare phylogenetic hypotheses from the plastid genome with independent phylogenetic estimates from numerous nuclear and mitochondrial loci. finally, while full plastid genome sequence data provide much power for resolving difficult phylogenetic relationships, it is not clear that they can resolve all plant relationships. theoretical work suggests that extremely large data sets may be necessary to resolve some relationships when the internal nodes are separated by very short branches, and recent analyses indicate that full plastid genomes are not sufficient to reject alternative topologies among monocots. indeed, the unresolved or conflicting parts of the green plant phylogeny in our analyses are generally associated with short internal branch lengths. thus, even if the model of evolution accurately reflects the true process of molecular evolution, and there is no systematic error, plastid genome data alone may not be sufficient to resolve all parts of the green plant tree of life. that is, the topology may not be identifiable with the plastid data alone. a recent analysis using a new diagnostic test for phylogenetic identifiability based on data cloning suggested that a backbone topology of angiosperms was identifiable from plastid sequence data using the gtr + Γ model, but the tree in this paper is much larger and the models more complex. in any case, it will be necessary to include perspectives from the nuclear genome and phenotypic data before we are confident about all deep-level relationships among green plants.  our diverse analyses provide a first approach to addressing some of the difficult issues associated with plastid phylogenetic analyses at this evolutionary depth and level of taxon sampling. the results of the analyses using different models, character-coding strategies, and character subsets suggest that much of the tree is robust to many different phylogenetic approaches, and they highlight regions of the tree that need more scrutiny. more sophisticated modelling approaches may more accurately characterize the heterogeneous processes of molecular evolution, but it is also crucial that the parameters of these complex models can be estimated by the data at hand. while it may be impossible for any model to reflect perfectly the complexities of molecular evolution, as we better characterize these processes it will be possible to examine through simulations their possible effects on phylogenetic analyses and to recognize phylogenetic error caused by model misspecification. 
we assembled protein-coding sequence data for genes from diverse green plant taxa with complete or nearly complete plastid genome sequences available from genbank. 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. next-generation sequencing has provided a wealth of plastid genome sequence data from an increasingly diverse set of green plants. here we explore the limits and challenges of inferring a comprehensive green plant phylogeny from available complete or nearly complete plastid genome sequence data. analyses of the plastid sequence data recovered a strongly supported framework of relationships for green plants. 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. although these data have helped resolve the phylogeny of numerous clades, their utility for inferring relationships across all green plants is uncertain. 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. 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. this framework includes: i) the placement of zygnematophyceace as sister to land plants, ii) a clade of extant gymnosperms with cycads + ginkgo sister to remaining extant gymnosperms and with gnetophytes sister to non-pinaceae conifers, and iii) within the monilophyte clade, equisetales + psilotales are sister to marattiales + leptosporangiate ferns.
 next-generation sequencing has provided a wealth of plastid genome sequence data from an increasingly diverse set of green plants. although these data have helped resolve the phylogeny of numerous clades, their utility for inferring relationships across all green plants is uncertain. viridiplantae originated 700- million years ago and may comprise as many as 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.  we assembled protein-coding sequence data for genes from 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.  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, ii) a clade of extant gymnosperms with cycads + ginkgo sister to remaining extant gymnosperms and with gnetophytes sister to non-pinaceae conifers, and iii) within the monilophyte clade, 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. composition biasphylogenomicsplastid genome sequencesplastomesry-codingviridiplantae 