Phylogenetic signal in nucleotide data from seed plants: implications for resolving the seed plant tree of life

Mitochondrial RNA editing sites affect the phylogenetic reconstruction of gymnosperms

•RNA editing sites may contain homoplasious signals that cause artifactual inferences in phylogenetic analyses.•Excluding RNA editing sites from gymnosperm mitochondrial genes restored the sister relationship of gnetophytes and Pinaceae.•Phylogenetic analysis based on mitochondrial genomic data should carefully evaluate the impact of RNA editing sites.

Phylogenomics and plastome evolution of Indigofera (Fabaceae)

Introduction

Indigofera L. is the third largest genus in Fabaceae and includes economically important species that are used for indigo dye-producing, medicinal, ornamental, and soil and water conservation. The genus is taxonomically difficult due to the high level of overlap in morphological characters of interspecies, fewer reliability states for classification, and extensive adaptive evolution. Previous characteristic-based taxonomy and nuclear ITS-based phylogenies have contributed to our understanding of Indigofera taxonomy and evolution. However, the lack of chloroplast genomic resources limits our comprehensive understanding of the phylogenetic relationships and evolutionary processes of Indigofera.

Methods

Here, we newly assembled 18 chloroplast genomes of Indigofera. We performed a series of analyses of genome structure, nucleotide diversity, phylogenetic analysis, species pairwise Ka/Ks ratios, and positive selection analysis by combining with allied species in Papilionoideae.

Results and discussion

The chloroplast genomes of Indigofera exhibited highly conserved structures and ranged in size from 157,918 to 160,040 bp, containing 83 protein-coding genes, 37 tRNA genes, and eight rRNA genes. Thirteen highly variable regions were identified, of which trnK-rbcL, ndhF-trnL, and ycf1 were considered as candidate DNA barcodes for species identification of Indigofera. Phylogenetic analysis using maximum likelihood (ML) and Bayesian inference (BI) methods based on complete chloroplast genome and protein-coding genes (PCGs) generated a well-resolved phylogeny of Indigofera and allied species. Indigofera monophyly was strongly supported, and four monophyletic lineages (i.e., the Pantropical, East Asian, Tethyan, and Palaeotropical clades) were resolved within the genus. The species pairwise Ka/Ks ratios showed values lower than 1, and 13 genes with significant posterior probabilities for codon sites were identified in the positive selection analysis using the branch-site model, eight of which were associated with photosynthesis. Positive selection of accD suggested that Indigofera species have experienced adaptive evolution to selection pressures imposed by their herbivores and pathogens. Our study provided insight into the structural variation of chloroplast genomes, phylogenetic relationships, and adaptive evolution in Indigofera. These results will facilitate future studies on species identification, interspecific and intraspecific delimitation, adaptive evolution, and the phylogenetic relationships of the genus Indigofera.

Cutting the long branches: Consilience as a path to unearth the evolutionary history of Gnetales

The Gnetales are one of the most fascinating groups within seed plants. Although the advent of molecular phylogenetics has generated some confidence in their phylogenetic placement of Gnetales within seed plants, their macroevolutionary history still presents many unknowns. Here, we review the reasons for such unknowns, and we focus the discussion on the presence of “long branches” both in their molecular and morphological history. The increased rate of molecular evolution and genome instability as well as the numerous unique traits (both reproductive and vegetative) in the Gnetales have been obstacles to a better understanding of their evolution. Moreover, the fossil record of the Gnetales, though relatively rich, has not yet been properly reviewed and investigated using a phylogenetic framework. Despite these apparent blocks to progress we identify new avenues to enable us to move forward. We suggest that a consilience approach, involving different disciplines such as developmental genetics, paleobotany, molecular phylogenetics, and traditional anatomy and morphology might help to “break” these long branches, leading to a deeper understanding of this mysterious group of plants.

Comparative Analyses of 3,654 Plastid Genomes Unravel Insights Into Evolutionary Dynamics and Phylogenetic Discordance of Green Plants

The plastid organelle is essential for many vital cellular processes and the growth and development of plants. The availability of a large number of complete plastid genomes could be effectively utilized to understand the evolution of the plastid genomes and phylogenetic relationships among plants. We comprehensively analyzed the plastid genomes of Viridiplantae comprising 3,654 taxa from 298 families and 111 orders and compared the genomic organizations in their plastid genomic DNA among major clades, which include gene gain/loss, gene copy number, GC content, and gene blocks. We discovered that some important genes that exhibit similar functions likely formed gene blocks, such as the psb family presumably showing co-occurrence and forming gene blocks in Viridiplantae. The inverted repeats (IRs) in plastid genomes have doubled in size across land plants, and their GC content is substantially higher than non-IR genes. By employing three different data sets [all nucleotide positions (nt123), only the first and second codon positions (nt12), and amino acids (AA)], our phylogenomic analyses revealed Chlorokybales + Mesostigmatales as the earliest-branching lineage of streptophytes. Hornworts, mosses, and liverworts forming a monophylum were identified as the sister lineage of tracheophytes. Based on nt12 and AA data sets, monocots, Chloranthales and magnoliids are successive sister lineages to the eudicots + Ceratophyllales clade. The comprehensive taxon sampling and analysis of different data sets from plastid genomes recovered well-supported relationships of green plants, thereby contributing to resolving some long-standing uncertainties in the plant phylogeny.

Loss of the IR region in conifer plastomes: Changes in the selection pressure and substitution rate of protein-coding genes

Plastid genomes (plastomes) have a quadripartite structure, but some species have drastically reduced or lost inverted repeat (IR) regions. IR regions are important for genome stability and the evolution rate. In the evolutionary process of gymnosperms, the typical IRs of conifers were lost, possibly affecting the evolutionary rate and selection pressure of genomic protein-coding genes. In this study, we selected 78 gymnosperm species (51 genera, 13 families) for evolutionary analysis. The selection pressure analysis results showed that negative selection effects were detected in all 50 common genes. Among them, six genes in conifers had higher ω values than non-conifers, and 12 genes had lower ω values. The evolutionary rate analysis results showed that 9 of 50 common genes differed between conifers and non-conifers. It is more obvious that in non-conifers, the rates of psbA (trst, trsv, ratio, dN, dS, and ω) were 2.6- to 3.1-fold of conifers. In conifers, trsv, ratio, dN, dS, and ω of ycf2 were 1.2- to 3.6-fold of non-conifers. In addition, the evolution rate of ycf2 in the IR was significantly reduced. psbA is undergoing dynamic change, with an abnormally high evolution rate as a small portion of it enters the IR region. Although conifers have lost the typical IR regions, we detected no change in the substitution rate or selection pressure of most protein-coding genes due to gene function, plant habitat, or newly acquired IRs.

The Evolution of euAPETALA2 Genes in Vascular Plants: From Plesiomorphic Roles in Sporangia to Acquired Functions in Ovules and Fruits

The field of evolutionary developmental biology can help address how morphological novelties evolve, a key question in evolutionary biology. In Arabidopsis thaliana, APETALA2 (AP2) plays a role in the development of key plant innovations including seeds, flowers, and fruits. AP2 belongs to the AP2/ETHYLENE RESPONSIVE ELEMENT BINDING FACTOR family which has members in all viridiplantae, making it one of the oldest and most diverse gene lineages. One key subclade, present across vascular plants is the euAPETALA2 (euAP2) clade, whose founding member is AP2. We reconstructed the evolution of the euAP2 gene lineage in vascular plants to better understand its impact on the morphological evolution of plants, identifying seven major duplication events. We also performed spatiotemporal expression analyses of euAP2/TOE3 genes focusing on less explored vascular plant lineages, including ferns, gymnosperms, early diverging angiosperms and early diverging eudicots. Altogether, our data suggest that euAP2 genes originally contributed to spore and sporangium development, and were subsequently recruited to ovule, fruit and floral organ development. Finally, euAP2 protein sequences are highly conserved; therefore, changes in the role of euAP2 homologs during development are most likely due to changes in regulatory regions.

Plastid phylogenomic analyses of Fagales reveal signatures of conflict and ancient chloroplast capture

Plastid phylogenomic analyses have shed light on many recalcitrant relationships across the angiosperm Tree of Life and continue to play an important role in plant phylogenetics alongside nuclear data sets given the utility of plastomes for revealing ancient and recent introgression. Here we conduct a plastid phylogenomic study of Fagales, aimed at exploring contentious relationships (e.g., the placement of Myricaceae and some intergeneric relationships in Betulaceae, Juglandaceae, and Fagaceae) and dissecting conflicting phylogenetic signals across the plastome. Combining 102 newly sequenced samples with publically available plastomes, we analyzed a dataset including 256 species and 32 of the 34 total genera of Fagales, representing the largest plastome-based study of the order to date. We find strong support for a sister relationship between Myricaceae and Juglandaceae, as well as strongly supported conflicting signal for alternative generic relationships in Betulaceae and Juglandaceae. These conflicts highlight the sensitivity of plastid phylogenomic analyses to genic composition, perhaps due to the prevalence of uninformative loci and heterogeneity in signal across different regions of the plastome. Phylogenetic relationships were geographically structured in subfamily Quercoideae, with Quercus being non-monophyletic and its sections forming clades with co-distributed Old World or New World genera of Quercoideae. Compared against studies based on nuclear genes, these results suggest extensive introgression and chloroplast capture in the early diversification of Quercus and Quercoideae. This study provides a critical plastome perspective on Fagales phylogeny, setting the stage for future studies employing more extensive data from the nuclear genome.

ASTRAL-Pro: Quartet-Based Species-Tree Inference despite Paralogy

Phylogenetic inference from genome-wide data (phylogenomics) has revolutionized the study of evolution because it enables accounting for discordance among evolutionary histories across the genome. To this end, summary methods have been developed to allow accurate and scalable inference of species trees from gene trees. However, most of these methods, including the widely used ASTRAL, can only handle single-copy gene trees and do not attempt to model gene duplication and gene loss. As a result, most phylogenomic studies have focused on single-copy genes and have discarded large parts of the data. Here, we first propose a measure of quartet similarity between single-copy and multicopy trees that accounts for orthology and paralogy. We then introduce a method called ASTRAL-Pro (ASTRAL for PaRalogs and Orthologs) to find the species tree that optimizes our quartet similarity measure using dynamic programing. By studying its performance on an extensive collection of simulated data sets and on real data sets, we show that ASTRAL-Pro is more accurate than alternative methods.

Measuring phylogenetic signal between categorical traits and phylogenies

MOTIVATION

Determining whether a trait and phylogeny share some degree of phylogenetic signal is a flagship goal in evolutionary biology. Signatures of phylogenetic signal can assist the resolution of a broad range of evolutionary questions regarding the tempo and mode of phenotypic evolution. However, despite the considerable number of strategies to measure it, few and limited approaches exist for categorical traits. Here, we used the concept of Shannon entropy and propose the δ statistic for evaluating the degree of phylogenetic signal between a phylogeny and categorical traits.

RESULTS

We validated δ as a measure of phylogenetic signal: the higher the δ-value the higher the degree of phylogenetic signal between a given tree and a trait. Based on simulated data we proposed a threshold-based classification test to pinpoint cases of phylogenetic signal. The assessment of the test's specificity and sensitivity suggested that the δ approach should only be applied to 20 or more species. We have further tested the performance of δ in scenarios of branch length and topology uncertainty, unbiased and biased trait evolution and trait saturation. Our results showed that δ may be applied in a wide range of phylogenetic contexts. Finally, we investigated our method in 14 360 mammalian gene trees and found that olfactory receptor genes are significantly associated with the mammalian activity patterns, a result that is congruent with expectations and experiments from the literature. Our application shows that δ can successfully detect molecular signatures of phenotypic evolution. We conclude that δ represents a useful measure of phylogenetic signal since many phenotypes can only be measured in categories.

AVAILABILITY AND IMPLEMENTATION

https://github.com/mrborges23/delta_statistic.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Forcing external constraints on tree inference using ASTRAL

Background

To account for genome-wide discordance among gene trees, several widely-used methods seek to find a species tree with the minimum distance to input gene trees. To efficiently explore the large space of species trees, some of these methods, including ASTRAL, use dynamic programming (DP). The DP paradigm can restrict the search space, and thus, ASTRAL and similar methods use heuristic methods to define a restricted search space. However, arbitrary constraints provided by the user on the output tree cannot be trivially incorporated into such restrictions. The ability to infer trees that honor user-defined constraints is needed for many phylogenetic analyses, but no solution currently exists for constraining the output of ASTRAL.

Results

We introduce methods that enable the ASTRAL dynamic programming to infer constrained trees in an effective and scalable manner. To do so, we adopt a recently developed tree completion algorithm and extend it to allow multifurcating input and output trees. In simulation studies, we show that the approach for honoring constraints is both effective and fast. On real data, we show that constrained searches can help interrogate branches not recovered in the optimal ASTRAL tree to reveal support for alternative hypotheses.

Conclusions

The new algorithm is added ASTRAL to all user-provided constraints on the species tree.

Reconstructing Krassilovia mongolica supports recognition of a new and unusual group of Mesozoic conifers

Previously unrecognized anatomical features of the cone scales of the enigmatic Early Cretaceous conifer Krassilovia mongolica include the presence of transversely oriented paracytic stomata, which is unusual for all other extinct and extant conifers. Identical stomata are present on co-occurring broad, linear, multiveined leaves assigned to Podozamites harrisii, providing evidence that K. mongolica and P. harrisii are the seed cones and leaves of the same extinct plant. Phylogenetic analyses of the relationships of the reconstructed Krassilovia plant place it in an informal clade that we name the Krassilovia Clade, which also includes Swedenborgia cryptomerioides-Podozamites schenkii, and Cycadocarpidium erdmanni-Podozamites schenkii. All three of these plants have linear leaves that are relatively broad compared to most living conifers, and that are also multiveined with transversely oriented paracytic stomata. We propose that these may be general features of the Krassilovia Clade. Paracytic stomata, and other features of this new group, recall features of extant and fossil Gnetales, raising questions about the phylogenetic homogeneity of the conifer clade similar to those raised by phylogenetic analyses of molecular data.

Pollination Drop Proteome and Reproductive Organ Transcriptome Comparison in Gnetum Reveals Entomophilous Adaptation

Gnetum possesses morphologically bisexual but functionally unisexual reproductive structures that exude sugary pollination drops to attract insects. Previous studies have revealed that the arborescent species (G. gnemon L.) and the lianoid species (G. luofuense C.Y.Cheng) possess different pollination syndromes. This study compared the proteome in the pollination drops of these two species using label-free quantitative techniques. The transcriptomes of fertile reproductive units (FRUs) and sterile reproductive units (SRUs) for each species were furthermore compared using Illumina Hiseq sequencing, and integrated proteomic and transcriptomic analyses were subsequently performed. Our results show that the differentially expressed proteins between FRUs and SRUs were involved in carbohydrate metabolism, the biosynthesis of amino acids and ovule defense. In addition, the differentially expressed genes between the FRUs and SRUs (e.g., MADS-box genes) were engaged in reproductive development and the formation of pollination drops. The integrated protein-transcript analyses revealed that FRUs and their exudates were relatively conservative while the SRUs and their exudates were more diverse, probably functioning as pollinator attractants. The evolution of reproductive organs appears to be synchronized with changes in the pollination drop proteome of Gnetum, suggesting that insect-pollinated adaptations are not restricted to angiosperms but also occur in gymnosperms.

Multiple historical processes obscure phylogenetic relationships in a taxonomically difficult group (Lobariaceae, Ascomycota)

In the age of next-generation sequencing, the number of loci available for phylogenetic analyses has increased by orders of magnitude. But despite this dramatic increase in the amount of data, some phylogenomic studies have revealed rampant gene-tree discordance that can be caused by many historical processes, such as rapid diversification, gene duplication, or reticulate evolution. We used a target enrichment approach to sample 400 single-copy nuclear genes and estimate the phylogenetic relationships of 13 genera in the lichen-forming family Lobariaceae to address the effect of data type (nucleotides and amino acids) and phylogenetic reconstruction method (concatenation and species tree approaches). Furthermore, we examined datasets for evidence of historical processes, such as rapid diversification and reticulate evolution. We found incongruence associated with sequence data types (nucleotide vs. amino acid sequences) and with different methods of phylogenetic reconstruction (species tree vs. concatenation). The resulting phylogenetic trees provided evidence for rapid and reticulate evolution based on extremely short branches in the backbone of the phylogenies. The observed rapid and reticulate diversifications may explain conflicts among gene trees and the challenges to resolving evolutionary relationships. Based on divergence times, the diversification at the backbone occurred near the Cretaceous-Paleogene (K-Pg) boundary (65 Mya) which is consistent with other rapid diversifications in the tree of life. Although some phylogenetic relationships within the Lobariaceae family remain with low support, even with our powerful phylogenomic dataset of up to 376 genes, our use of target-capturing data allowed for the novel exploration of the mechanisms underlying phylogenetic and systematic incongruence.

Significance of Photosynthetic Characters in the Evolution of Asian Gnetum (Gnetales)

Gnetum is a genus in the Gnetales that has a unique but ambiguous placement within seed plant phylogeny. Previous studies have shown that Gnetum has lower values of photosynthetic characters than those of other seed plants, but few Gnetum species have been studied, and those that have been studied are restricted to narrow taxonomic and geographic ranges. In addition, the mechanism underlying the lower values of photosynthetic characters in Gnetum remains poorly understood. Here, we investigated the photosynthetic characters of a Chinese lianoid species, i.e., Gnetum parvifolium, and co-occurring woody angiosperms growing in the wild, as well as seedlings of five Chinese Gnetum species cultivated in a greenhouse. The five Gnetum species had considerably lower values for photosynthesis parameters (net photosynthetic rate, transpiration rate, intercellular CO₂ concentration, and stomatal conductance) than those of other seed plant representatives. Interrelated analyses revealed that the low photosynthetic capacity may be an intrinsic property of Gnetum, and may be associated with its evolutionary history. Comparison of the chloroplast genomes (cpDNAs) of Gnetum with those of other seed plant representatives revealed that 17 coding genes are absent from the cpDNAs of all species of Gnetum. This lack of multiple functional genes from the cpDNAs probably leads to the low photosynthetic rates of Gnetum. Our results provide a new perspective on the evolutionary history of the Gnetales, and on the ecophysiological and genomic attributes of tropical biomes in general. These results could also be useful for the breeding and cultivation of Gnetum.

Origin of the Taxaceae aril: evolutionary implications of seed-cone teratologies in Pseudotaxus chienii

BACKGROUND AND AIMS

Fleshy structures that promote biotic dispersal by ingestion have evolved many times in seed plants. Within the yew family Taxaceae sensu lato (six genera, including Cephalotaxus), it remains controversial whether the characteristic fleshy structure surrounding the seed is interpreted as a novel outgrowth of the base of the ovule (i.e. an aril) or a fleshy seed coat that is entirely derived from the integument (i.e. a sarcotesta). This paper presents a detailed study of both wild-type and teratological seed cones of Pseudotaxus chienii, including morphology, anatomy and ontogeny.

METHODS

Wild-type and teratological seed cones were investigated with the classical paraffin technique and subsequent astrablue/safranin staining and scanning electron microscopy.

KEY RESULTS

The wild-type seed cone of Pseudotaxus possesses a fleshy white aril that is cup-like and not entirely fused to the seed. In the teratological seed cones investigated, the aril was bilobed and consisted of two free halves. In both wild-type and teratological cones, the aril was initiated as two lateral primordia in a transverse plane, but in wild-type cones the two primordia became extended into a ring primordium, which grew apically, leading to the cup-like shape. The teratological seed cones lacked a ring primordium and the two lateral aril lobes remained free throughout their entire ontogeny, alternating with the scale-like leaves inserted below them on the same branch; in some cases, these leaves also became fleshy.

CONCLUSIONS

Based on the ontogeny and arrangement of the two fleshy aril lobes in the teratological seed cones of Pseudotaxus, we suggest that the typical aril of Taxaceae could be readily interpreted as a fused pair of strongly swollen leaves rather than a modified integument. Our investigations of the cup-like aril of Pseudotaxus demonstrate a similarity not only with other Taxaceae but also with relatively distantly related conifers such as Phyllocladus (Podocarpaceae).

Insights on the persistence of pines (Pinus species) in the Late Cretaceous and their increasing dominance in the Anthropocene

Although gymnosperms were nearly swept away by the rise of the angiosperms in the Late Cretaceous, conifers, and pines (Pinus species) in particular, survived and regained their dominance in some habitats. Diversification of pines into fire-avoiding (subgenus Haploxylon) and fire-adapted (subgenus Diploxylon) species occurred in response to abiotic and biotic factors in the Late Cretaceous such as competition with emerging angiosperms and changing fire regimes. Adaptations/traits that evolved in response to angiosperm-fuelled fire regimes and stressful environments in the Late Cretaceous were key to pine success and are also contributing to a new "pine rise" in some areas in the Anthropocene. Human-mediated activities exert both positive and negative impacts of range size and expansion and invasions of pines. Large-scale afforestation with pines, human-mediated changes to fire regimes, and other ecosystem processes are other contributing factors. We discuss traits that evolved in response to angiosperm-mediated fires and stressful environments in the Cretaceous and that continue to contribute to pine persistence and dominance and the numerous ways in which human activities favor pines.

Ginkgo biloba's footprint of dynamic Pleistocene history dates back only 390,000 years ago

BACKGROUND

At the end of the Pliocene and the beginning of Pleistocene glaciation and deglaciation cycles Ginkgo biloba went extinct all over the world, and only few populations remained in China in relict areas serving as sanctuary for Tertiary relict trees. Yet the status of these regions as refuge areas with naturally existing populations has been proven not earlier than one decade ago. Herein we elaborated the hypothesis that during the Pleistocene cooling periods G. biloba expanded its distribution range in China repeatedly. Whole plastid genomes were sequenced, assembled and annotated, and sequence data was analyzed in a phylogenetic framework of the entire gymnosperms to establish a robust spatio-temporal framework for gymnosperms and in particular for G. biloba Pleistocene evolutionary history.

RESULTS

Using a phylogenetic approach, we identified that Ginkgoatae stem group age is about 325 million years, whereas crown group radiation of extant Ginkgo started not earlier than 390,000 years ago. During repeated warming phases, Gingko populations were separated and isolated by contraction of distribution range and retreated into mountainous regions serving as refuge for warm-temperate deciduous forests. Diversification and phylogenetic splits correlate with the onset of cooling phases when Ginkgo expanded its distribution range and gene pools merged.

CONCLUSIONS

Analysis of whole plastid genome sequence data representing the entire spatio-temporal genetic variation of wild extant Ginkgo populations revealed the deepest temporal footprint dating back to approximately 390,000 years ago. Present-day directional West-East admixture of genetic diversity is shown to be the result of pronounced effects of the last cooling period. Our evolutionary framework will serve as a conceptual roadmap for forthcoming genomic sequence data, which can then provide deep insights into the demographic history of Ginkgo.

Single-Copy Genes as Molecular Markers for Phylogenomic Studies in Seed Plants

Phylogenetic relationships among seed plant taxa, especially within the gymnosperms, remain contested. In contrast to angiosperms, for which several genomic, transcriptomic and phylogenetic resources are available, there are few, if any, molecular markers that allow broad comparisons among gymnosperm species. With few gymnosperm genomes available, recently obtained transcriptomes in gymnosperms are a great addition to identifying single-copy gene families as molecular markers for phylogenomic analysis in seed plants. Taking advantage of an increasing number of available genomes and transcriptomes, we identified single-copy genes in a broad collection of seed plants and used these to infer phylogenetic relationships between major seed plant taxa. This study aims at extending the current phylogenetic toolkit for seed plants, assessing its ability for resolving seed plant phylogeny, and discussing potential factors affecting phylogenetic reconstruction. In total, we identified 3,072 single-copy genes in 31 gymnosperms and 2,156 single-copy genes in 34 angiosperms. All studied seed plants shared 1,469 single-copy genes, which are generally involved in functions like DNA metabolism, cell cycle, and photosynthesis. A selected set of 106 single-copy genes provided good resolution for the seed plant phylogeny except for gnetophytes. Although some of our analyses support a sister relationship between gnetophytes and other gymnosperms, phylogenetic trees from concatenated alignments without 3rd codon positions and amino acid alignments under the CAT + GTR model, support gnetophytes as a sister group to Pinaceae. Our phylogenomic analyses demonstrate that, in general, single-copy genes can uncover both recent and deep divergences of seed plant phylogeny.

Dynamic Carboniferous tropical forests: new views of plant function and potential for physiological forcing of climate

Contents 1333 I. 1334 II. 1335 III. 1339 IV. 1344 V. 1347 VI. 1347 1348 1348 References 1348 SUMMARY: The Carboniferous, the time of Earth's penultimate icehouse and widespread coal formation, was dominated by extinct lineages of early-diverging vascular plants. Studies of nearest living relatives of key Carboniferous plants suggest that their physiologies and growth forms differed substantially from most types of modern vegetation, particularly forests. It remains a matter of debate precisely how differently and to what degree these long-extinct plants influenced the environment. Integrating biophysical analysis of stomatal and vascular conductivity with geochemical analysis of fossilized tissues and process-based ecosystem-scale modeling yields a dynamic and unique perspective on these paleoforests. This integrated approach indicates that key Carboniferous plants were capable of growth and transpiration rates that approach values found in extant crown-group angiosperms, differing greatly from comparatively modest rates found in their closest living relatives. Ecosystem modeling suggests that divergent stomatal conductance, leaf sizes and stem life span between dominant clades would have shifted the balance of soil-atmosphere water fluxes, and thus surface runoff flux, during repeated, climate-driven, vegetation turnovers. This synthesis highlights the importance of 'whole plant' physiological reconstruction of extinct plants and the potential of vascular plants to have influenced the Earth system hundreds of millions of years ago through vegetation-climate feedbacks.

β-Ginkgotides: Hyperdisulfide-constrained peptides from Ginkgo biloba

Hyperdisulfide-constrained peptides are distinguished by their high stability and diverse functions. Thus far, these peptides have been reported from animals only but their occurrence in plants are rare. Here, we report the discovery, synthesis and characterization of a hyperdisulfide-constrained peptides family of approximately 2 kDa, β-ginkgotides (β-gB1 and β-gB2) from Ginkgo biloba. Proteomic analysis showed β-ginkgotides contain 18‒20 amino acids, of which 16 residues form a conserved six-cysteine core with a highly clustered cysteine spacing of C‒CC‒C‒CC, an arrangement that has not been reported in cysteine-rich peptides. Disulfide mapping revealed a novel disulfide connectivity of CysI‒IV, CysII‒VI and CysIII‒V. Oxidative folding of synthetic β-gB1 to the native form was obtained in 70% yield. The synthetic β-gB1 displays a compact structure with no regular secondary structural elements, as determined by NMR spectroscopy. Transcriptomic analysis showed precursor βgb1 has a four-domain architecture and revealed an additional 76 β-ginkgotide-like peptides in 59 different gymnosperms, but none in angiosperms. Phylogenetic clustering analysis demonstrated β-ginkgotides belong to a new cysteine-rich peptide family. β-Ginkgotide is resistant to thermal, chemical and proteolytic degradation. Together, β-ginkgotides represent the first-in-class hyperdisulfide-constrained peptide family from plants with a novel scaffold that could be useful for engineering metabolically stable peptidyl therapeutics.

Diversification of Rosaceae since the Late Cretaceous based on plastid phylogenomics

Phylogenetic relationships in Rosaceae have long been problematic because of frequent hybridisation, apomixis and presumed rapid radiation, and their historical diversification has not been clarified. With 87 genera representing all subfamilies and tribes of Rosaceae and six of the other eight families of Rosales (outgroups), we analysed 130 newly sequenced plastomes together with 12 from GenBank in an attempt to reconstruct deep relationships and reveal temporal diversification of this family. Our results highlight the importance of improving sequence alignment and the use of appropriate substitution models in plastid phylogenomics. Three subfamilies and 16 tribes (as previously delimited) were strongly supported as monophyletic, and their relationships were fully resolved and strongly supported at most nodes. Rosaceae were estimated to have originated during the Late Cretaceous with evidence for rapid diversification events during several geological periods. The major lineages rapidly diversified in warm and wet habits during the Late Cretaceous, and the rapid diversification of genera from the early Oligocene onwards occurred in colder and drier environments. Plastid phylogenomics offers new and important insights into deep phylogenetic relationships and the diversification history of Rosaceae. The robust phylogenetic backbone and time estimates we provide establish a framework for future comparative studies on rosaceous evolution.

Ten genes and two topologies: an exploration of higher relationships in skipper butterflies (Hesperiidae)

Despite multiple attempts to infer the higher-level phylogenetic relationships of skipper butterflies (Family Hesperiidae), uncertainties in the deep clade relationships persist. The most recent phylogenetic analysis included fewer than 30% of known genera and data from three gene markers. Here we reconstruct the higher-level relationships with a rich sampling of ten nuclear and mitochondrial markers (7,726 bp) from 270 genera and find two distinct but equally plausible topologies among subfamilies at the base of the tree. In one set of analyses, the nuclear markers suggest two contrasting topologies, one of which is supported by the mitochondrial dataset. However, another set of analyses suggests mito-nuclear conflict as the reason for topological incongruence. Neither topology is strongly supported, and we conclude that there is insufficient phylogenetic evidence in the molecular dataset to resolve these relationships. Nevertheless, taking morphological characters into consideration, we suggest that one of the topologies is more likely.

The Complete Chloroplast Genome Sequence of Ampelopsis: Gene Organization, Comparative Analysis, and Phylogenetic Relationships to Other Angiosperms

Ampelopsis brevipedunculata is an economically important plant that belongs to the Vitaceae family of angiosperms. The phylogenetic placement of Vitaceae is still unresolved. Recent phylogenetic studies suggested that it should be placed in various alternative families including Caryophyllaceae, asteraceae, Saxifragaceae, Dilleniaceae, or with the rest of the rosid families. However, these analyses provided weak supportive results because they were based on only one of several genes. Accordingly, complete chloroplast genome sequences are required to resolve the phylogenetic relationships among angiosperms. Recent phylogenetic analyses based on the complete chloroplast genome sequence suggested strong support for the position of Vitaceae as the earliest diverging lineage of rosids and placed it as a sister to the remaining rosids. These studies also revealed relationships among several major lineages of angiosperms; however, they highlighted the significance of taxon sampling for obtaining accurate phylogenies. In the present study, we sequenced the complete chloroplast genome of A. brevipedunculata and used these data to assess the relationships among 32 angiosperms, including 18 taxa of rosids. The Ampelopsis chloroplast genome is 161,090 bp in length, and includes a pair of inverted repeats of 26,394 bp that are separated by small and large single copy regions of 19,036 bp and 89,266 bp, respectively. The gene content and order of Ampelopsis is identical to many other unrearranged angiosperm chloroplast genomes, including Vitis and tobacco. A phylogenetic tree constructed based on 70 protein-coding genes of 33 angiosperms showed that both Saxifragales and Vitaceae diverged from the rosid clade and formed two clades with 100% bootstrap value. The position of the Vitaceae is sister to Saxifragales, and both are the basal and earliest diverging lineages. Moreover, Saxifragales forms a sister clade to Vitaceae of rosids. Overall, the results of this study will contribute to better support of the evolution, molecular biology and genetic improvement of the plant Ampelopsis.

Molecular systematics of the freshwater stingrays (myliobatiformes: potamotrygonidae) of the Amazon, Orinoco, Magdalena, Esequibo, Caribbean, and Maracaibo basins (Colombia - Venezuela): evidence from three mitochondrial genes

Lack of adequate information about the taxonomic and evolutionary relationships, ecology, biology, and distribution of several species belonging to the family Potamotrygonidae makes these species vulnerable to anthropic activities, including commercial overexploitation for the ornamental fish market. The aim of this study was to investigate the systematic relationships among genera and species belonging to this family by analyses of three mitochondrial gene regions. Samples were collected from the main river basins in Colombia and Venezuela for four genera and seven species of the family, as well as for what appear to be unidentified species. Three mitochondrial molecular markers COI, Cytb, and ATP6 were amplified and sequenced. Maximum likelihood and Bayesian inference analysis were performed to obtain topologies for each marker and for a concatenated dataset including the three genes. Small dataset may compromise some methods estimations of sequence divergence in the ATP6 marker. Monophyly of the four genera in Potamotrygonidae was confirmed and phylogenetic relationships among members of the Potamotrygon genus were not clearly resolved. However, results obtained with the molecular marker Cytb appear to offer a good starting point to differentiate among genera and species as a tool that could be used for barcoding. The application of this gene as a barcode could be applied for management and regulation of extraction practices for these genera. Sequencing complete mitochondrial genomes would be the next step for testing evolutionary hypothesis among these genera. Population structure analyses should be undertaken for Paratrygon, Potamotrygon magdalenae and motoro.

Utility of characters evolving at diverse rates of evolution to resolve quartet trees with unequal branch lengths: analytical predictions of long-branch effects

BACKGROUND

The detection and avoidance of "long-branch effects" in phylogenetic inference represents a longstanding challenge for molecular phylogenetic investigations. A consequence of parallelism and convergence, long-branch effects arise in phylogenetic inference when there is unequal molecular divergence among lineages, and they can positively mislead inference based on parsimony especially, but also inference based on maximum likelihood and Bayesian approaches. Long-branch effects have been exhaustively examined by simulation studies that have compared the performance of different inference methods in specific model trees and branch length spaces.

RESULTS

In this paper, by generalizing the phylogenetic signal and noise analysis to quartets with uneven subtending branches, we quantify the utility of molecular characters for resolution of quartet phylogenies via parsimony. Our quantification incorporates contributions toward the correct tree from either signal or homoplasy (i.e. "the right result for either the right reason or the wrong reason"). We also characterize a highly conservative lower bound of utility that incorporates contributions to the correct tree only when they correspond to true, unobscured parsimony-informative sites (i.e. "the right result for the right reason"). We apply the generalized signal and noise analysis to classic quartet phylogenies in which long-branch effects can arise due to unequal rates of evolution or an asymmetrical topology. Application of the analysis leads to identification of branch length conditions in which inference will be inconsistent and reveals insights regarding how to improve sampling of molecular loci and taxa in order to correctly resolve phylogenies in which long-branch effects are hypothesized to exist.

CONCLUSIONS

The generalized signal and noise analysis provides analytical prediction of utility of characters evolving at diverse rates of evolution to resolve quartet phylogenies with unequal branch lengths. The analysis can be applied to identifying characters evolving at appropriate rates to resolve phylogenies in which long-branch effects are hypothesized to occur.

NDH expression marks major transitions in plant evolution and reveals coordinate intracellular gene loss

BACKGROUND

Key innovations have facilitated novel niche utilization, such as the movement of the algal predecessors of land plants into terrestrial habitats where drastic fluctuations in light intensity, ultraviolet radiation and water limitation required a number of adaptations. The NDH (NADH dehydrogenase-like) complex of Viridiplantae plastids participates in adapting the photosynthetic response to environmental stress, suggesting its involvement in the transition to terrestrial habitats. Although relatively rare, the loss or pseudogenization of plastid NDH genes is widely distributed across diverse lineages of photoautotrophic seed plants and mutants/transgenics lacking NDH function demonstrate little difference from wild type under non-stressed conditions. This study analyzes large transcriptomic and genomic datasets to evaluate the persistence and loss of NDH expression across plants.

RESULTS

Nuclear expression profiles showed accretion of the NDH gene complement at key transitions in land plant evolution, such as the transition to land and at the base of the angiosperm lineage. While detection of transcripts for a selection of non-NDH, photosynthesis related proteins was independent of the state of NDH, coordinate, lineage-specific loss of plastid NDH genes and expression of nuclear-encoded NDH subunits was documented in Pinaceae, gnetophytes, Orchidaceae and Geraniales confirming the independent and complete loss of NDH in these diverse seed plant taxa.

CONCLUSION

The broad phylogenetic distribution of NDH loss and the subtle phenotypes of mutants suggest that the NDH complex is of limited biological significance in contemporary plants. While NDH activity appears dispensable under favorable conditions, there were likely sufficiently frequent episodes of abiotic stress affecting terrestrial habitats to allow the retention of NDH activity. These findings reveal genetic factors influencing plant/environment interactions in a changing climate through 450 million years of land plant evolution.

Phylotranscriptomic analysis of the origin and early diversification of land plants

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

Phylogenetic signal detection from an ancient rapid radiation: Effects of noise reduction, long-branch attraction, and model selection in crown clade Apocynaceae

Crown clade Apocynaceae comprise seven primary lineages of lianas, shrubs, and herbs with a diversity of pollen aggregation morphologies including monads, tetrads, and pollinia, making them an ideal group for investigating the evolution and function of pollen packaging. Traditional molecular systematic approaches utilizing small amounts of sequence data have failed to resolve relationships along the spine of the crown clade, a likely ancient rapid radiation. The previous best estimate of the phylogeny was a five-way polytomy, leaving ambiguous the homology of aggregated pollen in two major lineages, the Periplocoideae, which possess pollen tetrads, and the milkweeds (Secamonoideae plus Asclepiadoideae), which possess pollinia. To assess whether greatly increased character sampling would resolve these relationships, a plastome sequence data matrix was assembled for 13 taxa of Apocynaceae, including nine newly generated complete plastomes, one partial new plastome, and three previously reported plastomes, collectively representing all primary crown clade lineages and outgroups. The effects of phylogenetic noise, long-branch attraction, and model selection (linked versus unlinked branch lengths among data partitions) were evaluated in a hypothesis-testing framework based on Shimodaira-Hasegawa tests. Discrimination among alternative crown clade resolutions was affected by all three factors. Exclusion of the noisiest alignment positions and topologies influenced by long-branch attraction resulted in a trichotomy along the spine of the crown clade consisting of Rhabdadenia+the Asian clade, Baisseeae+milkweeds, and Periplocoideae+the New World clade. Parsimony reconstruction on all optimal topologies after noise exclusion unambiguously supports parallel evolution of aggregated pollen in Periplocoideae (tetrads) and milkweeds (pollinia). Our phylogenomic approach has greatly advanced the resolution of one of the most perplexing radiations in Apocynaceae, providing the basis for study of convergent floral morphologies and their adaptive value.

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

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

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Phylogenomics and coalescent analyses resolve extant seed plant relationships

The extant seed plants include more than 260,000 species that belong to five main lineages: angiosperms, conifers, cycads, Ginkgo, and gnetophytes. Despite tremendous effort using molecular data, phylogenetic relationships among these five lineages remain uncertain. Here, we provide the first broad coalescent-based species tree estimation of seed plants using genome-scale nuclear and plastid data By incorporating 305 nuclear genes and 47 plastid genes from 14 species, we identify that i) extant gymnosperms (i.e., conifers, cycads, Ginkgo, and gnetophytes) are monophyletic, ii) gnetophytes exhibit discordant placements within conifers between their nuclear and plastid genomes, and iii) cycads plus Ginkgo form a clade that is sister to all remaining extant gymnosperms. We additionally observe that the placement of Ginkgo inferred from coalescent analyses is congruent across different nucleotide rate partitions. In contrast, the standard concatenation method produces strongly supported, but incongruent placements of Ginkgo between slow- and fast-evolving sites. Specifically, fast-evolving sites yield relationships in conflict with coalescent analyses. We hypothesize that this incongruence may be related to the way in which concatenation methods treat sites with elevated nucleotide substitution rates. More empirical and simulation investigations are needed to understand this potential weakness of concatenation methods.

Dancing together and separate again: gymnosperms exhibit frequent changes of fundamental 5S and 35S rRNA gene (rDNA) organisation

In higher eukaryotes, the 5S rRNA genes occur in tandem units and are arranged either separately (S-type arrangement) or linked to other repeated genes, in most cases to rDNA locus encoding 18S–5.8S–26S genes (L-type arrangement). Here we used Southern blot hybridisation, PCR and sequencing approaches to analyse genomic organisation of rRNA genes in all large gymnosperm groups, including Coniferales, Ginkgoales, Gnetales and Cycadales. The data are provided for 27 species (21 genera). The 5S units linked to the 35S rDNA units occur in some but not all Gnetales, Coniferales and in Ginkgo (∼30% of the species analysed), while the remaining exhibit separate organisation. The linked 5S rRNA genes may occur as single-copy insertions or as short tandems embedded in the 26S–18S rDNA intergenic spacer (IGS). The 5S transcript may be encoded by the same (Ginkgo, Ephedra) or opposite (Podocarpus) DNA strand as the 18S–5.8S–26S genes. In addition, pseudogenised 5S copies were also found in some IGS types. Both L- and S-type units have been largely homogenised across the genomes. Phylogenetic relationships based on the comparison of 5S coding sequences suggest that the 5S genes independently inserted IGS at least three times in the course of gymnosperm evolution. Frequent transpositions and rearrangements of basic units indicate relatively relaxed selection pressures imposed on genomic organisation of 5S genes in plants.

Chloroplast phylogenomics indicates that Ginkgo biloba is sister to cycads

Molecular phylogenetic studies have not yet reached a consensus on the placement of Ginkgoales, which is represented by the only living species, Ginkgo biloba (common name: ginkgo). At least six discrepant placements of ginkgo have been proposed. This study aimed to use the chloroplast phylogenomic approach to examine possible factors that lead to such disagreeing placements. We found the sequence types used in the analyses as the most critical factor in the conflicting placements of ginkgo. In addition, the placement of ginkgo varied in the trees inferred from nucleotide (NU) sequences, which notably depended on breadth of taxon sampling, tree-building methods, codon positions, positions of Gnetopsida (common name: gnetophytes), and including or excluding gnetophytes in data sets. In contrast, the trees inferred from amino acid (AA) sequences congruently supported the monophyly of a ginkgo and Cycadales (common name: cycads) clade, regardless of which factors were examined. Our site-stripping analysis further revealed that the high substitution saturation of NU sequences mainly derived from the third codon positions and contributed to the variable placements of ginkgo. In summary, the factors we surveyed did not affect results inferred from analyses of AA sequences. Congruent topologies in our AA trees give more confidence in supporting the ginkgo-cycad sister-group hypothesis.

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. [; ; ; ; ; .].

A functional phylogenomic view of the seed plants

A novel result of the current research is the development and implementation of a unique functional phylogenomic approach that explores the genomic origins of seed plant diversification. We first use 22,833 sets of orthologs from the nuclear genomes of 101 genera across land plants to reconstruct their phylogenetic relationships. One of the more salient results is the resolution of some enigmatic relationships in seed plant phylogeny, such as the placement of Gnetales as sister to the rest of the gymnosperms. In using this novel phylogenomic approach, we were also able to identify overrepresented functional gene ontology categories in genes that provide positive branch support for major nodes prompting new hypotheses for genes associated with the diversification of angiosperms. For example, RNA interference (RNAi) has played a significant role in the divergence of monocots from other angiosperms, which has experimental support in Arabidopsis and rice. This analysis also implied that the second largest subunit of RNA polymerase IV and V (NRPD2) played a prominent role in the divergence of gymnosperms. This hypothesis is supported by the lack of 24nt siRNA in conifers, the maternal control of small RNA in the seeds of flowering plants, and the emergence of double fertilization in angiosperms. Our approach takes advantage of genomic data to define orthologs, reconstruct relationships, and narrow down candidate genes involved in plant evolution within a phylogenomic view of species' diversification.

What is the phylogenetic signal limit from mitogenomes? The reconciliation between mitochondrial and nuclear data in the Insecta class phylogeny

BACKGROUND

Efforts to solve higher-level evolutionary relationships within the class Insecta by using mitochondrial genomic data are hindered due to fast sequence evolution of several groups, most notably Hymenoptera, Strepsiptera, Phthiraptera, Hemiptera and Thysanoptera. Accelerated rates of substitution on their sequences have been shown to have negative consequences in phylogenetic inference. In this study, we tested several methodological approaches to recover phylogenetic signal from whole mitochondrial genomes. As a model, we used two classical problems in insect phylogenetics: The relationships within Paraneoptera and within Holometabola. Moreover, we assessed the mitochondrial phylogenetic signal limits in the deeper Eumetabola dataset, and we studied the contribution of individual genes.

RESULTS

Long-branch attraction (LBA) artefacts were detected in all the datasets. Methods using Bayesian inference outperformed maximum likelihood approaches, and LBA was avoided in Paraneoptera and Holometabola when using protein sequences and the site-heterogeneous mixture model CAT. The better performance of this method was evidenced by resulting topologies matching generally accepted hypotheses based on nuclear and/or morphological data, and was confirmed by cross-validation and simulation analyses. Using the CAT model, the order Strepsiptera was recovered as sister to Coleoptera for the first time using mitochondrial sequences, in agreement with recent results based on large nuclear and morphological datasets. Also the Hymenoptera-Mecopterida association was obtained, leaving Coleoptera and Strepsiptera as the basal groups of the holometabolan insects, which coincides with one of the two main competing hypotheses. For the Paraneroptera, the currently accepted non-monophyly of Homoptera was documented as a phylogenetic novelty for mitochondrial data. However, results were not satisfactory when exploring the entire Eumetabola, revealing the limits of the phylogenetic signal that can be extracted from Insecta mitogenomes. Based on the combined use of the five best topology-performing genes we obtained comparable results to whole mitogenomes, highlighting the important role of data quality.

CONCLUSION

We show for the first time that mitogenomic data agrees with nuclear and morphological data for several of the most controversial insect evolutionary relationships, adding a new independent source of evidence to study relationships among insect orders. We propose that deeper divergences cannot be inferred with the current available methods due to sequence saturation and compositional bias inconsistencies. Our exploratory analysis indicates that the CAT model is the best dealing with LBA and it could be useful for other groups and datasets with similar phylogenetic difficulties.

An ILP solution for the gene duplication problem

Background

The gene duplication (GD) problem seeks a species tree that implies the fewest gene duplication events across a given collection of gene trees. Solving this problem makes it possible to use large gene families with complex histories of duplication and loss to infer phylogenetic trees. However, the GD problem is NP-hard, and therefore, most analyses use heuristics that lack any performance guarantee.

Results

We describe the first integer linear programming (ILP) formulation to solve instances of the gene duplication problem exactly. With simulations, we demonstrate that the ILP solution can solve problem instances with up to 14 taxa. Furthermore, we apply the new ILP solution to solve the gene duplication problem for the seed plant phylogeny using a 12-taxon, 6, 084-gene data set. The unique, optimal solution, which places Gnetales sister to the conifers, represents a new, large-scale genomic perspective on one of the most puzzling questions in plant systematics.

Conclusions

Although the GD problem is NP-hard, our novel ILP solution for it can solve instances with data sets consisting of as many as 14 taxa and 1, 000 genes in a few hours. These are the largest instances that have been solved to optimally to date. Thus, this work can provide large-scale genomic perspectives on phylogenetic questions that previously could only be addressed by heuristic estimates.

Comparative expression pattern analysis of WUSCHEL-related homeobox 2 (WOX2) and WOX8/9 in developing seeds and somatic embryos of the gymnosperm Picea abies

• In seed plants, current knowledge concerning embryonic pattern formation by polar auxin transport (PAT) and WUSCHEL-related homeobox (WOX) gene activity is primarily derived from studies on angiosperms, while less is known about these processes in gymnosperms. In view of the differences in their embryogeny, and the fact that somatic embryogenesis is used for mass propagation of conifers, a better understanding of embryo development is vital. • The expression patterns of PaWOX2 and PaWOX8/9 were followed with quantitative reverse transcription-polymerase chain reaction (qRT-PCR) and in situ hybridization (ISH) during seed and somatic embryo development in Norway spruce (Picea abies), and in somatic embryos treated with the PAT inhibitor N-1-naphthylphthalamic acid (NPA). • Both PaWOX2 and PaWOX8/9 were highly expressed at the early growth stages of zygotic and somatic embryos, and shared a similar expression pattern over the entire embryo. At later embryo stages, high expression of PaWOX8/9 became restricted to cotyledon primordia, epidermis, procambium and root apical meristem (RAM), which became most evident in NPA-treated somatic embryos, while expression of PaWOX2 was much lower. • Our results suggest an ancestral role of WOX in seed plant embryo development, and strengthen the proposed connection between PAT, PIN-FORMED (PIN) and WOX in the regulation of embryo patterning in seed plants.

Utility of a large, multigene plastid data set in inferring higher-order relationships in ferns and relatives (monilophytes)

PREMISE OF THE STUDY

The monilophytes (ferns and relatives)-the third largest group of land plants-exhibit a diverse array of vegetative and reproductive morphologies. Investigations into their early ecological and life-history diversification require accurate, well-corroborated phylogenetic estimates. We examined the utility of a large plastid-based data set in inferring backbone relationships for monilophytes. •

METHODS

We recovered 17 plastid genes for exemplar taxa using published and new primers. We compared results from maximum-likelihood and parsimony analyses, assessed the effects of removing rapidly evolving characters, and examined the extent to which our data corroborate or contradict the results of other studies, or resolve current ambiguities. •

KEY RESULTS

Considering multifamily clades, we found bootstrap support comparable to or better than that in published studies that used fewer genes from fewer or more taxa. We firmly establish filmy ferns (Hymenophyllales) as the sister group of all leptosporangiates except Osmundaceae, resolving the second deepest split in leptosporangiate-fern phylogeny. A clade comprising Ophioglossaceae and Psilotaceae is currently accepted as the sister group of other monilophytes, but we recover Equisetum in this position. We also recover marattioid and leptosporangiate ferns as sister groups. Maximum-likelihood rate-class estimates are somewhat skewed when a long-branch lineage (Selaginella) is included, negatively affecting bootstrap support for early branches. •

CONCLUSIONS

Our findings support the utility of this gene set in corroborating relationships found in previous studies, improving support, and resolving uncertainties in monilophyte phylogeny. Despite these advances, our results also underline the need for continued work on resolving the very earliest splits in monilophyte phylogeny.