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.

Molecular Traits of Long Non-protein Coding RNAs from Diverse Plant Species Show Little Evidence of Phylogenetic Relationships.

Long non-coding RNAs (lncRNAs) represent a diverse class of regulatory loci with roles in development and stress responses throughout all kingdoms of life. LncRNAs, however, remain under-studied in plants compared to animal systems. To address this deficiency, we applied a machine learning prediction tool, Classifying RNA by Ensemble Machine learning Algorithm (CREMA), to analyze RNAseq data from 11 plant species chosen to represent a wide range of evolutionary histories. Transcript sequences of all expressed and/or annotated loci from plants grown in unstressed (control) conditions were assembled and input into CREMA for comparative analyses. On average, 6.4% of the plant transcripts were identified by CREMA as encoding lncRNAs. Gene annotation associated with the transcripts showed that up to 99% of all predicted lncRNAs for Solanum tuberosum and Amborella trichopoda were missing from their reference annotations whereas the reference annotation for the genetic model plant Arabidopsis thaliana contains 96% of all predicted lncRNAs for this species. Thus a reliance on reference annotations for use in lncRNA research in less well-studied plants can be impeded by the near absence of annotations associated with these regulatory transcripts. Moreover, our work using phylogenetic signal analyses suggests that molecular traits of plant lncRNAs display different evolutionary patterns than all other transcripts in plants and have molecular traits that do not follow a classic evolutionary pattern. Specifically, GC content was the only tested trait of lncRNAs with consistently significant and high phylogenetic signal, contrary to high signal in all tested molecular traits for the other transcripts in our tested plant species.

Demographic expansion and genetic load of the halophyte model plant Eutrema salsugineum.

The halophyte model plant Eutrema salsugineum (Brassicaceae) disjunctly occurs in temperate to subarctic Asia and North America. This vast, yet extremely discontinuous distribution constitutes an ideal system to examine long-distance dispersal and the ensuing accumulation of deleterious mutations as expected in expanding populations of selfing plants. In this study, we resequenced individuals from 23 populations across the range of E. salsugineum. Our population genomic data indicate that E. salsugineum migrated "out of the Altai region" at least three times to colonize northern China, northeast Russia and western China. It then expanded its distribution into North America independently from northeast Russia and northern China, respectively. The species colonized northern China around 33.7 thousand years ago (kya) and underwent a considerable expansion in range size approximately 7-8 kya. The western China lineage is likely a hybrid derivative of the northern China and Altai lineages, originating approximately 25-30 kya. Deleterious alleles accumulated in a stepwise manner from (a) Altai to northern China and North America and (b) Altai to northeast Russia and North America. In summary, E. salsugineum dispersed from Asia to North America and deleterious mutations accumulated in a stepwise manner during the expansion of the species' distribution.

Prediction of plant lncRNA by ensemble machine learning classifiers.

BACKGROUND: In plants, long non-protein coding RNAs are believed to have essential roles in development and stress responses. However, relative to advances on discerning biological roles for long non-protein coding RNAs in animal systems, this RNA class in plants is largely understudied. With comparatively few validated plant long non-coding RNAs, research on this potentially critical class of RNA is hindered by a lack of appropriate prediction tools and databases. Supervised learning models trained on data sets of mostly non-validated, non-coding transcripts have been previously used to identify this enigmatic RNA class with applications largely focused on animal systems. Our approach uses a training set comprised only of empirically validated long non-protein coding RNAs from plant, animal, and viral sources to predict and rank candidate long non-protein coding gene products for future functional validation. RESULTS: Individual stochastic gradient boosting and random forest classifiers trained on only empirically validated long non-protein coding RNAs were constructed. In order to use the strengths of multiple classifiers, we combined multiple models into a single stacking meta-learner. This ensemble approach benefits from the diversity of several learners to effectively identify putative plant long non-coding RNAs from transcript sequence features. When the predicted genes identified by the ensemble classifier were compared to those listed in GreeNC, an established plant long non-coding RNA database, overlap for predicted genes from Arabidopsis thaliana, Oryza sativa and Eutrema salsugineum ranged from 51 to 83% with the highest agreement in Eutrema salsugineum. Most of the highest ranking predictions from Arabidopsis thaliana were annotated as potential natural antisense genes, pseudogenes, transposable elements, or simply computationally predicted hypothetical protein. Due to the nature of this tool, the model can be updated as new long non-protein coding transcripts are identified and functionally verified. CONCLUSIONS: This ensemble classifier is an accurate tool that can be used to rank long non-protein coding RNA predictions for use in conjunction with gene expression studies. Selection of plant transcripts with a high potential for regulatory roles as long non-protein coding RNAs will advance research in the elucidation of long non-protein coding RNA function.

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.

Exposure of two Eutrema salsugineum (Thellungiella salsuginea) accessions to water deficits reveals different coping strategies in response to drought.

Eutrema salsugineum is an extremophile related to Arabidopsis. Accessions from Yukon, Canada and Shandong, China, were evaluated for their tolerance to water deficits. Plants were exposed to two periods of water deficit separated by an interval of re-watering and recovery. All plants took the same time to wilt during the first drought exposure but Yukon plants took 1 day longer than Shandong plants following the second drought treatment. Following re-watering and turgor recovery, solute potentials of Shandong leaves returned to predrought values while those of Yukon leaves were lower than predrought levels consistent with having undergone osmotic adjustment. Polar metabolites profiled in re-watered plants showed that different metabolites are accumulated by Yukon and Shandong plants recovering from a water deficit with glucose more abundant in Yukon and fructose in Shandong leaves. The drought-responsive expression of dehydrin genes RAB18, ERD1, RD29A and RD22 showed greater changes in transcript abundance in Yukon relative to Shandong leaves during both water deficits and recovery with the greatest difference in expression appearing during the second drought. We propose that the initial exposure of Yukon plants to drought renders them more resilient to water loss during a subsequent water deficit leading to delayed wilting. Yukon plants also established a high leaf water content and increased specific leaf area during the second deficit. Shandong plants undergoing the same treatment regime do not show the same beneficial drought tolerance responses and likely use drought avoidance to cope with water deficits.

Investigation of intercellular salicylic acid accumulation during compatible and incompatible Arabidopsis-pseudomonas syringae interactions using a fast neutron-generated mutant allele of EDS5 identified by genetic mapping and whole-genome sequencing.

A whole-genome sequencing technique developed to identify fast neutron-induced deletion mutations revealed that iap1-1 is a new allele of EDS5 (eds5-5). RPS2-AvrRpt2-initiated effector-triggered immunity (ETI) was compromised in iap1-1/eds5-5 with respect to in planta bacterial levels and the hypersensitive response, while intra- and intercellular free salicylic acid (SA) accumulation was greatly reduced, suggesting that SA contributes as both an intracellular signaling molecule and an antimicrobial agent in the intercellular space during ETI. During the compatible interaction between wild-type Col-0 and virulent Pseudomonas syringae pv. tomato (Pst), little intercellular free SA accumulated, which led to the hypothesis that Pst suppresses intercellular SA accumulation. When Col-0 was inoculated with a coronatine-deficient strain of Pst, high levels of intercellular SA accumulation were observed, suggesting that Pst suppresses intercellular SA accumulation using its phytotoxin coronatine. This work suggests that accumulation of SA in the intercellular space is an important component of basal/PAMP-triggered immunity as well as ETI to pathogens that colonize the intercellular space.

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.

Thellungiella: an Arabidopsis-related model plant adapted to cold temperatures.

Thellungiella salsuginea, a wild crucifer that grows in subarctic Canada and is closely related to Arabidopsis thaliana, was examined for its suitability as a model plant for studies of tolerance to cold and freezing temperatures. Thellungiella completed its life cycle at 5 degrees C, demonstrating that temperature-sensitive processes such as seed germination and the production of pollen and seeds were resistant to cold temperatures. Moreover, the plant exhibited dramatically different vegetative and flowering phenotypes in response to growth at cold temperature and shifts to cold temperature. Northern analyses showed that genes induced by cold in Arabidopsis, including CBF1, the transcriptional activator for the cold-regulated (COR) genes COR15a and COR47, were also expressed in Thellungiella. Freezing tolerance, assayed by the regrowth of intact plants, increased from -13.0 to -18.5 degrees C after cold treatment. The plants lacked endogenous ice nucleation or anti-freeze activity, indicating a potential for supercooling. As a close relative to Arabidopsis, Thellungiella exhibits extreme cold tolerance and should be an important model system in the elucidation of stress tolerance mechanisms.

Transcriptional profiling implicates novel interactions between abiotic stress and hormonal responses in Thellungiella, a close relative of Arabidopsis.

Thellungiella, an Arabidopsis (Arabidopsis thaliana)-related halophyte, is an emerging model species for studies designed to elucidate molecular mechanisms of abiotic stress tolerance. Using a cDNA microarray containing 3,628 unique sequences derived from previously described libraries of stress-induced cDNAs of the Yukon ecotype of Thellungiella salsuginea, we obtained transcript profiles of its response to cold, salinity, simulated drought, and rewatering after simulated drought. A total of 154 transcripts were differentially regulated under the conditions studied. Only six of these genes responded to all three stresses of drought, cold, and salinity, indicating a divergence among the end responses triggered by each of these stresses. Unlike in Arabidopsis, there were relatively few transcript changes in response to high salinity in this halophyte. Furthermore, the gene products represented among drought-responsive transcripts in Thellungiella associate a down-regulation of defense-related transcripts with exposure to water deficits. This antagonistic interaction between drought and biotic stress response may demonstrate Thellungiella's ability to respond precisely to environmental stresses, thereby conserving energy and resources and maximizing its survival potential. Intriguingly, changes of transcript abundance in response to cold implicate the involvement of jasmonic acid. While transcripts associated with photosynthetic processes were repressed by cold, physiological responses in plants developed at low temperature suggest a novel mechanism for photosynthetic acclimation. Taken together, our results provide useful starting points for more in-depth analyses of Thellungiella's extreme stress tolerance.

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.