Integrated Genomic Analyses From Low-Depth Sequencing Help Resolve Phylogenetic Incongruence in the Bamboos (Poaceae: Bambusoideae).

The bamboos (Bambusoideae, Poaceae) comprise a major grass lineage with a complex evolutionary history involving ancient hybridization and allopolyploidy. About 1700 described species are classified into three tribes, Olyreae (herbaceous bamboos), Bambuseae (tropical woody bamboos), and Arundinarieae (temperate woody bamboos). Nuclear analyses strongly support monophyly of the woody tribes, whereas plastome analyses strongly support paraphyly, with Bambuseae sister to Olyreae. Our objectives were to clarify the origin(s) of the woody bamboo tribes and resolve the nuclear vs. plastid conflict using genomic tools. For the first time, plastid and nuclear genomic information from the same bamboo species were combined in a single study. We sampled 51 species of bamboos representing the three tribes, estimated their genome sizes and generated low-depth sample sequence data, from which plastomes were assembled and nuclear repeats were analyzed. The distribution of repeat families was found to agree with nuclear gene phylogenies, but also provides novel insights into nuclear evolutionary history. We infer two early, independent hybridization events, one between an Olyreae ancestor and a woody ancestor giving rise to the two Bambuseae lineages, and another between two woody ancestors giving rise to the Arundinarieae. Retention of the Olyreae plastome associated with differential dominance of nuclear genomes and subsequent diploidization in some lineages explains the paraphyly observed in plastome phylogenetic estimations. We confirm ancient hybridization and allopolyploidy in the origins of the extant woody bamboo lineages and propose biased fractionation and diploidization as important factors in their evolution.

Leaf shape and size track habitat transitions across forest-grassland boundaries in the grass family (Poaceae).

Grass leaf shape is a strong indicator of their habitat with linear leaves predominating in open areas and ovate leaves distinguishing forest-associated grasses. This pattern among extant species suggests that ancestral shifts between forest and open habitats may have coincided with changes in leaf shape or size. We tested relationships between habitat, climate, photosynthetic pathway, and leaf shape and size in a phylogenetic framework to evaluate drivers of leaf shape and size variation over the evolutionary history of the family. We also estimated the ancestral habitat of Poaceae and tested whether forest margins served as transitional zones for shifts between forests and grasslands. We found that grass leaf shape is converging toward different shape optima in the forest understory, forest margins, and open habitats. Leaf size also varies with habitat. Grasses have smaller leaves in open and drier areas, and in areas with high solar irradiance. Direct transitions between linear and ovate leaves are rare as are direct shifts between forest and open habitats. The most likely ancestral habitat of the family was the forest understory and forest margins along with an intermediate leaf shape served as important transitional habitat and morphology, respectively, for subsequent shifts across forest-grassland biome boundaries.

A 250 plastome phylogeny of the grass family (Poaceae): topological support under different data partitions.

The systematics of grasses has advanced through applications of plastome phylogenomics, although studies have been largely limited to subfamilies or other subgroups of Poaceae. Here we present a plastome phylogenomic analysis of 250 complete plastomes (179 genera) sampled from 44 of the 52 tribes of Poaceae. Plastome sequences were determined from high throughput sequencing libraries and the assemblies represent over 28.7 Mbases of sequence data. Phylogenetic signal was characterized in 14 partitions, including (1) complete plastomes; (2) protein coding regions; (3) noncoding regions; and (4) three loci commonly used in single and multi-gene studies of grasses. Each of the four main partitions was further refined, alternatively including or excluding positively selected codons and also the gaps introduced by the alignment. All 76 protein coding plastome loci were found to be predominantly under purifying selection, but specific codons were found to be under positive selection in 65 loci. The loci that have been widely used in multi-gene phylogenetic studies had among the highest proportions of positively selected codons, suggesting caution in the interpretation of these earlier results. Plastome phylogenomic analyses confirmed the backbone topology for Poaceae with maximum bootstrap support (BP). Among the 14 analyses, 82 clades out of 309 resolved were maximally supported in all trees. Analyses of newly sequenced plastomes were in agreement with current classifications. Five of seven partitions in which alignment gaps were removed retrieved Panicoideae as sister to the remaining PACMAD subfamilies. Alternative topologies were recovered in trees from partitions that included alignment gaps. This suggests that ambiguities in aligning these uncertain regions might introduce a false signal. Resolution of these and other critical branch points in the phylogeny of Poaceae will help to better understand the selective forces that drove the radiation of the BOP and PACMAD clades comprising more than 99.9% of grass diversity.

Grass plastomes reveal unexpected paraphyly with endemic species of Micrairoideae from India and new haplotype markers in Arundinoideae.

PREMISE OF THE STUDY: We investigated the little-studied Arundinoideae/Micrairoideae clade of grasses with an innovative plastome phylogenomic approach. This method gives robust results for taxa of uncertain phylogenetic placement. Arundinoideae comprise ∼45 species, although historically was much larger. Arundinoideae is notable for the widely invasive Phragmites australis. Micrairoideae comprise nine genera and ∼200 species. Some are threatened with extinction, including Hubbardia, some Isachne spp., and Limnopoa. Two micrairoid genera, Eriachne and Pheidochloa, exhibit C4 photosynthesis in this otherwise C3 subfamily and represent an independent origin of the C4 pathway among grasses. METHODS: Five new plastomes were sequenced with next-generation sequencing-by-synthesis methods. Plastomes were assembled by de novo methods and phylogenetically analyzed with eight other recently published arundinoid or micrairoid plastomes and 11 outgroup species. Stable carbon isotope ratios were determined for micrairoid and arundinoid species to investigate ambiguities in the proxy evidence for C4 photosynthesis. KEY RESULTS: Phylogenomic analyses showed strong support for ingroup nodes in the Arundinoideae/Micrairoideae subtree, including a paraphyletic clade of Hubbardieae with Isachneae. Anatomical, biochemical, and positively selected sites data are ambiguous with regard to the photosynthetic pathways in Micrairoideae. Species of Hubbardia, Isachne, and Limnopoa were definitively shown by δ13C measurements to be C3 and Eriachne to be C4. CONCLUSIONS: Our plastome phylogenomic analyses for Micrairoideae are the first phylogenetic results to indicate paraphyly between Isachneae and Hubbardieae. The definitive δ13C data for four genera of Micrairoideae indicates the breadth of variation possible in the proxy evidence for photosynthetic pathways of both C3 and C4 taxa.

Plastid phylogenomic study of species within the genus Zea: rates and patterns of three classes of microstructural changes.

This project examines the relationships within the genus Zea using complete plastid genomes (plastomes). While Zea mays has been well studied, congeneric species have yet to be as thoroughly examined. For this study four complete plastomes and a fifth nearly complete plastome were sequenced in the five species (Zea diploperennis, Zea perennis, Zea luxurians, Zea nicaraguensis, and Zea mays subsp. huehuetenangensis) by Sanger or next-generation methods. An analysis of the microstructural changes, such as inversions, insertion or deletion mutations (indels) and determination of their frequencies were performed for the complete plastomes. It was determined that 193 indels and 15 inversions occurred across the examined plastomes of Zea. Tandem repeat indels were the most common type of microstructural change observed. Divergence times were estimated using a noncorrelated relaxed clock method. Divergence dates for specific nodes relative to Zea were calculated to fall between 38,000 years before present (YBP) for the subspecies included in this study and 23,000 YBP for section Luxuriantes included in this study. The stem lineage of all Zea species was calculated to have diverged at 176,000 YBP. The calculated mutation rates for the genus fell within the range of 1.7E-8 to 3.5E-8 microstructural changes per site per year. These rates of change are not uniform, despite the close relationships of taxa in this study. Phylogenomic analyses using full plastome alignments were also conducted to compare tree topologies from different types of mutations. In most cases, the previous work examining Zea mitochondrial and nuclear data was confirmed.

Correction: The First Complete Plastid Genome from Joinvilleaceae (J. ascendens; Poales) Shows Unique and Unpredicted Rearrangements.

[This corrects the article DOI: 10.1371/journal.pone.0163218.].

The First Complete Plastid Genome from Joinvilleaceae (J. ascendens; Poales) Shows Unique and Unpredicted Rearrangements.

Joinvilleaceae is a family of tropical grass-like monocots that comprises only the genus Joinvillea. Previous studies have placed Joinvilleaceae in close phylogenetic proximity to the well-studied grass family. A full plastome sequence was determined and characterized for J. ascendens. The plastome was sequenced with next generation methods, fully assembled de novo and annotated. The assembly revealed two novel inversions specific to the Joinvilleaceae lineage and at least one novel plastid inversion in the Joinvilleaceae-Poaceae lineage. Two previously documented inversions in the Joinvilleaceae-Poaceae lineage and one previously documented inversion in the Poaceae lineage were also verified. Inversion events were identified visually and verified computationally by simulation mutations. Additionally, the loss and subsequent degradation of the accD gene in order Poales was explored extensively in Poaceae and J. ascendens. The two novel inversions along with changes in gene composition between families better delimited lineages in the Poales. The presence of large inversions and subsequent reversals in this small family suggested a high potential for large-scale rearrangements to occur in plastid genomes.

Evolutionary relationships in Panicoid grasses based on plastome phylogenomics (Panicoideae; Poaceae).

BACKGROUND: Panicoideae are the second largest subfamily in Poaceae (grass family), with 212 genera and approximately 3316 species. Previous studies have begun to reveal relationships within the subfamily, but largely lack resolution and/or robust support for certain tribal and subtribal groups. This study aims to resolve these relationships, as well as characterize a putative mitochondrial insert in one linage. RESULTS: 35 newly sequenced Panicoideae plastomes were combined in a phylogenomic study with 37 other species: 15 Panicoideae and 22 from outgroups. A robust Panicoideae topology largely congruent with previous studies was obtained, but with some incongruences with previously reported subtribal relationships. A mitochondrial DNA (mtDNA) to plastid DNA (ptDNA) transfer was discovered in the Paspalum lineage. CONCLUSIONS: The phylogenomic analysis returned a topology that largely supports previous studies. Five previously recognized subtribes appear on the topology to be non-monophyletic. Additionally, evidence for mtDNA to ptDNA transfer was identified in both Paspalum fimbriatum and P. dilatatum, and suggests a single rare event that took place in a common progenitor. Finally, the framework from this study can guide larger whole plastome sampling to discern the relationships in Cyperochloeae, Steyermarkochloeae, Gynerieae, and other incertae sedis taxa that are weakly supported or unresolved.

Phylogenetic estimation and morphological evolution of Arundinarieae (Bambusoideae: Poaceae) based on plastome phylogenomic analysis.

We explored phylogenetic relationships among the twelve lineages of the temperate woody bamboo clade (tribe Arundinarieae) based on plastid genome (plastome) sequence data. A representative sample of 28 taxa was used and maximum parsimony, maximum likelihood and Bayesian inference analyses were conducted to estimate the Arundinarieae phylogeny. All the previously recognized clades of Arundinarieae were supported, with Ampelocalamus calcareus (Clade XI) as sister to the rest of the temperate woody bamboos. Well supported sister relationships between Bergbambos tessellata (Clade I) and Thamnocalamus spathiflorus (Clade VII) and between Kuruna (Clade XII) and Chimonocalmus (Clade III) were revealed by the current study. The plastome topology was tested by taxon removal experiments and alternative hypothesis testing and the results supported the current plastome phylogeny as robust. Neighbor-net analyses showed few phylogenetic signal conflicts, but suggested some potentially complex relationships among these taxa. Analyses of morphological character evolution of rhizomes and reproductive structures revealed that pachymorph rhizomes were most likely the ancestral state in Arundinarieae. In contrast leptomorph rhizomes either evolved once with reversions to the pachymorph condition or multiple times in Arundinarieae. Further, pseudospikelets evolved independently at least twice in the Arundinarieae, but the ancestral state is ambiguous.

The floral transcriptomes of four bamboo species (Bambusoideae; Poaceae): support for common ancestry among woody bamboos.

BACKGROUND: Next-generation sequencing now allows for total RNA extracts to be sequenced in non-model organisms such as bamboos, an economically and ecologically important group of grasses. Bamboos are divided into three lineages, two of which are woody perennials with bisexual flowers, which undergo gregarious monocarpy. The third lineage, which are herbaceous perennials, possesses unisexual flowers that undergo annual flowering events. RESULTS: Transcriptomes were assembled using both reference-based and de novo methods. These two methods were tested by characterizing transcriptome content using sequence alignment to previously characterized reference proteomes and by identifying Pfam domains. Because of the striking differences in floral morphology and phenology between the herbaceous and woody bamboo lineages, MADS-box genes, transcription factors that control floral development and timing, were characterized and analyzed in this study. Transcripts were identified using phylogenetic methods and categorized as A, B, C, D or E-class genes, which control floral development, or SOC or SVP-like genes, which control the timing of flowering events. Putative nuclear orthologues were also identified in bamboos to use as phylogenetic markers. CONCLUSIONS: Instances of gene copies exhibiting topological patterns that correspond to shared phenotypes were observed in several gene families including floral development and timing genes. Alignments and phylogenetic trees were generated for 3,878 genes and for all genes in a concatenated analysis. Both the concatenated analysis and those of 2,412 separate gene trees supported monophyly among the woody bamboos, which is incongruent with previous phylogenetic studies using plastid markers.

Phylogenomics and Plastome Evolution of Tropical Forest Grasses (Leptaspis, Streptochaeta: Poaceae).

Studies of complete plastomes have proven informative for our understanding of the molecular evolution and phylogenomics of grasses. In this study, a plastome phylogenomic analysis sampled species from lineages of deeply diverging grasses including Streptochaeta spicata (Anomochlooideae), Leptaspis banksii, and L. zeylanica (both Pharoideae). Plastomes from next generation sequences for three species were assembled by de novo methods. The unambiguously aligned coding and non-coding sequences of the entire plastomes were aligned with those from 43 other grasses and the outgroup Joinvillea ascendens. Outgroup sampling of grasses has previously posed a challenge for plastome phylogenomic studies because of major rearrangements of the plastome. Here, over 81,000 bases of homologous sequence were aligned for phylogenomic and divergence estimation analyses. Rare genomic changes, including persistently long ψycf1 and ψycf2 loci, the loss of the rpoC1 intron, and a 21 base tandem repeat insert in the coding sequence for rps19 defined branch points in the grass phylogeny. Marked differences were seen in the topologies inferred from the complete plastome and two gene matrices, and mean maximum likelihood support values for the former were 10% higher. In the full plastome phylogenomic analyses, the two species of Anomochlooideae were monophyletic. Leptaspis and Pharus were found to be reciprocally monophyletic, with the estimated divergence of two Leptaspis species preceding those of Pharus by over 14 Ma, consistent with historical biogeography. Our estimates for deep divergences among grasses were older than previous such estimates, likely influenced by more complete taxonomic and molecular sampling and the use of recently available or previously unused fossil calibration points.

Resolving deep relationships of PACMAD grasses: a phylogenomic approach.

BACKGROUND: Plastome sequences for 18 species of the PACMAD grasses (subfamilies Panicoideae, Aristidoideae, Chloridoideae, Micrairoideae, Arundinoideae, Danthonioideae) were analyzed phylogenomically. Next generation sequencing methods were used to provide complete plastome sequences for 12 species. Sanger sequencing was performed to determine the plastome of one species, Hakonechloa macra, to provide a reference for annotation. These analyses were conducted to resolve deep subfamilial relationships within the clade. Divergence estimates were assessed to determine potential factors that led to the rapid radiation of this lineage and its dominance of warmer open habitats. RESULTS: New plastomes were completely sequenced and characterized for 13 PACMAD species. An autapomorphic ~1140 bp deletion was found in Hakonechloa macra putatively pseudogenizing rpl14 and eliminating rpl16 from this plastome. Phylogenomic analyses support Panicoideae as the sister group to the ACMAD clade. Complete plastome sequences provide greater support at deep nodes within the PACMAD clade. The initial diversification of PACMAD subfamilies was estimated to occur at 32.4 mya. CONCLUSIONS: Phylogenomic analyses of complete plastomes provides resolution for deep relationships of PACMAD grasses. The divergence estimate of 32.4 mya at the crown node of the PACMAD clade coincides with the Eocene-Oligocene Transition (EOT). The Eocene was a period of global cooling and drying, which led to forest fragmentation and the expansion of open habitats now dominated by these grasses. Understanding how these grasses are related and determining a cause for their rapid radiation allows for future predictions of grassland distribution in the face of a changing global climate.

Plastid phylogenomics of the cool-season grass subfamily: clarification of relationships among early-diverging tribes.

Whole plastid genomes are being sequenced rapidly from across the green plant tree of life, and phylogenetic analyses of these are increasing resolution and support for relationships that have varied among or been unresolved in earlier single- and multi-gene studies. Pooideae, the cool-season grass lineage, is the largest of the 12 grass subfamilies and includes important temperate cereals, turf grasses and forage species. Although numerous studies of the phylogeny of the subfamily have been undertaken, relationships among some 'early-diverging' tribes conflict among studies, and some relationships among subtribes of Poeae have not yet been resolved. To address these issues, we newly sequenced 25 whole plastomes, which showed rearrangements typical of Poaceae. These plastomes represent 9 tribes and 11 subtribes of Pooideae, and were analysed with 20 existing plastomes for the subfamily. Maximum likelihood (ML), maximum parsimony (MP) and Bayesian inference (BI) robustly resolve most deep relationships in the subfamily. Complete plastome data provide increased nodal support compared with protein-coding data alone at nodes that are not maximally supported. Following the divergence of Brachyelytrum, Phaenospermateae, Brylkinieae-Meliceae and Ampelodesmeae-Stipeae are the successive sister groups of the rest of the subfamily. Ampelodesmeae are nested within Stipeae in the plastome trees, consistent with its hybrid origin between a phaenospermatoid and a stipoid grass (the maternal parent). The core Pooideae are strongly supported and include Brachypodieae, a Bromeae-Triticeae clade and Poeae. Within Poeae, a novel sister group relationship between Phalaridinae and Torreyochloinae is found, and the relative branching order of this clade and Aveninae, with respect to an Agrostidinae-Brizinae clade, are discordant between MP and ML/BI trees. Maximum likelihood and Bayesian analyses strongly support Airinae and Holcinae as the successive sister groups of a Dactylidinae-Loliinae clade.

Evolution of the bamboos (Bambusoideae; Poaceae): a full plastome phylogenomic analysis.

BACKGROUND: Bambusoideae (Poaceae) comprise three distinct and well-supported lineages: tropical woody bamboos (Bambuseae), temperate woody bamboos (Arundinarieae) and herbaceous bamboos (Olyreae). Phylogenetic studies using chloroplast markers have generally supported a sister relationship between Bambuseae and Olyreae. This suggests either at least two origins of the woody bamboo syndrome in this subfamily or its loss in Olyreae. RESULTS: Here a full chloroplast genome (plastome) phylogenomic study is presented using the coding and noncoding regions of 13 complete plastomes from the Bambuseae, eight from Olyreae and 10 from Arundinarieae. Trees generated using full plastome sequences support the previously recovered monophyletic relationship between Bambuseae and Olyreae. In addition to these relationships, several unique plastome features are uncovered including the first mitogenome-to-plastome horizontal gene transfer observed in monocots. CONCLUSIONS: Phylogenomic agreement with previous published phylogenies reinforces the validity of these studies. Additionally, this study presents the first published plastomes from Neotropical woody bamboos and the first full plastome phylogenomic study performed within the herbaceous bamboos. Although the phylogenomic tree presented in this study is largely robust, additional studies using nuclear genes support monophyly in woody bamboos as well as hybridization among previous woody bamboo lineages. The evolutionary history of the Bambusoideae could be further clarified using transcriptomic techniques to increase sampling among nuclear orthologues and investigate the molecular genetics underlying the development of woody and floral tissues.