Assessing the value of DNA barcodes for molecular phylogenetics: effect of increased taxon sampling in lepidoptera

Background

A common perception is that DNA barcode datamatrices have limited phylogenetic signal due to the small number of characters available per taxon. However, another school of thought suggests that the massively increased taxon sampling afforded through the use of DNA barcodes may considerably increase the phylogenetic signal present in a datamatrix. Here I test this hypothesis using a large dataset of macrolepidopteran DNA barcodes.

Methodology/Principal Findings

Taxon sampling was systematically increased in datamatrices containing macrolepidopteran DNA barcodes. Sixteen family groups were designated as concordance groups and two quantitative measures; the taxon consistency index and the taxon retention index, were used to assess any changes in phylogenetic signal as a result of the increase in taxon sampling. DNA barcodes alone, even with maximal taxon sampling (500 species per family), were not sufficient to reconstruct monophyly of families and increased taxon sampling generally increased the number of clades formed per family. However, the scores indicated a similar level of taxon retention (species from a family clustering together) in the cladograms as the number of species included in the datamatrix was increased, suggesting substantial phylogenetic signal below the ‘family’ branch.

Conclusions/Significance

The development of supermatrix, supertree or constrained tree approaches could enable the exploitation of the massive taxon sampling afforded through DNA barcodes for phylogenetics, connecting the twigs resolved by barcodes to the deep branches resolved through phylogenomics.     

The Complete Mitochondrial Genome of Gossypium hirsutum and Evolutionary Analysis of Higher Plant Mitochondrial Genomes

Background

Mitochondria are the main manufacturers of cellular ATP in eukaryotes. The plant mitochondrial genome contains large number of foreign DNA and repeated sequences undergone frequently intramolecular recombination. Upland Cotton (Gossypium hirsutum L.) is one of the main natural fiber crops and also an important oil-producing plant in the world. Sequencing of the cotton mitochondrial (mt) genome could be helpful for the evolution research of plant mt genomes.

Methodology/Principal Findings

We utilized 454 technology for sequencing and combined with Fosmid library of the Gossypium hirsutum mt genome screening and positive clones sequencing and conducted a series of evolutionary analysis on Cycas taitungensis and 24 angiosperms mt genomes. After data assembling and contigs joining, the complete mitochondrial genome sequence of G. hirsutum was obtained. The completed G.hirsutum mt genome is 621,884 bp in length, and contained 68 genes, including 35 protein genes, four rRNA genes and 29 tRNA genes. Five gene clusters are found conserved in all plant mt genomes; one and four clusters are specifically conserved in monocots and dicots, respectively. Homologous sequences are distributed along the plant mt genomes and species closely related share the most homologous sequences. For species that have both mt and chloroplast genome sequences available, we checked the location of cp-like migration and found several fragments closely linked with mitochondrial genes.

Conclusion

The G. hirsutum mt genome possesses most of the common characters of higher plant mt genomes. The existence of syntenic gene clusters, as well as the conservation of some intergenic sequences and genic content among the plant mt genomes suggest that evolution of mt genomes is consistent with plant taxonomy but independent among different species.     

Morphological evolution in the variable resin-producing Detarieae (Fabaceae): do morphological characters retain a phylogenetic signal?

Background and Aims

Previous molecular phylogenetic studies disagree with the informal generic-level taxonomic groups based on morphology. In this study morphological characters in the caesalpinioid clade Detarieae are evaluated within a phylogenetic framework as a means of better understanding phylogenetic relationships and morphological evolution.

Methods

Morphological characters were observed and scored for representative species of Detarieae focusing on the resin-producing genera. Phylogenetic analyses were carried out with morphological characters alone and then combined with DNA sequences.

Key Results

Despite a high level of homoplasy, morphological data support several clades corresponding to those recovered in molecular phylogenetic analyses. The more strongly supported clades are each defined by at least one morphological synapomorphy. Several characters (e.g. apetaly) previously used to define informal generic groups evolved several times independently, leading to the differences observed with the molecular phylogenetic analyses. Although floral evolution is complex in Detarieae some patterns are recovered.

Conclusions

New informal taxonomic groupings are proposed based on the present findings. Floral evolution in the diverse Detarieae clade is characterized by a repeated tendency toward zygomorphy through the reduction of lateral petals and toward complete loss of petals.Key words: Caesalpinioideae, Detarieae, floral evolution, Leguminosae, morphology, phylogeny, resins, taxonomy     

New insights into Trimezieae (Iridaceae) phylogeny: what do molecular data tell us?

Background and Aims

The Neotropical tribe Trimezieae are taxonomically difficult. They are generally characterized by the absence of the features used to delimit their sister group Tigridieae. Delimiting the four genera that make up Trimezieae is also problematic. Previous family-level phylogenetic analyses have not examined the monophyly of the tribe or relationships within it. Reconstructing the phylogeny of Trimezieae will allow us to evaluate the status of the tribe and genera and to examine the suitability of characters traditionally used in their taxonomy.

Methods

Maximum parsimony and Bayesian phylogenetic analyses are presented for 37 species representing all four genera of Trimezieae. Analyses were based on nrITS sequences and a combined plastid dataset. Ancestral character state reconstructions were used to investigate the evolution of ten morphological characters previously considered taxonomically useful.

Key Results

Analyses of nrITS and plastid datasets strongly support the monophyly of Trimezieae and recover four principal clades with varying levels of support; these clades do not correspond to the currently recognized genera. Relationships within the four clades are not consistently resolved, although the conflicting resolutions are not strongly supported in individual analyses. Ancestral character state reconstructions suggest considerable homoplasy, especially in the floral characters used to delimit Pseudotrimezia.

Conclusions

The results strongly support recognition of Trimezieae as a tribe but suggest that both generic- and species-level taxonomy need revision. Further molecular analyses, with increased sampling of taxa and markers, are needed to support any revision. Such analyses will help determine the causes of discordance between the plastid and nuclear data and provide a framework for identifying potential morphological synapomorphies for infra-tribal groups. The results also suggest Trimezieae provide a promising model for evolutionary research.     

The role of seed traits in determining the phylogenetic structure of temperate plant communities

Background and Aims

Phylogenetic clustering of species within plant communities can be expected to result from environmental filtering acting on an evolutionary-conserved plant trait. One such a candidate trait is the embryo to seed-size ratio (E:S). A high E:S may allow faster germination immediately after imbibition, and is therefore assumed to be advantageous in dry habitats. In this study the hypothesis was tested that habitat filtering driven by soil moisture conditions and acting on seed germination and seedling establishment is an important ecological mechanism in structuring temperate plant communities.

Methods

Vegetation samplings were performed in three habitats located within 200 km of each other in western Europe: Ellenberg indicator values showed that the habitats selected differed substantially in terms of soil moisture and light availability. E.S ratio and seed mass data for all genera were obtained from literature. Data were analysed using recently developed phylogenetic methods.

Key Results

Genera with a similar E:S tend to co-occur, as low and high E:S genera dominate in moist and dry habitats, respectively. A phylogenetically clustered pattern of community structure was evident, and dispersion of E:S was positively related to phylogenetic dispersion.

Conclusions

The phenotypically and phylogenetically clustered pattern indicates that E:S-mediated habitat filtering is an important assembly process structuring the plant community of the temperate climate habitats studied.     

Martharaptor greenriverensis, a New Theropod Dinosaur from the Lower Cretaceous of Utah

Background

The Yellow Cat Member of the Cedar Mountain Formation (Early Cretaceous, Barremian?) of Utah has yielded a rich dinosaur fauna, including the basal therizinosauroid theropod Falcarius utahensis at its base. Recent excavation uncovered a new possible therizinosauroid taxon from a higher stratigraphic level in the Cedar Mountain Formation than F. utahensis.

Methodology/Principal Findings

Here we describe a fragmentary skeleton of the new theropod and perform a phylogenetic analysis to determine its phylogenetic position. The skeleton includes fragments of vertebrae, a scapula, forelimb and hindlimb bones, and an ischium. It also includes several well-preserved manual unguals. Manual and pedal morphology show that the specimen is distinct from other theropods from the Cedar Mountain Formation and from previously described therizinosauroids. It is here named as the holotype of a new genus and species, Martharaptor greenriverensis. Phylogenetic analysis places M. greenriverensis within Therizinosauroidea as the sister taxon to Alxasaurus + Therizinosauridae, although support for this placement is weak.

Conclusions/Significance

The new specimen adds to the known dinosaurian fauna of the Yellow Cat Member of the Cedar Mountain Formation. If the phylogenetic placement is correct, it also adds to the known diversity of Therizinosauroidea.     

Universal Artifacts Affect the Branching of Phylogenetic Trees,Not Universal Scaling Laws

Background

The superficial resemblance of phylogenetic trees to other branching structures allows searching for macroevolutionary patterns. However, such trees are just statistical inferences of particular historical events. Recent meta-analyses report finding regularities in the branching pattern of phylogenetic trees. But is this supported by evidence, or are such regularities just methodological artifacts? If so, is there any signal in a phylogeny?

Methodology

In order to evaluate the impact of polytomies and imbalance on tree shape, the distribution of all binary and polytomic trees of up to 7 taxa was assessed in tree-shape space. The relationship between the proportion of outgroups and the amount of imbalance introduced with them was assessed applying four different tree-building methods to 100 combinations from a set of 10 ingroup and 9 outgroup species, and performing covariance analyses. The relevance of this analysis was explored taking 61 published phylogenies, based on nucleic acid sequences and involving various taxa, taxonomic levels, and tree-building methods.

Principal Findings

All methods of phylogenetic inference are quite sensitive to the artifacts introduced by outgroups. However, published phylogenies appear to be subject to a rather effective, albeit rather intuitive control against such artifacts. The data and methods used to build phylogenetic trees are varied, so any meta-analysis is subject to pitfalls due to their uneven intrinsic merits, which translate into artifacts in tree shape. The binary branching pattern is an imposition of methods, and seldom reflects true relationships in intraspecific analyses, yielding artifactual polytomies in short trees. Above the species level, the departure of real trees from simplistic random models is caused at least by two natural factors –uneven speciation and extinction rates; and artifacts such as choice of taxa included in the analysis, and imbalance introduced by outgroups and basal paraphyletic taxa. This artifactual imbalance accounts for tree shape convergence of large trees.

Significance

There is no evidence for any universal scaling in the tree of life. Instead, there is a need for improved methods of tree analysis that can be used to discriminate the noise due to outgroups from the phylogenetic signal within the taxon of interest, and to evaluate realistic models of evolution, correcting the retrospective perspective and explicitly recognizing extinction as a driving force. Artifacts are pervasive, and can only be overcome through understanding the structure and biological meaning of phylogenetic trees.Catalan Abstract in Translation S1.     

Exploring the use of cytochrome oxidase c subunit 1 (COI) for DNA barcoding of free-living marine nematodes

Background

The identification of free-living marine nematodes is difficult because of the paucity of easily scorable diagnostic morphological characters. Consequently, molecular identification tools could solve this problem. Unfortunately, hitherto most of these tools relied on 18S rDNA and 28S rDNA sequences, which often lack sufficient resolution at the species level. In contrast, only a few mitochondrial COI data are available for free-living marine nematodes. Therefore, we investigate the amplification and sequencing success of two partitions of the COI gene, the M1-M6 barcoding region and the I3-M11 partition.

Methodology

Both partitions were analysed in 41 nematode species from a wide phylogenetic range. The taxon specific primers for the I3-M11 partition outperformed the universal M1-M6 primers in terms of amplification success (87.8% vs. 65.8%, respectively) and produced a higher number of bidirectional COI sequences (65.8% vs 39.0%, respectively). A threshold value of 5% K2P genetic divergence marked a clear DNA barcoding gap separating intra- and interspecific distances: 99.3% of all interspecific comparisons were >0.05, while 99.5% of all intraspecific comparisons were <0.05 K2P distance.

Conclusion

The I3-M11 partition reliably identifies a wide range of marine nematodes, and our data show the need for a strict scrutiny of the obtained sequences, since contamination, nuclear pseudogenes and endosymbionts may confuse nematode species identification by COI sequences.     

Phylogeny of Campanuloideae (Campanulaceae) with Emphasis on the Utility of Nuclear Pentatricopeptide Repeat (PPR) Genes

Background

The Campanuloideae (Campanulaceae) are a highly diverse clade of angiosperms found mostly in the Northern Hemisphere, with the highest diversity in temperate areas of the Old World. Chloroplast markers have greatly improved our understanding of this clade but many relationships remain unclear primarily due to low levels of molecular evolution and recent and rapid divergence. Furthermore, focusing solely on maternally inherited markers such as those from the chloroplast genome may obscure processes such as hybridization. In this study we explore the phylogenetic utility of two low-copy nuclear loci from the pentatricopeptide repeat gene family (PPR). Rapidly evolving nuclear loci may provide increased phylogenetic resolution in clades containing recently diverged or closely related taxa. We present results based on both chloroplast and low-copy nuclear loci and discuss the utility of such markers to resolve evolutionary relationships and infer hybridization events within the Campanuloideae clade.

Results

The inclusion of low-copy nuclear genes into the analyses provides increased phylogenetic resolution in two species-rich clades containing recently diverged taxa. We also obtain support for the placement of two early diverging lineages (Jasione and Musschia-Gadellia clades) that have previously been unresolved. Furthermore, phylogenetic analyses of PPR loci revealed potential hybridization events for a number of taxa (e.g., Campanula pelviformis and Legousia species). These loci offer greater overall topological support than obtained with plastid DNA alone.

Conclusion

This study represents the first inclusion of low-copy nuclear genes for phylogenetic reconstruction in Campanuloideae. The two PPR loci were easy to sequence, required no cloning, and the sequence alignments were straightforward across the entire Campanuloideae clade. Although potentially complicated by incomplete lineage sorting, these markers proved useful for understanding the processes of reticulate evolution and resolving relationships at a wide range of phylogenetic levels. Our results stress the importance of including multiple, independent loci in phylogenetic analyses.     

Gene Order Phylogeny of the Genus Prochlorococcus

Background

Using gene order as a phylogenetic character has the potential to resolve previously unresolved species relationships. This character was used to resolve the evolutionary history within the genus Prochlorococcus, a group of marine cyanobacteria.

Methodology/Principal Findings

Orthologous gene sets and their genomic positions were identified from 12 species of Prochlorococcus and 1 outgroup species of Synechococcus. From this data, inversion and breakpoint distance-based phylogenetic trees were computed by GRAPPA and FastME. Statistical support of the resulting topology was obtained by application of a 50% jackknife resampling technique. The result was consistent and congruent with nucleotide sequence-based and gene-content based trees. Also, a previously unresolved clade was resolved, that of MIT9211 and SS120.

Conclusions/Significance

This is the first study to use gene order data to resolve a bacterial phylogeny at the genus level. It suggests that the technique is useful in resolving the Tree of Life.     

Highly incomplete taxa can rescue phylogenetic analyses from the negative impacts of limited taxon sampling

Background

Phylogenies are essential to many areas of biology, but phylogenetic methods may give incorrect estimates under some conditions. A potentially common scenario of this type is when few taxa are sampled and terminal branches for the sampled taxa are relatively long. However, the best solution in such cases (i.e., sampling more taxa versus more characters) has been highly controversial. A widespread assumption in this debate is that added taxa must be complete (no missing data) in order to save analyses from the negative impacts of limited taxon sampling. Here, we evaluate whether incomplete taxa can also rescue analyses under these conditions (empirically testing predictions from an earlier simulation study).

Methodology/Principal Findings

We utilize DNA sequence data from 16 vertebrate species with well-established phylogenetic relationships. In each replicate, we randomly sample 4 species, estimate their phylogeny (using Bayesian, likelihood, and parsimony methods), and then evaluate whether adding in the remaining 12 species (which have 50, 75, or 90% of their data replaced with missing data cells) can improve phylogenetic accuracy relative to analyzing the 4 complete taxa alone. We find that in those cases where sampling few taxa yields an incorrect estimate, adding taxa with 50% or 75% missing data can frequently (>75% of relevant replicates) rescue Bayesian and likelihood analyses, recovering accurate phylogenies for the original 4 taxa. Even taxa with 90% missing data can sometimes be beneficial.

Conclusions

We show that adding taxa that are highly incomplete can improve phylogenetic accuracy in cases where analyses are misled by limited taxon sampling. These surprising empirical results confirm those from simulations, and show that the benefits of adding taxa may be obtained with unexpectedly small amounts of data. These findings have important implications for the debate on sampling taxa versus characters, and for studies attempting to resolve difficult phylogenetic problems.     

Phylogenetic analyses provide the first insights into the evolution of OVATE family proteins in land plants

Background and Aims

The OVATE gene encodes a nuclear-localized regulatory protein belonging to a distinct family of plant-specific proteins known as the OVATE family proteins (OFPs). OVATE was first identified as a key regulator of fruit shape in tomato, with nonsense mutants displaying pear-shaped fruits. However, the role of OFPs in plant development has been poorly characterized.

Methods

Public databases were searched and a total of 265 putative OVATE protein sequences were identified from 13 sequenced plant genomes that represent the major evolutionary lineages of land plants. A phylogenetic analysis was conducted based on the alignment of the conserved OVATE domain from these 13 selected plant genomes. The expression patterns of tomato SlOFP genes were analysed via quantitative real-time PCR. The pattern of OVATE gene duplication resulting in the expansion of the gene family was determined in arabidopsis, rice and tomato.

Key Results

Genes for OFPs were found to be present in all the sampled land plant genomes, including the early-diverged lineages, mosses and lycophytes. Phylogenetic analysis based on the amino acid sequences of the conserved OVATE domain defined 11 sub-groups of OFPs in angiosperms. Different evolutionary mechanisms are proposed for OVATE family evolution, namely conserved evolution and divergent expansion. Characterization of the AtOFP family in arabidopsis, the OsOFP family in rice and the SlOFP family in tomato provided further details regarding the evolutionary framework and revealed a major contribution of tandem and segmental duplications towards expansion of the OVATE gene family.

Conclusions

This first genome-wide survey on OFPs provides new insights into the evolution of the OVATE protein family and establishes a solid base for future functional genomics studies on this important but poorly characterized regulatory protein family in plants.     

Wood anatomical correlates with theoretical conductivity and wood density across China: evolutionary evidence of the functional differentiation of axial and radial parenchyma

Background and Aims

In recent years considerable effort has focused on linking wood anatomy and key ecological traits. Studies analysing large databases have described how these ecological traits vary as a function of wood anatomical traits related to conduction and support, but have not considered how these functions interact with cells involved in storage of water and carbohydrates (i.e. parenchyma cells).

Methods

We analyzed, in a phylogenetic context, the functional relationship between cell types performing each of the three xylem functions (conduction, support and storage) and wood density and theoretical conductivity using a sample of approx. 800 tree species from China.

Key Results

Axial parenchyma and rays had distinct evolutionary correlation patterns. An evolutionary link was found between high conduction capacity and larger amounts of axial parenchyma that is probably related to water storage capacity and embolism repair, while larger amounts of ray tissue have evolved with increased mechanical support and reduced hydraulic capacity. In a phylogenetic principal component analysis this association of axial parenchyma with increased conduction capacity and rays with wood density represented orthogonal axes of variation. In multivariate space, however, the proportion of rays might be positively associated with conductance and negatively with wood density, indicating flexibility in these axes in species with wide rays.

Conclusions

The findings suggest that parenchyma types may differ in function. The functional axes represented by different cell types were conserved across lineages, suggesting a significant role in the ecological strategies of the angiosperms.     

A New Baurusuchid (Crocodyliformes, Mesoeucrocodylia) from the Late Cretaceous of Brazil and the Phylogeny of Baurusuchidae

Background

Baurusuchidae is a group of extinct Crocodyliformes with peculiar, dog-faced skulls, hypertrophied canines, and terrestrial, cursorial limb morphologies. Their importance for crocodyliform evolution and biogeography is widely recognized, and many new taxa have been recently described. In most phylogenetic analyses of Mesoeucrocodylia, the entire clade is represented only by Baurusuchus pachecoi, and no work has attempted to study the internal relationships of the group or diagnose the clade and its members.

Methodology/Principal Findings

Based on a nearly complete skull and a referred partial skull and lower jaw, we describe a new baurusuchid from the Vale do Rio do Peixe Formation (Bauru Group), Late Cretaceous of Brazil. The taxon is diagnosed by a suite of characters that include: four maxillary teeth, supratemporal fenestra with equally developed medial and anterior rims, four laterally visible quadrate fenestrae, lateral Eustachian foramina larger than medial Eustachian foramen, deep depression on the dorsal surface of pterygoid wing. The new taxon was compared to all other baurusuchids and their internal relationships were examined based on the maximum parsimony analysis of a discrete morphological data matrix.

Conclusion

The monophyly of Baurusuchidae is supported by a large number of unique characters implying an equally large morphological gap between the clade and its immediate outgroups. A complex phylogeny of baurusuchids was recovered. The internal branch pattern suggests two main lineages, one with a relatively broad geographical range between Argentina and Brazil (Pissarrachampsinae), which includes the new taxon, and an endemic clade of the Bauru Group in Brazil (Baurusuchinae).     

Phylogenetic analysis of pelecaniformes (aves) based on osteological data: implications for waterbird phylogeny and fossil calibration studies

Background

Debate regarding the monophyly and relationships of the avian order Pelecaniformes represents a classic example of discord between morphological and molecular estimates of phylogeny. This lack of consensus hampers interpretation of the group''s fossil record, which has major implications for understanding patterns of character evolution (e.g., the evolution of wing-propelled diving) and temporal diversification (e.g., the origins of modern families). Relationships of the Pelecaniformes were inferred through parsimony analyses of an osteological dataset encompassing 59 taxa and 464 characters. The relationships of the Plotopteridae, an extinct family of wing-propelled divers, and several other fossil pelecaniforms (Limnofregata, Prophaethon, Lithoptila, ?Borvocarbo stoeffelensis) were also assessed. The antiquity of these taxa and their purported status as stem members of extant families makes them valuable for studies of higher-level avian diversification.

Methodology/Principal Findings

Pelecaniform monophyly is not recovered, with Phaethontidae recovered as distantly related to all other pelecaniforms, which are supported as a monophyletic Steganopodes. Some anatomical partitions of the dataset possess different phylogenetic signals, and partitioned analyses reveal that these discrepancies are localized outside of Steganopodes, and primarily due to a few labile taxa. The Plotopteridae are recovered as the sister taxon to Phalacrocoracoidea, and the relationships of other fossil pelecaniforms representing key calibration points are well supported, including Limnofregata (sister taxon to Fregatidae), Prophaethon and Lithoptila (successive sister taxa to Phaethontidae), and ?Borvocarbo stoeffelensis (sister taxon to Phalacrocoracidae). These relationships are invariant when ‘backbone’ constraints based on recent avian phylogenies are imposed.

Conclusions/Significance

Relationships of extant pelecaniforms inferred from morphology are more congruent with molecular phylogenies than previously assumed, though notable conflicts remain. The phylogenetic position of the Plotopteridae implies that wing-propelled diving evolved independently in plotopterids and penguins, representing a remarkable case of convergent evolution. Despite robust support for the placement of fossil taxa representing key calibration points, the successive outgroup relationships of several “stem fossil + crown family” clades are variable and poorly supported across recent studies of avian phylogeny. Thus, the impact these fossils have on inferred patterns of temporal diversification depends heavily on the resolution of deep nodes in avian phylogeny.     

Predicted protein-protein interactions in the moss Physcomitrella patens: a new bioinformatic resource

Background

Physcomitrella patens, a haploid dominant plant, is fast becoming a useful molecular genetics and bioinformatics tool due to its key phylogenetic position as a bryophyte in the post-genomic era. Genome sequences from select reference species were compared bioinformatically to Physcomitrella patens using reciprocal blasts with the InParanoid software package. A reference protein interaction database assembled using MySQL by compiling BioGrid, BIND, DIP, and Intact databases was queried for moss orthologs existing for both interacting partners. This method has been used to successfully predict interactions for a number of angiosperm plants.

Results

The first predicted protein-protein interactome for a bryophyte based on the interolog method contains 67,740 unique interactions from 5,695 different Physcomitrella patens proteins. Most conserved interactions among proteins were those associated with metabolic processes. Over-represented Gene Ontology categories are reported here.

Conclusion

Addition of moss, a plant representative 200 million years diverged from angiosperms to interactomic research greatly expands the possibility of conducting comparative analyses giving tremendous insight into network evolution of land plants. This work helps demonstrate the utility of “guilt-by-association” models for predicting protein interactions, providing provisional roadmaps that can be explored using experimental approaches. Included with this dataset is a method for characterizing subnetworks and investigating specific processes, such as the Calvin-Benson-Bassham cycle.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-015-0524-1) contains supplementary material, which is available to authorized users.