Transcriptome Analysis Reveals Roles of Anthocyanin- and Jasmonic Acid-Biosynthetic Pathways in Rapeseed in Response to High Light Stress.

Rapeseed (Brassica napus) is one of the major important oil crops worldwide and is largely cultivated in the Qinghai-Tibetan plateau (QTP), where long and strong solar-radiation is well-known. However, the molecular mechanisms underlying rapeseed's response to light stress are largely unknown. In the present study, the color of rapeseed seedlings changed from green to purple under high light (HL) stress conditions. Therefore, changes in anthocyanin metabolism and the transcriptome of rapeseed seedlings cultured under normal light (NL) and HL conditions were analyzed to dissect how rapeseed responds to HL at the molecular level. Results indicated that the contents of anthocyanins, especially glucosides of cyanidin, delphinidin, and petunidin, which were determined by liquid chromatography-mass spectrometry (LC-MS), increased by 9.6-, 4.2-, and 59.7-fold in rapeseed seedlings exposed to HL conditions, respectively. Next, RNA-sequencing analysis identified 7390 differentially expressed genes (DEGs), which included 4393 up-regulated and 2997 down-regulated genes. Among the up-regulated genes, many genes related to the anthocyanin-biosynthetic pathway were enriched. For example, genes encoding dihydroflavonol reductase (BnDFR) and anthocyanin synthase (BnANS) were especially induced by HL conditions, which was also confirmed by RT-qPCR analysis. In addition, two PRODUCTION OF ANTHOCYANIN PIGMENTATION 2 (BnPAP2) and GLABRA3 (BnGL3) genes encoding MYB-type and bHLH-type transcription factors, respectively, whose expression was also up-regulated by HL stress, were found to be associated with the changes in anthocyanin biosynthesis. Many genes involved in the jasmonic acid (JA)-biosynthetic pathway were also up-regulated under HL conditions. This finding, which is in agreement with the well-known positive regulatory role of JA in anthocyanin biosynthesis, suggests that the JA may also play a key role in the responses of rapeseed seedlings to HL. Collectively, these data indicate that anthocyanin biosynthesis-related and JA biosynthesis-related pathways mediate HL responses in rapeseed. These findings collectively provide mechanistic insights into the mechanisms involved in the response of rapeseed to HL stress, and the identified key genes may potentially be used to improve HL tolerance of rapeseed cultivars through genetic engineering or breeding strategies.

The East Asian Winter Monsoon Acts as a Major Selective Factor in the Intraspecific Differentiation of Drought-Tolerant Nitraria tangutorum in Northwest China.

The influence of Quaternary climate fluctuation on the geographical structure and genetic diversity of species distributed in the regions of the Qinghai-Tibet Plateau (QTP) has been well established. However, the underlying role of the East Asian monsoon system (EAMS) in shaping the genetic structure of the population and the demography of plants located in the arid northwest of China has not been explored. In the present study, Nitraria tangutorum, a drought-tolerant desert shrub that is distributed in the EAMS zone and has substantial ecological and economic value, was profiled to better understand the influence of EAMS evolution on its biogeographical patterns and demographic history. Thus, the phylogeographical structure and historical dynamics of this plant species were elucidated using its five chloroplast DNA (cpDNA) fragments. Hierarchical structure analysis revealed three distinct, divergent lineages: West, East-A, and East-B. The molecular dating was carried out using a Bayesian approach to estimate the time of intraspecies divergence. Notably, the eastern region, which included East-A and East-B lineages, was revealed to be the original center of distribution and was characterized by a high level of genetic diversity, with the intraspecific divergence time dated to be around 2.53 million years ago (Ma). These findings, combined with the data obtained by ecological niche modeling analysis, indicated that the East lineages have undergone population expansion and differentiation, which were closely correlated with the development of the EAMS, especially the East Asian winter monsoon (EAWM). The West lineage appears to have originated from the migration of N. tangutorum across the Hexi corridor at around 1.85 Ma, and subsequent colonization of the western region. These results suggest that the EAWM accelerated the population expansion of N. tangutorum and subsequent intraspecific differentiation. These findings collectively provide new information on the impact of the evolution of the EAMS on intraspecific diversification and population demography of drought-tolerant plant species in northwest China.

Transcriptome Analysis Reveals Potential Roles of Abscisic Acid and Polyphenols in Adaptation of Onobrychis viciifolia to Extreme Environmental Conditions in the Qinghai-Tibetan Plateau.

A detailed understanding of the molecular mechanisms of plant stress resistance in the face of ever-changing environmental stimuli will be helpful for promoting the growth and production of crop and forage plants. Investigations of plant responses to various single abiotic or biotic factors, or combined stresses, have been extensively reported. However, the molecular mechanisms of plants in responses to environmental stresses under natural conditions are not clearly understood. In this study, we carried out a transcriptome analysis using RNA-sequencing to decipher the underlying molecular mechanisms of Onobrychis viciifolia responding and adapting to the extreme natural environment in the Qinghai-Tibetan Plateau (QTP). The transcriptome data of plant samples collected from two different altitudes revealed a total of 8212 differentially expressed genes (DEGs), including 5387 up-regulated and 2825 down-regulated genes. Detailed analysis of the identified DEGs uncovered that up-regulation of genes potentially leading to changes in hormone homeostasis and signaling, particularly abscisic acid-related ones, and enhanced biosynthesis of polyphenols play vital roles in the adaptive processes of O. viciifolia. Interestingly, several DEGs encoding uridine diphosphate glycosyltransferases, which putatively regulate phytohormone homeostasis to resist environmental stresses, were also discovered. Furthermore, numerous DEGs encoding transcriptional factors, such as members of the myeloblastosis (MYB), homeodomain-leucine zipper (HD-ZIP), WRKY, and nam-ataf1,2-cuc2 (NAC) families, might be involved in the adaptive responses of O. viciifolia to the extreme natural environmental conditions. The DEGs identified in this study represent candidate targets for improving environmental stress resistance of O. viciifolia grown in higher altitudes of the QTP, and can provide deep insights into the molecular mechanisms underlying the responses of this plant species to the extreme natural environmental conditions of the QTP.

Natural Products, Traditional Uses and Pharmacological Activities of the Genus Biebersteinia (Biebersteiniaceae).

Medicinal plants have been known as a rich source of natural products (NPs). Due to their diverse chemical structures and remarkable pharmacological activities, NPs are regarded as important repertoires for drug discovery and development. Biebersteinia plant species belong to the Biebersteiniaceae family, and have been used in folk medicines in China and Iran for ages. However, the chemical properties, bioactivities and modes of action of the NPs produced by medicinal Biebersteinia species are poorly understood despite the fact that there are only four known Biebersteinia species worldwide. Here, we reviewed the chemical classifications and diversity of the various NPs found in the four known Biebersteinia species. We found that the major chemical categories in these plants include flavonoids, alkaloids, phenylpropanoids, terpenoids, essential oils and fatty acids. We also discussed the anti-inflammatory, analgesic, antibacterial, antioxidant, antihypertensive and hypoglycemic effects of the four Biebersteinia species. We believe that the present review will facilitate the exploration of traditional uses and pharmacological properties of Biebersteinia species, extraction of the NPs and elucidation of their molecular mechanisms, as well as the development of novel drugs based on the reported properties and mode-of-action.

Adaptive signals of flowering time pathways in wild barley from Israel over 28 generations.

Flowering time is one of the most critical traits for plants' life cycles, which is influenced by various environment changes, such as global warming. Previous studies have suggested that to guarantee reproductive success, plants have shifted flowering times to adapt to global warming. Although many studies focused on the molecular mechanisms of early flowering, little was supported by the repeated sampling at different time points through the changing climate. To fully dissect the temporal and spatial evolutionary genetics of flowering time, we investigated nucleotide variation in ten flowering time candidate genes and nine reference genes for the same ten wild-barley populations sampled 28 years apart (1980-2008). The overall genetic differentiation was significantly greater in the descendant populations (2008) compared with the ancestral populations (1980); however, local adaptation tests failed to detect any single-nucleotide polymorphism (SNP)/indel under spatial-diversifying selection at either time point. By contrast, the WFABC (Wright-Fisher ABC-based approach) that detected 54 SNPs/indels was under strong selection during the past 28 generations. Moreover, all these 54 alleles were segregated in the ancestral populations, but fixed in the descendent populations. Among the top ten SNPs/indels, seven were located in genes of FT1 (FLOWERING TIME LOCUS T 1), CO1 (CONSTANS-LIKE PROTEIN 1), and VRN-H2 (VERNALIZATION-H2), which have been documented to be associated with flowering time regulation in barley cultivars. This study might suggest that all ten populations have undergone parallel evolution over the past few decades in response to global warming, and even an overwhelming local adaptation and ecological differentiation.

Transcriptomes Divergence of Ricotia lunaria Between the Two Micro-Climatic Divergent Slopes at "Evolution Canyon" I, Israel.

As one of the hotspot regions for sympatric speciation studies, Evolution Canyon (EC) became an ideal place for its high level of microclimatic divergence interslopes. In this study, to highlight the genetic mechanisms of sympatric speciation, phenotypic variation on flowering time and transcriptomic divergence were investigated between two ecotypes of Ricotia lunaria, which inhabit the opposite temperate and tropical slopes of EC I (Lower Nahal Oren, Mount Carmel, Israel) separated by 100 m at the bottom of the slopes. Growth chamber results showed that flowering time of the ecotype from south-facing slope population # 3 (SFS 3) was significantly 3 months ahead of the north-facing slope population # 5 (NFS 5). At the same floral development stage, transcriptome analysis showed that 1,064 unigenes were differentially expressed between the two ecotypes, which enriched in the four main pathways involved in abiotic and/or biotic stresses responses, including flavonoid biosynthesis, α-linolenic acid metabolism, plant-pathogen interaction and linoleic acid metabolism. Furthermore, based on Ka/Ks analysis, nine genes were suggested to be involved in the ecological divergence between the two ecotypes, whose homologs functioned in RNA editing, ABA signaling, photoprotective response, chloroplasts protein-conducting channel, and carbohydrate metabolism in Arabidopsis thaliana. Among them, four genes, namely, SPDS1, FCLY, Tic21 and BGLU25, also showed adaptive divergence between R. lunaria and A. thaliana, suggesting that these genes could play an important role in plant speciation, at least in Brassicaceae. Based on results of both the phenotype of flowering time and comparative transcriptome, we hypothesize that, after long-time local adaptations to their interslope microclimatic environments, the molecular functions of these nine genes could have been diverged between the two ecotypes. They might differentially regulate the expression of the downstream genes and pathways that are involved in the interslope abiotic stresses, which could further diverge the flowering time between the two ecotypes, and finally induce the reproductive isolation establishment by natural selection overruling interslope gene flow, promoting sympatric speciation.

Transcriptional responses to phosphate starvation in Brachypodium distachyon roots.

Brachypodium distachyon is a model plant that has recently emerged in grass research. Although the growth and photochemical efficiency of this species respond strongly to phosphate (Pi) availability, its Pi starvation response mechanism, which controls the Pi homeostasis, remains largely unknown. This study presents the transcriptomic response profiles of Pi-deficient roots at growth stages during which the plants were starved but obvious growth defects were absent. The results identify several phosphate transporters (i.e., PHO1), purple acid phosphatases, and SYG1/PHO81/XPR1 (SPX) domain-containing proteins out of a total of 1740 differentially expressed genes (DEGs). In particular, the transcription factor ethylene response factor (ERF), basic helix-loop-helix (bHLH), and WRKY family genes were the three most abundant DEG groups and the latter was significantly enriched. Comparative transcriptome analysis of Brachypodium versus Arabidopsis and rice revealed the presence of several common components in response to Pi fluctuations. Most significantly, jasmonic acid (JA) signaling-related genes were overrepresented in gene ontology (GO) enrichment tests. The presence of a possible link between low Pi response, inositol polyphosphates, and JA signaling is therefore discussed.

Population dynamics of Agriophyllum squarrosum, a pioneer annual plant endemic to mobile sand dunes, in response to global climate change.

Climate change plays an important role in the transition of ecosystems. Stratigraphic investigations have suggested that the Asian interior experienced frequent transitions between grassland and desert ecosystems as a consequence of global climate change. Using maternally and bi-parentally inherited markers, we investigated the population dynamics of Agriophyllum squarrosum (Chenopodiaceae), an annual pioneer plant endemic to mobile sand dunes. Phylogeographic analysis revealed that A. squarrosum could originate from Gurbantunggut desert since ~1.6 Ma, and subsequently underwent three waves of colonisation into other deserts and sandy lands corresponding to several glaciations. The rapid population expansion and distribution range shifts of A. squarrosum from monsoonal climate zones suggested that the development of the monsoonal climate significantly enhanced the population growth and gene flow of A. squarrosum. These data also suggested that desertification of the fragile grassland ecosystems in the Qinghai-Tibetan Plateau was more ancient than previously suggested and will be aggravated under global warming in the future. This study provides new molecular phylogeographic insights into how pioneer annual plant species in desert ecosystems respond to global climate change, and facilitates evaluation of the ecological potential and genetic resources of future crops for non-arable dry lands to mitigate climate change.

The role of East Asian monsoon system in shaping population divergence and dynamics of a constructive desert shrub Reaumuria soongarica.

Both of the uplift of Qinghai-Tibet Plateau (QTP) and the development of East Asian monsoon system (EAMS) could have comprehensively impacted the formation and evolution of Arid Central Asia (ACA). To understand how desert plants endemic to ACA responded to these two factors, we profiled the historical population dynamics and distribution range shift of a constructive desert shrub Reaumuria soongarica (Tamaricaceae) based on species wide investigation of sequence variation of chloroplast DNA and nuclear ribosomal ITS. Phylogenetic analysis uncovered a deep divergence occurring at ca. 2.96 Mya between the western and eastern lineages of R. soongarica, and ecological niche modeling analysis strongly supported that the monsoonal climate could have fragmented its habitats in both glacial and interglacial periods and impelled its intraspecific divergence. Additionally, the population from the east monsoonal zone expanded rapidly, suggesting that the local monsoonal climate significantly impacted its population dynamics. The isolation by distance tests supported strong maternal gene flow along the direction of the East Asian winter monsoon, whose intensification induced the genetic admixture along the latitudinal populations of R. soongarica. Our results presented a new case that the development of EAMS had prominently impacted the intraspecific divergence and population dynamics of this desert plant.

Reference gene selection for quantitative real-time PCR normalization in Reaumuria soongorica.

Despite its superiority for evaluating gene expression, real-time quantitative polymerase chain reaction (qPCR) results can be significantly biased by the use of inappropriate reference genes under different experimental conditions. Reaumuria soongorica is a dominant species of desert ecosystems in arid central Asia. Given the increasing interest in ecological engineering and potential genetic resources for arid agronomy, it is important to analyze gene function. However, systematic evaluation of stable reference genes should be performed prior to such analyses. In this study, the stabilities of 10 candidate reference genes were analyzed under 4 kinds of abiotic stresses (drought, salt, dark, and heat) within 4 accessions (HG010, HG020, XGG030, and XGG040) from 2 different habitats using 3 algorithms (geNorm, NormFinder, and BestKeeper). After validation of the ribulose-1,5-bisphosphate carboxylase/oxygenase large unite (rbcL) expression pattern, our data suggested that histone H2A (H2A) and eukaryotic initiation factor 4A-2 (EIF4A2) were the most stable reference genes, cyclophilin (CYCL) was moderate, and elongation factor 1α (EF1α) was the worst choice. This first systematic analysis for stably expressed genes will facilitate future functional analyses and deep mining of genetic resources in R. soongorica and other species of the Reaumuria genus.

Transcriptomic analysis of a tertiary relict plant, extreme xerophyte Reaumuria soongorica to identify genes related to drought adaptation.

BACKGROUND: Reaumuria soongorica is an extreme xerophyte shrub widely distributed in the desert regions including sand dune, Gobi and marginal loess of central Asia which plays a crucial role to sustain and restore fragile desert ecosystems. However, due to the lacking of the genomic sequences, studies on R. soongorica had mainly limited in physiological responses to drought stress. Here, a deep transcriptomic sequencing of R. soongorica will facilitate molecular functional studies and pave the path to understand drought adaptation for a desert plant. METHODOLOGY/PRINCIPAL FINDINGS: A total of 53,193,660 clean paired-end reads was generated from the Illumina HiSeq™ 2000 platform. By assembly with Trinity, we got 173,700 contigs and 77,647 unigenes with mean length of 677 bp and N50 of 1109 bp. Over 55% (43,054) unigenes were successfully annotated based on sequence similarity against public databases as well as Rfam and Pfam database. Local BLAST and Kyoto Encyclopedia of Genes and Genomes (KEGG) maps were used to further exhausting seek for candidate genes related to drought adaptation and a set of 123 putative candidate genes were identified. Moreover, all the C4 photosynthesis genes existed and were active in R. soongorica, which has been regarded as a typical C3 plant. CONCLUSION/SIGNIFICANCE: The assembled unigenes in present work provide abundant genomic information for the functional assignments in an extreme xerophyte R. soongorica, and will help us exploit the genetic basis of how desert plants adapt to drought environment in the near future.

Differential salt tolerance in seedlings derived from dimorphic seeds of Atriplex centralasiatica: from physiology to molecular analysis.

Seed dimorphism provides plants with alternative strategies for survival in unfavorable environments. Here, we investigated the physiological responses and differential gene expression caused by salinity exposure in Atriplex centralasiatica plants grown from the two different seed morphs. Seedlings derived from yellow seeds (YS) showed a greater salt tolerance than those derived from brown seeds (BS). Salt treatment induced nitric oxide (NO) synthesis in roots, and seedlings derived from YS produced greater amounts of NO than did those from BS. Analyses of NO scavenging during salt stress revealed that NO contributed to the differential salt tolerance in seedlings derived from the two seed morphs by modulating antioxidative enzyme activity, hydrogen peroxide accumulation and the ion equilibrium. We also applied transcriptomics and subsequent microarray analysis to evaluate the differential gene expression during salt treatment. These genes encoded proteins related to osmotic and ionic homeostasis, redox equilibrium and signal transduction. A select group of genes including GH3.3, CAT1/2, TIP1, SIHP1 and EXP1 were further confirmed with RT-PCR analysis. These results revealed that the enhanced salt tolerance of seedlings from YS appeared to be governed by a superior ability to achieve ionic homeostasis and redox equilibrium, a rapid response to salt stress, and ultimately better growth potential. NO serves as a vital regulator in these processes.

Nitric oxide is associated with long-term zinc tolerance in Solanum nigrum.

Nitric oxide (NO) has been identified as a signal molecule that interplays with reactive oxygen species in response to heavy metal stresses. Roles of NO in regulating cadmium toxicity and iron deficiency have been proposed; however, the function of NO in zinc (Zn) tolerance in plants remains unclear. Here, we investigated NO accumulation and its role in plant Zn tolerance. Zn-induced NO production promoted an increase in reactive oxygen species accumulation in Solanum nigrum roots by modulating the expression and activity of antioxidative enzymes. Subsequently, programmed cell death (PCD) was observed in primary root tips. Inhibiting NO accumulation by 2-phenyl-4,4,5,5-tetramethyl-imidazoline-1-oxyl-3-oxide (a specific NO scavenger) or N(G)-nitro-l-arginine-methyl ester (a NO synthase inhibitor) prevented the increase of superoxide radical and hydrogen peroxide as well as the subsequent cell death in the root tips, supporting the role of NO in Zn-induced PCD in the root tips. Zn-induced NO production affected the length of primary roots, the number of lateral roots, and root hair growth and thereby modulated root system architecture and activity. Investigation of metal contents in Zn-treated roots suggests that NO is required for metal (especially iron) uptake and homeostasis in plants exposed to excess Zn. Taken together, our results indicate that NO production and the subsequent PCD in root tips exposed to excess Zn are favorable for the S. nigrum seedling response to long-term Zn toxicity by modulating root system architecture and subsequent adaptation to Zn stress.

Salt affects plant Cd-stress responses by modulating growth and Cd accumulation.

Cadmium contamination is a serious environmental problem for modern agriculture and human health. Salinity affects plant growth and development, and interactions between salt and cadmium have been reported. However, the molecular mechanisms of salinity-cadmium interactions are not fully understood. Here, we show that a low concentration of salt alleviates Cd-induced growth inhibition and increases Cd accumulation in Arabidopsis thaliana. Supplementation with low concentrations of salt reduced the reactive oxygen species level in Cd-stressed roots by increasing the contents of proline and glutathione and down-regulating the expression of RCD1, thereby protecting the plasma membrane integrity of roots under cadmium stress. Salt supplementation substantially reduces the Cd-induced elevation of IAA oxidase activity, thereby maintaining auxin levels in Cd-stressed plants, as indicated by DR5::GUS expression. Salt supply increased Cd absorption in roots and increased Cd accumulation in leaves, implying that salt enhances both Cd uptake in roots and the root-to-shoot translocation of Cd. The elevated Cd accumulation in plants in response to salt was found to be correlated with the elevated levels of phytochelatin the expression of heavy metal transporters AtHMA1-4, especially AtHMA4. Salt alleviated growth inhibition caused by Cd and increased Cd accumulation also was observed in Cd accumulator Solanum nigrum.

Protective effects of proline against cadmium toxicity in micropropagated hyperaccumulator, Solanum nigrum L.

Solanum nigrum is a newly discovered Cd-hyperaccumulator. In the present study, the protective effects of proline against cadmium toxicity of callus and regenerated shoots of S. nigrum are investigated based on a high frequency in vitro shoot regeneration system. Proline pretreatment reduces the reactive oxygen species levels and protects the plasma membrane integrity of callus under cadmium stress, and therefore improves the cadmium tolerance in S. nigrum. Inductively coupled plasma mass spectroscopy analysis shows that exogenous proline increases the cadmium accumulation in callus and regenerated shoots of S. nigrum. Further analysis indicates that the improvement of cadmium tolerance caused by proline pretreatment is correlated with an increase of superoxide dismutase and catalase activity and intracellular total glutathione content. The interaction between proline and enzymic or non-enzymic antioxidants is discussed.