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.

Differences in the genomic diversity, structure, and inbreeding patterns in wild and managed populations of Agave potatorum Zucc. used in the production of Tobalá mezcal in Southern Mexico.

Agave potatorum Zucc. locally known as Tobalá, is an important species for mezcal production. It is a perennial species that takes 10 to 15 years to reach reproductive age. Because of high demand of Tobalá mezcal and the slow maturation of the plants, its wild populations have been under intense anthropogenic pressure. The main objective of this study was to estimate the genome-wide diversity in A. potatorum and determine if the type of management has had any effect on its diversity, inbreeding and structure. We analyzed 174 individuals (105 wild, 42 cultivated and 27 from nurseries) from 34 sites with a reduced representation genomic method (ddRADseq), using 14,875 SNPs. The diversity measured as expected heterozygosity was higher in the nursery and wild plants than in cultivated samples. We did not find private alleles in the cultivated and nursery plants, which indicates that the individuals under management recently derived from wild populations, which was supported by higher gene flow estimated from wild populations to the managed plants. We found low but positive levels of inbreeding (FIS = 0.082), probably related to isolation of the populations. We detected low genetic differentiation among populations (FST = 0.0796), with positive and significant isolation by distance. The population genetic structure in the species seems to be related to elevation and ecology, with higher gene flow among populations in less fragmented areas. We detected an outlier locus related to the recognition of pollen, which is also relevant to self-incompatibility protein (SI). Due to seed harvest and long generation time, the loss of diversity in A. potatorum has been gradual and artificial selection and incipient management have not yet caused drastic differences between cultivated and wild plants. Also, we described an agroecological alternative to the uncontrolled extraction of wild individuals.

Historical, temporal, and geographic dynamism of the interaction between Agave and Leptonycteris nectar-feeding bats.

PREMISE: The interaction between ecological and evolutionary processes has been recognized as an important factor shaping the evolutionary history of species. Some authors have proposed different ecological and evolutionary hypotheses concerning the relationships between plants and their pollinators; a special case is the interaction and suspected coevolution among Agave spp. and their main pollinators, the Leptonycteris bats. Agave spp. have, in general, a pollination syndrome compatible with chiropterophily including floral shape and size, nocturnal nectar production, and nectar quality and sugar concentration. Our goal was to analyze the interaction Agave-Leptonycteris and its dynamics during three different climate scenarios. METHODS: We modeled the Agave-Leptonycteris interaction in its spatial and temporal components during the Pleistocene using Ecological Niche Models (ENMs) and three climate scenarios: Current, Last Glacial Maximum (LGM), and Last InterGlacial (LIG). Furthermore, we analyzed the geographic correlation between 96 Agave spp. and two of the Mexican Tequila bats, genus Leptonycteris. RESULTS: We found that Leptonycteris spp. interact with different Agave spp. over their migratory routes. We propose an interaction refuge in Metztitlán and Tehuacán-Cuicatlán areas, where Agave- Leptonycteris interaction has probably remained active. During the nonmigratory season, both bat species consume nectar of almost the same Agave spp., suggesting the possibility of a diffuse coevolution among Agave and Leptonycteris bats. CONCLUSIONS: We propose that in the areas related to migratory bat movements, each bat species interacts with different Agave spp., whereas in the areas occupied by nonmigrant individuals, both bat species consume nectar of almost the same Agave taxa.

Genomic analysis unveils reduced genetic variability but increased proportion of heterozygotic genotypes of the intensively managed mezcal agave, Agave angustifolia.

PREMISE: The central Oaxaca Basin has a century-long history of agave cultivation and is hypothesized to be the region of origin of other cultivated crops. Widely cultivated for mezcal production, the perennial crop known as "espadín" is putatively derived from wild Agave angustifolia. Nevertheless, little is known about its genetic relationship to the wild A. angustifolia or how the decades-long clonal propagation has affected its genetics. METHODS: Using restriction-site-associated DNA sequencing and over 8000 single-nucleotide polymorphisms, we studied aspects of the population genomics of wild and cultivated A. angustifolia in Puebla and Oaxaca, Mexico. We assessed patterns of genetic diversity, inbreeding, distribution of genetic variation, and differentiation among and within wild populations and plantations. RESULTS: Genetic differentiation between wild and cultivated plants was strong, and both gene pools harbored multiple unique alleles. Nevertheless, we found several cultivated individuals with high genetic affinity with wild samples. Higher heterozygosity was observed in the cultivated individuals, while in total, they harbored considerably fewer alleles and presented higher linkage disequilibrium compared to the wild plants. Independently of geographic distance among sampled plantations, the genetic relatedness of the cultivated plants was high, suggesting a common origin and prevalent role of clonal propagation. CONCLUSIONS: The considerable heterozygosity found in espadín is contained within a network of highly related individuals, displaying high linkage disequilibrium generated by decades of clonal propagation and possibly by the accumulation of somatic mutations. Wild A. angustifolia, on the other hand, represents a significant genetic diversity reservoir that should be carefully studied and conserved.

Conservation genomics of Agave tequilana Weber var. azul: low genetic differentiation and heterozygote excess in the tequila agave from Jalisco, Mexico.

Background: Genetic diversity is fundamental for the survival of species. In particular, in a climate change scenario, it is crucial that populations maintain genetic diversity so they can adapt to novel environmental conditions. Genetic diversity in wild agaves is usually high, with low genetic differentiation among populations, in part maintained by the agave pollinators such as the nectarivorous bats. In cultivated agaves, patterns of genetic diversity vary according to the intensity of use, management, and domestication stage. In Agave tequilana Weber var. azul (A. tequilana thereafter), the plant used for tequila production, clonal propagation has been strongly encouraged. These practices may lead to a reduction in genetic diversity. Methods: We studied the diversity patterns with genome-wide SNPs, using restriction site associated DNA sequencing in cultivated samples of A. tequilana from three sites of Jalisco, Mexico. For one locality, seeds were collected and germinated in a greenhouse. We compared the genomic diversity, levels of inbreeding, genetic differentiation, and connectivity among studied sites and between adults and juvenile plants. Results: Agave tequilana presented a genomic diversity of HT = 0.12. The observed heterozygosity was higher than the expected heterozygosity. Adults were more heterozygous than juveniles. This could be a consequence of heterosis or hybrid vigor. We found a shallow genetic structure (average paired FST = 0.0044). In the analysis of recent gene flow, we estimated an average migration rate among the different populations of m = 0.25. In particular, we found a population that was the primary source of gene flow and had greater genomic diversity (HE and HO ), so we propose that this population should continue to be monitored as a potential genetic reservoir. Discussion: Our results may be the consequence of more traditional management in the studied specific region of Jalisco. Also, the exchange of seeds or propagules by producers and the existence of gene flow due to occasional sexual reproduction may play an important role in maintaining diversity in A. tequilana. For populations to resist pests, to continue evolving and reduce their risk of extinction under a climate change scenario, it is necessary to maintain genetic diversity. Under this premise we encourage to continue acting in conservation programs for this species and its pollinators.

Demographic modelling helps track the rapid and recent divergence of a conifer species pair from Central Mexico.

Secondary contact of recently diverged species may have several outcomes, ranging from rampant hybridization to reinforced reproductive isolation. In plants, selfing tolerance and disjunct reproductive phenology may lead to reproductive isolation at contact zones. However, they may also evolve under both allopatric or parapatric frameworks and originate from adaptive and/or neutral forces. Inferring the historical demography of diverging taxa is thus a crucial step to identify factors that may have led to putative reproductive isolation. We explored various competing demographypotheses to account for the rapid divergence of a fir species complex (Abies flinckii-A. religiosa) distributed in "sky-islands" across central Mexico (i.e., along the Trans-Mexican Volcanic Belt; TMVB). Despite co-occurring in two independent sympatric regions (west and centre), these taxa rarely interbreed because of disjunct reproductive phenologies. We genotyped 1147 single nucleotide polymorphisms, generated by GBS (genotyping by sequencing), across 23 populations, and compared multiple scenarios based on the geological history of the TMVB. The best-fitting model revealed one of the most rapid and complete speciation cases for a conifer species-pair, dating back to ~1.2 million years ago. Coupled with the lack of support for stepwise colonization, our coalescent inferences point to an early cessation of interspecific gene flow under parapatric speciation; ancestral gene flow during divergence was asymmetrical (mostly from western firs into A. religiosa) and exclusive to the most ancient (i.e., central) contact zone. Factors promoting rapid reproductive isolation should be explored in other slowly evolving species complexes as they may account for the large tropical and subtropical diversity.

Genomic Analyses of Wild and Cultivated Bacanora Agave (Agave angustifolia var. pacifica) Reveal Inbreeding, Few Signs of Cultivation History and Shallow Population Structure.

Due to the recent increase in demand for agave-based beverages, many wild agave populations have experienced rapid decline and fragmentation, whereas cultivated plants are now managed at monocultural plantations, in some cases involving clonal propagation. We examined the relative effect of migration, genetic drift, natural selection and human activities on the genetic repertoire of Agave angustifolia var. pacifica, an agave used for bacanora (an alcoholic spirit similar to tequila) production in northwestern Mexico. We sampled 34 wild and cultivated sites and used over eleven thousand genome-wide SNPs. We found shallow genetic structure among wild samples, although we detected differentiation between coastal and inland sites. Surprisingly, no differentiation was found between cultivated and wild populations. Moreover, we detected moderate inbreeding (FIS ~ 0.13) and similar levels of genomic diversity in wild and cultivated agaves. Nevertheless, the cultivated plants had almost no private alleles and presented evidence of clonality. The overall low genetic structure in A. angustifolia var. pacifica is apparently the result of high dispersibility promoted by pollinators and the possibility of clonal reproduction. Incipient cultivation history and reliance on wild seeds and plants are probably responsible for the observed patterns of high genetic connectivity and considerable diversity in cultivated samples.

Uses, Knowledge and Extinction Risk Faced by Agave Species in Mexico.

We compiled an updated database of all Agave species found in Mexico and analyzed it with specific criteria according to their biological parameters to evaluate the conservation and knowledge status of each species. Analyzing the present status of all Agave species not only provides crucial information for each species, but also helps determine which ones require special protection, especially those which are heavily used or cultivated for the production of distilled beverages. We conducted an extensive literature review search and compiled the conservation status of each species using mainstream criteria by IUCN. The information gaps in the database indicate a lack of knowledge and research regarding specific Agave species and it validates the need to conduct more studies on this genus. In total, 168 Agave species were included in our study, from which 89 are in the subgenus Agave and 79 in the subgenus Littaea. Agave lurida and A. nizandensis, in the subgenus Agave and Littaea, respectively, are severely endangered, due to their endemism, lack of knowledge about pollinators and floral visitors, and their endangered status according to the IUCN Red List. Some species are at risk due to the loss of genetic diversity resulting from production practices (i.e., Agave tequilana), and others because of excessive and unchecked overharvesting of wild plants, such as A. guadalajarana, A. victoriae-reginae, A. kristenii, and others. Given the huge economic and ecological importance of plants in the genus Agave, our review will be a milestone to ensure their future and continued provision of ecosystem services for humans, as well as encouraging further research in Agave species in an effort to enhance awareness of their conservation needs and sustainable use, and the implementation of eco-friendly practices in the species management.

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.

Evolutionary ecology of Agave: distribution patterns, phylogeny, and coevolution (an homage to Howard S. Gentry).

With more than 200 species, the genus Agave is one of the most interesting and complex groups of plants in the world, considering for instance its great diversity and adaptations. The adaptations include the production of a single, massive inflorescence (the largest among plants) where after growing for many years, sometimes more than 30, the rosette dies shortly afterward, and the remarkable coevolution with their main pollinators, nectarivorous bats, in particular of the genus Leptonycteris. The physiological adaptations of Agave species include a photosynthetic metabolism that allows efficient use of water and a large degree of succulence, helping to store water and resources for their massive flowering event. Ecologically, the agaves are keystone species on which numerous animal species depend for their subsistence due to the large amounts of pollen and nectar they produce, that support many pollinators, including bats, perching birds, hummingbirds, moths, and bees. Moreover, in many regions of Mexico and in the southwestern United States, agaves are dominant species. We describe the contributions of H. S. Gentry to the understanding of agaves and review recent advances on the study of the ecology and evolution of the genus. We analyze the present and inferred past distribution patterns of different species in the genus, describing differences in their climatic niche and adaptations to dry conditions. We interpret these patterns using molecular clock data and phylogenetic analyses and information of their coevolving pollinators and from phylogeographic, morphological, and ecological studies and discuss the prospects for their future conservation and management.

Non-adaptive evolutionary processes governed the diversification of a temperate conifer lineage after its migration into the tropics.

Constructing phylogenetic relationships among closely related species is a recurrent challenge in evolutionary biology, particularly for long-lived taxa with large effective population sizes and uncomplete reproductive isolation, like conifers. Conifers further have slow evolutionary rates, which raises the question of whether adaptive or non/adaptive processes were predominantly involved when they rapidly diversified after migrating from temperate regions into the tropical mountains. Indeed, fine-scale phylogenetic relationships within several conifer genus remain under debate. Here, we studied the phylogenetic relationships of endemic firs (Abies, Pinaceae) discontinuously distributed in the montane forests from the Southwestern United States to Guatemala, and addressed several hypotheses related to adaptive and non-adaptive radiations. We derived over 80 K SNPs from genotyping by sequencing (GBS) for 45 individuals of nine Mesoamerican species to perform phylogenetic analyses. Both Maximum Likelihood and quartets-inference phylogenies resulted in a well-resolved topology, showing a single fir lineage divided in four subgroups that coincided with the main mountain ranges of Mesoamerica; thus having important taxonomic implications. Such subdivision fitted a North-South isolation by distance framework, in which non-adaptive allopatric processes seemed the rule. Interestingly, several reticulations were observed within subgroups, especially in the central-south region, which may explain past difficulties for generating infrageneric phylogenies. Further evidence for non-adaptive processes was obtained from analyses of 21 candidate-gene regions, which exhibited diminishing values of πa/πs and Ka/Ks with latitude, thus indicating reduced efficiency of purifying selection towards the Equator. Our study indicates that non-adaptive allopatric processes may be key generators of species diversity and endemism in the tropics.

The role of environment, local adaptation, and past climate fluctuation on the amount and distribution of genetic diversity in two subspecies of Mexican wild Zea mays.

PREMISE: Past climate fluctuations during the Holocene and Pleistocene shaped the distribution of several plant species in temperate areas over the world. Wild maize, commonly known as teosinte, is a good system to evaluate the effects of historical climate fluctuations on genetic diversity due to its wide distribution in Mexico with contrasting environmental conditions. We explored the influence of contemporary factors and historical environmental shifts on genetic diversity, including present and three historical periods using neutral markers. METHODS: We used 22 nuclear microsatellite loci to examine the genetic diversity of 14 populations of Zea mays subsp. parviglumis and 15 populations of Zea mays subsp. mexicana (527 individuals total). We implemented genetic structure analyses to evaluate genetic differentiation between and within subspecies. We applied coalescent-based demographic analysis and species distribution modeling to evaluate the effects of historical environmental shifts. RESULTS: We found 355 alleles in total for the two subspecies and variable levels of diversity in each (Z. mays subsp. parviglumis expected heterozygosity HE = 0.3646-0.7699; Z. mays subsp. mexicana HE = 0.5885-0.7671). We detected significant genetic structure among populations (DEST = 0.4332) with significant heterozygote deficiency (FIS = 0.1796), and variable selfing rates (sg2 = 0.0-0.3090). The Bayesian assignment analysis differentiated four genetic groups. Demographic and species distribution modeling analysis suggested that environmental shifts were influential in the amount of genetic diversity. CONCLUSIONS: Our analyses suggest that the current genetic diversity in teosinte is shaped by factors such as local adaptation and genetic isolation, along with historical environmental fluctuations.

Phylogeny, Diversification Rate, and Divergence Time of Agave sensu lato (Asparagaceae), a Group of Recent Origin in the Process of Diversification.

Agave sensu lato is one of the most diverse and complex genera of Asparagaceae, with more than 250 species. The morphological, ecological, and evolutionary diversity of the group has complicated its taxonomical study. We conducted phylogenetic analyses of DNA sequence data to reconstruct the phylogenetic relationships of the Agave genus. We included 107 species of the Asparagaceae family from which 83 correspond to the Agave sensu lato clade (Agave sensu stricto + Polianthes + Manfreda and Prochnyanthes, which together represent 30% of the genus) and as outgroups the genera Dasylirion, Hesperoyucca, Chlorogalum, Camassia, Hesperaloe, Yucca, Beschorneria, and Furcraea, in order to estimate the age and propose the history of their diversification. Previous studies postulated the relevance of the Miocene in the speciation rates of the agaves, as well as the relevance of the type of inflorescence in its diversification. However, these assertions have not been well supported. The analysis of chloroplast regions resulted in low resolution, which could be the consequence of the few variable sites. On the other hand, the internal transcribed spacer (ITS) implemented in our analysis ensued in higher resolution and better support values. Our phylogenetic analyses recovered five groups; one is the Striatae group, which is the sister group to Agave sensu stricto clade. Within this clade, we found three main groups with high support; these groups are not related with previous morphological proposals. We also analyzed the dates of origin and diversification rates. A Bayesian analysis of macroevolutionary mixtures indicated two significant shifts; the first was identified at 6.18 Ma, where the speciation rate increased to 4.10 species/Mya, this shift occurred during the late Miocene period, characterized by the emergence of arid biomes in North America. The second was identified at a stem age of 2.68 Ma where the speciation rate increased to 6.04 species/Mya. Concerning the ancestral reconstruction state of the inflorescence type in the Agave sensu stricto clade, the spike inflorescence character was predominant in the early-diverging groups, whereas the late-diverging groups present panicle inflorescences as the predominant character and higher speciation rates.

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.

Genomic consequences of dietary diversification and parallel evolution due to nectarivory in leaf-nosed bats.

BACKGROUND: The New World leaf-nosed bats (Phyllostomids) exhibit a diverse spectrum of feeding habits and innovations in their nutrient acquisition and foraging mechanisms. However, the genomic signatures associated with their distinct diets are unknown. RESULTS: We conducted a genomic comparative analysis to study the evolutionary dynamics related to dietary diversification and specialization. We sequenced, assembled, and annotated the genomes of five Phyllostomid species: one insect feeder (Macrotus waterhousii), one fruit feeder (Artibeus jamaicensis), and three nectar feeders from the Glossophaginae subfamily (Leptonycteris yerbabuenae, Leptonycteris nivalis, and Musonycteris harrisoni), also including the previously sequenced vampire Desmodus rotundus. Our phylogenomic analysis based on 22,388 gene families displayed differences in expansion and contraction events across the Phyllostomid lineages. Independently of diet, genes relevant for feeding strategies and food intake experienced multiple expansions and signatures of positive selection. We also found adaptation signatures associated with specialized diets: the vampire exhibited traits associated with a blood diet (i.e., coagulation mechanisms), whereas the nectarivore clade shares a group of positively selected genes involved in sugar, lipid, and iron metabolism. Interestingly, in fruit-nectar-feeding Phyllostomid and Pteropodids bats, we detected positive selection in two genes: AACS and ALKBH7, which are crucial in sugar and fat metabolism. Moreover, in these two proteins we found parallel amino acid substitutions in conserved positions exclusive to the tribe Glossophagini and to Pteropodids. CONCLUSIONS: Our findings illuminate the genomic and molecular shifts associated with the evolution of nectarivory and shed light on how nectar-feeding bats can avoid the adverse effects of diets with high glucose content.

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.

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.