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.

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.

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.

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.