Assessment of Absolute Substitution Model Fit Accommodating Time-Reversible and Non-Time-Reversible Evolutionary Processes.

The loss of information accompanying assessment of absolute fit of substitution models to phylogenetic data negatively affects the discriminatory power of previous methods and can make them insensitive to lineage-specific changes in the substitution process. As an alternative, I propose evaluating absolute fit of substitution models based on a novel statistic which describes the observed data without information loss and which is unlikely to become zero-inflated with increasing numbers of taxa. This method can accommodate gaps and is sensitive to lineage-specific shifts in the substitution process. In simulation experiments, it exhibits greater discriminatory power than previous methods. The method can be implemented in both Bayesian and Maximum Likelihood phylogenetic analyses, and used to screen any set of models. Recently, it has been suggested that model selection may be an unnecessary step in phylogenetic inference. However, results presented here emphasize the importance of model fit assessment for reliable phylogenetic inference. [Absolute model fit; model misspecification; origin of plastids; phylogenomics.].

The grapevine VvibZIPC22 transcription factor is involved in the regulation of flavonoid biosynthesis.

In grapevine, flavonoids constitute one of the most abundant subgroups of secondary metabolites, influencing the quality, health value, and typicity of wines. Their synthesis in many plant species is mainly regulated at the transcriptional level by modulation of flavonoid pathway genes either by single regulators or by complexes of different regulators. In particular, bZIP and MYB factors interact synergistically in the recognition of light response units present in the promoter of some genes of the pathway, thus mediating light-dependent flavonoid biosynthesis. We recently identified VvibZIPC22, a member of clade C of the grapevine bZIP family, in a quantitative trait locus (QTL) specifically associated with kaemperol content in mature berries. Here, to validate the involvement of this candidate gene in the fine regulation of flavonol biosynthesis, we characterized its function by in vitro and in vivo experiments. A role for this gene in the control of flavonol biosynthesis was indeed confirmed by its highest expression at flowering and during UV light-mediated induction, paralleled by accumulation of the flavonol synthase 1 transcript and flavonol compounds. The overexpression of VvibZIPC22 in tobacco caused a significant increase in several flavonoids in the flower, via induction of general and specific genes of the pathway. In agreement with this evidence, VvibZIPC22 was able to activate the promoters of specific genes of the flavonoid pathway, alone or together with other factors, as revealed by transient reporter assays. These findings, supported by in silico indications, allowed us to propose VvibZIPC22 as a new regulator of flavonoid biosynthesis in grapevine.

Streptophyte algae and the origin of land plants revisited using heterogeneous models with three new algal chloroplast genomes.

The phylogenetic branching order of the green algal groups that gave rise to land plants remains uncertain despite its fundamental importance to understanding plant evolution. Previous studies have demonstrated that land plants evolved from streptophyte algae, but different lineages of streptophytes have been suggested to be the sister group of land plants. To better understand the evolutionary history of land plants and to determine the potential effects of "long-branch attraction" in phylogenetic reconstruction, we analyzed a chloroplast genome data set including three new chloroplast genomes from streptophyte algae: Coleochaetae orbicularis (Coleochaetales), Nitella hookeri (Charales), and Spirogyra communis (Zygnematales). We further applied a site pattern sorting method together with site- and time-heterogeneous models to investigate the branching order among streptophytes and land plants. Our chloroplast phylogenomic analyses support previous hypotheses based on nuclear data in placing Zygnematales alone, or a clade consisting of Coleochaetales plus Zygnematales, as the closest living relatives of land plants.

Phylogenetic analysis of 47 chloroplast genomes clarifies the contribution of wild species to the domesticated apple maternal line.

Both the origin of domesticated apple and the overall phylogeny of the genus Malus are still not completely resolved. Having this as a target, we built a 134,553-position-long alignment including two previously published chloroplast DNAs (cpDNAs) and 45 de novo sequenced, fully colinear chloroplast genomes from cultivated apple varieties and wild apple species. The data produced are free from compositional heterogeneity and from substitutional saturation, which can adversely affect phylogeny reconstruction. Phylogenetic analyses based on this alignment recovered a branch, having the maximum bootstrap support, subtending a large group of the cultivated apple sorts together with all analyzed European wild apple (Malus sylvestris) accessions. One apple cultivar was embedded in a monophylum comprising wild M. sieversii accessions and other Asian apple species. The data demonstrate that M. sylvestris has contributed chloroplast genome to a substantial fraction of domesticated apple varieties, supporting the conclusion that different wild species should have contributed the organelle and nuclear genomes to the domesticated apple.

The evolutionary root of flowering plants.

Correct rooting of the angiosperm radiation is both challenging and necessary for understanding the origins and evolution of physiological and phenotypic traits in flowering plants. The problem is known to be difficult due to the large genetic distance separating flowering plants from other seed plants and the sparse taxon sampling among basal angiosperms. Here, we provide further evidence for concern over substitution model misspecification in analyses of chloroplast DNA sequences. We show that support for Amborella as the sole representative of the most basal angiosperm lineage is founded on sequence site patterns poorly described by time-reversible substitution models. Improving the fit between sequence data and substitution model identifies Trithuria, Nymphaeaceae, and Amborella as surviving relatives of the most basal lineage of flowering plants. This finding indicates that aquatic and herbaceous species dominate the earliest extant lineage of flowering plants. [; ; ; ; ; .].

The mitochondrial genome of Malus domestica and the import-driven hypothesis of mitochondrial genome expansion in seed plants.

Mitochondrial genomes of spermatophytes are the largest of all organellar genomes. Their large size has been attributed to various factors; however, the relative contribution of these factors to mitochondrial DNA (mtDNA) expansion remains undetermined. We estimated their relative contribution in Malus domestica (apple). The mitochondrial genome of apple has a size of 396 947 bp and a one to nine ratio of coding to non-coding DNA, close to the corresponding average values for angiosperms. We determined that 71.5% of the apple mtDNA sequence was highly similar to sequences of its nuclear DNA. Using nuclear gene exons, nuclear transposable elements and chloroplast DNA as markers of promiscuous DNA content in mtDNA, we estimated that approximately 20% of the apple mtDNA consisted of DNA sequences imported from other cell compartments, mostly from the nucleus. Similar marker-based estimates of promiscuous DNA content in the mitochondrial genomes of other species ranged between 21.2 and 25.3% of the total mtDNA length for grape, between 23.1 and 38.6% for rice, and between 47.1 and 78.4% for maize. All these estimates are conservative, because they underestimate the import of non-functional DNA. We propose that the import of promiscuous DNA is a core mechanism for mtDNA size expansion in seed plants. In apple, maize and grape this mechanism contributed far more to genome expansion than did homologous recombination. In rice the estimated contribution of both mechanisms was found to be similar.

A genome-wide phylogenetic reconstruction of family 1 UDP-glycosyltransferases revealed the expansion of the family during the adaptation of plants to life on land.

For almost a decade, our knowledge on the organisation of the family 1 UDP-glycosyltransferases (UGTs) has been limited to the model plant A. thaliana. The availability of other plant genomes represents an opportunity to obtain a broader view of the family in terms of evolution and organisation. Family 1 UGTs are known to glycosylate several classes of plant secondary metabolites. A phylogeny reconstruction study was performed to get an insight into the evolution of this multigene family during the adaptation of plants to life on land. The organisation of the UGTs in the different organisms was also investigated. More than 1500 putative UGTs were identified in 12 fully sequenced and assembled plant genomes based on the highly conserved PSPG motif. Analyses by maximum likelihood (ML) method were performed to reconstruct the phylogenetic relationships existing between the sequences. The results of this study clearly show that the UGT family expanded during the transition from algae to vascular plants and that in higher plants the clustering of UGTs into phylogenetic groups appears to be conserved, although gene loss and gene gain events seem to have occurred in certain lineages. Interestingly, two new phylogenetic groups, named O and P, that are not present in A. thaliana were discovered.

Mitochondrial DNA of Vitis vinifera and the issue of rampant horizontal gene transfer.

The mitochondrial genome of grape (Vitis vinifera), the largest organelle genome sequenced so far, is presented. The genome is 773,279 nt long and has the highest coding capacity among known angiosperm mitochondrial DNAs (mtDNAs). The proportion of promiscuous DNA of plastid origin in the genome is also the largest ever reported for an angiosperm mtDNA, both in absolute and relative terms. In all, 42.4% of chloroplast genome of Vitis has been incorporated into its mitochondrial genome. In order to test if horizontal gene transfer (HGT) has also contributed to the gene content of the grape mtDNA, we built phylogenetic trees with the coding sequences of mitochondrial genes of grape and their homologs from plant mitochondrial genomes. Many incongruent gene tree topologies were obtained. However, the extent of incongruence between these gene trees is not significantly greater than that observed among optimal trees for chloroplast genes, the common ancestry of which has never been in doubt. In both cases, we attribute this incongruence to artifacts of tree reconstruction, insufficient numbers of characters, and gene paralogy. This finding leads us to question the recent phylogenetic interpretation of Bergthorsson et al. (2003, 2004) and Richardson and Palmer (2007) that rampant HGT into the mtDNA of Amborella best explains phylogenetic incongruence between mitochondrial gene trees for angiosperms. The only evidence for HGT into the Vitis mtDNA found involves fragments of two coding sequences stemming from two closteroviruses that cause the leaf roll disease of this plant. We also report that analysis of sequences shared by both chloroplast and mitochondrial genomes provides evidence for a previously unknown gene transfer route from the mitochondrion to the chloroplast.

Automated removal of noisy data in phylogenomic analyses.

Noisy data, especially in combination with misalignment and model misspecification can have an adverse effect on phylogeny reconstruction; however, effective methods to identify such data are few. One particularly important class of noisy data is saturated positions. To avoid potential errors related to saturation in phylogenomic analyses, we present an automated procedure involving the step-wise removal of the most variable positions in a given data set coupled with a stopping criterion derived from correlation analyses of pairwise ML distances calculated from the deleted (saturated) and the remaining (conserved) subsets of the alignment. Through a comparison with existing methods, we demonstrate both the effectiveness of our proposed procedure for identifying noisy data and the effect of the removal of such data using a well-publicized case study involving placental mammals. At the least, our procedure will identify data sets requiring greater data exploration, and we recommend its use to investigate the effect on phylogenetic analyses of removing subsets of variable positions exhibiting weak or no correlation to the rest of the alignment. However, we would argue that this procedure, by identifying and removing noisy data, facilitates the construction of more accurate phylogenies by, for example, ameliorating potential long-branch attraction artefacts.

Phylogenetic analyses with systematic taxon sampling show that mitochondria branch within Alphaproteobacteria.

Though it is well accepted that mitochondria originated from an alphaproteobacteria-like ancestor, the phylogenetic relationship of the mitochondrial endosymbiont to extant Alphaproteobacteria is yet unresolved. The focus of much debate is whether the affinity between mitochondria and fast-evolving alphaproteobacterial lineages reflects true homology or artefacts. Approaches such as site exclusion have been claimed to mitigate compositional heterogeneity between taxa, but this comes at the cost of information loss, and the reliability of such methods is so far unproven. Here we demonstrate that site-exclusion methods produce erratic phylogenetic estimates of mitochondrial origin. Thus, previous phylogenetic hypotheses on the origin of mitochondria based on pretreated datasets should be re-evaluated. We applied alternative strategies to reduce phylogenetic noise by systematic taxon sampling while keeping site substitution information intact. Cross-validation based on a series of trees placed mitochondria robustly within Alphaproteobacteria, sharing an ancient common ancestor with Rickettsiales and currently unclassified marine lineages.

Removal of noisy characters from chloroplast genome-scale data suggests revision of phylogenetic placements of Amborella and Ceratophyllum.

It is widely appreciated that noisy, highly variable data can impede phylogeney reconstruction. Researchers have for a long time omitted problematic data from phylogenetic analyses, such as the third-codon positions and variable regions. In the analyses of the phylogenetic relations of the angiosperms; however, inclusion of complete gene sequences into genomic-scale alignments has become a common practice. Here we demonstrate that this practice can be misleading. We show that support of the basal-most position of Amborella trichopoda among the angiosperms in the chloroplast genomic data is based only on a tiny subset (< 1% of the total alignment length) of the most variable positions in alignment, exhibiting mean maximum likelihood (ML) distance among the angiosperm operational taxonomic units (OTUs) approximately 36 substitutions/site. Exclusion of these positions leads to disappearance of the basal Amborella branch. Likewise, the recently reported sister-group relationship of Ceratophyllum to the eudicots is based on the presence of 2% of the most variable positions in the genomic alignment, exhibiting, on average, 20 substitutions/site in comparison among the angiosperm OTUs. These observations highlight a need for excluding a certain proportion of saturated positions in alignment from phylogenomic analyses.

A Novel Test for Absolute Fit of Evolutionary Models Provides a Means to Correctly Identify the Substitution Model and the Model Tree.

A novel test is described that visualizes the absolute model-data fit of the substitution and tree components of an evolutionary model. The test utilizes statistics based on counts of character state matches and mismatches in alignments of observed and simulated sequences. This comparison is used to assess model-data fit. In simulations conducted to evaluate the performance of the test, the test estimator was able to identify both the correct tree topology and substitution model under conditions where the Goldman-Cox test-which tests the fit of a substitution model to sequence data and is also based on comparing simulated replicates with observed data-showed high error rates. The novel test was found to identify the correct tree topology within a wide range of DNA substitution model misspecifications, indicating the high discriminatory power of the test. Use of this test provides a practical approach for assessing absolute model-data fit when testing phylogenetic hypotheses.

The Rpv3-3 Haplotype and Stilbenoid Induction Mediate Downy Mildew Resistance in a Grapevine Interspecific Population.

The development of new resistant varieties to the oomycete Plasmopara viticola (Berk.& Curt) is a promising way to combat downy mildew (DM), one of the major diseases threatening the cultivated grapevine (Vitis vinifera L.). Taking advantage of a segregating population derived from "Merzling" (a mid-resistant hybrid) and "Teroldego" (a susceptible landrace), 136 F1 individuals were characterized by combining genetic, phenotypic, and gene expression data to elucidate the genetic basis of DM resistance and polyphenol biosynthesis upon P. viticola infection. An improved consensus linkage map was obtained by scoring 192 microsatellite markers. The progeny were screened for DM resistance and production of 42 polyphenols. QTL mapping showed that DM resistance is associated with the herein named Rpv3-3 specific haplotype and it identified 46 novel metabolic QTLs linked to 30 phenolics-related parameters. A list of the 95 most relevant candidate genes was generated by specifically exploring the stilbenoid-associated QTLs. Expression analysis of 11 genes in Rpv3-3 +/- genotypes displaying disparity in DM resistance level and stilbenoid accumulation revealed significant new candidates for the genetic control of stilbenoid biosynthesis and oligomerization. These overall findings emphasized that DM resistance is likely mediated by the major Rpv3-3 haplotype and stilbenoid induction.

Chloroplast genome phylogenetics: why we need independent approaches to plant molecular evolution.

The traditional approach to plant molecular phylogenetics involves amplifying, sequencing and analyzing one or a few genes from many species and is conducive to broad taxon sampling. An independent approach involves chloroplast genome sequencing, providing much larger amounts of data per taxon but for a smaller number of species. In principle, the two strategies can inform each other but in practice their results sometimes conflict for reasons that are currently debated. An Opinion article published in the October 2004 issue of Trends in Plant Science cautioned against the pursuit of genome-based phylogenies. Here, we provide a different perspective on issues at the heart of the current debate and defend the use of chloroplast genome phylogenetics for crucial species because it provides an independent test of hypotheses generated by the traditional approach.

A new test of phylogenetic model fitness addresses the issue of the basal angiosperm phylogeny.

We readdress the issue of phylogeny of the basal extant angiosperms employing a source previously not systematically investigated, specifically, the non-coding sequences of cpDNA. Comparison of trees with and without grasses or the outgroup (Pinus) in our analyses revealed no rearrangements in tree topology that might be expected if LBA were distorting the position of the magnoliids. For each model applied, irrespective of whether monocots or ANITA members appeared basally divergent, the orchid Phalaenopsis assumed the same position on the trees with the reduced taxon set as did the branch bearing the orchid plus the grasses in the full alignment. However, our new test of model fitness revealed a different flaw influencing the placement of monocots, which is related to model mis-specification. This flaw similarly affects the full alignment and the alignment with grasses removed. In both cases the models favouring a relatively derived position for the monocots and basal placement of the branch of Amborella plus Nymphaea provide better overall prediction of the observed data structure. In the view of apparent unsuitability of the bootstrap method for large data sets, our novel test provides a new means of exploring conflicts caused by systematic errors in phylogenetic analyses.

Genomes of Stigonematalean cyanobacteria (subsection V) and the evolution of oxygenic photosynthesis from prokaryotes to plastids.

Cyanobacteria forged two major evolutionary transitions with the invention of oxygenic photosynthesis and the bestowal of photosynthetic lifestyle upon eukaryotes through endosymbiosis. Information germane to understanding those transitions is imprinted in cyanobacterial genomes, but deciphering it is complicated by lateral gene transfer (LGT). Here, we report genome sequences for the morphologically most complex true-branching cyanobacteria, and for Scytonema hofmanni PCC 7110, which with 12,356 proteins is the most gene-rich prokaryote currently known. We investigated components of cyanobacterial evolution that have been vertically inherited, horizontally transferred, and donated to eukaryotes at plastid origin. The vertical component indicates a freshwater origin for water-splitting photosynthesis. Networks of the horizontal component reveal that 60% of cyanobacterial gene families have been affected by LGT. Plant nuclear genes acquired from cyanobacteria define a lower bound frequency of 611 multigene families that, in turn, specify diazotrophic cyanobacterial lineages as having a gene collection most similar to that possessed by the plastid ancestor.

Systematic error in seed plant phylogenomics.

Resolving the closest relatives of Gnetales has been an enigmatic problem in seed plant phylogeny. The problem is known to be difficult because of the extent of divergence between this diverse group of gymnosperms and their closest phylogenetic relatives. Here, we investigate the evolutionary properties of conifer chloroplast DNA sequences. To improve taxon sampling of Cupressophyta (non-Pinaceae conifers), we report sequences from three new chloroplast (cp) genomes of Southern Hemisphere conifers. We have applied a site pattern sorting criterion to study compositional heterogeneity, heterotachy, and the fit of conifer chloroplast genome sequences to a general time reversible + G substitution model. We show that non-time reversible properties of aligned sequence positions in the chloroplast genomes of Gnetales mislead phylogenetic reconstruction of these seed plants. When 2,250 of the most varied sites in our concatenated alignment are excluded, phylogenetic analyses favor a close evolutionary relationship between the Gnetales and Pinaceae-the Gnepine hypothesis. Our analytical protocol provides a useful approach for evaluating the robustness of phylogenomic inferences. Our findings highlight the importance of goodness of fit between substitution model and data for understanding seed plant phylogeny.

The genome of the domesticated apple (Malus × domestica Borkh.).

We report a high-quality draft genome sequence of the domesticated apple (Malus × domestica). We show that a relatively recent (>50 million years ago) genome-wide duplication (GWD) has resulted in the transition from nine ancestral chromosomes to 17 chromosomes in the Pyreae. Traces of older GWDs partly support the monophyly of the ancestral paleohexaploidy of eudicots. Phylogenetic reconstruction of Pyreae and the genus Malus, relative to major Rosaceae taxa, identified the progenitor of the cultivated apple as M. sieversii. Expansion of gene families reported to be involved in fruit development may explain formation of the pome, a Pyreae-specific false fruit that develops by proliferation of the basal part of the sepals, the receptacle. In apple, a subclade of MADS-box genes, normally involved in flower and fruit development, is expanded to include 15 members, as are other gene families involved in Rosaceae-specific metabolism, such as transport and assimilation of sorbitol.

The molecular phylogeny of Rebutia (Cactaceae) and its allies demonstrates the influence of paleogeography on the evolution of South American mountain cacti.

The tropical Andes harbor a major part of the world's plant biodiversity. The montane cacti of the tribes Browningieae, Cereeae, and Trichocereeae underwent extensive radiation and thus are well suited as a model group to study the diversification of Andean plants. We reconstructed their phylogeny employing three noncoding chloroplast regions and explained it in the context of the geological history of South America. We found that the clade of cephalia-bearing cacti with naked pericarpels is centered in northeastern Brazil, whereas almost all other clades comprise Andean species. The spatial split between the clades was probably caused by the Andean uplift and the concurrent formation of intracontinental marine basins in the Tertiary. The phylogenetic reconstructions based on parsimony and Bayesian approaches do not reflect the traditional delimitation of the tribes and of the large genera. Our results suggest that Rebutia s.l. and Echinopsis s.l. are not monophyletic and that Sulcorebutia, Weingartia, and Cintia should be united into one genus. Even though this "Weingartia-complex" and the genus Gymnocalycium are similar in size and morphological diversity, Gymnocalycium has a very high molecular divergence suggesting a comparably older radiation.