Variability among the Isolates of Broad Bean Mottle Virus and Encapsidation of Host RNAs.

Broad bean mottle bromovirus infects legume plants and is transmissible by insects. Several broad bean mottle virus (BBMV) isolates have been identified, including one in England (isolate Ba) and five in the Mediterranean countries: Libya (LyV), Morocco (MV), Syria (SV), Sudan (TU) and Tunisia (TV). Previously, we analyzed the nucleotide sequence of the Ba RNA and here we report on and compare it with another five Mediterranean variants. The RNA segments in the latter ones were extensively homologous, with some SNPs, single nucleotide deletions and insertions, while the number of mutations was higher in isolate Ba. Both the 5' and 3' untranslated terminal regions (UTRs) among the corresponding RNAs are highly conserved, reflecting their functionality in virus replication. The AUG initiation codons are within suboptimal contexts, possibly to adjust/regulate translation. The proteins 1a, 2a, 3a and coat protein (CP) are almost identical among the five isolates, but in Ba they have more amino acid (aa) substitutions. Phylogenetic analysis revealed the isolates from Morocco and Syria clustering with the isolate from England, while the variants from Libya, Tunisia and Sudan created a different clade. The BBMV isolates encapsidate a high content of host (ribosomal and messenger) RNAs. Our studies present BBMV as a useful model for bromoviruses infecting legumes.

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.

A 313 plastome phylogenomic analysis of Pooideae: Exploring relationships among the largest subfamily of grasses.

In this study, we analyzed 313 plastid genomes (plastomes) of Poaceae with a focus on expanding our current knowledge of relationships among the subfamily Pooideae, which represented over half the dataset (164 representatives). In total, 47 plastomes were sequenced and assembled for this study. This is the largest study of its kind to include plastome-level data, to not only increase sampling at both the taxonomic and molecular levels with the aim of resolving complex and reticulate relationships, but also to analyze the effects of alignment gaps in large-scale analyses, as well as explore divergences in the subfamily with an expanded set of 14 accepted grass fossils for more accurate calibrations and dating. Incorporating broad systematic assessments of Pooideae taxa conducted by authors within the last five years, we produced a robust phylogenomic reconstruction for the subfamily, which included all but two supergeneric taxa (Calothecinae and Duthieeae). We further explored how including alignment gaps in plastome analyses oftentimes can produce incorrect or misinterpretations of complex or reticulate relationships among taxa of Pooideae. This presented itself as consistently changing relationships at specific nodes for different stripping thresholds (percentage-based removal of gaps per alignment column). Our summary recommendation for large-scale genomic plastome datasets is to strip alignment columns of all gaps to increase pairwise identity and reduce errant signal from poly A/T bias. To do this we used the "mask alignment" tool in Geneious software. Finally, we determined an overall divergence age for Pooideae of roughly 84.8 Mya, which is in line with, but slightly older than most recent estimates.

Transcriptional responses to water stress and recovery in a drought-tolerant fescue wild grass (Festuca ovina; Poaceae).

Water stress associated with drought-like conditions is a major factor limiting plant growth and impacts productivity of natural plant communities and agricultural crops. Molecular responses of plants to water stress have been studied most extensively in model species and crops, few of which have evolved natural drought tolerance. In the current study, we examined physiological and transcriptomic responses at multiple timepoints during increasing water stress and following initial recovery from stress in a drought-tolerant C3 species, Festuca ovina. Results demonstrated non-linear transcriptomic changes during increasing stress, but largely linear declines in physiological measurements during this same period. Transcription factors represented approximately 12.7% of all differentially expressed genes. In total, 117 F. ovina homologs of previously identified and molecularly characterized drought-responsive plant genes were identified. This information will be valuable for further investigations of the molecular mechanisms involved in drought tolerance in C3 plants.

Plastome phylogenomics and phylogenetic diversity of endangered and threatened grassland species (Poaceae) in a North American tallgrass prairie.

Native grasslands are one of the most endangered ecosystems in North America. In this study, we examined the ecological and evolutionary roles of endangered and threatened (e/t) grasses by establishing robust evolutionary relationships with other nonthreatened native and introduced grass species of the community. We hypothesized that the phylogenomic distribution of e/t species of grasses in Illinois would be phylogenetically clustered because closely related species would be vulnerable to the same threats and have similar requirements for survival. This study presents the first time a phylogeny based on complete plastome DNA of Poaceae was analyzed by phylogenetic diversity analysis. To avoid the disturbance of e/t populations, DNA was extracted from herbarium specimens. Next-generation sequencing (NGS) techniques were used to sequence DNA of plastid genomes (plastomes). The resulting phylogenomic tree was analyzed by phylogenetic diversity metrics. The extracted DNA successfully produced complete plastomes demonstrating that herbarium material is a practical source of DNA for genomic studies. The phylogenomic tree was strongly supported and defined Dichanthelium as a separate clade from Panicum. The phylogenetic metrics revealed phylogenetic clustering of e/t species, confirming our hypothesis.

Zygnema circumcarinatum UTEX 1559 chloroplast and mitochondrial genomes provide insight into land plant evolution.

The complete chloroplast and mitochondrial genomes of Charophyta have shed new light on land plant terrestrialization. Here, we report the organellar genomes of the Zygnema circumcarinatum strain UTEX 1559, and a comparative genomics investigation of 33 plastomes and 18 mitogenomes of Chlorophyta, Charophyta (including UTEX 1559 and its conspecific relative SAG 698-1a), and Embryophyta. Gene presence/absence was determined across these plastomes and mitogenomes. A comparison between the plastomes of UTEX 1559 (157 548 bp) and SAG 698-1a (165 372 bp) revealed very similar gene contents, but substantial genome rearrangements. Surprisingly, the two plastomes share only 85.69% nucleotide sequence identity. The UTEX 1559 mitogenome size is 215 954 bp, the largest among all sequenced Charophyta. Interestingly, this large mitogenome contains a 50 kb region without homology to any other organellar genomes, which is flanked by two 86 bp direct repeats and contains 15 ORFs. These ORFs have significant homology to proteins from bacteria and plants with functions such as primase, RNA polymerase, and DNA polymerase. We conclude that (i) the previously published SAG 698-1a plastome is probably from a different Zygnema species, and (ii) the 50 kb region in the UTEX 1559 mitogenome might be recently acquired as a mobile element.

Continued Adaptation of C4 Photosynthesis After an Initial Burst of Changes in the Andropogoneae Grasses.

C$_{4}$ photosynthesis is a complex trait that sustains fast growth and high productivity in tropical and subtropical conditions and evolved repeatedly in flowering plants. One of the major C$_{4}$ lineages is Andropogoneae, a group of $\sim $1200 grass species that includes some of the world's most important crops and species dominating tropical and some temperate grasslands. Previous efforts to understand C$_{4}$ evolution in the group have compared a few model C$_{4}$ plants to distantly related C$_{3}$ species so that changes directly responsible for the transition to C$_{4}$ could not be distinguished from those that preceded or followed it. In this study, we analyze the genomes of 66 grass species, capturing the earliest diversification within Andropogoneae as well as their C$_{3}$ relatives. Phylogenomics combined with molecular dating and analyses of protein evolution show that many changes linked to the evolution of C$_{4}$ photosynthesis in Andropogoneae happened in the Early Miocene, between 21 and 18 Ma, after the split from its C$_{3}$ sister lineage, and before the diversification of the group. This initial burst of changes was followed by an extended period of modifications to leaf anatomy and biochemistry during the diversification of Andropogoneae, so that a single C$_{4}$ origin gave birth to a diversity of C$_{4}$ phenotypes during 18 million years of speciation events and migration across geographic and ecological spaces. Our comprehensive approach and broad sampling of the diversity in the group reveals that one key transition can lead to a plethora of phenotypes following sustained adaptation of the ancestral state. [Adaptive evolution; complex traits; herbarium genomics; Jansenelleae; leaf anatomy; Poaceae; phylogenomics.].

Cell Wall Enzymes in Zygnema circumcarinatum UTEX 1559 Respond to Osmotic Stress in a Plant-Like Fashion.

Previous analysis of charophyte green algal (CGA) genomes and transcriptomes for specific protein families revealed that numerous land plant characteristics had already evolved in CGA. In this study, we have sequenced and assembled the transcriptome of Zygnema circumcarinatum UTEX 1559, and combined its predicted protein sequences with those of 13 additional species [five embryophytes (Emb), eight charophytes (Cha), and two chlorophytes (Chl) as the outgroup] for a comprehensive comparative genomics analysis. In total 25,485 orthologous gene clusters (OGCs, equivalent to protein families) of the 14 species were classified into nine OGC groups. For example, the Cha+Emb group contains 4,174 OGCs found in both Cha and Emb but not Chl species, representing protein families that have evolved in the common ancestor of Cha and Emb. Different OGC groups were subjected to a Gene Ontology (GO) enrichment analysis with the Chl+Cha+Emb group (including 5,031 OGCs found in Chl and Cha and Emb) as the control. Interestingly, nine of the 20 top enriched GO terms in the Cha+Emb group are cell wall-related, such as biological processes involving celluloses, pectins, lignins, and xyloglucans. Furthermore, three glycosyltransferase families (GT2, 8, 43) were selected for in-depth phylogenetic analyses, which confirmed their presence in UTEX 1559. More importantly, of different CGA groups, only Zygnematophyceae has land plant cellulose synthase (CesA) orthologs, while other charophyte CesAs form a CGA-specific CesA-like (Csl) subfamily (likely also carries cellulose synthesis activity). Quantitative real-time-PCR experiments were performed on selected GT family genes in UTEX 1559. After osmotic stress treatment, significantly elevated expression was found for GT2 family genes ZcCesA, ZcCslC and ZcCslA-like (possibly mannan and xyloglucan synthases, respectively), as well as for GT8 family genes (possibly pectin synthases). All these suggest that the UTEX 1559 cell wall polysaccharide synthesis-related genes respond to osmotic stress in a manner that is similar to land plants.

Plastome phylogenomics and characterization of rare genomic changes as taxonomic markers in plastome groups 1 and 2 Poeae (Pooideae; Poaceae).

A phylogenomic analysis of 42 complete plastid genomes (plastomes), including 16 that were newly sequenced, was conducted. Plastomes were sampled from 19 subtribes of Pooideae, to investigate relationships within and between Chloroplast Group 1 (Aveneae) and Group 2 (Poeae) species. Two data partitions: complete plastomes, and a combined plastome and rare genomic change (RGC) data matrix, were analyzed. Overall, 156 non-ambiguous RGC were identified, of which homology was inferred for 38 RGC. Among the 38 RGC identified, six were synapomorphic among the Group 1 subtribes: Aveninae, Agrostidinae, and Anthoxanthinae, (Phalaridinae + Torreyochloinae), and 27 were synapomorphic among the Group 2 subtribes: Loliinae, (Ammochloinae + Parapholiinae + Dactylidinae), Parapholiinae, Dactylidinae, Poinae, and Coleanthinae. Four RGC were determined to be homoplasious in Groups 1 and 2. Two other RGC originated through intrastrand deletion events. The remaining RGC events likely originated through recombination given their size and lack of sequence evidence for other types of mutations. This study also determined that relationships between taxa, even those only weakly supported in previous studies, could be inferred with strong support when utilizing complete plastomes.

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.

Investigation of mitochondrial-derived plastome sequences in the Paspalum lineage (Panicoideae; Poaceae).

BACKGROUND: The grass family (Poaceae), ca. 12,075 species, is a focal point of many recent studies that aim to use complete plastomes to reveal and strengthen relationships within the family. The use of Next Generation Sequencing technology has revealed intricate details in many Poaceae plastomes; specifically the trnI - trnL intergenic spacer region. This study investigates this region and the putative mitochondrial inserts within it in complete plastomes of Paspalum and other Poaceae. RESULTS: Nine newly sequenced plastomes, seven of which contain an insert within the trnI - trnL intergenic spacer, were combined into plastome phylogenomic and divergence date analyses with 52 other species. A robust Paspalum topology was recovered, originating at 10.6 Ma, with the insert arising at 8.7 Ma. The alignment of the insert across Paspalum reveals 21 subregions with pairwise homology in 19. In an analysis of emergent self-organizing maps of tetranucleotide frequencies, the Paspalum insert grouped with mitochondrial DNA. CONCLUSIONS: A hypothetical ancestral insert, 17,685 bp in size, was found in the trnI - trnL intergenic spacer for the Paspalum lineage. A different insert, 2808 bp, was found in the same region for Paraneurachne muelleri. Seven different intrastrand deletion events were found within the Paspalum lineage, suggesting selective pressures to remove large portions of noncoding DNA. Finally, a tetranucleotide frequency analysis was used to determine that the origin of the insert in the Paspalum lineage is mitochondrial DNA.

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.

Microsatellite and major histocompatibility complex variation in an endangered rattlesnake, the Eastern Massasauga (Sistrurus catenatus).

Genetic diversity is fundamental to maintaining the long-term viability of populations, yet reduced genetic variation is often associated with small, isolated populations. To examine the relationship between demography and genetic variation, variation at hypervariable loci (e.g., microsatellite DNA loci) is often measured. However, these loci are selectively neutral (or near neutral) and may not accurately reflect genomewide variation. Variation at functional trait loci, such as the major histocompatibility complex (MHC), can provide a better assessment of adaptive genetic variation in fragmented populations. We compared patterns of microsatellite and MHC variation across three Eastern Massasauga (Sistrurus catenatus) populations representing a gradient of demographic histories to assess the relative roles of natural selection and genetic drift. Using 454 deep amplicon sequencing, we identified 24 putatively functional MHC IIB exon 2 alleles belonging to a minimum of six loci. Analysis of synonymous and nonsynonymous substitution rates provided evidence of historical positive selection at the nucleotide level, and Tajima's D provided support for balancing selection in each population. As predicted, estimates of microsatellite allelic richness, observed, heterozygosity, and expected heterozygosity varied among populations in a pattern qualitatively consistent with demographic history and abundance. While MHC allelic richness at the population and individual levels revealed similar trends, MHC nucleotide diversity was unexpectedly high in the smallest population. Overall, these results suggest that genetic variation in the Eastern Massasauga populations in Illinois has been shaped by multiple evolutionary mechanisms. Thus, conservation efforts should consider both neutral and functional genetic variation when managing captive and wild Eastern Massasauga populations.

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.