Dynamic response of RNA editing to temperature in grape by RNA deep sequencing

PlantC2U: deep learning of cross-species sequence landscapes predicts plastid C-to-U RNA editing in plants

In plants, C-to-U RNA editing mainly occurs in plastid and mitochondrial transcripts, which contributes to a complex transcriptional regulatory network. More evidence reveals that RNA editing plays critical roles in plant growth and development. However, accurate detection of RNA editing sites using transcriptome sequencing data alone is still challenging. In the present study, we develop PlantC2U, which is a convolutional neural network, to predict plastid C-to-U RNA editing based on the genomic sequence. PlantC2U achieves >95% sensitivity and 99% specificity, which outperforms the PREPACT tool, random forests, and support vector machines. PlantC2U not only further checks RNA editing sites from transcriptome data to reduce possible false positives, but also assesses the effect of different mutations on C-to-U RNA editing based on the flanking sequences. Moreover, we found the patterns of tissue-specific RNA editing in the mangrove plant Kandelia obovata, and observed reduced C-to-U RNA editing rates in the cold stress response of K. obovata, suggesting their potential regulatory roles in plant stress adaptation. In addition, we present RNAeditDB, available online at https://jasonxu.shinyapps.io/RNAeditDB/. Together, PlantC2U and RNAeditDB will help researchers explore the RNA editing events in plants and thus will be of broad utility for the plant research community.

Multi-scale analysis of heat stress acclimation in Arabidopsis seedlings highlights the primordial contribution of energy-transducing organelles

Much progress has been made in understanding the molecular mechanisms of plant adaptation to heat stress. However, the great diversity of models and stress conditions, and the fact that analyses are often limited to a small number of approaches, complicate the picture. We took advantage of a liquid culture system in which Arabidopsis seedlings are arrested in their development, thus avoiding interference with development and drought stress responses, to investigate through an integrative approach seedlings' global response to heat stress and acclimation. Seedlings perfectly tolerate a noxious heat shock (43°C) when subjected to a heat priming treatment at a lower temperature (38°C) the day before, displaying a thermotolerance comparable to that previously observed for Arabidopsis. A major effect of the pre-treatment was to partially protect energy metabolism under heat shock and favor its subsequent rapid recovery, which was correlated with the survival of seedlings. Rapid recovery of actin cytoskeleton and mitochondrial dynamics were another landmark of heat shock tolerance. The omics confirmed the role of the ubiquitous heat shock response actors but also revealed specific or overlapping responses to priming, heat shock, and their combination. Since only a few components or functions of chloroplast and mitochondria were highlighted in these analyses, the preservation and rapid recovery of their bioenergetic roles upon acute heat stress do not require extensive remodeling of the organelles. Protection of these organelles is rather integrated into the overall heat shock response, thus allowing them to provide the energy required to elaborate other cellular responses toward acclimation.

Comparative analysis of the organelle genomes of Aconitum carmichaelii revealed structural and sequence differences and phylogenetic relationships

In this study, we conducted an assembly and analysis of the organelle genomes of Aconitum carmichaelii. Our investigation encompassed the examination of organelle genome structures, gene transfer events, and the environmental selection pressures affecting A. carmichaelii. The results revealed distinct evolutionary patterns in the organelle genomes of A. carmichaelii. Especially, the plastome exhibited a more conserved structure but a higher nucleotide substitution rate (NSR), while the mitogenome displayed a more complex structure with a slower NSR. Through homology analysis, we identified several instances of unidirectional protein-coding genes (PCGs) transferring from the plastome to the mitogenome. However, we did not observe any events which genes moved from the mitogenome to the plastome. Additionally, we observed multiple transposable element (TE) fragments in the organelle genomes, with both organelles showing different preferences for the type of nuclear TE insertion. Divergence time estimation suggested that rapid differentiation occurred in Aconitum species approximately 7.96 million years ago (Mya). This divergence might be associated with the reduction in CO2 levels and the significant uplift of the Qinghai-Tibet Plateau (QTP) during the late Miocene. Selection pressure analysis indicated that the dN/dS values of both organelles were less than 1, suggested that organelle PCGs were subject to purification selection. However, we did not detect any positively selected genes (PSGs) in Subg. Aconitum and Subg. Lycoctonum. This observation further supports the idea that stronger negative selection pressure on organelle genes in Aconitum results in a more conserved amino acid sequence. In conclusion, this study contributes to a deeper understanding of organelle evolution in Aconitum species and provides a foundation for future research on the genetic mechanisms underlying the structure and function of the Aconitum plastome and mitogenome.

The emerging role of epitranscriptome in shaping stress responses in plants

KEY MESSAGE

RNA modifications and editing changes constitute 'epitranscriptome' and are crucial in regulating the development and stress response in plants. Exploration of the epitranscriptome and associated machinery would facilitate the engineering of stress tolerance in crops. RNA editing and modifications post-transcriptionally decorate almost all classes of cellular RNAs, including tRNAs, rRNAs, snRNAs, lncRNAs and mRNAs, with more than 170 known modifications, among which m6A, Ψ, m5C, 8-OHG and C-to-U editing are the most abundant. Together, these modifications constitute the "epitranscriptome", and contribute to changes in several RNA attributes, thus providing an additional structural and functional diversification to the "cellular messages" and adding another layer of gene regulation in organisms, including plants. Numerous evidences suggest that RNA modifications have a widespread impact on plant development as well as in regulating the response of plants to abiotic and biotic stresses. High-throughput sequencing studies demonstrate that the landscapes of m6A, m5C, Am, Cm, C-to-U, U-to-G, and A-to-I editing are remarkably dynamic during stress conditions in plants. GO analysis of transcripts enriched in Ψ, m6A and m5C modifications have identified bonafide components of stress regulatory pathways. Furthermore, significant alterations in the expression pattern of genes encoding writers, readers, and erasers of certain modifications have been documented when plants are grown in challenging environments. Notably, manipulating the expression levels of a few components of RNA editing machinery markedly influenced the stress tolerance in plants. We provide updated information on the current understanding on the contribution of RNA modifications in shaping the stress responses in plants. Unraveling of the epitranscriptome has opened new avenues for designing crops with enhanced productivity and stress resilience in view of global climate change.

A Genome-Wide Analysis of the Pentatricopeptide Repeat Protein Gene Family in Two Kiwifruit Species with an Emphasis on the Role of RNA Editing in Pathogen Stress

Kiwifruit is a perennial fruit tree with high nutritional and economic value; however, various pathogen stresses have resulted in reductions in its yield and quality. Pentatricopeptide repeat proteins (PPRs), characterized by tandem repetitions of 35 amino acid motifs, play roles in RNA editing, mRNA stability, and splicing. They may also regulate plant development and growth. Nevertheless, the roles of PPRs in plant development and disease resistance remain unclear. In this study, we focused on the roles of PPRs in the fruit development and pathogen stress of kiwifruit and conducted a series of analyses of the PPR gene family in two representative kiwifruit species (Actinidia chinensis (Ach) and Actinidia eriantha (Ace)) with markedly different degrees of disease resistance. A total of 497 and 499 PPRs were identified in Ach and Ace, respectively. All the kiwifruit PPRs could be phylogenetically divided into four subfamilies. There were about 40.68% PPRs predicted to be localized to mitochondria or chloroplasts. A synteny analysis showed that the expansion of the kiwifruit PPRs mainly originated from segmental duplication. Based on RNA-seq data from the fruit over 12 periods of development and maturity, a weighted correlation network analysis suggested that two PPRs, Actinidia20495.t1 and Actinidia15159.t1, may be involved in fruit development and maturation. In addition, we observed different responses with respect to the expression of PPRs and RNA editing between resistant and susceptible kiwifruits following infection with pathogenic bacteria, indicating the regulatory role of PPRs in the stress response via the modulation of RNA editing. The differentially expressed upstream transcription factors of the PPRs were further identified; they may regulate resistance adaption by modulating the expression of the PPRs. Collectively, these results suggest that PPRs play roles in the development and disease resistance of kiwifruit and provide candidate genes for further clarifying the resistance mechanisms in kiwifruits.

Chloroplast Pan-Genomes and Comparative Transcriptomics Reveal Genetic Variation and Temperature Adaptation in the Cucumber

Although whole genome sequencing, genetic variation mapping, and pan-genome studies have been done on a large group of cucumber nuclear genomes, organelle genome information is largely unclear. As an important component of the organelle genome, the chloroplast genome is highly conserved, which makes it a useful tool for studying plant phylogeny, crop domestication, and species adaptation. Here, we have constructed the first cucumber chloroplast pan-genome based on 121 cucumber germplasms, and investigated the genetic variations of the cucumber chloroplast genome through comparative genomic, phylogenetic, haplotype, and population genetic structure analysis. Meanwhile, we explored the changes in expression of cucumber chloroplast genes under high- and low-temperature stimulation via transcriptome analysis. As a result, a total of 50 complete chloroplast genomes were successfully assembled from 121 cucumber resequencing data, ranging in size from 156,616-157,641 bp. The 50 cucumber chloroplast genomes have typical quadripartite structures, consisting of a large single copy (LSC, 86,339-86,883 bp), a small single copy (SSC, 18,069-18,363 bp), and two inverted repeats (IRs, 25,166-25,797 bp). Comparative genomic, haplotype, and population genetic structure results showed that there is more genetic variation in Indian ecotype cucumbers compared to other cucumber cultivars, which means that many genetic resources remain to be explored in Indian ecotype cucumbers. Phylogenetic analysis showed that the 50 cucumber germplasms could be classified into 3 types: East Asian, Eurasian + Indian, and Xishuangbanna + Indian. The transcriptomic analysis showed that matK were significantly up-regulated under high- and low-temperature stresses, further demonstrating that cucumber chloroplasts respond to temperature adversity by regulating lipid metabolism and ribosome metabolism. Further, accD has higher editing efficiency under high-temperature stress, which may contribute to the heat tolerance. These studies provide useful insight into genetic variation in the chloroplast genome, and established the foundation for exploring the mechanisms of temperature-stimulated chloroplast adaptation.

Identification of genetic and biochemical mechanisms associated with heat shock and heat stress adaptation in grain amaranths

Heat stress is poised to become a major factor negatively affecting plant performance worldwide. In terms of world food security, increased ambient temperatures are poised to reduce yields in cereals and other economically important crops. Grain amaranths are known to be productive under poor and/or unfavorable growing conditions that significantly affect cereals and other crops. Several physiological and biochemical attributes have been recognized to contribute to this favorable property, including a high water-use efficiency and the activation of a carbon starvation response. This study reports the behavior of the three grain amaranth species to two different stress conditions: short-term exposure to heat shock (HS) conditions using young plants kept in a conditioned growth chamber or long-term cultivation under severe heat stress in greenhouse conditions. The latter involved exposing grain amaranth plants to daylight temperatures that hovered around 50°C, or above, for at least 4 h during the day and to higher than normal nocturnal temperatures for a complete growth cycle in the summer of 2022 in central Mexico. All grain amaranth species showed a high tolerance to HS, demonstrated by a high percentage of recovery after their return to optimal growing conditions. The tolerance observed coincided with increased expression levels of unknown function genes previously shown to be induced by other (a)biotic stress conditions. Included among them were genes coding for RNA-binding and RNA-editing proteins, respectively. HS tolerance was also in accordance with favorable changes in several biochemical parameters usually induced in plants in response to abiotic stresses. Conversely, exposure to a prolonged severe heat stress seriously affected the vegetative and reproductive development of all three grain amaranth species, which yielded little or no seed. The latter data suggested that the usually stress-tolerant grain amaranths are unable to overcome severe heat stress-related damage leading to reproductive failure.

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.

The PPR protein RARE1-mediated editing of chloroplast accD transcripts is required for fatty acid biosynthesis and heat tolerance in Arabidopsis

It has been reported that Arabidopsis chloroplast accD transcripts undergo RNA editing and that loss of accD-C794 RNA editing does not affect plant growth under normal conditions. To date, the exact biological role of accD-C794 editing has remained elusive. Here, we reveal an unexpected role for accD-C794 editing in response to heat stress. Loss of accD-C794 editing results in a yellow and dwarf phenotype with decreased chloroplast gene expression under heat stress, and artificial improvement of C794-edited accD gene expression enhances heat tolerance in Arabidopsis. These data suggest that accD-C794 editing confers heat tolerance in planta. We also found that treatment with the product of acetyl coenzyme A carboxylase (ACCase) could allay mutant phenotypic characteristics and showed that a mutation in the CAC3 gene for the α-subunit of ACCase was associated with dwarfism under heat stress. These observations indicate that defective accD-C794 editing may be intrinsic to reduced ACCase activity, thereby contributing to heat sensitivity. ACCase catalyzes the committed step of de novo fatty acid (FA) biosynthesis. FA content analysis revealed that unsaturated oleic (C18:1) and linoleic acids (C18:2) were low in the accD-C794 editing-defective mutant but high in the C794-edited accD-overexpressing plants compared with the wild type. Supplying exogenous C18:1 and C18:2 could rescue the mutant phenotype, suggesting that these FAs play an essential role in tolerance to heat stress. Transmission electron microscopy observations showed that heat stress seriously affected the membrane architecture in accD editing-defective mutants but not in accD-overexpressing plants. These results provide the first evidence that accD-C794 editing regulates FA biosynthesis for maintenance of membrane structural homeostasis under heat stress.

Genome-wide identification and expression analysis of peach multiple organellar RNA editing factors reveals the roles of RNA editing in plant immunity

BACKGROUND

Multiple organellar RNA editing factor (MORF) genes play key roles in chloroplast developmental processes by mediating RNA editing of Cytosine-to-Uracil conversion. However, the function of MORF genes in peach (Prunus persica), a perennial horticultural crop species of Rosaceae, is still not well known, particularly the resistance to biotic and abiotic stresses that threaten peach yield seriously.

RESULTS

In this study, to reveal the regulatory roles of RNA editing in plant immunity, we implemented genome-wide analysis of peach MORF (PpMORF) genes in response to biotic and abiotic stresses. The chromosomal and subcellular location analysis showed that the identified seven PpMORF genes distributed on three peach chromosomes were mainly localized in the mitochondria and chloroplast. All the PpMORF genes were classified into six groups and one pair of PpMORF genes was tandemly duplicated. Based on the meta-analysis of two types of public RNA-seq data under different treatments (biotic and abiotic stresses), we observed down-regulated expression of PpMORF genes and reduced chloroplast RNA editing, especially the different response of PpMORF2 and PpMORF9 to pathogens infection between resistant and susceptible peach varieties, indicating the roles of MORF genes in stress response by modulating the RNA editing extent in plant immunity. Three upstream transcription factors (MYB3R-1, ZAT10, HSFB3) were identified under both stresses, they may regulate resistance adaption by modulating the PpMORF gene expression.

CONCLUSION

These results provided the foundation for further analyses of the functions of MORF genes, in particular the roles of RNA editing in plant immunity. In addition, our findings will be conducive to clarifying the resistance mechanisms in peaches and open up avenues for breeding new cultivars with high resistance.

RNA Editing in Chloroplast: Advancements and Opportunities

Many eukaryotic and prokaryotic organisms employ RNA editing (insertion, deletion, or conversion) as a post-transcriptional modification mechanism. RNA editing events are common in these organelles of plants and have gained particular attention due to their role in the development and growth of plants, as well as their ability to cope with abiotic stress. Owing to rapid developments in sequencing technologies and data analysis methods, such editing sites are being accurately predicted, and many factors that influence RNA editing are being discovered. The mechanism and role of the pentatricopeptide repeat protein family of proteins in RNA editing are being uncovered with the growing realization of accessory proteins that might help these proteins. This review will discuss the role and type of RNA editing events in plants with an emphasis on chloroplast RNA editing, involved factors, gaps in knowledge, and future outlooks.

Genome-wide detection of RNA editing events during the hair follicles cycle of Tianzhu white yak

BACKGROUND

The hair coat is available for the yak to live in the harsh environment of the plateau. Besides, improving the hair production of yak is necessary for its textile industry development. Hair grows from hair follicles (HFs). The HFs undergo periodic growth after birth and are regulated by the complex gene regulatory network. However, the molecular mechanism of HFs regeneration in the Tianzhu white yak remains unclear. RNA editing is a post-transcriptional mechanism that regulates gene expression and produces new transcripts. Hence, we investigated the influence of the A-to-I RNA editing events on the HFs cycle of the Tianzhu white yak.

RESULTS

We finally identified 54,707 adenosine-to-inosine (A-to-I) RNA editing sites (RESs) from RNA sequencing data of the HFs cycle in the Tianzhu white yak. Annotation results showed RESs caused missense amino acid changes in 7 known genes. And 202 A-to-I editing sites altered 23 target genes of 140 microRNAs. A total of 1,722 differential RESs were identified during the HFs cycle of Tianzhu white yak. GO and KEGG enrichment analysis revealed several signaling pathways and GO terms involved skin development, hair growth, and HFs cycle. Such as genes with differential RNA editing levels were significantly enriched in the peroxisome, metabolic pathways, Notch signaling pathway, and PPAR signaling pathway. Besides, the editing sites in HFs development-related genes FAS, APCDD1, WWOX, MPZL3, RUNX1, KANK2, DCN, DSC2, LEPR, HEPHL1, and PTK2B were suggested as the potential RESs involving HFs development.

CONCLUSION

This study investigated the global A-to-I RNA editing events during the HFs cycle of yak skin tissue and expanded the knowledge of A-to-I RNA editing on the HFs cycle. Furthermore, this study revealed that RNA editing-influenced genes may regulate the HFs cycle by participating in the HFs development-related pathways. The findings might provide new insight into the regulation of RNA editing in hair growth.

Genome-Wide Identification and Characterization of RNA/DNA Differences Associated with Fusarium graminearum Infection in Wheat

RNA/DNA difference (RDD) is a post-transcriptional modification playing a crucial role in regulating diverse biological processes in eukaryotes. Although it has been extensively studied in plant chloroplast and mitochondria genomes, RDDs in plant nuclear genomes are not well studied at present. Here, we investigated the RDDs associated with fusarium head blight (FHB) through a novel method by comparing the RNA-seq data between Fusarium-infected and control samples of four wheat genotypes. A total of 187 high-confidence unique RDDs in 36 genes were identified, representing the first landscape of the FHB-responsive RDD in wheat. The majority (26) of these 36 RDD genes were correlated either positively or negatively with FHB levels. Effects of these RDDs on RNA and protein sequences have been identified, their editing frequency and the expression level of the corresponding genes provided, and the prediction of the effect on the minimum folding free energy of mRNA, miRNA binding, and colocation of RDDs with conserved domains presented. RDDs were predicted to induce modifications in the mRNA and protein structures of the corresponding genes. In two genes, TraesCS1B02G294300 and TraesCS3A02G263900, editing was predicted to enhance their affinity with tae-miR9661-5p and tae-miR9664-3p, respectively. To our knowledge, this study is the first report of the association between RDD and FHB in wheat; this will contribute to a better understanding of the molecular basis underlying FHB resistance, and potentially lead to novel strategies to improve wheat FHB resistance through epigenetic methods.

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.

Genome-Wide Analysis of Multiple Organellar RNA Editing Factor (MORF) Family in Kiwifruit (Actinidia chinensis) Reveals Its Roles in Chloroplast RNA Editing and Pathogens Stress

Kiwifruit (Actinidia chinensis) is well known for its high vitamin C content and good taste. Various diseases, especially bacterial canker, are a serious threat to the yield of kiwifruit. Multiple organellar RNA editing factor (MORF) genes are pivotal factors in the RNA editosome that mediates Cytosine-to-Uracil RNA editing, and they are also indispensable for the regulation of chloroplast development, plant growth, and response to stresses. Although the kiwifruit genome has been released, little is known about MORF genes in kiwifruit at the genome-wide level, especially those involved in the response to pathogens stress. In this study, we identified ten MORF genes in the kiwifruit genome. The genomic structures and chromosomal locations analysis indicated that all the MORF genes consisted of three conserved motifs, and they were distributed widely across the seven linkage groups and one contig of the kiwifruit genome. Based on the structural features of MORF proteins and the topology of the phylogenetic tree, the kiwifruit MORF gene family members were classified into six groups (Groups A-F). A synteny analysis indicated that two pairs of MORF genes were tandemly duplicated and five pairs of MORF genes were segmentally duplicated. Moreover, based on analysis of RNA-seq data from five tissues of kiwifruit, we found that both expressions of MORF genes and chloroplast RNA editing exhibited tissue-specific patterns. MORF2 and MORF9 were highly expressed in leaf and shoot, and may be responsible for chloroplast RNA editing, especially the ndhB genes. We also observed different MORF expression and chloroplast RNA editing profiles between resistant and susceptible kiwifruits after pathogen infection, indicating the roles of MORF genes in stress response by modulating the editing extend of mRNA. These results provide a solid foundation for further analyses of the functions and molecular evolution of MORF genes, in particular, for clarifying the resistance mechanisms in kiwifruits and breeding new cultivars with high resistance.

RNA Editing and Its Roles in Plant Organelles

RNA editing, a vital supplement to the central dogma, yields genetic information on RNA products that are different from their DNA templates. The conversion of C-to-U in mitochondria and plastids is the main kind of RNA editing in plants. Various factors have been demonstrated to be involved in RNA editing. In this minireview, we summarized the factors and mechanisms involved in RNA editing in plant organelles. Recently, the rapid development of deep sequencing has revealed many RNA editing events in plant organelles, and we further reviewed these events identified through deep sequencing data. Numerous studies have shown that RNA editing plays essential roles in diverse processes, such as the biogenesis of chloroplasts and mitochondria, seed development, and stress and hormone responses. Finally, we discussed the functions of RNA editing in plant organelles.

Pi-starvation induced transcriptional changes in barley revealed by a comprehensive RNA-Seq and degradome analyses

BACKGROUND

Small RNAs (sRNAs) are 20-30 nt regulatory elements which are responsible for plant development regulation and participate in many plant stress responses. Insufficient inorganic phosphate (Pi) concentration triggers plant responses to balance the internal Pi level.

RESULTS

In this study, we describe Pi-starvation-responsive small RNAs and transcriptome changes in barley (Hordeum vulgare L.) using Next-Generation Sequencing (NGS) RNA-Seq data derived from three different types of NGS libraries: (i) small RNAs, (ii) degraded RNAs, and (iii) functional mRNAs. We find that differentially and significantly expressed miRNAs (DEMs, Bonferroni adjusted p-value < 0.05) are represented by 15 molecules in shoot and 13 in root; mainly various miR399 and miR827 isomiRs. The remaining small RNAs (i.e., those without perfect match to reference sequences deposited in miRBase) are considered as differentially expressed other sRNAs (DESs, p-value Bonferroni correction < 0.05). In roots, a more abundant and diverse set of other sRNAs (DESs, 1796 unique sequences, 0.13% from the average of the unique small RNA expressed under low-Pi) contributes more to the compensation of low-Pi stress than that in shoots (DESs, 199 unique sequences, 0.01%). More than 80% of differentially expressed other sRNAs are up-regulated in both organs. Additionally, in barley shoots, up-regulation of small RNAs is accompanied by strong induction of two nucleases (S1/P1 endonuclease and 3'-5' exonuclease). This suggests that most small RNAs may be generated upon nucleolytic cleavage to increase the internal Pi pool. Transcriptomic profiling of Pi-starved barley shoots identifies 98 differentially expressed genes (DEGs). A majority of the DEGs possess characteristic Pi-responsive cis-regulatory elements (P1BS and/or PHO element), located mostly in the proximal promoter regions. GO analysis shows that the discovered DEGs primarily alter plant defense, plant stress response, nutrient mobilization, or pathways involved in the gathering and recycling of phosphorus from organic pools.

CONCLUSIONS

Our results provide comprehensive data to demonstrate complex responses at the RNA level in barley to maintain Pi homeostasis and indicate that barley adapts to Pi-starvation through elicitation of RNA degradation. Novel P-responsive genes were selected as putative candidates to overcome low-Pi stress in barley plants.

white panicle2 encoding thioredoxin z, regulates plastid RNA editing by interacting with multiple organellar RNA editing factors in rice

Thioredoxins (TRXs) occur in plant chloroplasts as complex disulphide oxidoreductases. Although many biological processes are regulated by thioredoxins, the regulatory mechanism of chloroplast TRXs are largely unknown. Here we report a rice white panicle2 mutant caused by a mutation in the thioredoxin z gene, an orthologue of AtTRX z in Arabidopsis. white panicle2 (wp2) seedlings exhibited a high-temperature-sensitive albinic phenotype. We found that plastid multiple organellar RNA editing factors (MORFs) were the regulatory targets of thioredoxin z. We showed that OsTRX z protein physically interacts with OsMORFs in a redox-dependent manner and that the redox state of a conserved cysteine in the MORF box is essential for MORF-MORF interactions. wp2 and OsTRX z knockout lines show reduced editing efficiencies in many plastidial-encoded genes especially under high-temperature conditions. An Arabidopsis trx z mutant also exhibited significantly reduced chloroplast RNA editing. Our combined results suggest that thioredoxin z regulates chloroplast RNA editing in plants by controlling the redox state of MORFs.

Heat-Responsive miRNAs Participate in the Regulation of Male Fertility Stability in Soybean CMS-Based F1 under High Temperature Stress

MicroRNAs (miRNAs), a class of noncoding small RNAs (sRNAs), are widely involved in the response to high temperature (HT) stress at both the seedling and flowering stages. To dissect the roles of miRNAs in regulating male fertility in soybean cytoplasmic male sterility (CMS)-based F1 under HT, sRNA sequencing was performed using flower buds from HT-tolerant and HT-sensitive CMS-based F1 combinations (NF1 and YF1, respectively). A total of 554 known miRNAs, 59 new members of known miRNAs, 712 novel miRNAs, and 1145 target genes of 580 differentially expressed miRNAs (DEMs) were identified under normal temperature and HT conditions. Further integrated analysis of sRNA and transcriptome sequencing found that 21 DEMs and 15 differentially expressed target genes, such as gma-miR397a/Laccase 2, gma-miR399a/Inorganic phosphate transporter 1-4, and gma-miR4413a/PPR proteins, mitochondrial-like, were negatively regulated under HT stress. Furthermore, all members of the gma-miR156 family were suppressed by HT stress in both NF1 and YF1, but were highly expressed in YF1 under HT condition. The negative correlation between gma-miR156b and its target gene squamosa promoter-binding protein-like 2b was confirmed by expression analysis, and overexpression of gma-miR156b in Arabidopsis led to male sterility under HT stress. With these results, we proposed that miRNAs play an important role in the regulation of male fertility stability in soybean CMS-based F1 under HT stress.