Genome-wide analysis of long noncoding RNA and mRNA profiles in PRRSV-infected porcine alveolar macrophages

Porcine reproductive and respiratory syndrome (PRRS), which is caused by PRRS virus (PRRSV), is one of the most globally devastating swine diseases. It is essential to develop new strategy to control PRRS via an understanding of mechanisms that PRRSV utilizes to interfere with the host's innate immunity. In this study, we deeply sequenced and analyzed long noncoding RNA (lncRNA) and mRNA expression profiles of the porcine alveolar macrophages (PAMs) after PRRSV infection. 126 lncRNAs and 753 mRNAs were differentially expressed between PRRSV-infected and control PAMs. The co-expressed genes of down-regulated lncRNAs were significantly enriched within NF-kappa B and toll-like receptor signaling pathways. Co-expression network analysis indicated that part of the dysregulated lncRNAs associated with the interferon-induced genes. These dysregulated lncRNAs may play an important role in the host's innate immune responses to PRRSV infection. However, further research is required to characterize the function of these lncRNAs.     

Genome-Wide Gene Expression Profiles in Lung Tissues of Pig Breeds Differing in Resistance to Porcine Reproductive and Respiratory Syndrome Virus

Porcine reproductive and respiratory syndrome (PRRS) caused by PRRS virus (PRRSV) is an infectious disease characterized by severe reproductive deficiency in pregnant sows, typical respiratory symptoms in piglets, and high mortality rate of piglets. In this study, we employed an Affymetrix microarray chip to compare the gene expression profiles of lung tissue samples from Dapulian (DPL) pigs (a Chinese indigenous pig breed) and Duroc×Landrace×Yorkshire (DLY) pigs after infection with PRRSV. During infection with PRRSV, the DLY pigs exhibited a range of clinical features that typify the disease, whereas the DPL pigs showed only mild signs of the disease. Overall, the DPL group had a lower percentage of CD4⁺ cells and lower CD4⁺/CD8⁺ratios than the DLY group (p<0.05). For both IL-10 and TNF-α, the DLY pigs had significantly higher levels than the DPL pigs (p<0.01). The DLY pigs have lower serum IFN-γ levels than the DPL pigs (p<0.01). The serum IgG levels increased slightly from 0 dpi to 7 dpi, and peaked at 14 dpi (p<0.0001). Microarray data analysis revealed 16 differentially expressed (DE) genes in the lung tissue samples from the DLY and DPL pigs (q≤5%), of which LOC100516029 and LOC100523005 were up-regulated in the PRRSV-infected DPL pigs, while the other 14 genes were down-regulated in the PRRSV-infected DPL pigs compared with the PRRSV-infected DLY pigs. The mRNA expression levels of 10 out of the 16 DE genes were validated by real-time quantitative RT-PCR and their fold change was consistent with the result of microarray data analysis. We further analyzed the mRNA expression level of 8 differentially expressed genes between the DPL and DLY pigs for both uninfected and infected groups, and found that TF and USP18 genes were important in underlying porcine resistance or susceptibility to PRRSV.     

参与调控意大利蜜蜂工蜂中肠基因表达的东方蜜蜂微孢子虫miRNA的组学解析及其调控网络

[目的]本研究结合前期已获得的miRNA和mRNA组学数据对东方蜜蜂微孢子虫Nosema ceranae 的差异表达 miRNA(differentially expressed miRNA,DEmiRNA)靶向意大利蜜蜂 Apis mellifera ligustica 工蜂中肠的 mRNA 和差异表达 mRNA(d...     

MicroRNA expression profile in RAW264.7 cells in response to Brucella melitensis infection

MicroRNA (miRNA) is small non-coding RNA with approximate 22 nt in length. Recent studies indicate that miRNAs play significant roles in pathogen-host interactions. Brucella organisms are Gram-negative facultative intracellular bacteria that cause Brucellosis. Brucella strains infect macrophages and establish chronic infection by altering host life activities including apoptosis and autophagy. Here, we report a comprehensive analysis of miRNA expression profiles in mock- and Brucella-infected RAW264.7 cells using high-throughput sequencing approach. In total, 344 unique miRNAs were co-expressed in the two libraries, in which 57 miRNAs were differentially expressed. Eight differentially expressed miRNAs with high abundance were subjected to further analysis. The GO enrichment analysis suggests that the putative target genes of these differentially expressed miRNAs are involved in apoptosis, autophagy and immune response. In particular, a total of 25 target genes are involved in regulating apoptosis and autophagy, indicating that these miRNAs may play important regulatory roles in the Brucella-host interactions. Furthermore, the interactions of miR-1981 and its target genes, Bcl-2 and Bid, were validated by luciferase assay. The results show that miR-1981 mimic up-regulated the luciferase activity of psiCHECK-2 Bcl-2 3' UTR, but the luciferase activity of psiCHECK-2 Bid 3' UTR was not changed significantly. Taken together, these data provide valuable framework on Brucella induced miRNA expression in RAW264.7 cells, and suggest that Brucella may establish chronic infection by regulating miRNA expression profile.     

Dysregulated genes and miRNAs in the apoptosis pathway in colorectal cancer patients

Apoptosis is genetically regulated and involves intrinsic and extrinsic pathways. We examined 133 genes within these pathways to identify whether they are expressed differently in colorectal carcinoma (CRC) and normal tissue (N?=?217) and if they are associated with similar differential miRNA expression. Gene expression data (RNA-Seq) and miRNA expression data (Agilent Human miRNA Microarray V19.0) were generated. We focused on dysregulated genes with a fold change (FC) of >?1.50 or <?0.67, that were significant after adjustment for multiple comparisons. miRNA:mRNA seed-region matches were determined. Twenty-three genes were significantly downregulated (FC?<?0.67) and 18 were significantly upregulated (FC?>?1.50). Of these 41 genes, 11 were significantly associated with miRNA differential expression. BIRC5 had the greatest number of miRNA associations (14) and the most miRNAs with a seed-region match (10). Four of these matches, miR-145-5p, miR-150-5p, miR-195-5p, and miR-650, had a negative beta coefficient. CSF2RB was associated with ten total miRNAs (five with a seed-region match, and one miRNA, miR-92a-3p, with a negative beta coefficient). Of the three miRNAs associated with CTSS, miR-20b-5p, and miR-501-3p, had a seed-region match and a negative beta coefficient between miRNA:mRNA pairs. Several miRNAs that were associated with dysregulated gene expression, seed-region matches, and negative beta coefficients also were associated with CRC-specific survival. Our data suggest that miRNAs could influence several apoptosis-related genes. BIRC5, CTSS, and CSF2R all had seed-region matches with miRNAs that would favor apoptosis. Our study identifies several miRNA associated with apoptosis-related genes, that if validated, could be important therapeutic targets.     

MiRNA expression profile and miRNA–mRNA integrated analysis (MMIA) during podocyte differentiation

The podocyte is a prominent cell type, which encases the capillaries of glomerulus. Podocyte-selective deletion of Dicer or Drosha was reported to induce proteinuria and glomerulosclerosis, suggesting the essential role of microRNA (miRNA) in podocytes for renal function. However, no comprehensive miRNA expression or miRNA–mRNA integrated analysis (MMIA) can be found during podocyte differentiation. Herein, miRNA and mRNA microarrays are presented, which were carried out in differentiated and undifferentiated mouse podocyte cell lines (MPC5). A total of 50 abnormal miRNAs (26 down-regulated and 24 up-regulated) were identified in differentiated and undifferentiated podocytes. Using MMIA, 80 of the 743 mRNAs (>twofold change) were predicted for potential crosstalk with 30 miRNAs of the 50 abnormal miRNAs. In addition, the gene ontology of mRNAs and the pathway analysis of miRNAs revealed a new potential-regulated network during podocyte differentiation. The expressions of three remarkably changed miRNAs (miR-34c, miR-200a and miR-467e) and four mRNAs (Runx1t1, Atp2a2, Glrp1, and Mmp15), were randomly chosen for further validation by the quantitative real-time polymerase chain reaction, and their expression trends were consistent with the microarray data. Reference searching was also conducted to confirm our data and to find potential new molecules and miRNA-target pairs involved in the podocyte differentiation. The dual luciferase reporter assay for miR-200a/GLRX and let-7b/ARL4D confirmed the prediction of MMIA. The results of this study provide a detailed integration of mRNA and miRNA during podocyte differentiation. The molecular integration mode will open up new perspectives for a better understanding of the mechanism during podocyte differentiation.     

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.     

意大利蜜蜂幼虫肠道在球囊菌胁迫后期的差异表达微小RNA及其靶基因分析

【目的】蜜蜂球囊菌Ascosphaera apis(简称球囊菌)是一种特异性侵染蜜蜂幼虫肠道的致死性真菌病原。微小RNA(microRNA,mi RNA)是一种在转录后水平对m RNA进行负调控的关键调控因子。本研究旨在对意大利蜜蜂Apis mellifera ligustica(简称意蜂)幼虫肠道在球囊菌胁迫后期的差异表达mi RNAs(DEmiRNAs)及其靶基因进行深入分析,为揭示DEmiRNAs在胁迫应答中的作用提供重要信息。【方法】利用small RNA-seq(s RNA-seq)技术对正常及球囊菌侵染的意蜂6日龄幼虫肠道(分别表示为Am CK和Am T)进行测序,通过相关生物信息学软件对DEmiRNAs进行预测和分析。利用TargetFinder软件预测DEmiRNAs的靶基因,然后利用Blast软件对靶基因进行GO和KEGG数据库的功能注释。通过Cytoscape软件构建DEmiRNAs与其靶m RNAs的调控网络。通过茎环实时荧光定量PCR(stem-loop RT-qPCR)验证测序数据的可靠性。【结果】Am CK和Am T的测序分别得到9 230 496和10 823 667条有效序列标签; Am CK vs Am T比较组中包含15个上调和6个下调mi RNAs,分别结合3 503和3 252个靶基因,它们可分别注释到40和39个GO terms以及104和99条代谢通路;进一步分析发现调控网络中的17个DEmiRNAs靶向结合116个与丝氨酸蛋白酶相关的m RNAs,14个DEmiRNAs靶向结合54个与泛素介导的蛋白水解相关的m RNAs。Stem-loop RT-qPCR验证结果显示随机选择的4个DEmiRNAs (mi R-251-x,mi R-9277-y,mi R-1672-x和mi R-4968-y)的表达量变化趋势与测序结果一致,表明测序数据真实可信。【结论】本研究率先对意蜂幼虫肠道在球囊菌胁迫后期的DEmiRNAs及其靶基因进行预测与分析,提供了mi RNAs的表达谱和差异表达信息。ame-miR-927b,mi R-429-y和mi R-8440-y等DEmiRNAs可能参与对丝氨酸蛋白酶的调控和对泛素介导的蛋白水解的调控,DEmiRNAs与靶基因之间存在复杂的调控网络关系。DEmiRNAs可参与到意蜂幼虫肠道与球囊菌间的互作。     

MicroRNA-target gene responses to root knot nematode (Meloidogyne incognita) infection in cotton (Gossypium hirsutum L.)

MicroRNAs (miRNAs) are a large class of small regulatory RNA molecules, however no study has been performed to elucidate the role of miRNAs in cotton (Gossypium hirsutum) response to the root knot nematode (RKN, Meloidogyne incognita) infection. We selected 28 miRNAs and 8 miRNA target genes to investigate the miRNA-target gene response to M. incognita infection. Our results show that RKN infection significantly affected the expression of several miRNAs and their targeted genes. After 10 days of RKN infection, expression fold changes on miRNA expressions ranged from down-regulated by 33% to upregulated by 406%; meanwhile the expression levels of miRNA target genes were 45.8% to 231%. Three miRNA-target pairs, miR159-MYB, miR319-TCP4 and miR167-ARF8, showed inverse expression patterns between gene targets and their corresponded miRNAs, suggesting miRNA-mediated gene regulation in cotton roots in response to RKN infection.     

星形细胞瘤生存预后相关miRNA-mRNA调控关系对筛选

利用GEO数据库中的芯片数据,筛选与星形细胞瘤生存预后相关的miRNA-mRNA调控关系对,为后续研究提供理论支持.下载芯片数据利用R语言进行差异表达分析,得到星形细胞瘤较正常组织表达显著改变的miRNA与mRNA;通过miRNA靶基因预测,将靶基因与差异表达mRNA取交集,明确mRNA与miRNA之间的关系;利用GE...     

Coordinated miRNA/mRNA expression profiles for understanding breed-specific metabolic characters of liver between Erhualian and large white pigs

MicroRNAs (miRNAs) are involved in the regulation of various metabolic processes in the liver, yet little is known on the breed-specific expression profiles of miRNAs in coordination with those of mRNAs. Here we used two breeds of male newborn piglets with distinct metabolic characteristics, Large White (LW) and Erhualian (EHL), to delineate the hepatic expression profiles of mRNA with microarray and miRNAs with both deep sequencing and microarray, and to analyze the functional relevance of integrated miRNA and mRNA expression in relation to the physiological and biochemical parameters. EHL had significantly lower body weight and liver weight at birth, but showed elevated serum levels of total cholesterol (TCH), high-density lipoprotein cholesterol (HDLC) and low-density lipoprotein cholesterol (LDLC), as well as higher liver content of cholesterol. Higher serum cortisol and lower serum insulin and leptin were also observed in EHL piglets. Compared to LW, 30 up-regulated and 18 down-regulated miRNAs were identified in the liver of EHL, together with 298 up-regulated and 510 down-regulated mRNAs (FDR<10%). RT-PCR validation of some differentially expressed miRNAs (DEMs) further confirmed the high-throughput data analysis. Using a target prediction algorithm, we found significant correlation between the up-regulated miRNAs and down-regulated mRNAs. Moreover, differentially expressed genes (DEGs), which are involved in proteolysis, were predicted to be mediated by DEMs. These findings provide new information on the miRNA and mRNA profiles in porcine liver, which would shed light on the molecular mechanisms underlying the breed-specific traits in the pig, and may serve as a basis for further investigation into the biological functions of miRNAs in porcine liver.     

MicroRNA profile analysis of a feline kidney cell line before and after infection with mink enteritis virus

MicroRNAs (miRNAs) are small regulatory RNAs that play a significant role in eukaryotes by targeting mRNAs for cleavage or translational repression. Recent studies have also shown them to be associated with cellular changes following viral infection. Mink enteritis virus (MEV) is one of the most important viral pathogens in the mink industry. To study the involvement of miRNAs in the MEV infection process, we used Illumina's ultrahigh throughput approach to sequencing miRNA libraries from the feline kidney (F81) cell line before and after infection with MEV. Using this bioinformatics approach we identified 196 known mammalian miRNA orthologs belonging to 152 miRNA families in F81 cells. Additionally, 97 miRNA*s of these miRNAs were detected. As well as known miRNAs, 384 and 398 novel miRNA precursor candidates were identified in uninfected and MEV-infected F81 cells respectively that have not been reported in other mammals. In MEV-infected cells 3 miRNAs were significantly down-regulated and 4 up-regulated including 3 significantly. The majority (12 of 16) of randomly selected miRNA expression profiles by qRT-PCR were consistent with those identified by deep sequencing. A total of 88 miRNAs were predicted to target interferon-associated genes; 6 appear to target the 3′UTR of MEV-specific receptor transferring receptor mRNAs; and 8 to target the MEV mRNA coding region. No miRNAs coded by MEV itself were detected.     

Analysis of miRNA-seq combined with gene expression profile reveals the complexity of salinity stress response in <Emphasis Type="Italic">Oryza sativa</Emphasis>

High salinity is a major abiotic stressor that affects crop productivity and quality. While proper seedling growth is critical for crop reproduction under high salinity stress. Nowadays, genes/miRNAs expression is used for studying salinity stress response in rice seedlings. However, analysis of miRNA combined with gene expression is rare. To this end, we used miRNA-seq and gene expression profile to ascertain 6335 genes (3276 genes up-regulated, 3059 genes down-regulated) and 126 miRNAs (47 miRNAs up-regulated, 79 miRNAs down-regulated) that respond to salinity stress in rice seedlings. We then used these 126 miRNAs (including the novel miRNA osa-Chr12_1506) to identify 121 differentially expressed predicted target genes. In addition, we identified 34 miRNA-target RNA pairs, consisting of 9 differentially expressed miRNAs with complementary expression patterns. Combined with previous studies, we proposed a simple model for the molecular mechanism of a 12-h salinity stress response in rice seedlings. The findings lead to a deeper understanding of the function of miRNAs and genes that respond to salinity, and contributed to the elucidation of the complex mechanisms activated by salinity stress.     

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.     

Comparison of microRNAs in the intramuscular adipose tissue from Jinhua and Landrace pigs

MicroRNA (miRNA) is critically involved in lipogenesis occurring in various body parts of humans and animals. In this study, to further investigate the role and distribution of miRNA in porcine intramuscular adipose tissue, small RNAs were extracted from Jinhua and Landrace pigs to identify the expression profiles of miRNAs. miRNA expression profiles revealed that 558 miRNAs including 287 known and 271 novel miRNAs were identified, and 220 of them showed differential expression in the pigs. Gene ontology and Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis suggested that the target genes of the differentially expressed miRNAs were involved in fatty metabolism. In conclusion, the current study reveals the active participation of miRNAs in the regulation of adipogenesis in the intramuscular adipose tissue of Jinhua and Landrace pigs.