miR482 regulation of NBS-LRR defense genes during fungal pathogen infection in cotton

Opportunities and Challenges in Studies of Host-Pathogen Interactions and Management of Verticillium dahliae in Tomatoes

Tomatoes (Solanum lycopersicum L.) are a valuable horticultural crop that are grown and consumed worldwide. Optimal production is hindered by several factors, among which Verticillium dahliae, the cause of Verticillium wilt, is considered a major biological constraint in temperate production regions. V. dahliae is difficult to mitigate because it is a vascular pathogen, has a broad host range and worldwide distribution, and can persist in soil for years. Understanding pathogen virulence and genetic diversity, host resistance, and plant-pathogen interactions could ultimately inform the development of integrated strategies to manage the disease. In recent years, considerable research has focused on providing new insights into these processes, as well as the development and integration of environment-friendly management approaches. Here, we discuss the current knowledge on the race and population structure of V. dahliae, including pathogenicity factors, host genes, proteins, enzymes involved in defense, and the emergent management strategies and future research directions for managing Verticillium wilt in tomatoes.

HDAC inhibitor affects soybean miRNA482bd expression under salt and osmotic stress

MicroRNAs (miRNAs) are small non-coding molecules that modulate gene expression through targeting mRNA by specific-sequence cleavage, translation inhibition, or transcriptional regulation. miRNAs are key molecules in regulatory networks in abiotic stresses such as salt stress and water deficit in plants. Throughout the world, soybean is a critical crop, the production of which is affected by environmental stress conditions. In this study, RNA-Seq libraries from leaves of soybean under salt treatment were analyzed. 17 miRNAs and 31 putative target genes were identified with inverse differential expression patterns, indicating miRNA-target interaction. The differential expression of six miRNAs, including miR482bd-5p, and their potential targets, were confirmed by RT-qPCR. The miR482bd-5p expression was repressed, while its potential HEC1 and BAK1 targets were increased. Polyethylene glycol experiment was used to simulate drought stress, and miR482bd-5p, HEC1, and BAK1 presented a similar expression pattern, as found in salt stress. Histone modifications occur in response to abiotic stress, where histone deacetylases (HDACs) can lead to gene repression and silencing. The miR482bd-5p epigenetic regulation by histone deacetylation was evaluated by using the SAHA-HDAC inhibitor. The miR482bd-5p was up-regulated, and HEC1 was down-regulated under SAHA-salt treatment. It suggests an epigenetic regulation, where the miRNA gene is repressed by HDAC under salt stress, reducing its transcription, with an associated increase in the HEC1 target expression.

Expansion of MIR482/2118 by a class-II transposable element in cotton

Some plant microRNA (miRNA) families contain multiple members generating identical or highly similar mature miRNA variants. Mechanisms underlying the expansion of miRNA families remain elusive, although tandem and/or segmental duplications have been proposed. In this study of two tetraploid cottons, Gossypium hirsutum and Gossypium barbadense, and their extant diploid progenitors, Gossypium arboreum and Gossypium raimondii, we investigated the gain and loss of members of the miR482/2118 superfamily, which modulates the expression of nucleotide-binding site leucine-rich repeat (NBS-LRR) disease resistance genes. We found significant expansion of MIR482/2118d in G. barbadense, G. hirsutum and G. raimondii, but not in G. arboreum. Several newly expanded MIR482/2118d loci have mutated to produce different miR482/2118 variants with altered target-gene specificity. Based on detailed analysis of sequences flanking these MIR482/2118 loci, we found that this expansion of MIR482/2118d and its derivatives resulted from an initial capture of an MIR482/2118d by a class-II DNA transposable element (TE) in G. raimondii prior to the tetraploidization event, followed by transposition to new genomic locations in G. barbadense, G. hirsutum and G. raimondii. The 'GosTE' involved in the capture and proliferation of MIR482/2118d and its derivatives belongs to the PIF/Harbinger superfamily, generating a 3-bp target site duplication upon insertion at new locations. All orthologous MIR482/2118 loci in the two diploids were retained in the two tetraploids, but mutation(s) in miR482/2118 were observed across all four species as well as in different cultivars of both G. barbadense and G. hirsutum, suggesting a dynamic co-evolution of miR482/2118 and its NBS-LRR targets. Our results provide fresh insights into the mechanisms contributing to MIRNA proliferation and enrich our knowledge on TEs.

Screening of optimal reference genes for qRT-PCR and preliminary exploration of cold resistance mechanisms in Prunus mume and Prunus sibirica varieties

Prunus sibirica and Prunus mume are closely related plant species that differ in cold tolerance. Hybrids of P. sibirica and true mume, belonging to the apricot mei group, inherited strong cold resistance from P. sibirica. These materials are favourable for research on the molecular mechanisms of cold resistance. However, no suitable reference genes have been identified for analysing gene expression patterns between P. sibirica and P. mume. Ten candidate reference genes were assessed, namely, actins (ACT2-1, ACT2-2, ACT2-3, ACT2-4), protein phosphatase 2A-1 (PP2A-1), ubiquitins (UBQ2, UBQ3), ubiquitin extension protein (UBQ1) and tubulins (TUB1, TUB2), with four distinct algorithms (geNorm, NormFinder, BestKeeper and RefFinder). UBQ2 was recognized as the best reference gene in stems and buds across materials (P. sibirica; 'Xiaohong Zhusha', 'Beijing Yudie', and 'Xiao Lve' for true mume; and 'Dan Fenghou', 'Fenghou', and 'Yanxing' for apricot mei) under cold stress. In addition, the temporal and spatial expression patterns of PmCBF6 and PmLEA10 among seven varieties during winter periods were analysed using UBQ2 as a reference gene. The expression differed significantly among cultivars, which may contribute to their differences in cold tolerance. This paper confirmed the strong cold tolerance of apricot mei. And the best internal reference gene suitable for seven varieties was selected: UBQ2. Based on the above results, the expression of PmCBF6 and PmLEA10 genes during wintering in seven varieties was analysed. The molecular mechanisms of cold resistance were found to be possibly different in different varieties of P. sibirica and P. mume.

High-Throughput Sequencing and Expression Analysis Suggest the Involvement of Pseudomonas putida RA-Responsive microRNAs in Growth and Development of Arabidopsis

Beneficial soil microorganisms largely comprise of plant growth-promoting rhizobacteria (PGPR), which adhere to plant roots and facilitate their growth and development. Pseudomonas putida (RA) strain MTCC5279 is one such PGPR that exhibits several characteristics of plant growth promotion, such as P-solubilization, and siderophores and IAA production. Plant–PGPR interactions are very complex phenomena, and essentially modulate the expression of numerous genes, consequently leading to changes in the physiological, biochemical, cellular and molecular responses of plants. Therefore, in order to understand the molecular bases of plant–PGPR interactions, we carried out the identification of microRNAs from the roots of Arabidopsis upon P. putida RA-inoculation, and analyses of their expression. MicroRNAs (miRNAs) are 20- to 24-nt non-coding small RNAs known to regulate the expression of their target genes. Small RNA sequencing led to the identification of 293 known and 67 putative novel miRNAs, from the control and RA-inoculated libraries. Among these, 15 known miRNAs showed differential expression upon RA-inoculation in comparison to the control, and their expressions were corroborated by stem-loop quantitative real-time PCR. Overall, 28,746 and 6931 mRNAs were expected to be the targets of the known and putative novel miRNAs, respectively, which take part in numerous biological, cellular and molecular processes. An inverse correlation between the expression of RA-responsive miRNAs and their target genes also strengthened the crucial role of RA in developmental regulation. Our results offer insights into the understanding of the RA-mediated modulation of miRNAs and their targets in Arabidopsis, and pave the way for the further exploitation and characterization of candidate RA-responsive miRNA(s) for various crop improvement strategies directed towards plant sustainable growth and development.

Water management and phenology influence the root-associated rice field microbiota

Microbial communities associated with plants are greatly influenced by water availability in soil. In flooded crops, such as rice, the impact of water management on microbial dynamics is not fully understood. Here, we present a comprehensive study of the rice microbiota investigated in an experimental field located in one of the most productive areas of northern Italy. The microbiota associated with paddy soil and root was investigated using 454 pyrosequencing of 16S, ITS and 18S rRNA gene amplicons under two different water managements, upland (non-flooded, aerobic) and lowland (traditional flooding, anaerobic), at three plant development stages. Results highlighted a major role of the soil water status in shaping microbial communities, while phenological stage had low impacts. Compositional shifts in prokaryotic and fungal communities upon water management consisted in significant abundance changes of Firmicutes, Methanobacteria, Chloroflexi, Sordariomycetes, Dothideomycetes and Glomeromycotina. A vicariance in plant beneficial microbes and between saprotrophs and pathotrophs was observed between lowland and upland. Moreover, through network analysis, we demonstrated different co-abundance dynamics between lowland and upland conditions with a major impact on microbial hubs (strongly interconnected microbes) that fully shifted to aerobic microbes in the absence of flooding.

Transcriptional regulation of MdmiR285N microRNA in apple (Malus x domestica) and the heterologous plant system Arabidopsis thaliana

Malus x domestica microRNA MdmiR285N is a potential key regulator of plant immunity, as it has been predicted to target 35 RNA transcripts coding for different disease resistance proteins involved in plant defense to pathogens. In this study, the promoter region of MdmiR285N was isolated from the apple genome and analyzed in silico to detect potential regulatory regions controlling its transcription. A complex network of putative regulatory elements involved in plant growth and development, and in response to different hormones and stress conditions, was identified. Activity of the β-Glucoronidase (GUS) reporter gene driven by the promoter of MdmiR285N was examined in transgenic apple, demonstrating that MdmiR285N was expressed during the vegetative growth phase. Similarly, in transgenic Arabidopsis thaliana, spatial and temporal patterns of GUS expression revealed that MdmiR285N was differentially regulated during seed germination, vegetative phase change, and reproductive development. To elucidate the role of MdmiR285N in plant immunity, MdmiR285N expression in wild-type apple plants and GUS activity in transgenic apple and Arabidopsis thaliana plants were monitored in response to Erwinia amylovora and Pseudomonas syringae pv. Tomato DC3000. A significant decrease of MdmiR285N levels and GUS expression was observed during host-pathogen infections. Overall, these data suggest that MdmiR285N is involved in the biotic stress response, plant growth, and reproductive development.

Identification of vital candidate microRNA/mRNA pairs regulating ovule development using high-throughput sequencing in hazel

BACKGROUND

Hazels (Corylus spp.) are economically important nut-producing species in which ovule development determines seed plumpness, one of the key parameters reflecting nut quality. microRNAs (miRNAs) play important roles in RNA silencing and the post-transcriptional regulation of gene expression. However, very little is currently known regarding the miRNAs involved in regulating ovule growth and development.

RESULTS

In this study, we accordingly sought to determine the important miRNAs involved in ovule development and growth in hazel. We examined ovules at four developmental stages, namely, ovule formation (Ov1), early ovule growth (Ov2), rapid ovule growth (Ov3), and ovule maturity (Ov4). On the basis of small RNA and mRNA sequencing using the Illumina sequencing platform, we identified 970 miRNAs in hazel, of which 766 and 204 were known and novel miRNAs, respectively. In Ov1-vs-Ov2, Ov1-vs-Ov3, Ov1-vs-Ov4, Ov2-vs-Ov3, Ov2-vs-Ov4, and Ov3-vs-Ov4 paired comparisons, 471 differentially expressed microRNAs (DEmiRNAs) and their 3117 target differentially expressed messenger RNAs (DEmRNAs) formed 11,199 DEmiRNA/DEmRNA pairs, with each DEmiRNA changing the expression of an average of 6.62 target mRNAs. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of all DEmRNAs revealed 29 significantly enriched KEGG pathways in the six paired comparisons, including protein export (ko03060), fatty acid elongation (ko00062), starch and sucrose metabolism (ko00500), fatty acid biosynthesis (ko00061), and amino sugar and nucleotide sugar metabolism (ko00520). Our results indicate that DEmiRNA/DEmRNA pairs showing opposite change trends were related to stress tolerance, embryo and seed development, cell proliferation, auxin transduction, and the biosynthesis of proteins, starch, and fats may participate in ovule growth and development.

CONCLUSIONS

These findings contribute to a better understanding of ovule development at the level of post-transcriptional regulation, and lay the foundation for further functional analyses of hazelnut ovule growth and development.

Genome-wide identification and characterization of microRNAs by small RNA sequencing for low nitrogen stress in potato

Nitrogen is an important nutrient for plant growth and tuber quality of potato. Since potato crop requires high dose of N, improving nitrogen use efficiency (NUE) of plant is an inevitable approach to minimize N fertilization. The aim of this study was to identify and characterize microRNAs (miRNAs) by small RNA sequencing in potato plants grown in aeroponic under two contrasting N (high and low) regimes. A total of 119 conserved miRNAs belonging to 41 miRNAs families, and 1002 putative novel miRNAs were identified. From total, 52 and 54 conserved miRNAs, and 404 and 628 putative novel miRNAs were differentially expressed in roots and shoots, respectively under low N stress. Of total 34,135 predicted targets, the gene ontology (GO) analysis indicated that maximum targets belong to biological process followed by molecular function and cellular component. Eexpression levels of the selected miRNAs and targets were validated by real time-quantitative polymerase chain reaction (RT-qPCR) analysis. Two predicted targets of potential miRNAs (miR397 and miR398) were validated by 5' RLM-RACE (RNA ligase mediated rapid amplification of cDNA ends). In general, predicted targets are associated with stress-related, kinase, transporters and transcription factors such as universal stress protein, heat shock protein, salt-tolerance protein, calmodulin binding protein, serine-threonine protein kinsae, Cdk10/11- cyclin dependent kinase, amino acid transporter, nitrate transporter, sugar transporter, transcription factor, F-box family protein, and zinc finger protein etc. Our study highlights that miR397 and miR398 play crucial role in potato during low N stress management. Moreover, study provides insights to modulate miRNAs and their predicted targets to develop N-use efficient potato using transgenic/genome-editing tools in future.

Small RNA Sequencing Analysis of miRNA Expression Reveals Novel Insihts into Root Formation under Root Restriction Cultivation in Grapevine (Vitis vinifera L.)

Root restriction cultivation (RRC) can influence plant root architecture, but its root phenotypic changes and molecular mechanisms are still unknown. In this study, phenotype observations of grapevine root under RRC and control cultivation (nRC) at 12 time points were conducted, and the root phenotype showed an increase of adventitious and lateral root numbers and root tip degeneration after RRC cultivation from 70 days after planting (DAP). The 70 and 125 DAP sampling of two different cultivations, named nR70, RR70, nR125, and RR125, were selected for small RNA sequencing. A total of 153 known miRNAs and 119 predicted novel miRNAs were obtained. Furthermore, BLAST was used to predict the novel miRNAs with miRBase databases using the default parameters; 96 of the 119 predicted novel miRNAs were similar to other species, and the remaining 23 grapevine-specific novel miRNAs were obtained. There were 26, 33, 26, and 32 miRNAs that were differentially expressed in different comparison groups (RR70 vs. nR70, RR125 vs. nR125, nR125 vs. nR70 and RR125 vs. RR70). Target genes prediction of differentially expressed miRNAs was annotated on a variety of biological processes, and 24 participated in root development. Moreover, multiple miRNAs were found to jointly regulate lateral root development under root restriction conditions. The miRNA expression pattern comparison between RRC and nRC may provide a framework for the future analysis of miRNAs associated with root development in grapevine.

Integrated small RNA and mRNA expression profiles reveal miRNAs and their target genes in response to Aspergillus flavus growth in peanut seeds

BACKGROUND

MicroRNAs are important gene expression regulators in plants immune system. Aspergillus flavus is the most common causal agents of aflatoxin contamination in peanuts, but information on the function of miRNA in peanut-A. flavus interaction is lacking. In this study, the resistant cultivar (GT-C20) and susceptible cultivar (Tifrunner) were used to investigate regulatory roles of miRNAs in response to A. flavus growth.

RESULTS

A total of 30 miRNAs, 447 genes and 21 potential miRNA/mRNA pairs were differentially expressed significantly when treated with A. flavus. A total of 62 miRNAs, 451 genes and 44 potential miRNA/mRNA pairs exhibited differential expression profiles between two peanut varieties. Gene Ontology (GO) analysis showed that metabolic-process related GO terms were enriched. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses further supported the GO results, in which many enriched pathways were related with biosynthesis and metabolism, such as biosynthesis of secondary metabolites and metabolic pathways. Correlation analysis of small RNA, transcriptome and degradome indicated that miR156/SPL pairs might regulate the accumulation of flavonoids in resistant and susceptible genotypes. The miR482/2118 family might regulate NBS-LRR gene which had the higher expression level in resistant genotype. These results provided useful information for further understanding the roles of miR156/157/SPL and miR482/2118/NBS-LRR pairs.

CONCLUSIONS

Integration analysis of the transcriptome, miRNAome and degradome of resistant and susceptible peanut varieties were performed in this study. The knowledge gained will help to understand the roles of miRNAs of peanut in response to A. flavus.

The Cotton miR477-CBP60A Module Participates in Plant Defense Against Verticillium dahlia

Previous reports have shown that, when Verticillium dahliae localizes at the root surface, many microRNAs (miRNAs) were identified at the early induction stage. Here, we constructed two groups from two timepoints of small RNA (sRNA) in cotton root responses to V. dahliae at the later induction stage, pathogen localizing in the interior of root tissue. We identified 71 known and 378 novel miRNAs from six libraries of the pathogen-induced and the control sRNAs. Combined with degradome and sRNA sequencing, 178 corresponding miRNA target genes were identified, in which 40 target genes from differentially expressed miRNAs were primarily associated with oxidation-reduction and stress responses. More importantly, we characterized the cotton miR477-CBP60A module in the later response of the plant to V. dahliae infection. A β-glucuronidase fusion reporter and cleavage site analysis showed that ghr-miR477 directly cleaved the messenger RNA of GhCBP60A in the posttranscriptional process. The ghr-miR477-silencing decreased plant resistance to this fungus, while the knockdown of GhCBP60A increased plant resistance, which regulated GhICS1 expression to determine salicylic acid level. Our data documented that numerous later-inducible miRNAs in the plant response to V. dahliae, suggesting that these miRNAs play important roles in plant resistance to vascular disease.

The ghr-miR164 and GhNAC100 modulate cotton plant resistance against Verticillium dahlia

MicroRNAs (miRNAs) participate in plant development and defence through post-transcriptional regulation of the target genes. However, few miRNAs were reported to regulate cotton plant disease resistance. Here, we characterized the cotton miR164-NAC100 module in the later induction stage response of the plant to Verticillium dahliae infection. The results of GUS fusing reporter and transcript identity showed that ghr-miR164 can directly cleave the mRNA of GhNAC100 in the post-transcriptional process. The ghr-miR164 positively regulated the cotton plant resistance to V. dahliae according to analyses of its over-expression and knockdown. In link with results, the knockdown of GhNAC100 increased the plant resistance to V. dahliae. Based on LUC reporter, expression analyses and yeast one-hybrid (Y1H) assays, GhNAC100 bound to the CGTA-box of GhPR3 promoter and repressed its expression, negatively regulating plant disease resistance. These results showed that the ghr-miR164 and GhNAC100 module fine-tunes plant defence through the post-transcriptional regulation, which documented that miRNAs play important roles in plant resistance to vascular disease.

Transcriptome-wide identification and profiling of miRNAs in a stress-tolerant conifer Sabina chinensis

miRNAs are important regulatory components involving in many biological processes, including plant development, vegetative and reproductive growth, and stress response. However, identification and characterization of miRNAs still remain limited for conifer species. In this study, with deep sequencing, we obtained 1,314,450 unique reads with 18-30 nt length from a stress-tolerant conifer, Sabina chinensis. We identified 37 conserved and 103 novel miRNAs, their unique characteristics were further analyzed, and 10 randomly selected were validated by qRT-PCR. Through miRNA target predictions and annotations, we found miRNA may have several targets as well a target could be regulated by several miRNAs, and a total of 2,397 mRNAs were predicted to be targets of the 140 miRNAs. These targets included not only important transcription factors such as auxin response factors, but also indispensable non-transcriptional factor proteins. Pathway-based analysis showed that S. chinensis miRNAs are involved in 172 metabolic pathways, of which 3 were discovered in adaptation-related pathways, indicating their possible relevance to the species' stress-tolerance characteristics. This study is expected to lay the foundation for exploring the regulative roles of miRNAs in development, growth, and response to environmental stresses of S. chinensis.

Colonisation of Oncidium orchid roots by the endophyte Piriformospora indica restricts Erwinia chrysanthemi infection, stimulates accumulation of NBS-LRR resistance gene transcripts and represses their targeting micro-RNAs in leaves

BACKGROUND

Erwinia chrysanthemi (Ec) is a destructive pathogen which causes soft-rot diseases in diverse plant species including orchids. We investigated whether colonization of Oncidium roots by the endophytic fungus Piriformospora indica (Pi) restricts Ec-induced disease development in leaves, and whether this might be related to the regulation of nucleotide binding site-leucine rich repeat (NBS-LRR) Resistance (R) genes.

RESULTS

Root colonization of Oncidium stackings by Pi restricts progression of Ec-induced disease development in the leaves. Since Pi does not inhibit Ec growth on agar plates, we tested whether NBS-LRR R gene transcripts and the levels of their potential target miRNAs in Oncidium leaves might be regulated by Pi. Using bioinformatic tools, we first identified NBS-LRR R gene sequences from Oncidium, which are predicted to be targets of miRNAs. Among them, the expression of two R genes was repressed and the accumulation of several regulatory miRNA stimulated by Ec in the leaves of Oncidium plants. This correlated with the progression of disease development, jasmonic and salicylic acid accumulation, ethylene synthesis and H₂O₂ production after Ec infection of Oncidium leaves. Interestingly, root colonization by Pi restricted disease development in the leaves, and this was accompanied by higher expression levels of several defense-related R genes and lower expression level of their target miRNA.

CONCLUSION

Based on these data we propose that Pi controls the levels of NBS-LRR R mRNAs and their target miRNAs in leaves. This regulatory circuit correlates with the protection of Oncidium plants against Ec infection, and molecular and biochemical investigations will demonstrate in the future whether, and if so, to what extent these two observations are related to each other.

A transcriptomic profile of topping responsive non-coding RNAs in tobacco roots (Nicotiana tabacum)

BACKGROUND

Non-coding RNAs (ncRNAs), including microRNAs (miRNAs), long ncRNAs (lncRNAs) and circular RNAs (circRNAs), accomplish remarkable variety of biological functions. However, the composition of ncRNAs and their interactions with coding RNAs in modulating and controlling of cellular process in plants is largely unknown. Using a diverse group of high-throughput sequencing strategies, the mRNA, miRNA, lncRNA and circRNA compositions of tobacco (Nicotiana tabacum) roots determined and their alteration and potential biological functions in response to topping treatment analyzed.

RESULTS

A total of 688 miRNAs, 7423 non-redundant lncRNAs and 12,414 circRNAs were identified, among which, some selected differentially expressed RNAs were verified by quantitative real-time PCR. Using the differentially expressed RNAs, a co-expression network was established that included all four types of RNAs. The number of circRNAs identified were higher than that of miRNAs and lncRNAs, but only two circRNAs were present in the co-expression network. LncRNAs appear to be the most active ncRNAs based on their numbers presented in the co-expression network, but none of them seems to be an eTM (endogenous Target Mimicry) of miRNAs. Integrated with analyses of sequence interaction, several mRNA-circRNA-miRNA interaction networks with a potential role in the regulation of nicotine biosynthesis were uncovered, including a QS-circQS-miR6024 interaction network. In this network miR6024 was significantly down-regulated, while the expression levels of its two targets, circQS and its host gene QS, were sharply increased following the topping treatment.

CONCLUSIONS

These results illustrated the transcriptomic profiles of tobacco roots, the organ responsible for nicotine biosynthesis. mRNAs always play the most important roles, while ncRNAs are also expressed extensively for topping treatment response, especially circRNAs are the most activated in the ncRNA pool. These studies also provided insights on the coordinated regulation module of coding and non-coding RNAs in a single plant biological sample. The findings reported here indicate that ncRNAs appear to form interaction complex for the regulation of stress response forming regulation networks with transcripts involved in nicotine biosynthesis in tobacco.

Identification of microRNAs regulated by tobacco curly shoot virus co-infection with its betasatellite in Nicotiana benthamiana

BACKGROUND

MicroRNAs (miRNAs) are a class of 21-24 nucleotide endogenous non-coding small RNAs that play important roles in plant development and defense responses to biotic and abiotic stresses. Tobacco curly shoot virus (TbCSV) is a monopartite begomovirus, cause leaf curling and plant stunting symptoms in many Solanaceae plants. The betasatellite of TbCSV (TbCSB) induces more severe symptoms and enhances virus accumulation when co-infect the plants with TbCSV.

METHODS

In this study, miRNAs regulated by TbCSV and TbCSB co-infection in Nicotiana benthamiana were characterized using high-throughput sequencing technology.

RESULTS

Small RNA sequencing analysis revealed that a total of 13 known miRNAs and 42 novel miRNAs were differentially expressed in TbCSV and TbCSB co-infected N. benthamiana plants. Several potential miRNA-targeted genes were identified through data mining and were involved in both catalytic and metabolic processes, in addition to plant defense mechanisms against virus infections according to Gene Ontology (GO) analyses. In addition, the expressions of several differentially expressed miRNAs and their miRNA-targeted gene were validated through quantitative real time polymerase chain reaction (qRT-PCR) approach.

CONCLUSIONS

A large number of miRNAs are identified, and their target genes, functional annotations also have been explored. Our results provide the information on N. benthamiana miRNAs and would be useful to further understand miRNA regulatory mechanisms after TbCSV and TbCSB co-infection.

Roles of Small RNAs in Virus-Plant Interactions

Small RNAs (sRNAs), including microRNAs (miRNAs) and short interfering RNAs (siRNAs), are non-coding but powerful RNA molecules of 20–30 nucleotides in length. sRNAs play crucial regulatory roles in diverse plant biological processes. Recently, many studies on sRNAs have been reported. We summarize new findings of sRNAs in virus-plant interactions to accelerate the function analysis of sRNAs. The main content of this review article includes three parts: virus-responsive sRNAs, function analysis of sRNAs in virus pathogenicity or host resistance, and some sRNAs-mediated underlying mechanisms in virus-plant interactions. New findings of sRNAs deepen our understanding about sRNAs’ roles, which might contribute to the design of novel control measures against plant viruses.

Role of MicroRNAs and small RNAs in regulation of developmental processes and agronomic traits in Gossypium species

Small RNAs (sRNAs) are short, non-coding, 17-24 nucleotides long RNA molecules that play vital roles in regulating gene expression in every known organism investigated to date including cotton (Gossypium ssp.). These tiny RNA molecules target diverse categories of genes from different bioliogical and metabolic processes and have been reported in the three domains of life. Small RNAs, including miRNAs, are involved in ovule and fiber development, biotic and abiotic stresses, fertility, and other biochemical processes in cotton species. Also, sRNAs are the critical components in RNA interference pathway. In this article, we have reviewed the research efforts related to the isolation and characterization of miRNAs using molecular and genomic approaches. The progress made in understanding the functional roles of miRNAs in regulation, alteration, and inactivation of fundamental plant processes and traits of importance in cotton are presented here.

Identification of Fusarium graminearum-responsive miRNAs and their targets in wheat by sRNA sequencing and degradome analysis

Fusarium head blight (FHB), a prevalent disease of bread wheat (Triticum aestivum L.) caused by Fusarium graminearum, leads to considerable losses of yield and quality in wheat production. MicroRNAs (miRNAs) are important regulators of plant defense responses. Here, to better understand the F. graminearum-responsive miRNAs, we constructed sRNA libraries for wheat cultivar Sumai 3 challenged with F. graminearum and sterile water, respectively. As a result, a total of 203 known miRNAs from 46 families and 68 novel miRNAs were identified. Among them, 18 known and six novel miRNAs were found to be differentially expressed between the F. graminearum-infected samples and the controls and thus were considered to be responsive to F. graminearum. The expression patterns of eight miRNAs were further validated by stem-loop qRT-PCR. Meanwhile, target genes were validated by degradome sequencing. Integrative analysis of the differentially expressed miRNAs and their targets revealed complex miRNA-mediated regulatory networks involved in the response of wheat to F. graminearum infection. Our findings are expected to facilitate a better understanding of the miRNA regulation in wheat-F. graminearum interaction.

Identification and profiling of upland cotton microRNAs at fiber initiation stage under exogenous IAA application

Background

Cotton is the most essential textile crop worldwide, and phytohormones are critical for cotton fiber development. One example is the role of auxin in fiber initiation, but we know little molecular basis. MicroRNAs (miRNAs) have a significant function in cotton development; nevertheless their role in fiber initiation remains unclear. Here, exogenous IAA was applied to cotton plant before anthesis. Utilizing small RNA sequencing, the mechanism underlying miRNA-mediated regulation of fiber initiation under exogenous IAA treatment was investigated.

Results

With exogenous IAA application, the endogenous IAA and GA contents of IAA treated (IT) ovules were higher than control (CK) ovules at the fiber initiation stage, while endogenous ABA content was lower in IT than CK. Using scanning electron microscopy, we found the fiber number and size were significantly promoted in IT at 0 DPA. Fiber quality analysis showed that fiber length, uniformity, strength, elongation, and micronaire of IT were higher than CK, though not statistically significant, while lint percent was significantly higher in IT. We generated six small RNA libraries using − 3, 0, and 3 DPA ovules of IT and CK, and identified 58 known miRNAs and 83 novel miRNAs together with the target genes. The differential expressed miRNAs number between IT and CK at − 3, 0, 3 DPA was 34, 16 and 24, respectively. Gene ontology and KEGG pathway enrichment analyses for the target genes of the miRNAs expressed in a differential manner showed that they were significantly enriched in 30 terms and 8 pathways. QRT-PCR for those identified miRNAs and the target genes related to phytohormones and fiber development was performed, and results suggested a potential role of these miRNAs in fiber initiation.

Conclusions

The exogenous IAA application affected the relative phytohormone contents in ovule and promoted fiber initiation in cotton. Identification and profiling of miRNAs and their targets at the fiber initiation stage provided insights for miRNAs’ regulation function of fiber initiation. These findings not only shed light on the regulatory network of fiber growth but also offer clues for cotton fiber amelioration strategies in cotton.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5760-8) contains supplementary material, which is available to authorized users.

Integration of mRNA and miRNA Analysis Reveals the Molecular Mechanism Underlying Salt and Alkali Stress Tolerance in Tobacco

Salinity is one of the most severe forms of abiotic stress and affects crop yields worldwide. Plants respond to salinity stress via a sophisticated mechanism at the physiological, transcriptional and metabolic levels. However, the molecular regulatory networks involved in salt and alkali tolerance have not yet been elucidated. We developed an RNA-seq technique to perform mRNA and small RNA (sRNA) sequencing of plants under salt (NaCl) and alkali (NaHCO3) stress in tobacco. Overall, 8064 differentially expressed genes (DEGs) and 33 differentially expressed microRNAs (DE miRNAs) were identified in response to salt and alkali stress. A total of 1578 overlapping DEGs, which exhibit the same expression patterns and are involved in ion channel, aquaporin (AQP) and antioxidant activities, were identified. Furthermore, genes involved in several biological processes, such as "photosynthesis" and "starch and sucrose metabolism," were specifically enriched under NaHCO3 treatment. We also identified 15 and 22 miRNAs that were differentially expressed in response to NaCl and NaHCO3, respectively. Analysis of inverse correlations between miRNAs and target mRNAs revealed 26 mRNA-miRNA interactions under NaCl treatment and 139 mRNA-miRNA interactions under NaHCO3 treatment. This study provides new insights into the molecular mechanisms underlying the response of tobacco to salinity stress.

Integrated transcriptome, small RNA and degradome sequencing approaches provide insights into Ascochyta blight resistance in chickpea

Ascochyta blight (AB) is one of the major biotic stresses known to limit the chickpea production worldwide. To dissect the complex mechanisms of AB resistance in chickpea, three approaches, namely, transcriptome, small RNA and degradome sequencing were used. The transcriptome sequencing of 20 samples including two resistant genotypes, two susceptible genotypes and one introgression line under control and stress conditions at two time points (3rd and 7th day post inoculation) identified a total of 6767 differentially expressed genes (DEGs). These DEGs were mainly related to pathogenesis-related proteins, disease resistance genes like NBS-LRR, cell wall biosynthesis and various secondary metabolite synthesis genes. The small RNA sequencing of the samples resulted in the identification of 651 miRNAs which included 478 known and 173 novel miRNAs. A total of 297 miRNAs were differentially expressed between different genotypes, conditions and time points. Using degradome sequencing and in silico approaches, 2131 targets were predicted for 629 miRNAs. The combined analysis of both small RNA and transcriptome datasets identified 12 miRNA-mRNA interaction pairs that exhibited contrasting expression in resistant and susceptible genotypes and also, a subset of genes that might be post-transcriptionally silenced during AB infection. The comprehensive integrated analysis in the study provides better insights into the transcriptome dynamics and regulatory network components associated with AB stress in chickpea and, also offers candidate genes for chickpea improvement.

Analysis of microRNAs, phased small interfering RNAs and their potential targets in Rosarugosa Thunb

BACKGROUND

MicroRNAs (miRNAs) are small non-coding RNAs that play important roles by regulating other genes. Rosa rugosa Thunb. is an important ornamental and edible plant, yet there are only a few studies on the miRNAs and their functions in R. rugosa.

RESULTS

We sequenced 10 samll RNA profiles from the roots, petals, pollens, stamens, and leaves and 4 RNA-seq profiles in leaves and petals to analysis miRNA, phasiRNAs and mRNAs in R. rugosa. In addition, we acquired a degradome sequencing profile from leaf of R. rugosa to identify miRNA and phasiRNA targets using the SeqTar algorithm. We have identified 321 conserved miRNA homologs including primary transcripts for 25 conserved miRNAs, and 22 novel miRNAs. We identified 592 putative targets of the conserved miRNAs or tasiRNAs that showed significant accumulations of degradome reads. We found differential expression patterns of conserved miRNAs in five different tissues of R. rugosa. We identified three hundred and thirty nine 21 nucleotide (nt) PHAS loci, and forty nine 24 nt PHAS loci, respectively. Our results suggest that miR482 triggers generations of phasiRNAs by targeting nucleotide-binding, leucine-rich repeat (NB-LRR) disease resistance genes in R. rugosa. Our results also suggest that the deregulated genes in leaves and petals are significantly enriched in GO terms and KEGG pathways related to metabolic processes and photosynthesis.

CONCLUSIONS

These results significantly enhanced our knowledge of the miRNAs and phasiRNAs, as well as their potential functions, in R. rugosa.

Identification of novel miRNAs and their target genes from Populus szechuanica infected with Melampsora larici-populina

Two novel miRNAs were selected from a pre-constructed RNA library of Populus szechuanica infected with the foliar rust fungus Melampsora larici-populina in order to detect the genes regulated as targets of the miRNAs novel_mir_11 and novel_mir_357. The novel miRNAs were identified from P. szechuanica using stem-loop methods and their precursors were able to fold into a complete stem loop structure. The predicted target genes of the novel miRNAs were verified with RNA ligase-mediated 5' rapid amplification of cDNA ends (RLM-5'RACE). The full-length sequences of target genes, RPM1 and RPS2/5, in P. szechuanica were obtained through rapid amplification of cDNA ends (RACE) and officially named PsRPM1 and PsRPS2/5. These genes contain nucleotide binding site-leucine-rich repeats (NBS-LRR) domains typical of resistance genes. The expression levels of miRNAs and their target genes in different periods post infection were analysed with quantitative real-time PCR (qRT-PCR). After infection with the foliar rust fungus, the expression levels of the novel miRNAs and their target genes were dynamic. Both novel_mir_11 and novel_mir_357 negatively regulated the expression of their target genes. In this study, the regulatory effects of two novel miRNAs through their target genes were characterized to provide further mechanistic information regarding the interaction between Populus and a foliar rust fungus. Results of this study improve our understanding of the defence response mechanisms of Populus and will stimulate future work to characterize strategies to prevent and control Populus diseases.

Characterization of genome-wide microRNAs and their roles in development and biotic stress in pear

Using a genome-wide analysis of miRNAs in 'Yali' pear (Pyrus bretschneideri) via the next-generation high-throughput sequencing of small RNAs with a bioinformatics analysis, we found that pbr-miR156, pbr-miR164, pbr-miR399, and pbr-miR482 and their target genes function in viral defense in 'Duli' and 'Hongbaoshi'. pbr-miR160, pbr-miR168, pbr-miR171, and pbr-miR319 and their targets function in auxin signaling pathways in 'Zhongai 4' and 'Zhongai 5'. Successful fruit production in pear (Pyrus spp.) depends on the use of optimal combinations of rootstocks and scions. Deciphering plant-pathogen defense mechanisms and hormone signaling pathways is an important step towards developing pear rootstocks and varieties with improved qualities. In the current study, we combined next-generation sequencing of small RNAs with a bioinformatics analysis to systematically identify and characterize 298 miRNAs in the pear scion cultivar 'Yali' (Pyrus bretschneideri). We also analyzed miRNAs in three rootstock varieties ('Duli', 'Zhongai 4', and 'Zhongai 5') and one scion cultivar ('Hongbaoshi'). We found that pbr-miR156, pbr-miR164, pbr-miR399, and pbr-miR482 are induced following infection with the pear virus Apple stem pitting virus (ASPV), and identified their target genes (pbRPS6, pbNAC, pbTLR, and pbRX-CC, respectively), which participate in viral defense pathways in 'Duli' and 'Hongbaoshi'. Furthermore, we identified pbr-miR160, pbr-miR168, pbr-miR171, and pbr-miR319, and found that the production of these miRNAs was suppressed under low levels of synthetic auxin. The targets of these miRNAs (pbARF, pbAEC, pbSCL, and pbTCP4) respond to auxin signaling pathways in 'Zhongai 4' and 'Zhongai 5'. Our results lay the foundation for breeding improved pear cultivars.

Tomato lncRNA23468 functions as a competing endogenous RNA to modulate NBS-LRR genes by decoying miR482b in the tomato-Phytophthora infestans interaction

Our previous studies indicated that tomato miR482b could negatively regulate the resistance of tomato to Phytophthora infestans and the expression of miR482b was decreased after inoculation with P. infestans. However, the mechanism by which the accumulation of miR482b is suppressed remains unclear. In this study, we wrote a program to identify 89 long noncoding RNA (lncRNA)-originated endogenous target mimics (eTMs) for 46 miRNAs from our RNA-Seq data. Three tomato lncRNAs, lncRNA23468, lncRNA01308 and lncRNA13262, contained conserved eTM sites for miR482b. When lncRNA23468 was overexpressed in tomato, miR482b expression was significantly decreased, and the expression of the target genes, NBS-LRRs, was significantly increased, resulting in enhanced resistance to P. infestans. Silencing lncRNA23468 in tomato led to the increased accumulation of miR482b and decreased accumulation of NBS-LRRs, as well as reduced resistance to P. infestans. In addition, the accumulation of both miR482b and NBS-LRRs was not significantly changed in tomato plants that overexpressed lncRNA23468 with a mutated eTM site. Based on the VIGS system, a target gene of miR482b, Solyc02g036270.2, was silenced. The disease symptoms of the VIGS-Solyc02g036270.2 tomato plants were in accordance with those of tomato plants in which lncRNA23468 was silenced after inoculation with P. infestans. More severe disease symptoms were found in the modified plants than in the control plants. Our results demonstrate that lncRNAs functioning as eTMs may modulate the effects of miRNAs in tomato and provide insight into how the lncRNA23468-miR482b-NBS-LRR module regulates tomato resistance to P. infestans.

RNA Interference: A Natural Immune System of Plants to Counteract Biotic Stressors

During plant-pathogen interactions, plants have to defend the living transposable elements from pathogens. In response to such elements, plants activate a variety of defense mechanisms to counteract the aggressiveness of biotic stressors. RNA interference (RNAi) is a key biological process in plants to inhibit gene expression both transcriptionally and post-transcriptionally, using three different groups of proteins to resist the virulence of pathogens. However, pathogens trigger an anti-silencing mechanism through the expression of suppressors to block host RNAi. The disruption of the silencing mechanism is a virulence strategy of pathogens to promote infection in the invaded hosts. In this review, we summarize the RNA silencing pathway, anti-silencing suppressors, and counter-defenses of plants to viral, fungal, and bacterial pathogens.

Comparative transcriptome analysis shows the defense response networks regulated by miR482b

The transcriptomic profile in the leaves of miR482b-overexpressing tomato plants revealed that miR482b may suppress alpha-linolenic acid metabolism, cysteine and methionine metabolism, plant-pathogen interaction, and the MAPK pathway to reduce resistance to Phytophthora infestans. Our previous study showed that tomato miR482b acted as a negative regulator during tomato resistance to Phytophthora infestans by silencing NBS-LRR genes. To investigate pathways related to miR482b, the transcriptomic profile of tomato plants that overexpressed miR482b was constructed. A total of 47,124,670 raw sequence reads from the leaves of miR482b-overexpressing tomato plants were generated by Illumina sequencing. A total of 746 genes in miR482b-overexpressing tomato plants were found to show significantly differential expression relative to those in wild-type tomato plants, including 132 up-regulated genes and 614 down-regulated genes. GO and KEGG enrichment analyses showed that plant-pathogen interaction, the MAPK pathway, and the pathways related to JA and ET biosynthesis were affected by miR482b in tomato. qRT-PCR results showed that all the enriched genes in these pathways were down-regulated in tomato plants that overexpressed miR482b and up-regulated in tomato plants that overexpressed an NBS-LRR gene (Soly02g036270.2, the target gene of miR482b). After P. infestans infection, the expression of the enriched genes showed a time-dependent response, and the genes played different roles between resistant tomato (Solanum pimpinellifolium L3708) and tomato susceptible to P. infestans (S. lycopersicum Zaofen No. 2). Our results have, therefore, demonstrated that miR482b is an important component of defense response network. This will also help to identify candidate genes involved in plant-pathogen interaction.

Chromosome-based survey sequencing reveals the genome organization of wild wheat progenitor Triticum dicoccoides

Wild emmer wheat (Triticum turgidum ssp. dicoccoides) is the progenitor of wheat. We performed chromosome-based survey sequencing of the 14 chromosomes, examining repetitive sequences, protein-coding genes, miRNA/target pairs and tRNA genes, as well as syntenic relationships with related grasses. We found considerable differences in the content and distribution of repetitive sequences between the A and B subgenomes. The gene contents of individual chromosomes varied widely, not necessarily correlating with chromosome size. We catalogued candidate agronomically important loci, along with new alleles and flanking sequences that can be used to design exome sequencing. Syntenic relationships and virtual gene orders revealed several small-scale evolutionary rearrangements, in addition to providing evidence for the 4AL-5AL-7BS translocation in wild emmer wheat. Chromosome-based sequence assemblies contained five novel miRNA families, among 59 families putatively encoded in the entire genome which provide insight into the domestication of wheat and an overview of the genome content and organization.

Genome-wide identification, characterization, and evolutionary analysis of NBS-encoding resistance genes in barley

In this study, a systematic analysis of Nucleotide-Binding Site (NBS) disease resistance (R) gene family in the barley, Hordeum vulgare L. cv. Bowman, genome was performed. Using multiple computational analyses, we could identify 96 regular NBS-encoding genes and characterize them on the bases of structural diversity, conserved protein signatures, genomic distribution, gene duplications, differential expression, selection pressure, codon usage, regulation by microRNAs and phylogenetic relationships. Depending on the presence or absence of CC and LRR domains; the identified NBS genes were assigned to four distinct groups; NBS-LRR (53.1%), CC-NBS-LRR (14.6%), NBS (26%), and CC-NBS (6.3%). NBS-associated domain analysis revealed the presence of signal peptides, zinc fingers, diverse kinases, and other structural features. Eighty-five of the identified NBS-encoding genes were mapped onto the seven barley chromosomes, revealing that 50% of them were located on chromosomes 7H, 2H, and 3H, with a tendency of NBS genes to be clustered in the distal telomeric regions of the barley chromosomes. Nine gene clusters, representing 22.35% of total mapped barley NBS-encoding genes, were found, suggesting that tandem duplication stands for an important mechanism in the expansion of this gene family in barley. Phylogenetic analysis determined 31 HvNBS orthologs from rice and Brachypodium. 87 out of 96 HvNBSs were supported by expression evidence, exhibiting various and quantitatively uneven expression patterns across distinct tissues, organs, and development stages. Fourteen potential miRNA-R gene target pairs were further identified, providing insight into the regulation of NBS genes expression. These findings offer candidate target genes to engineer disease-resistant barley genotypes, and promote our understanding of the evolution of NBS-encoding genes in Poaceae crops.

A Single-Nucleotide Polymorphism in the Promoter of a Hairpin RNA Contributes to Alternaria alternata Leaf Spot Resistance in Apple (Malus × domestica)

Apple leaf spot caused by the Alternaria alternata f. sp mali (ALT1) fungus is one of the most devastating diseases of apple (Malus × domestica). We identified a hairpin RNA (hpRNA) named MdhpRNA277 that produces small RNAs and is induced by ALT1 infection in 'Golden Delicious' apple. MdhpRNA277 produces mdm-siR277-1 and mdm-siR277-2, which target five resistance (R) genes that are expressed at high levels in resistant apple variety 'Hanfu' and at low levels in susceptible variety 'Golden Delicious' following ALT1 infection. MdhpRNA277 was strongly induced in 'Golden Delicious' but not 'Hanfu' following ALT1 inoculation. MdhpRNA277 promoter activity was much stronger in inoculated 'Golden Delicious' versus 'Hanfu'. We identified a single-nucleotide polymorphism (SNP) in the MdhpRNA277 promoter region between 'Golden Delicious' (pMdhpRNA277-GD) and 'Hanfu' (pMdhpRNA277-HF). The transcription factor MdWHy binds to pMdhpRNA277-GD, but not to pMdhpRNA277-HF Transgenic 'GL-3' apple expressing pMdhpRNA277-GD:MdhpRNA277 was more susceptible to ALT1 infection than plants expressing pMdhpRNA277-HF:MdhpRNA277 due to induced mdm-siR277 accumulation and reduced expression of the five target R genes. We confirmed that the SNP in pMdhpRNA277 is associated with A. alternata leaf spot resistance by crossing. This SNP could be used as a marker to distinguish between apple varieties that are resistant or susceptible to A. alternata leaf spot.

The island cotton NBS-LRR gene GbaNA1 confers resistance to the non-race 1 Verticillium dahliae isolate Vd991

Wilt caused by Verticillium dahliae significantly reduces cotton yields, as host resistance in commercially cultivated Gossypium species is lacking. Understanding the molecular basis of disease resistance in non-commercial Gossypium species could galvanize the development of Verticillium wilt resistance in cultivated species. Nucleotide-binding site leucine-rich repeat (NBS-LRR) proteins play a central role in plant defence against pathogens. In this study, we focused on the relationship between a locus enriched with eight NBS-LRR genes and Verticillium wilt resistance in G. barbadense. Independent virus-induced gene silencing of each of the eight NBS-LRR genes in G. barbadense cultivar Hai 7124 revealed that silencing of GbaNA1 alone compromised the resistance of G. barbadense to V. dahliae isolate Vd991. In cultivar Hai 7124, GbaNA1 could be induced by V. dahliae isolate Vd991 and by ethylene, jasmonic acid and salicylic acid. Nuclear protein localization of GbaNA1 was demonstrated by transient expression. Sequencing of the GbaNA1 orthologue in nine G. hirsutum accessions revealed that all carried a non-functional allele, caused by a premature peptide truncation. In addition, all 10 G. barbadense and nine G. hirsutum accessions tested carried a full-length (∼1140 amino acids) homologue of the V. dahliae race 1 resistance gene Gbve1, although some sequence polymorphisms were observed. Verticillium dahliae Vd991 is a non-race 1 isolate that lacks the Ave1 gene. Thus, the resistance imparted by GbaNA1 appears to be mediated by a mechanism distinct from recognition of the fungal effector Ave1.

Genome-wide identification and comprehensive analysis of microRNAs and phased small interfering RNAs in watermelon

BACKGROUND

MicroRNAs (miRNAs) are a class of endogenous small non-coding RNAs involved in the post-transcriptional gene regulation and play a critical role in plant growth, development and stress responses. Watermelon (Citrullus lanatus L.) is one of the important agricultural crops worldwide. However, the watermelon miRNAs and phasiRNAs and their functions are not well explored.

RESULTS

Here we carried out computational and experimental analysis of miRNAs and phased small interfering RNAs (phasiRNAs) in watermelon by analyzing 14 small RNA profiles from roots, leaves, androecium, petals, and fruits, and one published small RNA profile of mixed tissues. To identify the targets of miRNAs and phasiRNAs, we generated a degradome profile for watermelon leaf which is analyzed using the SeqTar algorithm. We identified 97 conserved pre-miRNAs, of which 58 have not been reported previously and 348 conserved mature miRNAs without precursors. We also found 9 novel pre-miRNAs encoding 18 mature miRNAs. One hundred and one 21 nucleotide (nt) PHAS loci, and two hundred and forty one 24 nt PHAS loci were also identified. We identified 127 conserved targets of the conserved miRNAs and TAS3-derived tasiRNAs by analyzing a degradome profile of watermelon leaf.

CONCLUSIONS

The presented results provide a comprehensive view of small regulatory RNAs and their targets in watermelon.

Function identification of miR482b, a negative regulator during tomato resistance to Phytophthora infestans

Tomato is an important horticultural and economic crop cultivated worldwide. As Phytophthora infestans becomes a huge threat to tomato production, it is necessary to study the resistance mechanisms of tomato against P. infestans. Our previous research has found that miR482 might be involved in tomato-P. infestans interaction. In this study, miR482b precursor was cloned from Solanum pimpinellifolium "L3708" and miR482b was shown to decrease in abundance in tomato following P. infestans infection. Compared to wild-type tomato plants, tomato plants that overexpressed miR482b displayed more serious disease symptoms after P. infestans infection, with more necrotic cells, longer lesion diameters, and increased P. infestans abundance. Meanwhile, silencing of miR482b was performed by short tandem target mimic (STTM), resulting in enhancement of tomato resistance to P. infestans. Using miRNA and degradome data sets, NBS-LRR disease-resistance genes targeted by miR482b were validated. Negative correlation between the expression of miR482b and its target genes was found in all miR482b-overexpressing and -silencing tomato plants. Our results provide insight into tomato miR482b involved in the response to P. infestans infection, and demonstrate that miR482b-NBS-LRR is an important component in the network of tomato-P. infestans interaction.

Genetic dissection of the fuzzless seed trait in Gossypium barbadense

Cotton fibres are single-celled trichomes arising from the epidermal cells of the seed coat and may be either long (lint) or very short (fuzz). The dominant fuzzless N1 of Gossypium hirsutum is a defective allele of the At-subgenome homoeolog of MYB25-like, but the genetic components underlying the recessive fuzzless trait from G. barbadense (Gb) are unknown. We have identified five genetic loci, including a major contributing locus containing MYB25-like_Dt, associated with Gb fuzzless seeds based on genotyping of fuzzy and fuzzless near isogenic lines (NILs) from an interspecies cross (G. barbadense × G. hirsutum). At 3 d post-anthesis when fuzz fibres are initiating, expression of MYB25-like_Dt was significantly lower in fuzzless NILs than in fuzzy seeded NILs, while higher MYB25-like_Dt expression was associated with more seed fuzz across different cotton genotypes. Phenotypic and genotypic analysis of MYB25-like homoeoalleles in cottons showing different fibre phenotypes and their crossing progeny indicated that both MYB25-like_At and MYB25-like_Dt are associated with lint development, and that fuzz development is mainly determined by the expression level of MYB25-like_Dt at ~3 d post-anthesis. Expression of Gb fuzzless seeds depends on genetic background and interactions amongst the multiple loci identified. MYB25-like_Dt is one of the best candidates for N2.

Plant Responses to Pathogen Attack: Small RNAs in Focus

Small RNAs (sRNA) are a significant group of gene expression regulators for multiple biological processes in eukaryotes. In plants, many sRNA silencing pathways produce extensive array of sRNAs with specialized roles. The evidence on record advocates for the functions of sRNAs during plant microbe interactions. Host sRNAs are reckoned as mandatory elements of plant defense. sRNAs involved in plant defense processes via different pathways include both short interfering RNA (siRNA) and microRNA (miRNA) that actively regulate immunity in response to pathogenic attack via tackling pathogen-associated molecular patterns (PAMPs) and other effectors. In response to pathogen attack, plants protect themselves with the help of sRNA-dependent immune systems. That sRNA-mediated plant defense responses play a role during infections is an established fact. However, the regulations of several sRNAs still need extensive research. In this review, we discussed the topical advancements and findings relevant to pathogen attack and plant defense mediated by sRNAs. We attempted to point out diverse sRNAs as key defenders in plant systems. It is hoped that sRNAs would be exploited as a mainstream player to achieve food security by tackling different plant diseases.

Phased secondary small interfering RNAs in Panaxnotoginseng

Background

Recent results demonstrated that either non-coding or coding genes generate phased secondary small interfering RNAs (phasiRNAs) guided by specific miRNAs. Till now, there is no studies for phasiRNAs in Panax notoginseng (Burk.) F.H. Chen (P. notoginseng), an important traditional Chinese herbal medicinal plant species.

Methods

Here we performed a genome-wide discovery of phasiRNAs and its host PHAS loci in P. notoginseng by analyzing small RNA sequencing profiles. Degradome sequencing profile was used to identify the trigger miRNAs of these phasiRNAs and potential targets of phasiRNAs. We also used RLM 5’-RACE to validate some of the identified phasiRNA targets.

Results

After analyzing 24 small RNA sequencing profiles of P. notoginseng, 204 and 90 PHAS loci that encoded 21 and 24 nucleotide (nt) phasiRNAs, respectively, were identified. Furthermore, we found that phasiRNAs produced from some pentatricopeptide repeat-contain (PPR) genes target another layer of PPR genes as validated by both the degradome sequencing profile and RLM 5’-RACE analysis. We also found that miR171 with 21 nt triggers the generations of 21 nt phasiRNAs from its conserved targets.

Conclusions

We validated that some phasiRNAs generated from PPRs and TASL genes are functional by targeting other PPRs in trans. These results provide the first set of PHAS loci and phasiRNAs in P. notoginseng, and enhance our understanding of PHAS in plants.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-4331-0) contains supplementary material, which is available to authorized users.

MiRNA-seq-based profiles of miRNAs in mulberry phloem sap provide insight into the pathogenic mechanisms of mulberry yellow dwarf disease

A wide range of miRNAs have been identified as phloem-mobile molecules that play important roles in coordinating plant development and physiology. Phytoplasmas are associated with hundreds of plant diseases, and the pathogenesis involved in the interactions between phytoplasmas and plants is still poorly understood. To analyse the molecular mechanisms of phytoplasma pathogenicity, the miRNAs profiles in mulberry phloem saps were examined in response to phytoplasma infection. A total of 86 conserved miRNAs and 19 novel miRNAs were identified, and 30 conserved miRNAs and 13 novel miRNAs were differentially expressed upon infection with phytoplasmas. The target genes of the differentially expressed miRNAs are involved in diverse signalling pathways showing the complex interactions between mulberry and phytoplasma. Interestingly, we found that mul-miR482a-5p was up-regulated in the infected phloem saps, and grafting experiments showed that it can be transported from scions to rootstock. Based on the results, the complexity and roles of the miRNAs in phloem sap and the potential molecular mechanisms of their changes were discussed. It is likely that the phytoplasma-responsive miRNAs in the phloem sap modulate multiple pathways and work cooperatively in response to phytoplasma infection, and their expression changes may be responsible for some symptoms in the infected plants.

Small RNAs, emerging regulators critical for the development of horticultural traits

Small RNAs (sRNAs) have been recently recognized as key genetic and epigenetic regulators in various organisms, ranging from the modification of DNA and histone methylations to the modulation of the abundance of coding or non-coding RNAs. In plants, major regulatory sRNAs are classified as respective microRNA (miRNA) and small interfering RNA (siRNA) species, with the former primarily engaging in posttranscriptional regulation while the latter in transcriptional one. Many of these characterized sRNAs are involved in regulation of diverse biological programs, processes, and pathways in response to developmental cues, environmental signals/stresses, pathogen infection, and pest attacks. Recently, sRNAs-mediated regulations have also been extensively investigated in horticultural plants, with many novel mechanisms unveiled, which display far more mechanistic complexity and unique regulatory features compared to those studied in model species. Here, we review the recent progress of sRNA research in horticultural plants, with emphasis on mechanistic aspects as well as their relevance to trait regulation. Given that major and pioneered sRNA research has been carried out in the model and other plants, we also discuss ongoing sRNA research on these plants. Because miRNAs and phased siRNAs (phasiRNAs) are the most studied sRNA regulators, this review focuses on their biogenesis, conservation, function, and targeted genes and traits as well as the mechanistic relation between them, aiming at providing readers comprehensive information instrumental for future sRNA research in horticulture crops.

Roles of Non-Coding RNA in Sugarcane-Microbe Interaction

Studies have highlighted the importance of non-coding RNA regulation in plant-microbe interaction. However, the roles of sugarcane microRNAs (miRNAs) in the regulation of disease responses have not been investigated. Firstly, we screened the sRNA transcriptome of sugarcane infected with Acidovorax avenae. Conserved and novel miRNAs were identified. Additionally, small interfering RNAs (siRNAs) were aligned to differentially expressed sequences from the sugarcane transcriptome. Interestingly, many siRNAs aligned to a transcript encoding a copper-transporter gene whose expression was induced in the presence of A. avenae, while the siRNAs were repressed in the presence of A. avenae. Moreover, a long intergenic non-coding RNA was identified as a potential target or decoy of miR408. To extend the bioinformatics analysis, we carried out independent inoculations and the expression patterns of six miRNAs were validated by quantitative reverse transcription-PCR (qRT-PCR). Among these miRNAs, miR408-a copper-microRNA-was downregulated. The cleavage of a putative miR408 target, a laccase, was confirmed by a modified 5'RACE (rapid amplification of cDNA ends) assay. MiR408 was also downregulated in samples infected with other pathogens, but it was upregulated in the presence of a beneficial diazotrophic bacteria. Our results suggest that regulation by miR408 is important in sugarcane sensing whether microorganisms are either pathogenic or beneficial, triggering specific miRNA-mediated regulatory mechanisms accordingly.

ghr-miR5272a-mediated regulation of GhMKK6 gene transcription contributes to the immune response in cotton

Fusarium wilt is a major biotic stress affecting the productivity of cotton (Gossypium hirsutum). Although mitogen-activated protein kinase (MAPK) cascades play critical roles in plant disease resistance, their intricate regulation under fungal stress remains unclear, especially with regards to microRNA-mediated regulation of MAPK gene expression. In this study, we report that the MAPK kinase gene GhMKK6 and ghr-miR5272a work together in cotton resistance to Fusarium wilt. Silencing GhMKK6 in cotton decreased resistance to F. oxysporum by repressing the expression of known disease-resistance genes. Furthermore, although GhMKK6 played a positive role in disease resistance, excessive GhMKK6 activation caused an excessive hypersensitive response. ghr-miR5272a, a major regulator, prevents this excessive response by regulating GhMKK6 expression. ghr-miR5272a targets the GhMKK6 3'-untranslated region in cotton. Overexpressing miR5272a decreased the expression of GhMKK6 and disease-resistance genes, and increased sensitivity to F. oxysporum, yielding a similar phenotype to GhMKK6-silenced cotton. Overall, these results demonstrate that the ghr-miR5272a-mediated regulation of GhMKK6 expression contributes to the immune response in cotton, and reveal a new feedback loop mechanism in plant disease response.

Resistance Genes in Piper colubrinum: In Silico Survey From Leaf Transcriptome and Expression Studies Upon Challenge Inoculation with Phytophthora capsici

The oomycetes, Phytophthora capsici, cause foot rot disease in black pepper. Piper colubrinum Link, a distant relative of cultivated black pepper, is highly resistant to this destructive pathogen. Identification of resistance (R) genes in P. colubrinum and the study of its expression profile during interaction with the pathogen can help in understanding the resistance mechanism involved. In the present study, 1289 R gene-related transcripts were mined from P. colubrinum transcriptome, clustered, and classified according to the conserved motifs and domains. Transcripts belonging to four major R gene classes were identified in P. colubrinum, but TIR-NBS-LRR-type R genes were absent. The relative expression of 12 selected R genes was studied using two virulent isolates of P. capsici, and these were found to be upregulated in the initial hours of plant pathogen interaction. The R genes studied were expressed even in aseptically maintained tissue-cultured plants and uninoculated greenhouse-grown plants at basal level suggesting that the plants are geared up with the R gene all the time and are under continuous surveillance for the pathogen and basal level of R gene expression do not require a pathogen trigger. ACT, ATUB, and EIF3E were identified as the most stable reference genes that can be used for real-time PCR study. The present study identified promising R genes in P. colubrinum which can be used in developing Phytophthora-resistant black pepper.

MicroRNA-Mediated Gene Silencing in Plant Defense and Viral Counter-Defense

MicroRNAs (miRNAs) are non-coding RNAs of approximately 20–24 nucleotides in length that serve as central regulators of eukaryotic gene expression by targeting mRNAs for cleavage or translational repression. In plants, miRNAs are associated with numerous regulatory pathways in growth and development processes, and defensive responses in plant–pathogen interactions. Recently, significant progress has been made in understanding miRNA-mediated gene silencing and how viruses counter this defense mechanism. Here, we summarize the current knowledge and recent advances in understanding the roles of miRNAs involved in the plant defense against viruses and viral counter-defense. We also document the application of miRNAs in plant antiviral defense. This review discusses the current understanding of the mechanisms of miRNA-mediated gene silencing and provides insights on the never-ending arms race between plants and viruses.

The use of high-throughput small RNA sequencing reveals differentially expressed microRNAs in response to aster yellows phytoplasma-infection in Vitis vinifera cv. 'Chardonnay'

Phytoplasmas are cell wall-less plant pathogenic bacteria responsible for major crop losses throughout the world. In grapevine they cause grapevine yellows, a detrimental disease associated with a variety of symptoms. The high economic impact of this disease has sparked considerable interest among researchers to understand molecular mechanisms related to pathogenesis. Increasing evidence exist that a class of small non-coding endogenous RNAs, known as microRNAs (miRNAs), play an important role in post-transcriptional gene regulation during plant development and responses to biotic and abiotic stresses. Thus, we aimed to dissect complex high-throughput small RNA sequencing data for the genome-wide identification of known and novel differentially expressed miRNAs, using read libraries constructed from healthy and phytoplasma-infected Chardonnay leaf material. Furthermore, we utilised computational resources to predict putative miRNA targets to explore the involvement of possible pathogen response pathways. We identified multiple known miRNA sequence variants (isomiRs), likely generated through post-transcriptional modifications. Sequences of 13 known, canonical miRNAs were shown to be differentially expressed. A total of 175 novel miRNA precursor sequences, each derived from a unique genomic location, were predicted, of which 23 were differentially expressed. A homology search revealed that some of these novel miRNAs shared high sequence similarity with conserved miRNAs from other plant species, as well as known grapevine miRNAs. The relative expression of randomly selected known and novel miRNAs was determined with real-time RT-qPCR analysis, thereby validating the trend of expression seen in the normalised small RNA sequencing read count data. Among the putative miRNA targets, we identified genes involved in plant morphology, hormone signalling, nutrient homeostasis, as well as plant stress. Our results may assist in understanding the role that miRNA pathways play during plant pathogenesis, and may be crucial in understanding disease symptom development in aster yellows phytoplasma-infected grapevines.