Phylogenetically diverse wild plant species use common biochemical strategies to thrive in the Atacama Desert.

The best ideotypes are under mounting pressure due to increased aridity. Understanding the conserved molecular mechanisms that evolve in wild plants adapted to harsh environments is crucial in developing new strategies for agriculture. Yet our knowledge of such mechanisms in wild species is scant. We performed metabolic pathway reconstruction using transcriptome information from 32 Atacama and phylogenetically related species that do not live in Atacama (sister species). We analyzed reaction enrichment to understand the commonalities and differences of Atacama plants. To gain insights into the mechanisms that ensure survival, we compared expressed gene isoform numbers and gene expression patterns between the annotated biochemical reactions from 32 Atacama and sister species. We found biochemical convergences characterized by reactions enriched in at least 50% of the Atacama species, pointing to potential advantages against drought and nitrogen starvation, for instance. These findings suggest that the adaptation in the Atacama Desert may result in part from shared genetic legacies governing the expression of key metabolic pathways to face harsh conditions. Enriched reactions corresponded to ubiquitous compounds common to extreme and agronomic species and were congruent with our previous metabolomic analyses. Convergent adaptive traits offer promising candidates for improving abiotic stress resilience in crop species.

Plant ecological genomics at the limits of life in the Atacama Desert.

The Atacama Desert in Chile-hyperarid and with high-ultraviolet irradiance levels-is one of the harshest environments on Earth. Yet, dozens of species grow there, including Atacama-endemic plants. Herein, we establish the Talabre-Lejía transect (TLT) in the Atacama as an unparalleled natural laboratory to study plant adaptation to extreme environmental conditions. We characterized climate, soil, plant, and soil-microbe diversity at 22 sites (every 100 m of altitude) along the TLT over a 10-y period. We quantified drought, nutrient deficiencies, large diurnal temperature oscillations, and pH gradients that define three distinct vegetational belts along the altitudinal cline. We deep-sequenced transcriptomes of 32 dominant plant species spanning the major plant clades, and assessed soil microbes by metabarcoding sequencing. The top-expressed genes in the 32 Atacama species are enriched in stress responses, metabolism, and energy production. Moreover, their root-associated soils are enriched in growth-promoting bacteria, including nitrogen fixers. To identify genes associated with plant adaptation to harsh environments, we compared 32 Atacama species with the 32 closest sequenced species, comprising 70 taxa and 1,686,950 proteins. To perform phylogenomic reconstruction, we concatenated 15,972 ortholog groups into a supermatrix of 8,599,764 amino acids. Using two codon-based methods, we identified 265 candidate positively selected genes (PSGs) in the Atacama plants, 64% of which are located in Pfam domains, supporting their functional relevance. For 59/184 PSGs with an Arabidopsis ortholog, we uncovered functional evidence linking them to plant resilience. As some Atacama plants are closely related to staple crops, these candidate PSGs are a "genetic goldmine" to engineer crop resilience to face climate change.

Shotgun proteomics of peach fruit reveals major metabolic pathways associated to ripening.

BACKGROUND: Fruit ripening in Prunus persica melting varieties involves several physiological changes that have a direct impact on the fruit organoleptic quality and storage potential. By studying the proteomic differences between the mesocarp of mature and ripe fruit, it would be possible to highlight critical molecular processes involved in the fruit ripening. RESULTS: To accomplish this goal, the proteome from mature and ripe fruit was assessed from the variety O'Henry through shotgun proteomics using 1D-gel (PAGE-SDS) as fractionation method followed by LC/MS-MS analysis. Data from the 131,435 spectra could be matched to 2740 proteins, using the peach genome reference v1. After data pre-treatment, 1663 proteins could be used for comparison with datasets assessed using transcriptomic approaches and for quantitative protein accumulation analysis. Close to 26% of the genes that code for the proteins assessed displayed higher expression at ripe fruit compared to other fruit developmental stages, based on published transcriptomic data. Differential accumulation analysis between mature and ripe fruit revealed that 15% of the proteins identified were modulated by the ripening process, with glycogen and isocitrate metabolism, and protein localization overrepresented in mature fruit, as well as cell wall modification in ripe fruit. Potential biomarkers for the ripening process, due to their differential accumulation and gene expression pattern, included a pectin methylesterase inhibitor, a gibbellerin 2-beta-dioxygenase, an omega-6 fatty acid desaturase, a homeobox-leucine zipper protein and an ACC oxidase. Transcription factors enriched in NAC and Myb protein domains would target preferentially the genes encoding proteins more abundant in mature and ripe fruit, respectively. CONCLUSIONS: Shotgun proteomics is an unbiased approach to get deeper into the proteome allowing to detect differences in protein abundance between samples. This technique provided a resolution so that individual gene products could be identified. Many proteins likely involved in cell wall and sugar metabolism, aroma and color, change their abundance during the transition from mature to ripe fruit.

A Prunus persica genome-wide RNA-seq approach uncovers major differences in the transcriptome among chilling injury sensitive and non-sensitive varieties.

Chilling injury represents a major constrain for crops productivity. Prunus persica, one of the most relevant rosacea crops, have early season varieties that are resistant to chilling injury, in contrast to late season varieties, which display chilling symptoms such as mealiness (dry, sandy fruit mesocarp) after prolonged storage at chilling temperatures. To uncover the molecular processes related to the ability of early varieties to withstand mealiness, postharvest and genome-wide RNA-seq assessments were performed in two early and two late varieties. Differences in juice content and ethylene biosynthesis were detected among early and late season fruits that became mealy after exposed to prolonged chilling. Principal component and data distribution analysis revealed that cold-stored late variety fruit displayed an exacerbated and unique transcriptome profile when compared to any other postharvest condition. A differential expression analysis performed using an empirical Bayes mixture modeling approach followed by co-expression and functional enrichment analysis uncover processes related to ethylene, lipids, cell wall, carotenoids and DNA metabolism, light response, and plastid homeostasis associated to the susceptibility or resistance of P. persica varieties to chilling stress. Several of the genes related to these processes are in quantitative trait loci (QTL) associated to mealiness in P. persica. Together, these analyses exemplify how P. persica can be used as a model for studying chilling stress in plants.

Genome-Wide Prediction of Metabolic Enzymes, Pathways, and Gene Clusters in Plants.

Plant metabolism underpins many traits of ecological and agronomic importance. Plants produce numerous compounds to cope with their environments but the biosynthetic pathways for most of these compounds have not yet been elucidated. To engineer and improve metabolic traits, we need comprehensive and accurate knowledge of the organization and regulation of plant metabolism at the genome scale. Here, we present a computational pipeline to identify metabolic enzymes, pathways, and gene clusters from a sequenced genome. Using this pipeline, we generated metabolic pathway databases for 22 species and identified metabolic gene clusters from 18 species. This unified resource can be used to conduct a wide array of comparative studies of plant metabolism. Using the resource, we discovered a widespread occurrence of metabolic gene clusters in plants: 11,969 clusters from 18 species. The prevalence of metabolic gene clusters offers an intriguing possibility of an untapped source for uncovering new metabolite biosynthesis pathways. For example, more than 1,700 clusters contain enzymes that could generate a specialized metabolite scaffold (signature enzymes) and enzymes that modify the scaffold (tailoring enzymes). In four species with sufficient gene expression data, we identified 43 highly coexpressed clusters that contain signature and tailoring enzymes, of which eight were characterized previously to be functional pathways. Finally, we identified patterns of genome organization that implicate local gene duplication and, to a lesser extent, single gene transposition as having played roles in the evolution of plant metabolic gene clusters.

Patterns of metabolite changes identified from large-scale gene perturbations in Arabidopsis using a genome-scale metabolic network.

Metabolomics enables quantitative evaluation of metabolic changes caused by genetic or environmental perturbations. However, little is known about how perturbing a single gene changes the metabolic system as a whole and which network and functional properties are involved in this response. To answer this question, we investigated the metabolite profiles from 136 mutants with single gene perturbations of functionally diverse Arabidopsis (Arabidopsis thaliana) genes. Fewer than 10 metabolites were changed significantly relative to the wild type in most of the mutants, indicating that the metabolic network was robust to perturbations of single metabolic genes. These changed metabolites were closer to each other in a genome-scale metabolic network than expected by chance, supporting the notion that the genetic perturbations changed the network more locally than globally. Surprisingly, the changed metabolites were close to the perturbed reactions in only 30% of the mutants of the well-characterized genes. To determine the factors that contributed to the distance between the observed metabolic changes and the perturbation site in the network, we examined nine network and functional properties of the perturbed genes. Only the isozyme number affected the distance between the perturbed reactions and changed metabolites. This study revealed patterns of metabolic changes from large-scale gene perturbations and relationships between characteristics of the perturbed genes and metabolic changes.

Identification of Metabolite and Lipid Profiles in a Segregating Peach Population Associated with Mealiness in Prunus persica (L.) Batsch.

The peach is the third most important temperate fruit crop considering fruit production and harvested area in the world. Exporting peaches represents a challenge due to the long-distance nature of export markets. This requires fruit to be placed in cold storage for a long time, which can induce a physiological disorder known as chilling injury (CI). The main symptom of CI is mealiness, which is perceived as non-juicy fruit by consumers. The purpose of this work was to identify and compare the metabolite and lipid profiles between two siblings from contrasting populations for juice content, at harvest and after 30 days at 0 °C. A total of 119 metabolites and 189 lipids were identified, which showed significant differences in abundance, mainly in amino acids, sugars and lipids. Metabolites displaying significant changes from the E1 to E3 stages corresponded to lipids such as phosphatidylglycerol (PG), monogalactosyldiacylglycerol (MGDG) and lysophosphatidylcholines (LPC), and sugars such as fructose 1 and 1-fructose-6 phosphate. These metabolites might be used as early stage biomarkers associated with mealiness at harvest and after cold storage.

Transcriptome analysis during ripening of table grape berry cv. Thompson Seedless.

Ripening is one of the key processes associated with the development of major organoleptic characteristics of the fruit. This process has been extensively characterized in climacteric fruit, in contrast with non-climacteric fruit such as grape, where the process is less understood. With the aim of studying changes in gene expression during ripening of non-climacteric fruit, an Illumina based RNA-Seq transcriptome analysis was performed on four developmental stages, between veraison and harvest, on table grapes berries cv Thompson Seedless. Functional analysis showed a transcriptional increase in genes related with degradation processes of chlorophyll, lipids, macromolecules recycling and nucleosomes organization; accompanied by a decrease in genes related with chloroplasts integrity and amino acid synthesis pathways. It was possible to identify several processes described during leaf senescence, particularly close to harvest. Before this point, the results suggest a high transcriptional activity associated with the regulation of gene expression, cytoskeletal organization and cell wall metabolism, which can be related to growth of berries and firmness loss characteristic to this stage of development. This high metabolic activity could be associated with an increase in the transcription of genes related with glycolysis and respiration, unexpected for a non-climacteric fruit ripening.

Genomic signatures of specialized metabolism in plants.

All plants synthesize basic metabolites needed for survival (primary metabolism), but different taxa produce distinct metabolites that are specialized for specific environmental interactions (specialized metabolism). Because evolutionary pressures on primary and specialized metabolism differ, we investigated differences in the emergence and maintenance of these processes across 16 species encompassing major plant lineages from algae to angiosperms. We found that, relative to their primary metabolic counterparts, genes coding for specialized metabolic functions have proliferated to a much greater degree and by different mechanisms and display lineage-specific patterns of physical clustering within the genome and coexpression. These properties illustrate the differential evolution of specialized metabolism in plants, and collectively they provide unique signatures for the potential discovery of novel specialized metabolic processes.