Small RNA and degradome deep sequencing reveals important roles of microRNAs in cotton (Gossypium hirsutum L.) response to root-knot nematode Meloidogyne incognita infection

Investigation of cotton response to nematode infection will allow us to better understand the cotton immune defense mechanism and design a better biotechnological approach for efficiently managing pest nematodes in cotton. In this study, we firstly treated cotton by root knot nematode (RKN, Meloidogyne incognita) infections, then we employed the high throughput deep sequencing technology to sequence and genome-widely identify all miRNAs in cotton; finally, we analyzed the functions of these miRNAs in cotton response to RKN infections. A total of 266 miRNAs, including 193 known and 73 novel miRNAs, were identified by deep sequencing technology, which belong to 67 conserved and 66 novel miRNA families, respectively. A majority of identified miRNA families only contain one miRNA; however, miR482 family contains 14 members and some others contain 2–13 members. Certain miRNAs were specifically expressed in RKN-infected cotton roots and others were completely inhibited by RKN infection. A total of 50 miRNAs were differentially expressed after RKN infection, in which 28 miRNAs were up-regulated and 22 were inhibited by RKN treatment. Based on degradome sequencing, 87 gene targets were identified to be targeted by 57 miRNAs. These miRNA-targeted genes are involved in the interaction of cotton plants and nematode infection. Based on GO (gene ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) analysis, 466 genes from all 636 miRNA targets were mapped to 6340 GO terms, 181 genes from 228 targets of differentially expressed miRNAs were mapped to 1588 GO terms. The GO terms were then categorized into the three main GO classes: biological processes, cellular components, and molecular functions. The targets of differentially expressed miRNAs were enriched in 43 GO terms, including 22 biological processes, 10 cellular components, and 11 molecular functions (p < 0.05). Many identified processes were associated with organism responses to the environmental stresses, including regulation of nematode larval development, response to nematode, and response to flooding. Our results will enhance the study and application of developing new cotton cultivars for nematode resistance.     

Discovery and profiling of novel and conserved microRNAs during flower development in Carya cathayensis via deep sequencing

Hickory (Carya cathayensis Sarg.) is an economically important woody plant in China, but its long juvenile phase delays yield. MicroRNAs (miRNAs) are critical regulators of genes and important for normal plant development and physiology, including flower development. We used Solexa technology to sequence two small RNA libraries from two floral differentiation stages in hickory to identify miRNAs related to flower development. We identified 39 conserved miRNA sequences from 114 loci belonging to 23 families as well as two novel and ten potential novel miRNAs belonging to nine families. Moreover, 35 conserved miRNA*s and two novel miRNA*s were detected. Twenty miRNA sequences from 49 loci belonging to 11 families were differentially expressed; all were up-regulated at the later stage of flower development in hickory. Quantitative real-time PCR of 12 conserved miRNA sequences, five novel miRNA families, and two novel miRNA*s validated that all were expressed during hickory flower development, and the expression patterns were similar to those detected with Solexa sequencing. Finally, a total of 146 targets of the novel and conserved miRNAs were predicted. This study identified a diverse set of miRNAs that were closely related to hickory flower development and that could help in plant floral induction.     

Identification and characterization of microRNAs and their target genes in Brassica oleracea

The microRNAs are a new class of small non-coding endogenous RNAs with lengths of approximately ~21 nt. MicroRNAs perform their biological function via the degradation of the target mRNAs or by inhibiting protein translation. Until recently, only limited numbers of miRNAs were identified in Brassica oleracea, a vegetable widely cultivated around the world. In present study, 193 potential miRNA candidates were identified from 17 expressed sequence tag (ESTs) and 152 genome survey sequences (GSSs) in B. oleracea. These miRNA candidates were classified into 70 families using a well-defined comparative genome-based computational analysis. Most miRNAs belong to the miRNA169, miR5021, miR156 and miR158 families. Of these, 36 miRNA families are firstly found in Brassica species. Around 1393 B. oleracea genes were predicted as candidate targets of 175 miRNAs. The mutual relationship between miRNAs and the candidate target genes was verified by checking differentially expression levels using quantitative real-time polymerase chain reaction (qRT-PCR) and 5' RLM-RACE analyses. These target genes participate in multiple biological and metabolic processes, including signal transduction, stress response, and plant development. Gene Ontology analysis shows that the 818, 514, and 265 target genes are involved in molecular functions, biological processes, and cellular component respectively. The Kyoto Encyclopedia of Genes and Genomes (KEGG) Pathway enrichment analysis suggests that these miRNAs might regulate 186 metabolic pathways, including those of lipid, energy, starch and sucrose, fatty acid and nitrogen.     

Identification and analyses of miRNA genes in allotetraploid Gossypium hirsutum fiber cells based on the sequenced diploid G. raimondii genome

The plant genome possesses a large number of microRNAs（miRNAs）mainly 21-24 nucleotides in length.They play a vital role in regulation of target gene expression at various stages throughout the whole plant life cycle.Here we sequenced and analyzed～10 million non-coding RNAs（ncRNAs）derived from fiber tissue of the allotetraploid cotton（Gossypium hirsutum）1 days post-anthesis using ncRNA-seq technology.In terms of distinct reads,24 nt ncRNA is by far the dominant species,followed by 21 nt and 23 nt ncRNAs. Using ab initio prediction,we identified and characterized a total of 562 candidate miRNA gene loci on the recently assembled D₅ genome of the diploid cotton G.raimondii.Of all the 562 predicted miRNAs,22 were previously discovered in cotton species and 187 had sequence conservation and homology to homologous miRNAs of other plant species.Nucleotide bias analysis showed that the 9th and 1 st positions were significantly conserved among different types of miRNA genes.Among the 463 putative miRNA target genes,most significant up/down-regulation occurred in 10-20 days post-anthesis,indicating that miRNAs played an important role during the elongation and secondary cell wall synthesis stages of cotton fiber development.The discovery of new miRNA genes will help understand the mechanisms of miRNA generation and regulation in cotton.     

Identification of soybean microRNAs involved in soybean cyst nematode infection by deep sequencing

Soybean cyst nematode (SCN), Heterodera glycines, is the most devastating pathogen of soybean worldwide. MicroRNAs (miRNAs) are a class of small, non-coding RNAs that are known to play important role in plant stress response. However, there are few reports profiling the miRNA expression patterns during pathogen stress. We sequenced four small RNA libraries from two soybean cultivar (Hairbin xiaoheidou, SCN race 3 resistant, Liaodou 10, SCN race 3 susceptible) that grown under un-inoculated and SCN-inoculated soil. Small RNAs were mapped to soybean genome sequence, 364 known soybean miRNA genes were identified in total. In addition, 21 potential miRNA candidates were identified. Comparative analysis of miRNA profiling indicated 101 miRNAs belong to 40 families were SCN-responsive. We also found 20 miRNAs with different express pattern even between two cultivars of the same species. These findings suggest that miRNA paly important role in soybean response to SCN and have important implications for further identification of miRNAs under pathogen stress.