A Systematic Analysis of miRNA Transcriptome in Marek’s Disease Virus-Induced Lymphoma Reveals Novel and Differentially Expressed miRNAs

Marek’s disease is a lymphoproliferative neoplastic disease of the chicken, which poses a serious threat to poultry health. Marek’s disease virus (MDV)-induced T-cell lymphoma is also an excellent biomedical model for neoplasia research. Recently, miRNAs have been demonstrated to play crucial roles in mediating neoplastic transformation. To investigate host miRNA expression profiles in the tumor transformation phase of MDV infection, we performed deep sequencing in two MDV-infected samples (tumorous spleen and MD lymphoma from liver), and two non-infected controls (non-infected spleen and lymphocytes). In total, 187 and 16 known miRNAs were identified in chicken and MDV, respectively, and 17 novel chicken miRNAs were further confirmed by qPCR. We identified 28 down-regulated miRNAs and 11 up-regulated miRNAs in MDV-infected samples by bioinformatic analysis. Of nine further tested by qPCR, seven were verified. The gga-miR-181a, gga-miR-26a, gga-miR-221, gga-miR-222, gga-miR-199*, and gga-miR-140* were down-regulated, and gga-miR-146c was up-regulated in MDV-infected tumorous spleens and MD lymphomas. In addition, 189 putative target genes for seven differentially expressed miRNAs were predicted. The luciferase reporter gene assay showed interactions of gga-miR-181a with MYBL1, gga-miR-181a with IGF2BP3, and gga-miR-26a with EIF3A. Differential expression of miRNAs and the predicted targets strongly suggest that they contribute to MDV-induced lymphomagenesis.     

Analysis of deep sequencing exosome-microRNA expression profile derived from CP-II reveals potential role of gga-miRNA-451 in inflammation

Mycoplasma gallisepticum (MG) can cause chronic respiratory disease (CRD) in chickens. While several studies have reported the inflammatory functions of microRNAs during MG infection, the mechanism by which exosomal miRNAs regulate MG-induced inflammation remains to be elucidated. The expression of exosome-microRNA derived from MG-infected chicken type II pneumocytes (CP-II) was screened, and the target genes and function of differentially expressed miRNAs (DEGs) were predicted. To verify the role of exosomal gga-miR-451, Western blot, ELISA and RT-qPCR were used in this study. The results showed that a total of 722 miRNAs were identified from the two exosomal small RNA (sRNA) libraries, and 30 miRNAs (9 up-regulated and 21 down-regulated) were significantly differentially expressed. The target miRNAs were significantly enriched in the treatment group, such as cell cycle, Toll-like receptor signalling pathway and MAPK signalling pathway. The results have also confirmed that gga-miR-451-absent exosomes derived from MG-infected CP-II cells increased inflammatory cytokine production in chicken fibroblast cells (DF-1), and wild-type CP-II cell–derived exosomes displayed protective effects. Collectively, our work suggests that exosomes from MG-infected CP-II cells alter the dynamics of the DF-1 cells, and may contribute to pathology of the MG infection via exosomal gga-miR-451 targeting YWHAZ involving in inflammation.     

携带与非携带松材线虫的松墨天牛miRNA表达谱比较分析

【目的】松墨天牛是松树的重要蛀干害虫,也是林业重大外来入侵种松材线虫的媒介昆虫。虽然松墨天牛和松材线虫互作的化学生态和分子进化机制受到人们的广泛关注,但miRNA等表观遗传因子在天牛—线虫互作中的作用未见报道。【方法】使用illumina HiSeq 2000平台进行miRNA高通量测序,得到4个携带线虫的天牛miRNA库和4个未携带线虫的天牛miRNA库。此外,对鉴定出的miRNA进行了差异表达分析,并对这些miRNA的靶基因进行了GO注释和KEGG通路富集分析。【结果】在携带线虫的天牛表皮、脂肪体、中肠和气管样本中分别鉴定出780、802、617和762个miRNA;在未携带松材线虫的天牛的不同组织样本中分别鉴定出784、723、713和837个miRNA。在携带松材线虫的松墨天牛中,某些已知miRNA表达量会显著升高,如miR-14、miR279和miR-312等。差异表达miRNA的功能大多指向代谢、免疫等方面。【结论】miRNA在松墨天牛和松材线虫互作中起着重要的调控作用。     

中华蜜蜂幼虫肠道响应球囊菌胁迫的microRNA应答分析

【目的】蜜蜂球囊菌(Ascosphaera apis,简称球囊菌)是一种能够侵染中华蜜蜂(Apis cerana cerana,简称中蜂)幼虫的致死性真菌病原。微小RNA(microRNA,miRNA)可通过在转录后水平靶向抑制或降解mRNA而参与宿主与病原互作过程。本研究旨在对球囊菌胁迫的中蜂6日龄幼虫肠道的差异表达miRNA(DEmiRNA)及其靶基因进行深入分析,进而揭示DEmiRNA在中蜂响应球囊菌胁迫应答过程中的作用。【方法】利用Illumina MiSeq平台对正常及球囊菌胁迫的中蜂6日龄幼虫肠道(AcCK和AcT)进行测序,通过相关生物信息学软件预测DEmiRNA及其靶基因。通过Blast将靶基因注释到GO和KEGG数据库。利用Cytoscape软件构建DEmiRNA与其靶mRNA的调控网络。通过Stem-loop RT-PCR和qPCR验证测序数据的可靠性。【结果】本研究共预测出537个miRNA,其长度分布介于16–35 nt之间,且不同长度的miRNA首位碱基偏向性差异明显。通过Stem-loop RT-PCR证实了10个novel miRNA的表达。AcCK vs AcT比较组共有54个DEmiRNA,包含31个上调和23个下调miRNA,可分别靶向结合6170和8199个靶基因。GO分类结果显示上调和下调miRNA的靶基因分别涉及47和47个条目,富集基因数最多的皆为结合细胞进程和催化活性。KEGG代谢通路(pathway)富集分析结果表明上调和下调miRNA的靶基因分别富集在134和126条pathway,富集基因数最多的均为内吞作用和内质网中的蛋白质加工。调控网络分析结果表明,DEmiRNA及其靶mRNA形成十分复杂的调控关系;31个DEmiRNA可靶向结合51个与泛素介导的蛋白水解相关的mRNA,18个DEmiRNA可靶向结合14个与Jak-STAT信号通路相关的mRNA;miR-1277-x、miR-26-x、miR-27-y、miR-30-x、miR-6052-x等16个miRNA共同参与了上述两条免疫通路的调控。最后,随机挑选3个DEmiRNA进行qPCR验证,结果证明了测序数据的可靠性。【结论】本研究提供了中蜂幼虫肠道在球囊菌胁迫后期的miRNA的表达谱和差异表达信息,揭示了球囊菌与宿主之间在miRNA组学水平存在复杂的互作。miR-6052-x和miR-1277-x作为调控网络的核心可能通过影响细胞凋亡参与宿主的免疫防御,miR-26-x和miR-30-x可能通过调控Jak-STAT信号通路参与宿主的胁迫应答。本研究筛选出的关键DEmiRNA有望作为治疗白垩病的分子靶标。     

鸡microRNA的研究进展

microRNA (miRNA)是参与基因转录后调控的一类重要的非编码小RNA分子.以下简要总结鸡miRNA的数量与染色体分布,同时概述了鸡miRNA对免疫、胚胎发育和病毒感染的调控作用,最后对鸡miRNA的应用前景也进行了简单探讨,以期为miRNA在禽业生产中的深入研究和应用提供参考.     

Chicken gga-miR-130a targets HOXA3 and MDFIC and inhibits Marek’s disease lymphoma cell proliferation and migration

Marek’s disease (MD) is an infectious disease of chickens caused by MD virus (MDV), which is a herpesvirus that initiates tumor formation. Studies have indicated that microRNAs (miRNAs) are linked with the development of cancers or tumors. Previously, gga-miR-130a was discovered downregulated in MDV-infected tissues. Here, we aimed to explore the further function of gga-miR-130a in MD. The expression of gga-miR-130a in MDV-infected and uninfected spleens was detected by quantitative real-time PCR (qRT-PCR). Subsequently, proliferation and migration assays of MDV-transformed lymphoid cells (MSB1) were carried out by transfecting gga-miR-130a. The target genes of gga-miR-130a were predicted using TargetScan and miRDB and clustered through Gene Ontology analysis. The target genes were validated by western blot, qRT-PCR, and a dual luciferase reporter assay. Our results show that the expression of gga-miR-130a was reduced in MDV-infected spleens. Gga-miR-130a showed an inhibitory effect on MSB1 cell proliferation and migration. Two target genes, homeobox A3 (HOXA3) and MyoD family inhibitor domain containing (MDFIC), were predicted and clustered to cell proliferation. Results indicate that gga-miR-130a regulates HOXA3 and MDFIC at the protein level but not at the mRNA level. Moreover, the gga-miR-130a binding sites of two target genes have been confirmed. We conclude that gga-miR-130a can arrest MSB1 cell proliferation and migration, and target HOXA3 and MDFIC, which are both involved in the regulation of cell proliferation. Collectively, gga-miR-130a plays a critical role in the tumorigenesis associated with chicken MD.     

Expression patterns and miRNA regulation of DNA methyltransferases in chicken primordial germ cells

DNA methylation is widespread in most species, from bacteria to mammals, and is crucial for genomic imprinting, gene expression, and embryogenesis. DNA methylation occurs via two major classes of enzymatic reactions: maintenance-type methylation catalyzed by DNA (cytosine-5-)-methyltransferase (DNMT) 1, and de novo methylation catalyzed by DNMT 3 alpha (DNMT3A) and -beta (DNMT3B). The expression pattern and regulation of DNMT genes in primordial germ cells (PGCs) and germ line cells has not been sufficiently established in birds. Therefore, we employed bioinformatics, RT-PCR, real-time PCR, and in situ hybridization analyses to examine the structural conservation and conserved expression patterns of chicken DNMT family genes. We further examined the regulation of a candidate de novo DNA methyltransferase gene, cDNMT3B by cotransfection of cDNMT3B 3'UTR- and cDNMT3B 3'UTR-specific miRNAs through a dual fluorescence reporter assay. All cDNMT family members were differentially detected during early embryonic development. Of interest, cDNMT3B expression was highly detected in early embryos and in PGCs. During germ line development and sexual maturation, cDNMT3B expression was reestablished in a female germ cell-specific manner. In the dual fluorescence reporter assay, cDNMT3B expression was significantly downregulated by four miRNAs: gga-miR-15c (25.82%), gga-miR-29b (30.01%), gga-miR-383 (30.0%), and gga-miR-222 (31.28%). Our data highlight the structural conservation and conserved expression patterns of chicken DNMTs. The miRNAs investigated in this study may induce downregulation of gene expression in chicken PGCs and germ cells.     

固始鸡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调控鸡下丘脑发育的分子机制提供了有益线索.     

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.     

gga-miR-375 Plays a Key Role in Tumorigenesis Post Subgroup J Avian Leukosis Virus Infection

Avian leukosis is a neoplastic disease caused in part by subgroup J avian leukosis virus J (ALV-J). Micro ribonucleic acids (miRNAs) play pivotal oncogenic and tumour-suppressor roles in tumour development and progression. However, little is known about the potential role of miRNAs in avian leukosis tumours. We have found a novel tumour-suppressor miRNA, gga-miR-375, associated with avian leukosis tumorigenesis by miRNA microarray in a previous report. We have also previously studied the biological function of gga-miR-375; Overexpression of gga-miR-375 significantly inhibited DF-1 cell proliferation, and significantly reduced the expression of yes-associated protein 1 (YAP1) by repressing the activity of a luciferase reporter carrying the 3′-untranslated region of YAP1. This indicates that gga-miR-375 is frequently downregulated in avian leukosis by inhibiting cell proliferation through YAP1 oncogene targeting. Overexpression of gga-miR-375 markedly promoted serum starvation induced apoptosis, and there may be the reason why the tumour cycle is so long in the infected chickens. In vivo assays, gga-miR-375 was significantly downregulated in chicken livers 20 days after infection with ALV-J, and YAP1 was significantly upregulated 20 days after ALV-J infection (P<0.05). We also found that expression of cyclin E, an important regulator of cell cycle progression, was significantly upregulated (P<0.05). Drosophila inhibitor of apoptosis protein 1 (DIAP1), which is related to caspase-dependent apoptosis, was also significantly upregulated after infection. Our data suggests that gga-miR-375 may function as a tumour suppressor thereby regulating cancer cell proliferation and it plays a key role in avian leukosis tumorigenesis.     

侵染中华蜜蜂6日龄幼虫的蜜蜂球囊菌的微小RNA差异表达谱及调控网络

【目的】蜜蜂球囊菌(Ascosphaeraapis,简称球囊菌)是专性侵染蜜蜂幼虫的致死性真菌病原。MicroRNA(miRNA)作为一类重要的基因表达调控因子,能够广泛参与真菌及其宿主的相互作用过程。本研究通过比较分析球囊菌孢子(AaCK)和侵染中华蜜蜂(Apis cerana cerana,简称中蜂) 6日龄幼虫肠道内的球囊菌(AaT)的smallRNA(sRNA)组学数据对球囊菌的差异表达miRNA(differentiallyexpressed miRNA,DEmiRNA)、靶mRNA及二者间的调控网络进行全面解析,旨在揭示miRNA介导的球囊菌对中蜂幼虫的侵染机制。【方法】对于球囊菌侵染的中蜂6日龄幼虫肠道的small RNA-seq (sRNA-seq)数据,利用BLAST工具连续比对东方蜜蜂(Apiscerana)和球囊菌的参考基因组筛滤得到AaT的sRNA组学数据。分别将AaCK和AaT的sRNA组学数据比对miRBase数据库,对球囊菌侵染宿主前后miRNA的数量和结构特征进行分析。联用RNAhybrid+svm_light、Miranda和Tar...     

The Oncogenic MicroRNA OncomiR-21 Overexpressed during Marek's Disease Lymphomagenesis Is Transactivated by the Viral Oncoprotein Meq

Gallid herpesvirus 2 (GaHV-2) is an oncogenic herpesvirus that causes T lymphoma in chicken. GaHV-2 encodes a basic leucine zipper (bZIP) protein of the AP-1 family, Meq. Upon formation of homo- or heterodimers with c-Jun, Meq may modulate the expression of viral and cellular genes involved in lymphomagenesis. GaHV-2 also encodes viral microRNAs (miRNAs) involved in latency and apoptosis escape. However, little is known about cellular miRNA deregulation during the development of GaHV-2-associated lymphoma. We determined the cellular miRNA expression profiles of chickens infected with a very virulent strain (RB-1B) or a vaccine strain (CVI988) or noninfected. Among the most deregulated cellular miRNAs, we focused our efforts on gga-miR-21, which is upregulated during GaHV-2 infection. We mapped the gga-miR-21 promoter to the 10th intron of the TMEM49 gene and found it to be driven by AP-1- and Ets-responsive elements. We show here that the viral oncoprotein Meq binds to this promoter, thereby transactivating gga-miR-21 expression. We confirmed that this miRNA targets chicken programmed death cell 4 (PDCD4) and promotes tumor cell growth and apoptosis escape. Finally, gga-miR-21 was overexpressed only during infection with a very virulent strain (RB-1B) and not during infection with a nononcogenic strain (CVI988), providing further evidence for its role in GaHV-2 lymphomagenesis. Our data therefore suggest an additional role for Meq in GaHV-2-mediated lymphomagenesis through the induction of miR-21 expression.     

Novel Biomarkers for Non-functioning Invasive Pituitary Adenomas were Identified by Using Analysis of microRNAs Expression Profile

The microRNAs (miRNAs) are involved in multiple pathological processes among various types of tumors. However, the functions of miRNAs in benign brain tumors are largely unexplored. In order to explore the pathogenesis of the invasiveness in non-functional pituitary adenoma (NFPA), the miRNAs expression profile was analyzed between invasive and non-invasive non-functional pituitary adenoma by miRNAs microarray. Six most significant differentially expressed miRNAs were identified including four upregulated miRNAs hsa-miR-181b-5p, hsa-miR-181d, hsa-miR-191-3p, and hsa-miR-598 and two downregulated miRNAs hsa-miR-3676-5p and hsa-miR-383. The functions and corresponding signaling pathways of differentially expressed miRNAs were investigated by bioinformatics techniques, including Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis. The result of GO analysis indicates regulation of voltage-gated potassium channel activity, positive regulation of sodium ion transport, positive regulation of GTPase activity, negative regulation of Notch signaling pathway, etc. KEGG pathway reveals a series of biological processes, including prolactin signaling pathway, endocrine and other factor-regulated calcium reabsorption, fatty acid metabolism, neuroactive ligand-receptor interaction, etc. The miRNAs hsa-miR-181a-5p was verified by quantitative real-time PCR, and the expression level was in accordance with the microarray result. Our result can provide the evidence on featured miRNAs which play a prominent role in pituitary adenoma as effective biomarkers and therapeutic targets in the future.     

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.     

感染GCRV的草鱼miRNA测序与比较分析

为探索基于高通量筛选抗性主效miRNA功能分子的方法应用于抗性草鱼(Ctenopharyngodon idellus)的选育,研究借助高通量测序鉴定分析感染GCRV前后的草鱼头肾组织中miRNAs.测序结果共鉴定出821个成熟miRNAs,其中118个在GCRV攻毒组特异表达,82个为未攻毒组特有;差异分析结果显示,G...     

蜜蜂球囊菌的microRNA鉴定及其调控网络分析

【目的】本研究利用small RNA-seq技术对球囊菌的纯培养进行测序,对球囊菌的micro RNAs miRNAs)进行预测、鉴定和分析,进而构建miRNAs-mRNAs的调控网络。【方法】利用Illumina Hiseq Xten平台对球囊菌菌丝与孢子进行测序,通过相关生物信息学软件对球囊菌的miRNAs进行预测和分析,通过茎环(Stem-loop)PCR对部分miRNAs进行鉴定,利用Cytoskype软件构建miRNAs-mRNAs的调控网络。【结果】本研究共获得48268696条clean reads,预测出118个球囊菌的miRNAs,它们的长度分布介于18–25 nt之间,不同长度的mi RNA的首位碱基偏好性差异明显。Stem-loop PCR验证结果显示共有10个miRNAs能够扩增出符合预期的目的片段,说明多数miRNAs可能真实存在。共预测出6529个球囊菌miRNAs的靶基因,其中5725个能够注释到Nr、Swissprot、KOG、GO和KEGG数据库。进一步分析结果显示有24个靶基因注释在MAPK信号通路。Cytoskype软件分析结果显示球囊菌的miRNAs与mRNAs之间存在复杂的调控网络,绝大多数的miRNAs处于调控网络的内部且同时结合多个mRNAs。【结论】本研究率先对球囊菌的miRNAs及miRNAs-mRNAs调控网络进行全面分析,研究结果丰富了对球囊菌miRNAs的认识,为其基础生物学信息提供了有益补充,也为阐明球囊菌致病的分子机理打下了一定基础。