The 'queen of the Andes' (Puya raimondii) is genetically fragile and fragmented: a consequence of long generation time and semelparity?

Understanding how life history shapes genetic diversity is a fundamental issue in evolutionary biology, with important consequences for conservation. However, we still have an incomplete picture of the impact of life history on genome-wide patterns of diversity, especially in long-lived semelparous plants. Puya raimondii is a high-altitude semelparous species from the Andes that flowers at 40-100 years of age. We sequenced the whole genome and estimated the nucleotide diversity of 200 individuals sampled from nine populations. Coalescent-based approaches were then used to infer past population dynamics. Finally, these results were compared with results obtained for the iteroparous species, Puya macrura. The nine populations of P. raimondii were highly divergent, highly inbred, and carried an exceptionally high genetic load. They are genetically depauperate, although, locally in the genome, balancing selection contributed to the maintenance of genetic polymorphism. While both P. raimondii and P. macrura went through a severe bottleneck during the Pleistocene, P. raimondii did not recover from it and continuously declined, while P. macrura managed to bounce back. Our results demonstrate the importance of life history, in particular generation time and reproductive strategy, in affecting population dynamics and genomic variation, and illustrate the genetic fragility of long-lived semelparous plants.

A plastid phylogenomic framework for the palm family (Arecaceae).

BACKGROUND: Over the past decade, phylogenomics has greatly advanced our knowledge of angiosperm evolution. However, phylogenomic studies of large angiosperm families with complete species or genus-level sampling are still lacking. The palms, Arecaceae, are a large family with ca. 181 genera and 2600 species and are important components of tropical rainforests bearing great cultural and economic significance. Taxonomy and phylogeny of the family have been extensively investigated by a series of molecular phylogenetic studies in the last two decades. Nevertheless, some phylogenetic relationships within the family are not yet well-resolved, especially at the tribal and generic levels, with consequent impacts for downstream research. RESULTS: Plastomes of 182 palm species representing 111 genera were newly sequenced. Combining these with previously published plastid DNA data, we were able to sample 98% of palm genera and conduct a plastid phylogenomic investigation of the family. Maximum likelihood analyses yielded a robustly supported phylogenetic hypothesis. Phylogenetic relationships among all five palm subfamilies and 28 tribes were well-resolved, and most inter-generic phylogenetic relationships were also resolved with strong support. CONCLUSIONS: The inclusion of nearly complete generic-level sampling coupled with nearly complete plastid genomes strengthened our understanding of plastid-based relationships of the palms. This comprehensive plastid genome dataset complements a growing body of nuclear genomic data. Together, these datasets form a novel phylogenomic baseline for the palms and an increasingly robust framework for future comparative biological studies of this exceptionally important plant family.

Genetic markers in Andean Puya species (Bromeliaceae) with implications on plastome evolution and phylogeny.

The Andean plant endemic Puya is a striking example of recent and rapid diversification from central Chile to the northern Andes, tracking mountain uplift. This study generated 12 complete plastomes representing nine Puya species and compared them to five published plastomes for their features, genomic evolution, and phylogeny. The total size of the Puya plastomes ranged from 159,542 to 159,839 bp with 37.3%-37.4% GC content. The Puya plastomes were highly conserved in organization and structure with a typical quadripartite genome structure. Each of the 17 consensus plastomes harbored 133 genes, including 87 protein-coding genes, 38 tRNA (transfer RNA) genes, and eight rRNA (ribosomal RNA) genes; we found 69-78 tandem repeats, 45-60 SSRs (simple sequence repeats), and 8-22 repeat structures among 13 species. Four protein-coding genes were identified under positive site-specific selection in Puya. The complete plastomes and hypervariable regions collectively provided pronounced species discrimination in Puya and a practical tool for future phylogenetic studies. The reconstructed phylogeny and estimated divergence time for the lineage suggest that the diversification of Puya is related to Andean orogeny and Pleistocene climatic oscillations. This study provides plastome resources for species delimitation and novel phylogenetic and biogeographic studies.

Fingerprints of climatic changes through the late Cenozoic in southern Asian flora: Magnolia section Michelia (Magnoliaceae).

BACKGROUND AND AIMS: Ongoing global warming is a challenge for humankind. A series of drastic climatic changes have been proven to have occurred throughout the Cenozoic based on a variety of geological evidence, which helps to better understand our planet's future climate. Notably, extant biomes have recorded drastic environmental shifts. The climate in southern Asia, which hosts high biodiversity, is deeply impacted by the Asian monsoon. The origins and evolutionary dynamics of biomes occurring between the tropics and sub-tropics in southern Asia have probably been deeply impacted by climatic changes; however, these aspects remain poorly studied. We tested whether the evolutionary dynamics of the above biomes have recorded the drastic, late Cenozoic environmental shifts, by focusing on Magnolia section Michelia of the family Magnoliaceae. METHODS: We established a fine time-calibrated phylogeny of M. section Michelia based on complete plastid genomes and inferred its ancestral ranges. Finally, we estimated the evolutionary dynamics of this section through time, determining its diversification rate and the dispersal events that occurred between tropical and sub-tropical areas. KEY RESULTS: The tropical origin of M. section Michelia was dated to the late Oligocene; however, the diversification of its core group (i.e. M. section Michelia subsection Michelia) has occurred mainly from the late Miocene onward. Two key evolutionary shifts (dated approx. 8 and approx. 3 million years ago, respectively) were identified, each of them probably in response to drastic climatic changes. CONCLUSION: Here, we inferred the underlying evolutionary dynamics of biomes in southern Asia, which probably reflect late Cenozoic climatic changes. The occurrence of modern Asian monsoons was probably fundamental for the origin of M. section Michelia; moreover, the occurrence of asymmetric dispersal events between the tropics and sub-tropics hint at an adaptation strategy of M. section Michelia to global cooling, in agreement with the tropical conservatism hypothesis.

The utility of DNA barcodes to confirm the identification of palm collections in botanical gardens.

The palm family (Arecaceae) is of high ecological and economic value, yet identification in the family remains a challenge for both taxonomists and horticulturalists. The family consists of approximately 2600 species across 181 genera and DNA barcoding may be a useful tool for species identification within the group. However, there have been few systematic evaluations of DNA barcodes for the palm family. In the present study, five DNA barcodes (rbcL, matK, trnH-psbA, ITS, ITS2) were evaluated for species identification ability across 669 samples representing 314 species and 100 genera in the Arecaceae, employing four analytical methods. The ITS gene region was found to not be a suitable barcode for the palm family, due in part, to low recovery rates and paralogous gene copies. Among the four analyses used, species resolution for ITS2 was much higher than that achieved with the plastid barcodes alone (rbcL, matK, trnH-psbA), and the barcode combination ITS2 + matK + rbcL gave the highest resolution among all single barcodes and their combinations, followed by ITS2 + matK. Among 669 palm samples analyzed, 110 samples (16.3%) were found to be misidentified. The 2992 DNA barcode sequences generated in this study greatly enriches the existing identification toolbox available to plant taxonomists that are interested in researching genetic relationships among palm taxa as well as for horticulturalists that need to confirm palm collections for botanical garden curation and horticultural applications. Our results indicate that the use of the ITS2 DNA barcode gene region provides a useful and cost-effective tool to confirm the identity of taxa in the Palm family.

Extensive reticulate evolution within Fargesia (s.l.) (Bambusoideae: Poaceae) and its allies: Evidence from multiple nuclear markers.

Reticulate evolution resulting from hybridization and introgression has been recognized as a creative source of species and diversification in bamboos. Previous phylogenetic studies revealed that Fargesia (s.l.) (Fargesia and Yushania) was divided into the Fargesia spathe clade and the non-spathe clade. Interestingly, the Fargesia spathe clade may have originated from hybridization among other clades within Fargesia (s.l.). Understanding the hybrid origin of this clade requires a robust phylogenetic framework in which major clades within Fargesia (s.l.) are resolved. Here, we used three nuclear genes to reconstruct the evolutionary history of Fargesia (s.l.) and its allies to identify putative patterns in the origin of the Fargesia spathe clade and to examine the extent to which reticulate evolution has occurred at the interspecific level in bamboos. Bashania species form a clade with Fargesia (s.l.), which is further divided into Group I and Group II. The Fargesia spathe clade, the Alpine Bashania clade, and Fargesia yajiangensis comprise Group I, while the Bashania fargesii clade and the remaining Fargesia (s.l.) species form Group II. Incongruence between the current nuclear-based and previous plastid phylogenies demonstrate several possible hybridization events among Fargesia (s.l.) species and related taxa, which have given rise to the Fargesia spathe clade, the Phyllostachys clade, and the Ampelocalamus clade. We also detected several putative hybrid species of Fargesia (s.l.). Our results show that reticulate evolution has played a prominent role in Fargesia (s.l.) evolution, which could, in part, account for the taxonomic difficulty associated with Fargesia (s.l.) and the alpine bamboos. The study also underscores the importance of hybridization in the evolution of bamboos, at both intergeneric and intrageneric levels.

Straight From the Plastome: Molecular Phylogeny and Morphological Evolution of Fargesia (Bambusoideae: Poaceae).

Fargesia is ecologically and economically important in mountainous forests. Many Fargesia species are also important sources of food for some endangered animals such as the giant panda. Recent molecular phylogenetic analyses have revealed Fargesia as a polyphyletic group despite some unclear lineage affinities. In the present study, we reconstructed the phylogeny of Fargesia and its allies, including Thamnocalamus, Arundinaria (incl. Bashania), Yushania, Indocalamus, Ampelocalamus and Phyllostachys, from a plastome sequence matrix that contained 20 Fargesia and five Yushania species as ingroups, 16 species from nine other bamboo genera plus Oryza sativa and Zea mays as outgroups. Fargesia and its allies were broken into eight clades. Several Fargesia species were clustered into the Thamnocalamus clade and the Drepanostachyum + Himalayacalamus clade that rendered the polyphyly of Fargesia. The remaining six clades, including the Fargesia spathe clade, the Phyllostachys clade, Arundinaria fargesii, the Ampelocalamus clade, the Fargesia grossa clade, and the Fargesia macclureana clade, were identified. Molecular phylogenetic analyses supported that Yushania should be included in Fargesia (s.l.) which had synapomorphy of expanded leaf sheaths in varying degree at the basis of inflorescences, and further divided into the Fargesia spathe clade, the Fargesia grossa clade, and the Fargesia macclureana clade. All sampled species of Yushania were nested within the Fargesia grossa clade. The probable model of the origin of the species in the Fargesia spathe clade with spathe-like leaf sheath syndrome was proposed. Moreover, the formation of the spathe-like leaf sheath syndrome may be correlated with cold climatic conditions in Quaternary. Our results provide new sight into the phylogenetic relationship within Fargesia.

Microsatellite markers for the endangered Puya raimondii in Peru.

PREMISE: Microsatellite primers were developed for Puya raimondii (Bromeliaceae), an endangered species distributed in the Andean Mountains of Bolivia and Peru. METHODS AND RESULTS: Genome skimming of P. raimondii, P. macrura, and P. hutchisonii resulted in the selection of 46 pairs of cross-species microsatellite markers. Of these, 12 microsatellite primer pairs produced clear and polymorphic bands in P. raimondii. These primer sets were then used for the detection of potential polymorphisms in 84 P. raimondii individuals collected from four populations in Peru. The number of alleles per locus ranged from one to six, and the observed and expected levels of heterozygosity ranged from 0.000 to 0.8929 and from 0.000 to 0.7662, respectively. CONCLUSIONS: The microsatellite markers developed in this study will be useful for future population genetic analyses and breeding system studies in P. raimondii.

Morphological variation in leaf dissection of Rheum palmatum complex (Polygonaceae).

AIMS: Rheum palmatum complex comprises all taxa within section Palmata in the genus Rheum, including R. officinale, R. palmatum, R. tanguticum, R. tanguticum var. liupanshanense and R. laciniatum. The identification of the taxa in section Palmata is based primarily on the degree of leaf blade dissection and the shape of the lobes; however, difficulties in species identification may arise from their significant variation. The aim of this study is to analyze the patterns of variation in leaf blade characteristics within and among populations through population-based sampling covering the entire distribution range of R. palmatum complex. METHODS: Samples were taken from 2340 leaves from 780 individuals and 44 populations representing the four species, and the degree of leaf blade dissection and the shape of the lobe were measured to yield a set of quantitative data. Furthermore, those data were statistically analyzed. IMPORTANT FINDINGS: The statistical analysis showed that the degree of leaf blade dissection is continuous from lobed to parted, and the shape of the lobe is also continuous from broadly triangular to lanceolate both within and between populations. We suggested that taxa in section Palmata should be considered as one species. Based on the research on the R. palmatum complex, we considered that the quantitative characteristics were greatly influenced by the environment. Therefore, it is not reliable to delimitate the species according to the continuously quantitative vegetative characteristics.

Merging Fargesia dracocephala into Fargesia decurvata (Bambusoideae, Poaceae): implications from morphological and ITS sequence analyses.

AIMS: Fargesia decurvata is closely allied with F. dracocephala and differs in 5 major characters (i.e. the culm sheath blade base shape, the width of the culm sheath blade base, the auricle shape, and the lower surface of leaf blade) in Fargesia. It is difficult to distinguish these two species because of existing of transitional statements of characters. The aims of this paper are to (i) investigate whether the variation of the characters is continuous or not; (ii) reveal whether the publishment of F. dracocephala was the result of discontinuous sampling of F. decurvata or not. METHODS: Ten populations of F. decurvata and F. dracocephala were investigated in their entire distribution (including type localities). The statements of 5 major characters were measured from 693 annual and 693 perennial culms of 231 individuals in 10 populations, and analyzed at population, individual and culm levels. UPGMA cluster analysis was carried out based on 29 characters from 10 populations of F. decurvata and F. dracocephala and 2 populations of F. qinlingensis as outgroup. The ITS sequences were also sequenced and analyzed. IMPORTANT FINDINGS: Five major characters exhibited great variation not only at population level, but at individual level within a population, even the culm level within an individual and in different parts of the same culm. Cluster analyses showed that 10 populations of F. decurvata and F. dracocephala were not divided into two species, but they were well separated with outgroup. There was no difference in floral organ between F. decurvata and F. dracocephala. MP and NJ trees based on ITS sequences showed the same results with the cluster analysis on morphological characters. All the facts indicated that the publishment of F. dracocephala was the result of discontinuous sampling of F. decurvata, and F. dracocephala should be treated as the synonym of F. decurvata.

Genetic diversity of the endemic and medicinally important plant Rheum officinale as revealed by Inter-Simpe Sequence Repeat (ISSR) Markers.

Rheum officinale Baill., an important but endangered medicinal herb, is endemic to China. Inter-simple sequence repeat (ISSR) markers were employed to investigate the genetic diversity and differentiation of 12 populations of R. officinale. Thirteen selected primers yielded 189 bright and discernible bands, with an average of 14.54 per primer. The genetic diversity was low at the population level, but pretty high at the species level (H = 0.1008, I = 0.1505, PPB = 28.95% vs. H = 0.3341, I = 0.5000, PPB = 95.24%, respectively) by POPGENE analysis. Analysis of molecular variance (AMOVA) showed that the genetic variation was found mainly among populations (74.38%), in line with the limited gene flow (N(m) = 0.2766) among populations. Mantel test revealed a significant correlation between genetic and geographic distances (r = 0.5381, P = 0.002), indicating the role of geographic isolation in shaping the present population genetic structure. Both Bayesian analysis and UPGMA cluster analysis demonstrated the similar results. Our results imply that the conservation efforts should aim to preserve all the extant populations of this endangered species, and cultivation is proposed in this study.