Designing filters for fast-known NcRNA identification

Detecting members of known noncoding RNA (ncRNA) families in genomic DNA is an important part of sequence annotation. However, the most widely used tool for modeling ncRNA families, the covariance model (CM), incurs a high-computational cost when used for genome-wide search. This cost can be reduced by using a filter to exclude sequences that are unlikely to contain the ncRNA of interest, applying the CM only where it is likely to match strongly. Despite recent advances, designing an efficient filter that can detect ncRNA instances lacking strong conservation while excluding most irrelevant sequences remains challenging. In this work, we design three types of filters based on multiple secondary structure profiles (SSPs). An SSP augments a regular profile (i.e., a position weight matrix) with secondary structure information but can still be efficiently scanned against long sequences. Multi-SSPbased filters combine evidence from multiple SSP matches and can achieve high sensitivity and specificity. Our SSP-based filters are extensively tested in BRAliBase III data set, Rfam 9.0, and a published soil metagenomic data set. In addition, we compare the SSPbased filters with several other ncRNA search tools including Infernal (with profile HMMs as filters), ERPIN, and tRNAscan-SE. Our experiments demonstrate that carefully designed SSP filters can achieve significant speedup over unfiltered CM search while maintaining high sensitivity for various ncRNA families. The designed filters and filter-scanning programs are available at our website: www.cse.msu.edu/~yannisun/ssp/.     

小立碗藓（Physcomitrella patens）植物microRNA结合靶位预测

利用857条植物miRNA序列对27546条小立碗藓mRNA序列进行搜索,预测出162个植物miRNA家族在小立碗藓中存在结合靶位。miRNA结合靶位数目和miRNA协同作用网络分析结果同时显示,miR482和miR1168在小立碗藓中结合靶位多、协同作用广,提示它们对于小立碗藓可能具有重要生物学功能。52个菜茵衣藻特有的miRNA被预测在小立碗藓中存在结合靶位,显示小立碗藓在从藻类向种子植物进化过程中处在独特演化位置。     

Plant noncoding RNA gene discovery by "single-genome comparative genomics"

Plant genomes have undergone multiple rounds of duplications that contributed massively to the growth of gene families. The structure of resulting families has been studied in depth for protein-coding genes. However, little is known about the impact of duplications on noncoding RNA (ncRNA) genes. Here we perform a systematic analysis of duplicated regions in the rice genome in search of such ncRNA repeats. We observe that, just like their protein counterparts, most ncRNA genes have undergone multiple duplications that left visible sequence conservation footprints. The extent of ncRNA gene duplication in plants is such that these sequence footprints can be exploited for the discovery of novel ncRNA gene families on a large scale. We developed an SVM model that is able to retrieve likely ncRNA candidates among the 100,000+ repeat families in the rice genome, with a reasonably low false-positive discovery rate. Among the nearly 4000 ncRNA families predicted by this means, only 90 correspond to putative snoRNA or miRNA families. About half of the remaining families are classified as structured RNAs. New candidate ncRNAs are particularly enriched in UTR and intronic regions. Interestingly, 89% of the putative ncRNA families do not produce a detectable signal when their sequences are compared to another grass genome such as maize. Our results show that a large fraction of rice ncRNA genes are present in multiple copies and are species-specific or of recent origin. Intragenome comparison is a unique and potent source for the computational annotation of this major class of ncRNA.     

Identification of novel microRNAs in rice (Oryza sativa) based on the cleavage signals in precursors

MicroRNAs (miRNAs) are crucial regulators of gene expression in plants and a growing number of novel miRNA genes have been cloned in rice in recent years. However, there is no evidence that all miRNAs have been discovered, especially for those low expression ones which are difficult to be found by conventional methods. By taking advantage of the finding that DCL1-mediated cleavage signals for the processing of the miRNA precursors could be used as the clues for novel miRNAs’ discovery, a genome-wide search for rice miRNA candidates was carried out. As a result, 51 previously validated miRNAs and 24 novel miRNA candidates were discovered. After target prediction and degradome sequencing data-based validation, coupled with reverse approach retest, 10 miRNA candidate–mRNA target pairs were further identified, providing a basis for in-depth functional analysis of these miRNA candidates. Besides, some isomiRs found in this study showed more likely to be the real miRNAs. We also found an exceptional example which did not obey the rule that 22-nt miRNAs have the ability to trigger the phased siRNAs production from the cleaved targets.     

Global analysis of non-coding small RNAs in Arabidopsis in response to jasmonate treatment by deep sequencing technology

In plants, non-coding small RNAs play a vital role in plant development and stress responses. To explore the possible role of non-coding small RNAs in the regulation of the jasmonate (JA) pathway, we compared the non-coding small RNAs between the JA-deficient aos mutant and the JA-treated wild type Arabidopsis via high-throughput sequencing. Thirty new miRNAs and 27 new miRNA candidates were identified through bioinformatics approach. Forty-nine known miRNAs (belonging to 24 families), 15 new miRNAs and new miRNA candidates (belonging to 11 families) and 3 tasiRNA families were induced by JA, whereas 1 new miRNA, 1 tasiRNA family and 22 known miRNAs (belonging to 9 families) were repressed by JA.     

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 of putative miRNAs from the deep-branching unicellular flagellates

MicroRNAs (miRNAs) are a class of extensively studied RNAi-associated small RNAs that play a critical role in eukaryotic gene regulation. However, knowledge on the miRNA and its regulation in unicellular eukaryotes is very limited. In order to obtain a better understanding on the origin of miRNA regulation system, we used deep-sequencing technology to investigate the miRNA expression pattern in four deep-branching unicellular flagellates: Giardia lamblia, Trichomonas vaginalis, Tritrichomonas foetus, and Pentatrichomonas hominis. In addition to the known miRNAs that have been described in G. lamblia and T. vaginalis, we identified 14 ancient animal miRNA families and 13 plant-specific families. Bioinformatics analysis also identified four novel miRNA candidates with reliable precursor structures derived from mature tRNAs. Our results indicated that miRNAs are likely to be a general feature for gene regulation throughout unicellular and multicellular eukaryotes and some of them may derive from unconventional ncRNAs such as snoRNA and tRNA.     

Sorting the wheat from the chaff: identifying miRNAs in genomic survey sequences of Triticum aestivum chromosome 1AL

Individual chromosome-based studies of bread wheat are beginning to provide valuable structural and functional information about one of the world's most important crops. As new genome sequences become available, identifying miRNA coding sequences is arguably as important a task as annotating protein coding sequences, but one that is not as well developed. We compared conservation-based identification of conserved miRNAs in 1.5× coverage survey sequences of wheat chromosome 1AL with a predictive method based on pre-miRNA hairpin structure alone. In total, 42 sequences expected to encode conserved miRNAs were identified on chromosome 1AL, including members of several miRNA families that have not previously been reported to be expressed in T. aestivum. In addition, we demonstrate that a number of sequences previously annotated as novel wheat miRNAs are closely related to transposable elements, particularly Miniature Inverted Terminal repeat Elements (MITEs). Some of these TE-miRNAs may well have a functional role, but separating true miRNA coding sequences from TEs in genomic sequences is far from straightforward. We propose a strategy for annotation to minimize the risk of mis-identifying TE sequences as miRNAs.     

Cloning, structural analysis, and mapping of the B30 and B7 multigenic families to the major histocompatibility complex (MHC) and other chromosomal regions

 We present the cloning, structural analysis, and mapping of new members belonging to two multigenic families, the B30-RING finger family and the B7.1-B7.2 family, as well as two genes derived by exon shuffling from members of these families. Eight new members were found and three of them map to the human major histocompatibilitiy complex (MHC) region. Phylogenic and physical mapping analysis allowed us to decipher the evolutionary story of these two multigenic families and to shed light on the evolution of the MHC region. We also show that a deductive analysis can be used to predict the existence of a given gene. Received: 17 February 1997 / Revised: 17 March 1997