Abstract
The continually growing human population creates a concomitantly increasing demand for nutritious crops with high yields. Advances in high throughput sequencing technologies have revealed the genetic architecture of major crops. This includes extensive information enabling comprehensive genetic markers for breeding selection, new gene discoveries, and novel gene regulatory strategies for crop editing. RNA structure is an important type of genetic feature, essential for post-transcriptional regulation of gene expression. Here, we summarize recent advances in genome-wide RNA structure studies in crops and review the associated RNA structure-mediated regulation of gene expression. We also discuss the functional importance of those single nucleotide variations that induce large RNA structure disparities. Lastly, we discuss the potential role of RNA structure in crop molecular breeding.
Introduction
Crop breeding to achieve high productivity is hugely important to meet the growing demands of food sustainability and security (Hickey et al., ). Crop breeding strategies are based on the selection of enhanced traits among cultivars including crop yield, quality, ease of cultivation, tolerance to environmental stresses, and resistance against pathogens etc. (Breseghello and Coelho, ; Tian et al., ). Recent advances in high-throughput sequencing technologies offer a greater understanding of the genetic basis of phenotypic variations enabling precise breeding guided by the genomics (Bevan et al., ). Over recent decades, tremendous effort has generated vast genomics resources comprising comprehensive maps of genome variations in most major crops including rice, wheat, sorghum and tomato (Bevan et al., ). These genomic maps have determined millions of single-nucleotide variations (SNVs) across diverse cultivars facilitating association studies with acquired phenotypes (Huang and Han, ). The majority of these identified SNVs are in non-coding regions (e.g., ~97% in rice) rather than in coding regions, raising a key limitation in understanding the functions of these non-coding SNVs and thereby utilizing them to increase the genetic diversity in generating crop lines with diverse traits (Huang and Han, ). This poses the challenges in understanding the functions those non-coding SNVs, assisting the enhancement of genetic diversity.
Emerging evidence has shown that SNVs have the potential to alter RNA structure, offering scope for understanding the functional role of the majority of SNVs (Solem et al., ). Besides sequence content, RNA structure is the other key genetic property that serves an essential role in post-transcriptional regulation of processes, such as RNA degradation, translation, and splicing (Zhang and Ding, ). RNA structure may offer another way to interpret the results of genome-wide association studies (GWAS) (Solem et al., ). Recent advances in high-throughput RNA structure profiling have revealed RNA structure landscapes in crops such as rice and wheat (Deng et al., ; Su et al., ; Yang et al., ), not only opening novel insights into the unique RNA structural complexities in regulating gene expression in crops but also providing the potential of RNA structure-guided molecular breeding.
Here, we provide an overview of these recent discoveries of RNA structure functionalities in crops, focusing on rice and wheat. We summarize the distinctive features of crop RNA structuromes and their involvement in the post-transcriptional regulation of gene expression. We also review how crops adopted RNA structure in response to changing environmental conditions including temperature and nutrients. We emphasize the identification of SNVs in crops, known as riboSNitches, that induce large RNA structure disparities, as well as the corresponding functional impact of these SNVs. We promote discussion of how to consider adding RNA structure as a new perspective in crop breeding design.
Unique RNA Structure Features in Crops
Recent advances in in vivo RNA structure profiling have generated RNA structure information over tens of thousands of genes at single nucleotide resolution (Ding et al., ; Rouskin et al., ; Spitale et al., ; Deng et al., ; Yang et al., ). Two types of chemical probing methods have been successfully applied in crops (Deng et al., ; Su et al., ; Yang et al., ). One is the dimethyl sulfate (DMS)-based method whereby single-stranded A and C are methylated by DMS (Deng et al., ; Su et al., ), whilst the other method is the SHAPE (Selective 2'-Hydroxyl Acylation analyzed by Primer Extension)-based method, whereby the single-strandedness of all four nucleotides are acetylated (Yang et al., ). Both in vivo RNA structuromes in rice and wheat captured the RNA structures for over half of their transcriptomes, providing a comprehensive view of general RNA structure features (Deng et al., ; Yang et al., ). Notably, most mRNAs in vivo did not fold into in silico-predicted structures that are thermostable RNA secondary structures, indicating in vivo RNA structures in crops maintained their flexibility for folding (Deng et al., ; Yang et al., ). Global structure features in rice and wheat are similar to those in Arabidopsis with a single-stranded region upstream of the start codon in the 5'UTR and a triplet periodic trend of RNA structure patterns across CDS regions (Ding et al., ; Deng et al., ; Yang et al., ).
Most major crops such as rice and wheat have high GC-content genomes and transcriptomes (Sorrells et al., ). Considering that RNAs with high GC content tend to fold into strong structures, one would expect that RNA structures in both rice and wheat should be very strong. Surprisingly, the high GC content in mRNAs did not lead to strong RNA structures in vivo (Deng et al., ; Yang et al., ). Notably, the translatome analysis in both rice and wheat revealed that mRNAs with high GC content tend to be highly translated, suggesting that RNA structures in crops might have evolved to maintain some flexibility even with high GC content to ensure active and efficient RNA biological processes such as translation (Deng et al., ; Yang et al., ).
Wheat, as a typical allopolyploid crop, contains two or more divergent genomes (subgenomes) (Yang et al., ). A recent study on the tetraploid durum wheat cultivar, Kronos (2n = 4x = 28, BBAA), found that up to 39.5% of the homoeologs between A and B subgenomes showed significantly different RNA structure features (Yang et al., ). Among them, 50.9% of homoeologs displayed stronger structures in the A subgenome over those in the B subgenome, while 49.1% of homoeologous pairs exhibited weaker structures in the A subgenome compared to the B subgenome (Yang et al., ). This asymmetry at the RNA structure level offers another layer of subgenomic diversity.
Apart from RNA secondary structure, a recent study using selective 2′-hydroxyl acylation with lithium ion-based primer extension with high-throughput sequencing (SHALiPE-seq) identified the global existence of RNA G-quadruplex structure (RG4), a specific RNA tertiary structure in rice (Yang et al., ). SinceRG4 plays an important role in regulating translation, it is suggested that RG4 may be one type of riboregulators across the rice transcriptome (Yang et al., ). Notably, rice has more RG4s with stronger folding status than Arabidopsis which is likely due to the high GC-content transcriptome of rice (Yang et al., ). Interestingly, 121 orthologous gene pairs containing RG4s that are conserved between Arabidopsis and rice suggest that RG4 may be strongly selected during evolution (Yang et al., ).
RNA Structure-Mediated Translational Regulations of Gene Expression in Crops
With the unveiling of RNA structure landscapes in crops, new insights have been exploited in understanding the RNA structure-mediated post-transcriptional regulation of gene expression (Zhang and Ding, ). The translation process is one of the key post-transcriptional regulations of gene expression (Zhang and Ding, ). The in vivo RNA structure landscapes in both rice and wheat revealed that the single-stranded region upstream of the start codon in the 5′UTR was significantly associated with high translation efficiency (Deng et al., ; Yang et al., ). This single-strandedness may promote ribosome binding and/or enhance ribosome initiation. Moreover, stronger triplet periodic trends of RNA structure patterns across CDS regions were observed in those mRNAs with higher translation efficiency, indicating that the pattern might facilitate ribosome elongation (Deng et al., ; Yang et al., ). Notably, in wheat there were significant anti-correlations between average base pair probability (BPP) and translation efficiency in both 5'UTR and CDS regions, but not in 3'UTR region, suggesting that weak structures in these two genic regions enhance translation (Yang et al., ). This significant effect of RNA structure in regulating translation provides potential in crop editing by altering the RNA structure to obtain preferable translational levels.
Recent polysome profiling in tetraploid durum wheat discovered translational subgenome asymmetry (Yang et al., ). Since both A and B subgenomes share a high degree of sequence similarity (~95.42%), the differences of both GC content and codon preference between A and B subgenomes are very subtle, thereby not significantly contributing to subgenomic translational difference (Yang et al., ). In contrast, the subgenomic RNA structure difference was significantly associated with the subgenomic translational difference (Yang et al., ). Notably, this relationship was much closer in those homoeologs with biological functions such as abiotic stress response, biotic stress response, metal ion response, phytohormone signaling, transcriptional and translational regulations (Figure 1A) (Yang et al., ). Therefore, compared to GC content and codon preference, RNA structure may play a more prevalent role in regulating translational subgenome asymmetry (Yang et al., ).
Figure 1
Apart from transcriptome-wide studies, an antisense-mediated translational regulation has been discovered in rice. The antisense RNA (cis-NATpho1.2) interacts with its sense RNA, the phosphate transporter PHOSPHATE1.2 (PHO1.2) at a high GC region across the PHO1.2 start codon (Reis et al.,
RNA Structure-Mediated Degradation in Crops
The folding status of RNA is crucial in determining its thermodynamic properties, thereby its stability (Bevilacqua et al.,
SNVs Are Capable of Changing RNA Structure-Mediated Regulations of Gene Expression in Crops
RNA folding status largely relies on sequence content (Mathews et al.,
One of the studies in rice found that a 1,236 nt long non-coding RNA, referred to as long-day-specific male-fertility–associated RNA (LDMAR) regulated photoperiod-sensitive male sterility (Ding et al.,
Four genes (VRN1-VRN3 and VRN-D4) were previously identified for regulating vernalization in wheat (Kippes et al.,
In addition to these individual examples of the importance of SNV-induced RNA structure alteration, a recent genome-wide study on the in vivo RNA structure landscape in the tetraploid durum wheat cultivar, Kronos, identified 3,564 SNVs which induced large structure disparities, known as riboSNitches between the A and B subgenomes (Yang et al.,
Durum wheat (T. turgidum ssp. durum) is a major cereal grain used for pasta and couscous production, evolved from domesticated emmer wheat (T. turgidum ssp. Dicoccum). Domesticated emmer wheat was derived from wild emmer wheat (T. turgidum ssp. Dicoccoides) in the Fertile Crescent about 10,000 years ago (Maccaferri et al.,
Future Perspectives
These recent studies provide evidence that RNA structure could serve as an important genetic property to guide crop breeding. In contrast to other cis-regulatory motifs such as transcription factor binding sites, influencing transcriptional levels, RNA structure features mostly affect post-transcriptional regulation of gene expression, directly impacting protein synthesis (Gebauer and Hentze,
Figure 2

Schematic of the diversity among phenotypes, genome sequences, RNA structures, and translational levels across natural varieties. Red squares represent the sequence variations. Small blue circles represent small ribosomal subunits, whilst big blue circles indicate large ribosomal subunits.
It was proposed that the subgenome asymmetry in wheat was evolved in control of both vital genes and genes involved in environmental adaptations (Feldman et al.,
The discovery of riboSNitches selected during wheat domestication also offers scope for innovative crop breeding. In the riboSNitch example illustrated above (Figure 1B), the sequences of the homoeologous pair in the wild emmers are the same between A and B subgenomes that lead to the same subgenomic RNA structures, thereby the same subgenomic translational efficiencies (Figure 1B). In contrast, the sequence divergences between A and B subgenomes in the domesticated cultivars of durum wheat populations cause differentiated RNA structures, and subsequent different translational efficiencies (Figure 1B). This riboSNitch could serve as a molecular marker during the selection of breeding lines. Additionally, if this subgenomic divergence was selected during domestication, nucleotide-editing CRISPR technology could precisely modify the RNA structures to obtain greater divergences between A and B subgenomes, leading to a greater differentiation of translational subgenome expression asymmetry, that may offer a strategy for accelerating the domestication process.
Adaptation to changing environmental conditions such as abiotic stress and biotic stress is essential to crop growth. Most environmental factors such as temperature, heavy metals, nutrients, salt, and light affect RNA folding (Chung et al.,
Funding
HZ was supported by the National Natural Science Foundation of China (32170229) and the National Key Research and Development Program of China (2021YFF1000900).
Publisher's Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Statements
Author contributions
WS, LD, and HZ wrote the manuscript together. All authors contributed to the article and approved the submitted version.
Acknowledgments
We thank Dr. Yiliang Ding (John Innes Centre) and Dr. Azahara C. MartÃn (John Innes Centre) for their suggestions.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
References
1
BevanM. W.UauyC.WulffB. B.ZhouJ.KrasilevaK.ClarkM. D. (2017). Genomic innovation for crop improvement. Nature543, 346–354. 10.1038/nature22011
2
BevilacquaP. C.RitcheyL. E.SuZ.AssmannS. M. (2016). Genome-wide analysis of RNA secondary structure. Annu. Rev. Genet.50, 235–266. 10.1146/annurev-genet-120215-035034
3
BreseghelloF.CoelhoA. S. (2013). Traditional and modern plant breeding methods with examples in rice (Oryza sativa L.). J. Agric. Food Chem.61, 8277–8286. 10.1021/jf305531j
4
ChungB. Y. W.BalcerowiczM.Di AntonioM.JaegerK. E.GengF.FranaszekK.et al. (2020). An RNA thermoswitch regulates daytime growth in Arabidopsis. Nat. Plants6, 522–532. 10.1038/s41477-020-0633-3
5
CorleyM.SolemA.PhillipsG.LackeyL.ZiehrB.VincentH. A.et al. (2017). An RNA structure-mediated, posttranscriptional model of human alpha-1-antitrypsin expression. Proc. Natl. Acad. Sci. U. S. A.114, E10244–E53. 10.1073/pnas.1706539114
6
DengH.CheemaJ.ZhangH.WoolfendenH.NorrisM.LiuZ.et al. (2018). Rice in vivo RNA structurome reveals RNA secondary structure conservation and divergence in plants. Mol. Plant. 11, 607–622. 10.1016/j.molp.2018.01.008
7
DingJ.LuQ.OuyangY.MaoH.ZhangP.YaoJ.et al. (2012). A long noncoding RNA regulates photoperiod-sensitive male sterility, an essential component of hybrid rice. Proc. Natl. Acad. Sci. U. S. A.109, 2654–2659. 10.1073/pnas.1121374109
8
DingY.TangY.KwokC. K.ZhangY.BevilacquaP. C.AssmannS. M. (2014). In vivo genome-wide profiling of RNA secondary structure reveals novel regulatory features. Nature505, 696–700. 10.1038/nature12756
9
FeldmanM.LevyA. A.FahimaT.KorolA. (2012). Genomic asymmetry in allopolyploid plants: wheat as a model. J. Exp. Bot.63, 5045–5059. 10.1093/jxb/ers192
10
GebauerF.HentzeM. W. (2004). Molecular mechanisms of translational control. Nat. Rev. Mol. Cell Biol.5, 827–835. 10.1038/nrm1488
11
HickeyL. T.HafeezA. N.RobinsonH.JacksonS. A.Leal-BertioliS. C. M.TesterM.et al. (2019). Breeding crops to feed 10 billion. Nat. Biotechnol.37, 744–754. 10.1038/s41587-019-0152-9
12
HiraoI.EllingtonA. D. (1995). Re-creating the RNA world. Curr. Biol.5, 1017–1022. 10.1016/S0960-9822(95)00205-3
13
HuangX.HanB. (2014). Natural variations and genome-wide association studies in crop plants. Annu. Rev. Plant Biol.65, 531–551. 10.1146/annurev-arplant-050213-035715
14
JabnouneM.SeccoD.LecampionC.RobagliaC.ShuQ.PoirierY. (2013). A rice cis-natural antisense RNA acts as a translational enhancer for its cognate mRNA and contributes to phosphate homeostasis and plant fitness. Plant Cell25, 4166–4182. 10.1105/tpc.113.116251
15
KippesN.DebernardiJ. M.Vasquez-GrossH. A.AkpinarB. A.BudakH.KatoK.et al. (2015). Identification of the VERNALIZATION 4 gene reveals the origin of spring growth habit in ancient wheats from South Asia. Proc. Natl. Acad. Sci. U. S. A.112, E5401–E5410. 10.1073/pnas.1514883112
16
KrajewskiS. S.NarberhausF. (2014). Temperature-driven differential gene expression by RNA thermosensors. Biochim. Biophys. Acta1839, 978–988. 10.1016/j.bbagrm.2014.03.006
17
LeskC.RowhaniP.RamankuttyN. (2016). Influence of extreme weather disasters on global crop production. Nature529, 84–87. 10.1038/nature16467
18
MaccaferriM.HarrisN. S.TwardziokS. O.PasamR. K.GundlachH.SpannaglM.et al. (2019). Durum wheat genome highlights past domestication signatures and future improvement targets. Nat. Genet.51, 885–895. 10.1038/s41588-019-0381-3
19
MathewsD. H.MossW. N.TurnerD. H. (2010). Folding and finding RNA secondary structure. Cold Spring Harb. Perspect. Biol.2, a003665. 10.1101/cshperspect.a003665
20
ReisR. S.DeforgesJ.SchmidtR. R.SchippersJ. H. M.PoirierY. (2021). An antisense noncoding RNA enhances translation via localized structural rearrangements of its cognate mRNA. Plant Cell33, 1381–1397. 10.1093/plcell/koab010
21
RobertsonM. P.JoyceG. F. (2012). The origins of the RNA world. Cold Spring Harb. Perspect. Biol.4, a003608. 10.1101/cshperspect.a003608
22
RouskinS.ZubradtM.WashietlS.KellisM.WeissmanJ. S. (2014). Genome-wide probing of RNA structure reveals active unfolding of mRNA structures in vivo. Nature505, 701–705. 10.1038/nature12894
23
SolemA. C.HalvorsenM.RamosS. B.LaederachA. (2015). The potential of the riboSNitch in personalized medicine. Wiley Interdiscip. Rev. RNA6, 517–532. 10.1002/wrna.1291
24
SorrellsM. E.La RotaM.Bermudez-KandianisC. E.GreeneR. A.KantetyR.MunkvoldJ. D.et al. (2003). Comparative DNA sequence analysis of wheat and rice genomes. Genome Res.13, 1818–1827. 10.1101/gr.1113003
25
SpitaleR. C.FlynnR. A.ZhangQ. C.CrisalliP.LeeB.JungJ. W.et al. (2015). Structural imprints in vivo decode RNA regulatory mechanisms. Nature519, 486–490. 10.1038/nature14263
26
StreitnerC.DanismanS.WehrleF.SchoningJ. C.AlfanoJ. R.StaigerD. (2008). The small glycine-rich RNA binding protein AtGRP7 promotes floral transition in Arabidopsis thaliana. Plant J.56, 239–250. 10.1111/j.1365-313X.2008.03591.x
27
SuZ.TangY.RitcheyL. E.TackD. C.ZhuM.BevilacquaP. C.et al. (2018). Genome-wide RNA structurome reprogramming by acute heat shock globally regulates mRNA abundance. Proc. Natl. Acad. Sci. U. S. A.115, 12170–12175. 10.1073/pnas.1807988115
28
TianZ.WangJ. W.LiJ.HanB. (2021). Designing future crops: challenges and strategies for sustainable agriculture. Plant J.105, 1165–1178. 10.1111/tpj.15107
29
WanY.QuK.ZhangQ. C.FlynnR. A.ManorO.OuyangZ.et al. (2014). Landscape and variation of RNA secondary structure across the human transcriptome. Nature505, 706–709. 10.1038/nature12946
30
XiaoJ.XuS.LiC.XuY.XingL.NiuY.et al. (2014). O-GlcNAc-mediated interaction between VER2 and TaGRP2 elicits TaVRN1 mRNA accumulation during vernalization in winter wheat. Nat. Commun.5, 4572. 10.1038/ncomms5572
31
YangX.CheemaJ.ZhangY.DengH.DuncanS.UmarM. I.et al. (2020). RNA G-quadruplex structures exist and function in vivo in plants. Genome Biol.21, 226. 10.1186/s13059-020-02142-9
32
YangX.YuH.SunW.DingL.LiJ.CheemaJ.et al. (2021). Wheat in vivo RNA structure landscape reveals a prevalent role of RNA structure in modulating translational subgenome expression asymmetry. Genome Biol.22, 326. 10.1186/s13059-021-02549-y
33
ZhangH.DingY. (2021). Novel insights into the pervasive role of RNA structure in post-transcriptional regulation of gene expression in plants. Biochem. Soc. Trans.49, 1829–1839. 10.1042/BST20210318
34
ZhangY.MasselK.GodwinI. D.GaoC. (2018). Applications and potential of genome editing in crop improvement. Genome Biol.19, 210. 10.1186/s13059-018-1586-y
35
ZhuJ.LiC.PengX.ZhangX. (2021). RNA architecture influences plant biology. J. Exp. Bot.72, 4144–4160. 10.1093/jxb/erab030
Summary
Keywords
RNA structure, translation, molecular breeding, crop improvement, riboSNitches
Citation
Sun W, Ding L and Zhang H (2022) The Potential Role of RNA Structure in Crop Molecular Breeding. Front. Plant Sci. 13:868771. doi: 10.3389/fpls.2022.868771
Received
03 February 2022
Accepted
21 March 2022
Published
02 May 2022
Volume
13 - 2022
Edited by
Paula Duque, Gulbenkian Institute of Science (IGC), Portugal
Reviewed by
Zhao Su, The Pennsylvania State University (PSU), United States
Updates

Check for updates
Copyright
© 2022 Sun, Ding and Zhang.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Huakun Zhang zhanghk045@nenu.edu.cn
†These authors have contributed equally to this work
This article was submitted to Plant Cell Biology, a section of the journal Frontiers in Plant Science
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.