Unraveling the genetic complexity of a cultivated breeding population of "yerba mate" (Ilex paraguariensis St. Hil.).

We report the first comprehensive multilocus molecular characterization of cultivated Ilex paraguariensis plants belonging to a breeding program. Using nuclear and homologous chloroplastidic molecular markers, we have genotyped 158 plants from four plantation sites. Analyses of the nuclear data (187 variable dominant loci) allowed detecting high diversity (0.569), the occurrence of four distinct genetic clusters, and a low but significant differentiation among sites. Additionally, 20 chloroplastidic alleles were identified applying five microsatellite polymorphic markers, and a high chloroplastidic diversity was recognized (0.505); two haplogroups were distinguished amongst the 63 haplotypes detected. Our results from both nuclear and plastidic markers indicate that most genetic variation reside within plantations sites (≥ 95%), and that these plantations were established on highly variable materials (either as seeds or plantlets) derived from, at least, 63 maternal lineages. Moreover, our study suggests that the genetic structure within each plantation site was most likely shaped by past admixture favored by farmers´ practices during the establishment of each plantation. Also, subsequent constraints in gene flow and/or a low level of shared polymorphism among plantations could have contributed to current structure.

Transcript and metabolic adjustments triggered by drought in Ilex paraguariensis leaves.

MAIN CONCLUSION: Abscisic acid is involved in the drought response of Ilex paraguariensis. Acclimation includes root growth stimulation, stomatal closure, osmotic adjustment, photoprotection, and regulation of nonstructural carbohydrates and amino acid metabolisms. Ilex paraguariensis (yerba mate) is cultivated in the subtropical region of South America, where the occurrence of drought episodes limit yield. To explore the mechanisms that allow I. paraguariensis to overcome dehydration, we investigated (1) how gene expression varied between water-stressed and non-stressed plants and (2) in what way the modulation of gene expression was linked to physiological status and metabolite composition. A total of 4920 differentially expressed transcripts were obtained through RNA-Seq after water deprivation. Drought induced the expression of several transcripts involved in the ABA-signalling pathway. Stomatal closure and leaf osmotic adjustments were promoted to minimize water loss, and these responses were accompanied by a high transcriptional remodeling of stress perception, signalling and transcriptional regulation, the photoprotective and antioxidant systems, and other stress-responsive genes. Simultaneously, significant changes in metabolite contents were detected. Glutamine, phenylalanine, isomaltose, fucose, and malate levels were shown to be positively correlated with dehydration. Principal component analysis showed differences in the metabolic profiles of control and stressed leaves. These results provide a comprehensive overview of how I. paraguariensis responds to dehydration at transcriptional and metabolomic levels and provide further characterization of the molecular mechanisms associated with drought response in perennial subtropical species.

Yerba mate (Ilex paraguariensis, A. St.-Hil.) de novo transcriptome assembly based on tissue specific genomic expression profiles.

BACKGROUND: The most common infusion in southern Latin-American countries is prepared with dried leaves of Ilex paraguariensis A. St.-Hil., an aboriginal ancestral beverage known for its high polyphenols concentration currently consumed in > 90% of homes in Argentina, in Paraguay and Uruguay. The economy of entire provinces heavily relies on the production, collection and manufacture of Ilex paraguariensis, the fifth plant species with highest antioxidant activity. Polyphenols are associated to relevant health benefits including strong antioxidant properties. Despite its regional relevance and potential biotechnological applications, little is known about functional genomics and genetics underlying phenotypic variation of relevant traits. By generating tissue specific transcriptomic profiles, we aimed to comprehensively annotate genes in the Ilex paraguariensis phenylpropanoid pathway and to evaluate differential expression profiles. RESULTS: In this study we generated a reliable transcriptome assembly based on a collection of 15 RNA-Seq libraries from different tissues of Ilex paraguariensis. A total of 554 million RNA-Seq reads were assembled into 193,897 transcripts, where 24,612 annotated full-length transcripts had complete ORF. We assessed the transcriptome assembly quality, completeness and accuracy using BUSCO and TransRate; consistency was also evaluated by experimentally validating 11 predicted genes by PCR and sequencing. Functional annotation against KEGG Pathway database identified 1395 unigenes involved in biosynthesis of secondary metabolites, 531 annotated transcripts corresponded to the phenylpropanoid pathway. The top 30 differentially expressed genes among tissue revealed genes involved in photosynthesis and stress response. These significant differences were then validated by qRT-PCR. CONCLUSIONS: Our study is the first to provide data from whole genome gene expression profiles in different Ilex paraguariensis tissues, experimentally validating in-silico predicted genes key to the phenylpropanoid (antioxidant) pathway. Our results provide essential genomic data of potential use in breeding programs for polyphenol content. Further studies are necessary to assess if the observed expression variation in the phenylpropanoid pathway annotated genes is related to variations in leaves' polyphenol content at the population scale. These results set the current reference for Ilex paraguariensis genomic studies and provide a substantial contribution to research and biotechnological applications of phenylpropanoid secondary metabolites.

Genetic and Phytochemical Analysis to Evaluate the Diversity and Relationships of Mate (Ilex paraguariensis A.St.-Hil.) Elite Genetic Resources in a Germplasm Collection.

The aim of this study was to evaluate the phytochemical and genetic diversity, relationships and identification of mate (Ilex paraguariensis A.St.-Hil.) elite genetic resources belonging to the Brazilian germplasm collection and mate breeding program. Mate has been studied due to the presence of phytochemical compounds, especially methylxanthines and phenolic compounds. The samples were collected from the leaves of 76 mate elite genetic resources (16 progenies × 5 localities). Total DNA was extracted from mate leaves and 20 random primers were used for DNA amplification. Methylxanthines (caffeine and theobromine) and phenolic compounds (chlorogenic, neochlorogenic, and criptochlorogenic acids) were quantified by HPLC. The genetic divergence estimated was higher within (92%) than among (8%) the different populations. Analysis of genetic distance between origins provided the formation of two groups by UPGMA cluster analysis, with higher polymorphism (94.9%). The average content of caffeine ranged from 0.01 to 1.38% and theobromine of 0.10 - 0.85% (w/w). The caffeoylquinic acids concentrations (1.43 - 5.38%) showed a gradient 3-CQA > 5-CQA > 4-CQA. The coefficient of genetic variation (CVg) was of low magnitude for all mono-caffeoylquinics acids. Significant correlations (positive and negative) were observed between the phytochemical compounds. Genetic diversity analysis performed by RAPD markers showed a greater intra-populational diversity; genetic resources with low caffeine and higher theobromine content were identified and can be used in breeding programs; the correlation between methylxanthines and phenolic compounds can be used as a good predictor in future studies.

Genetic diversity of wild germplasm of "yerba mate" (Ilex paraguariensis St. Hil.) from Uruguay.

The "yerba mate" tree, Ilex paraguariensis St. Hil., is a crop native to subtropical South America, marketed for the elaboration of the highly popular "mate" beverage. The Uruguayan germplasm occupies the southernmost area of the species distribution range and carries adaptations to environments that considerably differ from the current production area. We characterized the genetic variability of the germplasm from this unexplored area by jointly analyzing individuals from the diversification center (ABP, Argentina, Brazil and Paraguay) with 19 nuclear and 11 plastidic microsatellite markers. For the Uruguayan germplasm, we registered 55 alleles (18 % private), and 80 genotypes (44 % exclusive), whereas 63 alleles (28.6 % private) and 81 genotypes (42 % exclusive) were recorded for individuals from ABP. Only two plastidic haplotypes were detected. Distance-based and multilocus genotype analyses showed that individuals from ABP intermingle and that the Uruguayan germplasm is differentiated in three gene-pools. Significant positive correlations between genetic and geographic distances were detected. Our results concur in that ABP individuals harbor greater genetic variation than those from the tail of the distribution, as to the number of alleles (1.15-fold), He (1.19-fold), Rs (1.39-fold), and the between-group genetic distances (1.16-fold). Also the shape of the genetic landscape interpolation analysis suggests that the genetic variation decays southward towards the Uruguayan territory. We showed that Uruguayan germplasm hosts a combination of nuclear alleles not present in the central region, constituting a valuable breeding resource. Future conservation efforts should concentrate in collecting numerous individuals of "yerba mate" per site to gather the existent variation.

Exploring the genes of yerba mate (Ilex paraguariensis A. St.-Hil.) by NGS and de novo transcriptome assembly.

Yerba mate (Ilex paraguariensis A. St.-Hil.) is an important subtropical tree crop cultivated on 326,000 ha in Argentina, Brazil and Paraguay, with a total yield production of more than 1,000,000 t. Yerba mate presents a strong limitation regarding sequence information. The NCBI GenBank lacks an EST database of yerba mate and depicts only 80 DNA sequences, mostly uncharacterized. In this scenario, in order to elucidate the yerba mate gene landscape by means of NGS, we explored and discovered a vast collection of I. paraguariensis transcripts. Total RNA from I. paraguariensis was sequenced by Illumina HiSeq-2000 obtaining 72,031,388 pair-end 100 bp sequences. High quality reads were de novo assembled into 44,907 transcripts encompassing 40 million bases with an estimated coverage of 180X. Multiple sequence analysis allowed us to predict that yerba mate contains ∼ 32,355 genes and 12,551 gene variants or isoforms. We identified and categorized members of more than 100 metabolic pathways. Overall, we have identified ∼ 1,000 putative transcription factors, genes involved in heat and oxidative stress, pathogen response, as well as disease resistance and hormone response. We have also identified, based in sequence homology searches, novel transcripts related to osmotic, drought, salinity and cold stress, senescence and early flowering. We have also pinpointed several members of the gene silencing pathway, and characterized the silencing effector Argonaute1. We predicted a diverse supply of putative microRNA precursors involved in developmental processes. We present here the first draft of the transcribed genomes of the yerba mate chloroplast and mitochondrion. The putative sequence and predicted structure of the caffeine synthase of yerba mate is presented. Moreover, we provide a collection of over 10,800 SSR accessible to the scientific community interested in yerba mate genetic improvement. This contribution broadly expands the limited knowledge of yerba mate genes, and is presented as the first genomic resource of this important crop.

A succinate dehydrogenase flavoprotein subunit-like transcript is upregulated in Ilex paraguariensis leaves in response to water deficit and abscisic acid.

Ilex paraguariensis plants were subjected to progressive soil water deficit, and differential display (DD) was used to analyse gene expression in leaves to characterise physiological responses to mild and severe water deficits. A cDNA fragment showing strong homology with the flavoprotein subunit (SDH1) of succinate:ubiquinone oxidoreductase (succinate dehydrogenase, SDH, EC 1.3.5.1) was upregulated in plants exposed to drought. Quantitative real-time PCR revealed that the SDH1-like transcript level began to increase when the leaf relative water content (RWC) decreased to 78% and peaked when the RWC dropped to 57%. A correlation between abscisic acid (ABA) concentration and variations in transcript levels was assessed by GC-SIM. After rehydration, SDH1 mRNA and ABA returned to their initial levels. In stressed leaves sprayed with ABA SDH1 mRNA accumulated in greater levels compared to stressed leaves that did not receive ABA. Moreover, the enzymatic activity of succinate dehydrogenase increased 1.5-fold in the mature leaves of ABA-treated plants. This physiological response may be related to the tendency of this species to minimise water losses through stomatal closure in the early stages of dehydration to avoid tissue desiccation. As the leaf water potential diminished due to an increase in water restriction, I. paraguariensis leaf tissues reacted by making osmotic adjustments to sustain tissue metabolic activity, which enables the recovery of photosynthesis upon re-watering. These results provide new insights concerning the linkage between plant respiration and photosynthetic metabolism that could be potentially further used in breeding programs aiming water tolerant genotypes.

Analysis of free-radical scavenging of Yerba Mate (Ilex paraguriensis) using electron spin resonance and radical-induced DNA damage.

Mate (MT) is a popular South American beverage that has been used as a traditional medicine for centuries, spurring recent interest in its nutraceutical properties. MT is prepared as an infusion of leaves from the Yerba Mate (llex paraguriensis) tree. MT has been reported to have antioxidant properties in vitro and in vivo, but these have not been fully characterized in terms of effects against specific radicals. Accordingly, we examined the antioxidant effects of an MT infusion against hydroxyl and superoxide radicals in both chemical and cell culture assays. MT infusions were prepared at 3.10 g/L in boiling water and diluted to experimental dilutions from this stock. Electron spin resonance (ESR) experiments indicated that MT scavenged hydroxyl radicals (produced via the Fenton reaction) and superoxide radicals (produced via the xanthine/xanthine oxidase enzymatic reaction) at all concentrations tested (P < 0.05). Further controls indicated that superoxide radical scavenging was not due to xanthine oxidase inhibition. MT scavenged hydroxyl radicals and decreased cellular oxygen consumption in a dose-dependent manner in Cr(VI)-stimulated RAW 264.7 cells, based on ESR and oxygraph measurements (P < 0.05). Similarly, MT also inhibited hydroxyl-radical-induced lipid peroxidation and DNA damage in a dose-dependent manner in RAW 264.7 cells, based on malondialdehyde and Comet assay data (P < 0.05). This study indicates that MT possesses potent antioxidant effects against hydroxyl and superoxide radicals in both chemical and cell culture systems, as well as DNA-protective properties. These data further clarify the reported antioxidant effects of Yerba Mate infusions.

Genomic screening in dioecious "yerba mate" tree (Ilex paraguariensis A. St. Hill., Aquifoliaceae) through representational difference analysis.

The "yerba mate" tree, Ilex paraguariensis, is a functionally dioecious crop species with economic relevance in several South American countries. We report a genomic screening accomplished through representational difference analysis (RDA) in male and female I. paraguariensis trees. The aim of the present paper was to investigate the occurrence of sex-related genomic differences in order to develop an early gender detection molecular method that could help reducing energy inputs during the "yerba mate" processing and that could be suitable for breeding programs. An intra-experiment redundancy was detected via SSCP analysis and sequence characterization. Taking together both reciprocal RDA assays, fragments isolated can be discriminated into three main categories. The first category of fragments shows spurious affinities with available deposited sequences and could be considered as specific to I. paraguariensis. The second category comprises sequences identified as organellar or ribosomal plant DNA. Sequences grouped in the third category involve clones akin to conserved domains of retrotransposons (RNaseH, integrases and/or chromodomains) from at least two distinct lineages of Ty3/Gypsy retrotransposons and one from Ty1/Copia retroelements, which in addition are associated to sex determination regions of the Solanaceae, Caricaceae and Salicaceae. A contig sequence was assembled that codes for an integrase core domain and a chromodomain. The phylogenetic analysis of the so-called IPRE (for I. paraguariensis retroelement) integrase domain indicates that it belongs to the Del lineage of the Chromoviridae. This is the first report of mobile elements isolated and detected from the "yerba mate" tree. Although RDA derived fragments, so far tested, have been retrieved from both sexes with similar sequences, association to sex related regions cannot be completely discarded. Implications of present results are further discussed.