马蔺根系响应Cd胁迫的miRNA高通量测序分析

为了解Cd胁迫下马蔺﹝Iris lactea var. chinensis ( Fisch.) Koidz.〕根系miRNA的表达模式,采用高通量测序法对100μmol·L-1 Cd胁迫0(CK)和24 h(Cd)后马蔺根系的sRNA文库(分别为CK和Cd文库)进行分析,筛选出显著差异表达的miRNA,并对这些miRNA的靶基因功能进行预测;在此基础上,采用qRT-PCR技术对部分miRNA及其靶基因的表达模式进行验证。结果表明：在CK和Cd文库中,未注释的sRNA序列较多,分别占各自sRNA特异序列总数的86.4%和80.5%;在已注释的sRNA序列中,miRNA所占比例最低(分别为0.3%和0.5%),而rRNA所占比例最高(分别为9.4%和11.8%);2个文库中的sRNA长度主要为21~24 nt,且均以21 nt为最多。从Cd胁迫下马蔺根系sRNA中共筛选出32个显著差异表达的miRNA,其中20个miRNA表达量下调(分别属于miR165、miR166、miR167、miR168、miR390和miR396家族),12个miRNA表达量上调。功能预测结果表明：这些miRNA靶基因的功能主要集中在生物学过程、细胞组分和分子功能3个方面;而从KEGG通路富集分析看,富集在核糖体、氨基酸生物合成和碳代谢3个通路上的差异表达miRNA的靶基因数分别为122、88和82个。 qRT-PCR验证结果表明：在CK和Cd文库间,11个差异表达miRNA及8个靶基因的相对表达量均有显著差异(P<0.05);其中,11个miRNA相对表达量的上调和下调趋势与上述筛选结果一致,并且miRNA相对表达量上调时,其靶基因的相对表达量下调,反之亦然,说明Cd胁迫条件下马蔺根系的miRNA负调控其靶基因的表达,并且这些靶基因主要参与编码转录因子、HD-ZIP蛋白和信号蛋白等过程。     

Identification and comparative analysis of the miRNA expression profiles from four tissues of Micropterus salmoides using deep sequencing

In the present study, four small RNA libraries were constructed from an M. salmoides population and sequenced using deep sequencing technology. A total of 9,888,822; 8,519,365; 20,566,198; and 15,762,254 raw reads representing 666,097; 755,711; 978,923; and 840,175 unique sequences were obtained from the spleen, liver, kidney, and muscle libraries, respectively. As a result, 509 known miRNAs belonging to 143 families and 1157 novel miRNAs were identified. The miRNAs displayed diverse expression levels among the four libraries, among which most of the known miRNAs were expressed at higher levels than the novel miRNAs. Furthermore, stem-loop qRT-PCR was applied to validate and profile the expression of the differentially expressed miRNAs in the four different tissues, which revealed that some miRNAs showed tissue specific expression. The identification of miRNAs in M. salmoides will provide new information and enhance our understanding of the functions of miRNAs in regulating biological processes.     

Subtyping glioblastoma by combining miRNA and mRNA expression data using compressed sensing-based approach

In the clinical practice, many diseases such as glioblastoma, leukemia, diabetes, and prostates have multiple subtypes. Classifying subtypes accurately using genomic data will provide individualized treatments to target-specific disease subtypes. However, it is often difficult to obtain satisfactory classification accuracy using only one type of data, because the subtypes of a disease can exhibit similar patterns in one data type. Fortunately, multiple types of genomic data are often available due to the rapid development of genomic techniques. This raises the question on whether the classification performance can significantly be improved by combining multiple types of genomic data. In this article, we classified four subtypes of glioblastoma multiforme (GBM) with multiple types of genome-wide data (e.g., mRNA and miRNA expression) from The Cancer Genome Atlas (TCGA) project. We proposed a multi-class compressed sensing-based detector (MCSD) for this study. The MCSD was trained with data from TCGA and then applied to subtype GBM patients using an independent testing data. We performed the classification on the same patient subjects with three data types, i.e., miRNA expression data, mRNA (or gene expression) data, and their combinations. The classification accuracy is 69.1% with the miRNA expression data, 52.7% with mRNA expression data, and 90.9% with the combination of both mRNA and miRNA expression data. In addition, some biomarkers identified by the integrated approaches have been confirmed with results from the published literatures. These results indicate that the combined analysis can significantly improve the accuracy of classifying GBM subtypes and identify potential biomarkers for disease diagnosis.     

Integrated analysis of mutations,miRNA and mRNA expression in glioblastoma

Background  

Glioblastoma arises from complex interactions between a variety of genetic alterations and environmental perturbations. Little attention has been paid to understanding how genetic variations, altered gene expression and microRNA (miRNA) expression are integrated into networks which act together to alter regulation and finally lead to the emergence of complex phenotypes and glioblastoma.     

Integrated analyses of zebrafish miRNA and mRNA expression profiles identify miR-29b and miR-223 as potential regulators of optic nerve regeneration

Background

Unlike mammals, zebrafish have the ability to regenerate damaged parts of their central nervous system (CNS) and regain functionality of the affected area. A better understanding of the molecular mechanisms involved in zebrafish regeneration may therefore provide insight into how CNS repair might be induced in mammals. Although many studies have described differences in gene expression in zebrafish during CNS regeneration, the regulatory mechanisms underpinning the differential expression of these genes have not been examined.

Results

We used microarrays to analyse and integrate the mRNA and microRNA (miRNA) expression profiles of zebrafish retina after optic nerve crush to identify potential regulatory mechanisms that underpin central nerve regeneration. Bioinformatic analysis identified 3 miRNAs and 657 mRNAs that were differentially expressed after injury. We then combined inverse correlations between our miRNA expression and mRNA expression, and integrated these findings with target predictions from TargetScan Fish to identify putative miRNA-gene target pairs. We focused on two over-expressed miRNAs (miR-29b and miR-223), and functionally validated seven of their predicted gene targets using RT-qPCR and luciferase assays to confirm miRNA-mRNA binding. Gene ontology analysis placed the miRNA-regulated genes (eva1a, layna, nefmb, ina, si:ch211-51a6.2, smoc1, sb:cb252) in key biological processes that included cell survival/apoptosis, ECM-cytoskeleton signaling, and heparan sulfate proteoglycan binding,

Conclusion

Our results suggest a key role for miR-29b and miR-223 in zebrafish regeneration. The identification of miRNA regulation in a zebrafish injury model provides a framework for future studies in which to investigate not only the cellular processes required for CNS regeneration, but also how these mechanisms might be regulated to promote successful repair and return of function in the injured mammalian brain.

Electronic supplementary material

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

Naive and primed murine pluripotent stem cells have distinct miRNA expression profiles

Over the last years, the microRNA (miRNA) pathway has emerged as a key component of the regulatory network of pluripotency. Although clearly distinct states of pluripotency have been described in vivo and ex vivo, differences in miRNA expression profiles associated with the developmental modulation of pluripotency have not been extensively studied so far. Here, we performed deep sequencing to profile miRNA expression in naive (embryonic stem cell [ESC]) and primed (epiblast stem cell [EpiSC]) pluripotent stem cells derived from mouse embryos of identical genetic background. We developed a graphical representation method allowing the rapid identification of miRNAs with an atypical profile including mirtrons, a small nucleolar RNA (snoRNA)-derived miRNA, and miRNAs whose biogenesis may differ between ESC and EpiSC. Comparison of mature miRNA profiles revealed that ESCs and EpiSCs exhibit very different miRNA signatures with one third of miRNAs being differentially expressed between the two cell types. Notably, differential expression of several clusters, including miR290-295, miR17-92, miR302/367, and a large repetitive cluster on chromosome 2, was observed. Our analysis also showed that differentiation priming of EpiSC compared to ESC is evidenced by changes in miRNA expression. These dynamic changes in miRNAs signature are likely to reflect both redundant and specific roles of miRNAs in the fine-tuning of pluripotency during development.     

Integrating miRNA and mRNA expression profiles in response to heat stress-induced injury in rat small intestine

The molecular mechanisms underlying the pathophysiology of heat stress in the small intestine remain undefined. Furthermore, little information is available concerning changes in microRNA (miRNA) expression following heat stress. The present study sought to evaluate miRNA and mRNA expression profiles in the rat small intestine in response to heat stress. Male Sprague-Dawley rats were subjected to 2?h of heat stress daily for ten consecutive days. Rats were sacrificed at specific time points immediately following heat treatment, and morphological changes in the small intestine were determined. The miRNA and mRNA expression profiles from sample of small intestine were evaluated by microarray analysis. Heat stress caused pronounced morphological damage in the rat small intestine, most severe within the jejunum after 3?days of heat treatment. A mRNA microarray analysis found 270 genes to be up-regulated and 122 genes down-regulated (P?≤?0.01, ≥2.0-fold change) in the jejunum after heat treatment. A miRNA microarray analysis found 18 miRNAs to be up-regulated and 11 down-regulated in the jejunum after heat treatment (P?≤?0.05). Subsequent bioinformatic analyses of the differentially expressed mRNAs and miRNAs were carried out to integrate miRNA and mRNA expression and revealed that alterations in mRNA following heat stress were negatively correlated with miRNA expression. These findings significantly advance our understanding of the regulatory mechanisms underlying the pathophysiology of heat stress-induced injury in the small intestine, specifically with regard to miRNAs.     

Computational analysis of drought stress-associated miRNAs and miRNA co-regulation network in Physcomitrella patens

miRNAs are non-coding small RNAs that involve diverse biological processes. Until now, little is known about their roles in plant drought resistance. Physcomitrella patens is highly tolerant to drought; however, it is not clear about the basic biology of the traits that contribute P. patens this important character. In this work, we discovered 16 drought stress-associated miRNA (DsAmR) families in P. patens through computational analysis. Due to the possible discrepancy of expression periods and tissue distributions between potential DsAmRs and their targeting genes, and the existence of false positive results in computational identification, the prediction results should be examined with further experimental validation. We also constructed an miRNA co-regulation network, and identified two network hubs, miR902a-5p and miR414, which may play important roles in regulating drought-resistance traits. We distributed our results through an online database named ppt-miRBase, which can be accessed at http://bioinfor.cnu.edu.cn/ppt_miRBase/index.php. Our methods in finding DsAmR and miRNA co-regulation network showed a new direction for identifying miRNA functions.     

Cd胁迫对黄菖蒲幼苗4种抗氧化酶活性的影响

采用水培法对Cd胁迫下黄菖蒲（Iris pseudacorus L．）幼苗叶片和根系中超氧化物歧化酶（SOD）、过氧化物酶（POD）、过氧化氢酶（CAT）及抗坏血酸过氧化物酶（APX）的活性进行了研究。结果表明,10、40和120mg·L^-1Cd胁迫下,黄菖蒲幼苗叶片和根系中4种酶活性的变化不同。10和40mg·L^-1 cd胁迫下,黄菖蒲幼苗叶片和根系中的POD及APX活性、叶片中的SOD活性及根系中的CAT均明显高于对照;在120mg·L^-1 Cd胁迫下,叶片中的POD活性及根系中的POD和CAT活性均高于对照;各处理组根系中的SOD活性均低于对照。随处理时间的延长,40和120mg·L^-1Cd胁迫处理组叶片的CAT活性和120mg·L^-1Cd胁迫处理组根系的APX活性逐渐降低,其他处理组不同酶的活性逐渐升高或先升后降。黄菖蒲叶片及根系中的4种酶对Cd胁迫的响应能力有差异,其中POD可能是黄菖蒲耐Cd胁迫的主要抗性诱导酶。     

Integrated analysis of miRNA and mRNA expression in childhood medulloblastoma compared with neural stem cells

Medulloblastoma (MB) is the most common malignant brain tumor in children and a leading cause of cancer-related mortality and morbidity. Several molecular sub-types of MB have been identified, suggesting they may arise from distinct cells of origin. Data from animal models indicate that some MB sub-types arise from multipotent cerebellar neural stem cells (NSCs). Hence, microRNA (miRNA) expression profiles of primary MB samples were compared to CD133+ NSCs, aiming to identify deregulated miRNAs involved in MB pathogenesis. Expression profiling of 662 miRNAs in primary MB specimens, MB cell lines, and human CD133+ NSCs and CD133- neural progenitor cells was performed by qRT-PCR. Clustering analysis identified two distinct sub-types of MB primary specimens, reminiscent of sub-types obtained from their mRNA profiles. 21 significantly up-regulated and 12 significantly down-regulated miRNAs were identified in MB primary specimens relative to CD133+ NSCs (p<0.01). The majority of up-regulated miRNAs mapped to chromosomal regions 14q32 and 17q. Integration of the predicted targets of deregulated miRNAs with mRNA expression data from the same specimens revealed enrichment of pathways regulating neuronal migration, nervous system development and cell proliferation. Transient over-expression of a down-regulated miRNA, miR-935, resulted in significant down-regulation of three of the seven predicted miR-935 target genes at the mRNA level in a MB cell line, confirming the validity of this approach. This study represents the first integrated analysis of MB miRNA and mRNA expression profiles and is the first to compare MB miRNA expression profiles to those of CD133+ NSCs. We identified several differentially expressed miRNAs that potentially target networks of genes and signaling pathways that may be involved in the transformation of normal NSCs to brain tumor stem cells. Based on this integrative approach, our data provide an important platform for future investigations aimed at characterizing the role of specific miRNAs in MB pathogenesis.     

Differential expression profiles of miRNA in the serum of sarcopenic rats

As the geriatric population and life expectancy increase, the interest in preventing geriatric diseases, such as sarcopenia, is increasing. However, the causes of sarcopenia are unclear, and current diagnostic methods for sarcopenia are unreliable. We hypothesized that the changes in the expression of certain miRNAs may be associated with the pathophysiology of sarcopenia. Herein, we analyzed the miRNA expression profiles in the blood of young (3-months-old) healthy rats, old sarcopenic (17-months-old) rats, and age-matched (17-months-old) control rats. The changes in miRNA expression levels were analyzed using Bowtie 2 software. A total of 523 miRNAs were detected in the rat serum. Using scatter plots and clustering heatmap data, we found 130 miRNAs that were differentially expressed in sarcopenic rats (>2-fold change) compared to the expression in young healthy and age-matched control rats. With a threshold of >5-fold change, we identified 14 upregulated miRNAs, including rno-miR-133b-3p, rno-miR-133a-3p, rno-miR-133c, rno-miR-208a-3p, and rno-miR434-5p among others in the serum of sarcopenic rats. A protein network map based on these 14 miRNAs identified the genes involved in skeletal muscle differentiation, among which Notch1, Egr2, and Myocd represented major nodes. The data obtained in this study are potentially useful for the early diagnosis of sarcopenia and for the identification of novel therapeutic targets for the treatment and/or prevention of sarcopenia.     

Massive analysis of cDNA Ends (MACE) and miRNA expression profiling identifies proatherogenic pathways in chronic kidney disease

Epigenetic dysregulation contributes to the high cardiovascular disease burden in chronic kidney disease (CKD) patients. Although microRNAs (miRNAs) are central epigenetic regulators, which substantially affect the development and progression of cardiovascular disease (CVD), no data on miRNA dysregulation in CKD-associated CVD are available until now. We now performed high-throughput miRNA sequencing of peripheral blood mononuclear cells from ten clinically stable hemodialysis (HD) patients and ten healthy controls, which allowed us to identify 182 differentially expressed miRNAs (e.g., miR-21, miR-26b, miR-146b, miR-155). To test biological relevance, we aimed to connect miRNA dysregulation to differential gene expression. Genome-wide gene expression profiling by MACE (Massive Analysis of cDNA Ends) identified 80 genes to be differentially expressed between HD patients and controls, which could be linked to cardiovascular disease (e.g., KLF6, DUSP6, KLF4), to infection / immune disease (e.g., ZFP36, SOCS3, JUND), and to distinct proatherogenic pathways such as the Toll-like receptor signaling pathway (e.g., IL1B, MYD88, TICAM2), the MAPK signaling pathway (e.g., DUSP1, FOS, HSPA1A), and the chemokine signaling pathway (e.g., RHOA, PAK1, CXCL5). Formal interaction network analysis proved biological relevance of miRNA dysregulation, as 68 differentially expressed miRNAs could be connected to 47 reciprocally expressed target genes. Our study is the first comprehensive miRNA analysis in CKD that links dysregulated miRNA expression with differential expression of genes connected to inflammation and CVD. After recent animal data suggested that targeting miRNAs is beneficial in experimental CVD, our data may now spur further research in the field of CKD-associated human CVD.     

mRNA Targeting to Endoplasmic Reticulum Precedes Ago Protein Interaction and MicroRNA (miRNA)-mediated Translation Repression in Mammalian Cells

MicroRNA (miRNA) binds to the 3′-UTR of its target mRNAs to repress protein synthesis. Extensive research was done to understand the mechanism of miRNA-mediated repression in animal cells. Considering the progress in understanding the mechanism, information about the subcellular sites of miRNA-mediated repression is surprisingly limited. In this study, using an inducible expression system for an miRNA target message, we have delineated how a target mRNA passes through polysome association and Ago2 interaction steps on rough endoplasmic reticulum (ER) before the miRNA-mediated repression sets in. From this study, de novo formed target mRNA localization to the ER-bound polysomes manifested as the earliest event, which is followed by Ago2 micro-ribonucleoprotein binding, and translation repression of target message. Compartmentalization of this process to rough ER membrane ensures enrichment of miRNA-targeted messages and micro-ribonucleoprotein components on ER upon reaching a steady state.