sRNAanno-a database repository of uniformly annotated small RNAs in plants

The genomic insights of intertidal adaptation in Bryopsis corticulans

Many marine green algae thrive in intertidal zones, adapting to complex light environments that fluctuate between low underwater light and intense sunlight. Exploring their genomic bases could help to comprehend the diversity of adaptation strategies in response to environmental pressures. Here, we developed a novel and practical strategy to assemble high-confidence algal genomes and sequenced a high-quality genome of Bryopsis corticulans, an intertidal zone macroalga in the Bryopsidales order of Chlorophyta that originated 678 million years ago. Comparative genomic analyses revealed a previously overlooked whole genome duplication event in a closely related species, Caulerpa lentillifera. A total of 100 genes were acquired through horizontal gene transfer, including a homolog of the cryptochrome photoreceptor CRY gene. We also found that all four species studied in Bryopsidales lack key photoprotective genes (LHCSR, PsbS, CYP97A3, and VDE) involved in the xanthophyll cycle and energy-dependent quenching processes. We elucidated that the expansion of light-harvesting antenna genes and the biosynthesis pathways for siphonein and siphonaxanthin in B. corticulans likely contribute to its adaptation to intertidal light conditions. Our study unraveled the underlying special genetic basis of Bryopsis' adaptation to intertidal environments, advancing our understanding of plant adaptive evolution.

Analysis of the Genes from Gibberellin, Jasmonate, and Auxin Signaling Under Drought Stress: A Genome-Wide Approach in Castor Bean (Ricinus communis L.)

Castor bean (Ricinus communis L.) can tolerate long periods of dehydration, allowing the investigation of gene circuits involved in drought tolerance. Genes from gibberellins, jasmonates, and auxin signaling are important for crosstalk in the developmental and environmental adaptation process to drought conditions. However, the genes related to these signals, as well as their transcription profiles under drought, remain poorly characterized in the castor bean. In the present work, genes from gibberellins, jasmonates, and auxin signaling were identified and molecularly characterized. These analyses allowed us to identify genes encoding receptors, inhibitory proteins, and transcription factors from each signaling pathway in the castor bean genome. Chromosomal distribution, gene structure, evolutionary relationships, and conserved motif analyses were performed. Expression analysis through RNA-seq and RT-qPCR revealed that gibberellins, jasmonates, and auxin signaling were modulated at multiple levels under drought, with notable changes in specific genes. The gibberellin receptor RcGID1c was downregulated in response to drought, and RcDELLA3 was strongly repressed, whereas its homologues were not, reinforcing the suggestion of a nuanced regulation of gibberellin signaling during drought. Considering jasmonate signaling, the downregulation of the transcription factor RcMYC2 aligned with the drought tolerance observed in mutants lacking this gene. Altogether, these analyses have provided insights into hormone signaling in the castor bean, unveiling transcriptional responses that enhance our understanding of high drought tolerance in this plant. This knowledge opens avenues for identifying potential candidate genes suitable for genetic manipulation in biotechnological approaches.

Evolution of Plant Conserved microRNAs After Whole-Genome Duplications

MicroRNAs (miRNAs) are a specialized class of small silencing RNAs that regulate gene expression in numerous biological processes in eukaryotes. While the evolutionary dynamics of protein-coding genes after plant whole-genome duplications (WGDs) has been extensively studied, the patterns of evolution for conserved MIRNAs (miRNA genes) post-WGDs are less understood. In this study, we systematically investigated miRNAs and their targets in 6 plant species with varying WGD histories. Our findings reveal that WGDs significantly contribute to the expansion of conserved miRNA families. Notably, through homologous analyses of conserved miRNA families, we discovered that beyond the loci derived from WGDs and other duplication events, some conserved miRNA families have independently gained new loci and/or lost syntenic loci in specific lineage or species through evolution. Additionally, our analyses of sequence divergence in conserved miRNAs showed that the mature sequences of miRNA duplicates gradually diverge following WGDs, with this sequence divergence being correlated with that of their adjacent protein-coding genes after recent WGDs. Furthermore, expression and functional divergence analyses of duplicated targets in different miRNA-target interaction scenarios suggest that conserved miRNAs may play crucial roles in regulating the expression of duplicated genes and related regulatory networks following WGDs. In summary, our analyses reveal universal evolutionary patterns of plant conserved miRNAs following WGDs and provide evidence that some miRNA copies in conserved families originated independently during evolution.

CsMPDB 1.0: An interactive web application for visualizing and exploring the microRNAs and phasiRNAs of tea plant (Camellia sinensis var. sinensis 'Shuchazao')

The microRNAs and phasiRNAs of plant are small non-coding RNAs with important functions through regulating gene expression at the post-transcriptional level. However, identifying miRNAs, phasiRNAs and their target genes from numerous sequencing raw data requires multiple software and command-line operations, which are time-consuming and labor-intensive for non-model plants. Therefore, we present CsMPDB (miRNAs and phasiRNAs database of Camellia sinensis), an interactive web application with multiple analysis modules developed to visualize and explore miRNA and phasiRNA in tea plants based on 259 sRNA-seq samples and 24 degradome-seq samples in NCBI. The source code for the CsMPDB was written in R/shiny. It is compatible, extendable, and portable to be easily set up on different operating systems, and can be accessed at http://myshiny.cpolar.io/CsMPDB. This application plays an important role in accelerating the functional study of sRNAs in the transcriptional regulation of tea plants, and has important reference value for the development of sRNA databases of other species.

Comprehensive genome annotation of the model ciliate Tetrahymena thermophila by in-depth epigenetic and transcriptomic profiling

The ciliate Tetrahymena thermophila is a well-established unicellular model eukaryote, contributing significantly to foundational biological discoveries. Despite its acknowledged importance, current studies on Tetrahymena biology face challenges due to gene annotation inaccuracy, particularly the notable absence of untranslated regions (UTRs). To comprehensively annotate the Tetrahymena macronuclear genome, we collected extensive transcriptomic data spanning various cell stages. To ascertain transcript orientation and transcription start/end sites, we incorporated data on epigenetic marks displaying enrichment towards the 5' end of gene bodies, including H3 lysine 4 tri-methylation (H3K4me3), histone variant H2A.Z, nucleosome positioning and N6-methyldeoxyadenine (6mA). Cap-seq data was subsequently applied to validate the accuracy of identified transcription start sites. Additionally, we integrated Nanopore direct RNA sequencing (DRS), strand-specific RNA sequencing (RNA-seq) and assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) data. Using a newly developed bioinformatic pipeline, coupled with manual curation and experimental validation, our work yielded substantial improvements to the current gene models, including the addition of 2,481 new genes, updates to 23,936 existing genes, and the incorporation of 8,339 alternatively spliced isoforms. Furthermore, novel UTR information was annotated for 26,687 high-confidence genes. Intriguingly, 20% of protein-coding genes were identified to have natural antisense transcripts characterized by high diversity in alternative splicing, thus offering insights into understanding transcriptional regulation. Our work will enhance the utility of Tetrahymena as a robust genetic toolkit for advancing biological research, and provides a promising framework for genome annotation in other eukaryotes.

In silico identification of chilli genome encoded MicroRNAs targeting the 16S rRNA and secA genes of "Candidatus phytoplasma trifolii"

Phytoplasma, a potentially hazardous pathogen associated with witches' broom, is an economically harmful disease-producing bacteria that damages chilli cultivation. Phytoplasma-infected plants display various symptoms that indicate significant disruptions in normal plant physiology and behaviour. Diseases caused by phytoplasma are widespread and have a major economic impact on crop quality and yield. This work focuses on identifying and examining chilli microRNAs (miRNAs) as potential targets against the 16S rRNA and secA gene of "Candidatus Phytoplasma trifolii" ("Ca. P. trifolii") through plant miRNA prediction algorithms. Mature chilli miRNAs (CA-miRNAs) were collected and used to hybridise the 16S rRNA and secA genes. A total of four common CA-miRNAs were picked according to genetic consensus. Three algorithms applied in the present study suggested that the physiologically relevant, top-ranked miR169b_2 has a possibly specific site at nucleotide position 1,006 for targeting the 'Ca. P. trifolii' 16S rRNA gene. The circos algorithm was then utilised to create the miRNA-mRNA regulatory network. The free energy between the miRNA:mRNA duplex was also computed, and the best value of -17.46 kcal/mol was obtained for CA-miR166c_2. Currently, there are no suitable commercial 'Ca. P. trifolii'-resistant chilli crops. As a result, the expected biological data provide useful evidence for developing 'Ca. P. trifolii'-resistant chilli plants.

Almond Grafting for Plum Pox Virus Resistance Triggers Significant Transcriptomic and Epigenetic Shifts in Peaches

Sharka disease, caused by the plum pox virus (PPV), negatively impacts stone fruit production, resulting in economic losses. It has been demonstrated that grafting the almond (Prunus dulcis (Miller) D.A. Webb) variety 'Garrigues' into susceptible peach (Prunus persica (L.) Batsch) rootstocks can result in PPV resistance. The molecular circuits related to grafting in Prunus species, however, have not been fully investigated. In this study, susceptible peach rootstocks 'GF305' were either heterografted with 'Garrigues' almond or homografted with the same cultivar. Peach samples were collected at two stages of scion development, with ungrafted plants utilized as controls. Profiles of transcripts, small RNAs (sRNAs), and DNA methylation were obtained and analyzed on a genome-wide scale. Homografting and heterografting significantly altered the transcriptome and methylome of peach rootstocks, with these modifications being more pronounced during the early stages of scion development. The profiles of sRNAs were significantly more impacted when almonds were used as a scion as opposed to peaches, likely due to the transmission of PPV-unrelated viral sequences. Gene expression differences resulting from DNA methylation alterations are more thoroughly documented at the promoter sequences of genes than within their bodies. This study suggests that the 'Garrigues' almond variety triggers a complex defense response in the peach rootstock, potentially involving the interplay of epigenetic modifications and small RNA-mediated priming of antiviral defenses, which ultimately may contribute to PPV resistance.

Genome-Wide Identification and Evolutionary Analysis of Functional BBM-like Genes in Plant Species

Background/Objectives: BABY BOOM (BBM), a transcription factor from the APETALA2 (AP2) protein family, plays a critical role in somatic embryo induction and apomixis. BBM has now been widely applied to induce apomixis or enhance plant transformation and regeneration efficiency through overexpression or ectopic expression. However, the structural and functional evolutionary history of BBM genes in plants is still not well understood. Methods: The protein sequences of 10 selected plant species were used to locate the branch of BBM-Like by key domain identification and phylogenetic tree construction. The identified BBML genes were used for further conserved motif identification, gene structural analysis, miRNA binding site prediction, cis-acting element prediction, collinear analysis, protein-protein interaction network construction, three-dimensional structure modeling, molecular docking, and expression pattern analysis. Results: A total of 24 BBML proteins were identified from 10 representative plant species. Phylogenetic relationship analysis displayed that BBML proteins from eudicots and monocots were divided into two clusters, with monocots exhibiting a higher number of BBMLs. Gene duplication events indicated that whole genome/segmental duplication were the primary drivers of BBML genes' evolution in the tested species, with purifying selection playing a key role during evolution processes. Comparative analysis of motif, domains, and gene structures revealed that most BBMLs were highly evolutionarily conserved. The expression patterns of BBML genes revealed significant tissue specificity, particularly in the root and embryo. We also constructed protein-protein interaction networks and molecular docking models to identify functional pathways and key amino acid residues of BBML proteins. The functions of BBMLs may differ between monocots and eudicots, as suggested by the functional enrichment of interacting proteins. Conclusions: Our research delved into the molecular mechanism, evolutionary relationships, functional differentiation, and expression patterns of BBML genes across plants, laying the groundwork for further investigations into the molecular properties and biological roles of BBMLs.

Genome-wide analysis of miR172-mediated response to heavy metal stress in chickpea (Cicer arietinum L.): physiological, biochemical, and molecular insights

BACKGROUND

Chickpea (Cicer arietinum L.), a critical diploid legume in the Fabaceae family, is a rich source of protein, vitamins, and minerals. However, heavy metal toxicity severely affects its growth, yield, and quality. MicroRNAs (miRNAs) play a crucial role in regulating plant responses to both abiotic and biotic stress, including heavy metal exposure, by suppressing the expression of target genes. Plants respond to heavy metal stress through miRNA-mediated regulatory mechanisms at multiple physiological, biochemical, and molecular levels. Although the Fabaceae family is well represented in miRNA studies, chickpeas have been notably underrepresented. This study aimed to investigate the effects of heavy metal-induced stress, particularly from 100 µM concentrations of cadmium (Cd), chromium (Cr), nickel (Ni), lead (Pb), and 30 µM arsenic (As), on two chickpea varieties: ILC 482 (sensitive) and Azkan (tolerant). The assessment focused on physiological, biochemical, and molecular parameters. Furthermore, a systematic characterization of the miR172 gene family in the chickpea genome was conducted to better understand the plant's molecular response to heavy metal stress.

RESULTS

Variance analysis indicated significant effects of genotype (G), treatment (T), and genotype-by-treatment (GxT) interactions on plant growth, physiological, and biochemical parameters. Heavy metal stress negatively impacted plant growth in chickpea genotypes ILC 482 and Azkan. A reduction in chlorophyll content and relative leaf water content was observed, along with increased cell membrane damage. In ILC 482, the highest hydrogen peroxide (H₂O₂) levels in shoot tissue were recorded under As, Cd, and Ni treatments, while in Azkan, peak levels were observed with Pb treatment. Malondialdehyde (MDA) levels in root tissue were highest in ILC 482 under Cd and Ni exposure and in Azkan under As, Cr, and Cd treatments. Antioxidant enzyme activities, including superoxide dismutase (SOD), catalase (CAT), peroxidase (POD), and ascorbate peroxidase (APX), were significantly elevated under heavy metal stress in both genotypes. Gene expression analysis revealed upregulation of essential antioxidant enzyme genes, such as SOD, CAT, and APX, with APX showing notable increases in both shoot and root tissues compared to the control. Additionally, seven miR172 genes (miR172a, miR172b, miR172c, miR172d, miR172e, miR172f, and miR172g) were identified in the chickpea genome, distributed across five chromosomes. All genes exhibited conserved hairpin structures essential for miRNA functionality. Phylogenetic analysis grouped these miR172 genes into three clades, suggesting strong evolutionary conservation with other plant species. The expression analysis of miR172 and its target genes under heavy metal stress showed varied expression patterns, indicating their role in enhancing heavy metal tolerance in chickpea.

CONCLUSIONS

Heavy metal stress significantly impaired plant growth and physiological and biochemical parameters in chickpea genotypes, except for cell membrane damage. The findings underscore the critical role of miR172 and its target genes in modulating chickpea's response to heavy metal stress. These insights provide a foundational understanding for developing stress-tolerant chickpea varieties through miRNA-based genetic engineering approaches.

Involvement of MicroRNAs in the Hypersensitive Response of Capsicum Plants to the Capsicum Chlorosis Virus at Elevated Temperatures

The orthotospovirus capsicum chlorosis virus (CaCV) is an important pathogen affecting capsicum plants. Elevated temperatures may affect disease progression and pose a potential challenge to capsicum production. To date, CaCV-resistant capsicum breeding lines have been established; however, the impact of an elevated temperature of 35 °C on this genetic resistance remains unexplored. Thus, this study aimed to investigate how high temperature (HT) influences the response of CaCV-resistant capsicum to the virus. Phenotypic analysis revealed a compromised resistance in capsicum plants grown at HT, with systemic necrotic spots appearing in 8 out of 14 CaCV-infected plants. Molecular analysis through next-generation sequencing identified 105 known and 83 novel microRNAs (miRNAs) in CaCV-resistant capsicum plants. Gene ontology revealed that phenylpropanoid and lignin metabolic processes, regulated by Can-miR408a and Can- miR397, are likely involved in elevated-temperature-mediated resistance-breaking responses. Additionally, real-time PCR validated an upregulation of Can-miR408a and Can-miR397 by CaCV infection at HT; however, only the Laccase 4 transcript, targeted by Can-miR397, showed a tendency of negative correlation with this miRNA. Overall, this study provides the first molecular insights into how elevated temperature affects CaCV resistance in capsicum plants and reveals the potential role of miRNA in temperature-sensitive tospovirus resistance.

SynGAP: a synteny-based toolkit for gene structure annotation polishing

Genome sequencing has become a routine task for biologists, but the challenge of gene structure annotation persists, impeding accurate genomic and genetic research. Here, we present a bioinformatics toolkit, SynGAP (Synteny-based Gene structure Annotation Polisher), which uses gene synteny information to accomplish precise and automated polishing of gene structure annotation of genomes. SynGAP offers exceptional capabilities in the improvement of gene structure annotation quality and the profiling of integrative gene synteny between species. Furthermore, an expression variation index is designed for comparative transcriptomics analysis to explore candidate genes responsible for the development of distinct traits observed in phylogenetically related species.

Allopolyploidization from two dioecious ancestors leads to recurrent evolution of sex chromosomes

Polyploidization presents an unusual challenge for species with sex chromosomes, as it can lead to complex combinations of sex chromosomes that disrupt reproductive development. This is particularly true for allopolyploidization between species with different sex chromosome systems. Here, we assemble haplotype-resolved chromosome-level genomes of a female allotetraploid weeping willow (Salix babylonica) and a male diploid S. dunnii. We show that weeping willow arose from crosses between a female ancestor from the Salix-clade, which has XY sex chromosomes on chromosome 7, and a male ancestor from the Vetrix-clade, which has ancestral XY sex chromosomes on chromosome 15. We find that weeping willow has one pair of sex chromosomes, ZW on chromosome 15, that derived from the ancestral XY sex chromosomes in the male ancestor of the Vetrix-clade. Moreover, the ancestral 7X chromosomes from the female ancestor of the Salix-clade have reverted to autosomal inheritance. Duplicated intact ARR17-like genes on the four homologous chromosomes 19 likely have contributed to the maintenance of dioecy during polyploidization and sex chromosome turnover. Taken together, our results suggest the rapid evolution and reversion of sex chromosomes following allopolyploidization in weeping willow.

SapBase: A central portal for functional and comparative genomics of Sapindaceae species

The Sapindaceae family, encompassing a wide range of plant forms such as herbs, vines, shrubs, and trees, is widely distributed across tropical and subtropical regions. This family includes economically important crops like litchi, longan, rambutan, and ackee. With the wide application of genomic technologies in recent years, several Sapindaceae plant genomes have been decoded, leading to an accumulation of substantial omics data in this field. This surge in data highlights the pressing need for a unified genomic data center capable of storing, sharing, and analyzing these data. Here, we introduced SapBase, that is, the Sapindaceae Genome Database. SapBase houses seven published plant genomes alongside their corresponding gene structure and functional annotations, small RNA annotations, gene expression profiles, gene pathways, and synteny block information. It offers user-friendly features for gene information mining, co-expression analysis, and inter-species comparative genomic analysis. Furthermore, we showcased SapBase's extensive capacities through a detailed bioinformatic analysis of a MYB gene in litchi. Thus, SapBase could serve as an integrative genomic resource and analysis platform for the scientific exploration of Sapinaceae species and their comparative studies with other plants.

Functional characterization of polyphenol oxidase OfPPO2 supports its involvement in parallel biosynthetic pathways of acteoside

Acteoside is a bioactive phenylethanoid glycoside widely distributed throughout the plant kingdom. Because of its two catechol moieties, acteoside displays a variety of beneficial activities. The biosynthetic pathway of acteoside has been largely elucidated, but the assembly logic of two catechol moieties in acteoside remains unclear. Here, we identified a novel polyphenol oxidase OfPPO2 from Osmanthus fragrans, which could hydroxylate various monophenolic substrates, including tyrosine, tyrosol, tyramine, 4-hydroxyphenylacetaldehyde, salidroside, and osmanthuside A, leading to the formation of corresponding catechol-containing intermediates for acteoside biosynthesis. OfPPO2 could also convert osmanthuside B into acteoside, creating catechol moieties directly via post-modification of the acteoside skeleton. The reverse transcription-quantitative polymerase chain reaction (RT-qPCR) analysis and subcellular localization assay further support the involvement of OfPPO2 in acteoside biosynthesis in planta. These findings suggest that the biosynthesis of acteoside in O. fragrans may follow "parallel routes" rather than the conventionally considered linear route. In support of this hypothesis, the glycosyltransferase OfUGT and the acyltransferase OfAT could direct the flux of diphenolic intermediates generated by OfPPO2 into acteoside. Significantly, OfPPO2 and its orthologs constitute a functionally conserved enzyme family that evolved independently from other known biosynthetic enzymes of acteoside, implying that the substrate promiscuity of this PPO family may offer acteoside-producing plants alternative ways to synthesize acteoside. Overall, this work expands our understanding of parallel pathways plants may employ to efficiently synthesize acteoside, a strategy that may contribute to plants' adaptation to environmental challenges.

Identification of Arabidopsis thaliana small RNAs responsive to the fungal pathogen Botrytis cinerea at an early stage of interaction

In plants, small RNAs (sRNAs), mainly microRNAs (miRNAs) and small interfering RNAs (siRNAs), have been described as key regulators of plant development, growth, and abiotic and biotic responses. Despite reports indicating the involvement of certain sRNAs in regulating the interaction between Botrytis cinerea (a major necrotrophic fungal phytopathogen) and host plants, there remains a lack of analysis regarding the potential regulatory roles of plant sRNAs during early stages of the interaction despite early immune responses observed then during infection. We present the first transcriptome-wide analysis of small RNA expression on the early interaction between the necrotrophic fungus Botrytis cinerea and the model plant Arabidopsis thaliana. We found that evolutionary conserved A. thaliana miRNAs were the sRNAs that accumulated the most in the presence of B. cinerea. The upregulation of miR167, miR159 and miR319 was of particular interest because these, together with their target transcripts, are involved in the fine regulation of the plant hormone signaling pathways. We also describe that miR173, which triggers the production of secondary siRNAs from TAS1 and TAS2 loci, as well as secondary siRNAs derived from these loci, is upregulated in response to B. cinerea. Thus, at an early stage of the interaction there are transcriptional changes of sRNA-guided silencing pathway genes and of a subset of sRNAs that targeted genes from the PPR gene superfamily, and these may be important mechanisms regulating the interaction between A. thaliana and B. cinerea. This work provides the basis for a better understanding of the regulation mediated by sRNAs during early B. cinerea-plant interaction and may help in the development of more effective strategies for its control.

Multiple strategies, including 6mA methylation, affecting plant alternative splicing in allopolyploid peanut

Alternative splicing (AS), an important post-transcriptional regulation mechanism in eukaryotes, can significantly increase transcript diversity and contribute to gene expression regulation and many other complicated developmental processes. While plant gene AS events are well described, few studies have investigated the comprehensive regulation machinery of plant AS. Here, we use multi-omics to analyse peanut AS events. Using long-read isoform sequencing, 146 464 full-length non-chimeric transcripts were obtained, resulting in annotation corrections for 1782 genes and the identification of 4653 new loci. Using Iso-Seq RNA sequences, 271 776 unique splice junctions were identified, 82.49% of which were supported by transcriptome data. We characterized 50 977 polyadenylation sites for 23 262 genes, 12 369 of which had alternative polyadenylation sites. AS allows differential regulation of the same gene by miRNAs at the isoform level coupled with polyadenylation. In addition, we identified many long non-coding RNAs and fusion transcripts. There is a suppressed effect of 6mA on AS and gene expression. By analysis of chromatin structures, the genes located in the boundaries of topologically associated domains, proximal chromosomal telomere regions, inter- or intra-chromosomal loops were found to have more unique splice isoforms, higher expression, lower 6mA and more transposable elements (TEs) in their gene bodies than the other genes, indicating that chromatin interaction, 6mA and TEs play important roles in AS and gene expression. These results greatly refine the peanut genome annotation and contribute to the study of gene expression and regulation in peanuts. This work also showed AS is associated with multiple strategies for gene regulation.

A novel ilarvirus protein CP-RT is expressed via stop codon readthrough and suppresses RDR6-dependent RNA silencing

Ilarviruses are a relatively understudied but important group of plant RNA viruses that includes a number of crop pathogens. Their genomes comprise three RNA segments encoding two replicase subunits, movement protein, coat protein (CP), and (in some ilarvirus subgroups) a protein that suppresses RNA silencing. Here we report that, in many ilarviruses, RNA3 encodes an additional protein (termed CP-RT) as a result of ribosomal readthrough of the CP stop codon into a short downstream readthrough (RT) ORF. Using asparagus virus 2 as a model, we find that CP-RT is expressed in planta where it functions as a weak suppressor of RNA silencing. CP-RT expression is essential for persistent systemic infection in leaves and shoot apical meristem. CP-RT function is dependent on a putative zinc-finger motif within RT. Replacing the asparagus virus 2 RT with the RT of an ilarvirus from a different subgroup restored the ability to establish persistent infection. These findings open up a new avenue for research on ilarvirus silencing suppression, persistent meristem invasion and vertical transmission.

The miR166 targets CsHDZ3 genes to negatively regulate drought tolerance in tea plant (Camellia sinensis)

Drought is the stressor with a significant adverse impact on the yield stability of tea plants. HD-ZIP III transcription factors (TFs) play important regulatory roles in plant growth, development, and stress responses. However, whether and how HD-ZIP III TFs are involved in drought response and tolerance in tea plants remains unclear. Here, we identified seven HD-ZIP III genes (CsHDZ3-1 to CsHDZ3-7) in tea plant genome. The evolutionary analysis demonstrated that CsHDZ3 members were subjected to purify selection. Subcellular localization analysis revealed that all seven CsHDZ3s located in the nucleus. Yeast self-activation and dual-luciferase reporter assays demonstrated that CsHDZ3-1 to CsHDZ3-4 have trans-activation ability whereas CsHDZ3-5 to CsHDZ3-7 served as transcriptional inhibitors. The qRT-PCR assay showed that all seven CsHDZ3 genes could respond to simulated natural drought stress and polyethylene glycol treatment. Further assays verified that all CsHDZ3 genes can be cleaved by csn-miR166. Overexpression of csn-miR166 inhibited the expression of seven CsHDZ3 genes and weakened drought tolerance of tea leaves. In contrast, suppression of csn-miR166 promoted the expression of seven CsHDZ3 genes and enhanced drought tolerance of tea leaves. These findings established the foundation for further understanding the mechanism of CsHDZ3-miR166 modules' participation in drought responses and tolerance.

MicroRNA482/2118 is lineage-specifically involved in gibberellin signalling via the regulation of GID1 expression by targeting noncoding PHAS genes and subsequently instigated phasiRNAs

MicroRNA482/2118 (miR482/2118) is a 22-nt miRNA superfamily, with conserved functions in disease resistance and plant development. It usually instigates the production of phased small interfering RNAs (phasiRNAs) from its targets to expand or reinforce its silencing effect. Using a new high-quality reference genome sequence and comprehensive small RNA profiling, we characterized a newly evolved regulatory pathway of miR482/2118 in litchi. In this pathway, miR482/2118 cleaved a novel noncoding trans-acting gene (LcTASL1) and triggered phasiRNAs to regulate the expression of gibberellin (GA) receptor gene GIBBERELLIN INSENSITIVE DWARF1 (GID1) in trans; another trans-acting gene LcTASL2, targeted by LcTASL1-derived phasiRNAs, produced phasiRNAs as well to target LcGID1 to reinforce the silencing effect of LcTASL1. We found this miR482/2118-TASL-GID1 pathway was likely involved in fruit development, especially the seed development in litchi. In vivo construction of the miR482a-TASL-GID1 pathway in Arabidopsis could lead to defects in flower and silique development, analogous to the phenotype of gid1 mutants. Finally, we found that a GA-responsive transcription factor, LcGAMYB33, could regulate LcMIR482/2118 as a feedback mechanism of the sRNA-silencing pathway. Our results deciphered a lineage-specifically evolved regulatory module of miR482/2118, demonstrating the high dynamics of miR482/2118 function in plants.

sRNAminer: A multifunctional toolkit for next-generation sequencing small RNA data mining in plants

Small RNAs (sRNAs), found extensively in plants, play an essential role in plant growth and development. Although various sRNA analysis tools have been developed for plants, the use of most of them depends on programming and command-line environments, which is a challenge for many wet-lab biologists. Furthermore, current sRNA analysis tools mostly focus on the analysis of certain type of sRNAs and are resource-intensive, normally demanding an immense amount of time and effort to learn the use of numerous tools or scripts and assemble them into a workable pipeline to get the final results. Here, we present sRNAminer, a powerful stand-alone toolkit with a user-friendly interface that integrates all common functions for the analysis of three major types of plant sRNAs: microRNAs (miRNAs), phased small interfering RNAs (phasiRNAs), and heterochromatic siRNAs (hc-siRNAs). We constructed a curated or "golden" set of MIRNA and PHAS loci, which was used to assess the performance of sRNAminer in comparison to other existing tools. The results showed that sRNAminer outperformed these tools in multiple aspects, highlighting its functionality. In addition, to enable an efficient evaluation of sRNA annotation results, we developed Integrative Genomics Viewer (IGV)-sRNA, a modified genome browser optimized from IGV and we incorporated it as a functional module in sRNAminer. IGV-sRNA can display a wealth of sRNA-specific features, enabling a more comprehensive understanding of sRNA data. sRNAminer and IGV-sRNA are both platform-independent software that can be run under all operating systems. They are now freely available at https://github.com/kli28/sRNAminer and https://gitee.com/CJchen/IGV-sRNA.

Plant non-coding RNAs: The new frontier for the regulation of plant development and adaptation to stress

Most plant transcriptomes constitute functional non-coding RNAs (ncRNAs) that lack the ability to encode proteins. In recent years, more research has demonstrated that ncRNAs play important regulatory roles in almost all plant biological processes by modulating gene expression. Thus, it is important to study the biogenesis and function of ncRNAs, particularly in plant growth and development and stress tolerance. In this review, we systematically explore the process of formation and regulatory mechanisms of ncRNAs, particularly those of microRNAs (miRNAs), small interfering RNAs (siRNAs), long non-coding RNAs (lncRNAs), and circular RNAs (circRNAs). Additionally, we provide a comprehensive overview of the recent advancements in ncRNAs research, including their regulation of plant growth and development (seed germination, root growth, leaf morphogenesis, floral development, and fruit and seed development) and responses to abiotic and biotic stress (drought, heat, cold, salinity, pathogens and insects). We also discuss research challenges and provide recommendations to advance the understanding of the roles of ncRNAs in agronomic applications.

Global hypermethylation of the N6-methyladenosine RNA modification associated with apple heterografting

Grafting can facilitate better scion performance and is widely used in plants. Numerous studies have studied the involvement of mRNAs, small RNAs, and epigenetic regulations in the grafting process. However, it remains unclear whether the mRNA N6-methyladenosine (m6A) modification participates in the apple (Malus x domestica Borkh.) grafting process. Here, we decoded the landscape of m6A modification profiles in 'Golden delicious' (a cultivar, Gd) and Malus prunifolia 'Fupingqiuzi' (a unique rootstock with resistance to environmental stresses, Mp), as well as their heterografted and self-grafted plants. Interestingly, global hypermethylation of m6A occurred in both heterografted scion and rootstock compared with their self-grafting controls. Gene Ontology (GO) term enrichment analysis showed that grafting-induced differentially m6A-modified genes were mainly involved in RNA processing, epigenetic regulation, stress response, and development. Differentially m6A-modified genes harboring expression alterations were mainly involved in various stress responses and fatty acid metabolism. Furthermore, grafting-induced mobile mRNAs with m6A and gene expression alterations mainly participated in ABA synthesis and transport (e.g. carotenoid cleavage dioxygenase 1 [CCD1] and ATP-binding cassette G22 [ABCG22]) and abiotic and biotic stress responses, which might contribute to the better performance of heterografted plants. Additionally, the DNA methylome analysis also demonstrated the DNA methylation alterations during grafting. Downregulated expression of m6A methyltransferase gene MdMTA (ortholog of METTL3) in apples induced the global m6A hypomethylation and distinctly activated the expression level of DNA demethylase gene MdROS1 (REPRESSOR OF SILENCING 1) showing the possible association between m6A and 5mC methylation in apples. Our results reveal the m6A modification profiles in the apple grafting process and enhance our understanding of the m6A regulatory mechanism in plant biological processes.

Role of miRNAs in sucrose stress response, reactive oxygen species, and anthocyanin biosynthesis in Arabidopsis thaliana

In plants, sucrose is the main transported disaccharide that is the primary product of photosynthesis and controls a multitude of aspects of the plant life cycle including structure, growth, development, and stress response. Sucrose is a signaling molecule facilitating various stress adaptations by crosstalk with other hormones, but the molecular mechanisms are not well understood. Accumulation of high sucrose concentrations is a hallmark of many abiotic and biotic stresses, resulting in the accumulation of reactive oxygen species and secondary metabolite anthocyanins that have antioxidant properties. Previous studies have shown that several MYeloBlastosis family/MYB transcription factors are positive and negative regulators of sucrose-induced anthocyanin accumulation and subject to microRNA (miRNA)-mediated post-transcriptional silencing, consistent with the notion that miRNAs may be "nodes" in crosstalk signaling by virtue of their sequence-guided targeting of different homologous family members. In this study, we endeavored to uncover by deep sequencing small RNA and mRNA transcriptomes the effects of exogenous high sucrose stress on miRNA abundances and their validated target transcripts in Arabidopsis. We focused on genotype-by-treatment effects of high sucrose stress in Production of Anthocyanin Pigment 1-Dominant/pap1-D, an activation-tagged dominant allele of MYB75 transcription factor, a positive effector of secondary metabolite anthocyanin pathway. In the process, we discovered links to reactive oxygen species signaling through miR158/161/173-targeted Pentatrico Peptide Repeat genes and two novel non-canonical targets of high sucrose-induced miR408 and miR398b*(star), relevant to carbon metabolic fluxes: Flavonoid 3'-Hydroxlase (F3'H), an important enzyme in determining the B-ring hydroxylation pattern of flavonoids, and ORANGE a post-translational regulator of Phytoene Synthase expression, respectively. Taken together, our results contribute to understanding the molecular mechanisms of carbon flux shifts from primary to secondary metabolites in response to high sugar stress.

Manual correction of genome annotation improved alternative splicing identification of Artemisia annua

Gene annotation is essential for genome-based studies. However, algorithm-based genome annotation is difficult to fully and correctly reveal genomic information, especially for species with complex genomes. Artemisia annua L. is the only commercial resource of artemisinin production though the content of artemisinin is still to be improved. Genome-based genetic modification and breeding are useful strategies to boost artemisinin content and therefore, ensure the supply of artemisinin and reduce costs, but better gene annotation is urgently needed. In this study, we manually corrected the newly released genome annotation of A. annua using second- and third-generation transcriptome data. We found that incorrect gene information may lead to differences in structural, functional, and expression levels compared to the original expectations. We also identified alternative splicing events and found that genome annotation information impacted identifying alternative splicing genes. We further demonstrated that genome annotation information and alternative splicing could affect gene expression estimation and gene function prediction. Finally, we provided a valuable version of A. annua genome annotation and demonstrated the importance of gene annotation in future research.

Overexpression of Sly-miR398b Compromises Disease Resistance against Botrytis cinerea through Regulating ROS Homeostasis and JA-Related Defense Genes in Tomato

MicroRNAs (miRNAs) have been shown to be critical components in plant immunity. MicroRNA398 (miR398) is a highly conserved miRNA in all land plants and plays crucial roles in diverse biotic stress responses. However, the role of miR398 has not yet been characterized in tomato resistance against Botrytis cinerea. In this report, the transcript levels of sly-miR398b were strongly decreased in B. cinerea-infected leaves and the overexpression of sly-miR398b resulted in enhanced susceptibility. The attenuated expression of cytosol Cu/Zn-SOD (CSD1), chloroplast Cu/Zn-SOD (CSD2), and guaiacol peroxidase (GPOD), as well as the decreased activities of superoxide dismutase (SOD) and GPOD, collectively led to increased hydrogen peroxide (H₂O₂) accumulation in sly-miR398b overexpressing plants. Furthermore, sly-miR398b was induced by methyl jasmonate (MeJA) treatment. The overexpression of sly-miR398b suppressed the expression of TomLoxD, LapA, and PR-STH2 in response to B. cinerea and MeJA treatment. Our data demonstrate that sly-miR398b overexpression negatively regulates the resistance to B. cinerea in tomato by inducing the accumulation of reactive oxygen species (ROS) and downregulating the expression of MeJA-responsive defense genes.

Integrative time-scale and multi-omics analysis of host responses to viroid infection

Viroids are circular RNAs of minimal complexity compelled to subvert plant-regulatory networks to accomplish their infectious process. Studies focused on the response to viroid-infection have mostly addressed specific regulatory levels and considered specifics infection-times. Thus, much remains to be done to understand the temporal evolution and complex nature of viroid-host interactions. Here we present an integrative analysis of the temporal evolution of the genome-wide alterations in cucumber plants infected with hop stunt viroid (HSVd) by integrating differential host transcriptome, sRNAnome and methylome. Our results support that HSVd promotes the redesign of the cucumber regulatory-pathways predominantly affecting specific regulatory layers at different infection-phases. The initial response was characterised by a reconfiguration of the host-transcriptome by differential exon-usage, followed by a progressive transcriptional downregulation modulated by epigenetic changes. Regarding endogenous small RNAs, the alterations were limited and mainly occurred at the late stage. Significant host-alterations were predominantly related to the downregulation of transcripts involved in plant-defence mechanisms, the restriction of pathogen-movement and the systemic spreading of defence signals. We expect that these data constituting the first comprehensive temporal-map of the plant-regulatory alterations associated with HSVd infection could contribute to elucidate the molecular basis of the yet poorly known host-response to viroid-induced pathogenesis.

Small RNAs >26 nt in length associate with AGO1 and are upregulated by nutrient deprivation in the alga Chlamydomonas

Small RNAs (sRNAs) associate with ARGONAUTE (AGO) proteins forming effector complexes with key roles in gene regulation and defense responses against molecular parasites. In multicellular eukaryotes, extensive duplication and diversification of RNA interference (RNAi) components have resulted in intricate pathways for epigenetic control of gene expression. The unicellular alga Chlamydomonas reinhardtii also has a complex RNAi machinery, including 3 AGOs and 3 DICER-like proteins. However, little is known about the biogenesis and function of most endogenous sRNAs. We demonstrate here that Chlamydomonas contains uncommonly long (>26 nt) sRNAs that associate preferentially with AGO1. Somewhat reminiscent of animal PIWI-interacting RNAs, these >26 nt sRNAs are derived from moderately repetitive genomic clusters and their biogenesis is DICER-independent. Interestingly, the sequences generating these >26-nt sRNAs have been conserved and amplified in several Chlamydomonas species. Moreover, expression of these longer sRNAs increases substantially under nitrogen or sulfur deprivation, concurrently with the downregulation of predicted target transcripts. We hypothesize that the transposon-like sequences from which >26-nt sRNAs are produced might have been ancestrally targeted for silencing by the RNAi machinery but, during evolution, certain sRNAs might have fortuitously acquired endogenous target genes and become integrated into gene regulatory networks.

Delineating the tissue-mediated drought stress governed tuning of conserved miR408 and its targets in rice

Engineering drought tolerance in rice needs to focus on regulators that enhance tolerance while boosting plant growth and vigor. The present study delineated the concealed function and tissue-mediated interplay of the miR408/target module in imparting drought stress tolerance in rice. The plant miR408 family comprises three dominant mature forms (21 nt), including a distinct monocot variant (F-7 with 5' C) and is divided into six groups. miR408 majorly cleaves genes belonging to the blue copper protein in addition to several other species-specific targets in plants. Comparative sequence analysis in 4726 rice accessions identified 22 sequence variants (SNP and InDELs) in its promoter (15) and pre-miR408 region. Haplotype analysis of the sequence variants indicated eight haplotypes (three: Japonica-specific and five: Indica-specific) of the miR408 promoter. In drought-tolerant Nagina 22, miR408 follows flag leaf preferential expression. Under drought conditions, its levels are upregulated in flag leaf and roots which seems to be regulated by a differential fraction of methylated cytosines (mCs) in the precursor region. The active pool of miR408 regulated targets under control and drought conditions is impacted by the tissue type. Comparative expression analysis of the miR408/target module under different sets of conditions features 83 targets exhibiting antagonistic expression in rice, out of which 12 genes, including four PLANTACYANINS (OsUCL6, 7, 9 and 30), PIRIN, OsLPR1, OsCHUP1, OsDOF12, OsBGLU1, glycine-rich cell wall gene, OsDUT, and OsERF7, are among the high confidence targets. Further, overexpression of MIR408 in drought-sensitive rice cultivar (PB1) leads to the massive enhancement of vegetative growth in rice with improved ETR and Y(II) and enhanced dehydration stress tolerance. The above results suggest that miR408 is likely to act as a positive regulator of growth and vigor, as well as dehydration stress, making it a potential candidate for engineering drought tolerance in rice.

Long non-coding RNA and microRNA landscape of two major domesticated cotton species

Allotetraploid cotton plants Gossypium hirsutum and Gossypium barbadense have been widely cultivated for their natural, renewable textile fibres. Even though ncRNAs in domesticated cotton species have been extensively studied, systematic identification and annotation of lncRNAs and miRNAs expressed in various tissues and developmental stages under various biological contexts are limited. This influences the comprehension of their functions and future research on these cotton species. Here, we report high confidence lncRNAs and miRNA collection from G. hirsutum accession and G. barbadense accession using large-scale RNA-seq and small RNA-seq datasets incorporated into a user-friendly database, CoNCRAtlas. This database provides a wide range and depth of lncRNA and miRNA annotation based on the systematic integration of extensive annotations such as expression patterns derived from transcriptome data analysis in thousands of samples, as well as multi-omics annotations. We assume this comprehensive resource will accelerate evolutionary and functional studies in ncRNAs and inform future breeding programs for cotton improvement. CoNCRAtlas is accessible at http://www.nipgr.ac.in/CoNCRAtlas/.

Small RNAs and Karma methylation in Elaeis guineensis mother palms are linked to high clonal mantling

The mantled phenotype is an abnormal somaclonal variant arising from the oil palm cloning process and severe phenotypes lead to oil yield losses. Hypomethylation of the Karma retrotransposon within the B-type MADS-box EgDEF1 gene has been associated with this phenotype. While abnormal Karma-EgDEF1 hypomethylation was detected in mantled clones, we examined the methylation state of Karma in ortets that gave rise to high mantling rates in their clones. Small RNAs (sRNAs) were proposed to play a role in Karma hypomethylation as part of the RNA-directed DNA methylation process, hence differential expression analysis of sRNAs between the ortet groups was conducted. While no sRNA was differentially expressed at the Karma-EgDEF1 region, three sRNA clusters were differentially regulated in high-mantling ortets. The first two down-regulated clusters were possibly derived from long non-coding RNAs while the third up-regulated cluster was derived from the intron of a DnaJ chaperone gene. Several predicted mRNA targets for the first two sRNA clusters conversely displayed increased expression in high-mantling relative to low-mantling ortets. These predicted mRNA targets may be associated with defense or pathogenesis response. In addition, several differentially methylated regions (DMRs) were identified in Karma and its surrounding regions, mainly comprising subtle CHH hypomethylation in high-mantling ortets. Four of the 12 DMRs were located in a region corresponding to hypomethylated areas at the 3'end of Karma previously reported in mantled clones. Further investigations on these sRNAs and DMRs may indicate the predisposition of certain ortets towards mantled somaclonal variation.

Taxon-specific, phased siRNAs underlie a speciation locus in monkeyflowers

Taxon-specific small RNA loci are widespread in eukaryotic genomes, yet their role in lineage-specific adaptation, phenotypic diversification, and speciation is poorly understood. Here, we report that a speciation locus in monkeyflowers (Mimulus), YELLOW UPPER (YUP), contains an inverted repeat region that produces small interfering RNAs (siRNAs) in a phased pattern. Although the inverted repeat is derived from a partial duplication of a protein-coding gene that is not involved in flower pigmentation, one of the siRNAs targets and represses a master regulator of floral carotenoid pigmentation. YUP emerged with two protein-coding genes that control other aspects of flower coloration as a "superlocus" in a subclade of Mimulus and has contributed to subsequent phenotypic diversification and pollinator-mediated speciation in the descendant species.

An efficient deep learning based predictor for identifying miRNA-triggered phasiRNA loci in plant

Phasic small interfering RNAs are plant secondary small interference RNAs that typically generated by the convergence of miRNAs and polyadenylated mRNAs. A growing number of studies have shown that miRNA-initiated phasiRNA plays crucial roles in regulating plant growth and stress responses. Experimental verification of miRNA-initiated phasiRNA loci may take considerable time, energy and labor. Therefore, computational methods capable of processing high throughput data have been proposed one by one. In this work, we proposed a predictor (DIGITAL) for identifying miRNA-initiated phasiRNAs in plant, which combined a multi-scale residual network with a bi-directional long-short term memory network. The negative dataset was constructed based on positive data, through replacing 60% of nucleotides randomly in each positive sample. Our predictor achieved the accuracy of 98.48% and 94.02% respectively on two independent test datasets with different sequence length. These independent testing results indicate the effectiveness of our model. Furthermore, DIGITAL is of robustness and generalization ability, and thus can be easily extended and applied for miRNA target recognition of other species. We provide the source code of DIGITAL, which is freely available at https://github.com/yuanyuanbu/DIGITAL.

A wheat integrative regulatory network from large-scale complementary functional datasets enables trait-associated gene discovery for crop improvement

Gene regulation is central to all aspects of organism growth, and understanding it using large-scale functional datasets can provide a whole view of biological processes controlling complex phenotypic traits in crops. However, the connection between massive functional datasets and trait-associated gene discovery for crop improvement is still lacking. In this study, we constructed a wheat integrative gene regulatory network (wGRN) by combining an updated genome annotation and diverse complementary functional datasets, including gene expression, sequence motif, transcription factor (TF) binding, chromatin accessibility, and evolutionarily conserved regulation. wGRN contains 7.2 million genome-wide interactions covering 5947 TFs and 127 439 target genes, which were further verified using known regulatory relationships, condition-specific expression, gene functional information, and experiments. We used wGRN to assign genome-wide genes to 3891 specific biological pathways and accurately prioritize candidate genes associated with complex phenotypic traits in genome-wide association studies. In addition, wGRN was used to enhance the interpretation of a spike temporal transcriptome dataset to construct high-resolution networks. We further unveiled novel regulators that enhance the power of spike phenotypic trait prediction using machine learning and contribute to the spike phenotypic differences among modern wheat accessions. Finally, we developed an interactive webserver, wGRN (http://wheat.cau.edu.cn/wGRN), for the community to explore gene regulation and discover trait-associated genes. Collectively, this community resource establishes the foundation for using large-scale functional datasets to guide trait-associated gene discovery for crop improvement.

Mobile ARGONAUTE 1d binds 22-nt miRNAs to generate phasiRNAs important for low-temperature male fertility in rice

Phased small interfering RNAs (phasiRNAs) are abundantly expressed in anthers and linked to environment-related male fertility in grasses, yet how they function under different environmental conditions remains unclear. Here, we identified a rice (Oryza sativa) low temperature-induced Argonaute (AGO) protein, OsAGO1d, that is responsible for generating phasiRNAs and preserving male fertility at low temperature. Loss of OsAGO1d function causes low-temperature male sterility associated with delayed programmed cell death of tapetal cells during anther development. OsAGO1d binds miR2118 and miR2275 family members and triggers phasiRNA biogenesis; it also binds 21-nt phasiRNAs with a 5' terminal U. In total, phasiRNAs from 972 loci are OsAGO1d-dependent. OsAGO1d protein moves from anther wall cells into meiocytes, where it loads miR2275 to produce 24-nt phasiRNAs. Together, our results show that OsAGO1d acts as a mobile signal to fine-tune phasiRNA production and this function is important for male fertility at low temperature.

SimiR396d targets SiGRF1 to regulate drought tolerance and root growth in foxtail millet

MicroRNAs play critical roles in growth, development and abiotic stress responses. SimR396d is a miRNA whose expression level is much higher in foxtail millet roots than other tissues. Whether SimR396d is involved in foxtail millet root growth and response to abiotic stress is still unknown. Here, we demonstrate that SimiR396d modulates both drought response and root growth in foxtail millet. The expression of SimiR396d is induced by PEG treatment. Overexpression of SimiR396d enhances drought tolerance and root length, while knockdown SimiR396d expression using target mimics of SimiR396d (MIM396) resulted in reduced drought tolerance and shortened root length. Furthermore, we identified and confirmed a plant-specific transcription factor, growth-regulating factor 1 (SiGRF1), as a direct target of SimiR396d. Overexpression of SiGRF1 in foxtail millet resulted in suppressed root growth and reduced sensitivity to drought stress. Moreover, ethylene signaling is necessary for SimiR396d and SiGRF1 to participate in the regulation of plant root growth. These results revealed a pivotal role of SimiR396d in drought tolerance and root growth in foxtail millet. SimiR396d-SiGRF1 regulatory module provides a strategy to improve drought-stress resistance of crop.

The roles of microRNAs in horticultural plant disease resistance

The development of the horticultural industry is largely limited by disease and excessive pesticide application. MicroRNAs constitute a major portion of the transcriptomes of eukaryotes. Various microRNAs have been recognized as important regulators of the expression of genes involved in essential biological processes throughout the whole life cycle of plants. Recently, small RNA sequencing has been applied to study gene regulation in horticultural plants. In this review, we summarize the current understanding of the biogenesis and contributions of microRNAs in horticultural plant disease resistance. These microRNAs may potentially be used as genetic resources for improving disease resistance and for molecular breeding. The challenges in understanding horticultural plant microRNA biology and the possibilities to make better use of these horticultural plant gene resources in the future are discussed in this review.

Ascorbate-Glutathione Cycle Genes Families in Euphorbiaceae: Characterization and Evolutionary Analysis

Ascorbate peroxidase (APX), Monodehydroascorbate Reductase (MDAR), Dehydroascorbate Reductase (DHAR) and Glutathione Reductase (GR) enzymes participate in the ascorbate-glutathione cycle, which exerts a central role in the antioxidant metabolism in plants. Despite the importance of this antioxidant system in different signal transduction networks related to development and response to environmental stresses, the pathway has not yet been comprehensively characterized in many crop plants. Among different eudicotyledons, the Euphorbiaceae family is particularly diverse with some species highly tolerant to drought. Here the APX, MDAR, DHAR, and GR genes in Ricinus communis, Jatropha curcas, Manihot esculenta, and Hevea brasiliensis were identified and characterized. The comprehensive phylogenetic and genomic analyses allowed the classification of the genes into different classes, equivalent to cytosolic, peroxisomal, chloroplastic, and mitochondrial enzymes, and revealed the duplication events that contribute to the expansion of these families within plant genomes. Due to the high drought stress tolerance of Ricinus communis, the expression patterns of ascorbate-glutathione cycle genes in response to drought were also analyzed in leaves and roots, indicating a differential expression during the stress. Altogether, these data contributed to the characterization of the expression pattern and evolutionary analysis of these genes, filling the gap in the proposed functions of core components of the antioxidant mechanism during stress response in an economically relevant group of plants.

Gma-miR408 Enhances Soybean Cyst Nematode Susceptibility by Suppressing Reactive Oxygen Species Accumulation

Soybean cyst nematode (SCN, Heterodera glycine) is a serious damaging disease in soybean worldwide, thus resulting in severe yield losses. MicroRNA408 (miR408) is an ancient and highly conserved miRNA involved in regulating plant growth, development, biotic and abiotic stress response. Here, we analyzed the evolution of miR408 in plants and verified four miR408 members in Glycine max. In the current research, highly upregulated gma-miR408 expressing was detected during nematode migration and syncytium formation response to soybean cyst nematode infection. Overexpressing and silencing miR408 vectors were transformed to soybean to confirm its potential role in plant and nematode interaction. Significant variations were observed in the MAPK signaling pathway with low OXI1, PR1, and wounding of the overexpressing lines. Overexpressing miR408 could negatively regulate soybean resistance to SCN by suppressing reactive oxygen species accumulation. Conversely, silencing miR408 positively regulates soybean resistance to SCN. Overall, gma-miR408 enhances soybean cyst nematode susceptibility by suppressing reactive oxygen species accumulation.

Multiple Functions of MiRNAs in Brassica napus L

The worldwide climate changes every year due to global warming, waterlogging, drought, salinity, pests, and pathogens, impeding crop productivity. Brassica napus is one of the most important oil crops in the world, and rapeseed oil is considered one of the most health-beneficial edible vegetable oils. Recently, miRNAs have been found and confirmed to control the expression of targets under disruptive environmental conditions. The mechanism is through the formation of the silencing complex that mediates post-transcriptional gene silencing, which pairs the target mRNA and target cleavage and/or translation inhibition. However, the functional role of miRNAs and targets in B. napus is still not clarified. This review focuses on the current knowledge of miRNAs concerning development regulation and biotic and abiotic stress responses in B. napus. Moreover, more strategies for miRNA manipulation in plants are discussed, along with future perspectives, and the enormous amount of transcriptome data available provides cues for miRNA functions in B. napus. Finally, the construction of the miRNA regulatory network can lead to the significant development of climate change-tolerant B. napus through miRNA manipulation.

RiceNCexp: a rice non-coding RNA co-expression atlas based on massive RNA-seq and small-RNA seq data

Non-coding RNAs (ncRNAs) play important roles in regulating expression of protein-coding genes. Although gene expression databases have emerged in a timely manner, a comprehensive expression database for ncRNAs is still lacking. Herein, we constructed a rice ncRNA co-expression atlas (RiceNCexp), based on 491 RNA-seq and 274 small RNA (sRNA)-seq datasets. RiceNCexp hosts four types of ncRNAs, namely lncRNAs, PHAS genes, miRNAs, and phasiRNAs. RiceNCexp provides comprehensive expression information for rice ncRNAs in 22 tissues/organs, an efficient tau-based mining tool for tissue-specific ncRNAs, and the robust co-expression analysis among ncRNAs or between ncRNAs and protein-coding genes, based on 116 pairs of RNA-seq and sRNA-seq libraries from the same experiments. In summary, RiceNCexp is a user-friendly and comprehensive rice ncRNA co-expression atlas and can be freely accessed at https://cbi.njau.edu.cn/RiceNCexp/.

Tomato MicroRNAs and Their Functions

MicroRNAs (miRNAs) define an essential class of non-coding small RNAs that function as posttranscriptional modulators of gene expression. They are coded by MIR genes, several hundreds of which exist in the genomes of Arabidopsis and rice model plants. The functional analysis of Arabidopsis and rice miRNAs indicate that their miRNAs regulate a wide range of processes including development, reproduction, metabolism, and stress. Tomato serves as a major model crop for the study of fleshy fruit development and ripening but until recently, information on the identity of its MIR genes and their coded miRNAs was limited and occasionally contradictory. As a result, the majority of tomato miRNAs remained uncharacterized. Recently, a comprehensive annotation of tomato MIR genes has been carried out by several labs and us. In this review, we curate and organize the resulting partially overlapping MIR annotations into an exhaustive and non-redundant atlas of tomato MIR genes. There are 538 candidate and validated MIR genes in the atlas, of which, 169, 18, and 351 code for highly conserved, Solanaceae-specific, and tomato-specific miRNAs, respectively. Furthermore, a critical review of functional studies on tomato miRNAs is presented, highlighting validated and possible functions, creating a useful resource for future tomato miRNA research.

TVIR: a comprehensive vegetable information resource database for comparative and functional genomic studies

Vegetables are an indispensable part of the daily diet of humans. Therefore, it is vital to systematically study the genomic data of vegetables and build a platform for data sharing and analysis. In this study, a comprehensive platform for vegetables with a user-friendly Web interface-The Vegetable Information Resource (TVIR, http://tvir.bio2db.com)-was built based on the genomes of 59 vegetables. TVIR database contains numerous important functional genes, including 5215 auxin genes, 2437 anthocyanin genes, 15 002 flowering genes, 79 830 resistance genes, and 2639 glucosinolate genes of 59 vegetables. In addition, 2597 N6-methyladenosine (m6A) genes were identified, including 513 writers, 1058 erasers, and 1026 readers. A total of 2 101 501 specific clustered regularly interspaced short palindromic repeat (CRISPR) guide sequences and 17 377 miRNAs were detected and deposited in TVIR database. Information on gene synteny, duplication, and orthologs is also provided for 59 vegetable species. TVIR database contains 2 346 850 gene annotations by the Swiss-Prot, TrEMBL, Gene Ontology (GO), Pfam, and Non-redundant (Nr) databases. Synteny, Primer Design, Blast, and JBrowse tools are provided to facilitate users in conducting comparative genomic analyses. This is the first large-scale collection of vegetable genomic data and bioinformatic analysis. All genome and gene sequences, annotations, and bioinformatic results can be easily downloaded from TVIR. Furthermore, transcriptome data of 98 vegetables have been collected and collated, and can be searched by species, tissues, or different growth stages. TVIR is expected to become a key hub for vegetable research globally. The database will be updated with newly assembled vegetable genomes and comparative genomic studies in the future.

Identification and characterization of miRNAs and PHAS loci related to the early development of the embryo and endosperm in Fragaria × ananassa

BACKGROUND

The strawberry fleshy fruit is actually enlarged receptacle tissue, and the successful development of the embryo and endosperm is essential for receptacle fruit set. MicroRNAs (miRNAs) and phased small interfering RNAs (phasiRNAs) play indispensable regulatory roles in plant growth and development. However, miRNAs and phasiRNAs participating in the regulation of strawberry embryo and endosperm development have yet to be explored.

RESULTS

Here, we performed genome-wide identification of miRNA and phasiRNA-producing loci (PHAS) in strawberry seeds with a focus on those involved in the development of the early embryo and endosperm. We found that embryos and endosperm have different levels of small RNAs. After bioinformatics analysis, the results showed that a total of 404 miRNAs (352 known and 52 novel) and 156 PHAS genes (81 21-nt and 75 24-nt genes) could be found in strawberry seed-related tissues, of which four and nine conserved miRNA families displayed conserved expression in the endosperm and embryo, respectively. Based on refined putative annotation of PHAS loci, some auxin signal-related genes, such as CM3, TAR2, AFB2, ASA1, NAC and TAS3, were found, which demonstrates that IAA biosynthesis is important for endosperm and embryo development during early fruit growth. Additionally, some auxin signal-related conserved (miR390-TAS3) and novel (miR156-ASA1) trigger-PHAS pairs were identified.

CONCLUSIONS

Taken together, these results expand our understanding of sRNAs in strawberry embryo and endosperm development and provide a genomic resource for early-stage fruit development.

MicroRNAs: emerging regulators in horticultural crops

Horticulture is one of the oldest agricultural practices with great popularity throughout the world. Horticultural crops include fruits, vegetables, ornamental plants, as well as medicinal and beverage plants. They are cultivated for food, specific nutrition, and medical use, or for aesthetic pleasure. MicroRNAs (miRNAs), which constitute a major class of endogenous small RNAs in plants, affect a multitude of developmental and physiological processes by imparting sequence specificity to gene regulation. Over the past decade, tens of thousands of miRNAs have been identified in more than 100 horticultural crops and their critical roles in regulating quality development of diverse horticultural crops have been demonstrated. Here, we review how miRNAs have emerged as important regulators and promising tools for horticultural crop improvement.

Rapid sequence and functional diversification of a miRNA superfamily targeting calcium signaling components in seed plants

MicroRNA (miRNA)-directed posttranscriptional gene silencing (miR-PTGS) is an integral component of gene regulatory networks governing plant development and responses to the environment. The sequence homology between Sly-miR4376, a miRNA common to Solanaceae and reported to target autoinhibited Ca²⁺ -ATPase 10 (ACA10) messenger RNA (mRNA) in tomato, and Arabidopsis miR391 (Ath-miR391), previously annotated as a nonconserved member of the deeply conserved miR390 family, has prompted us to revisit the function of Ath-miR391, as well as its regulatory conservation. A combination of genetic, molecular, and bioinformatic analyses revealed a hidden conservation for miR-PTGS of ACA10 homologs in spermatophytes. We found that the Arabidopsis ACA10 mRNA undergoes miR391-directed cleavage in vivo. Furthermore, transgenic overexpression of miR391 recapitulated the compact inflorescence (cif) phenotypes characteristic of ACA10 loss-of-function mutants, due to miR391-directed PTGS of ACA10. Significantly, comprehensive data mining revealed robust evidence for widespread PTGS of ACA10 homologs directed by a superfamily of related miRNAs sharing a conserved sequence core. Intriguingly, the ACA-targeting miRNAs in Poaceae also direct PTGS for calmodulin-like proteins which are putative Ca²⁺ sensors. The PTGS of ACA10 homologs is therefore directed by a miRNA superfamily that is of ancient origin and has undergone rapid sequence diversification associated with functional innovation.

Genome-Wide Analysis of miR159 Gene Family and Predicted Target Genes Associated with Environmental Stress in Dendrobium officinale: A Bioinformatics Study

Dendrobium officinale (D. officinale) is a widely used traditional Chinese medicine with high economic value. MicroR159 (miR159) is an ancient and conserved microRNA (miRNA) family in land plants, playing roles in the progress of growth and development, as well as the stress response. In order to find out functions of miR159 in D. officinale, multiple bioinformatic approaches were employed and 10 MIR159 genes were found, localizing on seven chromosomes and an unanchored segment of the D. officinale genome. All of the precursor sequences of Dof-miR159 could form a stable stem-loop structure. The phylogenetic analysis revealed that the MIR159 genes of D. officinale were divided into five clades. Furthermore, the conservation analysis suggested that the 2 to 20 nt region of miR159 mature sequences were highly conserved among family members. The promoter analysis of MIR159s showed that the majority of the predicted cis-elements were related to environmental stress or hormones. In total, five classes of genes were predicted to be the target genes of Dof-miR159s, including GAMYB transcription factors, which had been confirmed in many other land plants. The expression patterns of predicted target genes revealed their potential roles in the growth and development of D. officinale, as well as in cold and drought stress responses. In conclusion, our results illustrated the stress-related miR159-targeted genes in D. officinale, which could provide candidate genes for resistance breeding in the future.

Interactive Web-based Annotation of Plant MicroRNAs with iwa-miRNA

MicroRNAs (miRNAs) are important regulators of gene expression. The large-scale detection and profiling of miRNAs have been accelerated with the development of high-throughput small RNA sequencing (sRNA-Seq) techniques and bioinformatics tools. However, generating high-quality comprehensive miRNA annotations remains challenging due to the intrinsic complexity of sRNA-Seq data and inherent limitations of existing miRNA prediction tools. Here, we present iwa-miRNA, a Galaxy-based framework that can facilitate miRNA annotation in plant species by combining computational analysis and manual curation. iwa-miRNA is specifically designed to generate a comprehensive list of miRNA candidates, bridging the gap between already annotated miRNAs provided by public miRNA databases and new predictions from sRNA-Seq datasets. It can also assist users in selecting promising miRNA candidates in an interactive mode, contributing to the accessibility and reproducibility of genome-wide miRNA annotation. iwa-miRNA is user-friendly and can be easily deployed as a web application for researchers without programming experience. With flexible, interactive, and easy-to-use features, iwa-miRNA is a valuable tool for the annotation of miRNAs in plant species with reference genomes. We also illustrate the application of iwa-miRNA for miRNA annotation using data from plant species with varying genomic complexity. The source codes and web server of iwa-miRNA are freely accessible at http://iwa-miRNA.omicstudio.cloud/.

Dynamically expressed small RNAs, substantially driven by genomic structural variants, contribute to transcriptomic changes during tomato domestication

Tomato has undergone extensive selections during domestication. Recent progress has shown that genomic structural variants (SVs) have contributed to gene expression dynamics during tomato domestication, resulting in changes of important traits. Here, we performed comprehensive analyses of small RNAs (sRNAs) from nine representative tomato accessions. We demonstrate that SVs substantially contribute to the dynamic expression of the three major classes of plant sRNAs: microRNAs (miRNAs), phased secondary short interfering RNAs (phasiRNAs), and 24-nucleotide heterochromatic siRNAs (hc-siRNAs). Changes in the abundance of phasiRNAs and 24-nucleotide hc-siRNAs likely contribute to the alteration of mRNA gene expression in cis during tomato domestication, particularly for genes associated with biotic and abiotic stress tolerance. We also observe that miRNA expression dynamics are associated with imprecise processing, alternative miRNA-miRNA* selections, and SVs. SVs mainly affect the expression of less-conserved miRNAs that do not have established regulatory functions or low abundant members in highly expressed miRNA families. Our data highlight different selection pressures on miRNAs compared to phasiRNAs and 24-nucleotide hc-siRNAs. Our findings provide insights into plant sRNA evolution as well as SV-based gene regulation during crop domestication. Furthermore, our dataset provides a rich resource for mining the sRNA regulatory network in tomato.

Integrated Degradome and Srna Sequencing Revealed miRNA-mRNA Regulatory Networks between the Phloem and Developing Xylem of Poplar

Lignin and cellulose are the most abundant natural organic polymers in nature. MiRNAs are a class of regulatory RNAs discovered in mammals, plants, viruses, and bacteria. Studies have shown that miRNAs play a role in lignin and cellulose biosynthesis by targeting key enzymes. However, the specific miRNAs functioning in the phloem and developing xylem of Populus deltoides are still unknown. In this study, a total of 134 miRNAs were identified via high-throughput small RNA sequencing, including 132 known and two novel miRNAs, six of which were only expressed in the phloem. A total of 58 differentially expressed miRNAs (DEmiRNAs) were identified between the developing xylem and the phloem. Among these miRNAs, 21 were significantly upregulated in the developing xylem in contrast to the phloem and 37 were significantly downregulated. A total of 2431 target genes of 134 miRNAs were obtained via high-throughput degradome sequencing. Most target genes of these miRNAs were transcription factors, including AP2, ARF, bHLH, bZIP, GRAS, GRF, MYB, NAC, TCP, and WRKY genes. Furthermore, 13 and nine miRNAs were involved in lignin and cellulose biosynthesis, respectively, and we validated the miRNAs via qRT-PCR. Our study explores these miRNAs and their regulatory networks in the phloem and developing xylem of P.deltoides and provides new insight into wood formation.

miR390-tasiRNA3-ARF4 pathway is involved in regulating flowering time in woodland strawberry

miR390-tasiRNA3-ARF4 pathway was identified in woodland strawberry. FvemiR390 was involved in the regulation of flowering time, and miR390-tasiRNA3-ARF4 regulated flowering time through FveAP1/FveFUL in woodland strawberry. miRNA is an important type of regulator, and widely involved in plant growth, development and stress response. As a conserved miRNA family, the function of miR390 has been studied in many species, but poorly understood in woodland strawberry. In this study, we found that the members of miR390 family were highly conservative, and FvemiR390a and FvemiR390b have the same mature sequence. Therefore, we chose FveMIR390a to generate FvemiR390 mature sequence for functional studies. Subsequently, the result of transient gene expression assay proved that FvemiR390 negatively regulates FveARF4 through miR390-tasiRNA3-ARF4 pathway. Using transgenic plants, we discovered that the overexpression of FveMIR390a delayed flowering in woodland strawberry. Further studies revealed that the expressions of FveAP1 and FveFUL were lower in transgenic plants, which indicates miR390-tasiRNA3-ARF4 pathway delays flowering time through the FveAP1/FveFUL in woodland strawberry. Moreover, the expression of FvemiR390 responded to exogenous hormones, which also provides a reference for the application of exogenous hormones in regulating the flowering time of woodland strawberry.