High-throughput sequence analysis of small RNAs in skotomorphogenic seedlings of Brassica rapa ssp. rapa

Skotomorphogenic development is the process by which seedlings adapt to a stressful dark environment. Such metabolic responses to abiotic stresses in plants are known to be regulated in part by microRNAs (miRNAs); however, little is known about the involvement of miRNAs in the regulation of skotomorphogenesis. To identify miRNAs at the genome-wide level in skotomorphogenic seedlings of turnip (Brassica rapa subsp. rapa), an important worldwide root vegetable, we used Solexa sequencing to sequence a small RNA library from seedlings grown in the dark for 4 days. Deep sequencing showed that the small RNAs (sRNAs) were predominantly 21 to 24 nucleotides long. Specifically, 13,319,035 reads produced 359,531 unique sRNAs including rRNA, tRNA, miRNA, small nuclear RNA (snRNA), small nucleolar RNA (snoRNA), and unannotated sRNAs. Sequence analysis identified 96 conserved miRNAs belonging to 36 miRNA families and 576 novel miRNAs. qRT-PCR confirmed that the miRNAs were expressed during skotomorphogenesis similar to the trends shown by the Solexa sequencing results. A total of 2013 potential targets were predicted, and the targets of BrmiR157, BrmiR159 and BrmiR160 were proved to be regulated by miRNA-guided cleavage. These results show that specific regulatory miRNAs are present in skotomorphogenic seedlings of turnip and may play important roles in growth, development, and response to dark environment.     

Diversity and Expression of MicroRNAs in the Filarial Parasite,Brugia malayi

Human filarial parasites infect an estimated 120 million people in 80 countries worldwide causing blindness and the gross disfigurement of limbs and genitals. An understanding of RNA-mediated regulatory pathways in these parasites may open new avenues for treatment. Toward this goal, small RNAs from Brugia malayi adult females, males and microfilariae were cloned for deep-sequencing. From ∼30 million sequencing reads, 145 miRNAs were identified in the B. malayi genome. Some microRNAs were validated using the p19 RNA binding protein and qPCR. B. malayi miRNAs segregate into 99 families each defined by a unique seed sequence. Sixty-one of the miRNA families are highly conserved with homologues in arthropods, vertebrates and helminths. Of those miRNAs not highly conserved, homologues of 20 B. malayi miRNA families were found in vertebrates. Nine B. malayi miRNA families appear to be filarial-specific as orthologues were not found in other organisms. The miR-2 family is the largest in B. malayi with 11 members. Analysis of the sequences shows that six members result from a recent expansion of the family. Library comparisons found that 1/3 of the B. malayi miRNAs are differentially expressed. For example, miR-71 is 5–7X more highly expressed in microfilariae than adults. Studies suggest that in C.elegans, miR-71 may enhance longevity by targeting the DAF-2 pathway. Characterization of B. malayi miRNAs and their targets will enhance our understanding of their regulatory pathways in filariads and aid in the search for novel therapeutics.     

Genome-wide identification of novel microRNAs from genome sequences using computational approach in the mudskipper (<Emphasis Type="Italic">Boleophthalmus pectinirostris</Emphasis>)

MicroRNAs (miRNAs), approximately 22 nucleotides (nt) long, are small, non-coding RNA molecules with important regulatory functions in gene expression. They are mostly conserved among the organisms and this conservation makes them a good source for the identification of novel miRNAs by computational genomic homology. The miRNA repertoire of a major aquaculture species, Boleophthalmus pectinirostris, has been unknown until recently. Currently, the B. pectinirostris whole-genome sequences have been completed, making it more convenient for us to focus on computational prediction for novel miRNA homologs. Following a range of strict filtering criteria, a total of 62 potential miRNAs were identified for the first time; they belong to 39 different miRNA families. All these miRNAs were observed in the stem portion of the stable stem–loop structures. The minimum free energy (MFE) of the predicted miRNAs ranged from ?21.6 to ?62.7 kcal/mol with an average of ?39.2 kcal/mol. The A + U ranged from 32.5 to 69.1% with an average value of 52.2%. The phylogenetic analysis of predicted miRNAs revealed that miR-23a-3p, miR-184-3p, miR-214-5p, and miR-338-3p from B. pectinirostris are evolutionary highly conserved showing more similarity with other fish species. To verify the predicted miRNAs, selected miRNAs representing 16 of the 39 families were confirmed by stem–loop RT-PCR, indicating that the computational approach that we used to identify the miRNAs is a highly efficient and affordable alternative method. Taken together, these findings provide a reference point for further research on miRNAs identification in fish species, meanwhile, our study also will be useful for further insight into biological functions of miRNAs and improved understanding of genome in B. pectinirostris.     

Identification of microRNAs from Plutella xylostella larvae associated with parasitization by Diadegma semiclausum

MicroRNAs (miRNAs) as small non-coding RNAs play important roles in many biological processes such as development, cell signaling and immune response. Small RNA deep sequencing technology provided an opportunity for a thorough survey of miRNAs in a global key pest Plutella xylostella as well as comparative analysis of miRNA expression profile of the insect in association with parasitization by Diadegma semiclausum. Combining the deep sequencing data and bioinformatics, 235 miRNAs were identified from P. xylostella. Differential expression of host cellular miRNAs in response to parasitism was examined by making small RNA libraries from parasitized and naive second instar larvae of P. xylostella. Bantam, miR-276*, miR-10, miR-31 and miR-184 were detected as five most abundant miRNAs in both libraries and 96 miRNAs were identified that were differentially expressed after parasitization. Bantam*, miR-184 and miR-281* were significantly down-regulated and two miRNAs miR-279b and miR-2944b* were highly induced in parasitized larvae. Interestingly, high copy numbers and differential expression of several miRNA passenger strands (miRNA*) suggest their potential roles in host-parasitoid interaction. In conclusion, expression profiling of miRNAs provided insights into their possible involvement in insect immune response to parasitism and offer an important resource for further studies.     

Deep Sequencing and Microarray Hybridization Identify Conserved and Species-Specific MicroRNAs during Somatic Embryogenesis in Hybrid Yellow Poplar

Background

To date, several studies have indicated a major role for microRNAs (miRNAs) in regulating plant development, but miRNA-mediated regulation of the developing somatic embryo is poorly understood, especially during early stages of somatic embryogenesis in hardwood plants. In this study, Solexa sequencing and miRNA microfluidic chips were used to discover conserved and species-specific miRNAs during somatic embryogenesis of hybrid yellow poplar (Liriodendron tulipifera×L. chinense).

Methodology/Principal Findings

A total of 17,214,153 reads representing 7,421,623 distinct sequences were obtained from a short RNA library generated from small RNAs extracted from all stages of somatic embryos. Through a combination of deep sequencing and bioinformatic analyses, we discovered 83 sequences with perfect matches to known miRNAs from 33 conserved miRNA families and 273 species-specific candidate miRNAs. MicroRNA microarray results demonstrated that many conserved and species-specific miRNAs were expressed in hybrid yellow poplar embryos. In addition, the microarray also detected another 149 potential miRNAs, belonging to 29 conserved families, which were not discovered by deep sequencing analysis. The biological processes and molecular functions of the targets of these miRNAs were predicted by carrying out BLAST search against Arabidopsis thaliana GenBank sequences and then analyzing the results with Gene Ontology.

Conclusions

Solexa sequencing and microarray hybridization were used to discover 232 candidate conserved miRNAs from 61 miRNA families and 273 candidate species-specific miRNAs in hybrid yellow poplar. In these predicted miRNAs, 64 conserved miRNAs and 177 species-specific miRNAs were detected by both sequencing and microarray hybridization. Our results suggest that miRNAs have wide-ranging characteristics and important roles during all stages of somatic embryogenesis in this economically important species.     

Revealing liver specific microRNAs linked with carbohydrate metabolism of farmed carp,Labeo rohita (Hamilton, 1822)

The role of microRNA in gene regulation during developmental biology has been well depicted in several organisms. The present study was performed to investigate miRNAs role in the liver tissues during carbohydrate metabolism and their targets in the farmed carp rohu, Labeo rohita, which is economically important species in aquaculture. Using Illumina-HiSeq technology, a total of 22,612,316; 44,316,046 and 13,338,434 clean reads were obtained from three small-RNA libraries. We have identified 138 conserved and 161 novel miRNAs and studies revealed that miR-22, miR-122, miR-365, miR-200, and miR-146 are involved in carbohydrate metabolism. Further analysis depicted mature miRNA and their predicted target sites in genes that were involved in developmental biology, cellular activities, transportation, etc. This is the first report of the presence of miRNAs in liver tissue of rohu and their comparative profile linked with metabolism serves as a vital resource as a biomarker.     

Normalization of microRNA expression levels in quantitative RT-PCR assays: identification of suitable reference RNA targets in normal and cancerous human solid tissues

Proper normalization is a critical but often an underappreciated aspect of quantitative gene expression analysis. This study describes the identification and characterization of appropriate reference RNA targets for the normalization of microRNA (miRNA) quantitative RT-PCR data. miRNA microarray data from dozens of normal and disease human tissues revealed ubiquitous and stably expressed normalization candidates for evaluation by qRT-PCR. miR-191 and miR-103, among others, were found to be highly consistent in their expression across 13 normal tissues and five pair of distinct tumor/normal adjacent tissues. These miRNAs were statistically superior to the most commonly used reference RNAs used in miRNA qRT-PCR experiments, such as 5S rRNA, U6 snRNA, or total RNA. The most stable normalizers were also highly conserved across flash-frozen and formalin-fixed paraffin-embedded lung cancer tumor/NAT sample sets, resulting in the confirmation of one well-documented oncomir (let-7a), as well as the identification of novel oncomirs. These findings constitute the first report describing the rigorous normalization of miRNA qRT-PCR data and have important implications for proper experimental design and accurate data interpretation.     

Identification and comparative analysis of the pearl oyster <Emphasis Type="Italic">Pinctada fucata</Emphasis> hemocytes microRNAs in response to <Emphasis Type="Italic">Vibrio alginolyticus</Emphasis> infection

MicroRNAs (miRNAs) are a class of noncoding RNA molecules that function as negative regulators of gene expression and play important roles in a wide spectrum of biological processes, including in immune response. However, the physiological regulation function of Pinctada fucata miRNAs, specially their immunomodulation has not been explored yet. Here, two small RNA libraries from hemocytes of P. fucata with or without Vibrio alginolyticus infection were constructed and sequenced using the high-throughput Illumina deep sequencing technology. In total, 11,939,992 and 11,083,327 raw reads, corresponding to 10,993,546 and 9,988,179 clean reads, were respectively obtained in the control and infected libraries. A total of 276 miRNAs, including 225 known miRNAs and 51 putative novel miRNAs, were identified by bioinformatic analysis. By using pairwise comparison between two libraries, 93 miRNAs were found to be significantly differentially expressed, with 42 and 51 miRNAs exhibiting up-regulation and down-regulation, respectively. Thereinto, some known miRNAs were considered to be immune-related. Real-time PCR were implemented for 6 miRNAs co-expressed in the control and infected samples, and agreement was confirmed between the high-throughput sequencing and real-time PCR data. After miRNA targets were predicted, GO and KEGG pathway enrichment analysis were performed, and the results indicated that ten of the differentially expressed miRNAs were involved in immune-related pathways, and might participate in the host immune response to V. alginolyticus. These results of identification and comparative analysis of miRNAs might deepen our understanding of host-pathogen interactions and immune defense mechanisms in P. fucata.