Transcriptomic analysis of salt stress responsive genes in Rhazya stricta

Comparative analysis of the complete mitogenome of Geoffroea decorticans: a native tree surviving in the Atacama Desert

Chañar (Geoffroea decorticans (Gill., ex Hook. & Arn.) Burkart) has been highly significant for indigenous people in the Atacama Desert for over 3,000 years. Through evolutionary processes, the G. decorticans mitogenome likely underwent changes facilitating its adaptation to the extreme conditions of the Atacama Desert. Here, we compare the mitochondrial genome of G. decorticans with those of other Papilionoideae family species. The complete mitogenome of G. decorticans was sequenced and assembled, making it the first in the genus Geoffroea. The mitogenome contained 383,963 base pairs, consisting of 33 protein coding genes, 21 transfer RNA genes, and 3 ribosomal RNA genes. The Chañar mitogenome is relatively compact, and has two intact genes (sdh4 and nad1) which were not observed in most other species. Additionally, Chañar possessed the highest amount of mitochondrial DNA of plastid origin among angiosperm species. The phylogenetic analysis of the mitogenomes of Chañar and 12 other taxa displayed a high level of consistency in taxonomic classification, when compared to those of the plastid genome. Atp8 was subjected to positive selection, while the ccmFc and rps1 were subjected to neutral selection. This study provides valuable information regarding its ability to survive the extreme environmental conditions of the Atacama Desert.

Deciphering the enigma of RNA editing in the ATP1_alpha subunit of ATP synthase in Triticum aestivum

There is evidence that RNA editing is related to plant cellular stress as well as electron transport organelles, such as mitochondria. The mitochondrial atp1 gene encodes the alpha-subunit of Atp synthase. Control as well as two periods of drought stress treatments were analyzed in the cDNAs generated from the mitochondrial atp1 gene of two cultivars of Triticum aestivum [Giza 168 (G168) and Gemmiza 10 (GM10)]. Following RNA-seq data assembly, atp1 cDNAs from the control (acc. no. OQ129415), 2-hour (acc. no. OQ129416), and 12-hour (acc. no. OQ129417) time points of the T. aestivum cultivar G168 were obtained. Control (acc. no. OQ129419), 2-hour (acc. no. OQ129420), and 12-hour (acc. no. OQ129421) samples all included reconstructed atp1 transcripts from Gemmiza 10. Atp1 transcripts were assembled using the wheat atp1 gene (acc. no. NC_036024). RNA-seq raw data was utilized to identify 11 RNA editing sites in atp1 in the tolerant cultivar Giza168 and 6 in the sensitive cultivar Gemmiza10. The significant difference in RNA editing observed between control and drought stress conditions in sites led to synonymous amino acids. This led to no change in tertiary structure between tolerant and sensitive cultivars. But the change was focused between produced protein and its correspondence sequence on DNA.

Effect of salinity on ccmfn gene RNA editing of mitochondria in wild barley and uncommon types of RNA editing

The primary function of mitochondria is cellular respiration and energy production. Cytochrome C complex is an essential complex that transports electrons in the respiratory chain between complex III and complex IV. One of this complex's main subunits is CcmFN, which is believed to be crucial for holocytochrome assembly. In wild-type plant Hordeum vulgare subsp. spontaneum, four ccmfn cDNAs are subjected to high salt stress (500 mM salinity), 0 h (or control) (GenBank accession no. ON764850), after 2 h (GenBank accession no. ON7648515), after 12 h (GenBank accession no. ON764852), and after 24 h (GenBank accession no. ON764853) and mtDNA of ccmfn gene (GenBank accession no. ON764854). Using raw data from RNA-seq, 47 sites with nucleotide and amino acid modifications were detected. There were ten different RNA editing types, with most of them are C to U. Unusual editing types in plants have also been found, such as A to C, C to A, A to G, A to U, T to A, T to C, C to G, G to C, and T to G. High levels of editing were observed in control as well as treatments of salinity stress. Amino acid changes were found in 43 sites; nearly all showed hydrophilic to hydrophilic alterations. Only C749 showed regulation under salinity stress.

Distribution characteristics of selenium, cadmium and arsenic in rice grains and their genetic dissection by genome-wide association study

High selenium (Se) and low cadmium (Cd) and arsenic (As) contents in rice grains were good for human health. The genetic basis and relationship of Se, Cd and As concentrations in rice grains are still largely unknown. In the present study, large variations were observed in Se, Cd and As concentrations in brown and milled rice in normal and Se treatment conditions in 307 rice accessions from 3K Rice Genomes Project. Se fertilizer treatment greatly increased Se concentrations but had no obvious changes in concentrations of Cd and As both in brown and milled rice. Total of 237 QTL were identified for Se, Cd and As concentrations in brown and milled rice in normal and Se treatment conditions as well as ratio of concentrations under Se treatment to normal conditions. Only 19 QTL (13.4%) were mapped for concentrations of Se and Cd, Se and As, and Se, Cd and As in the same or adjacent regions, indicating that most Se concentration QTL are independent of Cd and As concentration QTL. Forty-three favorable alleles were identified for 40 candidate genes by gene-based association study and haplotype analysis in 14 important QTL regions. Se-enriched rice variety will be developed by pyramiding favorable alleles at different Se QTL and excluding undesirable alleles at Cd and As QTL, or combining favorable alleles at Se QTL with the alleles at Se-sensitive QTL by marker-assisted selection.

Identification and expression analysis of the grape pentatricopeptide repeat (PPR) gene family in abiotic stress

The pentatricopeptide repeat (PPR) is one of the largest gene family in plants, and play important role in regulating plant growth, development and abiotic stress response. However, PPR genes have been poorly studied in grapes. In this study, based on the grape genome database, bioinformatics methods and quantitative real-time PCR (qRT-PCR) were used to identify the VvPPR family and the response to abiotic stress. A total of 181 PPR genes were identified in grape and divided into two subfamilies. Subcellular localization predicted that this gene family mainly functions in chloroplasts, nucleus, and mitochondria. Protein-protein interaction prediction indicated that there may be interaction between VvPPR44,53 and VvPPR44. The promoter region of VvPPR gene family contained various cis-acting elements, which were related to light and hormone. Expression pattern analysis showed that the VvPPR gene family was highly expressed in grape leaves, buds and carpel tissues. qRT-PCR results showed that the expression of VvPPR genes in roots was higher than stems and leaves under NAA, SA, ABA, MeJA and GA3 treatments. VvPPR8 was significantly up-regulated after GA3 and MeJA treatment for 24 h, VvPPR53 was significantly up-regulated after SA, NAA, ABA and MeJA treatment. In addition, In grape leaves, VvPPR53 was up-regulated under PEG, Nacl and 4 ℃ treatments. These data indicate that VvPPR gene family members are responsive to hormones and abiotic stresses, and that there are some differences in the degree of response and expression in different grape tissues. This study provides a certain theoretical basis for grape resistance breeding.

Decoding RNA Editing Sites Through Transcriptome Analysis in Rice Under Alkaline Stress

Ribonucleic acid editing (RE) is a post-transcriptional process that altered the genetics of RNA which provide the extra level of gene expression through insertion, deletions, and substitutions. In animals, it converts nucleotide residues C-U. Similarly in plants, the role of RNA editing sites (RES) in rice under alkaline stress is not fully studied. Rice is a staple food for most of the world population. Alkaline stress cause reduction in yield. Here, we explored the effect of alkaline stress on RES in the whole mRNA from rice chloroplast and mitochondria. Ribonucleic acid editing sites in both genomes (3336 RESs) including chloroplast (345 RESs) and mitochondria (2991 RESs) with average RES efficiency ∼55% were predicted. Our findings showed that majority of editing events found in non-synonymous codon changes and change trend in amino acids was hydrophobic. Four types of RNA editing A-G (A-I), C-T (C-U), G-A, and T-C were identified in treated and untreated samples. Overall, RNA editing efficiency was increased in the treated samples. Analysis of Gene Ontology revealed that mapped genes were engaged in many biological functions and molecular processes. We also checked the expression of pentatricopeptide repeat (PPR), organelle zinc-finger (OZI), and multiple organellar RNA editing factors/RNA editing factor interacting proteins genes in control and treatment, results revealed upregulation of PPR and OZ1 genes in treated samples. This induction showed the role of these genes in RNA editing. The current findings report that RNA editing increased under alkaline stress which may contribute in adaptation for rice by changing amino acids in edited genes (88 genes). These findings will provide basis for identification of RES in other crops and also will be useful in alkaline tolerance development in rice.

Evaluating Variation in Germination and Growth of Landraces of Barley (Hordeum vulgare L.) Under Salinity Stress

Ongoing climate change is resulting in increasing areas of salinity affected soils, rising saline groundwater and droughts resulting in irrigation with brackish water. This leads to increased salinity stress in crops that are already grown on marginal agricultural lands, such as barley. Tolerance to salinity stress is limited in the elite barley cultivar pools, but landraces of barley hold potential sources of tolerance due to their continuous selection on marginal lands. This study analyzed 140 heritage cultivars and landrace lines of barley, including 37 Scottish Bere lines that were selected from coastal regions, to screen for tolerance to salinity stress. Tolerance to salinity stress was screened by looking at the germination speed and the early root growth during germination, and the pre-maturity biomass accumulation during early growth stages. Results showed that most lines increased germination time, and decreased shoot biomass and early root growth with greater salinity stress. Elite cultivars showed increased response to the salinity, compared to the landrace lines. Individual Bere and landrace lines showed little to no effect of increased salinity in one or more experiments, one line showed high salinity tolerance in all experiments-Bere 49 A 27 Shetland. A Genome Wide Association Screening identified a number of genomic regions associated with increased tolerance to salinity stress. Two chromosomal regions were found, one associated with shoot biomass on 5HL, and another associated with early root growth, in each of the salinities, on 3HS. Within these regions a number of promising candidate genes were identified. Further analysis of these new regions and candidate genes should be undertaken, along with field trials, to identify targets for future breeding for salinity tolerance.

The first report of RNA U to C or G editing in the mitochondrial NADH dehydrogenase subunit 5 (Nad5) transcript of wild barley

BACKGROUND

Nad dehydrogenase complex in mtDNA has a significant role in cellular respiration. One of the largest subunits in the complex is subunit 5 (Nad5).

METHODS AND RESULTS

Four cDNAs of the Hordeum vulgare subsp. spontaneum nad5 gene have been characterized and subjected to four phases of 0.5 M salinity, at 0 h (control, accession no. MT235236), after 2 h (acc. no. MT235237), after 12 h (acc. no. MT235238) and after 24 h (acc. no. MT235239). Utilizing raw data from RNA-seq, ten RNA editing sites were reported. Seven sites have common editing from C to U in positions (C1490, C1859, C1895, C1900, C1901, C1916, C1918). A rare editing event U to C was detected in two positions (U1650 and U1652) and a novel editing event U to G was for the first time in positions nad5-U231. The highest editing level was shown in 2 and 12 h after salinity exposure. After 24 h, these edits were disrupted, possibly due to the launch of the programed cell death mechanism. However, the RNA editing in positions U1650, U1652 and U231 was fixed at all exposure times.

CONCLUSIONS

Although study clarified the role of salinity stress in nad5 RNA editing sites, the main achievements are first report of U to G RNA editing in plants at position U231 and first report of U to C editing in the nad5 gene at U1650 and U1652.

Chlorophyll Metabolism and Gene Expression in Response to Submergence Stress and Subsequent Recovery in Perennial Ryegrass Accessions Differing in Growth Habits

Submergence-induced leaf senescence may alter chlorophyll metabolism. The objective of this study was to characterize chlorophyll biosynthesis and degradation in contrasting perennial ryegrass (Lolium perenne) in response to submergence stress and recovery. The light-green and fast-growing accession PI238938 and the darker-green and slow-growing cultivar BrightStar SLT were exposed to 0, 6 h, 1-, 3-, and 7-d of submergence stress and 1- and 5- d of de-submergence, respectively. Plant growth of PI238938 were more severely inhibited by submergence stress and recovery. Both accessions showed increased leaf malondialdehyde under stress and recovery, but reduced chlorophyll (Chl) concentrations were observed at 3- and 7-d of stress and at recovery. The reduction in Chl was more severe in BrightStar SLT at 7 d of stress. The concentration of 5-aminolevulenic acid was unaffected by stress but increased at 1d of recovery. Activities of 5-aminolevulinic acid dehydratase (ALAD) involved in Chl biosynthesis remained unchanged under stress and recovery, while the activities of Chl degrading enzymes chlorophyllase (CHL) and pheophytinase (PPH) increased at 3 d or 7 d of stress, and returned to the control level after recovery in both accessions. The downregulation of Chl-biosynthetic genes CHLI, POR, and CHLP and the upregulation of Chl-degrading genes CLH, PPH, and SGR were observed in both accessions under most of the stress periods. BrightStar SLT exhibited much lower expressions of the Chl-biosynthetic genes PBGD, CHS, and CHID under stress, while PI238938 had remarkably higher expressions of genes involved in Chl breakdown including CLH, PPH, PAO, RCCR, and SGR, and the expressions of these genes remained at a higher level at 1 d of recovery. The results indicated that submergence-induced leaf senescence and declines in Chl were associated with downregulation of more Chl-biosynthetic genes in slow-growing genotype and upregulation of more Chl-degrading genes in fast-growing genotype of perennial ryegrass.

Protective Roles of Cytosolic and Plastidal Proteasomes on Abiotic Stress and Pathogen Invasion

Protein malfunction is typically caused by abiotic stressors. To ensure cell survival during conditions of stress, it is important for plant cells to maintain proteins in their respective functional conformation. Self-compartmentalizing proteases, such as ATP-dependent Clp proteases and proteasomes are designed to act in the crowded cellular environment, and they are responsible for degradation of misfolded or damaged proteins within the cell. During different types of stress conditions, the levels of misfolded or orphaned proteins that are degraded by the 26S proteasome in the cytosol and nucleus and by the Clp proteases in the mitochondria and chloroplasts increase. This allows cells to uphold feedback regulations to cellular-level signals and adjust to altered environmental conditions. In this review, we summarize recent findings on plant proteolytic complexes with respect to their protective functions against abiotic and biotic stressors.

Salinity effects on nad3 gene RNA editing of wild barley mitochondria

Nad complex plays a very important role during cellular respiration. nad3 (nad dehydrogenase subunit 3) is one of the biggest subunits in this complex. Four cDNAs of nad3 gene were characterized in Hordeum vulgare subsp. spontaneum at exposed to four periods of 500 mM salinity, 0 h or control (accession no. MN066165), after 2 h (accession no. MN066166), after 12 h (accession no. MN066167) and after 24 h (accession no. MN066168) using RNA-seq raw data. Seventeen RNA editing sites were found in positions (or nucleotide nos. C5, C39, C44, C61, C62, C79, C80, C147, C185, C190, C191, C208, C209, C275, C317, C344, C349) within the nad3 coding region. These alterations represent differential editing at four exposure times. The maximum editing rate was revealed 2 and 12 h after salinity exposure. However, these edits were disrupted after 24 h probably due to the initiation of program cell death machinery. We found that RNA editing not only improved protein function but also may improve codon bias by altering the nucleotide without any change in amino acid. Characterization of pentatricopeptide repeat-containing protein At4g13650 (PPRSp1) in wild barley helped us to understand the behavior of editing sites C190 and C191 under salinity. Position - 6 in cis-element upstream editing sites of C155, C190 and C191 may be vital to the editing process in these sites by PPRSp1 protein. The differential editing of this gene under salinity led to a relationship between RNA editing and cellular respiration regulation.

A genome-wide analysis of pentatricopeptide repeat (PPR) protein-encoding genes in four Gossypium species with an emphasis on their expression in floral buds, ovules, and fibers in upland cotton

Cotton is the most important natural fiber used in textiles. Breeding for "three-lines", i.e., cytoplasmic male sterility (CMS)-based sterile (A), maintainer (B), and restorer (R) line, is a promising approach to harness hybrid vigor in cotton. Pentatricopeptide repeat (PPR) protein-encoding genes play an important role in plant growth and development including restoration of CMS plants to male fertility. However, PPRs, especially those contributing to CMS and fiber development, remain largely unknown in cotton. In this study, a genome-wide identification and characterization of PPR gene family in four Gossypium species with genome sequences (G. arboreum, G. raimondii, G. hirsutum, and G. barbadense) were performed, and expressed PPR genes in developing floral buds, ovules, and fibers were compared to identify possible PPRs related to CMS restoration and fiber development. A total of 539, 558, 1032, and 1055 PPRs were predicted in the above four species, respectively, which were further mapped to chromosomes for a synteny analysis. Through an RNA-seq analysis, 86% (882) PPRs were expressed in flowering buds of upland cotton (G. hirsutum); however, only 11 and 6 were differentially expressed (DE) between restorer R and its near-isogenic (NI) B and between R and its NI A line, respectively. Another RNA-seq analysis identified the expression of only 54% (556) PPRs in 0 and 3 day(s) post-anthesis (DPA) ovules and 24% (247) PPRs in 10 DPA fibers; however, only 59, 6, and 27 PPRs were DE in 0 and 3 DPA ovules, and 10 DPA fibers between two backcross inbred lines (BILs) with differing fiber length, respectively. Only 2 PPRs were DE between Xuzhou 142 and its fiberless and fuzzless mutant. Quantitative RT-PCR analysis confirmed the validity of the RNA-seq results for the gene expression pattern. Therefore, only a very small number of PPRs may be associated with fertility restoration of CMS and genetic differences in fiber initiation and elongation. These results lay a foundation for understanding the roles of PPR genes in cotton, and will be useful in the prioritization of candidate PPR gene functional validation for cotton CMS restoration and fiber development.

The Promise of Molecular and Genomic Techniques for Biodiversity Research and DNA Barcoding of the Arabian Peninsula Flora

The Arabian Peninsula is known to have a comprehensive and rich endowment of unique and genetically diverse plant genetic resources. Analysis and conservation of biological diversity is a crucial issue to the whole Arabian Peninsula. The rapid and accurate delimitation and identification of a species is crucial to genetic diversity analysis and the first critical step in the assessment of distribution, population abundance and threats related to a particular target species. During the last two decades, classical strategies of evaluating genetic variability, such as morphology and physiology, have been greatly complemented by phylogenetic, taxonomic, genetic diversity and breeding research molecular studies. At present, initiatives are taking place around the world to generate DNA barcode libraries for vascular plant flora and to make these data available in order to better understand, conserve and utilize biodiversity. The number of herbarium collection-based plant evolutionary genetics and genomics studies being conducted has been increasing worldwide. The herbaria provide a rich resource of already preserved and identified material, and these as well as freshly collected samples from the wild can be used for creating a reference DNA barcode library for the vascular plant flora of a region. This review discusses the main molecular and genomic techniques used in plant identification and biodiversity analysis. Hence, we highlight studies emphasizing various molecular techniques undertaken during the last 10 years to study the plant biodiversity of the Arabian Peninsula. Special emphasis on the role of DNA barcoding as a powerful tool for plant biodiversity analysis is provided, along with the crucial role of herbaria in creating a DNA barcode library.

Molecular and Functional Diversity of RNA Editing in Plant Mitochondria

RNA editing is a fundamental biochemical process relating to the modification of nucleotides in messenger RNAs of functional genes in cells. RNA editing leads to re-establishment of conserved amino acid residues for functional proteins in nuclei, chloroplasts, and mitochondria. Identification of RNA editing factors that contributes to target site recognition increases our understanding of RNA editing mechanisms. Significant progress has been made in recent years in RNA editing studies for both animal and plant cells. RNA editing in nuclei and mitochondria of animal cells and in chloroplast of plant cells has been extensively documented and reviewed. RNA editing has been also extensively documented on plant mitochondria. However, functional diversity of RNA editing factors in plant mitochondria is not overviewed. Here, we review the biological significance of RNA editing, recent progress on the molecular mechanisms of RNA editing process, and function diversity of editing factors in plant mitochondrial research. We will focus on: (1) pentatricopeptide repeat proteins in Arabidopsis and in crop plants; (2) the progress of RNA editing process in plant mitochondria; (3) RNA editing-related RNA splicing; (4) RNA editing associated flower development; (5) RNA editing modulated male sterile; (6) RNA editing-regulated cell signaling; and (7) RNA editing involving abiotic stress. Advances described in this review will be valuable in expanding our understanding in RNA editing. The diverse functions of RNA editing in plant mitochondria will shed light on the investigation of molecular mechanisms that underlies plant development and abiotic stress tolerance.

Stepwise response of MeJA-induced genes and pathways in leaves of C. roseus

Catharanthus roseus is a perennial herb known for the production of important terpenoid indole alkaloids (TIAs) in addition to a variety of phenolic compounds. The goal of the present work was to detect the prolonged effects of MeJA (6 uM) treatment across time (up to 24 days) in order to detect the stepwise response of MeJA-induced genes and pathways in leaves of C. rouses. Prolonged exposure of plants to MeJA (6 uM) treatment for different time points (6, 12 and 24 days) indicated that genes in the indole alkaloid biosynthesis pathway and upstream pathways were triggered earlier (e.g., 6 days) than those in the anthocyanin biosynthesis pathway and its upstream pathways (e.g., 12 days). Three enzymes, e.g., T16H, OMT, and D4H, in the six-step vindoline biosynthesis and two enzymes, e.g., TDC and STR, acting consecutively in the conversion of tryptophan to strictosidine, were activated after 6 days of MeJA treatment. Two other key enzymes, e.g., TRP and CYP72A1, acting concurrently upstream of the TIA biosynthesis pathway were upregulated after 6 days. The genes encoding TDC and STR might concurrently act as a master switch of the TIA pathway towards the production of the indole alkaloids. On the other hand, we speculate that the gene encoding PAL enzyme also acts as the master switch of phenylpropanoid biosynthesis and the downstream flavonoid biosynthesis and anthocyanin biosynthesis pathways towards the production of several phenolic compounds. PAL and the downstream enzymes were activated 12 days after treatment. Cluster analysis confirmed the concordant activities of the flower- and silique-specific bHLH25 transcription factor and the key enzyme in the TIA biosynthesis pathway, e.g., STR. Due to the stepwise response of the two sets of pathways, we speculate that enzymes activated earlier likely make TIA biosynthesis pathway a more favourable target in C. roseus than anthocyanin biosynthesis pathway.

Strengthening desert plant biotechnology research in the United Arab Emirates: a viewpoint

The biotechnology of desert plants is a vast subject. The main applications in this broad field of study comprises of plant tissue culture, genetic engineering, molecular markers and others. Biotechnology applications have the potential to address biodiversity conservation as well as agricultural, medicinal, and environmental issues. There is a need to increase our knowledge of the genetic diversity through the use of molecular genetics and biotechnological approaches in desert plants in the Arabian Gulf region including those in the United Arab Emirates (UAE). This article provides a prospective research for the study of UAE desert plant diversity through DNA fingerprinting as well as understanding the mechanisms of both abiotic stress resistance (including salinity, drought and heat stresses) and biotic stress resistance (including disease and insect resistance). Special attention is given to the desert halophytes and their utilization to alleviate the salinity stress, which is one of the major challenges in agriculture. In addition, symbioses with microorganisms are thought to be hypothesized as important components of desert plant survival under stressful environmental conditions. Thus, factors shaping the diversity and functionality of plant microbiomes in desert ecosystems are also emphasized in this article. It is important to establish a critical mass for biotechnology research and applications while strengthening the channels for collaboration among research/academic institutions in the area of desert plant biotechnology.

Transcriptomic studies reveal a key metabolic pathway contributing to a well-maintained photosynthetic system under drought stress in foxtail millet (Setaria italica L.)

Drought stress is one of the most important abiotic factors limiting crop productivity. A better understanding of the effects of drought on millet (Setaria italica L.) production, a model crop for studying drought tolerance, and the underlying molecular mechanisms responsible for drought stress responses is vital to improvement of agricultural production. In this study, we exposed the drought resistant F1 hybrid, M79, and its parental lines E1 and H1 to drought stress. Subsequent physiological analysis demonstrated that M79 showed higher photosynthetic energy conversion efficiency and drought tolerance than its parents. A transcriptomic study using leaves collected six days after drought treatment, when the soil water content was about ∼20%, identified 3066, 1895, and 2148 differentially expressed genes (DEGs) in M79, E1 and H1 compared to the respective untreated controls, respectively. Further analysis revealed 17 Gene Ontology (GO) enrichments and 14 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways in M79, including photosystem II (PSII) oxygen-evolving complex, peroxidase (POD) activity, plant hormone signal transduction, and chlorophyll biosynthesis. Co-regulation analysis suggested that these DEGs in M79 contributed to the formation of a regulatory network involving multiple biological processes and pathways including photosynthesis, signal transduction, transcriptional regulation, redox regulation, hormonal signaling, and osmotic regulation. RNA-seq analysis also showed that some photosynthesis-related DEGs were highly expressed in M79 compared to its parental lines under drought stress. These results indicate that various molecular pathways, including photosynthesis, respond to drought stress in M79, and provide abundant molecular information for further analysis of the underlying mechanism responding to this stress.

Transcriptomic studies reveal a key metabolic pathway contributing to a well-maintained photosynthetic system under drought stress in foxtail millet (Setaria italica L.)

Drought stress is one of the most important abiotic factors limiting crop productivity. A better understanding of the effects of drought on millet (Setaria italica L.) production, a model crop for studying drought tolerance, and the underlying molecular mechanisms responsible for drought stress responses is vital to improvement of agricultural production. In this study, we exposed the drought resistant F₁ hybrid, M79, and its parental lines E1 and H1 to drought stress. Subsequent physiological analysis demonstrated that M79 showed higher photosynthetic energy conversion efficiency and drought tolerance than its parents. A transcriptomic study using leaves collected six days after drought treatment, when the soil water content was about ∼20%, identified 3066, 1895, and 2148 differentially expressed genes (DEGs) in M79, E1 and H1 compared to the respective untreated controls, respectively. Further analysis revealed 17 Gene Ontology (GO) enrichments and 14 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways in M79, including photosystem II (PSII) oxygen-evolving complex, peroxidase (POD) activity, plant hormone signal transduction, and chlorophyll biosynthesis. Co-regulation analysis suggested that these DEGs in M79 contributed to the formation of a regulatory network involving multiple biological processes and pathways including photosynthesis, signal transduction, transcriptional regulation, redox regulation, hormonal signaling, and osmotic regulation. RNA-seq analysis also showed that some photosynthesis-related DEGs were highly expressed in M79 compared to its parental lines under drought stress. These results indicate that various molecular pathways, including photosynthesis, respond to drought stress in M79, and provide abundant molecular information for further analysis of the underlying mechanism responding to this stress.