Caught in the Act: Incipient Speciation at the Southern limit of Viburnum in the Central Andes.

A fundamental objective of evolutionary biology is to understand the origin of independently evolving species. Phylogenetic studies of species radiations rarely are able to document ongoing speciation; instead, modes of speciation, entailing geographic separation and/or ecological differentiation, are posited retrospectively. The Oreinotinus clade of Viburnum has radiated recently from north to south through the cloud forests of Mexico and Central America to the Central Andes. Our analyses support a hypothesis of incipient speciation in Oreinotinus at the southern edge of its geographic range, from central Peru to northern Argentina. Although several species and infraspecific taxa of have been recognized in this area, multiple lines of evidence and analytical approaches (including analyses of phylogenetic relationships, genetic structure, leaf morphology, and climatic envelopes) favor the recognition of just a single species, V. seemenii. We show that what has previously been recognized as V. seemenii f. minor has recently occupied the drier Tucuman-Bolivian forest region from Samaipata in Bolivia to Salta in northern Argentina. Plants in these populations form a well-supported clade with a distinctive genetic signature and they have evolved smaller, narrower leaves. We interpret this as the beginning of a within-species divergence process that has elsewhere in the neotropics resulted repeatedly in Viburnum species with a particular set of leaf ecomorphs. Specifically, the southern populations are in the process of evolving the small, glabrous, and entire leaf ecomorph that has evolved in four other montane areas of endemism. As predicted based on our studies of leaf ecomorphs in Chiapas, Mexico, these southern populations experience generally drier conditions, with large diurnal temperature fluctuations. In a central portion of the range of V. seemenii, characterized by wetter climatic conditions, we also document what may be the initial differentiation of the leaf ecomorph with larger, pubescent, and toothy leaves. The emergence of these ecomorphs thus appears to be driven by adaptation to subtly different climatic conditions in separate geographic regions, as opposed to parapatric differentiation along elevational gradients as suggested by Viburnum species distributions in other parts of the neotropics.

Flora diversity survey and establishment of a plant DNA barcode database of Lomas ecosystems in Peru.

Lomas formations or "fog oases" are islands of vegetation in the desert belt of the west coast of South America, with a unique vegetation composition among the world's deserts. However, plant diversity and conservation studies have long been neglected, and there exists a severe gap in plant DNA sequence information. To address the lack of DNA information, we conducted field collections and laboratory DNA sequencing to establish a DNA barcode reference library of Lomas plants from Peru. This database provides 1,207 plant specimens and 3,129 DNA barcodes data corresponding with collections from 16 Lomas locations in Peru, during 2017 and 2018. This database will facilitate both rapid species identification and basic studies on plant diversity, thereby enhancing our understanding of Lomas flora's composition and temporal variation, and providing valuable resources for conserving plant diversity and maintaining the stability of the fragile Lomas ecosystems.

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.

Genetic diversity and genetic structure of Puya raimondii (Bromeliaceae) for its conservation in the Peruvian Andes

Resumen Puya raimondii es una especie endémica de los altos Andes de Perú y Bolivia. En el Perú se distribuye desde 8.068501°S, 16.170280°W hasta 16.180580° S, 70.658873° W, entre los 3600 y 4800 m de altitud, viviendo en condiciones climáticas extremas propias de la Puna, donde juega un papel ecológico importante. Pese a la amplia distribución de las poblaciones de P. raimondii en el Perú, aparentemente son bastante uniformes morfológicamente; por lo que surgen las siguientes preguntas: ¿Podrán las actuales herramientas moleculares mostrar diferencias entre las numerosas poblaciones? ¿Son suficientes las áreas de conservación establecidas para P. raimondii ya que albergan la variabilidad existente? Para responder a estas interrogantes, este trabajo tuvo como objetivo evaluar la diversidad genética y estructura genética en una población del norte del país, Pachapaqui (departamento de Ancash), una población del centro, Yanacancha (Junín), y una población del sur, Lampa - sector Choconchaca (Puno), utilizando marcadores microsatélites (SSR) específicos para la especie. Los parámetros de diversidad genética utilizados incluyeron número de alelos (A), alelos exclusivos (RA), heterocigosidad observada (Ho), heterocigosidad esperada (He) e índice de contenido polimórfico (PIC). Los resultados mostraron que el número total de A varió de 2 ‒ 13, los valores de He fueron 0 ‒ 0.723 y Ho 0 ‒ 0.929, con un He promedio de 0.217, indicando una diversidad genética moderada a alta, siendo la población de Lampa-sector Choconchaca, la que presentó mayor diversidad alélica y mayor diversidad genética. La prueba de Hardy-Weinberg mostró que las poblaciones se encuentran en desequilibrio HW, el análisis estadístico indica un 65% de variación genética a nivel poblacional y valores de FST (0.426) y RST (0.650) que indican alta diferenciación genética entre poblaciones, con dos grupos genéticos (K=2) que corresponden a las poblaciones del centro-norte y sur del Perú. Los resultados brindan información útil para establecer estrategias de conservación para P. raimondii, que conduzcan a la creación de una área de conservación adicional para proteger a las poblaciones del sur del Perú.

Abstract Puya raimondii is an endemic species from the high Andes of Peru and Bolivia. In Peru it is distributed from 8.068501°S, 16.170280°W to 16.180580°S, 70.658873°W, between 3600 and 4800 m, living in extreme climatic conditions typical of the Puna, where it plays an important ecological role. Despite the wide distribution of P. raimondii populations in Peru, they appear to be fairly uniform morphologically. The following questions arise: Will the current molecular tools be able to show differences between the numerous populations? Are the conservation areas established for P. raimondii sufficient since they harbor the existing variability? To answer these questions, this work aimed to evaluate the genetic diversity and genetic structure in a northern population, Pachapaqui (Ancash department), a central population, Yanacancha (Junin), and a southern population, Lampa - Choconchaca sector (Puno), using microsatellite markers (SSR) specific for the species. The genetic diversity parameters used included number of alleles (A), exclusive alleles (RA), observed heterozygosity (Ho), expected heterozygosity (He), and polymorphic content index (PIC). The results showed that the total number of A varied from 2 - 13, the He values were 0 ‒ 0.723 and Ho 0 ‒ 0.929, with an average He of 0.217, indicating a moderate to high genetic diversity, being the population of Lampa-Choconchaca sector, the one that presented the greatest allelic diversity and the greatest genetic diversity. The Hardy-Weinberg test showed that the populations are in HW disequilibrium, the statistical analysis indicates 65% of the genetic variation at the population level and values of FST (0.426) and RST (0.650) that indicate high genetic differentiation among populations, with two genetic groups (K=2) that correspond to the populations of northern-central and southern Peru. The results provide useful information to establish conservation strategies for P. raimondii, which lead to the creation of an additional conservation area to protect the populations in southern Peru.

Draft genome of Puya raimondii (Bromeliaceae), the Queen of the Andes.

Puya raimondii, the Queen of the Andes, is an endangered high Andean species in the Bromeliaceae family. Here, we report its first genome to promote its conservation and evolutionary study. Comparative genomics showed P. raimondii diverged from Ananas comosus about 14.8 million years ago, and the long terminal repeats were likely to contribute to the genus diversification in last 3.5 million years. The gene families related to plant reproductive development and stress responses significantly expanded in the genome. At the same time, gene families involved in disease defense, photosynthesis and carbohydrate metabolism significantly contracted, which may be an evolutionary strategy to adapt to the harsh conditions in high Andes. The demographic history analysis revealed the P. raimondii population size sharply declined in the Pleistocene and then increased in the Holocene. We also designed and tested 46 pairs of universal primers for amplifying orthologous single-copy nuclear genes in Puya species.

Quaternary diversification of a columnar cactus in the driest place on earth.

PREMISE: The cactus family (Cactaceae) is a speciose lineage with an almost entirely New World distribution. The genus Eulychnia with eight currently recognized species is endemic to the Atacama and Peruvian Deserts. Here we investigated the phylogeny of this group based on a complete taxon sampling to elucidate species delimitation and biogeographic history of the genus. METHODS: A family-wide Bayesian molecular clock dating based on plastid sequence data was conducted to estimate the age of Eulychnia and its divergence from its sister genus Austrocactus. A second data set obtained from genotyping by sequencing (GBS) was analyzed, using the family-wide age estimate as a secondary calibration to date the GBS phylogeny using a penalized likelihood approach. Ancestral ranges were inferred employing the dispersal extinction cladogenesis approach. RESULTS: Our GBS phylogeny of Eulychnia was fully resolved with high support values nearly throughout the phylogeny. The split from Austrocactus occurred in the late Miocene, and Eulychnia diversified during the early Quaternary. Three lineages were retrieved: Eulychnia ritteri from Peru is sister to all Chilean species, which in turn fall into two sister clades of three and four species, respectively. Diversification in the Chilean clades started in the early Pleistocene. Eulychnia likely originated at the coastal range of its distribution and colonized inland locations several times. CONCLUSIONS: Diversification of Eulychnia during the Pleistocene coincides with long periods of hyperaridity alternated with pluvial phases. Hyperaridity caused habitat fragmentation, ultimately leading to speciation and resulting in the current allopatric distribution of taxa.

Cinchona L. "Árbol de la Quina": repoblamiento y reforestación en el Perú / Cinchona L. "Quina Tree": Repopulation and reforestation in Peru

Resumen Este trabajo se centra en las recientes decisiones tomadas por el Ministerio de Agricultura del Perú sobre la recolección, propagación y siembra de especies de quina. Señalamos las deficiencias de estas decisiones a la luz de los avances en el conocimiento de la taxonomía y sistemática de Cinchona. Destacamos el papel de la ciencia analizando e informando sobre las decisiones en el uso sostenible de los recursos naturales.

Abstract This work focuses about the recent decisions made by the Ministry of Agriculture of Peru on the collection, propagation, and sowing of Cinchona species are reasons to point out the deficiencies of these decisions in light of advances in the knowledge of Cinchona taxonomy and systematics. We emphasize the role of the science analysing and informing about decisions in the sustainable use of natural resources.

Characterization of 16 microsatellite markers for the Oreinotinus clade of Viburnum (Adoxaceae).

PREMISE OF THE STUDY: Microsatellite loci were isolated from four species of Viburnum (Adoxaceae) to study population structure and assess species boundaries among morphologically similar South American Viburnum species of the Oreinotinus clade. METHODS AND RESULTS: Using a microsatellite-enriched library and mining next-generation sequence data, 16 microsatellites were developed. Each locus was tested on two populations of V. triphyllum and one population of V. pichinchense. For nuclear loci, one to 13 alleles were recovered, expected heterozygosity ranged from 0 to 0.8975, Simpson diversity index ranged from 0.0167 to 1.000, and Shannon diversity index ranged from 0 to 2.3670 in a given population. For the mitochondrial locus, three to six alleles were recovered and unbiased haploid diversity values ranged from 0.756 to 0.853 in a given population. CONCLUSIONS: The 16 microsatellite loci developed for the Oreinotinus clade (Viburnum, Adoxaceae) will inform investigations of population structure and species boundaries within this group.

Genetic enablers underlying the clustered evolutionary origins of C4 photosynthesis in angiosperms.

The evolutionary accessibility of novel adaptations varies among lineages, depending in part on the genetic elements present in each group. However, the factors determining the evolutionary potential of closely related genes remain largely unknown. In plants, CO2-concentrating mechanisms such as C4 and crassulacean acid metabolism (CAM) photosynthesis have evolved numerous times in distantly related groups of species, and constitute excellent systems to study constraints and enablers of evolution. It has been previously shown for multiple proteins that grasses preferentially co-opted the same gene lineage for C4 photosynthesis, when multiple copies were present. In this work, we use comparative transcriptomics to show that this bias also exists within Caryophyllales, a distantly related group with multiple C4 origins. However, the bias is not the same as in grasses and, when all angiosperms are considered jointly, the number of distinct gene lineages co-opted is not smaller than that expected by chance. These results show that most gene lineages present in the common ancestor of monocots and eudicots produced gene descendants that were recruited into C4 photosynthesis, but that C4-suitability changed during the diversification of angiosperms. When selective pressures drove C4 evolution, some copies were preferentially co-opted, probably because they already possessed C4-like expression patterns. However, the identity of these C4-suitable genes varies among clades of angiosperms, and C4 phenotypes in distant angiosperm groups thus represent genuinely independent realizations, based on different genetic precursors.

A chloroplast tree for Viburnum (Adoxaceae) and its implications for phylogenetic classification and character evolution.

• Premise of the study: Despite recent progress, significant uncertainties remain concerning relationships among early-branching lineages within Viburnum (Adoxaceae), prohibiting a new classification and hindering studies of character evolution and the increasing use of Viburnum in addressing a wide range of ecological and evolutionary questions. We hoped to resolve these issues by sequencing whole plastid genomes for representative species and combining these with molecular data previously obtained from an expanded taxon sample.• Methods: We performed paired-end Illumina sequencing of plastid genomes of 22 Viburnum species and combined these data with a 10-gene data set to infer phylogenetic relationships for 113 species. We used the results to devise a comprehensive phylogenetic classification and to analyze the evolution of eight morphological characters that vary among early-branching lineages.• Key results: With greatly increased levels of confidence in most of the early branches, we propose a phylogenetic classification of Viburnum, providing formal phylogenetic definitions for 30 clades, including 13 with names recognized under the International Code of Nomenclature for Algae, Fungi, and Plants, eight with previously proposed informal names, and nine newly proposed names for major branches. Our parsimony reconstructions of bud structure, leaf margins, inflorescence form, ruminate endosperm, extrafloral nectaries, glandular trichomes, palisade anatomy, and pollen exine showed varying levels of homoplasy, but collectively provided morphological support for some, though not all, of the major clades.• Conclusions: Our study demonstrates the value of next-generation plastid sequencing, the ease of creating a formal phylogenetic classification, and the utility of such a system in describing patterns of character evolution.

Shared origins of a key enzyme during the evolution of C4 and CAM metabolism.

CAM and C4 photosynthesis are two key plant adaptations that have evolved independently multiple times, and are especially prevalent in particular groups of plants, including the Caryophyllales. We investigate the origin of photosynthetic PEPC, a key enzyme of both the CAM and C4 pathways. We combine phylogenetic analyses of genes encoding PEPC with analyses of RNA sequence data of Portulaca, the only plants known to perform both CAM and C4 photosynthesis. Three distinct gene lineages encoding PEPC exist in eudicots (namely ppc-1E1, ppc-1E2 and ppc-2), one of which (ppc-1E1) was recurrently recruited for use in both CAM and C4 photosynthesis within the Caryophyllales. This gene is present in multiple copies in the cacti and relatives, including Portulaca. The PEPC involved in the CAM and C4 cycles of Portulaca are encoded by closely related yet distinct genes. The CAM-specific gene is similar to genes from related CAM taxa, suggesting that CAM has evolved before C4 in these species. The similar origin of PEPC and other genes involved in the CAM and C4 cycles highlights the shared early steps of evolutionary trajectories towards CAM and C4, which probably diverged irreversibly only during the optimization of CAM and C4 phenotypes.

Genetic variability of an unusual apomictic triploid cactus--Haageocereus tenuis Ritter--from the Coast of Central Peru.

Haageocereus tenuis is a prostrate cactus restricted to a small area of 2 km(2) near the city of Lima, Peru. The species is triploid and propagates mainly through stem fragmentation. In addition, propagation via agamospermy is documented and adventitious embryony is also inferred as a mechanism. Although seedling recruitment has not been observed in nature, we have shown that asexually produced seeds are viable. About 45 adult individuals, plus 9 individuals obtained from seeds, were sampled and 5 microsatellite markers were used to assess genetic variability. Microsatellite analysis confirms that individuals from the only existing population are genetically identical and that the population likely represents a single clone. The absence of mutations in any individual, even in highly variable microsatellite loci, may indicate that the species is also of recent origin. Other prostrate species of Haageocereus are suspected to be occasional apomicts. This phenomenon has significant implications for the evolutionary biology and ecology of Haageocereus and other clonal Cactaceae.

Contemporaneous and recent radiations of the world's major succulent plant lineages.

The cacti are one of the most celebrated radiations of succulent plants. There has been much speculation about their age, but progress in dating cactus origins has been hindered by the lack of fossil data for cacti or their close relatives. Using a hybrid phylogenomic approach, we estimated that the cactus lineage diverged from its closest relatives ≈35 million years ago (Ma). However, major diversification events in cacti were more recent, with most species-rich clades originating in the late Miocene, ≈10-5 Ma. Diversification rates of several cactus lineages rival other estimates of extremely rapid speciation in plants. Major cactus radiations were contemporaneous with those of South African ice plants and North American agaves, revealing a simultaneous diversification of several of the world's major succulent plant lineages across multiple continents. This short geological time period also harbored the majority of origins of C(4) photosynthesis and the global rise of C(4) grasslands. A global expansion of arid environments during this time could have provided new ecological opportunity for both succulent and C(4) plant syndromes. Alternatively, recent work has identified a substantial decline in atmospheric CO(2) ≈15-8 Ma, which would have strongly favored C(4) evolution and expansion of C(4)-dominated grasslands. Lowered atmospheric CO(2) would also substantially exacerbate plant water stress in marginally arid environments, providing preadapted succulent plants with a sharp advantage in a broader set of ecological conditions and promoting their rapid diversification across the landscape.

C(4) eudicots are not younger than C(4) monocots.

C(4) photosynthesis is a plant adaptation to high levels of photorespiration. Physiological models predict that atmospheric CO(2) concentration selected for C(4) grasses only after it dropped below a critical threshold during the Oligocene (∼30 Ma), a hypothesis supported by phylogenetic and molecular dating analyses. However the same models predict that CO(2) should have reached much lower levels before selecting for C(4) eudicots, making C(4) eudicots younger than C(4) grasses. In this study, different phylogenetic datasets were combined in order to conduct the first comparative analysis of the age of C(4) origins in eudicots. Our results suggested that all lineages of C(4) eudicots arose during the last 30 million years, with the earliest before 22 Ma in Chenopodiaceae and Aizoaceae, and the latest probably after 2 Ma in Flaveria. C(4) eudicots are thus not globally younger than C(4) monocots. All lineages of C(4) plants evolved in a similar low CO(2) atmosphere that predominated during the last 30 million years. Independent C(4) origins were probably driven by different combinations of specific factors, including local ecological characteristics such as habitat openness, aridity, and salinity, as well as the speciation and dispersal history of each clade. Neither the lower number of C(4) species nor the frequency of C(3)-C(4) intermediates in eudicots can be attributed to a more recent origin, but probably result from variation in diversification and evolutionary rates among the different groups that evolved the C(4) pathway.

Isolation and characterization of novel microsatellite markers for Arctium minus (Compositae).

PREMISE OF THE STUDY: Microsatellite primers were developed for the invasive plant Arctium minus to investigate the effects of facultative self-pollination and the biannual habit on population genetic structure, as well as the colonization of the Americas by this Eurasian species. • METHODS AND RESULTS: Sixteen di- and trinucleotide microsatellite loci were identified in six populations. The number of alleles per locus ranged from one to 10, observed heterozygosities ranged from 0.000 to 0.897, and the mean value of F(IS) was 0.316. • CONCLUSIONS: These results indicate the utility of these loci in future studies of population genetics in A. minus.

Characterization of polymorphic microsatellite loci in Haageocereus (Trichocereeae, Cactaceae).

PREMISE OF THE STUDY: Microsatellites were isolated from two species of the genus Haageocereus (H. tenuis and H. pseudomelanostele) to be applied in studies of genetic diversity and population structure. • METHODS AND RESULTS: Five loci were employed in a preliminary study of genetic diversity and population differentiation in two rare (H. tenuis and H. repens) and two widespread (H. acranthus and H. pseudomelanostele) species, yielding between one and 44 alleles per locus. All five loci were polymorphic, with overall levels of observed and expected heterozygosities ranging from 0.478 to 0.871 and from 0.564 to 0.956, respectively. Three additional loci were scored in H. pseudomelanostele. These eight plus the remaining 11 loci were amplified from putative parents of three hybrids involving Haageocereus and Espostoa. • CONCLUSIONS: These markers will facilitate analysis of genetic diversity, hybridization, and population differentiation throughout Haageocereus and Espostoa.

Composición florística de la cuenca del río Ilo-Moquegua y Lomas de Ilo, Moquegua, Perú / Floristic composition of river Ilo-Moquegua basin and Lomas of Ilo, Moquegua, Peru

Se presenta un estudio de la composición y riqueza de la flora vascular en la cuenca del Río Ilo-Moquegua y las Lomas de Ilo, a lo largo de una gradiente altitudinal que va desde los 0 hasta 4600m basado en colectas intensivas, revisión de ejemplares de herbario y consulta bibliográfica. De acuerdo al presente estudio, la flora vascular de la cuenca del río Ilo: Moquegua y las Lomas de Ilo esta integrada por 63 familias, 233 géneros y 394 especies. Las magnoliopsida representan el 83 por ciento de las especies y las liliopsida el 15 por ciento. Las familias con mayor número de géneros y especies son Asteraceae (41 géneros y 60 especies), Poaceace (28 y 44), Solanaceae 811 y 329, Fabaceae (17 y 26), malvaceae (11 y 21) Brassicaceae (10 y 15) Boriginaceae (9 y 15) y Cactaceae 810 y 32), Fabacear 817 y 26), Malvaceae 811 y 21), Brasicaceae (10 y 15) Boraginaceae 89 y 15) y Cactaceae (10 y 14). Las formas biológicas dominantes son las hierbas (72 por ciento), seguidas por los arbustos (21 por ciento), plantas que pueden ser hierbas o arbustos (5 por ciento), árboles (2 por ciento) y plantas parásitas (menos del 1 por ciento). De las 55 especies endémicas del Perú reportadas en el presente estudio, 10 están restringidas al departamento de Moquegua. Se presenta un total de 176 nuevas adiciones a la flora del departamento.