Short Internal Branches Research Articles

Phylogenomics of the adaptive radiation of Triturus newts supports gradual ecological niche expansion towards an incrementally aquatic lifestyle

Newts of the genus Triturus (marbled and crested newts) exhibit substantial variation in the number of trunk vertebrae (NTV) and a higher NTV corresponds to a longer annual aquatic period. Because the Triturus phylogeny has thwarted resolution to date, the evolutionary history of NTV, annual aquatic period, and their potential coevolution has remained unclear. To resolve the phylogeny of Triturus, we generated a c. 6000 transcriptome-derived marker data set using a custom target enrichment probe set, and conducted phylogenetic analyses using: (1) data concatenation with RAxML, (2) gene-tree summary with ASTRAL, and (3) species-tree estimation with SNAPP. All analyses produce the same, highly supported topology, despite cladogenesis having occurred over a short timeframe, resulting in short internal branch lengths. Our new phylogenetic hypothesis is consistent with the minimal number of inferred changes in NTV count necessary to explain the diversity in NTV observed today. Although a causal relationship between NTV, body form, and aquatic ecology has yet to be experimentally established, our phylogeny indicates that these features have evolved together, and suggest that they may underlie the adaptive radiation that characterizes Triturus.

Multispecies coalescent analysis confirms standing phylogenetic instability in Hexapoda.

The multispecies coalescent (MSC) has been increasingly used in phylogenomic analyses due to the accommodation of gene tree topological heterogeneity by taking into account population-level processes, such as incomplete lineage sorting. In this sense, the phylogeny of insect species, which are characterized by their large effective population sizes, is suitable for a coalescent-based analysis. Furthermore, studies so far recovered short internal branches at early divergences of the insect tree of life, indicating fast evolutionary radiations that increase the probability of incomplete lineage sorting in deep time. Here, we investigated the performance of the MSC for a phylogenomic data set of hexapods compiled by Misof etal. (2014, Science 346:763). Our analysis recovered the monophyly of most insect orders, and major phylogenetic relationships were in agreement with current insect systematics. We identified, however, some evolutionary associations that were consistently problematic. Most noticeable, Hexapod monophyly was disrupted by the sister group relationship between the remiped crustacean and Insecta. Additionally, the interordinal relationships within Polyneoptera and Neuropteroidea were found to be phylogenetically unstable. We show that these controversial phylogenetic arrangements were also poorly supported by previous analyses, and therefore, we evaluated their robustness to stochastic errors from sampling sites and terminals, confirming standing problems in hexapod phylogeny in the genomics age.

Plastome sequences and exploration of tree-space help to resolve the phylogeny of riceflowers (Thymelaeaceae: Pimelea)

Data sets comprising small numbers of genetic markers are not always able to resolve phylogenetic relationships. This has frequently been the case in molecular systematic studies of plants, with many analyses being based on sequence data from only two or three chloroplast genes. An example of this comes from the riceflowers Pimelea Banks & Sol. ex Gaertn. (Thymelaeaceae), a large genus of flowering plants predominantly distributed in Australia. Despite the considerable morphological variation in the genus, low sequence divergence in chloroplast markers has led to the phylogeny of Pimelea remaining largely uncertain. In this study, we resolve the backbone of the phylogeny of Pimelea in comprehensive Bayesian and maximum-likelihood analyses of plastome sequences from 41 taxa. However, some relationships received only moderate to poor support, and the Pimelea clade contained extremely short internal branches. By using topology-clustering analyses, we demonstrate that conflicting phylogenetic signals can be found across the trees estimated from individual chloroplast protein-coding genes. A relaxed-clock dating analysis reveals that Pimelea arose in the mid-Miocene, with most divergences within the genus occurring during a subsequent rapid diversification. Our new phylogenetic estimate offers better resolution and is more strongly supported than previous estimates, providing a platform for future taxonomic revisions of both Pimelea and the broader subfamily. Our study has demonstrated the substantial improvements in phylogenetic resolution that can be achieved using plastome-scale data sets in plant molecular systematics.

Genome rearrangements and phylogeny reconstruction in Yersinia pestis

Genome rearrangements have played an important role in the evolution of Yersinia pestis from its progenitor Yersinia pseudotuberculosis. Traditional phylogenetic trees for Y. pestis based on sequence comparison have short internal branches and low bootstrap supports as only a small number of nucleotide substitutions have occurred. On the other hand, even a small number of genome rearrangements may resolve topological ambiguities in a phylogenetic tree. We reconstructed phylogenetic trees based on genome rearrangements using several popular approaches such as Maximum likelihood for Gene Order and the Bayesian model of genome rearrangements by inversions. We also reconciled phylogenetic trees for each of the three CRISPR loci to obtain an integrated scenario of the CRISPR cassette evolution. Analysis of contradictions between the obtained evolutionary trees yielded numerous parallel inversions and gain/loss events. Our data indicate that an integrated analysis of sequence-based and inversion-based trees enhances the resolution of phylogenetic reconstruction. In contrast, reconstructions of strain relationships based on solely CRISPR loci may not be reliable, as the history is obscured by large deletions, obliterating the order of spacer gains. Similarly, numerous parallel gene losses preclude reconstruction of phylogeny based on gene content.

Coalescent-Based Analyses of Genomic Sequence Data Provide a Robust Resolution of Phylogenetic Relationships among Major Groups of Gibbons.

The phylogenetic relationships among extant gibbon species remain unresolved despite numerous efforts using morphological, behavorial, and genetic data and the sequencing of whole genomes. A major challenge in reconstructing the gibbon phylogeny is the radiative speciation process, which resulted in extremely short internal branches in the species phylogeny and extensive incomplete lineage sorting with extensive gene-tree heterogeneity across the genome. Here, we analyze two genomic-scale data sets, with ∼10,000 putative noncoding and exonic loci, respectively, to estimate the species tree for the major groups of gibbons. We used the Bayesian full-likelihood method bpp under the multispecies coalescent model, which naturally accommodates incomplete lineage sorting and uncertainties in the gene trees. For comparison, we included three heuristic coalescent-based methods (mp-est, SVDQuartets, and astral) as well as concatenation. From both data sets, we infer the phylogeny for the four extant gibbon genera to be (Hylobates, (Nomascus, (Hoolock, Symphalangus))). We used simulation guided by the real data to evaluate the accuracy of the methods used. Astral, while not as efficient as bpp, performed well in estimation of the species tree even in presence of excessive incomplete lineage sorting. Concatenation, mp-est and SVDQuartets were unreliable when the species tree contains very short internal branches. Likelihood ratio test of gene flow suggests a small amount of migration from Hylobates moloch to H. pileatus, while cross-genera migration is absent or rare. Our results highlight the utility of coalescent-based methods in addressing challenging species tree problems characterized by short internal branches and rampant gene tree-species tree discordance.

A pilot study applying the plant Anchored Hybrid Enrichment method to New World sages (Salvia subgenus Calosphace; Lamiaceae)

We conducted a pilot study using Anchored Hybrid Enrichment to resolve relationships among a mostly Neotropical sage lineage that may have undergone a recent evolutionary radiation. Conventional markers (ITS, trnL-trnF and trnH-psbA) have not been able to resolve the relationships among species nor within portions of the backbone of the lineage. We sampled 12 representative species of subgenus Calosphace and included one species of Salvia’s s.l. closest relative, Lepechinia, as outgroup. Hybrid enrichment and sequencing were successful, yielding 448 alignments of individual loci with an average length of 704bp. The performance of the phylogenomic data in phylogenetic reconstruction was superior to that of conventional markers, increasing both support and resolution. Because the captured loci vary in the amount of net phylogenetic informativeness at different phylogenetic depths, these data are promising in phylogenetic reconstruction of this group and likely other lineages within Lamiales. However, special attention should be placed on the amount of phylogenetic noise that the data could potentially contain. A prior exploration step using phylogenetic informativeness profiles to detect loci with sites with disproportionately high substitution rates (showing “phantom” spikes) and, if required, the ensuing filtering of the problematic data is recommended. In our dataset, filtering resulted in increased support and resolution for the shallow nodes in maximum likelihood phylogenetic trees resulting from concatenated analyses of all the loci. Additionally, it is expected that an increase in sampling (loci and taxa) will aid in resolving weakly supported, short deep internal branches.

Diversification of Bromelioideae (Bromeliaceae) in the Brazilian Atlantic rainforest: A case study in Aechmea subgenus Ortgiesia

Aechmea subgenus Ortgiesia comprises ca. 20 species distributed in Brazil, Argentina, Paraguay, and Uruguay, with a center of diversity in the Brazilian Atlantic rainforest. We examined interspecific relationships of Ortgiesia based on Amplified Fragment Length Polymorphisms (AFLP). Ninety-six accessions belonging to 14 species of Ortgiesia were sampled, and genotyped with 11 AFLP primer combinations. The neighbor joining (NJ) tree depicted two main genetic groups within Aechmea subgenus Ortgiesia, and four subgroups. The NJ tree showed short internal branches, indicating an overall shallow genetic divergence among Ortgiesia species as expected for the recently radiated subfamily Bromelioideae. Our results suggest that hybridization and/or incomplete lineage sorting may have hampered the reconstruction of interspecific relationships in Aechmea subgenus Ortgiesia. The mapping of petal color (yellow, blue, pink, or white), inflorescence type (simple or compound), and inflorescence shape (ellipsoid, subcylindric, cylindric, or pyramidal) against the NJ tree indicated that these characters are of limited taxonomic use in Aechmea subgenus Ortgiesia due to homoplasy. An analysis of the current distribution of Ortgiesia identified the southern region of the Brazilian Atlantic rainforest, between latitudes of 26° and 27°S, as the center of diversity for the subgenus.

Phylogenomic analyses resolve an ancient trichotomy at the base of Ischyropsalidoidea (Arachnida, Opiliones) despite high levels of gene tree conflict and unequal minority resolution frequencies

Phylogenetic resolution of ancient rapid radiations has remained problematic despite major advances in statistical approaches and DNA sequencing technologies. Here we report on a combined phylogenetic approach utilizing transcriptome data in conjunction with Sanger sequence data to investigate a tandem of ancient divergences in the harvestmen superfamily Ischyropsalidoidea (Arachnida, Opiliones, Dyspnoi). We rely on Sanger sequences to resolve nodes within and between closely related genera, and use RNA-seq data from a subset of taxa to resolve a short and ancient internal branch. We use several analytical approaches to explore this succession of ancient diversification events, including concatenated and coalescent-based analyses and maximum likelihood gene trees for each locus. We evaluate the robustness of phylogenetic inferences using a randomized locus sub-sampling approach, and find congruence across these methods despite considerable incongruence across gene trees. Incongruent gene trees are not recovered in frequencies expected from a simple multispecies coalescent model, and we reject incomplete lineage sorting as the sole contributor to gene tree conflict. Using these approaches we attain robust support for higher-level phylogenetic relationships within Ischyropsalidoidea.

The Challenges of Resolving a Rapid, Recent Radiation: Empirical and Simulated Phylogenomics of Philippine Shrews.

Phylogenetic relationships in recent, rapid radiations can be difficult to resolve due to incomplete lineage sorting and reliance on genetic markers that evolve slowly relative to the rate of speciation. By incorporating hundreds to thousands of unlinked loci, phylogenomic analyses have the potential to mitigate these difficulties. Here, we attempt to resolve phylogenetic relationships among eight shrew species (genus Crocidura) from the Philippines, a phylogenetic problem that has proven intractable with small (< 10 loci) data sets. We sequenced hundreds of ultraconserved elements and whole mitochondrial genomes in these species and estimated phylogenies using concatenation, summary coalescent, and hierarchical coalescent methods. The concatenated approach recovered a maximally supported and fully resolved tree. In contrast, the coalescent-based approaches produced similar topologies, but each had several poorly supported nodes. Using simulations, we demonstrate that the concatenated tree could be positively misleading. Our simulations also show that the tree shape we tend to infer, which involves a series of short internal branches, is difficult to resolve, even if substitution models are known and multiple individuals per species are sampled. As such, the low support we obtained for backbone relationships in our coalescent-based inferences reflects a real and appropriate lack of certainty. Our results illuminate the challenges of estimating a bifurcating tree in a rapid and recent radiation, providing a rare empirical example of a nearly simultaneous series of speciation events in a terrestrial animal lineage as it spreads across an oceanic archipelago.

Coalescent methods are robust to the simultaneous effects of long branches and incomplete lineage sorting.

It is well known that species with elevated substitution rates can give rise to disproportionately long branches in the species tree. This combination of long and short branches can contribute to long-branch artifacts (LBA). Despite efforts to remedy LBA via increased taxon sampling and methodological improvements in gene tree estimation, it remains unclear how long and short branches affect species tree estimation in the presence of incomplete lineage sorting (ILS). Here, we examine the combined influence of long external and short internal branches on concatenation and coalescent methods using both simulated and empirical data. Our results demonstrate that the presence of long and short branches alone does not obviously confound the consistency of concatenation and coalescent methods. However, when long external and short internal branches occur simultaneously with high ILS, concatenation methods can be misled, especially when two of these long branches are sister lineages. In contrast, coalescent methods are more robust under these circumstances. This is particularly relevant because this topological pattern also characterizes numerous ancient rapid radiations across the tree of life. Because short internal branches can increase the potential for ILS and gene tree discordance, our results collectively suggest that coalescent methods are more likely to infer the correct species tree in cases of ancient rapid radiations where long external and short internal branches are in close phylogenetic proximity.

Coalescent versus concatenation methods and the placement of Amborella as sister to water lilies.

The molecular era has fundamentally reshaped our knowledge of the evolution and diversification of angiosperms. One outstanding question is the phylogenetic placement of Amborella trichopoda Baill., commonly thought to represent the first lineage of extant angiosperms. Here, we leverage publicly available data and provide a broad coalescent-based species tree estimation of 45 seed plants. By incorporating 310 nuclear genes, our coalescent analyses strongly support a clade containing Amborella plus water lilies (i.e., Nymphaeales) that is sister to all other angiosperms across different nucleotide rate partitions. Our results also show that commonly applied concatenation methods produce strongly supported, but incongruent placements of Amborella: slow-evolving nucleotide sites corroborate results from coalescent analyses, whereas fast-evolving sites place Amborella alone as the first lineage of extant angiosperms. We further explored the performance of coalescent versus concatenation methods using nucleotide sequences simulated on (i) the two alternate placements of Amborella with branch lengths and substitution model parameters estimated from each of the 310 nuclear genes and (ii) three hypothetical species trees that are topologically identical except with respect to the degree of deep coalescence and branch lengths. Our results collectively suggest that the Amborella alone placement inferred using concatenation methods is likely misled by fast-evolving sites. This appears to be exacerbated by the combination of long branches in stem group angiosperms, Amborella, and Nymphaeales with the short internal branch separating Amborella and Nymphaeales. In contrast, coalescent methods appear to be more robust to elevated substitution rates.

&lt;span lang="EN-GB"&gt;Phylogenetic relationships of Discyphus scopulariae (Orchidaceae, Cranichideae) inferred from plastid and nuclear DNA sequences: evidence supporting recognition of a new subtribe, Discyphinae&lt;/span&gt;

The monospecific genus Discyphus, previously considered a member of Spiranthinae (Orchidoideae: Cranichideae), displays both vegetative and floral morphological peculiarities that are out of place in that subtribe. These include a single, sessile, cordate leaf that clasps the base of the inflorescence and lies flat on the substrate, petals that are long-decurrent on the column, labellum margins free from sides of the column and a column provided with two separate, cup-shaped stigmatic areas. Because of its morphological uniqueness, the phylogenetic relationships of Discyphus have been considered obscure. In this study, we analyse nucleotide sequences of plastid and nuclear DNA under maximum parsimony and maximum likelihood criteria with the aim of clarifying its systematic position and discussing its peculiar morphology in an explicit phylogenetic context. Our analyses failed to support inclusion of Discyphus in Spiranthinae, signifying instead that this genus represents an additional isolated lineage of “core spiranthids.” The notable morphological disparity among such major lineages, as compared with the short internal branches subtending them in the molecular trees, would support the hypothesis that Discyphus represents a relict from an early radiation that also gave rise to Cranichidinae and Spiranthinae, putatively driven by adaptation to different pollinators given the morphological differences in floral morphology among these taxa.

Structure secondaire de l’ARN 18S et phylogénie basale du genre Carabus L., 1758 (Coleoptera : Carabidae)

SummaryThe phylogeny of the Carabinae (genera Cychrus, Pamborus, Ceroglossus, Calosoma, Carabus) is still not well resolved in its deepest nodes, despite several attempts in recent years. This group appeared during the Palaeogene, but phylogenetic analyses are difficult, which is correlated to the topology of the conflicting trees obtained, with long terminal and short internal branches. To circumvent this classical problem, we studied a non-coding and rather conservative nuclear marker (18SrRNA), then we focussed on its most variable regions, looking for some ‘molecular signatures’ originating from deep inter-nodes. For this, we compared the secondary structure of the helices H23-3 and H43. We found that the genus Carabus is monophyletic, the clade Calosoma being its sister group. Within Carabus, the Arcifera is the sister group of a diversified clade, named Eucarabi, grouping the Spinulati, the Crenolimbi and the Heterocarabi. Unexpectedly, the Heterocarabi is a very homogenous clade considering 18SrRNA, despite a high morphological diversity. The close phylogenetic relationship between Spinulati and Crenolimbi is emphasised. A clade grouping the genera Pamborus and Ceroglossus is a probable hypothesis but is not yet robustly demonstrated.

Discordance of Species Trees with Their Most Likely Gene Trees: A Unifying Principle

A labeled gene tree topology that disagrees with a labeled species tree topology is said to be anomalous if it is more probable under a coalescent model for gene lineage evolution than the labeled gene tree topology that matches the species tree. It has previously been shown that as a consequence of short internal branches of the species tree, for every labeled species tree topology with five or more taxa, and for asymmetric four-taxon species tree topologies, an assignment of species tree branch lengths can be made which gives rise to anomalous gene trees (AGTs). Here, I offer an alternative characterization of this result—a labeled species tree topology produces AGTs if and only if it contains two consecutive internal branches in an ancestor–descendant relationship—and I provide a proof that follows from the change in perspective. The reformulation and alternative proof of the existence result for AGTs provide the insight that it is not merely short internal branches that generate AGTs, but instead, short internal branches that are arranged consecutively.

Incomplete Lineage Sorting Is Common in Extant Gibbon Genera

We sequenced reduced representation libraries by means of Illumina technology to generate over 1.5 Mb of orthologous sequence from a representative of each of the four extant gibbon genera (Nomascus, Hylobates, Symphalangus, and Hoolock). We used these data to assess the evolutionary relationships between the genera by evaluating the likelihoods of all possible bifurcating trees involving the four taxa. Our analyses provide weak support for a tree with Nomascus and Hylobates as sister taxa and with Hoolock and Symphalangus as sister taxa, though bootstrap resampling suggests that other phylogenetic scenarios are also possible. This uncertainty is due to short internal branch lengths and extensive incomplete lineage sorting across taxa. The true phylogenetic relationships among gibbon genera will likely require a more extensive whole-genome sequence analysis.

Phylogenetic analysis of Primula section Primula reveals rampant non-monophyly among morphologically distinct species

The type section of Primula (Primulaceae), here considered to include seven species, is phylogenetically quite isolated in its genus. Although its species are popular ornamentals, traditional medicinal plants and model organisms for the study of heterostyly, the section has not yet been studied from a phylogenetic or evolutionary perspective. Using phylogenetic analysis of nuclear ITS and plastid data from all species and subspecies, we find that widespread Primula elatior is genetically heterogeneous and non-monophyletic to most if not all of the other ingroup taxa. The Genealogical Sorting Index (GSI) indicates that the assumption of all currently accepted species being independent lineages is consistent with the data. It is possible that P. elatior in its current circumscription may represents the disjointed remnant of an ancestral species from which the other recognized species diverged. However, based on available data, the alternative possibility of introgression explaining the non-monophyly of this species cannot be excluded. Species trees show P. elatior and P. veris as sister species. Primula vulgaris and P. juliae are closely related, while, in contrast to previous assumptions, P. renifolia does not appear to be a close relative of P. megaseifolia. With the section’s isolation from the rest of the genus and very short internal branches, our dataset also presents a case study of the confounding effects of different branch length priors on the Bayesian estimation of resulting branch length estimates. Experimental runs using different priors confirm the problem of resulting estimates varying by orders of magnitude, while topology and relative branch lengths seem unaffected.

Multiple gene sequences resolve phylogenetic relationships in the shorebird suborder Scolopaci (Aves: Charadriiformes)

Shorebirds (Charadriiformes) are a diverse assemblage of species renowned for their variation in behavior, morphology, and life-history traits, but comparative studies of trait variation remain limited by the lack of a well-supported phylogeny based on DNA sequences. In this study we build upon previous shorebird phylogenies to construct the first sequence-based species-level phylogeny for the Scolopaci, one of three shorebird suborders. We sampled 84 species in the Scolopaci, and collected data for five genes (one nuclear and four mitochondrial) via PCR and sequencing or from GenBank. The phylogeny was estimated using Bayesian inference on a partitioned dataset of 6365 aligned base pairs, and was well-supported except for the radiations within Tringa and Calidris. The shanks and phalaropes are sister to the snipes, woodcocks and dowitchers, which in turn are sister to the sandpipers. The godwits and curlews are successive sister-groups to these clades, and the morphologically disparate taxa (jacanas, painted snipes, seedsnipes, and the Plains-wanderer) are the basal sister-group in the Scolopaci. We show that Tringa, Gallinago, and Calidris are paraphyletic assemblages, and thus are in need of taxonomic revision. The clade of Calidridine sandpipers has very short internal branches indicative of a relatively recent rapid radiation, and will require a gene tree/species tree approach to resolve relationships among species.

The Evolution of Defense Mechanisms Correlate with the Explosive Diversification of Autodigesting Coprinellus Mushrooms (Agaricales, Fungi)

Bursts of diversification are known to have contributed significantly to the extant morphological and species diversity, but evidence for many of the theoretical predictions about adaptive radiations have remained contentious. Despite their tremendous diversity, patterns of evolutionary diversification and the contribution of explosive episodes in fungi are largely unknown. Here, using the genus Coprinellus (Psathyrellaceae, Agaricales) as a model, we report the first explosive fungal radiation and infer that the onset of the radiation correlates with a change from a multilayered to a much simpler defense structure on the fruiting bodies. We hypothesize that this change constitutes a key innovation, probably relaxing constraints on diversification imposed by nutritional investment into the development of protective tissues of fruiting bodies. Fossil calibration suggests that Coprinellus mushrooms radiated during the Miocene coinciding with global radiation of large grazing mammals following expansion of dry open grasslands. In addition to diversification rate-based methods, we test the hard polytomy hypothesis, by analyzing the resolvability of internal nodes of the backbone of the putative radiation using Reversible-Jump MCMC. We discuss potential applications and pitfalls of this approach as well as how biologically meaningful polytomies can be distinguished from alignment shortcomings. Our data provide insights into the nature of adaptive radiations in general by revealing a deceleration of morphological diversification through time. The dynamics of morphological diversification was approximated by obtaining the temporal distribution of state changes in discrete traits along the trees and comparing it with the tempo of lineage accumulation. We found that the number of state changes correlate with the number of lineages, even in parts of the tree with short internal branches, and peaks around the onset of the explosive radiation followed by a slowdown, most likely because of the decrease in available niches.

A new genus of miniaturized and pug-nosed gecko from South America (Sphaerodactylidae: Gekkota).

Sphaerodactyl geckos comprise five genera distributed across Central and South America and the Caribbean. We estimated phylogenetic relationships among sphaerodactyl genera using both separate and combined analyses of seven nuclear genes. Relationships among genera were incongruent at different loci and phylogenies were characterized by short, in some cases zero length, internal branches and poor phylogenetic support at most nodes. We recovered a polyphyletic Coleodactylus, with Coleodactylus amazonicus being deeply divergent from the remaining Coleodactylus species sampled. The C. amazonicus lineage possessed unique codon deletions in the genes PTPN12 and RBMX while the remaining Coleodactylus species had unique codon deletions in RAG1. Topology tests could not reject a monophyletic Coleodactylus, but we show that short internal branch lengths decreased the accuracy of topology tests because there were not enough data along short branches to support one phylogenetic hypothesis over another. Morphological data corroborated results of the molecular phylogeny, with Coleodactylus exhibiting substantial morphological heterogeneity. We identified a suite of unique craniofacial features that differentiate C. amazonicus not only from other Coleodactylus species, but also from all other geckos. We describe this novel sphaerodactyl lineage as a new genus, Chatogekko gen. nov. We present a detailed osteology of Chatogekko, characterizing osteological correlates of miniaturization that provide a framework for future studies in sphaerodactyl systematics and biology.

The phylogeny of the BEP clade in grasses revisited: Evidence from the whole-genome sequences of chloroplasts

Despite the considerable efforts to reconstruct the phylogeny of grasses, the relationships among the subfamilies Bambusoideae, Pooideae and Ehrhartoideae in the BEP clade remain unresolved. Here we completely sequenced three chloroplast genomes of representative species from Bambusoideae and Ehrhartoideae and obtained 19 additional chloroplast genome sequences of other grasses from GenBank. Using sequences of 76 chloroplast protein-coding genes from the 22 grass species, we fully resolved the phylogeny of the BEP clade. Our results strongly supported the (B,P)E hypothesis, i.e., Bambusoideae and Pooideae are more closely related than Ehrhartoideae. This result was not biased by systematic or sampling errors and was impervious to phylogenetic methods or model specification. The divergence time estimate suggests that the initial diversification of the BEP clade into three subfamilies happened within a short time period (∼4MY). The presence of these short internal branches may explain the inability of previous studies to achieve a confident resolution of the BEP clade. The combination of the sequences of the entire chloroplast genomes provided sufficient phylogenetic information to resolve the BEP phylogeny fully. These results provide a valuable evolutionary framework for comparative and functional genomic studies using the grass family as a model system.