应用生物信息学寻找山羊新的microRNA分子及其实验验证

摘要: microRNA(miRNA)是一类长约22个碱基的非编码RNA分子, 在转录后水平调节基因的表达及其在细胞的增殖、分化、凋亡等过程中起着重要的调控作用。根据miRNA分子具有一定的保守性, 文章将人、小鼠、牛、猪和狗5种哺乳动物已知的miRNA分子与NCBI公布的与山羊具有极高同源性的绵羊基因组序列对比, 获得11条miRNA候选分子, 然后通过逆转录聚合酶链反应(RT-PCR)验证, 发现在山羊脑组织中这11条分子均有表达, 肝脏组织中有5条分子表达, 初步确定为山羊新的miRNA分子, 为寻找山羊miRNA提供了新的 思路。     

webFOG: A web tool to map genomic features onto genes

A large number of new genomic features are being discovered using high throughput techniques. The next challenge is to automatically map them to the reference genome for further analysis and functional annotation. We have developed a tool that can be used to map important genomic features to the latest version of the human genome and also to annotate new features. These genomic features could be of many different source types, including miRNAs, microarray primers or probes, Chip-on-Chip data, CpG islands and SNPs to name a few. A standalone version and web interface for the tool can be accessed through: http://populationhealth.qimr.edu.au/cgi-bin/webFOG/index.cgi. The project details and source code is also available at http://www.bioinformatics.org/webfog.     

Distribution of miRNA genes in the pig genome

Background

Recent completion of swine genome may simplify the production of swine as a large biomedical model. Here we studied sequence and location of known swine miRNA genes, key regulators of protein-coding genes at the level of RNA, and compared them to human and mouse data to prioritize future molecular studies.

Results

Distribution of miRNA genes in pig genome shows no particular relation to different genomic features including protein coding genes - proportions of miRNA genes in intergenic regions, introns and exons roughly agree with the size of these regions in the pig genome. Our analyses indicate that host genes harbouring intragenic miRNAs are longer from other protein-coding genes, however, no important GO enrichment was found. Swine mature miRNAs show high sequence similarity to their human and mouse orthologues. Location of miRNA genes relative to protein-coding genes is also similar among studied species, however, there are differences in the precise position in particular intergenic regions and within particular hosts. The most prominent difference between pig and human miRNAs is a large group of pig-specific sequences (53% of swine miRNAs). We found no evidence that this group of evolutionary new pig miRNAs is different from old miRNAs genes with respect to genomic location except that they are less likely to be clustered.

Conclusions

There are differences in precise location of orthologues miRNA genes in particular intergenic regions and within particular hosts, and their meaning for coexpression with protein-coding genes deserves experimental studies. Functional studies of a large group of pig-specific sequences in future may reveal limits of the pig as a model organism to study human gene expression.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-015-0166-3) contains supplementary material, which is available to authorized users.     

水稻MicroRNA的预测及实验验证

根据已报道水稻pre-miRNA的序列与结构信息,利用支持向量机（support vector machine, SVM）方法在miRNA前体上预测成熟区,产生一个模型——mature-SVM.它预测水稻成熟区的敏感性和特异性分别为86.7% 和100%;然后,用这个模型对从水稻基因组中筛选出的46.501条pre-miRNA进行成熟链预测,此外再根据miRNA的作用原理用blast程序所进一步的筛选,得到了127条pre-miRNA及成熟miRNA;除去其中已知的21条,最后得到106条候选的新的水稻miRNA. 从中随机挑取10条进行Northern验证,结果有4条miRNA得到确认.     

CID-miRNA: a web server for prediction of novel miRNA precursors in human genome

microRNAs (miRNA) are a class of non-protein coding functional RNAs that are thought to regulate expression of target genes by direct interaction with mRNAs. miRNAs have been identified through both experimental and computational methods in a variety of eukaryotic organisms. Though these approaches have been partially successful, there is a need to develop more tools for detection of these RNAs as they are also thought to be present in abundance in many genomes. In this report we describe a tool and a web server, named CID-miRNA, for identification of miRNA precursors in a given DNA sequence, utilising secondary structure-based filtering systems and an algorithm based on stochastic context free grammar trained on human miRNAs. CID-miRNA analyses a given sequence using a web interface, for presence of putative miRNA precursors and the generated output lists all the potential regions that can form miRNA-like structures. It can also scan large genomic sequences for the presence of potential miRNA precursors in its stand-alone form. The web server can be accessed at http://mirna.jnu.ac.in/cidmirna/.     

人源microRNA分子表达库的构建

microRNA(miRNA)是一类长度为22nt左右的单链非编码小RNA分子,通过与靶mRNA分子结合而沉默其表达.目前,虽然在多种生物中发现了大量的miRNA,但对它们的功能还知之甚少.为了深入研究miRNA的功能,构建了一个包括170多种人源miRNA表达载体的miRNA分子表达库,并对部分表达载体采用RNA印迹及双荧光素酶分析技术进行验证.实验证明:这些miRNA表达载体在HEK-293细胞内可以高水平表达miRNA前体和成熟的miRNA,并且能抑制含有相应靶位点的报告基因的表达.这些结果表明:该miRNA表达库可以表达功能miRNA,并可用于miRNA功能的筛选和研究.     

Improved annotation of C. elegans microRNAs by deep sequencing reveals structures associated with processing by Drosha and Dicer

MicroRNAs (miRNAs) are small regulatory RNAs that are essential in all studied metazoans. Research has focused on the prediction and identification of novel miRNAs, while little has been done to validate, annotate, and characterize identified miRNAs. Using Illumina sequencing, ～20 million small RNA sequences were obtained from Caenorhabditis elegans. Of the 175 miRNAs listed on the miRBase database, 106 were validated as deriving from a stem-loop precursor with hallmark characteristics of miRNAs. This result suggests that not all sequences identified as miRNAs belong in this category of small RNAs. Our large data set of validated miRNAs facilitated the determination of general sequence and structural characteristics of miRNAs and miRNA precursors. In contrast to previous observations, we did not observe a preference for the 5' nucleotide of the miRNA to be unpaired compared to the 5' nucleotide of the miRNA*, nor a preference for the miRNA to be on either the 5' or 3' arm of the miRNA precursor stem-loop. We observed that steady-state pools of miRNAs have fairly homogeneous termini, especially at their 5' end. Nearly all mature miRNA-miRNA* duplexes had two nucleotide 3' overhangs, and there was a preference for a uracil in the first and ninth position of the mature miRNA. Finally, we observed that specific nucleotides and structural distortions were overrepresented at certain positions adjacent to Drosha and Dicer cleavage sites. Our study offers a comprehensive data set of C. elegans miRNAs and their precursors that significantly decreases the uncertainty associated with the identity of these molecules in existing databases.     

SynBrowse: a synteny browser for comparative sequence analysis

MOTIVATION: The recent efforts of various sequence projects to sequence deeply into various phylogenies provide great resources for comparative sequence analysis. A generic and portable tool is essential for scientists to visualize and analyze sequence comparisons. RESULTS: We have developed SynBrowse, a synteny browser for visualizing and analyzing genome alignments both within and between species. It is intended to help scientists study macrosynteny, microsynteny and homologous genes between sequences. It can also aid with the identification of uncharacterized genes, putative regulatory elements and novel structural features of a species. SynBrowse is a GBrowse (the Generic Genome Browser) family software tool that runs on top of the open source BioPerl modules. It consists of two components: a web-based front end and a set of relational database back ends. Each database stores pre-computed alignments from a focus sequence to reference sequences in addition to the genome annotations of the focus sequence. The user interface lets end users select a key comparative alignment type and search for syntenic blocks between two sequences and zoom in to view the relationships among the corresponding genome annotations in detail. SynBrowse is portable with simple installation, flexible configuration, convenient data input and easy integration with other components of a model organism system. AVAILABILITY: The software is available at http://www.gmod.org CONTACT: vbrendel@iastate.edu     

Alignment of BLAST high-scoring segment pairs based on the longest increasing subsequence algorithm

MOTIVATION:The popular BLAST algorithm is based on a local similarity search strategy, so its high-scoring segment pairs (HSPs) do not have global alignment information. When scientists use BLAST to search for a target protein or DNA sequence in a huge database like the human genome map, the existence of repeated fragments, homologues or pseudogenes in the genome often makes the BLAST result filled with redundant HSPs. Therefore, we need a computational strategy to alleviate this problem. RESULTS: In the gene discovery group of Celera Genomics, I developed a two-step method, i.e. a BLAST step plus an LIS step, to align thousands of cDNA and protein sequences into the human genome map. The LIS step is based on a mature computational algorithm, Longest Increasing Subsequence (LIS) algorithm. The idea is to use the LIS algorithm to find the longest series of consecutive HSPs in the BLAST output. Such a BLAST+LIS strategy can be used as an independent alignment tool or as a complementary tool for other alignment programs like Sim4 and GenWise. It can also work as a general purpose BLAST result processor in all sorts of BLAST searches. Two examples from Celera were shown in this paper.     

miR-503 represses human cell proliferation and directly targets the oncogene DDHD2 by non-canonical target pairing

Background

The pathways regulating the transition of mammalian cells from quiescence to proliferation are mediated by multiple miRNAs. Despite significant improvements in our understanding of miRNA targeting, the majority of miRNA regulatory networks are still largely unknown and require experimental validation.

Results

Here we identified miR-503, miR-103, and miR-494 as negative regulators of proliferation in primary human cells. We experimentally determined their genome wide target profiles using RNA-induced silencing complex (RISC) immunoprecipitations and gene expression profiling. Analysis of the genome wide target profiles revealed evidence of extensive regulation of gene expression through non-canonical target pairing by miR-503. We identified the proto-oncogene DDHD2 as a target of miR-503 that requires pairing outside of the canonical 5′ seed region of miR-503, representing a novel mode of miRNA-target pairing. Further bioinformatics analysis implicated miR-503 and DDHD2 in breast cancer tumorigenesis.

Conclusions

Our results provide an extensive genome wide set of targets for miR-503, miR-103, and miR-494, and suggest that miR-503 may act as a tumor suppressor in breast cancer by its direct non-canonical targeting of DDHD2.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1279-9) contains supplementary material, which is available to authorized users.     

COBALT: constraint-based alignment tool for multiple protein sequences

MOTIVATION: A tool that simultaneously aligns multiple protein sequences, automatically utilizes information about protein domains, and has a good compromise between speed and accuracy will have practical advantages over current tools. RESULTS: We describe COBALT, a constraint based alignment tool that implements a general framework for multiple alignment of protein sequences. COBALT finds a collection of pairwise constraints derived from database searches, sequence similarity and user input, combines these pairwise constraints, and then incorporates them into a progressive multiple alignment. We show that using constraints derived from the conserved domain database (CDD) and PROSITE protein-motif database improves COBALT's alignment quality. We also show that COBALT has reasonable runtime performance and alignment accuracy comparable to or exceeding that of other tools for a broad range of problems. AVAILABILITY: COBALT is included in the NCBI C++ toolkit. A Linux executable for COBALT, and CDD and PROSITE data used is available at: ftp://ftp.ncbi.nlm.nih.gov/pub/agarwala/cobalt     

病毒编码的miRNA:基因表达新的调控因子

微小RNA(microRNA,miRNA)是一类长约22个核苷酸的RNA,在数量、序列、结构、表达和功能上具有多样性。目前,通过生物信息学手段和分子克隆方法,已发现了3518种miRNA,在控制细胞的生长发育、分化、凋亡等过程中发挥着十分重要的作用。最近研究发现疱疹病毒、多瘤病毒、逆转录病毒的某些病毒基因组也能够编码miRNA,这些miRNA在调控病毒基因自身表达以及病毒与宿主相互作用方面可能起重要的作用。某些病毒甚至能够利用宿主体内的miRNA调控其自身表达。找出病毒可能编码的miRNA,探索其对病毒感染、复制、表达的作用,有助于病毒分子生物学的研究,也会为研发防治病毒的新方法和新途径提供新的思路。     

Discovery of novel microRNAs in rat kidney using next generation sequencing and microarray validation

MicroRNAs (miRNAs) are small non-coding RNAs that regulate a variety of biological processes. The latest version of the miRBase database (Release 18) includes 1,157 mouse and 680 rat mature miRNAs. Only one new rat mature miRNA was added to the rat miRNA database from version 16 to version 18 of miRBase, suggesting that many rat miRNAs remain to be discovered. Given the importance of rat as a model organism, discovery of the completed set of rat miRNAs is necessary for understanding rat miRNA regulation. In this study, next generation sequencing (NGS), microarray analysis and bioinformatics technologies were applied to discover novel miRNAs in rat kidneys. MiRanalyzer was utilized to analyze the sequences of the small RNAs generated from NGS analysis of rat kidney samples. Hundreds of novel miRNA candidates were examined according to the mappings of their reads to the rat genome, presence of sequences that can form a miRNA hairpin structure around the mapped locations, Dicer cleavage patterns, and the levels of their expression determined by both NGS and microarray analyses. Nine novel rat hairpin precursor miRNAs (pre-miRNA) were discovered with high confidence. Five of the novel pre-miRNAs are also reported in other species while four of them are rat specific. In summary, 9 novel pre-miRNAs (14 novel mature miRNAs) were identified via combination of NGS, microarray and bioinformatics high-throughput technologies.     

Computational identification of miRNAs in medicinal plant Senecio vulgaris (Groundsel)

RNAs Interference plays a very important role in gene silencing. In vitro identification of miRNAs is a slow process as it is difficult to isolate them. Nucleotide sequences of miRNAs are highly conserved among the plants and, this form the key feature behind the identification of miRNAs in plant species by homology alignment. In silico identification of miRNAs from EST database is emerging as a novel, faster and reliable approach. Here EST sequences of Senecio vulgaris (Groundsel) were searched against known miRNA sequences by using BLASTN tool. A total of 10 miRNAs were identified from 1956 EST sequences and 115 GSS sequences. The most stable miRNA identified is svu-mir-1. This approach will accelerate advance research in regulation of gene expression in Groundsel by interfering RNAs.