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

The performance of outgroup-free rooting under evolutionary radiations

Tree rooting implies a temporal dimension to phylogenies. Only after defining the position of the root node is that the ancestral-descendant relationship between branches can be fully deduced. Rooting has been usually carried out by employing evolutionarily close outgroup lineages, which is a drawback when these lineages are unavailable or unknown. Alternatively, outgroup-free rooting methods were proposed, which rely on the constancy of evolutionary rates to varying degrees. In this work we analyzed the performance of two of these methods, the midpoint rooting (MPR) and the minimal ancestor deviation (MAD), in rooting topologies evolved under challenging scenarios of fast evolutionary radiations derived from empirical data, characterized by short internal branches near the crown node. Considering all branch length combinations investigated, both methods exhibited average success rates below 50%, although MAD slightly outperformed MPR. Moreover, tree balance significantly impacted the relative performance of the methods. We found that, in four-taxa unrooted trees, the outcome of whether both methodologies will correctly root the tree can be roughly predicted by two simple dimensionless metrics: the coefficient of variation of the external branch lengths, and the ratio between the internal branch length to the total sum of branch lengths, which were employed to devise a general linear model that allowed calculating the probability of correct placing the root node for any four-taxa tree. We predicted that the performance of both outgroup-free rooting methods on loci representing the placental mammal radiation ranged between 50% and 75%.

Exploring the Robustness of the Parsimonious Reconciliation Method in Host-Symbiont Cophylogeny

The aim of this paper is to explore the robustness of the parsimonious host-symbiont tree reconciliation method under editing or small perturbations of the input. The editing involves making different choices of unique symbiont mapping to a host in the case where multiple associations exist. This is made necessary by the fact that the tree reconciliation model is currently unable to handle such associations. The analysis performed could however also address the problem of errors. The perturbations are re-rootings of the symbiont tree to deal with a possibly wrong placement of the root specially in the case of fast-evolving species. In order to do this robustness analysis, we introduce a simulation scheme specifically designed for the host-symbiont cophylogeny context, as well as a measure to compare sets of tree reconciliations, both of which are of interest by themselves.

The morphophysiological dormancy in Amborella trichopoda seeds is a pleisiomorphic trait in angiosperms

BACKGROUND AND AIMS

Recent parsimony-based reconstructions suggest that seeds of early angiosperms had either morphophysiological or physiological dormancy, with the former considered as more probable. The aim of this study was to determine the class of seed dormancy present in Amborella trichopoda , the sole living representative of the most basal angiosperm lineage Amborellales, with a view to resolving fully the class of dormancy present at the base of the angiosperm clade.

METHODS

Drupes of A. trichopoda without fleshy parts were germinated and dissected to observe their structure and embryo growth. Pre-treatments including acid scarification, gibberellin treatment and seed excision were tested to determine their influence on dormancy breakage and germination. Character-state mapping by maximum parsimony, incorporating data from the present work and published sources, was then used to determine the likely class of dormancy present in early angiosperms.

KEY RESULTS

Germination in A. trichopoda requires a warm stratification period of at least approx. 90 d, which is followed by endosperm swelling, causing the water-permeable pericarp-mesocarp envelope to split open. The embryo then grows rapidly within the seed, to radicle emergence some 17 d later and cotyledon emergence after an additional 24 d. Gibberellin treatment, acid scarification and excision of seeds from the surrounding drupe tissues all promoted germination by shortening the initial phase of dormancy, prior to embryo growth.

CONCLUSIONS

Seeds of A. trichopoda have non-deep simple morphophysiological dormancy, in which mechanical resistance of the pericarp-mesocarp envelope plays a key role in the initial physiological phase. Maximum parsimony analyses, including data obtained in the present work, indicate that morphophysiological dormancy is likely to be a pleisiomorphic trait in flowering plants. The significance of this conclusion for studies of early angiosperm evolution is discussed.

Dioecy in Amborella trichopoda: evidence for genetically based sex determination and its consequences for inferences of the breeding system in early angiosperms

BACKGROUND AND AIMS

This work aimed to gain insight into the breeding system at the base of living angiosperms through both character state reconstructions and the study of sex ratios and phenotypes in the likely sister to all other living angiosperms, Amborella trichopoda .

METHODS

Sex phenotypes were mapped onto a phylogeny of basally diverging angiosperms using maximum parsimony. In parallel, sex ratios and phenotypes were studied over two consecutive flowering seasons in an ex situ population of A. trichopoda , while the sex ratio of an in situ population was also assessed.

KEY RESULTS

Parsimony analyses failed to resolve the breeding system present at the base of living angiosperms, but indicated the importance of A. trichopoda for the future elucidation of this question. The ex situ A. trichopoda population studied showed a primary sex ratio close to 1:1, though sex ratio bias was found in the in situ population studied. Instances of sexual instability were quantified in both populations.

CONCLUSIONS

Sex ratio data support the presence of genetic sex determination in A. trichopoda , whose further elucidation may guide inferences on the breeding system at the base of living angiosperms. Sexual instability in A. trichopoda suggests the operation of epigenetic mechanisms, and the evolution of dioecy via a gynodioecious intermediate.

Rooting Trees, Methods for

Identifying the root of a phylogenetic tree is important because incorrectly rooted phylogenetic trees may mislead evolutionary and taxonomic inferences. Many techniques for inferring the root have been proposed, but each has shortcomings that may make it inappropriate for any particular dataset. Here we outline the various ways to root phylogenetic trees, which include: outgroup, midpoint rooting, molecular clock rooting, and Bayesian molecular clock rooting. In addition, we discuss the pros and cons and also list software availability for each of the rooting methods.

Time-Scaled Evolutionary Analysis of the Transmission and Antibiotic Resistance Dynamics of Staphylococcus aureus Clonal Complex 398

Staphylococcus aureus clonal complex 398 (CC398) is associated with disease in humans and livestock, and its origins and transmission have generated considerable interest. We performed a time-scaled phylogenetic analysis of CC398, including sequenced isolates from the United Kingdom (Scotland), along with publicly available genomes. Using state-of-the-art methods for mapping traits onto phylogenies, we quantified transitions between host species to identify sink and source populations for CC398 and employed a novel approach to investigate the gain and loss of antibiotic resistance in CC398 over time. We identified distinct human- and livestock-associated CC398 clades and observed multiple transmissions of CC398 from livestock to humans and between countries, lending quantitative support to previous reports. Of note, we identified a subclade within the livestock-associated clade comprised of isolates from hospital environments and newborn babies, suggesting that livestock-associated CC398 is capable of onward transmission in hospitals. In addition, our analysis revealed significant differences in the dynamics of resistance to methicillin and tetracycline related to contrasting historical patterns of antibiotic usage between the livestock industry and human medicine. We also identified significant differences in patterns of gain and loss of different tetracycline resistance determinants, which we ascribe to epistatic interactions between the resistance genes and/or differences in the modes of inheritance of the resistance determinants.

Maximum likelihood inference of small trees in the presence of long branches

The statistical basis of maximum likelihood (ML), its robustness, and the fact that it appears to suffer less from biases lead to it being one of the most popular methods for tree reconstruction. Despite its popularity, very few analytical solutions for ML exist, so biases suffered by ML are not well understood. One possible bias is long branch attraction (LBA), a regularly cited term generally used to describe a propensity for long branches to be joined together in estimated trees. Although initially mentioned in connection with inconsistency of parsimony, LBA has been claimed to affect all major phylogenetic reconstruction methods, including ML. Despite the widespread use of this term in the literature, exactly what LBA is and what may be causing it is poorly understood, even for simple evolutionary models and small model trees. Studies looking at LBA have focused on the effect of two long branches on tree reconstruction. However, to understand the effect of two long branches it is also important to understand the effect of just one long branch. If ML struggles to reconstruct one long branch, then this may have an impact on LBA. In this study, we look at the effect of one long branch on three-taxon tree reconstruction. We show that, counterintuitively, long branches are preferentially placed at the tips of the tree. This can be understood through the use of analytical solutions to the ML equation and distance matrix methods. We go on to look at the placement of two long branches on four-taxon trees, showing that there is no attraction between long branches, but that for extreme branch lengths long branches are joined together disproportionally often. These results illustrate that even small model trees are still interesting to help understand how ML phylogenetic reconstruction works, and that LBA is a complicated phenomenon that deserves further study.

Long-branch attraction and the phylogeny of true water bugs (Hemiptera: Nepomorpha) as estimated from mitochondrial genomes

BACKGROUND

Most previous studies of morphological and molecular data have consistently supported the monophyly of the true water bugs (Hemiptera: Nepomorpha). An exception is a recent study by Hua et al. (BMC Evol Biol 9: 134, 2009) based on nine nepomorphan mitochondrial genomes. In the analysis of Hua et al. (BMC Evol Biol 9: 134, 2009), the water bugs in the group Pleoidea formed the sister group to a clade that consisted of Nepomorpha (the remaining true water bugs) + Leptopodomorpha (shore bugs) + Cimicomorpha (assassin bugs and relatives) + Pentatomomorpha (stink bugs and relatives), thereby suggesting that fully aquatic hemipterans evolved independently at least twice. Based on these results, Hua et al. (BMC Evol Biol 9: 134, 2009) elevated the Pleoidea to a new infraorder, the Plemorpha.

RESULTS

Our reanalysis suggests that the lack of support for the monophyly of the true water bugs (including Pleoidea) by Hua et al. (BMC Evol Biol 9: 134, 2009) likely resulted from inadequate taxon sampling. In particular, long-branch attraction (LBA) between the distant outgroup taxa and Pleoidea, as well as LBA among taxa in the ingroup, made Nepomorpha appear to be polyphyletic. We used three complementary strategies to test and alleviate the effects of LBA: (1) the removal of distant outgroups from the analysis; (2) the addition of closely related outgroups; and (3) the addition of a mitochondrial genome from a second family of Pleoidea. We also performed likelihood-ratio tests to examine the support for monophyly of Nepomorpha with different combinations of taxa included in the analysis. Furthermore, we found that specimens of Helotrephes sp. were misidentified as Paraplea frontalis (Fieber, 1844) by Hua et al. (BMC Evol Biol 9: 134, 2009).

CONCLUSIONS

All analyses that included the addition of more taxa significantly and consistently supported the placement of Pleoidea within the Nepomorpha (i.e., supported the monophyly of the traditional true water bugs). Our analyses further support a close relationship between Notonectoidea and Pleoidea within Nepomorpha, and the superfamilies Nepoidea, Ochteroidea, Naucoroidea, and Pleoidea are resolved as monophyletic in all trees with strong support. Our results also confirmed that monophyly of Nepomorpha clearly is not refuted by the mitochondrial genome data.

Cabomba as a model for studies of early angiosperm evolution

BACKGROUND

The angiosperms, or flowering plants, diversified in the Cretaceous to dominate almost all terrestrial environments. Molecular phylogenetic studies indicate that the orders Amborellales, Nymphaeales and Austrobaileyales, collectively termed the ANA grade, diverged as separate lineages from a remaining angiosperm clade at a very early stage in flowering plant evolution. By comparing these early diverging lineages, it is possible to infer the possible morphology and ecology of the last common ancestor of the extant angiosperms, and this analysis can now be extended to try to deduce the developmental mechanisms that were present in early flowering plants. However, not all species in the ANA grade form convenient molecular-genetic models.

SCOPE

The present study reviews the genus Cabomba (Nymphaeales), which shows a range of features that make it potentially useful as a genetic model. We focus on characters that have probably been conserved since the last common ancestor of the extant flowering plants. To facilitate the use of Cabomba as a molecular model, we describe methods for its cultivation to flowering in the laboratory, a novel Cabomba flower expressed sequence tag database, a well-adapted in situ hybridization protocol and a measurement of the nuclear genome size of C. caroliniana. We discuss the features required for species to become tractable models, and discuss the relative merits of Cabomba and other ANA-grade angiosperms in molecular-genetic studies aimed at understanding the origin of the flowering plants.

Parallel structural evolution of auxin response factors in the angiosperms

Here we analyze the structural evolution of the paralogous transcription factors ETTIN (ETT/ARF3) and AUXIN RESPONSE FACTOR 4 (ARF4), which control the development of floral organs and leaves in the model angiosperm Arabidopsis. ETT is truncated at its C terminus, and consequently lacks two regulatory domains present in most other ARFs, including ARF4. Our analysis indicates ETT and ARF4 to have been generated by the duplication of a non-truncated ARF gene prior to the radiation of the extant angiosperms. We furthermore show that either ETT or ARF4 orthologs have become modified to encode truncated ARF proteins, lacking C-terminal regulatory domains, in representatives of three groups that separated early in angiosperm evolution: Amborellales, Nymphaeales and the remaining angiosperm clade. Interestingly, the production of truncated ARF4 transcripts in Amborellales occurs through an alternative splicing mechanism, rather than through a permanent truncation, as in the other groups studied. To gain insight into the potential functional significance of truncations to ETT and ARF4, we tested the capacity of native, truncated and chimeric coding sequences of these genes to restore a wild-type phenotype to Arabidopsis ett mutants. We discuss the results of this analysis in the context of the structural evolution of ARF genes in the angiosperms.

Different gymnosperm outgroups have (mostly) congruent signal regarding the root of flowering plant phylogeny

We examined multiple plastid genes from a diversity of gymnosperm lineages to explore the consistency of signal among different outgroups for rooting flowering plant phylogeny. For maximum parsimony (MP), most outgroups attach on a branch of the underlying ingroup tree that leads to Amborella. Maximum likelihood (ML) analyses either root angiosperms on a nearby branch or find split support for these neighboring root placements, depending on the outgroup. The inclusion of two species of Hydatellaceae, recently recognized as an ancient line of angiosperms, does not aid in inference of the root. Cost profiles for placing the root in suboptimal locations are highly correlated across most outgroup comparisons, even comparing MP and ML profiles. Those for Gnetales are the most deviant of all those considered. This divergent outgroup either attaches on a long eudicot branch with moderate bootstrap support in MP analyses or supports no particular root location in ML analysis. Removing the most rapidly evolving sites in rate classifications based on two divergent angiosperm root placements with Gnetales yields strongly conflicting root placements in MP analysis, despite substantial overlap in the estimated sets of conservative sites. However, the generally high consistency in rooting signal among distantly related gymnosperm clades suggests that the long branch connecting angiosperms to their extant relatives may not interfere substantially with inference of the angiosperm root.

Floral variation and floral genetics in basal angiosperms

Recent advances in phylogeny reconstruction and floral genetics set the stage for new investigations of the origin and diversification of the flower. We review the current state of angiosperm phylogeny, with an emphasis on basal lineages. With the surprising inclusion of Hydatellaceae with Nymphaeales, recent studies support the topology of Amborella sister to all other extant angiosperms, with Nymphaeales and then Austrobaileyales as subsequent sisters to all remaining angiosperms. Notable modifications from most recent analyses are the sister relationships of Chloranthaceae with the magnoliids and of Ceratophyllaceae with eudicots. We review "trends" in floral morphology and contrast historical, intuitive interpretations with explicit character-state reconstructions using molecular-based trees, focusing on (1) the size, number, and organization of floral organs; (2) the evolution of the perianth; (3) floral symmetry; and (4) floral synorganization. We provide summaries of those genes known to affect floral features that contribute to much of floral diversity. Although most floral genes have not been investigated outside of a few model systems, sufficient information is emerging to identify candidate genes for testing specific hypotheses in nonmodel plants. We conclude with a set of evo-devo case studies in which floral genetics have been linked to variation in floral morphology.

An empirical test of the midpoint rooting method

The outgroup method is widely used to root phylogenetic trees. An accurate root indication, however, strongly depends on the availability of a proper outgroup. An alternate rooting method is the midpoint rooting (MPR). In this case, the root is set at the midpoint between the two most divergent operational taxonomic units. Although the midpoint rooting algorithm has been extensively used, the efficiency of this method in retrieving the correct root remains untested. In the present study, we empirically tested the success rate of the MPR in obtaining the outgroup root for a given phylogenetic tree. This was carried out by eliminating outgroups in 50 selected data sets from 33 papers and rooting the trees with the midpoint method. We were thus able to compare the root position retrieved by each method. Data sets were separated into three categories with different root consistencies: data sets with a single outgroup taxon (54% success rate for MPR), data sets with multiple outgroup taxa that showed inconsistency in root position (82% success rate), and data sets with multiple outgroup taxa in which root position was consistent (94% success rate). Interestingly, the more consistent the outgroup root is, the more successful MPR appears to be. This is a strong indication that the MPR method is valuable, particularly for cases where a proper outgroup is unavailable.

Phylogenetic and evolutionary implications of complete chloroplast genome sequences of four early-diverging angiosperms: Buxus (Buxaceae), Chloranthus (Chloranthaceae), Dioscorea (Dioscoreaceae), and Illicium (Schisandraceae)

We have determined the complete chloroplast genome sequences of four early-diverging lineages of angiosperms, Buxus (Buxaceae), Chloranthus (Chloranthaceae), Dioscorea (Dioscoreaceae), and Illicium (Schisandraceae), to examine the organization and evolution of plastid genomes and to estimate phylogenetic relationships among angiosperms. For the most part, the organization of these plastid genomes is quite similar to the ancestral angiosperm plastid genome with a few notable exceptions. Dioscorea has lost one protein-coding gene, rps16; this gene loss has also happened independently in four other land plant lineages, liverworts, conifers, Populus, and legumes. There has also been a small expansion of the inverted repeat (IR) in Dioscorea that has duplicated trnH-GUG. This event has also occurred multiple times in angiosperms, including in monocots, and in the two basal angiosperms Nuphar and Drimys. The Illicium chloroplast genome is unusual by having a 10 kb contraction of the IR. The four taxa sequenced represent key groups in resolving phylogenetic relationships among angiosperms. Illicium is one of the basal angiosperms in the Austrobaileyales, Chloranthus (Chloranthales) remains unplaced in angiosperm classifications, and Buxus and Dioscorea are early-diverging eudicots and monocots, respectively. We have used sequences for 61 shared protein-coding genes from these four genomes and combined them with sequences from 35 other genomes to estimate phylogenetic relationships using parsimony, likelihood, and Bayesian methods. There is strong congruence among the trees generated by the three methods, and most nodes have high levels of support. The results indicate that Amborella alone is sister to the remaining angiosperms; the Nymphaeales represent the next-diverging clade followed by Illicium; Chloranthus is sister to the magnoliids and together this group is sister to a large clade that includes eudicots and monocots; and Dioscorea represents an early-diverging lineage of monocots just internal to Acorus.

Pollen tube growth in association with a dry-type stigmatic transmitting tissue and extragynoecial compitum in the basal angiosperm Kadsura longipedunculata (Schisandraceae)

Spatial features of pollen tube growth and the composition of the extracellular matrix (ECM) of transmitting tissue in carpels of Kadsura longipedunculata, a member of the basal angiosperm taxon Schisandraceae, were characterized to identify features of transmitting tissue that might have been important for pollen-carpel interactions during the early history of angiosperms. In addition to growing extracellularly along epidermal cells that make up stigmatic crests of individual carpels, pollen tubes grow on abaxial carpel epidermal cells between unfused carpels along an extragynoecial compitum to subsequently enter an adjacent carpel, a feature important for enhancing seed set in apocarpous species. Histo- and immunochemical data indicated that transmitting tissue ECM is not freely flowing as previously hypothesized. Rather, the ECM is similar to that of a dry-type stigma whereby a cuticular boundary with associated esterase activity confines a matrix containing methyl-esterified homogalacturonans. The Schisandraceae joins an increasing number of basal angiosperm taxa that have a transmitting tissue ECM similar to a dry-type stigma, thereby challenging traditional views that the ancestral pollen tube pathway was similar to a wet-type stigma covered with a freely flowing exudate. Dry-type stigmas are posited to provide tighter control over pollen capture, retention, and germination than wet-type stigmas.

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

Background

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

Results

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

Conclusion

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

An Empirical Assessment of Long-Branch Attraction Artefacts in Deep Eukaryotic Phylogenomics

In the context of exponential growing molecular databases, it becomes increasingly easy to assemble large multigene data sets for phylogenomic studies. The expected increase of resolution due to the reduction of the sampling (stochastic) error is becoming a reality. However, the impact of systematic biases will also become more apparent or even dominant. We have chosen to study the case of the long-branch attraction artefact (LBA) using real instead of simulated sequences. Two fast-evolving eukaryotic lineages, whose evolutionary positions are well established, microsporidia and the nucleomorph of cryptophytes, were chosen as model species. A large data set was assembled (44 species, 133 genes, and 24,294 amino acid positions) and the resulting rooted eukaryotic phylogeny (using a distant archaeal outgroup) is positively misled by an LBA artefact despite the use of a maximum likelihood-based tree reconstruction method with a complex model of sequence evolution. When the fastest evolving proteins from the fast lineages are progressively removed (up to 90%), the bootstrap support for the apparently artefactual basal placement decreases to virtually 0%, and conversely only the expected placement, among all the possible locations of the fast-evolving species, receives increasing support that eventually converges to 100%. The percentage of removal of the fastest evolving proteins constitutes a reliable estimate of the sensitivity of phylogenetic inference to LBA. This protocol confirms that both a rich species sampling (especially the presence of a species that is closely related to the fast-evolving lineage) and a probabilistic method with a complex model are important to overcome the LBA artefact. Finally, we observed that phylogenetic inference methods perform strikingly better with simulated as opposed to real data, and suggest that testing the reliability of phylogenetic inference methods with simulated data leads to overconfidence in their performance. Although phylogenomic studies can be affected by systematic biases, the possibility of discarding a large amount of data containing most of the nonphylogenetic signal allows recovering a phylogeny that is less affected by systematic biases, while maintaining a high statistical support.

The {beta}-amylase genes of grasses and a phylogenetic analysis of the Triticeae (Poaceae)

There are two forms of β-amylase in the Triticeae crop plants wheat, barley, and rye: an endosperm-specific form encoded by two or three closely linked genes, and a tissue-ubiquitous form encoded by a single gene. Both rice and corn have one ubiquitously expressed form encoded by a single gene. This study focuses on two phylogenetic analyses of β-amylase gene sequences. First, a phylogenetic analysis of coding sequences from wheat, barley, rye, rice, and corn was expected to clarify the relationship between the endosperm-specific and tissue-ubiquitous forms of the protein. Instead, it illustrates possible effects of distant outgroups, based on conflicting patterns of character state variation consistent with different root positions. Next, a broad sample of the monogenomic Triticeae was included in a phylogenetic analysis based on sequences from a portion of the tissue-ubiquitous β-amylase gene. The results were compared to existing Triticeae gene trees, among which extensive conflict had been noted in the past. One additional gene tree has not completely clarified the complexity of the group, but has shed additional light on reticulate phylogenetic patterns within the tribe, including relationships involving Eremopyrum, Thinopyrum, and the Triticum/Aegilops group.

Archaeobatrachian paraphyly and pangaean diversification of crown-group frogs

Current models for the early diversification of living frogs inferred from morphological, ontogenetic, or DNA sequence data invoke very different scenarios of character evolution and biogeography. To explore central controversies on the phylogeny of Anura, we analyzed nearly 4000 base pairs of mitochondrial and nuclear DNA for the major frog lineages. Likelihood-based analyses of this data set are congruent with morphological evidence in supporting a paraphyletic arrangement of archaeobatrachian frogs, with an (Ascaphus + Leiopelma) clade as the sister-group of all other living anurans. The stability of this outcome is reinforced by screening for phylogenetic bias resulting from site-specific rate variation, homoplasy, or the obligatory use of distantly related outgroups. Twenty-one alternative branching and rooting hypotheses were evaluated using a nonparametric multicomparison test and parametric bootstrapping. Relaxed molecular clock estimates situate the emergence of crown-group anurans in the Triassic, approximately 55 million years prior to their first appearance in the fossil record. The existence of at least four extant frog lineages on the supercontinent Pangaea before its breakup gains support from the estimation that three early splits between Laurasia- and Gondwana-associated families coincide with the initial rifting of these landmasses. This observation outlines the potential significance of this breakup event in the formation of separate Mesozoic faunal assemblages in both hemispheres.

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.

Evidence that CRABS CLAW and TOUSLED have conserved their roles in carpel development since the ancestor of the extant angiosperms

The carpel is the female reproductive organ specific to flowering plants. We aim to define the genes that controlled carpel development in the common ancestor of this group as a step toward determining the molecular events that were responsible for the evolution of the carpel. CRABS CLAW (CRC) and TOUSLED (TSL) control important aspects of carpel development in the model plant, Arabidopsis thaliana. The basal angiosperm species Amborella trichopoda and Cabomba aquatica very likely represent the two most early diverging groups of flowering plants. We have identified putative orthologues of CRC and TSL from A. trichopoda and C. aquatica, respectively. We demonstrate the expression patterns of these genes in carpels to be very highly conserved, both spatially and temporally, with those of their Arabidopsis orthologues. We argue that CRC and TSL in Arabidopsis are likely to have conserved their respective roles in carpel development since the common ancestor of the living flowering plants. We conclude that a divergent role shown for the CRC orthologue in rice, DROOPING LEAF, most probably arose specifically in the monocot lineage. We show that, in addition to its expression in carpels, the TSL orthologue of C. aquatica is expressed in tissues that contribute to buoyancy and argue that its role in these tissues may have arisen later than its role in carpel development.

Effect of teosinte cytoplasmic genomes on maize phenotype

Determining the contribution of organelle genes to plant phenotype is hampered by several factors, including the paucity of variation in the plastid and mitochondrial genomes. To circumvent this problem, evolutionary divergence between maize (Zea mays ssp. mays) and the teosintes, its closest relatives, was utilized as a source of cytoplasmic genetic variation. Maize lines in which the maize organelle genomes were replaced through serial backcrossing by those representing the entire genus, yielding alloplasmic sublines, or cytolines were created. To avoid the confounding effects of segregating nuclear alleles, an inbred maize line was utilized. Cytolines with Z. mays teosinte cytoplasms were generally indistinguishable from maize. However, cytolines with cytoplasm from the more distantly related Z. luxurians, Z. diploperennis, or Z. perennis exhibited a plethora of differences in growth, development, morphology, and function. Significant differences were observed for 56 of the 58 characters studied. Each cytoline was significantly different from the inbred line for most characters. For a given character, variation was often greater among cytolines having cytoplasms from the same species than among those from different species. The characters differed largely independently of each other. These results suggest that the cytoplasm contributes significantly to a large proportion of plant traits and that many of the organelle genes are phenotypically important.

Phylogenetic inference in Rafflesiales: the influence of rate heterogeneity and horizontal gene transfer

Background

The phylogenetic relationships among the holoparasites of Rafflesiales have remained enigmatic for over a century. Recent molecular phylogenetic studies using the mitochondrial matR gene placed Rafflesia, Rhizanthes and Sapria (Rafflesiaceae s. str.) in the angiosperm order Malpighiales and Mitrastema (Mitrastemonaceae) in Ericales. These phylogenetic studies did not, however, sample two additional groups traditionally classified within Rafflesiales (Apodantheaceae and Cytinaceae). Here we provide molecular phylogenetic evidence using DNA sequence data from mitochondrial and nuclear genes for representatives of all genera in Rafflesiales.

Results

Our analyses indicate that the phylogenetic affinities of the large-flowered clade and Mitrastema, ascertained using mitochondrial matR, are congruent with results from nuclear SSU rDNA when these data are analyzed using maximum likelihood and Bayesian methods. The relationship of Cytinaceae to Malvales was recovered in all analyses. Relationships between Apodanthaceae and photosynthetic angiosperms varied depending upon the data partition: Malvales (3-gene), Cucurbitales (matR) or Fabales (atp1). The latter incongruencies suggest that horizontal gene transfer (HGT) may be affecting the mitochondrial gene topologies. The lack of association between Mitrastema and Ericales using atp1 is suggestive of HGT, but greater sampling within eudicots is needed to test this hypothesis further.

Conclusions

Rafflesiales are not monophyletic but composed of three or four independent lineages (families): Rafflesiaceae, Mitrastemonaceae, Apodanthaceae and Cytinaceae. Long-branch attraction appears to be misleading parsimony analyses of nuclear small-subunit rDNA data, but model-based methods (maximum likelihood and Bayesian analyses) recover a topology that is congruent with the mitochondrial matR gene tree, thus providing compelling evidence for organismal relationships. Horizontal gene transfer appears to be influencing only some taxa and some mitochondrial genes, thus indicating that the process is acting at the single gene (not whole genome) level.

Long branch attraction, taxon sampling, and the earliest angiosperms: Amborella or monocots?

Background

Numerous studies, using in aggregate some 28 genes, have achieved a consensus in recognizing three groups of plants, including Amborella, as comprising the basal-most grade of all other angiosperms. A major exception is the recent study by Goremykin et al. (2003; Mol. Biol. Evol. 20:1499–1505), whose analyses of 61 genes from 13 sequenced chloroplast genomes of land plants nearly always found 100% support for monocots as the deepest angiosperms relative to Amborella, Calycanthus, and eudicots. We hypothesized that this conflict reflects a misrooting of angiosperms resulting from inadequate taxon sampling, inappropriate phylogenetic methodology, and rapid evolution in the grass lineage used to represent monocots.

Results

We used two main approaches to test this hypothesis. First, we sequenced a large number of chloroplast genes from the monocot Acorus and added these plus previously sequenced Acorus genes to the Goremykin et al. (2003) dataset in order to explore the effects of altered monocot sampling under the same analytical conditions used in their study. With Acorus alone representing monocots, strongly supported Amborella-sister trees were obtained in all maximum likelihood and parsimony analyses, and in some distance-based analyses. Trees with both Acorus and grasses gave either a well-supported Amborella-sister topology or else a highly unlikely topology with 100% support for grasses-sister and paraphyly of monocots (i.e., Acorus sister to "dicots" rather than to grasses). Second, we reanalyzed the Goremykin et al. (2003) dataset focusing on methods designed to account for rate heterogeneity. These analyses supported an Amborella-sister hypothesis, with bootstrap support values often conflicting strongly with cognate analyses performed without allowing for rate heterogeneity. In addition, we carried out a limited set of analyses that included the chloroplast genome of Nymphaea, whose position as a basal angiosperm was also, and very recently, challenged.

Conclusions

These analyses show that Amborella (or Amborella plus Nymphaea), but not monocots, is the sister group of all other angiosperms among this limited set of taxa and that the grasses-sister topology is a long-branch-attraction artifact leading to incorrect rooting of angiosperms. These results highlight the danger of having lots of characters but too few and, especially, molecularly divergent taxa, a situation long recognized as potentially producing strongly misleading molecular trees. They also emphasize the importance in phylogenetic analysis of using appropriate evolutionary models.

The mitochondrial DNA of land plants: peculiarities in phylogenetic perspective

Land plants exhibit a significant evolutionary plasticity in their mitochondrial DNA (mtDNA), which contrasts with the more conservative evolution of their chloroplast genomes. Frequent genomic rearrangements, the incorporation of foreign DNA from the nuclear and chloroplast genomes, an ongoing transfer of genes to the nucleus in recent evolutionary times and the disruption of gene continuity in introns or exons are the hallmarks of plant mtDNA, at least in flowering plants. Peculiarities of gene expression, most notably RNA editing and trans-splicing, are significantly more pronounced in land plant mitochondria than in chloroplasts. At the same time, mtDNA is generally the most slowly evolving of the three plant cell genomes on the sequence level, with unique exceptions in only some plant lineages. The slow sequence evolution and a variable occurrence of introns in plant mtDNA provide an attractive reservoir of phylogenetic information to trace the phylogeny of older land plant clades, which is as yet not fully resolved. This review attempts to summarize the unique aspects of land plant mitochondrial evolution from a phylogenetic perspective.

The Causes of self-sterility in natural populations of the relictual angiosperm, Illicium floridanum (Illiciaceae)

BACKGROUND AND AIMS

Illicium floridanum, a species belonging to the basal extant angiosperm taxon Illiciaceae, reportedly exhibits self-incompatibility (SI). To date, the site and timing of SI within the carpel of this species remains unidentified. Thus, the objective of this research was to determine the cellular and temporal aspects of SI in I. floridanum.

METHODS

Following controlled application of cross- and self-pollen in natural populations of I. floridanum, embryo sac development and temporal aspects of stigma receptivity, as well as pollen tube growth, fertilization, and embryo and endosperm development, were investigated with the aid of light and fluorescence microscopy.

KEY RESULTS

Flowers of I. floridanum exhibited complete dichogamy whereby stigmas only supported cross- and self-pollen tube growth prior to anther dehiscence. In contrast to earlier reports of SI in this species, a prezygotic SI resulting in rejection of self-pollen tube growth at the stigma was absent and there were no significant differences between cross- versus self-pollen germination and pollen tube growth within the style and ovary during the first 5 d after pollination. Structural development of the four-celled embryo sac was not differentially influenced by pollen type as noted to occur in other angiosperms with late-acting ovarian SI. The ovule micropyle and embryo sac were penetrated equally by cross- and self-pollen tubes. In addition, there were no statistically significant differences in cross- versus self-fertilization. A resting zygote and multicellular endosperm at a variety of developmental stages was present by 30 d after application of cross- or self-pollen.

CONCLUSIONS

In the clear absence of a prezygotic SI that was previously reported to result in differential self-pollen tube growth at the stigma, self- sterility in I. floridanum is likely due to early-acting inbreeding depression, although late-acting post-zygotic ovarian SI cannot be ruled out.

Mitochondrial DNA sequences reveal the photosynthetic relatives of Rafflesia, the world's largest flower

All parasites are thought to have evolved from free-living ancestors. However, the ancestral conditions facilitating the shift to parasitism are unclear, particularly in plants because the phylogenetic position of many parasites is unknown. This is especially true for Rafflesia, an endophytic holoparasite that produces the largest flowers in the world and has defied confident phylogenetic placement since its discovery >180 years ago. Here we present results of a phylogenetic analysis of 95 species of seed plants designed to infer the position of Rafflesia in an evolutionary context using the mitochondrial gene matR (1,806 aligned base pairs). Overall, the estimated phylogenetic tree is highly congruent with independent analyses and provides a strongly supported placement of Rafflesia with the order Malpighiales, which includes poinsettias, violets, and passionflowers. Furthermore, the phylogenetic placement of Mitrastema, another enigmatic, holoparasitic angiosperm with the order Ericales (which includes blueberries and persimmons), was obtained with these data. Although traditionally classified together, Rafflesia and Mitrastema are only distantly related, implying that their endoparasitic habits result from convergent evolution. Our results indicate that the previous significant difficulties associated with phylogenetic placement of holoparasitic plants may be overcome by using mitochondrial DNA so that a broader understanding of the origins and evolution of parasitism may emerge.

Evolution of the APETALA3 and PISTILLATA lineages of MADS-box-containing genes in the basal angiosperms

The B class genes, including homologs of the Arabidopsis loci APETALA3 (AP3) and PISTILLATA (PI ), appear to play a conserved role in the determination of petal and stamen identity across core eudicot angiosperms. Understanding how and when these functions evolved is a critical component of elucidating the evolution of flowers, particularly the appearance of petaloid perianth organs. Before comparisons of gene expression patterns or functions can be made, however, it is necessary to establish the orthology of AP3 and PI homologs from basal angiosperms. Here, we report the identification and analysis of 29 new representatives of the B gene lineage from basal ANITA and magnoliid dicot angiosperms. These studies indicate that gene duplications have occurred at every phylogenetic level, both before and after the duplication that produced the separate AP3 and PI lineages. Comparison of genomic structure among PI homologs indicates that a 12-nucleotide deletion that had been considered synapomorphic for the whole PI lineage actually arose within the ANITA grade, after the split of the Nymphaeales but before the separation of the Austrobaileyales. Evidence for alternative splicing of the Nymphaea AP3 homolog is also presented. The implications of these findings for angiosperm systematics, the conservation of AP3 and PI gene function, and the evolution of the ABC program are discussed.

Computational analysis of the evolution of the structure and function of 1-deoxy-D-xylulose-5-phosphate synthase, a key regulator of the mevalonate-independent pathway in plants

We investigated molecular evolution of 1-deoxy-D-xylulose-5-phosphate synthase (DXS), an important regulatory enzyme of the mevalonate-independent pathway involved in terpenoid biosynthesis. Sequence alignment showed that some regions, likely to be functionally important, were highly conserved among all of the plant DXS sequences analysed. Phylogenetic trees were inferred using DXS sequences from 11 species of Angiosperms and showed the division of the sequences into two classes. Clustering of DXS sequences did not correspond to phylogenetic relationships among the plant species studied. There was no consistency in the similarity of the variable regions in the secondary structure of the DXS functional protein except for Capsicum and Lycopersicon, both members of the Solanaceae. Hydrophobicity plots for the functional region of DXS revealed great similarity in their hydrophobic structure, consistent with the phylogenetic trees inferred, and with eight prominent hydrophilic and hydrophobic peaks. We also observed a consistent set of features common to the DXS transit peptides studied. These features were the same hydrophobic slope, a hydrophobic region in residues 35-45, and, in eight of 12 sequences, a Pro-Pro-Thr sequence at the C-terminal end. The transit sequences are likely bipartite and contain features that suggest the DXS protein is not only targeted to the chloroplast, but also to the thylakoid. To our knowledge this is the first suggestion that DXS is located specifically in the chloroplast thylakoid.

Rooting Phylogenetic Trees with Distant Outgroups: A Case Study from the Commelinoid Monocots

Phylogenetic rooting experiments demonstrate that two chloroplast genes from commelinoid monocot taxa that represent the closest living relatives of the pickerelweed family, Pontederiaceae, retain measurable signals regarding the position of that family's root. The rooting preferences of the chloroplast sequences were compared with those for artificial sequences that correspond to outgroups so divergent that their signal has been lost completely. These random sequences prefer the three longest branches in the unrooted ingroup topology and do not preferentially root on the branches favored by real outgroup sequences. However, the rooting behavior of the artificial sequences is not a simple function of branch length. The random outgroups preferentially root on long terminal ingroup branches, but many ingroup branches comparable in length to those favored by random sequences attract no or few hits. Nonterminal ingroup branches are generally avoided, regardless of their length. Comparisons of the ease of forcing sequences onto suboptimal roots indicate that real outgroups require a substantially greater rooting penalty than random outgroups for around half of the least-parsimonious candidate roots. Although this supports the existence of nonrandomized signal in the real outgroups, it also indicates that there is little power to choose among the optimal and nearly optimal rooting possibilities. A likelihood-based test rejects the hypothesis that all rootings of the subtree using real outgroup sequences are equally good explanations of the data and also eliminates around half of the least optimal candidate roots. Adding genes or outgroups can improve the ability to discriminate among different root locations. Rooting discriminatory power is shown to be stronger, in general, for more closely related outgroups and is highly correlated among different real outgroups, genes, and optimality criteria.

The root of the angiosperms revisited

Most recent phylogenetic analyses of basal angiosperms have converged on the placement of Amborella as sister to all other extant angiosperms. However, certain recent studies suggest that Amborella and Nymphaeales (water lilies) form a clade sister to all remaining angiosperms or that Nymphaeales alone are the sister to the remaining angiosperms. We report here (i) maximum parsimony, maximum likelihood, and Bayesian phylogenetic analyses of 11 genes (>15,000 bp per taxon) for 16 taxa, (ii) maximum parsimony analysis for a subset of these genes for 104 taxa, and (iii) tests of alternative rootings with the nonparametric bootstrap and the likelihood ratio test with the parametric bootstrap. In addition, we use simulation analyses to examine the amount of bias that may be present in our methods of phylogeny estimation. Amborella continues to receive strong bootstrap support as the sister to all other extant angiosperms, and three of four tests reject alternative hypotheses of the angiosperm root. Although we cannot conclusively choose between Amborella vs. Amborella + Nymphaeales as sister to all other angiosperms, most analyses favor the former rooting.