Expression profiling and regulatory network of cucumber microRNAs and their putative target genes in response to cucumber green mottle mosaic virus infection

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.

Meta-Analysis of Common and Differential Transcriptomic Responses to Biotic and Abiotic Stresses in Arabidopsis thaliana

Environmental stresses adversely affect crop growth and yield, resulting in major losses to plants. These stresses occur simultaneously in nature, and we therefore conducted a meta-analysis in this study to identify differential and shared genes, pathways, and transcriptomic mechanisms involved in Arabidopsis response to biotic and abiotic stresses. The results showed a total of 436/21 significant up-/downregulated differentially expressed genes (DEGs) in response to biotic stresses, while 476 and 71 significant DEGs were respectively up- and downregulated in response to abiotic stresses in Arabidopsis thaliana. In addition, 21 DEGs (2.09%) were commonly regulated in response to biotic and abiotic stresses. Except for WRKY45 and ATXTH22, which were respectively up-/down- and down-/upregulated in response to biotic and abiotic stresses, other common DEGs were upregulated in response to all biotic and abiotic treatments. Moreover, the transcription factors (TFs) bHLH, MYB, and WRKY were the common TFs in response to biotic and abiotic stresses. In addition, ath-miR414 and ath-miR5658 were identified to be commonly expressed in response to both biotic and abiotic stresses. The identified common genes and pathways during biotic and abiotic stresses may provide potential candidate targets for the development of stress resistance breeding programs and for the genetic manipulation of crop plants.

RNA interference and crop protection against biotic stresses

RNA interference (RNAi) is a universal phenomenon of RNA silencing or gene silencing with broader implications in important physiological and developmental processes of most eukaryotes, including plants. Small RNA (sRNA) are the critical drivers of the RNAi machinery that ensures down-regulation of the target genes in a homology-dependent manner and includes small-interfering RNAs (siRNAs) and micro RNAs (miRNAs). Plant researchers across the globe have exploited the powerful technique of RNAi to execute targeted suppression of desired genes in important crop plants, with an intent to improve crop protection against pathogens and pests for sustainable crop production. Biotic stresses cause severe losses to the agricultural productivity leading to food insecurity for future generations. RNAi has majorly contributed towards the development of designer crops that are resilient towards the various biotic stresses such as viruses, bacteria, fungi, insect pests, and nematodes. This review summarizes the recent progress made in the RNAi-mediated strategies against these biotic stresses, along with new insights on the future directions in research involving RNAi for crop protection.

New early phenotypic markers for cucumber green mottle mosaic virus disease in cucumbers exposed to fluctuating extreme temperatures

Studies of early stages of cucumber green mottle mosaic virus (CGMMV) disease have been recently focused on plant molecular responses. However, extreme diurnal environmental temperatures, characteristic of global climate changes, could affect plant susceptibility and disease phenotype progression. Our studies of CGMMV disease progression, under simulated extreme temperature waves, have revealed two new disease initiation phenotypes that developed gradually, preceding severe symptom manifestations of post-recovery CGMMV systemic infections. 'Early post-recovery stage' bright yellow islands (BYIs) with defined boundaries amid asymptomatic leaf blades were first emerging followed by 'late post-recovery stage' BYIs with diffused boundaries. A deduced CGMMV disease progression scheme, postulating BYI symptom occurrence time-windows, revealed BYIs in field grown cucumber plants exposed to extreme diurnal temperatures. Profiling ontology of cucumber differentially expressed genes in BYIs vs the associated dark-green surrounding tissues disclosed activation of jasmonic acid (JA) pathway in 'early post-recovery stage' BYIs. JA signaling was inactivated in 'late post-recovery stage' BYIs concomitant with increasing expressions of JA signaling inhibitors and downregulation of JA responsive phenylpropanoid pathway. Our results disclosed a new phenotypic description of CGMMV disease initiation, characteristic of cucumbers grown under extreme environmental temperature fluctuations. The BYI phenotypes could define a time-window for CGMMV disease management applications.

Genome-Wide Characterization and Expression Analysis of the SBP-Box Gene Family in Sweet Orange (Citrus sinensis)

SBP-box is an important plant-specific transcription factor family and is involved in diverse biological processes. Here, we identified a total of 15 SBP-BOX genes in the important fruit crop sweet orange (Citrus sinensis) and characterized their gene structures, conserved domain and motif, chromosomal location, and cis-acting regulatory elements. SBP genes were classified into four subfamilies based on the amino acid sequence homology, and the classification is equally strongly supported by the gene and protein structures. Our analysis revealed that segmental duplication events were the main driving force in the evolution of CsSBP genes, and gene pairs might undergo extensive purifying selection. Further synteny analysis of the SBP members among sweet orange and other plant species provides valuable information for clarifying the CsSBP family evolutionary relationship. According to publicly available RNA-seq data and qRT-PCR analysis from various sweet orange tissues, CsSBP genes may be expressed in different tissues and developmental stages. Gene expression analysis showed variable expression profiles of CsSBP genes under various abiotic stresses, such as high and low-temperature, salt, and wound treatments, demonstrating the potential role of SBP members in sweet orange response to abiotic stress. Noticeably, all CsSBP genes were also downregulated in sweet orange upon the infection of an important fungal pathogen Diaporthe citri. Our results provide valuable information for exploring the role of SBP-Box in sweet orange.

Long non‑coding RNA HCG11 suppresses the malignant phenotype of non‑small cell lung cancer cells by targeting a miR‑875/SATB2 axis

Long non‑coding RNAs (lncRNAs) are involved in the development and progression of a variety of diseases. However, the role of the lncRNA HLA complex group 11 (HCG11) in non‑small cell lung cancer (NSCLC) remains unclear. The present study showed that the expression levels of HCG11 were reduced in tumor tissues compared with adjacent normal tissues, and similar results were obtained in experiments using lung cancer cell lines. Additionally, patients with high HCG11 expression had an increased survival rate compared with patients with low HCG11 expression. Further studies have shown that overexpression of HCG11 inhibited NSCLC cell proliferation in vitro and in vivo. Interestingly, it was observed that HCG11 expression was negatively associated with the expression levels of oncogenic microRNA‑875 (miR‑875) in patient specimens. Specifically, HCG11 served as a sponge of miR‑875. Notably, it was determined that special AT‑rich sequence‑binding protein 2 (SATB2) was a direct target gene of miR‑875, and overexpression of miR‑875 largely abrogated the effects of HCG11 in NSCLC cells. In conclusion, HCG11 was shown to suppress the malignant properties of NSCLC cells by targeting a miR‑875/SATB2 axis, and may therefore be a promising target for the treatment of NSCLC.

Cataloguing the small RNA content of honey using next generation sequencing

•

Plant miRNAs are present in Australian polyfloral and Leptospermum scoparium honey.

•

Sequencing shows that honey contains a diverse range of small, non-coding RNAs.

•

Honey RNA comes from different phylogenies including invertebrates and prokaryotes.

•

Unique small RNA profiles can provide insight into honey production conditions.

Honey adulteration is a problem that effects the global honey industry and specifically, has been discovered in the Australian market. Common methods of adulteration include dilution with sugar syrup substitutes and the mislabelling of the floral and geographic origin(s) of honey. Current authentication tools rely on the molecular variability between different honeys, identifying unique chemical profiles and/or DNA signatures characteristic of a particular honey.

Honey is known to contain plant miRNAs derived from its floral source. To explore the composition and variability of honey RNA molecules, this is the first study to catalogue the small RNA content of Australian polyfloral table honey and New Zealand Leptospermum scoparium honey using next generation sequencing. The data shows that in addition to miRNAs, honey contains a variety of small non-coding RNAs including tRNA-derived fragments. Moreover, the honey small RNAs are derived from a range of phylogenetic sources, including from plant, invertebrate, and prokaryotic species. The data indicates that different honeys contain unique small RNA profiles, which suggests a novel avenue in developing molecular-based honey authentication tools.

A tobacco ringspot virus-based vector system for gene and microRNA function studies in cucurbits

Cucurbits are economically important crops worldwide. The genomic data of many cucurbits are now available. However, functional analyses of cucurbit genes and noncoding RNAs have been impeded because genetic transformation is difficult for many cucurbitaceous plants. Here, we developed a set of tobacco ringspot virus (TRSV)-based vectors for gene and microRNA (miRNA) function studies in cucurbits. A TRSV-based expression vector could simultaneously express GREEN FLUORESCENT PROTEIN (GFP) and heterologous viral suppressors of RNA silencing in TRSV-infected plants, while a TRSV-based gene silencing vector could knock down endogenous genes exemplified by PHYTOENE DESATURASE (PDS) in Cucumis melo, Citrullus lanatus, Cucumis sativus, and Nicotiana benthamiana plants. We also developed a TRSV-based miRNA silencing vector to dissect the functions of endogenous miRNAs. Four representative miRNAs, namely, miR159, miR166, miR172, and miR319, from different cucurbits were inserted into the TRSV vector using a short tandem target mimic strategy and induced characteristic phenotypes in TRSV-miRNA-infected plants. This TRSV-based vector system will facilitate functional genomic studies in cucurbits.

Genome-Wide Changes of Regulatory Non-Coding RNAs Reveal Pollen Development Initiated at Ecodormancy in Peach

Bud dormancy is under the regulation of complex mechanisms including genetic and epigenetic factors. To study the function of regulatory non-coding RNAs in winter dormancy release, we analyzed the small RNA and long non-coding RNA (lncRNA) expression from peach (Prunus persica) floral buds in endodormancy, ecodormancy and bud break stages. Small RNAs underwent a major shift in expression primarily between dormancy and flowering with specific pairs of microRNAs and their mRNA target genes undergoing coordinated differential expression. From endodormancy to ecodormancy, ppe-miR6285 was significantly upregulated while its target gene, an ASPARAGINE-RICH PROTEIN involved in the regulation of abscisic acid signaling, was downregulated. At ecodormancy, ppe-miR2275, a homolog of meiosis-specific miR2275 across angiosperms, was significantly upregulated, supporting microsporogenesis in anthers at a late stage of dormancy. The expression of 785 lncRNAs, unlike the overall expression pattern in the small RNAs, demonstrated distinctive expression signatures across all dormancy and flowering stages. We predicted that a subset of lncRNAs were targets of microRNAs and found 18 lncRNA/microRNA target pairs with both differentially expressed across time points. The genome-wide differential expression and network analysis of non-coding RNAs and mRNAs from the same tissues provide new candidate loci for dormancy regulation and suggest complex noncoding RNA interactions control transcriptional regulation across these key developmental time points.

Differential miRNA profiles in South African cassava mosaic virus-infected cassava landraces reveal clues to susceptibility and tolerance to cassava mosaic disease

Specific miRNA families are involved in susceptibility or antiviral immunity in plants. Manihot esculenta Crantz (cassava) is a perennial plant that is an important food security crop in sub-Saharan Africa. Cassava is susceptible to several begomoviruses that cause cassava mosaic disease (CMD). In this study, we investigated the leaf miRNAome response in a tolerant (TME3) and susceptible (T200) cassava landrace challenged with South African cassava mosaic virus. RNAseq was performed on leaf samples at 12, 32 and 67 days post infection (dpi), representing early, symptomatic and late persistent stages of CMD infection. Significantly, distinct profiles of conserved miRNA family expression between the T200 and TME3 landraces at the three infection stages were observed. Notably at 12 days post SACMV infection, TME3 exhibited significant downregulation (log2fold<2.0) of 42 %, compared to 9% in T200, of the conserved miRNA families. This demonstrates an overall early response to SACMV in TME3 prior to symptom appearance not observed in T200, and expression of a large cohort of miRNA-regulated genes. Notably, at early infection, downregulation of mes-miR162 and 168 that target antiviral posttransriptional gene silencing (PTGS) regulators DCL1 and AGO1, respectively, was observed in TME3, and AGO1 and DCL1 expression was higher compared to T200 post infection. Early rapid responses prior to symptom development, including RNA silencing, may be key to establishing the tolerance/recovery phenotype exhibited by TME3 landrace later on at 67 dpi. At recovery, TME3 was hallmarked by a highly significant down-regulation of mes-miR167. MiR167 targets an auxin responsive factor which plays a role in auxin signaling and adaptive responses to stress, suggesting the importance of the auxin signaling in recovery of SACMV-induced symptoms. The gene targets of these miRNAs and their associated networks may provide clues to the molecular basis of CMD tolerance in perennial hosts such as cassava.

Transcriptome analysis provides insights into the responses of sweet potato to sweet potato virus disease (SPVD)

•

Transcriptome responses in sweet potato infected with SPCSV and/or SPFMV were studied.

•

Numerous genes, miRNAs and phasiRNAs were responsive mainly to the dual infection.

•

Salicylic acid-mediated pathways play important roles in antiviral defense responses.

Sweet potato (Ipomoea batatas) ranks among the most important crops in the world and provides nutritional and economic sustainability for subsistence farmers in sub-Saharan Africa. Its production is mainly constrained by sweet potato virus disease (SPVD) caused by the coinfection of two positive-sense single-stranded RNA viruses, sweet potato chlorotic stunt virus (SPCSV) and sweet potato feathery mottle virus (SPFMV). Current understanding of sweet potato responses to SPCSV and SPFMV at the molecular level remains very limited. In this study, we performed deep sequencing of both messenger RNA (mRNA) and small RNA (sRNA) populations in an SPVD-susceptible cultivar ‘Beauregard’ upon viral infection, to identify biological pathways that contribute to both general and specific host responses to these important viral pathogens. We found that pathways related to stress response and signaling were significantly affected by viral infection. sRNA components of these pathways were predominantly affected in late stages of the coinfection by SPCSV and SPFMV. We identified several novel microRNAs that were responsive to viral infection, some of which were predicted to target nucleotide-binding site leucine-rich repeat (NBS-LRR) disease resistance genes. The downregulation of the salicylic acid-mediated defense response pathway in particular seems to be a result of the viral infection process, and can in part explain the susceptible nature of the ‘Beauregard’ cultivar.

Up-regulation of microRNA targets correlates with symptom severity in Citrus sinensis plants infected with two different isolates of citrus psorosis virus

MAIN CONCLUSION: miRNA targets from Citrus sinensis are predicted and validated using degradome data. They show an up-regulation upon infection with CPsV, with a positive correlation between target expression and symptom severity. Sweet orange (Citrus sinensis) may suffer from disease symptoms induced by virus infections, thus resulting in drastic economic losses. Infection of sweet orange plants with two isolates of citrus psorosis virus (CPsV), expressing different symptomatologies, alters the accumulation of a set of endogenous microRNAs (miRNAs). Here, we predicted ten putative targets from four down-regulated miRNAs: three belonging to the CCAAT-binding transcription factor family (CBFAs); an Ethylene-responsive transcription factor (RAP2-7); an Integrase-type DNA-binding superfamily protein (AP2B); Transport inhibitor response 1 (TIR1); GRR1-like protein 1-related (GRR1); Argonaute 2-related (AGO2), Argonaute 7 (AGO7), and a long non-coding RNA (ncRNA). We validated six of them through analysis of leaf degradome data. Expressions of the validated targets increase in infected samples compared to healthy tissue, showing a more striking up-regulation those samples with higher symptom severity. This study contributes to the understanding of the miRNA-mediated regulation of important transcripts in Citrus sinensis through target validation and shed light in the manner a virus can alter host regulatory mechanisms leading to symptom expression.

In silico determination of transposon-derived miRNAs and targets in Aegilops species

Transposable elements (TEs) are found almost in all living organism, shaping organisms' genomes. miRNAs are noncoding RNA types which are especially important in gene expression regulations. Many previously determined plant miRNAs are identical/homologous to transposons (TE-MIR). The aim of this study was computational characterization of novel TE-related miRNAs and their targets in Aegilops genome by using stringent criteria. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed by BLAST2GO. Seventeen novel TE-related miRNAs in Aegilops genome were identified for the first time. GO analyses indicated that 40 targets played different roles in biological processes, cellular components and molecular functions. Moreover, these genes were involved in 10 metabolic pathways such as purine metabolism, nitrogen metabolism, oxidative phosphorylation, etc. as a result of KEGG analyses. Identification of miRNAs and their targets are significant to understand miRNA-TEs relationships and even how TEs affect plant growth and development. Obtaining results of this study are expected to provide possible new insight into Aegilops and its related species, wheat, with respect to miRNAs evolution and domestication.Communicated by Ramaswamy H. Sarma.