Coding and long non-coding RNAs provide evidence of distinct transcriptional reprogramming for two ecotypes of the extremophile plant Eutrema salsugineum undergoing water deficit stress.

BACKGROUND: The severity and frequency of drought has increased around the globe, creating challenges in ensuring food security for a growing world population. As a consequence, improving water use efficiency by crops has become an important objective for crop improvement. Some wild crop relatives have adapted to extreme osmotic stresses and can provide valuable insights into traits and genetic signatures that can guide efforts to improve crop tolerance to water deficits. Eutrema salsugineum, a close relative of many cruciferous crops, is a halophytic plant and extremophyte model for abiotic stress research. RESULTS: Using comparative transcriptomics, we show that two E. salsugineum ecotypes display significantly different transcriptional responses towards a two-stage drought treatment. Even before visibly wilting, water deficit led to the differential expression of almost 1,100 genes for an ecotype from the semi-arid, sub-arctic Yukon, Canada, but only 63 genes for an ecotype from the semi-tropical, monsoonal, Shandong, China. After recovery and a second drought treatment, about 5,000 differentially expressed genes were detected in Shandong plants versus 1,900 genes in Yukon plants. Only 13 genes displayed similar drought-responsive patterns for both ecotypes. We detected 1,007 long non-protein coding RNAs (lncRNAs), 8% were only expressed in stress-treated plants, a surprising outcome given the documented association between lncRNA expression and stress. Co-expression network analysis of the transcriptomes identified eight gene clusters where at least half of the genes in each cluster were differentially expressed. While many gene clusters were correlated to drought treatments, only a single cluster significantly correlated to drought exposure in both ecotypes. CONCLUSION: Extensive, ecotype-specific transcriptional reprogramming with drought was unexpected given that both ecotypes are adapted to saline habitats providing persistent exposure to osmotic stress. This ecotype-specific response would have escaped notice had we used a single exposure to water deficit. Finally, the apparent capacity to improve tolerance and growth after a drought episode represents an important adaptive trait for a plant that thrives under semi-arid Yukon conditions, and may be similarly advantageous for crop species experiencing stresses attributed to climate change.

Evidence that tolerance of Eutrema salsugineum to low phosphate conditions is hard-wired by constitutive metabolic and root-associated adaptations.

MAIN CONCLUSION: The extremophyte Eutrema salsugineum (Yukon ecotype) has adapted to an environment low in available phosphate through metabolic and root-associated traits that enables it to efficiently retrieve, use, and recycle phosphorus. Efficient phosphate (Pi) use by plants would increase crop productivity under Pi-limiting conditions and reduce our reliance on Pi applied as fertilizer. An ecotype of Eutrema salsugineum originating from the Yukon, Canada, shows no evidence of decreased relative growth rate or biomass under low Pi conditions and, as such, offers a promising model for identifying mechanisms to improve Pi use by crops. We evaluated traits associated with efficient Pi use by Eutrema (Yukon ecotype) seedlings and 4-week-old plants, including acquisition, remobilization, and the operation of metabolic bypasses. Relative to Arabidopsis, Eutrema was slower to remobilize phosphorus (P) from senescing leaves, primary and lateral roots showed a lower capacity for rhizosphere acidification, and root acid phosphatase activity was more broadly distributed and not Pi responsive. Both species produced long root hairs on low Pi media, whereas Arabidopsis root hairs were well endowed with phosphatase activity. This capacity was largely absent in Eutrema. In contrast to Arabidopsis, maximal in vitro rates of pyrophosphate-dependent phosphofructokinase and phosphoenolpyruvate carboxylase activities were not responsive to low Pi conditions suggesting that Eutrema has a constitutive and likely preferential capacity to use glycolytic bypass enzymes. Rhizosphere acidification, exudation of acid phosphatases, and rapid remobilization of leaf P are unlikely strategies used by Eutrema for coping with low Pi. Rather, equipping an entire root system for Pi acquisition and utilizing a metabolic strategy suited to deficient Pi conditions offer better explanations for how Eutrema has adapted to thrive on alkaline, highly saline soil that is naturally low in available Pi.

Acclimation of the crucifer Eutrema salsugineum to phosphate limitation is associated with constitutively high expression of phosphate-starvation genes.

Eutrema salsugineum, a halophytic relative of Arabidopsis thaliana, was subjected to varying phosphate (Pi) treatments. Arabidopsis seedlings grown on 0.05 mm Pi displayed shortened primary roots, higher lateral root density and reduced shoot biomass allocation relative to those on 0.5 mm Pi, whereas Eutrema seedlings showed no difference in lateral root density and shoot biomass allocation. While a low Fe concentration mitigated the Pi deficiency response for Arabidopsis, Eutrema root architecture was unaltered, but adding NaCl increased Eutrema lateral root density almost 2-fold. Eutrema and Arabidopsis plants grown on soil without added Pi for 4 weeks had low shoot and root Pi content. Pi-deprived, soil-grown Arabidopsis plants were stunted with senescing older leaves, whereas Eutrema plants were visually indistinguishable from 2.5 mm Pi-supplemented plants. Genes associated with Pi starvation were analysed by RT-qPCR. EsIPS2, EsPHT1;4 and EsPAP17 showed up-regulated expression in Pi-deprived Eutrema plants, while EsPHR1, EsWRKY75 and EsRNS1 showed no induction. Absolute quantification of transcripts indicated that PHR1, WRKY75 and RNS1 were expressed at higher levels in Eutrema plants relative to those in Arabidopsis regardless of external Pi. The low phenotypic plasticity Eutrema displays to Pi supply is consistent with adaptation to chronic Pi deprivation in its extreme natural habitat.

RNA-Seq effectively monitors gene expression in Eutrema salsugineum plants growing in an extreme natural habitat and in controlled growth cabinet conditions.

BACKGROUND: The investigation of extremophile plant species growing in their natural environment offers certain advantages, chiefly that plants adapted to severe habitats have a repertoire of stress tolerance genes that are regulated to maximize plant performance under physiologically challenging conditions. Accordingly, transcriptome sequencing offers a powerful approach to address questions concerning the influence of natural habitat on the physiology of an organism. We used RNA sequencing of Eutrema salsugineum, an extremophile relative of Arabidopsis thaliana, to investigate the extent to which genetic variation and controlled versus natural environments contribute to differences between transcript profiles. RESULTS: Using 10 million cDNA reads, we compared transcriptomes from two natural Eutrema accessions (originating from Yukon Territory, Canada and Shandong Province, China) grown under controlled conditions in cabinets and those from Yukon plants collected at a Yukon field site. We assessed the genetic heterogeneity between individuals using single-nucleotide polymorphisms (SNPs) and the expression patterns of 27,016 genes. Over 39,000 SNPs distinguish the Yukon from the Shandong accessions but only 4,475 SNPs differentiated transcriptomes of Yukon field plants from an inbred Yukon line. We found 2,989 genes that were differentially expressed between the three sample groups and multivariate statistical analyses showed that transcriptomes of individual plants from a Yukon field site were as reproducible as those from inbred plants grown under controlled conditions. Predicted functions based upon gene ontology classifications show that the transcriptomes of field plants were enriched by the differential expression of light- and stress-related genes, an observation consistent with the habitat where the plants were found. CONCLUSION: Our expectation that comparative RNA-Seq analysis of transcriptomes from plants originating in natural habitats would be confounded by uncontrolled genetic and environmental factors was not borne out. Moreover, the transcriptome data shows little genetic variation between laboratory Yukon Eutrema plants and those found at a field site. Transcriptomes were reproducible and biological associations meaningful whether plants were grown in cabinets or found in the field. Thus RNA-Seq is a valuable approach to study native plants in natural environments and this technology can be exploited to discover new gene targets for improved crop performance under adverse conditions.

Transcriptomic and metabolomic analysis of Yukon Thellungiella plants grown in cabinets and their natural habitat show phenotypic plasticity.

BACKGROUND: Thellungiella salsuginea is an important model plant due to its natural tolerance to abiotic stresses including salt, cold, and water deficits. Microarray and metabolite profiling have shown that Thellungiella undergoes stress-responsive changes in transcript and organic solute abundance when grown under controlled environmental conditions. However, few reports assess the capacity of plants to display stress-responsive traits in natural habitats where concurrent stresses are the norm. RESULTS: To determine whether stress-responsive changes observed in cabinet-grown plants are recapitulated in the field, we analyzed leaf transcript and metabolic profiles of Thellungiella growing in its native Yukon habitat during two years of contrasting meteorological conditions. We found 673 genes showing differential expression between field and unstressed, chamber-grown plants. There were comparatively few overlaps between genes expressed under field and cabinet treatment-specific conditions. Only 20 of 99 drought-responsive genes were expressed both in the field during a year of low precipitation and in plants subjected to drought treatments in cabinets. There was also a general pattern of lower abundance among metabolites found in field plants relative to control or stress-treated plants in growth cabinets. Nutrient availability may explain some of the observed differences. For example, proline accumulated to high levels in cold and salt-stressed cabinet-grown plants but proline content was, by comparison, negligible in plants at a saline Yukon field site. We show that proline accumulated in a stress-responsive manner in Thellungiella plants salinized in growth cabinets and in salt-stressed seedlings when nitrogen was provided at 1.0 mM. In seedlings grown on 0.1 mM nitrogen medium, the proline content was low while carbohydrates increased. The relatively higher content of sugar-like compounds in field plants and seedlings on low nitrogen media suggests that Thellungiella shows metabolic plasticity in response to environmental stress and that resource availability can influence the expression of stress tolerance traits under field conditions. CONCLUSION: Comparisons between Thellungiella plants responding to stress in cabinets and in their natural habitats showed differences but also overlap between transcript and metabolite profiles. The traits in common offer potential targets for improving crops that must respond appropriately to multiple, concurrent stresses.

Identification of phosphomethylethanolamine N-methyltransferase from Arabidopsis and its role in choline and phospholipid metabolism.

Three sequential methylations of phosphoethanolamine (PEA) are required for the synthesis of phosphocholine (PCho) in plants. A cDNA encoding an N-methyltransferase that catalyzes the last two methylation steps was cloned from Arabidopsis by heterologous complementation of a Saccharomyces cerevisiae cho2, opi3 mutant. The cDNA encodes phosphomethylethanolamine N-methyltransferase (PMEAMT), a polypeptide of 475 amino acids that is organized as two tandem methyltransferase domains. PMEAMT shows 87% amino acid identity to a related enzyme, phosphoethanolamine N-methyltransferase, an enzyme in plants that catalyzes all three methylations of PEA to PCho. PMEAMT cannot use PEA as a substrate, but assays using phosphomethylethanolamine as a substrate result in both phosphodimethylethanolamine and PCho as products. PMEAMT is inhibited by the reaction products PCho and S-adenosyl-l-homocysteine, a property reported for phosphoethanolamine N-methyltransferase from various plants. An Arabidopsis mutant with a T-DNA insertion associated with locus At1g48600 showed no transcripts encoding PMEAMT. Shotgun lipidomic analyses of leaves of atpmeamt and wild-type plants generated phospholipid profiles showing the content of phosphatidylmethylethanolamine to be altered relative to wild type with the content of a 34:3 lipid molecular species 2-fold higher in mutant plants. In S. cerevisiae, an increase in PtdMEA in membranes is associated with reduced viability. This raises a question regarding the role of PMEAMT in plants and whether it serves to prevent the accumulation of PtdMEA to potentially deleterious levels.

Adenosine kinase deficiency is associated with developmental abnormalities and reduced transmethylation.

Adenosine (Ado) kinase (ADK; ATP:Ado 5' phosphotransferase, EC 2.7.1.20) catalyzes the salvage synthesis of adenine monophosphate from Ado and ATP. In Arabidopsis, ADK is encoded by two cDNAs that share 89% nucleotide identity and are constitutively, yet differentially, expressed in leaves, stems, roots, and flowers. To investigate the role of ADK in plant metabolism, lines deficient in this enzyme activity have been created by sense and antisense expression of the ADK1 cDNA. The levels of ADK activity in these lines range from 7% to 70% of the activity found in wild-type Arabidopsis. Transgenic plants with 50% or more of the wild-type activity have a normal morphology. In contrast, plants with less than 10% ADK activity are small with rounded, wavy leaves and a compact, bushy appearance. Because of the lack of elongation of the primary shoot, the siliques extend in a cluster from the rosette. Fertility is decreased because the stamen filaments do not elongate normally; hypocotyl and root elongation are reduced also. The hydrolysis of S-adenosyl-L-homo-cysteine (SAH) produced from S-adenosyl-L-methionine (SAM)-dependent methylation reactions is a key source of Ado in plants. The lack of Ado salvage in the ADK-deficient lines leads to an increase in the SAH level and results in the inhibition of SAM-dependent transmethylation. There is a direct correlation between ADK activity and the level of methylesterified pectin in seed mucilage, as monitored by staining with ruthenium red, immunofluorescence labeling, or direct assay. These results indicate that Ado must be steadily removed by ADK to prevent feedback inhibition of SAH hydrolase and maintain SAM utilization and recycling.