Evolution of NBS-LRR gene copies among Dicot plants and its regulation by members of the miR482/2118 superfamily of miRNAs

The role of microRNAs in NBS-LRR gene expression and its implications for plant immunity and crop development

Plants evolved, over millions of years, complex defense systems against pathogens. Once infected, the interaction between pathogen effector molecules and host receptors triggers plant immune responses, which include apoptosis, systemic immune response, among others. An important protein family responsible for pathogen effector recognition is the nucleotide binding site-leucine repeat rich (NBS-LRR) proteins. The NBS-LRR gene family is the largest disease resistance gene class in plants. These proteins are widely distributed in vascular plants and have a complex multigenic cluster distribution in plant genomes. To counteract the genetic load of such a large gene family on fitness cost, plants evolved a mechanism using post transcriptional gene silencing induced by small RNAs, particularly microRNAs. For the NBS-LRR gene family, the small RNAs involved in this silencing mechanism are mainly the microRNA482/2118 superfamily. This suppression mechanism is relieved upon pathogen infection, thus allowing increased NBS-LRR expression and triggering plant immunity. In this review, we will discuss the biogenesis of microRNAs and secondary RNAs involved in this silencing mechanism, biochemical and structural features of NBS-LRR proteins in response to pathogen effectors and the evolution of microRNA-based silencing mechanism with a focus on the miR482/2118 family. Furthermore, the biotechnological manipulation of microRNA expression, using both transgenic or genome editing approaches to improve cultivated plants will be discussed, with a focus on the miR482/2118 family in soybean.

Impact of biotic stresses on the Brassicaceae family and opportunities for crop improvement by exploiting genotyping traits

MAIN CONCLUSION

Utilizing RNAi, miRNA, siRNA, lncRNA and exploiting genotyping traits can help safeguard the food supply from illnesses and pest damage to Brassicas, as well as reduce yield losses caused by plant pathogens and insect pests. In the natural environment, plants face significant challenges in the form of biotic stress, due to various living organisms, leading to biological stress and a sharp decline in crop yields. To cope with these effects, plants have evolved specialized mechanisms to mitigate these challenges. Plant stress tolerance and resistance are influenced by genes associated with stress-responsive pathogens that interact with various stress-related signaling pathway components. Plants employ diverse strategies and mechanisms to combat biological stress, involving a complex network of transcription factors that interact with specific cis-elements to regulate gene expression. Understanding both plant developmental and pathogenic disease resistance mechanisms can allow us to develop stress-tolerant and -resistant crops. Brassica genus includes commercially important crops, e.g., broccoli, cabbage, cauliflower, kale, and rapeseed, cultivated worldwide, with several applications, e.g., oil production, consumption, condiments, fodder, as well as medicinal ones. Indeed, in 2020, global production of vegetable Brassica reached 96.4 million tones, a 10.6% rise from the previous decade. Taking into account their commercial importance, coupled to the impact that pathogens can have in Brassica productivity, yield losses up to 60%, this work complies the major diseases caused due to fungal, bacterial, viral, and insects in Brassica species. The review is structured into three parts. In the first part, an overview is provided of the various pathogens affecting Brassica species, including fungi, bacteria, viruses, and insects. The second part delves into the exploration of defense mechanisms that Brassica plants encounter against these pathogens including secondary metabolites, duplicated genes, RNA interference (RNAi), miRNA (micro-RNA), siRNA (small interfering RNA), and lncRNA (long non-coding RNA). The final part comprehensively outlines the current applications of CRISPR/Cas9 technology aimed at enhancing crop quality. Taken collectively, this review will contribute to our enhanced understanding of these mechanisms and their role in the development of resistance in Brassica plants, thus supporting strategies to protect this crucial crop.

Regulation of plant immunity via small RNA-mediated control of NLR expression

Plants use different receptors to detect potential pathogens: membrane-anchored pattern recognition receptors (PRRs) activated upon perception of pathogen-associated molecular patterns (PAMPs) that elicit pattern-triggered immunity (PTI); and intracellular nucleotide-binding leucine-rich repeat proteins (NLRs) activated by detection of pathogen-derived effectors, activating effector-triggered immunity (ETI). The interconnections between PTI and ETI responses have been increasingly reported. Elevated NLR levels may cause autoimmunity, with symptoms ranging from fitness cost to developmental arrest, sometimes combined with run-away cell death, making accurate control of NLR dosage key for plant survival. Small RNA-mediated gene regulation has emerged as a major mechanism of control of NLR dosage. Twenty-two nucleotide miRNAs with the unique ability to trigger secondary siRNA production from target transcripts are particularly prevalent in NLR regulation. They enhance repression of the primary NLR target, but also bring about repression of NLRs only complementary to secondary siRNAs. We summarize current knowledge on miRNAs and siRNAs in the regulation of NLR expression with an emphasis on 22 nt miRNAs and propose that miRNA and siRNA regulation of NLR levels provides additional links between PTI and NLR defense pathways to increase plant responsiveness against a broad spectrum of pathogens and control an efficient deployment of defenses.

Molecular perspectives on age-related resistance of plants to (viral) pathogens

Age-related resistance to microbe invasion is a commonly accepted concept in plant pathology. However, the impact of such age-dependent interactive phenomena is perhaps not yet sufficiently recognized by the broader plant science community. Toward cataloging an understanding of underlying mechanisms, this review explores recent molecular studies and their relevance to the concept. Examples describe differences in genetic background, transcriptomics, hormonal balances, protein-mediated events, and the contribution by short RNA-controlled gene silencing events. Throughout, recent findings with viral systems are highlighted as an illustration of the complexity of the interactions. It will become apparent that instead of uncovering a unifying explanation, we unveiled only trends. Nevertheless, with a degree of confidence, we propose that the process of plant age-related defenses is actively regulated at multiple levels. The overarching goal of this control for plants is to avoid a constitutive waste of resources, especially at crucial metabolically draining early developmental stages.

The Role of Viruses in Identifying and Analyzing RNA Silencing

Adaptive antiviral immunity in plants is an RNA-based mechanism in which small RNAs derived from both strands of the viral RNA are guides for an Argonaute (AGO) nuclease. The primed AGO specifically targets and silences the viral RNA. In plants this system has diversified to involve mobile small interfering RNAs (siRNAs), an amplification system involving secondary siRNAs and targeting mechanisms involving DNA methylation. Most, if not all, plant viruses encode multifunctional proteins that are suppressors of RNA silencing that may also influence the innate immune system and fine-tune the virus-host interaction. Animal viruses similarly trigger RNA silencing, although it may be masked in differentiated cells by the interferon system and by the action of the virus-encoded suppressor proteins. There is huge potential for RNA silencing to combat viral disease in crops, farm animals, and people, although there are complications associated with the various strategies for siRNA delivery including transgenesis. Alternative approaches could include using breeding or small molecule treatment to enhance the inherent antiviral capacity of infected cells.

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

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

Roles of RNA silencing in viral and non-viral plant immunity and in the crosstalk between disease resistance systems

RNA silencing is a well-established antiviral immunity system in plants, in which small RNAs guide Argonaute proteins to targets in viral RNA or DNA, resulting in virus repression. Virus-encoded suppressors of silencing counteract this defence system. In this Review, we discuss recent findings about antiviral RNA silencing, including the movement of RNA through plasmodesmata and the differentiation between plant self and viral RNAs. We also discuss the emerging role of RNA silencing in plant immunity against non-viral pathogens. This immunity is mediated by transkingdom movement of RNA into and out of the infected plant cells in vesicles or as extracellular nucleoproteins and, like antiviral immunity, is influenced by the silencing suppressors encoded in the pathogens' genomes. Another effect of RNA silencing on general immunity involves host-encoded small RNAs, including microRNAs, that regulate NOD-like receptors and defence signalling pathways in the innate immunity system of plants. These RNA silencing pathways form a network of processes with both positive and negative effects on the immune systems of plants.

Evolutionary footprint of plant immunity

There are pieces of evidence from genomic footprints and fossil records indicating that plants have co-evolved with microbes after terrestrialization for more than 407 million years. Therefore, to truly comprehend plant evolution, we need to understand the co-evolutionary process and history between plants and microbes. Recent developments in genomes and transcriptomes of a vast number of plant species as well as microbes have greatly expanded our knowledge of the evolution of the plant immune system. In this review, we summarize recent advances in the co-evolution between plants and microbes with emphasis on the plant side and point out future research needed for understanding plant-microbial co-evolution. Knowledge of the evolution and variation of the plant immune system will better equip us on designing crops with boosted performance in agricultural fields.

MicroRNA482/2118, a miRNA superfamily essential for both disease resistance and plant development

MicroRNAs (miRNAs) are a class of 21-24 nucleotides (nt) noncoding small RNAs ubiquitously distributed across the plant kingdom. miR482/2118, one of the conserved miRNA superfamilies originating from gymnosperms, has divergent main functions in core-angiosperms. It mainly regulates NUCLEOTIDE BINDING SITE-LEUCINE-RICH REPEAT (NBS-LRR) genes in eudicots, functioning as an essential component in plant disease resistance; in contrast, it predominantly targets numerous long noncoding RNAs (lncRNAs) in monocot grasses, which are vital for plant reproduction. Usually, miR482/2118 is 22-nt in length, which can trigger the production of phased small interfering RNAs (phasiRNAs) after directed cleavage. PhasiRNAs instigated from target genes of miR482/2118 enhance their roles in corresponding biological processes by cis-regulation on cognate genes and expands their function to other pathways via trans activity on different genes. This review summarizes the origin, biogenesis, conservation, and evolutionary characteristics of the miR482/2118 superfamily and delineates its diverse functions in disease resistance, plant development, stress responses, etc.

Identification of Micro Ribonucleic Acids and Their Targets in Response to Plasmodiophora brassicae Infection in Brassica napus

Clubroot disease, which is caused by the soil-borne pathogen Plasmodiophora brassicae War (P. brassicae), is one of the oldest and most destructive diseases of Brassica and cruciferous crops in the world. Plant microRNAs [micro ribonucleic acids (miRNAs)] play important regulatory roles in several developmental processes. Although the role of plant miRNAs in plant-microbe interaction has been extensively studied, there are only few reports on the specific functions of miRNAs in response to P. brassicae. This study investigated the roles of miRNAs and their targets during P. brassicae infection in a pair of Brassica napus near-isogenic lines (NILs), namely clubroot-resistant line 409R and clubroot-susceptible line 409S. Small RNA sequencing (sRNA-seq) and degradome-seq were performed on root samples of 409R and 409S with or without P. brassicae inoculation. sRNA-seq identified a total of 48 conserved and 72 novel miRNAs, among which 18 had a significant differential expression in the root of 409R, while only one miRNA was differentially expressed in the root of 409S after P. brassicae inoculation. The degradome-seq analysis identified 938 miRNA target transcripts, which are transcription factors, enzymes, and proteins involved in multiple biological processes and most significantly enriched in the plant hormone signal transduction pathway. Between 409R and 409S, we found eight different degradation pathways in response to P. brassicae infection, such as those related to fatty acids. By combining published transcriptome data, we identified a total of six antagonistic miRNA-target pairs in 409R that are responsive to P. brassicae infection and involved in pathways associated with root development, hypersensitive cell death, and chloroplast metabolic synthesis. Our results reveal that P. brassicae infection leads to great changes in miRNA pool and target transcripts. More interestingly, these changes are different between 409R and 409S. Clarification of the crosstalk between miRNAs and their targets may shed new light on the possible mechanisms underlying the pathogen resistance against P. brassicae.

NB-LRR genes: characteristics in three Solanum species and transcriptional response to Ralstonia solanacearum in tomato

NB-LRR genes in the three Solanum species showed specific constitution characteristics and evolved multiple clusters and duplicates. Some genes could respond to biotic stresses such as tomato bacterial wilt. Nucleotide-binding and leucine-rich repeat (NB-LRR, NLR) is a largest resistance gene family in plants, which plays a key role in response to biotic stresses. In this study, NB-LRR genes in cultivated tomato Solanum lycopersicum (Sl) and its wild relatives S. pennellii (Spe) and S. pimpinellifolium (Spi) were analyzed using bioinformatics approaches. In total, 238, 202 and 217 NB-LRR genes of 8 different types were found in Sl, Spe and Spi, respectively. The three species showed similar genomic characteristics. The NB-LRR genes were mainly distributed on chromosomes 4, 5 and 11 and located at the distal zones, forming multiple clusters and tandem duplicates. A large number of homologs appeared through gene expansion, with most Ka/Ks values being less than 1, indicating that purifying selection had occurred in evolution. These genes were mainly expressed in root and could respond to different biotic stresses. RT-qPCR analysis revealed that SlNLR genes could respond to tomato bacterial wilt, with SlNLR1 probably involved in the resistance response, whereas others being the opposite. The transcription factors (TFs) and interaction proteins that regulate target genes were mainly Dof, NAC and MYB families and kinases. The results provide a basis for the isolation and application of related genes in plant disease resistance breeding.

MiR1885 Regulates Disease Tolerance Genes in Brassica rapa during Early Infection with Plasmodiophora brassicae

Clubroot caused by Plasmodiophora brassicae is a severe disease of cruciferous crops that decreases crop quality and productivity. Several clubroot resistance-related quantitative trait loci and candidate genes have been identified. However, the underlying regulatory mechanism, the interrelationships among genes, and how genes are regulated remain unexplored. MicroRNAs (miRNAs) are attracting attention as regulators of gene expression, including during biotic stress responses. The main objective of this study was to understand how miRNAs regulate clubroot resistance-related genes in P. brassicae-infected Brassica rapa. Two Brassica miRNAs, Bra-miR1885a and Bra-miR1885b, were revealed to target TIR-NBS genes. In non-infected plants, both miRNAs were expressed at low levels to maintain the balance between plant development and basal immunity. However, their expression levels increased in P. brassicae-infected plants. Both miRNAs down-regulated the expression of the TIR-NBS genes Bra019412 and Bra019410, which are located at a clubroot resistance-related quantitative trait locus. The Bra-miR1885-mediated down-regulation of both genes was detected for up to 15 days post-inoculation in the clubroot-resistant line CR Shinki and in the clubroot-susceptible line 94SK. A qRT-PCR analysis revealed Bra019412 expression was negatively regulated by miR1885. Both Bra019412 and Bra019410 were more highly expressed in CR Shinki than in 94SK; the same expression pattern was detected in multiple clubroot-resistant and clubroot-susceptible inbred lines. A 5′ rapid amplification of cDNA ends analysis confirmed the cleavage of Bra019412 by Bra-miR1885b. Thus, miR1885s potentially regulate TIR-NBS gene expression during P. brassicae infections of B. rapa.

Roles of small RNAs in crop disease resistance

Small RNAs (sRNAs) are a class of short, non-coding regulatory RNAs that have emerged as critical components of defense regulatory networks across plant kingdoms. Many sRNA-based technologies, such as host-induced gene silencing (HIGS), spray-induced gene silencing (SIGS), virus-induced gene silencing (VIGS), artificial microRNA (amiRNA) and synthetic trans-acting siRNA (syn-tasiRNA)-mediated RNA interference (RNAi), have been developed as disease control strategies in both monocot and dicot plants, particularly in crops. This review aims to highlight our current understanding of the roles of sRNAs including miRNAs, heterochromatic siRNAs (hc-siRNAs), phased, secondary siRNAs (phasiRNAs) and natural antisense siRNAs (nat-siRNAs) in disease resistance, and sRNAs-mediated trade-offs between defense and growth in crops. In particular, we focus on the diverse functions of sRNAs in defense responses to bacterial and fungal pathogens, oomycete and virus in crops. Further, we highlight the application of sRNA-based technologies in protecting crops from pathogens. Further research perspectives are proposed to develop new sRNAs-based efficient strategies to breed non-genetically modified (GMO), disease-tolerant crops for sustainable agriculture.

Evolutionary and Characteristic Analysis of RING-DUF1117 E3 Ubiquitin Ligase Genes in Gossypium Discerning the Role of GhRDUF4D in Verticillium dahliae Resistance

Verticillium wilt, primarily induced by the soil-borne fungus Verticillium dahliae, is a serious threat to cotton fiber production. There are a large number of really interesting new gene (RING) domain-containing E3 ubiquitin ligases in Arabidopsis, of which three (At2g39720 (AtRHC2A), At3g46620 (AtRDUF1), and At5g59550 (AtRDUF2)) have a domain of unknown function (DUF) 1117 domain in their C-terminal regions. This study aimed to detect and characterize the RDUF members in cotton, to gain an insight into their roles in cotton’s adaptation to environmental stressors. In this study, a total of 6, 7, 14, and 14 RDUF (RING-DUF1117) genes were detected in Gossypium arboretum, G. raimondii, G. hirsutum, and G. barbadense, respectively. These RDUF genes were classified into three groups. The genes in each group were highly conserved based on gene structure and domain analysis. Gene duplication analysis revealed that segmental duplication occurred during cotton evolution. Expression analysis revealed that the GhRDUF genes were widely expressed during cotton growth and under abiotic stresses. Many cis-elements related to hormone response and environment stressors were identified in GhRDUF promoters. The predicted target miRNAs and transcription factors implied that GhRDUFs might be regulated by gra-miR482c, as well as by transcription factors, including MYB, C2H2, and Dof. The GhRDUF genes responded to cold, drought, and salt stress and were sensitive to jasmonic acid, salicylic acid, and ethylene signals. Meanwhile, GhRDUF4D expression levels were enhanced after V. dahliae infection. Subsequently, GhRDUF4D was verified by overexpression in Arabidopsis and virus-induced gene silencing treatment in upland cotton. We observed that V. dahliae resistance was significantly enhanced in transgenic Arabidopsis, and weakened in GhRDUF4D silenced plants. This study conducted a comprehensive analysis of the RDUF genes in Gossypium, hereby providing basic information for further functional studies.

miRNA Mediated Regulation and Interaction between Plants and Pathogens

Plants have evolved diverse molecular mechanisms that enable them to respond to a wide range of pathogens. It has become clear that microRNAs, a class of short single-stranded RNA molecules that regulate gene expression at the transcriptional or post-translational level, play a crucial role in coordinating plant-pathogen interactions. Specifically, miRNAs have been shown to be involved in the regulation of phytohormone signals, reactive oxygen species, and NBS-LRR gene expression, thereby modulating the arms race between hosts and pathogens. Adding another level of complexity, it has recently been shown that specific lncRNAs (ceRNAs) can act as decoys that interact with and modulate the activity of miRNAs. Here we review recent findings regarding the roles of miRNA in plant defense, with a focus on the regulatory modes of miRNAs and their possible applications in breeding pathogen-resistance plants including crops and trees. Special emphasis is placed on discussing the role of miRNA in the arms race between hosts and pathogens, and the interaction between disease-related miRNAs and lncRNAs.

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.

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.

Suppression of NB-LRR genes by miRNAs promotes nitrogen-fixing nodule development in Medicago truncatula

Plant genomes contain two major classes of innate immune receptors to recognize different pathogens. The pattern recognition receptors perceive conserved pathogen-associated molecular patterns and the resistance genes with nucleotide-binding (NB) and leucine-rich repeat (LRR) domains recognize specific pathogen effectors. The precise regulation of resistance genes is important since the unregulated expression of NB-LRR genes can inhibit growth and may result in autoimmunity in the absence of pathogen infection. It was shown that a subset of miRNAs could target NB-LRR genes and act as an important regulator of plant immunity in the absence of pathogens. Plants not only interact with pathogens, but they can also establish symbiotic interactions with microbes. Nitrogen-fixing symbiotic interaction and nodule formation of legumes may also require the suppression of host defence to prevent immune responses. We found that upon symbiotic interactions, miRNAs repressing NB-LRR expression are upregulated in the developing nodules of Medicago truncatula. Furthermore, we show that the suppression of the activity of the NB-LRR genes targeted by these miRNAs is important during nodule development. Our results suggest that the downregulation of NB-LRR resistance genes in the developing nodule produces a suitable niche that facilitates bacterial colonization and the development of an N-fixing nodule.

Evolution of Plant NLRs: From Natural History to Precise Modifications

Nucleotide-binding leucine-rich repeat receptors (NLRs) monitor the plant intracellular environment for signs of pathogen infection. Several mechanisms of NLR-mediated immunity arose independently across multiple species. These include the functional specialization of NLRs into sensors and helpers, the independent emergence of direct and indirect recognition within NLR subfamilies, the regulation of NLRs by small RNAs, and the formation of NLR networks. Understanding the evolutionary history of NLRs can shed light on both the origin of pathogen recognition and the common constraints on the plant immune system. Attempts to engineer disease resistance have been sparse and rarely informed by evolutionary knowledge. In this review, we discuss the evolution of NLRs, give an overview of previous engineering attempts, and propose how to use evolutionary knowledge to advance future research in the generation of novel disease-recognition capabilities.

Genome-wide identification of miRNAs and their targets during early somatic embryogenesis in Dimocarpus longan Lour

miRNAs are endogenous regulatory factors that play pivotal roles in post-transcriptional regulation. However, their specific roles in early somatic embryogenesis (SE) remain unclear. Study of the SE system is fundamental for clarifying the molecular mechanisms in Dimocarpus longan. We identified 289 known miRNAs from 106 different miRNA families and 1087 novel miRNAs during early longan SE, including embryogenic callus (EC), incomplete pro-embryogenic culture (ICpEC), globular embryo (GE), and non-embryogenic callus (NEC). The abundances of known miRNAs were concentrated in GE. The differentially expression (DE) miRNAs showed five expression patterns during early SE. Largely miRNAs were expressed highly and specially in EC, ICpEC, and GE, respectively. Some miRNAs and putative target genes were enriched in lignin metabolism. Most potential targets were related to the pathways of plant hormone signal transduction, alternative splicing, tyrosine metabolism and sulfur metabolism in early longan SE. The regulatory relationships between dlo-miR166a-3p and DlHD-zip8, dlo-miR397a and DlLAC7, dlo-miR408-3p and DlLAC12 were confirmed by RNA ligase-mediated rapid amplification of cDNA ends. The expression patterns of eight DE miRNAs detected by qRT-PCR were consistent with RNA-seq. Finally, the miRNA regulatory network in early SE was constructed, which provided new insight into molecular mechanism of early SE in longan.

A Brassica miRNA Regulates Plant Growth and Immunity through Distinct Modes of Action

In plants, high disease resistance often results in a reduction of yield. Therefore, breeding crops with balanced yield and disease resistance has become a major challenge. Recently, microRNA (miRNA)-mediated R gene turnover has been shown to be a protective mechanism used by plants to prevent autoimmunity in the absence of pathogens. However, whether these miRNAs play a role in plant growth and how miRNA-mediated R gene turnover responds to pathogen infection have rarely been explored. Here, we found that a Brassica miRNA, miR1885, targets both an immune receptor gene and a development-related gene for negative regulation through distinct modes of action. MiR1885 directly silences the TIR-NBS-LRR class of R gene BraTNL1 but represses the expression of the photosynthesis-related gene BraCP24 by targeting the Trans-Acting Silencing (TAS) gene BraTIR1 for trans-acting small interfering RNAs (tasiRNAs)-mediated silencing. We found that, under natural conditions, miR1885 was kept at low levels to maintain normal development and basal immunity but peaked during the floral transition to promote flowering. Interestingly, upon Turnip mosaic virus (TuMV) infection, miR1885-dependent trans-acting silencing of BraCP24 was enhanced to speed up the floral transition, whereas miR1885-mediated R gene turnover was overcome by TuMV-induced BraTNL1 expression, reflecting precise regulation of the arms race between plants and pathogens. Collectively, our results demonstrate that a single Brassica miRNA dynamically regulates both innate immunity and plant growth and responds to viral infection, revealing that Brassica plants have developed a sophisticated mechanism in modulating the interplay between growth, immunity, and pathogen infection.

Potato Virus Y Infection Alters Small RNA Metabolism and Immune Response in Tomato

Potato virus Y (PVY) isolate PVY^C-to induces growth reduction and foliar symptoms in tomato, but new vegetation displays symptom recovery at a later stage. In order to investigate the role of micro(mi)RNA and secondary small(s)RNA-regulated mechanisms in tomato defenses against PVY, we performed sRNA sequencing from healthy and PVY^C-to infected tomato plants at 21 and 30 days post-inoculation (dpi). A total of 792 miRNA sequences were obtained, among which were 123 canonical miRNA sequences, many isomiR variants, and 30 novel miRNAs. MiRNAs were mostly overexpressed in infected vs. healthy plants, whereas only a few miRNAs were underexpressed. Increased accumulation of isomiRs was correlated with viral infection. Among miRNA targets, enriched functional categories included resistance (R) gene families, transcription and hormone factors, and RNA silencing genes. Several 22-nt miRNAs were shown to target R genes and trigger the production of 21-nt phased sRNAs (phasiRNAs). Next, 500 phasiRNA-generating loci were identified, and were shown to be mostly active in PVY-infected tissues and at 21 dpi. These data demonstrate that sRNA-regulated host responses, encompassing miRNA alteration, diversification within miRNA families, and phasiRNA accumulation, regulate R and disease-responsive genes. The dynamic regulation of miRNAs and secondary sRNAs over time suggests a functional role of sRNA-mediated defenses in the recovery phenotype.

Identification and analysis of oxygen responsive microRNAs in the root of wild tomato (S. habrochaites)

BACKGROUND

MicroRNA (miRNA) are key players in regulating expression of target genes at post-transcriptional level. A number of miRNAs are implicated in modulating tolerance to various abiotic stresses. Waterlogging is an abiotic stress that deters plant growth and productivity by hypoxia. Dozens of reports mention about the miRNAs expressed in response to waterlogging and hypoxia. Despite the fact that tomato is a model vegetable but waterlogging sensitive crop, the role of miRNAs in hypoxia tolerance is poorly understood in tomato.

RESULTS

In this study, we investigated the differentially expressed miRNAs between hypoxia-treated and untreated wild tomato root by using high-throughput sequencing technology. A total of 33 known miRNAs were lowly expressed, whereas only 3 miRNAs showed higher expression in hypoxia-treated wild tomato root compared with untreated wild tomato root. Then two conserved and lowly expressed miRNAs, miR171 and miR390, were deactivated by Short Tandem Target Mimic (STTM) technology in Arabidopsis. As the results, the number and length of lateral roots were more in STTM171 and STTM390 transgenic lines compared with that of wild type plant, which partly phenocopy the increase root number and shortening the root length in hypoxia-treated wild tomato root.

CONCLUSIONS

The differentially expressed miRNAs between hypoxia-treated wild tomato and control root, which contribute to the auxin homeostasis, morphologic change, and stress response, might result in reduction in the biomass and length of the root in hypoxiated conditions.

Genome-wide comparative analysis in Solanaceous species reveals evolution of microRNAs targeting defense genes in Capsicum spp

MicroRNAs (miRNAs) play roles in various biological processes in plants including growth, development, and disease resistance. Previous studies revealed that some plant miRNAs produce secondary small interfering RNAs (siRNAs) such as phased, secondary siRNAs (phasiRNAs), and they regulate a cascade of gene expression. We performed a genome-wide comparative analysis of miRNAs in Solanaceous species (pepper, tomato, and potato), from an evolutionary perspective. Microsynteny of miRNAs was analysed based on the genomic loci and their flanking genes and most of the well-conserved miRNA genes maintained microsynteny in Solanaceae. We identified target genes of the miRNAs via degradome analysis and found that several miRNAs target many genes encoding nucleotide-binding leucine-rich repeat (NLR) or receptor-like proteins (RLPs), which are known to be major players in defense responses. In addition, disease-resistance-associated miRNAs trigger phasiRNA production in pepper, indicating amplification of the regulation of disease-resistance gene families. Among these, miR-n033a-3p, whose target NLRs have been duplicated in pepper, targets more NLRs belonging to specific subgroup in pepper than those in potato. miRNAs targeting resistance genes might have evolved to regulate numerous targets in Solanaceae, following expansion of target resistance genes. This study provides an insight into evolutionary relationship between miRNAs and their target defense genes in plants.

Enhanced resistance to bacterial and oomycete pathogens by short tandem target mimic RNAs in tomato

Plants, like animals, have complex disease resistance systems in which receptors interact directly or indirectly with effectors of disease produced by pests and pathogens. To minimize the fitness cost of these systems to the plant, there are miRNAs that target the mRNAs of a family of receptor proteins required for disease resistance. Target site mimics of these miRNAs confer enhanced quantitative resistance in tomato against an oomycete and a bacterium. These findings are consistent with a role of the receptor proteins in quantitative disease resistance and show how blocking these miRNAs could be a useful approach in crop protection.

Nucleotide binding site leucine-rich repeat (NLR) proteins of the plant innate immune system are negatively regulated by the miR482/2118 family miRNAs that are in a distinct 22-nt class of miRNAs with a double mode of action. First, they cleave the target RNA, as with the canonical 21-nt miRNAs, and second, they trigger secondary siRNA production using the target RNA as a template. Here, we address the extent to which the miR482/2118 family affects expression of NLR mRNAs and disease resistance. We show that structural differences of miR482/2118 family members in tomato (Solanum lycopersicum) are functionally significant. The predicted target of the miR482 subfamily is a conserved motif in multiple NLR mRNAs, whereas for miR2118b, it is a noncoding RNA target formed by rearrangement of several different NLR genes. From RNA sequencing and degradome data in lines expressing short tandem target mimic (STTM) RNAs of miR482/2118, we confirm the different targets of these miRNAs. The effect on NLR mRNA accumulation is slight, but nevertheless, the tomato STTM lines display enhanced resistance to infection with the oomycete and bacterial pathogens. These data implicate an RNA cascade of miRNAs and secondary siRNAs in the regulation of NLR RNAs and show that the encoded NLR proteins have a role in quantitative disease resistance in addition to dominant gene resistance that has been well characterized elsewhere. We also illustrate the use of STTM RNA in a biotechnological approach for enhancing quantitative disease resistance in highly bred cultivars.

Molecular mechanisms that limit the costs of NLR-mediated resistance in plants

Crop diseases cause significant yield losses, and the use of resistant cultivars can effectively mitigate these losses and control many plant diseases. Most plant resistance (R) genes encode immune receptors composed of nucleotide-binding and leucine-rich repeat (NLR) domains. These proteins mediate the specific recognition of pathogen avirulence effectors to induce defence responses. However, NLR-triggered immunity can be associated with a reduction in growth and yield, so-called 'fitness costs'. Recent data have shown that plants use an elaborate interplay of different mechanisms to control NLR gene transcript levels, as well as NLR protein abundance and activity, to avoid the associated cost of resistance in the absence of a pathogen. In this review, we discuss the different levels of NLR regulation (transcriptional, post-transcriptional and at the protein level). We address the apparent need for plants to maintain diverse modes of regulation. A recent model suggesting an equilibrium 'ON/OFF state' of NLR proteins, in the absence of a pathogen, provides the context for our discussion.

Evolutionary Divergence of TNL Disease-Resistant Proteins in Soybean (Glycine max) and Common Bean (Phaseolus vulgaris)

Disease-resistant genes (R genes) encode proteins that are involved in protecting plants from their pathogens and pests. Availability of complete genome sequences from soybean and common bean allowed us to perform a genome-wide identification and analysis of the Toll interleukin-1 receptor-like nucleotide-binding site leucine-rich repeat (TNL) proteins. Hidden Markov model (HMM) profiling of all protein sequences resulted in the identification of 117 and 77 regular TNL genes in soybean and common bean, respectively. We also identified TNL gene homologs with unique domains, and signal peptides as well as nuclear localization signals. The TNL genes in soybean formed 28 clusters located on 10 of the 20 chromosomes, with the majority found on chromosome 3, 6 and 16. Similarly, the TNL genes in common bean formed 14 clusters located on five of the 11 chromosomes, with the majority found on chromosome 10. Phylogenetic analyses of the TNL genes from Arabidopsis, soybean and common bean revealed less divergence within legumes relative to the divergence between legumes and Arabidopsis. Syntenic blocks were found between chromosomes Pv10 and Gm03, Pv07 and Gm10, as well as Pv01 and Gm14. The gene expression data revealed basal level expression and tissue specificity, while analysis of available microRNA data showed 37 predicted microRNA families involved in targeting the identified TNL genes in soybean and common bean.

Biogenesis of a 22-nt microRNA in Phaseoleae species by precursor-programmed uridylation

Phased, secondary siRNAs (phasiRNAs) represent a class of small RNAs in plants generated via distinct biogenesis pathways, predominantly dependent on the activity of 22-nt miRNAs. Most 22-nt miRNAs are processed by DCL1 from miRNA precursors containing an asymmetric bulge, yielding a 22/21-nt miRNA/miRNA* duplex. Here we show that miR1510, a soybean miRNA capable of triggering phasiRNA production from numerous nucleotide-binding leucine-rich repeat (NB-LRRs), previously described as 21 nt in its mature form, primarily accumulates as a 22-nt isoform via monouridylation. We demonstrate that, in Arabidopsis, this uridylation is performed by HESO1. Biochemical experiments showed that the 3' terminus of miR1510 is only partially 2'-O-methylated because of the terminal mispairing in the miR1510/miR1510* duplex that inhibits HEN1 activity in soybean. miR1510 emerged in the Phaseoleae ∼41-42 million years ago with a conserved precursor structure yielding a 22-nt monouridylated form, yet a variant in mung bean is processed directly in a 22-nt mature form. This analysis of miR1510 yields two observations: (i) plants can utilize postprocessing modification to generate abundant 22-nt miRNA isoforms to more efficiently regulate target mRNA abundances; and (ii) comparative analysis demonstrates an example of selective optimization of precursor processing of a young plant miRNA.

Biogenesis of a 22-nt microRNA in Phaseoleae species by precursor-programmed uridylation

Phased, secondary siRNAs (phasiRNAs) represent a class of small RNAs in plants generated via distinct biogenesis pathways, predominantly dependent on the activity of 22-nt miRNAs. Most 22-nt miRNAs are processed by DCL1 from miRNA precursors containing an asymmetric bulge, yielding a 22/21-nt miRNA/miRNA* duplex. Here we show that miR1510, a soybean miRNA capable of triggering phasiRNA production from numerous nucleotide-binding leucine-rich repeat (NB-LRRs), previously described as 21 nt in its mature form, primarily accumulates as a 22-nt isoform via monouridylation. We demonstrate that, in Arabidopsis, this uridylation is performed by HESO1. Biochemical experiments showed that the 3' terminus of miR1510 is only partially 2'-O-methylated because of the terminal mispairing in the miR1510/miR1510* duplex that inhibits HEN1 activity in soybean. miR1510 emerged in the Phaseoleae ∼41-42 million years ago with a conserved precursor structure yielding a 22-nt monouridylated form, yet a variant in mung bean is processed directly in a 22-nt mature form. This analysis of miR1510 yields two observations: (i) plants can utilize postprocessing modification to generate abundant 22-nt miRNA isoforms to more efficiently regulate target mRNA abundances; and (ii) comparative analysis demonstrates an example of selective optimization of precursor processing of a young plant miRNA.

miRVIT: A Novel miRNA Database and Its Application to Uncover Vitis Responses to Flavescence dorée Infection

Micro(mi)RNAs play crucial roles in plant developmental processes and in defense responses to biotic and abiotic stresses. In the last years, many works on small RNAs in grapevine (Vitis spp.) were published, and several conserved and putative novel grapevine-specific miRNAs were identified. In order to reorganize the high quantity of available data, we produced "miRVIT," the first database of all novel grapevine miRNA candidates characterized so far, and still not deposited in miRBase. To this aim, each miRNA accession was renamed, repositioned in the last version of the grapevine genome, and compared with all the novel and conserved miRNAs detected in grapevine. Conserved and novel miRNAs cataloged in miRVIT were then used for analyzing Vitis vinifera plants infected by Flavescence dorée (FD), one of the most severe phytoplasma diseases affecting grapevine. The analysis of small RNAs from healthy, recovered (plants showing spontaneous and stable remission of symptoms), and FD-infected "Barbera" grapevines showed that FD altered the expression profiles of several miRNAs, including those involved in cell development and photosynthesis, jasmonate signaling, and disease resistance response. The application of miRVIT in a biological context confirmed the effectiveness of the followed approach, especially for the identification of novel miRNA candidates in grapevine. miRVIT database is available at http://mirvit.ipsp.cnr.it. Highlights: The application of the newly produced database of grapevine novel miRNAs to the analysis of plants infected by Flavescence dorée reveals key roles of miRNAs in photosynthesis and jasmonate signaling.

Bacillus amyloliquefaciens FZB42 represses plant miR846 to induce systemic resistance via a jasmonic acid-dependent signalling pathway

Bacillus amyloliquefaciens FZB42 is a type of plant growth-promoting rhizobacterium (PGPR) which activates induced systemic resistance (ISR) in Arabidopsis. Blocking of the synthesis of cyclic lipopeptides and 2,3-butanediol by FZB42, which have been demonstrated to be involved in the priming of ISR, results in the abolishment of the plant defence responses. To further clarify the ISR activated by PGPRs at the microRNA (miRNA) level, small RNA (sRNA) libraries from Arabidopsis leaves after root irrigation with FZB42, FZB42ΔsfpΔalsS and control were constructed and sequenced. After fold change selection, promoter analysis and target prediction, miR846-5p and miR846-3p from the same precursor were selected as candidate ISR-associated miRNAs. miR846 belongs to the non-conserved miRNAs, specifically exists in Arabidopsis and its function in the plant defence response remains unclear. The disease severity of transgenic Arabidopsis overexpressing miR846 (OEmiR846) or knockdown miR846 (STTM846) against Pseudomonas syringae DC3000 suggests that the miR846 expression level in Arabidopsis is negatively correlated with disease resistance. Moreover, miR846 in Arabidopsis Col-0 is repressed after methyl jasmonate treatment. In addition, jasmonic acid (JA) signalling-related genes are up-regulated in STTM846, and the stomatal apertures of STTM846 are also less than those in Arabidopsis Col-0 after methyl jasmonate treatment. Furthermore, the disease resistance of STTM846 transgenic Arabidopsis against Pseudomonas syringae pv. tomato DC3000 (Pst DC3000) is blocked by the addition of the JA biosynthetic inhibitor diethyldiethiocarbamic acid (DIECA). Taken together, our results suggest that B. amyloliquefaciens FZB42 inoculation suppresses miR846 expression to induce Arabidopsis systemic resistance via a JA-dependent signalling pathway.

Genomic and epigenomic immunity in common bean: the unusual features of NB-LRR gene family

In plants, a key class of genes comprising most of disease resistance (R) genes encodes Nucleotide-binding leucine-rich repeat (NL) proteins. Access to common bean (Phaseolus vulgaris) genome sequence provides unparalleled insight into the organization and evolution of this large gene family (∼400 NL) in this important crop. As observed in other plant species, most common bean NL are organized in cluster of genes. However, a particularity of common bean is that these clusters are often located in subtelomeric regions close to terminal knobs containing the satellite DNA khipu. Phylogenetically related NL are spread between different chromosome ends, suggesting frequent exchanges between non-homologous chromosomes. NL peculiar location, in proximity to heterochromatic regions, led us to study their DNA methylation status using a whole-genome cytosine methylation map. In common bean, NL genes displayed an unusual body methylation pattern since half of them are methylated in the three contexts, reminiscent of the DNA methylation pattern of repeated sequences. Moreover, 90 NL were also abundantly targeted by 24 nt siRNA, with 90% corresponding to methylated NL genes. This suggests the existence of a transcriptional gene silencing mechanism of NL through the RdDM (RNA-directed DNA methylation) pathway in common bean that has not been described in other plant species.

Variability among Cucurbitaceae species (melon, cucumber and watermelon) in a genomic region containing a cluster of NBS-LRR genes

BACKGROUND

Cucurbitaceae species contain a significantly lower number of genes coding for proteins with similarity to plant resistance genes belonging to the NBS-LRR family than other plant species of similar genome size. A large proportion of these genes are organized in clusters that appear to be hotspots of variability. The genomes of the Cucurbitaceae species measured until now are intermediate in size (between 350 and 450 Mb) and they apparently have not undergone any genome duplications beside those at the origin of eudicots. The cluster containing the largest number of NBS-LRR genes has previously been analyzed in melon and related species and showed a high degree of interspecific and intraspecific variability. It was of interest to study whether similar behavior occurred in other cluster of the same family of genes.

RESULTS

The cluster of NBS-LRR genes located in melon chromosome 9 was analyzed and compared with the syntenic regions in other cucurbit genomes. This is the second cluster in number within this species and it contains nine sequences with a NBS-LRR annotation including two genes, Fom1 and Prv, providing resistance against Fusarium and Ppapaya ring-spot virus (PRSV). The variability within the melon species appears to consist essentially of single nucleotide polymorphisms. Clusters of similar genes are present in the syntenic regions of the two species of Cucurbitaceae that were sequenced, cucumber and watermelon. Most of the genes in the syntenic clusters can be aligned between species and a hypothesis of generation of the cluster is proposed. The number of genes in the watermelon cluster is similar to that in melon while a higher number of genes (12) is present in cucumber, a species with a smaller genome than melon. After comparing genome resequencing data of 115 cucumber varieties, deletion of a group of genes is observed in a group of varieties of Indian origin.

CONCLUSIONS

Clusters of genes coding for NBS-LRR proteins in cucurbits appear to have specific variability in different regions of the genome and between different species. This observation is in favour of considering that the adaptation of plant species to changing environments is based upon the variability that may occur at any location in the genome and that has been produced by specific mechanisms of sequence variation acting on plant genomes. This information could be useful both to understand the evolution of species and for plant breeding.

Integrated RNA-seq and sRNA-seq analysis reveals miRNA effects on secondary metabolism in Solanum tuberosum L

Light is a major environmental factor that affects metabolic pathways and stimulates the production of secondary metabolites in potato. However, adaptive changes in potato metabolic pathways and physiological functions triggered by light are partly explained by gene expression changes. Regulation of secondary metabolic pathways in potato has been extensively studied at transcriptional level, but little is known about the mechanisms of post-transcriptional regulation by miRNAs. To identify light-responsive miRNAs/mRNAs and construct putative metabolism pathways regulated by the miRNA-mRNA pairs, an integrated omics (sRNAome and transcriptome) analysis was performed to potato under light stimulus. A total of 31 and 48 miRNAs were identified to be differentially expressed in the leaves and tubers, respectively. Among the DEGs, 1353 genes in the leaves and 1841 genes in the tubers were upregulated, while 1595 genes in the leaves and 897 genes in the tubers were downregulated by light. Mapman enrichment analyses showed that genes related to MVA pathway, alkaloids-like, phenylpropanoids, flavonoids, and carotenoids metabolism were significantly upregulated, while genes associated with major CHO metabolism were repressed in the leaves and tubers. Integrated miRNA and mRNA profiles revealed that light-responsive miRNAs are important regulators in alkaloids metabolism, UMP salvage, lipid biosynthesis, and cellulose catabolism. Moreover, several miRNAs may participate in glycoalkaloids metabolism via JA signaling pathway, UDP-glucose biosynthesis and hydroxylation reaction. This study provides a global view of miRNA and mRNA expression profiles in potato response to light, our results suggest that miRNAs might play important roles in secondary metabolic pathways, especially in glycoalkaloid biosynthesis. The findings will enlighten us on the genetic regulation of secondary metabolite pathways and pave the way for future application of genetically engineered potato.

The Diversification of Plant NBS-LRR Defense Genes Directs the Evolution of MicroRNAs That Target Them

High expression of plant nucleotide binding site leucine-rich repeat (NBS-LRR) defense genes is often lethal to plant cells, a phenotype perhaps associated with fitness costs. Plants implement several mechanisms to control the transcript level of NBS-LRR defense genes. As negative transcriptional regulators, diverse miRNAs target NBS-LRRs in eudicots and gymnosperms. To understand the evolutionary benefits of this miRNA-NBS-LRR regulatory system, we investigated the NBS-LRRs of 70 land plants, coupling this analysis with extensive small RNA data. A tight association between the diversity of NBS-LRRs and miRNAs was found. The miRNAs typically target highly duplicated NBS-LRRs In comparison, families of heterogeneous NBS-LRRs were rarely targeted by miRNAs in Poaceae and Brassicaceae genomes. We observed that duplicated NBS-LRRs from different gene families periodically gave birth to new miRNAs. Most of these newly emerged miRNAs target the same conserved, encoded protein motif of NBS-LRRs, consistent with a model of convergent evolution for these miRNAs. By assessing the interactions between miRNAs and NBS-LRRs, we found nucleotide diversity in the wobble position of the codons in the target site drives the diversification of miRNAs. Taken together, we propose a co-evolutionary model of plant NBS-LRRs and miRNAs hypothesizing how plants balance the benefits and costs of NBS-LRR defense genes.

A network analysis of miRNA mediated gene regulation of rice: crosstalk among biological processes

To understand the network architecture of miRNA mediated regulations at the genomic and functional levels of rice, we have made an unambiguous annotation of the experimentally verified miRNAs, predicted their targets and the possible biological functions they can affect. Some functions, namely translational and protein modifications and photosynthesis are targeted by higher percentage of miRNA. Using transformation procedures, we constructed a genome scale miRNA-miRNA functional synergistic network (MFSN). The analysis of MFSN modules help to identify miRNAs co-regulating target genes having several interrelated biological processes. Some of these target genes are also co-expressed under particular conditions. For example, the genes co-expressed under drought conditions as well as those targeted by miRNAs present in a MFSN module have interdependent biological processes namely, photosynthesis, cell-wall biogenesis, root development and xylan synthesis. The stress-induced miRNAs and their distributions, and the presence of transcription factors in the target set of MFSN modules were also analyzed.

Phylogenetic tree-informed microRNAome analysis uncovers conserved and lineage-specific miRNAs in Camellia during floral organ development

In plants, miRNAs are endogenous small RNAs derived from single-stranded precursors with hairpin structures. The evolution of miRNAs and their targets represents one of the most dynamic circuits directing gene expression, which may play fundamental roles in shaping the development of distinct plant organs. Here we performed high-throughput small RNA sequencing in five organ types of Camellia azalea to capture the spatial profile of small non-coding RNA. In total we obtained >227 million high-quality reads and identified 175 miRNAs with mature and precursor sequences. We aligned the miRNAs to known miRNA databases and revealed some conserved as well as 'newly evolved' miRNA genes. Twelve miRNAs were identified to be specific in the genus Camellia, supporting the lineage-specific manner of expansion of 'young' miRNAs. Through differential expression analysis, we showed that many miRNAs were preferentially abundant in certain organ types. Moreover, hierarchical clustering analysis revealed distinctive expression patterns of tissue-specific miRNAs. Gene Ontology enrichment analysis of targets of stamen- and carpel-specific miRNA subclusters showed that miRNA-target regulatory circuits were involved in many important biological processes, enabling their proper specification and organogenesis, such as 'DNA integration' and 'fruit development'. Further, quantitative PCR of key miRNAs and their target genes revealed anti-correlated patterns, and uncovered the functions of key miRNA-target pairs in different floral organs. Taken together, this work yielded valuable information on miRNA-target regulation in the control of floral organ development and sheds light on the evolution of lineage-specific miRNAs in Camellia.

Development of molecular markers linked to disease resistance genes in common bean based on whole genome sequence

Common bean (Phaseolus vulgaris) is the most important grain legume for direct human consumption in the world, particularly in developing countries where it constitutes the main source of protein. Unfortunately, common bean yield stability is constrained by a number of pests and diseases. As use of resistant genotypes is the most economic and ecologically safe means for controlling plant diseases, efforts have been made to genetically characterize resistance genes (R genes) in common bean. Despite its agronomic importance, genomic resources available in common bean were limited until the recent sequencing of common bean genome (Andean genotype G19833). Besides allowing the annotation of Nucleotide Binding-Leucine Rich Repeat (NB-LRR) encoding gene family, which is the prevalent class of disease R genes in plants, access to the whole genome sequence of common bean can be of great help for intense selection to increase the overall efficiency of crop improvement programs using marker-assisted selection (MAS). This review presents the state of the art of common bean NB-LRR gene clusters, their peculiar location in subtelomeres and correlation with genetically characterized monogenic R genes, as well as how the availability of the whole genome sequence can boost the development of molecular markers for MAS.

Evolutionarily Dynamic, but Robust, Targeting of Resistance Genes by the miR482/2118 Gene Family in the Solanaceae

Plants are exposed to pathogens around the clock. A common resistance response in plants upon pathogen detection is localized cell death. Given the irreversible nature of this response, multiple layers of negative regulation are present to prevent the untimely or misexpression of resistance genes. One layer of negative regulation is provided by a recently discovered microRNA (miRNA) gene family, miR482/2118. This family targets the transcripts of resistance genes in plants. We investigated the evolutionary history and specificity of this miRNA gene family within the Solanaceae. This plant family includes many important crop species, providing a set of well-defined resistance gene repertoires. Across 14 species from the Solanaceae, we identified eight distinct miR482/2118 gene family members. Our studies show conservation of miRNA type and number in the group of wild tomatoes and, to a lesser extent, throughout the Solanaceae. The eight orthologous miRNA gene clusters evolved under different evolutionary constraints, allowing for individual subfunctionalization of the miRNAs. Despite differences in the predicted targeting behavior of each miRNA, the miRNA-R-gene network is robust due to its high degree of interconnectivity and redundant targeting. Our data suggest that the miR482/2118 gene family acts as an evolutionary buffer for R-gene sequence diversity.

Interspecific and intraspecific gene variability in a 1-Mb region containing the highest density of NBS-LRR genes found in the melon genome

BACKGROUND

Plant NBS-LRR -resistance genes tend to be found in clusters, which have been shown to be hot spots of genome variability. In melon, half of the 81 predicted NBS-LRR genes group in nine clusters, and a 1 Mb region on linkage group V contains the highest density of R-genes and presence/absence gene polymorphisms found in the melon genome. This region is known to contain the locus of Vat, an agronomically important gene that confers resistance to aphids. However, the presence of duplications makes the sequencing and annotation of R-gene clusters difficult, usually resulting in multi-gapped sequences with higher than average errors.

RESULTS

A 1-Mb sequence that contains the largest NBS-LRR gene cluster found in melon was improved using a strategy that combines Illumina paired-end mapping and PCR-based gap closing. Unknown sequence was decreased by 70% while about 3,000 SNPs and small indels were corrected. As a result, the annotations of 18 of a total of 23 NBS-LRR genes found in this region were modified, including additional coding sequences, amino acid changes, correction of splicing boundaries, or fussion of ORFs in common transcription units. A phylogeny analysis of the R-genes and their comparison with syntenic sequences in other cucurbits point to a pattern of local gene amplifications since the diversification of cucurbits from other families, and through speciation within the family. A candidate Vat gene is proposed based on the sequence similarity between a reported Vat gene from a Korean melon cultivar and a sequence fragment previously absent in the unrefined sequence.

CONCLUSIONS

A sequence refinement strategy allowed substantial improvement of a 1 Mb fragment of the melon genome and the re-annotation of the largest cluster of NBS-LRR gene homologues found in melon. Analysis of the cluster revealed that resistance genes have been produced by sequence duplication in adjacent genome locations since the divergence of cucurbits from other close families, and through the process of speciation within the family a candidate Vat gene was also identified using sequence previously unavailable, which demonstrates the advantages of genome assembly refinements when analyzing complex regions such as those containing clusters of highly similar genes.