MicroRNA-target gene responses to root knot nematode (Meloidogyne incognita) infection in cotton (Gossypium hirsutum L.)

MicroRNAs (miRNAs) are a large class of small regulatory RNA molecules, however no study has been performed to elucidate the role of miRNAs in cotton (Gossypium hirsutum) response to the root knot nematode (RKN, Meloidogyne incognita) infection. We selected 28 miRNAs and 8 miRNA target genes to investigate the miRNA-target gene response to M. incognita infection. Our results show that RKN infection significantly affected the expression of several miRNAs and their targeted genes. After 10 days of RKN infection, expression fold changes on miRNA expressions ranged from down-regulated by 33% to upregulated by 406%; meanwhile the expression levels of miRNA target genes were 45.8% to 231%. Three miRNA-target pairs, miR159-MYB, miR319-TCP4 and miR167-ARF8, showed inverse expression patterns between gene targets and their corresponded miRNAs, suggesting miRNA-mediated gene regulation in cotton roots in response to RKN infection.     

Phylogeny and Molecular Evolution of miR820 and miR396 microRNA Families in <Emphasis Type="Italic">Oryza</Emphasis> AA Genomes

The phylogeny and evolution of the microRNA families, miR820 and miR396, was analysed across the AA genomes of the Oryza species, the close relatives of domesticated rice. A highly dynamic evolution of the miR820 family was revealed. The number of copies of MIR820 genes, their chromosomal location and the mature microRNA sequence varied greatly with a total of 16 novel miR820 variants being identified. The phylogeny of pre-MIR820 sequences revealed that MIR820 genes of recently evolved Oryza AA genomes may have derived from sequence divergence of one or a few ancestral genes found in wild Australian perennial rice populations, Taxon B (jpn2)-MIR820 genes. Genomic scale duplication played an important role in the evolution of some of the miR396 family genes in AA genome Oryza species. miR396 family contained a MIR396 gene cluster (MIR396a and MIR396c) which was conserved across the cereal genomes. Nucleotide diversity analysis at these two MIR396 loci revealed that domesticated rice has retained less than 10% of the total diversity present in wild species. In contrast, the nucleotide sequence of four MIR396 loci remained almost conserved across domesticated and wild rices, indicating that they were under extreme functional constraint and may be involved in regulating some fundamental processes which are important both for wild and domesticated rices. Expression analysis demonstrated that miR820 variants were expressed in O. glaberrima O. barthi and O. longistaminata genome. These findings pose new challenges to explain the possible role of miR820 variants identified.     

Analysis of gradient-like expression of miR167 in Arabidopsis thaliana embryonic tissue

MicroRNAs (miRNAs) constitute a class of small non-coding endogenous RNAs that regulate gene expression by targeting mRNAs for degradation or translational repression. MiRNAs are intensively investigated and they have been found to be a pivotal component of developmental regulation processes. Recent studies showed the non-cell autonomous function of several miRNAs. We analyzed the accumulation pattern of selected miRNAs in Arabidopsis thaliana embryonic tissues. The majority of the investigated miRNAs showed uniform accumulation across the embryo suggesting their possible role at this developmental stage. In the case of miR167 however, we detected a gradient-like expression profile which in earlier studies has been considered to be the hallmark of the non-cell autonomous activity of miRNAs. Using reporter assay we analyzed the expression patterns of the four MIR167 precursor genes. We found that two of the precursor genes, MIR167A and MIR167B, also showed an overlapping gradient-like expression patterns in the embryo. These data indicate that in addition to non-cell autonomous activity of some miRNAs, the gradient-like expression patterns can be generated also by the specific expression characteristic of miRNA precursor genes.     

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.     

The evolutionary-developmental analysis of plant microRNAs

MicroRNAs (miRNAs) control many important aspects of plant development, suggesting these molecules may also have played key roles in the evolution of developmental processes in plants. However, evolutionary-developmental (evo-devo) studies of miRNAs have been held back by technical difficulties in gene identification. To help solve this problem, we have developed a two-step procedure for the efficient identification of miRNA genes in any plant species. As a test case, we have studied the evolution of the MIR164 family in the angiosperms. We have identified novel MIR164 genes in three species occupying key phylogenetic positions and used these, together with published sequence data, to partially reconstruct the evolution of the MIR164 family since the last common ancestor of the extant flowering plants. We use our evolutionary reconstruction to discuss potential roles for MIR164 genes in the evolution of leaf shape and carpel closure in the angiosperms. The techniques we describe may be applied to any miRNA family and should thus enable plant evo-devo to begin to investigate the contributions miRNAs have made to the evolution of plant development.     

Redundant and specific roles of individual MIR172 genes in plant development

Evolutionarily conserved microRNAs (miRNAs) usually have high copy numbers in the genome. The redundant and specific roles of each member of a multimember miRNA gene family are poorly understood. Previous studies have shown that the miR156-SPL-miR172 axis constitutes a signaling cascade in regulating plant developmental transitions. Here, we report the feasibility and utility of CRISPR-Cas9 technology to investigate the functions of all 5 MIR172 family members in Arabidopsis. We show that an Arabidopsis plant devoid of miR172 is viable, although it displays pleiotropic morphological defects. MIR172 family members exhibit distinct expression pattern and exert functional specificity in regulating meristem size, trichome initiation, stem elongation, shoot branching, and floral competence. In particular, we find that the miR156-SPL-miR172 cascade is bifurcated into specific flowering responses by matching pairs of coexpressed SPL and MIR172 genes in different tissues. Our results thus highlight the spatiotemporal changes in gene expression that underlie evolutionary novelties of a miRNA gene family in nature. The expansion of MIR172 genes in the Arabidopsis genome provides molecular substrates for the integration of diverse floral inductive cues, which ensures that plants flower at the optimal time to maximize seed yields.

This study uses CRISPR-Cas9 technology to investigate the functions of all five miR172 genes in Arabidopsis, finding that miRNA172 family members exhibit distinct expression pattern and exert functional specificity in regulating meristem size, trichome initiation, stem elongation, shoot branching and floral competence.     

The Regulatory Roles of microRNAs and Associated Target Genes during Early Somatic Embryogenesis in <i>Liriodendron Sino-Americanum</i>

Somatic cells respond to considerable stress, and go through a series of phytohormone pathways, then forming an embryo. The developmental process is recorded as somatic embryogenesis (SE). One of the key components regulating SE are the microRNAs (miRNAs). Despite previous studies, it is still not clear exactly how miRNAs exert their function of regulating targets during conditionally activated early SE. Here, we use Liriodendron sino-americanum as a model system and perform a combined analysis of microfluidic chips and degradome sequencing to study this process. We identified a total of 386 conserved miRNAs and 153 novel miRNAs during early SE. According to the ANOVA test, 239 miRNAs showed 12 distinct expression patterns. Through degradome sequencing, 419 targets and 198 targets were identified for 136 known miRNAs and 37 novel miRNAs, respectively. Gene Ontology (GO) and metabolism pathway enrichment analysis revealed that these targets were significantly involved in oxidation-reduction processes, calmodulin-mediated signal transduction pathways and carbohydrate metabolism. The genes that were related to stress responses, phytohormone pathways and plant metabolism were identified within the targets of miR319, miR395, miR408, miR472, miR482, miR390, miR2055, miR156, miR157, miR171, miR396, miR397, miR529, miR535 and miR159. According to promoter analysis, various cis-acting elements related to plant growth and development, phytohormones response and stress response were present in the promoter of the miRNAs. The differential expression patterns of 11 miRNA-target modules were confirmed by real-time quantitative PCR. The study demonstrated that the miRNA plays an important role in the early SE process by regulating its target and then participating in carbohydrate metabolism and stress response. It also provided a valuable resource for further research in determining the genetic mechanism of SE, and then facilitating breeding programs on plants.     

Identification and expression analysis of miRNAs in germination and seedling growth of Tibetan hulless barley

Germination and seedling growth are crucial for plant development and agricultural production. While, the regulatory mechanisms during these processes in Tibetan hulless barley (Hordeum vulgare L. var. nudum) are not well understood. Given the regulatory roles of microRNAs (miRNAs) in crop plants and the irreplaceability of barley in the highland area of China, we herein presented a genome-wide survey of miRNAs to reveal a potential regulatory network in the early developmental stages of two Tibetan hulless barleys, from which a total of 156 miRNAs was identified including 35 known and 121 novel ones. Six of the identified novel miRNAs were further experimentally validated. According to the evolutionary analysis, miR156, miR166, miR168, and miR171 were conserved across Tibetan hulless barleys and eight other seed plants. Expression profiles of ten known miRNAs showed that they were involved in phytohormone signaling, carbohydrate and lipid metabolism, as well as juvenile-adult transition during barley development. Moreover, a total of 1280 genes targeted by 101 miRNAs were predicted from both barley libraries. Three genes (PLN03212, MATE eukaryotic, and GRAS) were validated via the RNA ligase-mediated 5′-rapid amplification of cDNA ends (RLM-5' RACE) to be the targets of hvu-miR159a, hvu-miR166a, and hvu-miR171-3p, respectively. Based on the Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis of putative targets, the most abundant pathways were related to “metabolism”.These results revealed that miRNA-target pairs participating in the regulation of multigene expression and the embryonic development of Tibetan hulless barleys were controlled by complex mechanisms involving the concordant expression of different miRNAs and feedback loops among miRNAs as well as their targets. The study provides insight into the regulatory network of barley miRNAs for better understanding of miRNA functions during germination and seedling growth.     

miR482 Regulation of NBS-LRR Defense Genes during Fungal Pathogen Infection in Cotton

In this study, we characterized the miR482 family in cotton using existing small RNA datasets and the recently released draft genome sequence of Gossypium raimondii, a diploid cotton species whose progenitor is the putative contributor of the D_t (representing the D genome of tetraploid) genome of the cultivated tetraploid cotton species G. hirsutum and G. barbadense. Of the three ghr-miR482 members reported in G. hirsutum, ghr-miR482a has no homolog in G. raimondii, ghr-miR482b and ghr-miR482c each has a single homolog in G. raimondii. Gra-miR482d has five homologous loci (gra-miR482d, f-i) in G. raimondii and also exists in G. hirsutum (ghr-miR482d). A variant, miR482.2 that is a homolog of miR2118 in other species, is produced from several GHR-MIR482 loci in G. hirsutum. Approximately 12% of the G. raimondii NBS-LRR genes were predicted targets of various members of the gra-miR482 family. Based on the rationale that the regulatory relationship between miR482 and NBS-LRR genes will be conserved in G. raimondii and G. hirsutum, we investigated this relationship using G. hirsutum miR482 and G. raimondii NBS-LRR genes, which are not currently available in G. hirsutum. Ghr-miR482/miR482.2-mediated cleavage was confirmed for three of the four NBS-LRR genes analysed. As in tomato, miR482-mediated cleavage of NBS-LRR genes triggered production of phased secondary small RNAs in cotton. In seedlings of the susceptible cultivar Sicot71 (G. hirsutum) infected with the fungal pathogen Verticillium dahliae, the expression levels of ghr-miR482b/miR482b.2, ghr-miR482c and ghr-miR482d.2 were down-regulated, and several NBS-LRR targets of ghr-miR482c and ghr-miR482d were up-regulated. These results imply that, like tomato plants infected with viruses or bacteria, cotton plants are able to induce expression of NBS-LRR defence genes by suppression of the miRNA-mediated gene silencing pathway upon fungal pathogen attack.