Sex specific expression and distribution of small RNAs in papaya

Background

Regulatory function of small non-coding RNAs (sRNA) in response to environmental and developmental cues has been established. Additionally, sRNA, also plays an important role in maintaining the heterochromatin and centromere structures of the chromosome. Papaya, a trioecious species with recently evolved sex chromosomes, has emerged as an excellent model system to study sex determination and sex chromosome evolution in plants. However, role of small RNA in papaya sex determination is yet to be explored.

Results

We analyzed the high throughput sRNAs reads in the Illumina libraries prepared from male, female, and hermaphrodite flowers of papaya. Using the sRNA reads, we identified 29 miRNAs that were not previously reported from papaya. Including this and two previous studies, a total of 90 miRNAs has been identified in papaya. We analyzed the expression of these miRNAs in each sex types. A total of 65 miRNAs, including 31 conserved and 34 novel mirNA, were detected in at least one library. Fourteen of the 65 miRNAs were differentially expressed among different sex types. Most of the miRNA expressed higher in male flowers were related to the auxin signaling pathways, whereas the miRNAs expressed higher in female flowers were the potential regulators of the apical meristem identity genes. Aligning the sRNA reads identified the sRNA hotspots adjacent to the gaps of the X and Y chromosomes. The X and Y chromosomes sRNA hotspots has a 7.8 and 4.4 folds higher expression of sRNA, respectively, relative to the chromosome wide average. Approximately 75% of the reads aligned to the X chromosome hotspot was identical to that of the Y chromosome hotspot.

Conclusion

By analyzing the large-scale sRNA sequences from three sex types, we identified the sRNA hotspots flanking the gaps of papaya X, Y, and Y^h chromosome. The sRNAs expression patterns in these regions were reminiscent of the pericentromeric region indicating that the only remaining gap in each of these chromosomes is likely the centromere. We also identified 14 differentially expressed miRNAs in male, female and hermaphrodite flowers of papaya. Our results provide valuable information toward understanding the papaya sex determination.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-20) contains supplementary material, which is available to authorized users.     

Identification and bioinformatics analysis of microRNAs associated with stress and immune response in serum of heat-stressed and normal Holstein cows

MicroRNAs (miRNAs) are small single-stranded non-coding RNAs that have an important regulatory function in animal growth and developmental processes. However, the differential expression of miRNA and the role of these miRNAs in heat-stressed Holstein cows are still unknown. In this study, the profile of differentially expressed miRNAs and the target genes analysis in the serum of heat-stressed and normal Holstein cows were investigated by a Solexa deep-sequencing approach and bioinformatics. The data identified 52 differentially expressed miRNAs in 486 known miRNAs which were changed significantly between heat-stressed and normal Holstein cows (fold change >2, P < 0.001). Target genes analysis showed that at least 7 miRNAs (miR-19a, miR-19b, miR-146a, miR-30a-5p, miR-345-3p, miR-199a-3p, and miR-1246) were involved in the response to stress, oxidative stress, development of the immune system, and immune response among the identified 52 differentially expressed miRNAs. Five miRNAs (miR-27b, miR-181a, miR-181b, miR-26a, and miR-146b) were involved in stress and immune responses and the expression of five miRNAs was striking (P < 0.001). In addition, RT-qPCR and deep-sequencing methods showed that 8 miRNAs among the 12 selected miRNAs (miR-19a, miR-19b, miR-27b, miR-30a-5p, miR-181a, miR-181b, miR-345-3p, and miR-1246) were highly expressed in the serum of heat-stressed Holstein cows. GO and KEGG pathway analysis showed that these differentially expressed miRNAs were involved in a pathway that may differentially regulate the expression of stress response and immune response genes. Our study provides an overview of miRNAs expression profile and the interaction between miRNAs and their target genes, which will lead to further understanding of the important roles of miRNAs in heat-stressed Holstein cows.     

固始鸡1日龄和36周龄下丘脑高丰度差异性MiRNA的鉴定及功能预测分析

为鉴定鸡下丘脑发育相关特异性表达miRNA,基于固始鸡1日龄和36周龄下丘脑小RNA的Solexa测序数据,共鉴定到266种2个发育阶段共表达的miRNA,其中157种miRNA的表达水平被显著下调,22种被显著上调.聚类分析显示,鸡下丘脑高丰度差异性miRNA主要集中于let-7、mir-181、mir-30、mir-99、mir-1和mir-17等基因家族.另外,预测了10种高丰度差异性miRNA的靶基因,并进行了相应的GO分析和KEGG通路分析.结果显示,预测靶基因在发育过程、代谢过程、细胞过程和生物学过程调节等4个生物学过程以及细胞周期、粘着斑、TGF-beta信号通路和MAPK信号通路等通路中显著富集.研究结果为进一步揭示miRNA调控鸡下丘脑发育的分子机制提供了有益线索.     

Genome-wide identification and characterization of perirenal adipose tissue microRNAs in rabbits fed a high-fat diet

MicroRNAs (miRNAs) are a class of endogenous single-stranded RNA molecules that play an important role in gene regulation in animals by pairing with target gene mRNA. Extensive evidence shows that miRNAs are key players in metabolic regulation and the development of obesity. However, the systemic understanding of miRNAs in the adipogenesis of obese rabbits need further investigation. Here, seven small RNA libraries from rabbits fed either a standard normal diet (SND; n=3) or high-fat diet (HFD; n=4) were constructed and sequenced. Differentially expressed (DE) miRNAs were identified using the edgeR data analysis package from R. Software miRanda and RNAhybrid were used to predict the target genes of miRNAs. To further explore the functions of DE miRNAs, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis were performed. A total of 81449996 clean reads were obtained from the seven libraries, of which, 52 known DE miRNAs (24 up-regulated, 28 down-regulated) and 31 novel DE miRNAs (14 up-regulated, 17 down-regulated) were identified. GO enrichment analysis revealed that the DE miRNAs target genes were involved in intermediate filament cytoskeleton organization, intermediate filament-based process, and α-tubulin binding. DE miRNAs were involved in p53 signaling, linoleic acid metabolism, and other adipogenesis-related KEGG pathways. Our study further elucidates the possible functions of DE miRNAs in rabbit adipogenesis, contributing to the understanding of rabbit obesity.     

A comprehensive study of construction and analysis of competitive endogenous RNA networks in lung adenocarcinoma

BackgroundLong noncoding RNAs (lncRNAs) have gain increasing attention in lung adenocarcinoma. In this study, we aimed at constructing and analyzing the lncRNAs and the related proteins based competitive endogenous RNA (ceRNA) network.MethodsRNA expression data of lung adenocarcinoma were extracted from the TCGA database. Differentially expressed (DE) lncRNAs, messenger RNAs (mRNAs) and microRNAs (miRNAs) were identified and then a DElncRNA-DEmiRNA-DEmRNA ceRNA network was constructed for lung adenocarcinoma. We also analyzed the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment of the DEgenes. Kaplan-Meier survival curves were also been further utilized for exploring the prognostic factors.ResultsAfter compared and calculated lncRNA, mRNA and miRNA expression profiles between lung adenocarcinoma and normal samples, 1709 differential expressed lncRNAs, 2554 differential expressed mRNAs and 116 differential expressed miRNAs were finally identified. Afterwards, a lncRNA mediated ceRNA network was constructed, according to the interactions among 544 pairs of DElncRNA-DEmiRNA relationships and 47 pairs of DEmiRNA-DEmRNA relationships. As for the survival analyses, we found 10 DElncRNAs, 25 DEmRNAs and 7 miRNAs have statistically prognostic significance for overall survival, respectively.ConclusionsThis study provides meaningful information for deeper understanding the underlying molecular mechanism of lung adenocarcinoma and for evaluating prognosis, which could monitor recurrence, guide clinical treatment drugs and subsequent related researches.     

黄瓜AGAMOUS同源基因的分离和表达(英文)

黄瓜花发育的早期均有雌蕊和雄蕊原基的分化,但在发育过程中,由于雌蕊或雄蕊的发育受到阻滞,导致雄花和雌花的形成。近年来在拟南芥和金鱼草等植物中遗传学的研究表明,花器官的特征是由同源异形基因决定的。在拟南芥中,由于ＡＧ在决定雄蕊和雌蕊特征方面起重要作用,本研究利用 ＲＴ－ＰＣＲ技术,从黄瓜的雌、雄蕊中分离出 ＡＧ的同源基因,并对其在花发育过程中的表达和可能作用进行了分析。首先,根据ＡＧ同源基因的保守区域设计５’简并引物５’－ＧＡ（Ａ／Ｇ）ＡＴ（Ｔ／Ｃ／Ａ）ＡＡ（Ｔ／Ｃ／Ａ）ＡＡ（Ｇ／Ａ）（Ａ／Ｃ）Ｇ（Ｇ／Ｔ／Ｃ）ＡＴＣＧＡ（Ｃ／Ａ）ＡＡＣ－３’,然后进行３’,ＲＡ ＣＥ ＰＣＲ,扩增出约１ｋｂ大小的片段,序列分析表明该片段含有非常保守的ＭＡＤＳ ｂｏｘ。进而,利用５’ ＲＡＣＥ ＰＣＲ得到全长度的ｃＤＮＡ。该ｃＤＮＡ的核苷酸序列与ＣＵＭ１（黄瓜 ＭＡＤＳ ｂｏｘ ｇｅｎｅ１）同源性高达９７％。 ＣＵＭ１在接牵牛中过量表达可引起花萼变为心皮状和花瓣变为雄蕊,说明ＣＵＭ１为ＡＧ的同源基因。基于该基因与ＣＵＭ１序列上的高度同源,我们认为其为黄瓜的ＡＧ同源基因。该基因命名为ＣＭＢ１,基因银行登记号为ＡＦ２８６６４９。Ｓｏｕｔｈｅｒ     

Functional similarity analysis of human virus-encoded miRNAs

miRNAs are a class of small RNAs that regulate gene expression via RNA silencing machinery. Some viruses also encode miRNAs, contributing to the complex virus-host interactions. A better understanding of viral miRNA functions would be useful in designing new preventive strategies for treating diseases induced by viruses. To meet the challenge for how viruses module host gene expression by their encoded miRNAs, we measured the functional similarities among human viral miRNAs by using a method we reported previously. Higher order functions regulated by viral miRNAs were also identified by KEGG pathway analysis on their targets. Our study demonstrated the biological processes involved in virus-host interactions via viral miRNAs. Phylogenetic analysis suggested that viral miRNAs have distinct evolution rates compared with their corresponding genome.     

A dictionary on microRNAs and their putative target pathways

While in the last decade mRNA expression profiling was among the most popular research areas, over the past years the study of non-coding RNAs, especially microRNAs (miRNAs), has gained increasing interest. For almost 900 known human miRNAs hundreds of pretended targets are known. However, there is only limited knowledge about putative systemic effects of changes in the expression of miRNAs and their regulatory influence. We determined for each known miRNA the biochemical pathways in the KEGG and TRANSPATH database and the Gene Ontology categories that are enriched with respect to its target genes. We refer to these pathways and categories as target pathways of the corresponding miRNA. Investigating target pathways of miRNAs we found a strong relation to disease-related regulatory pathways, including mitogen-activated protein kinase (MAPK) signaling cascade, Transforming growth factor (TGF)-beta signaling pathway or the p53 network. Performing a sophisticated analysis of differentially expressed genes of 13 cancer data sets extracted from gene expression omnibus (GEO) showed that targets of specific miRNAs were significantly deregulated in these sets. The respective miRNA target analysis is also a novel part of our gene set analysis pipeline GeneTrail. Our study represents a comprehensive theoretical analysis of the relationship between miRNAs and their predicted target pathways. Our target pathways analysis provides a ‘miRNA-target pathway’ dictionary, which enables researchers to identify target pathways of differentially regulated miRNAs.     

Rapid and Efficient Isolation of High-Quality Small RNAs from Recalcitrant Plant Species Rich in Polyphenols and Polysaccharides

Small RNAs, including microRNAs (miRNAs) and small interfering RNAs (siRNAs), are important regulators of plant development and gene expression. The acquisition of high-quality small RNAs is the first step in the study of its expression and function analysis, yet the extraction method of small RNAs in recalcitrant plant tissues with various secondary metabolites is not well established, especially for tropical and subtropical plant species rich in polysaccharides and polyphenols. Here, we developed a simple and efficient method for high quality small RNAs extraction from recalcitrant plant species. Prior to RNA isolation, a precursory step with a CTAB-PVPP buffer system could efficiently remove compounds and secondary metabolites interfering with RNAs from homogenized lysates. Then, total RNAs were extracted by Trizol reagents followed by a differential precipitation of high-molecular-weight (HMW) RNAs using polyethylene glycol (PEG) 8000. Finally, small RNAs could be easily recovered from supernatant by ethanol precipitation without extra elimination steps. The isolated small RNAs from papaya showed high quality through a clear background on gel and a distinct northern blotting signal with miR159a probe, compared with other published protocols. Additionally, the small RNAs extracted from papaya were successfully used for validation of both predicted miRNAs and the putative conserved tasiARFs. Furthermore, the extraction method described here was also tested with several other subtropical and tropical plant tissues. The purity of the isolated small RNAs was sufficient for such applications as end-point stem-loop RT-PCR and northern blotting analysis, respectively. The simple and feasible extraction method reported here is expected to have excellent potential for isolation of small RNAs from recalcitrant plant tissues rich in polyphenols and polysaccharides.     

急性髓系白血病mRNA与非编码RNA差异表达谱及CeRNA调控网络分析

利用GEO数据库(gene expression omnibus database)通过生物信息学分析方法探讨急性髓系白血病(acute myelogenous leukemia,AML)的发病机制。检索GEO数据库中AML相关芯片数据集GSE142698、GSE142699和GSE96535。利用GEO2R分析得到差异mRNAs、miRNAs以及差异lncRNAs。利用在线生物信息学分析工具DAVID对差异mRNAs进行GO富集分析和KEGG通路分析。利用miRWalk数据库预测AML相关miRNAs的靶向mRNAs,利用Spongescan数据库预测AML相关miRNAs的靶向lncRNAs,构建lncRNA-miRNA-mRNA竞争性内源RNA （competing endogenous RNA,ceRNA）调控网络。共筛选出29个显著差异mRNAs、70个显著差异miRNAs和20 005个显著差异lncRNAs。GO富集分析和KEGG通路分析显示,差异表达基因主要涉及蛋白磷酸化、细胞分裂、细胞增殖的负调控、基因表达的正向调节、周期蛋白依赖的丝氨酸/苏氨酸激酶活性的调节等生物过程以及细胞周期、细胞衰老、癌症通路、PI3K-Akt通路等信号通路。将miRWalk数据库预测的靶向mRNAs与差异mRNAs取交集,Spongescan数据库预测的靶向lncRNAs与差异lncRNAs取交集,分别确定了25个mRNAs、6个lncRNAs参与AML相关ceRNA调控网络的构建。结果表明,lncRNAs可能作为关键的ceRNA,通过调控miRNA和相关靶基因参与AML的发生与发展,研究结果为AML诊断和治疗的分子生物学研究提供了新的依据。     

肾透明细胞癌预后相关miRNAs的生物信息学分析

目的寻找可作为肾透明细胞癌（ccRCC）生物标志物的miRNA,以及ccRCC与正常组织间miRNA差异表达情况。 方法利用TCGA数据库下载ccRCC中miRNA表达数据,分析肿瘤与正常组织间差异表达miRNA。使用Kaplan-Meier曲线对患者进行生存分析,筛选出表达情况与临床预后相关的miRNA。通过生物信息学对miRNA的靶基因进行预测,然后运用FunRich软件和ClueGO对靶基因进行GO和KEGG富集分析。 结果通过TCGA数据库分析发现,ccRCC较正常组织差异表达miRNA共54个,其中上调33个,下调21个。通过生存分析发现hsa-miR-21和hsa-miR-155与患者预后相关,P≤0.05。进一步通过Perl软件在Targetscan、miRDB、miRTarBase、miRPath这四个数据库中预测miRNA靶基因并将结果取交集,共发现129个靶基因。GO和KEGG分析结果表明,这些靶基因主要与转录因子活性、信号转导以及FoxO、TNF等信号通路密切相关。 结论通过生物信息学分析发现了ccRCC与正常组织的差异表达miRNA;其中hsa-miR-21和hsa-miR-155与患者总体生存率相关,并通过调控靶基因参与相关的信号通路进而影响ccRCC的发生发展进程,提示hsa-miR-21和hsa-miR-155可能是ccRCC潜在的生物标志物。     

Genome‐wide identification and characterization of microRNAs responding to ABA and GA in maize embryos during seed germination

• MicroRNAs (miRNAs) are an important class of non‐coding small RNAs that regulate the expression of target genes through mRNA cleavage or translational inhibition. Previous studies have revealed their roles in regulating seed dormancy and germination in model plants such as Arabidopsis thaliana, rice (Oryza sativa) and maize (Zea mays). However, the miRNA response to exogenous gibberellic acid (GA) and abscisic acid (ABA) during seed germination in maize has yet to be explored.
• In this study, small RNA libraries were generated and sequenced from maize embryos treated with GA, ABA or double‐distilled water as control.
• A total of 247 miRNAs (104 known and 143 novel) were identified, of which 45 known and 53 novel miRNAs were differentially expressed in embryos in the different treatment groups. In total, 74 (37 up‐regulated and 37 down‐regulated) and 55 (23 up‐regulated and 32 down‐regulated) miRNAs were expressed in response to GA and to ABA, respectively, and a total of 18 known and 38 novel miRNAs displayed differential expression between the GA‐ and ABA‐treated groups. Using bioinformatics tools, we predicted the target genes of the differentially expressed miRNAs. Using GO enrichment and KEGG pathway analysis of these targets, we showed that miRNAs differentially expressed in our samples affect genes encoding proteins involved in the peroxisome, ribosome and plant hormonal signalling pathways.
• Our results indicate that miRNA‐mediated gene expression influences the GA and ABA signalling pathways during seed germination.

     

Deregulated microRNAs in CD4+ T cells from individuals with latent tuberculosis versus active tuberculosis

The mechanisms of latent tuberculosis (TB) infection remain elusive. Roles of microRNA (miRNA) have been highlighted in pathogen–host interactions recently. To identify miRNAs involved in the immune response to TB, expression profiles of miRNAs in CD4⁺ T cells from patients with latent TB, active TB and healthy controls were investigated by microarray assay and validated by RT‐qPCR. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were used to analyse the significant functions and involvement in signalling pathways of the differentially expressed miRNAs. To identify potential target genes for miR‐29, interferon‐γ (IFN‐γ) mRNA expression was measured by RT‐qPCR. Our results showed that 27 miRNAs were deregulated among the three groups. RT‐qPCR results were generally consistent with the microarray data. We observed an inverse correlation between miR‐29 level and IFN‐γ mRNA expression in CD4⁺ T cells. GO and KEGG pathway analysis showed that the possible target genes of deregulated miRNAs were significantly enriched in mitogen‐activated protein kinase signalling pathway, focal adhesion and extracellular matrix receptor interaction, which might be involved in the transition from latent to active TB. In all, for the first time, our study revealed that some miRNAs in CD4⁺ T cells were altered in latent and active TB. Function and pathway analysis highlighted the possible involvement of miRNA‐deregulated mRNAs in TB. The study might help to improve understanding of the relationship between miRNAs in CD4⁺ T cells and TB, and laid an important foundation for further identification of the underlying mechanisms of latent TB infection and its reactivation.