Transcriptome analysis reveals molecular mechanisms underlying chloroplast biogenesis in albino Agave angustifolia plantlets.

Albino plants display partial or complete loss of photosynthetic pigments and defective thylakoid membrane development, consequently impairing plastid function and development. These distinctive attributes render albino plants excellent models for investigating chloroplast biogenesis. Despite their potential, limited exploration has been conducted regarding the molecular alterations underlying these phenotypes, extending beyond photosynthetic metabolism. In this study, we present a novel de novo transcriptome assembly of an albino somaclonal variant of Agave angustifolia Haw., which spontaneously emerged during the micropropagation of green plantlets. Additionally, RT-qPCR analysis was employed to validate the expression of genes associated with chloroplast biogenesis, and plastome copy numbers were quantified. This research aims to gain insight into the molecular disruptions affecting chloroplast development and ascertain whether the expression of critical genes involved in plastid development and differentiation is compromised in albino tissues of A. angustifolia. Our transcriptomic findings suggest that albino Agave plastids exhibit high proliferation, activation of the protein import machinery, altered transcription directed by PEP and NEP, dysregulation of plastome expression genes, reduced expression of photosynthesis-associated nuclear genes, disruption in the tetrapyrrole and carotenoid biosynthesis pathway, alterations in the plastid ribosome, and an increased number of plastome copies, among other alterations.

A newly identified glycosyltransferase AsRCOM provides resistance to purple curl leaf disease in agave.

BACKGROUND: Purple curl leaf disease brings a significant threat to the development of agave industry, the underlying mechanism of disease-resistant Agave sisalana. hybrid 11648 (A. H11648R) is still unknown. RESULTS: To excavate the crucial disease-resistant genes against purple curl leaf disease, we performed an RNA-seq analysis for A.H11648R and A.H11648 during different stages of purple curl leaf disease. The DEGs (differentially expressed genes) were mainly enriched in linolenic acid metabolism, starch and sucrose mechanism, phenylpropanoid biosynthesis, hypersensitive response (HR) and systemic acquired resistance. Further analysis suggested that eight candidate genes (4'OMT2, ACLY, NCS1, GTE10, SMO2, FLS2, SQE1 and RCOM) identified by WGCNA (weighted gene co-expression network analysis) may mediate the resistance to agave purple curl disease by participating the biosynthesis of benzylisoquinoline alkaloids, steroid, sterols and flavonoids, and the regulation of plant innate immunity and systemic acquired resistance. After qPCR verification, we found that AsRCOM, coding a glycosyltransferase and relevant to the regulation of plant innate immunity and systemic acquired resistance, may be the most critical disease-resistant gene. Finally, the overexpression of AsRCOM gene in agave could significantly enhance the resistance to purple curl disease with abundant reactive oxygen species (ROS) accumulations. CONCLUSIONS: Integrative RNA-seq analysis found that HR may be an important pathway affecting the resistance to purple curl leaf disease in agave, and identified glycosyltransferase AsRCOM as the crucial gene that could significantly enhance the resistance to purple curl leaf disease in agave, with obvious ROS accumulations.

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.

Phylogenomics and gene selection in Aspergillus welwitschiae: Possible implications in the pathogenicity in Agave sisalana.

Aspergillus welwitschiae causes bole rot disease in sisal (Agave sisalana and related species) which affects the production of natural fibers in Brazil, the main worldwide producer of sisal fibers. This fungus is a saprotroph with a broad host range. Previous research established A. welwitschiae as the only causative agent of bole rot in the field, but little is known about the evolution of this species and its strains. In this work, we performed a comparative genomics analysis of 40 Aspergillus strains. We show the conflicting molecular identity of this species, with one sisal-infecting strain sharing its last common ancestor with Aspergillus niger, having diverged only 833 thousand years ago. Furthermore, our analysis of positive selection reveals sites under selection in genes coding for siderophore transporters, Sodium­calcium exchangers, and Phosphatidylethanolamine-binding proteins (PEBPs). Herein, we discuss the possible impacts of these gene functions on the pathogenicity in sisal.

Revealing floral metabolite network in tuberose that underpins scent volatiles synthesis, storage and emission.

KEY MESSAGE: The role of central carbon metabolism in the synthesis and emission of scent volatiles in tuberose flowers was revealed through measurement of changes in transcripts and metabolites levels. Tuberose or Agave amica (Medikus) Thiede & Govaerts is a widely cultivated ornamental plant in several subtropical countries. Little is known about metabolite networking involved in biosynthesis of specialized metabolites utilizing primary metabolites. In this study, metabolite profiling and gene expression analyses were carried out from six stages of maturation throughout floral lifespan. Multivariate analysis indicated distinction between early and late maturation stages. Further, the roles of sugars viz. sucrose, glucose and fructose in synthesis, glycosylation and emission of floral scent volatiles were studied. Transcript levels of an ABC G family transporter (picked up from the floral transcriptome) was in synchronization with terpene volatiles emission during the anthesis stage. A diversion from phenylpropanoid/benzenoid to flavonoid metabolism was observed as flowers mature. Further, it was suggested that this metabolic shift could be mediated by isoforms of 4-Coumarate-CoA ligase along with Myb308 transcription factor. Maximum glycosylation of floral scent volatiles was shown to occur at the late mature stage when emission declined, facilitating both storage and export from the floral tissues. Thus, this study provides an insight into floral scent volatiles synthesis, storage and emission by measuring changes at transcripts and metabolites levels in tuberose throughout floral lifespan.

Transcriptional profiling of the CAM plant Agave salmiana reveals conservation of a genetic program for regeneration.

In plants, the best characterized plant regeneration process is de novo organogenesis. This type of regeneration is characterized by the formation of a multicellular structure called callus. Calli are induced via phytohormone treatment of plant sections. The callus formation in plants like Agave species with Crassulacean Acid Metabolism (CAM) is poorly studied. In this study, we induced callus formation from Agave salmiana leaves and describe cell arrangement in this tissue. Moreover, we determined and analyzed the transcriptional program of calli, as well as those of differentiated root and leaf tissues, by using RNA-seq. We were able to reconstruct 170,844 transcripts of which 40,644 have a full Open Reading Frame (ORF). The global profile obtained by Next Generation Sequencing (NGS) reveals that several callus-enriched protein coding transcripts are orthologs of previously reported factors highly expressed in Arabidopsis calli. At least 62 genes were differentially expressed in Agave calli, 50 of which were up-regulated. Several of these are actively involved in the perception of, and response to, auxin and cytokinin. Not only are these the first results for the A. salmiana callus, but they provide novel data from roots and leaves of this Agave species, one of the largest non-tree plants in nature.

Transcriptome analysis of bolting in A. tequilana reveals roles for florigen, MADS, fructans and gibberellins.

BACKGROUND: Reliable indicators for the onset of flowering are not available for most perennial monocarpic species, representing a drawback for crops such as bamboo, agave and banana. The ability to predict and control the transition to the reproductive stage in A. tequilana would represent an advantage for field management of agaves for tequila production and for the development of a laboratory model for agave species. RESULTS: Consistent morphological features could not be determined for the vegetative to reproductive transition in A. tequilana. However, changes in carbohydrate metabolism where sucrose decreased and fructans of higher degree of polymerization increased in leaves before and after the vegetative to reproductive transition were observed. At the molecular level, transcriptome analysis from leaf and shoot apical meristem tissue of A. tequilana plants from different developmental stages identified OASES as the most effective assembly program and revealed evidence for incomplete transcript processing in the highly redundant assembly obtained. Gene ontology analysis uncovered enrichment for terms associated with carbohydrate and hormone metabolism and detailed analysis of expression patterns for individual genes revealed roles for specific Flowering locus T (florigen), MADS box proteins, gibberellins and fructans in the transition to flowering. CONCLUSIONS: Based on the data obtained, a preliminary model was developed to describe the regulatory mechanisms underlying the initiation of flowering in A. tequilana. Identification of specific promoter and repressor Flowering Locus T and MADS box genes facilitates functional analysis and the development of strategies to modulate the vegetative to reproductive transition in A. tequilana.

Evaluating the efficacy of sample collection approaches and DNA metabarcoding for identifying the diversity of plants utilized by nectivorous bats.

In this study, we evaluated the efficacy of sample collection approaches and DNA metabarcoding to identify plants utilized by nectivorous bats. Samples included guano collected from beneath bat roosts and pollen-swabs from bat fur, both of which were subjected to DNA metabarcoding and visual identification of pollen (microscopy) to measure plant diversity. Our objectives were to determine whether DNA metabarcoding could detect likely food plants of nectivorous bats, whether sample types would produce different estimates of plant diversity, and to compare results of DNA metabarcoding to visual identification. Visual identification found that 99% of pollen was from Agave, which is thought to be the bats' main food source. The dominant taxon found by metabarcoding was also Agavoideae, but a broader diversity of plant species was also detected, many of which are likely "by-catch" from the broader environment. Metabarcoding outcomes differed between sample types, likely because pollen-swabs measured the plant species visited by bats and guano samples measured all items consumed in the bat's diet, even those that were not pollen or nectar. Overall, metabarcoding is a powerful, high-throughput tool to understand bat ecology and species interactions, but careful analysis of results is necessary to derive accurate ecological conclusions.

Diel rewiring and positive selection of ancient plant proteins enabled evolution of CAM photosynthesis in Agave.

BACKGROUND: Crassulacean acid metabolism (CAM) enhances plant water-use efficiency through an inverse day/night pattern of stomatal closure/opening that facilitates nocturnal CO2 uptake. CAM has evolved independently in over 35 plant lineages, accounting for ~ 6% of all higher plants. Agave species are highly heat- and drought-tolerant, and have been domesticated as model CAM crops for beverage, fiber, and biofuel production in semi-arid and arid regions. However, the genomic basis of evolutionary innovation of CAM in genus Agave is largely unknown. RESULTS: Using an approach that integrated genomics, gene co-expression networks, comparative genomics and protein structure analyses, we investigated the molecular evolution of CAM as exemplified in Agave. Comparative genomics analyses among C3, C4 and CAM species revealed that core metabolic components required for CAM have ancient genomic origins traceable to non-vascular plants while regulatory proteins required for diel re-programming of metabolism have a more recent origin shared among C3, C4 and CAM species. We showed that accelerated evolution of key functional domains in proteins responsible for primary metabolism and signaling, together with a diel re-programming of the transcription of genes involved in carbon fixation, carbohydrate processing, redox homeostasis, and circadian control is required for the evolution of CAM in Agave. Furthermore, we highlighted the potential candidates contributing to the adaptation of CAM functional modules. CONCLUSIONS: This work provides evidence of adaptive evolution of CAM related pathways. We showed that the core metabolic components required for CAM are shared by non-vascular plants, but regulatory proteins involved in re-reprogramming of carbon fixation and metabolite transportation appeared more recently. We propose that the accelerated evolution of key proteins together with a diel re-programming of gene expression were required for CAM evolution from C3 ancestors in Agave.

Retroelements and DNA Methylation Could Contribute to Diversity of 5S rDNA in Agave L.

Agave L. is a genus of economic importance, and many of the 166 species in the American plant genus Agave L. inhabit high-stress environments, which makes the genus promising for facing global climate change. However, sustainable use of economically important species without interfering with their ecology and evolution requires generating knowledge about the factors responsible for their genetic variation and diversity and, on this basis, their adaptation and speciation. Few genetic studies exploring the evolutionary relationships, speciation processes, genetic variability and diversity within species of Agave are currently available. Analyses of rDNA loci have been performed with the purpose of determining the genetic variability and diversity of the genus Agave, and these loci have been used as genetic markers of ploidy. However, the factors involved in the diversity of 5S rDNA regions in Agave have not yet been studied in depth. Our study explored the possible mechanisms of genetic (retroelements) and epigenetic (DNA methylation) diversity in 5S rDNA regions in Agave. We characterized the 5S rDNA gene tandem in species of the genus with different ploidy numbers and determined the levels of methylation in 13 haplotypes of 5S rDNA and in four non-transcribed spacers (NTS). Our results showed highly dynamic methylation with a high percentage in haplotypes and NTS of 5S rDNA regions in Agave. The characterization of the 5S rDNA tandem array in Agave revealed vestigial remains of the Cassandra terminal-repeat retrotransposon in miniature (TRIM). Our analysis supported previous results suggesting that in species of Agave L., regulation and diversity of 5S rDNA regions are controlled by coordinated genetic and epigenetic events, which will vary according to the species and the level of ploidy. The artificial pressure to which some agave crops are subjected may affect the mechanisms of evolution of gene 5S rDNA.

Population genetic analysis and bioclimatic modeling in Agave striata in the Chihuahuan Desert indicate higher genetic variation and lower differentiation in drier and more variable environments.

PREMISE OF THE STUDY: Is there an association between bioclimatic variables and genetic variation within species? This question can be approached by a detailed analysis of population genetics parameters along environmental gradients in recently originated species (so genetic drift does not further obscure the patterns). The genus Agave, with more than 200 recent species encompassing a diversity of morphologies and distributional patterns, is an adequate system for such analyses. We studied Agave striata, a widely distributed species from the Chihuahuan Desert, with a distinctive iteroparous reproductive ecology and two recognized subspecies with clear morphological differences. We used population genetic analyses along with bioclimatic studies to understand the effect of environment on the genetic variation and differentiation of this species. METHODS: We analyzed six populations of the subspecies A. striata subsp. striata, with a southern distribution, and six populations of A. striata subsp. falcata, with a northern distribution, using 48 ISSR loci and a total of 541 individuals (averaging 45 individuals per population). We assessed correlations between population genetics parameters (the levels of genetic variation and differentiation) and the bioclimatic variables of each population. We modeled each subspecies distribution and used linear correlations and multifactorial analysis of variance. KEY RESULTS: Genetic variation (measured as expected heterozygosity) increased at higher latitudes. Higher levels of genetic variation in populations were associated with a higher variation in environmental temperature and lower precipitation. Stronger population differentiation was associated with wetter and more variable precipitation in the southern distribution of the species. The two subspecies have genetic differences, which coincide with their climatic differences and potential distributions. CONCLUSIONS: Differences in genetic variation among populations and the genetic differentiation between A. striata subsp. striata and A. striata subsp. falcata is correlated with differences in environmental climatic variables along their distribution. We found two distinct gene pools that suggest active differentiation and perhaps incipient speciation. The detected association between genetic variation and environment variables indicates that climatic variables are playing an important role in the differentiation of A. striata.

Physiological differences and changes in global DNA methylation levels in Agave angustifolia Haw. albino variant somaclones during the micropropagation process.

KEY MESSAGE: Global DNA methylation changes caused by in vitro conditions are associated with the subculturing and phenotypic variation in Agave angustifolia Haw. While the relationship between the development of albinism and in vitro culture is well documented, the role of epigenetic processes in this development leaves some important questions unanswered. During the micropropagation of Agave angustifolia Haw., we found three different phenotypes, green (G), variegated (V) and albino (A). To understand the physiological and epigenetic differences among the somaclones, we analyzed several morphophysiological parameters and changes in the DNA methylation patterns in the three phenotypes during their in vitro development. We found that under in vitro conditions, the V plantlets maintained their CAM photosynthetic capacity, while the A variant showed no pigments and lost its CAM photosynthetic ability. Epigenetic analysis revealed that global DNA methylation increased in the G phenotype during the first two subcultures. However, after that time, DNA methylation levels declined. This hypomethylation correlated with the appearance of V shoots in the G plantlets. A similar correlation occurred in the V phenotype, where an increase of 2 % in the global DNA methylation levels was correlated with the generation of A shoots in the V plantlets. This suggests that an "epigenetic stress memory" during in vitro conditions causes a chromatin shift that favors the generation of variegated and albino shoots.

Functionally different PIN proteins control auxin flux during bulbil development in Agave tequilana.

In Agave tequilana, reproductive failure or inadequate flower development stimulates the formation of vegetative bulbils at the bracteoles, ensuring survival in a hostile environment. Little is known about the signals that trigger this probably unique phenomenon in agave species. Here we report that auxin plays a central role in bulbil development and show that the localization of PIN1-related proteins is consistent with altered auxin transport during this process. Analysis of agave transcriptome data led to the identification of the A. tequilana orthologue of PIN1 (denoted AtqPIN1) and a second closely related gene from a distinct clade reported as 'Sister of PIN1' (denoted AtqSoPIN1). Quantitative real-time reverse transcription-PCR (RT-qPCR) analysis showed different patterns of expression for each gene during bulbil formation, and heterologous expression of the A. tequilana PIN1 and SoPIN1 genes in Arabidopsis thaliana confirmed functional differences between these genes. Although no free auxin was detected in induced pedicel samples, changes in the levels of auxin precursors were observed. Taken as a whole, the data support the model that AtqPIN1 and AtqSoPIN1 have co-ordinated but distinct functions in relation to auxin transport during the initial stages of bulbil formation.

Agave biotechnology: an overview.

Agaves are plants of importance both in Mexican culture and economy and in other Latin-American countries. Mexico is reported to be the place of Agave origin, where today, scientists are looking for different industrial applications without compromising its sustainability and preserving the environment. To make it possible, a deep knowledge of all aspects involved in production process, agro-ecological management and plant biochemistry and physiology is required. Agave biotechnology research has been focusing on bio-fuels, beverages, foods, fibers, saponins among others. In this review, we present the advances and challenges of Agave biotechnology.

Genetic and morphological contrasts between wild and anthropogenic populations of Agave parryi var. huachucensis in south-eastern Arizona.

BACKGROUND AND AIMS: At least seven species of Agave, including A. parryi, were cultivated prehistorically in Arizona, serving as important sources of food and fibre. Many relict populations from ancient cultivation remain in the modern landscape, offering a unique opportunity to study pre-Columbian plant manipulation practices. This study examined genetic and morphological variation in six A. p. var. huachucensis populations of unknown origin to compare them with previous work on A. parryi populations of known origin, to infer their cultivation history and to determine whether artificial selection is evident in populations potentially managed by early agriculturalists. METHODS: Six A. p. var. huachucensis and 17 A. parryi populations were sampled, and morphometric, allozyme and microsatellite data were used to compare morphology and genetic structure in purportedly anthropogenic and wild populations, as well as in the two taxa. Analysis of molecular variance and Bayesian clustering were performed to partition variation associated with taxonomic identity and hypothesized evolutionary history, to highlight patterns of similarity among populations and to identify potential wild sources for the planting stock. KEY RESULTS: A p. var. huachucensis and A. parryi populations differed significantly both morphologically and genetically. Like A. parryi, wild A. p. var. huachucensis populations were more genetically diverse than the inferred anthropogenic populations, with greater expected heterozygosity, percentage of polymorphic loci and number of alleles. Inferred anthropogenic populations exhibited many traits indicative of past active cultivation: greater morphological uniformity, fixed heterozygosity for several loci (non-existent in wild populations), fewer multilocus genotypes and strong differentiation among populations. CONCLUSIONS: Where archaeological information is lacking, the genetic signature of many Agave populations in Arizona can be used to infer their evolutionary history and to identify potentially fruitful sites for archaeological investigation of ancient settlements and cultivation practices. The same approach can clearly be adopted for other species in similar situations.

Full-length agave transcriptome reveals candidate glycosyltransferase genes involved in hemicellulose biosynthesis.

Agave species are typical crassulacean acid metabolism (CAM) plants commonly cultivated to produce beverages, fibers, and medicines. To date, few studies have examined hemicellulose biosynthesis in Agave H11648, which is the primary cultivar used for fiber production. We conducted PacBio sequencing to obtain full-length transcriptome of five agave tissues: leaves, shoots, roots, flowers, and fruits. A total of 41,807 genes were generated, with a mean length of 2394 bp and an annotation rate of 97.12 % using public databases. We identified 42 glycosyltransferase genes related to hemicellulose biosynthesis, including mixed-linkage glucan (1), glucomannan (5), xyloglucan (16), and xylan (20). Their expression patterns were examined during leaf development and fungal infection, together with hemicellulose content. The results revealed four candidate glycosyltransferase genes involved in xyloglucan and xylan biosynthesis, including glucan synthase (CSLC), xylosyl transferase (XXT), xylan glucuronyltransferase (GUX), and xylan α-1,3-arabinosyltransferase (XAT). These genes can be potential targets for manipulating xyloglucan and xylan traits in agaves, and can also be used as candidate enzymatic tools for enzyme engineering. We have provided the first full-length transcriptome of agave, which will be a useful resource for gene identification and characterization in agave species. We also elucidated the hemicellulose biosynthesis machinery, which will benefit future studies on hemicellulose traits in agave.

De novo transcriptome assembly of drought tolerant CAM plants, Agave deserti and Agave tequilana.

BACKGROUND: Agaves are succulent monocotyledonous plants native to xeric environments of North America. Because of their adaptations to their environment, including crassulacean acid metabolism (CAM, a water-efficient form of photosynthesis), and existing technologies for ethanol production, agaves have gained attention both as potential lignocellulosic bioenergy feedstocks and models for exploring plant responses to abiotic stress. However, the lack of comprehensive Agave sequence datasets limits the scope of investigations into the molecular-genetic basis of Agave traits. RESULTS: Here, we present comprehensive, high quality de novo transcriptome assemblies of two Agave species, A. tequilana and A. deserti, built from short-read RNA-seq data. Our analyses support completeness and accuracy of the de novo transcriptome assemblies, with each species having a minimum of approximately 35,000 protein-coding genes. Comparison of agave proteomes to those of additional plant species identifies biological functions of gene families displaying sequence divergence in agave species. Additionally, a focus on the transcriptomics of the A. deserti juvenile leaf confirms evolutionary conservation of monocotyledonous leaf physiology and development along the proximal-distal axis. CONCLUSIONS: Our work presents a comprehensive transcriptome resource for two Agave species and provides insight into their biology and physiology. These resources are a foundation for further investigation of agave biology and their improvement for bioenergy development.

Molecular variability among isolates of Fusarium oxysporum associated with root rot disease of Agave tequilana.

In this study, 115 isolates of Fusarium oxysporum from roots of Agave tequilana Weber cv azul plants and soil in commercial plantations in western Mexico were characterized using morphological and molecular methods. Genetic analyses of monosporic isolates included restriction enzyme analysis of rDNA (ARDRA) using HaeIII and HinfI, and genetic diversity was determined using Box-PCR molecular markers. Box-PCR analysis generated 14 groups. The groups correlated highly with the geographic location of the isolate and sample type. These results demonstrate the usefulness of ARDRA and Box-PCR techniques in the molecular characterization of the Fusarium genus for the discrimination of pathogenic isolates.

Differential Repeat Accumulation in the Bimodal Karyotype of Agave L.

The genus Agave presents a bimodal karyotype with x = 30 (5L, large, +25S, small chromosomes). Bimodality within this genus is generally attributed to allopolyploidy in the ancestral form of Agavoideae. However, alternative mechanisms, such as the preferential accumulation of repetitive elements at the macrochromosomes, could also be important. Aiming to understand the role of repetitive DNA within the bimodal karyotype of Agave, genomic DNA from the commercial hybrid 11648 (2n = 2x = 60, 6.31 Gbp) was sequenced at low coverage, and the repetitive fraction was characterized. In silico analysis showed that ~67.6% of the genome is mainly composed of different LTR retrotransposon lineages and one satellite DNA family (AgSAT171). The satellite DNA localized at the centromeric regions of all chromosomes; however, stronger signals were observed for 20 of the macro- and microchromosomes. All transposable elements showed a dispersed distribution, but not uniform across the length of the chromosomes. Different distribution patterns were observed for different TE lineages, with larger accumulation at the macrochromosomes. The data indicate the differential accumulation of LTR retrotransposon lineages at the macrochromosomes, probably contributing to the bimodality. Nevertheless, the differential accumulation of the satDNA in one group of macro- and microchromosomes possibly reflects the hybrid origin of this Agave accession.

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.