Underlying and proximate drivers of biodiversity changes in Mesoamerican biosphere reserves.

Protected areas are of paramount relevance to conserving wildlife and ecosystem contributions to people. Yet, their conservation success is increasingly threatened by human activities including habitat loss, climate change, pollution, and species overexploitation. Thus, understanding the underlying and proximate drivers of anthropogenic threats is urgently needed to improve protected areas' effectiveness, especially in the biodiversity-rich tropics. We addressed this issue by analyzing expert-provided data on long-term biodiversity change (last three decades) over 14 biosphere reserves from the Mesoamerican Biodiversity Hotspot. Using multivariate analyses and structural equation modeling, we tested the influence of major socioeconomic drivers (demographic, economic, and political factors), spatial indicators of human activities (agriculture expansion and road extension), and forest landscape modifications (forest loss and isolation) as drivers of biodiversity change. We uncovered a significant proliferation of disturbance-tolerant guilds and the loss or decline of disturbance-sensitive guilds within reserves causing a "winner and loser" species replacement over time. Guild change was directly related to forest spatial changes promoted by the expansion of agriculture and roads within reserves. High human population density and low nonfarming occupation were identified as the main underlying drivers of biodiversity change. Our findings suggest that to mitigate anthropogenic threats to biodiversity within biosphere reserves, fostering human population well-being via sustainable, nonfarming livelihood opportunities around reserves is imperative.

Population Genomics of Domesticated Cucurbita ficifolia Reveals a Recent Bottleneck and Low Gene Flow with Wild Relatives.

Cucurbita ficifolia is a squash grown from Mexico to Bolivia. Its ancestor is unknown, but it has limited compatibility with wild xerophytic Cucurbita from Mexico's highlands. We assembled the reference genome of C. ficifolia and assessed the genetic diversity and historical demography of the crop in Mexico with 2524 nuclear single nucleotide polymorphisms (SNPs). We also evaluated the gene flow between C. ficifolia and xerophytic taxa with 6292 nuclear and 440 plastome SNPs from 142 individuals sampled in 58 sites across their area of sympatry. Demographic modelling of C. ficifolia supports an eight-fold decrease in effective population size at about 2409 generations ago (95% CI = 464-12,393), whereas plastome SNPs support the expansion of maternal lineages ca. 1906-3635 years ago. Our results suggest a recent spread of C. ficifolia in Mexico, with high genetic diversity (π = 0.225, FST = 0.074) and inbreeding (FIS = 0.233). Coalescent models suggest low rates of gene flow with C. radicans and C. pedatifolia, whereas ABBA-BABA tests did not detect significant gene flow with wild taxa. Despite the ecogeographic proximity of C. ficifolia and its relatives, this crop persists as a highly isolated lineage of puzzling origin.

Influence of plant domestication on plant-pollinator interactions: Floral attributes and floral visitor communities in wild and cultivated squash plants.

PREMISE: Domestication of plant species results in phenotypic modifications and changes in biotic interactions. Most studies have compared antagonistic plant-herbivore interactions of domesticated plants and their wild relatives, but little attention has been given to how domestication influences plant-pollinator interactions. Floral attributes and interactions of floral visitors were compared between sister taxa of the genus Cucurbita (Cucurbitaceae), the domesticated C. moschata, C. argyrosperma ssp. argyrosperma and its wild progenitor C. argyrosperma ssp. sororia in the place of origin. METHODS: We conducted univariate and multivariate analyses to compare floral morphological traits and analyzed floral reward (nectar and pollen) quantity and quality between flowers of wild and domesticated Cucurbita taxa. Staminate and pistillate flowers of all three taxa were video recorded, and visitation and behavior of floral visitors were registered and analyzed. RESULTS: Most floral morphological characteristics of flowers of domesticated taxa were larger in both staminate and pistillate flowers. Staminate and pistillate flowers presented distinct correlations between floral traits and integration indices between domesticated and wild species. Additionally, pollen quantity and protein to lipid ratio were greater in domesticated species. Cucurbit pollen specialists, Eucera spp., had the highest probability of visit for all Cucurbita taxa. CONCLUSIONS: We provide evidence that floral traits of domesticated and wild Cucurbita species experienced different selection pressures. Domesticated Cucurbita species may have more resources invested towards floral traits, thereby increasing attractiveness to pollinators and potentially plant reproductive success. Wild ancestor plant populations should be conserved in their centers of origin to preserve plant-pollinator interactions.

Incorporating evolutionary and threat processes into crop wild relatives conservation.

Crop wild relatives (CWR) intra- and interspecific diversity is essential for crop breeding and food security. However, intraspecific genetic diversity, which is central given the idiosyncratic threats to species in landscapes, is usually not considered in planning frameworks. Here, we introduce an approach to develop proxies of genetic differentiation to identify conservation areas, applying systematic conservation planning tools that produce hierarchical prioritizations of the landscape. It accounts for: (i) evolutionary processes, including historical and environmental drivers of genetic diversity, and (ii) threat processes, considering taxa-specific tolerance to human-modified habitats, and their extinction risk status. Our analyses can be used as inputs for developing national action plans for the conservation and use of CWR. Our results also inform public policy to mitigate threat processes to CWR (like crops living modified organisms or agriculture subsidies), and could advise future research (e.g. for potential germplasm collecting). Although we focus on Mesoamerican CWR within Mexico, our methodology offers opportunities to effectively guide conservation and monitoring strategies to safeguard the evolutionary resilience of any taxa, including in regions of complex evolutionary histories and mosaic landscapes.

Potential Distribution of Cedrela odorata L. in Mexico according to Its Optimal Thermal Range for Seed Germination under Different Climate Change Scenarios.

Cedrela odorata is a native tree of economic importance, as its wood is highly demanded in the international market. In this work, the current and future distributions of C. odorata in Mexico under climate change scenarios were analyzed according to their optimal temperature ranges for seed germination. For the present distribution, 256 localities of the species' presence were obtained from the Global Biodiversity Information Facility (GBIF) database and modelled with MaxEnt. For the potential distribution, the National Center for Atmospheric Research model (CCSM4) was used under conservative and drastic scenarios (RCP2.6 and RCP8.5 Watts/m2, respectively) for the intermediate future (2050) and far future (2070). Potential distribution models were built from occurrence data within the optimum germination temperature range of the species. The potential distribution expanded by 5 and 7.8% in the intermediate and far future, respectively, compared with the current distribution. With the increase in temperature, adequate environmental conditions for the species distribution should be met in the central Mexican state of Guanajuato. The states of Chihuahua, Mexico, Morelos, Guerrero, and Durango presented a negative trend in potential distribution. Additionally, in the far future, the state of Chihuahua it is likely to not have adequate conditions for the presence of the species. For the prediction of the models, the precipitation variable during the driest month presented the greatest contribution. When the humidity is not limiting, the thermal climatic variables are the most important ones. Models based on its thermal niche for seed germination allowed for the identification of areas where temperature will positively affect seed germination, which will help maximize the establishment of plant populations and adaptation to different climate change scenarios.

Thermal Niche for Seed Germination and Species Distribution Modelling of Swietenia macrophylla King (Mahogany) under Climate Change Scenarios.

Swietenia macrophylla is an economically important tree species propagated by seeds that lose their viability in a short time, making seed germination a key stage for the species recruitment. The objective of this study was to determine the cardinal temperatures and thermal time for seed germination of S. macrophylla; and its potential distribution under different climate change scenarios. Seeds were placed in germination chambers at constant temperatures from 5 to 45 °C and their thermal responses modelled using a thermal time approach. In addition, the potential biogeographic distribution was projected according to the Community Climate System Model version 4 (CCSM4). Germination rate reached its maximum at 37.3 ± 1.3 °C (To); seed germination decreased to near zero at 52.7 ± 2.2 °C (ceiling temperature, Tc) and at 12.8 ± 2.4 °C (base temperature, Tb). The suboptimal thermal time θ150 needed for 50% germination was ca. 190 °Cd, which in the current scenario is accumulated in 20 days. The CCSM4 model estimates an increase of the potential distribution of the species of 12.3 to 18.3% compared to the current scenario. The temperature had an important effect on the physiological processes of the seeds. With the increase in temperature, the thermal needs for germination are completed in less time, so the species will not be affected in its distribution. Although the distribution of the species may not be affected, it is crucial to generate sustainable management strategies to ensure its long-term conservation.

The domestication of Cucurbita argyrosperma as revealed by the genome of its wild relative.

Despite their economic importance and well-characterized domestication syndrome, the genomic impact of domestication and the identification of variants underlying the domestication traits in Cucurbita species (pumpkins and squashes) is currently lacking. Cucurbita argyrosperma, also known as cushaw pumpkin or silver-seed gourd, is a Mexican crop consumed primarily for its seeds rather than fruit flesh. This makes it a good model to study Cucurbita domestication, as seeds were an essential component of early Mesoamerican diet and likely the first targets of human-guided selection in pumpkins and squashes. We obtained population-level data using tunable Genotype by Sequencing libraries for 192 individuals of the wild and domesticated subspecies of C. argyrosperma across Mexico. We also assembled the first high-quality wild Cucurbita genome. Comparative genomic analyses revealed several structural variants and presence/absence of genes related to domestication. Our results indicate a monophyletic origin of this domesticated crop in the lowlands of Jalisco. We found evidence of gene flow between the domesticated and wild subspecies, which likely alleviated the effects of the domestication bottleneck. We uncovered candidate domestication genes that are involved in the regulation of growth hormones, plant defense mechanisms, seed development, and germination. The presence of shared selected alleles with the closely related species Cucurbita moschata suggests domestication-related introgression between both taxa.

Phylogeography and Genetic Diversity in a Southern North American Desert: Agave kerchovei From the Tehuacán-Cuicatlán Valley, Mexico.

The Tehuacán-Cuicatlán Valley, located at the southeast of the state of Puebla and the northeast of the state of Oaxaca in Central Mexico, south of the Trans-Mexican Volcanic Belt (TMVB), is of particular interest for understanding the evolutionary dynamics of arid and semi-arid environments, being one of the main reservoirs of biological diversity for the arid zones of North America, including the highest diversity of Agavaceae worldwide and high levels of endemism. Studying in detail the phylogeography, environmental history and population genetics of representative species will hopefully shed light on the evolutionary and ecological dynamics that generated the tremendous biodiversity and endemism of this important region in Mexico. We sequenced three non-coding regions of chloroplast genome of Agave kerchovei, a representative species of the Tehuacán Valley, generating 2,188 bp from 128 individuals sampled from eight populations throughout the species range. We used this data set to (i) characterize the levels of genetic diversity and genetic structure in A. kerchovei; (ii) predict the distribution of A. kerchovei for the present day, and to reconstruct the past geographical history of the species by constructing ecological niche models (ENM); and (iii) compare the levels of diversity in this species with those estimated for the widely distributed Agave lechuguilla. Agave kerchovei has high levels of total chloroplast genetic variation (Hd = 0.718), especially considering that it is a species with a very restricted distribution. However, intrapopulation diversity is low (zero in some populations), and genetic structure is high (F ST = 0.928, G ST = 0.824), which can be expected for endemic species with isolated populations. Our data suggest that Pleistocene glacial cycles have played an important role in the distribution of A. kerchovei, where the climatic variability of the region - likely associated with its topographic complexity - had a significant effect on the levels of genetic diversity and population dynamics, while the potential distribution of the species seems to be stable since the middle Holocene (6 kya). We conclude that in A. kerchovei there is a core group of populations in the Tehuacán Valley, and peripheric populations that appear to be evolving independently and thus the species is fundamentally an endemic species from the Tehuacán Valley while the populations outside the Valley appear to be in the process of incipient speciation.

Medicinal plants of Oaxaca, Mexico: ethnobotany and antibacterial activity / Plantas medicinales de Oaxaca, México: etnobotánica y actividad antibacteriana

Santiago Quiotepec, one of the oldest communities of the Tehuacán-Cuicatlán Valley (México), has a great tradition using medicinal plants. The aim of this study was to make an inventory of the medicinal species used by the inhabitants of Santiago Quiotepec and evaluate the antibacterial activity. An ethnobotanical study of medicinal plants was carried out, 60 informants mentioned that 66 species of plants are being used in the treatment of different diseases. Fifteen species were selected to evaluate the antibacterial activity in possible bacterial originated diseases treatment. The lowest values were presented in the hexane extract of Plumbago pulchella, with a MIC of 0.25 mg/mL over Staphylococcus aureus and S. epidermidis as well as the hexanic extract of Echinopterys eglandulosa showed a MIC of 0.25 mg/mL over Pseudomona aeruginosa.

Santiago Quiotepec es una de las comunidades más antiguas del valle de Tehuacán-Cuicatlán (México), y tiene una gran tradición en el uso de plantas medicinales. El objetivo de este estudio fue realizar un inventario de las especies medicinales utilizadas por los habitantes de Santiago Quiotepec y evaluar la actividad antibacteriana. Se realizó un estudio etnobotánico de plantas medicinales, 60 informantes mencionaron 66 especies de plantas utilizadas en el tratamiento de diferentes enfermedades. Quince especies utilizadas en la comunidad para tratar enfermedades de posible origen bacteriano fueron seleccionadas para evaluar la actividad antibacteriana. Los valores más bajos se presentaron en el extracto hexánico de Plumbago pulchella, con una CMI de 0.25 mg/ml sobre Staphylococcus aureus y S. epidermidis, así como el extracto hexánico de Echinopterys eglandulosa mostró una CMI de 0.25 mg/mL sobre Pseudomona aeruginosa.

Phylogeographic and population genetic analyses of Cucurbita moschata reveal divergence of two mitochondrial lineages linked to an elevational gradient.

PREMISE: Domestication usually involves local adaptation to environmental conditions. Cucurbita species are a promising model for studying these processes. Cucurbita moschata is the third major crop in the genus because of its economic value and because it displays high landrace diversity, but research about its genetic diversity, population structure, and phylogeography is limited. We aimed at understanding how geography and elevation shape the distribution of genetic diversity in C. moschata landraces in Mexico. METHODS: We sampled fruits from 24 localities throughout Mexico. We assessed 11 nuclear microsatellite loci, one mtDNA region, and three cpDNA regions but found no variation in cpDNA. We explored genetic structure with cluster analysis, and phylogeographic relationships with haplotype network analysis. RESULTS: Mitochondrial genetic diversity was high, and nuclear genetic differentiation among localities was intermediate compared to other domesticated Cucurbita. We found high levels of inbreeding. We recovered two mitochondrial lineages: highland (associated with the Trans-Mexican Volcanic Belt) and lowland. Nuclear microsatellites show that localities from the Yucatan Peninsula constitute a well-differentiated group. CONCLUSIONS: Mexico is an area of high diversity for C. moschata, and these landraces represent important plant genetic resources. In Mexico this species is characterized by divergence processes linked to an elevational gradient, which could be related to adaptation and may be of value for applications in agriculture. The Isthmus of Tehuantepec may be a partial barrier to gene flow. Morphological variation, agricultural management, and cultural differences may be related to this pattern of genetic structure, but further studies are needed.

Evolutionary Dynamics of Transferred Sequences Between Organellar Genomes in Cucurbita.

Twenty-nine DNA regions of plastid origin have been previously identified in the mitochondrial genome of Cucurbita pepo (pumpkin; Cucurbitaceae). Four of these regions harbor homolog sequences of rbcL, matK, rpl20-rps12 and trnL-trnF, which are widely used as molecular markers for phylogenetic and phylogeographic studies. We extracted the mitochondrial copies of these regions based on the mitochondrial genome of C. pepo and, along with published sequences for these plastome markers from 13 Cucurbita taxa, we performed phylogenetic molecular analyses to identify inter-organellar transfer events in the Cucurbita phylogeny and changes in their nucleotide substitution rates. Phylogenetic reconstruction and tree selection tests suggest that rpl20 and rbcL mitochondrial paralogs arose before Cucurbita diversification whereas the mitochondrial matK and trnL-trnF paralogs emerged most probably later, in the mesophytic Cucurbita clade. Nucleotide substitution rates increased one order of magnitude in all the mitochondrial paralogs compared to their original plastid sequences. Additionally, mitochondrial trnL-trnF sequences obtained by PCR from nine Cucurbita taxa revealed higher nucleotide diversity in the mitochondrial than in the plastid copies, likely related to the higher nucleotide substitution rates in the mitochondrial region and loss of functional constraints in its tRNA genes.

Tracing back the origin of pumpkins (Cucurbita pepo ssp. pepo L.) in Mexico.

Cucurbita pepo is an economically important crop, which consists of cultivated C. pepo ssp. pepo, and two wild taxa (C. pepo ssp. fraterna and C. pepo ssp. ovifera). We aimed at understanding the domestication and the diversity of C. pepo in Mexico. We used two chloroplast regions and nine nuclear microsatellite loci to assess the levels of genetic variation and structure for C. pepo ssp. pepo's landraces sampled in 13 locations in Mexico, five improved varieties, one C. pepo ssp. fraterna population and ornamental C. pepo ssp. ovifera. We tested four hypotheses regarding the origin of C. pepo ssp. pepo's ancestor through approximate Bayesian computation: C. pepo ssp. ovifera as the ancestor; C. pepo ssp. fraterna as the ancestor; an unknown extinct lineage as the ancestor; and C. pepo ssp. pepo as hybrid from C. pepo ssp. ovifera and C. pepo ssp. fraterna ancestors. Cucurbita pepo ssp. pepo showed high genetic variation and low genetic differentiation. Cucurbita pepo ssp. fraterna and C. pepo ssp. pepo shared two chloroplast haplotypes. The three subspecies were well differentiated for microsatellite loci. Cucurbita pepo ssp. fraterna was probably C. pepo ssp. pepo's wild ancestor, but subsequent hybridization between taxa complicate defining C. pepo ssp. pepo's ancestor.

The Genome of Cucurbita argyrosperma (Silver-Seed Gourd) Reveals Faster Rates of Protein-Coding Gene and Long Noncoding RNA Turnover and Neofunctionalization within Cucurbita.

Whole-genome duplications are an important source of evolutionary novelties that change the mode and tempo at which genetic elements evolve within a genome. The Cucurbita genus experienced a whole-genome duplication around 30 million years ago, although the evolutionary dynamics of the coding and noncoding genes in this genus have not yet been scrutinized. Here, we analyzed the genomes of four Cucurbita species, including a newly assembled genome of Cucurbita argyrosperma, and compared the gene contents of these species with those of five other members of the Cucurbitaceae family to assess the evolutionary dynamics of protein-coding and long intergenic noncoding RNA (lincRNA) genes after the genome duplication. We report that Cucurbita genomes have a higher protein-coding gene birth-death rate compared with the genomes of the other members of the Cucurbitaceae family. C. argyrosperma gene families associated with pollination and transmembrane transport had significantly faster evolutionary rates. lincRNA families showed high levels of gene turnover throughout the phylogeny, and 67.7% of the lincRNA families in Cucurbita showed evidence of birth from the neofunctionalization of previously existing protein-coding genes. Collectively, our results suggest that the whole-genome duplication in Cucurbita resulted in faster rates of gene family evolution through the neofunctionalization of duplicated genes.

Historical biogeography and phylogeny of Cucurbita: Insights from ancestral area reconstruction and niche evolution.

Knowledge of the role of geographical and ecological events associated to the divergence process of wild progenitors is important to understand the process of domestication. We analysed the temporal, spatial and ecological patterns of the diversification of Cucurbita, an American genus of worldwide economic importance. We conducted a phylogenetic analysis based on six chloroplast regions (5907 bp) to estimate diversification rates and dates of divergence between taxa. This is the first phylogenetic study to include C. radicans, a wild species that is endemic to the Trans Mexican Volcanic Belt. We performed analysis of ancestral area reconstruction and paleoreconstructions of species distribution models to understand shifts in wild species ranges. We used principal component analysis (PCA) and multivariate analysis of variance (MANOVA) to evaluate the environmental differentiation among taxa within each clade. The phylogenetic analyses showed good support for at least six independent domestication events in Cucurbita. The genus Cucurbita showed a time of divergence of 11.24 Ma (6.88-17 Ma 95% HDP), and the dates of divergence between taxa within each group ranged from 0.35 to 6.58 Ma, being the divergence between C. lundelliana and C. okeechobeensis subsp. martinezii the most recent. The diversification rate of the genus was constant through time. The diversification of most wild taxa occurred during the Pleistocene, and its date of divergence is concordant with the dates of divergence reported for specialized bees of the genera Xenoglossa and Peponapis, suggesting a process of coevolution between Cucurbita and their main pollinators that should be further investigated. Tests of environmental differentiation together with ancestral area reconstruction and species distribution models past projections suggest that divergence was promoted by the onset of geographic barriers and secondary range contraction and by expansion related to glacial-interglacial cycles.

Genetic Resources in the "Calabaza Pipiana" Squash (Cucurbita argyrosperma) in Mexico: Genetic Diversity, Genetic Differentiation and Distribution Models.

Analyses of genetic variation allow understanding the origin, diversification and genetic resources of cultivated plants. Domesticated taxa and their wild relatives are ideal systems for studying genetic processes of plant domestication and their joint is important to evaluate the distribution of their genetic resources. Such is the case of the domesticated subspecies C. argyrosperma ssp. argyrosperma, known in Mexico as calabaza pipiana, and its wild relative C. argyrosperma ssp. sororia. The main aim of this study was to use molecular data (microsatellites) to assess the levels of genetic variation and genetic differentiation within and among populations of domesticated argyrosperma across its distribution in Mexico in comparison to its wild relative, sororia, and to identify environmental suitability in previously proposed centers of domestication. We analyzed nine unlinked nuclear microsatellite loci to assess levels of diversity and distribution of genetic variation within and among populations in 440 individuals from 19 populations of cultivated landraces of argyrosperma and from six wild populations of sororia, in order to conduct a first systematic analysis of their genetic resources. We also used species distribution models (SDMs) for sororia to identify changes in this wild subspecies' distribution from the Holocene (∼6,000 years ago) to the present, and to assess the presence of suitable environmental conditions in previously proposed domestication sites. Genetic variation was similar among subspecies (HE = 0.428 in sororia, and HE = 0.410 in argyrosperma). Nine argyrosperma populations showed significant levels of inbreeding. Both subspecies are well differentiated, and genetic differentiation (FST) among populations within each subspecies ranged from 0.152 to 0.652. Within argyrosperma we found three genetic groups (Northern Mexico, Yucatan Peninsula, including Michoacan and Veracruz, and Pacific coast plus Durango). We detected low levels of gene flow among populations at a regional scale (<0.01), except for the Yucatan Peninsula, and the northern portion of the Pacific Coast. Our analyses suggested that the Isthmus of Tehuantepec is an effective barrier isolating southern populations. Our SDM results indicate that environmental characteristics in the Balsas-Jalisco region, a potential center of domestication, were suitable for the presence of sororia during the Holocene.