Genome-Wide Analysis of Multiple Organellar RNA Editing Factor Family in Poplar Reveals Evolution and Roles in Drought Stress

Comprehensive identification, characterization, and expression analysis of the MORF gene family in Brassica napus

BACKGROUND

RNA editing in chloroplast and mitochondrion transcripts of plants is an important type of post-transcriptional RNA modification in which members of the multiple organellar RNA editing factor gene family (MORF) play a crucial role. However, a systematic identification and characterization of MORF members in Brassica napus is still lacking.

RESULTS

In this study, a total of 43 MORF genes were identified from the genome of the Brassica napus cultivar "Zhongshuang 11". The Brassica napus MORF (BnMORF) family members were divided into three groups through phylogenetic analysis. BnMORF genes distributed on 14 chromosomes and expanded due to segmental duplication and whole genome duplication repetitions. The majority of BnMORF proteins were predicted to be localized to mitochondria and chloroplasts. The promoter cis-regulatory element analysis, spatial-temporal expression profiling, and co-expression network of BnMORF genes indicated the involvement of BnMORF genes in stress and phytohormone responses, as well as growth and development.

CONCLUSION

This study provides a comprehensive analysis of BnMORF genes and lays a foundation for further exploring their physiological functions in Brassica napus.

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.

Evolution of the 14-3-3 gene family in monocotyledons and dicotyledons and validation of MdGRF13 function in transgenic Arabidopsis thaliana

KEY MESSAGE

The 14-3-3 family is more highly conserved among monocotyledons, and overexpression of MdGRF13 improved drought and salt tolerance in transgenic Arabidopsis thaliana. The 14-3-3 are highly conserved regulatory proteins found in eukaryotes and play an essential role in plant growth, development and stress response. However, the 14-3-3 gene family evolution in monocotyledons and dicotyledons and the biological functions of the MdGRF13 under abiotic stress remain unknown. In our study, 195 members of the 14-3-3 family were identified from 12 species and divided into ε group and the Non-ε group. Synteny analysis within the 14-3-3 family indicated that segmental duplication events contributed to the expansion of the family. Selective pressure analysis indicated that purifying selection was a vital force in the 14-3-3 genes evolution, and monocotyledons had a lower million years ago (Mya) mean values than dicotyledons. Meanwhile, the codon adaptation index (CAI) and frequency of optical codons (FOP) are higher and the effective number of codons (Nc) is lower in monocotyledons 14-3-3 genes compared to dicotyledons. Moreover, the yeast two-hybrid (Y2H) demonstrated that MdGRF13 interacts with MdRD22, MdLHP1a and MdMORF1. Significantly, the malondialdehyde (MDA) content and relative conductivity were decreased, while the superoxide dismutase (SOD), peroxidase (POD), and catalase (CAT) activities were increased in transgenic Arabidopsis compared to the wild type (WT) under drought and salt stress. These results suggest that overexpression of MdGRF13 significantly improved the tolerance to drought and salt stress in transgenic Arabidopsis. Thus, our results provide a theoretical basis for exploring the evolution and function of the 14-3-3 gene family in monocotyledons and dicotyledons.

Genome-Scale Analysis of the Grapevine KCS Genes Reveals Its Potential Role in Male Sterility

Very long-chain fatty acid (VLCFA) synthesis in plants, is primarily rate-limited by the enzyme 3-ketoacyl CoA synthase (KCS), which also controls the rate and carbon chain length of VLCFA synthesis. Disruption of VLCFA during pollen development, may affect the pollen wall formation and ultimately lead to male sterility. Our study identified 24 grapevine KCS (VvKCS) genes and provided new names based on their relative chromosome distribution. Based on sequence alignment and phylogenetic investigation, these genes were grouped into seven subgroups, members of the same subgroup having similar motif structures. Synteny analysis of VvKCS genes, showed that the segmental duplication events played an important role in expanding this gene family. Expression profiles obtained from the transcriptome data showed different expression patterns of VvKCS genes in different tissues. Comparison of transcriptome and RT-qPCR data of the male sterile grape ‘Y−14’ and its fertile parent ‘Shine Muscat’, revealed that 10 VvKCS genes were significantly differentially expressed at the meiosis stage, which is a critical period of pollen wall formation. Further, joint analysis by weighted gene co-expression network analysis (WGCNA) and Kyoto Encyclopedia of Genes and Genomes (KEGG), revealed that five of these VvKCS (VvKCS6/15/19/20/24) genes were involved in the fatty acid elongation pathway, which may ultimately affect the structural integrity of the pollen wall in ‘Y−14’. This systematic analysis provided a foundation for further functional characterization of VvKCS genes, with the aim of grapevine precision breeding improvement.

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.

Characteristics and molecular identification of glyceraldehyde-3-phosphate dehydrogenases in poplar

Glyceraldehyde-3-phosphate dehydrogenase (GAPDH), an essential enzyme of the glycolysis metabolic pathway, plays a vital role in carbon metabolism, plant development, and stress resistance. As a kind of woody plant, poplars are widely cultivated for afforestation. Although the whole genome data of poplars have been published, little information is known about the GAPDH family of genes in poplar. This study performed a genome-wide identification of the poplar GAPDH family, and 13 determined PtGAPDH genes were identified from poplar genome. Phylogenetic tree showed that the PtGAPDH members were divided into PtGAPA/B, PtGAPC, PtGAPCp, and PtGAPN groups. A total of 13 PtGAPDH genes were distributed on eight chromosomes, 13 gene pairs belonging to segmented replication events were detected in poplar, and 23 collinearity gene pairs were determined between poplar and willow. The PtGAPDHcis-acting elements associated with growth and development as well as stress resistance revealed that PtGAPDHs might be involved in these processes. The phosphoglycerate kinase (PGK) and triose-phosphate isomerase (TPI) were predicted as the putative interaction proteins of PtGAPDHs. Gene ontology (GO) analysis showed that PtGAPDHs play a crucial role in the oxidation and reduction processes. PtGAPDH expression levels were induced by NaCl and PEG treatments, which implied that PtGAPDHs might be involved in stress response. Overexpression of PtGAPC1 significantly changed the contents of lipid and carbohydrate metabolites, which indicated that PtGAPC1 plays an essential role in metabolic regulation. This study highlights the characterizations and profiles of PtGAPDHs and reveals that PtGAPC1 is involved in the loop of lipid and carbohydrate metabolisms.

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.

QTL mapping of drought-related traits in the hybrids of Populus deltoides 'Danhong'×Populus simonii 'Tongliao1'

BACKGROUND

Poplar trees provide a large amount of wood material, but many parts of the world are arid or semi-arid areas because of insufficient annual precipitation, which seriously affects the growth of poplar trees. Populus simonii 'Tongliao1' shows strong tolerance to stress environments, and Populus deltoides 'Danhong' shows a stronger growth rate in a suitable environment. To identify drought tolerance-related QTLs and genes, an F₁ population derived from the cross between the 'Danhong' and 'Tongliao 1' Populus was assessed under drought stress.

RESULTS

We measured drought-related traits such as the relative height growth, relative diameter growth, leaf senescence number, specific leaf area, and leaf relative water content in the population under control and drought environments. The results showed that drought stress reduced the plant height relative growth, ground diameter relative growth, specific leaf area and leaf relative water content and increased the number of leaf drops. A total of 208 QTLs were identified by QTL mapping analysis, and they consisted of 92, 63 and 53 QTLs under control, drought stress treatment and drought index conditions, respectively. A molecular identification marker for drought tolerance, np2841, which was associated with a QTL (qDLRWC-LG10-1) for relative leaf water content, was initially developed. We mined 187 candidate genes for QTL regions of five traits under a drought environment. The reference genome annotation for Populus trichocarpa and a homologous gene analysis of Arabidopsis thaliana identified two candidate genes, Potri.003G171300 and Potri.012G123900, with significant functions in response to drought stress. We identified five key regulatory genes (Potri.006G273500, Potri.007G111500, Potri.007G111600, Potri.007G111700, and Potri.007G111800) related to drought tolerance through the poplar coexpression network.

CONCLUSION

In this study, our results indicate that the QTLs can effectively enhance the drought tolerance of poplar. It is a step closer towards unravelling the genetic basis of poplar drought tolerance-related traits, and to providing validated candidate genes and molecular markers for future genetic improvement.

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.

Genome-wide investigation and expression profiling of polyphenol oxidase (PPO) family genes uncover likely functions in organ development and stress responses in Populus trichocarpa

Background

Trees such as Populus are planted extensively for reforestation and afforestation. However, their successful establishment greatly depends upon ambient environmental conditions and their relative resistance to abiotic and biotic stresses. Polyphenol oxidase (PPO) is a ubiquitous metalloproteinase in plants, which plays crucial roles in mediating plant resistance against biotic and abiotic stresses. Although the whole genome sequence of Populus trichocarpa has long been published, little is known about the PPO genes in Populus, especially those related to drought stress, mechanical damage, and insect feeding. Additionally, there is a paucity of information regarding hormonal responses at the whole genome level.

Results

A genome-wide analysis of the poplar PPO family was performed in the present study, and 18 PtrPPO genes were identified. Bioinformatics and qRT-PCR were then used to analyze the gene structure, phylogeny, chromosomal localization, gene replication, cis-elements, and expression patterns of PtrPPOs. Sequence analysis revealed that two-thirds of the PtrPPO genes lacked intronic sequences. Phylogenetic analysis showed that all PPO genes were categorized into 11 groups, and woody plants harbored many PPO genes. Eighteen PtrPPO genes were disproportionally localized on 19 chromosomes, and 3 pairs of segmented replication genes and 4 tandem repeat genomes were detected in poplars. Cis-acting element analysis identified numerous growth and developmental elements, secondary metabolism processes, and stress-related elements in the promoters of different PPO members. Furthermore, PtrPPO genes were expressed preferentially in the tissues and fruits of young plants. In addition, the expression of some PtrPPOs could be significantly induced by polyethylene glycol, abscisic acid, and methyl jasmonate, thereby revealing their potential role in regulating the stress response. Currently, we identified potential upstream TFs of PtrPPOs using bioinformatics.

Conclusions

Comprehensive analysis is helpful for selecting candidate PPO genes for follow-up studies on biological function, and progress in understanding the molecular genetic basis of stress resistance in forest trees might lead to the development of genetic resources.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08028-9.

OsMORF9 is necessary for chloroplast development and seedling survival in rice

Chloroplasts are closely associated with the growth and development of higher plants. Accumulating evidence has revealed that the multiple organellar RNA editing factors (MORF) family of proteins influences plastidic and mitochondrial development through post-transcriptional regulation. However, the role of MORFs in regulating the development of chloroplasts in rice is still unclear. The OsMORF9 gene belongs to a small family of 7 genes in rice and is highly expressed in young leaves. We used the CRISPR/Cas9 system to mutate OsMORF9. The resulting knockout lines osmorf9-1 and osmorf9-2 exhibited an albino seedling lethal phenotype. Besides, the expression of many plastid-encoded genes involved in photosynthesis, the biogenesis of plastidic ribosomes and the editing and splicing of specific plastidic RNA molecules were severely affected in these two OsMORF9 mutants. Furthermore, yeast two-hybrid analysis revealed that OsMORF9 could interact with OsSLA4 and DUA1 which are members of the pentatricopeptide repeat (PPR) family of proteins. Analysis of subcellular localization of OsMORF9 also suggested that it might function in chloroplasts. The findings from the present study demonstrated the critical role of OsMORF9 in the biogenesis of chloroplast ribosomes, chloroplast development and seedling survival. This therefore provides new insights on the function of MORF proteins in rice.

Analysis of Populus glycosyl hydrolase family I members and their potential role in the ABA treatment and drought stress response

Glycoside hydrolase family 1 (GH1) β-glucosidases (BGLUs) are encoded by a large number of genes and are involved in many developmental processes and stress responses in plants. Due to their importance in plant growth and development, genome-wide analyses have been conducted in the model plant species Arabidopsis thaliana, rice and maize but not in woody plant species, which have important economic and ecological value. In this study, we systematically analyzed Populus BGLUs (PtBGLUs) and demonstrated the involvement of several genes under stress conditions. Forty-four PtBGLUs were identified in Populus databases; these genes were located on 11 chromosomes, and the proteins of several PtBGLU genes were highly similar. More than 90% of PtBGLUs contain three conserved motifs. Collinearity results showed that 44 PtBGLU genes resulted from 12 tandem and 5 segmental duplication events. Phylogenetic analysis revealed that 128 BGLU genes from Populus trichocarpa, A. thaliana and Oryza sativa could be classified into 4 subgroups and subgroup Ⅱ and Ⅳ were differently having PtBGLUs and AtBGLUs. We further investigated whether several PtBGLUs responded to drought stress and ABA treatment, and the results showed that most of the selected BGLU genes were expressed in response to stress, which is consistent with previous studies involving rice and Arabidopsis homologous genes. Large numbers of stress-, hormone-, and development-related elements in the PtBGLU promoters suggest that BGLU genes may be involved in complex networks. Taken together, our results provide valuable information for an improved understanding of β-glucosidase function in woody plants.

Genome-Wide Investigation of the NF-X1 Gene Family in Populus trichocarpa Expression Profiles during Development and Stress

Poplar are planted extensively in reforestation and afforestation. However, their successful establishment largely depends on the environmental conditions of the newly established plantation and their resistance to abiotic as well as biotic stresses. NF-X1, a widespread transcription factor in plants, plays an irreplaceable role in plant growth, development, and stress tolerance. Although the whole genome sequence of Populus trichocarpa has been published for a long time, little is known about the NF-X1 genes in poplar, especially those related to drought stress, mechanical damage, insect feeding, and hormone response at the whole genome level. In this study, whole genome analysis of the poplar NF-X1 family was performed, and 4 PtrNF-X1 genes were identified. Then, bioinformatics analysis and qRT-PCR were applied to analyze the gene structure, phylogeny, chromosomal localization, gene replication, Cis-elements, and expression patterns of PtrNF-X1genes. Sequence analysis revealed that one-quarter of the PtrNF-X1 genes did not contain introns. Phylogenetic analysis revealed that all NF-X1 genes were split into three subfamilies. The number of two pairs of segmented replication genes were detected in poplars. Cis-acting element analysis identified a large number of elements of growth and development and stress-related elements on the promoters of different NF-X1 members. In addition, some PtrNF-X1 could be significantly induced by polyethylene glycol (PEG) and abscisic acid (ABA), thus revealing their potential role in regulating stress response. Comprehensive analysis is helpful in selecting candidate NF-X1 genes for the follow-up study of the biological function, and molecular genetic progress of stress resistance in forest trees provides genetic resources.

Characterization, Expression, and Interaction Analyses of OsMORF Gene Family in Rice

The multiple organellar RNA editing factors (MORF) gene family plays a key role in organelle RNA editing in flowering plants. MORF genes expressions are also affected by abiotic stress. Although seven OsMORF genes have been identified in rice, few reports have been published on their expression patterns in different tissues and under abiotic stress, and OsMORF–OsMORF interactions. In this study, we analyzed the gene structure of OsMORF family genes. The MORF family members were divided into six subgroups in different plants based on phylogenetic analysis. Seven OsMORF genes were highly expressed in leaves. Six and seven OsMORF genes expressions were affected by cold and salt stresses, respectively. OsMORF–OsMORF interaction analysis indicated that OsMORF1, OsMORF8a, and OsMORF8b could each interact with themselves to form homomers. Moreover, five OsMORF proteins were shown to be able to interact with each other, such as OsMORF8a and OsMORF8b interacting with OsMORF1 and OsMORF2b, respectively, to form heteromers. These results provide information for further study of OsMORF gene function.