Long intergenic non-coding RNAs (lincRNAs) are transcribed from non-coding DNA sequences. Studies have revealed that aberrant expressions of lincRNAs are associated with various types of cancers and neurological disorders. Marek's disease (MD) is a highly contagious T-cell lymphoid neoplasia of chicken induced by Marek's disease virus (MDV). In this study, we first identified and validated linc-GALMD3 highly expressed in MDV-infected CD4+ T cells by RNA-Seq and qRT-PCR. By RNA-Seq analysis in MDCC-MSB1 cells after loss of function of linc-GALMD3 by shRNA, we found that linc-GALMD3 could positively cis-regulate its downstream gga-miR-223 gene expression. In contrast, it could trans-regulate the 748 differentially expressed genes (FDR < 0.01) that were mainly enriched into mitochondrial structure and cell cycle processes using GO analysis. Of these, the most significantly expressed gene EPYC might cause iris lesion in MD. The other eight genes, NDUFA4, NDUFB6, NDUFV1, NDUFS8, SDHB, UQCRC1, UQCRC2, and COX7A2, actively participated in oxidative phosphorylation in mitochondrial dysfunction and cell death. Most importantly, we found that the MDV replication was repressed when linc-GALMD3 was knocked down in CEF cells. Our results suggested that linc-GALMD3 might be a critical regulator in chicken MD and could be used as a candidate-promising mark for MD prevention, diagnosis, and treatment.

To explore the role of circular RNA (circRNA) circZFR in tumorigenic capacity of lung cancer (LC).

RNA-seq technology has become an important tool for quantifying the gene and transcript expression in transcriptome study. The two major difficulties for the gene and transcript expression quantification are the read mapping ambiguity and the overdispersion of the read distribution along reference sequence. Many approaches have been proposed to deal with these difficulties. A number of existing methods use Poisson distribution to model the read counts and this easily splits the counts into the contributions from multiple transcripts. Meanwhile, various solutions were put forward to account for the overdispersion in the Poisson models. By checking the similarities among the variation patterns of read counts for individual genes, we found that the count variation is exon-specific and has the conserved pattern across the samples for each individual gene. We introduce Gamma-distributed latent variables to model the read sequencing preference for each exon. These variables are embedded to the rate parameter of a Poisson model to account for the overdispersion of read distribution. The model is tractable since the Gamma priors can be integrated out in the maximum likelihood estimation. We evaluate the proposed approach, PGseq, using four real datasets and one simulated dataset, and compare its performance with other popular methods. Results show that PGseq presents competitive performance compared to other alternatives in terms of accuracy in the gene and transcript expression calculation and in the downstream differential expression analysis. Especially, we show the advantage of our method in the analysis of low expression.

This study aimed to investigate the regulatory role of differentially-expressed circular RNAs (circRNAs) in mouse cardiomyocytes during doxorubicin (DOX)-induced cardiotoxicity.

The quality and quantity of animal meat are closely related to the development of skeletal muscle, which, in turn, is determined by myogenic cells, including myoblasts and skeletal muscle satellite cells (SMSCs). Circular RNA, an endogenous RNA derivative formed through specific reverse splicing in mRNA precursors, has the potential to influence muscle development by binding to miRNAs or regulating gene expression involved in muscular growth at the transcriptional level. Previous high-throughput sequencing of circRNA in chicken liver tissue revealed a circular transcript, circIGF2BP3, derived from the gene encoding insulin-like growth factor 2 mRNA binding protein 3 (IGF2BP3). In this study, we confirmed the presence of the natural circular molecule of circIGF2BP3 through an RNase R enzyme tolerance assay. RT-qPCR results showed high circIGF2BP3 expression in the pectoral and thigh muscles of Yuexi frizzled feather chickens at embryonic ages 14 and 18, as well as at 7 weeks post-hatch. Notably, its expression increased during embryonic development, followed by a rapid decrease after birth. As well as using RT-qPCR, Edu, CCK-8, immunofluorescence, and Western blot techniques, we demonstrated that overexpressing circIGF2BP3 could promote the proliferation and differentiation of chicken primary myoblasts through upregulating genes such as proliferating cell nuclear antigen (PCNA), cyclin D1 (CCND1), cyclin E1 (CCNE1), cyclin dependent kinase 2 (CDK2), myosin heavy chain (MyHC), myoblast-determining 1 (MyoD1), myogenin (MyoG), and Myomaker. In conclusion, circIGF2BP3 promotes the proliferation and differentiation of myoblasts in chickens. This study establishes a foundation for further investigation into the biological functions and mechanisms of circIGF2BP3 in myoblasts proliferation and differentiation.

Natural antisense transcripts (NATs) are widely present in mammalian genomes and act as pivotal regulator molecules of gene expression. However, studies on NATs in the chicken are relatively rare. We identified a novel antisense transcript in the chicken, designated GHR-AS-EST, transcribed from the growth hormone receptor (GHR) locus, which encodes a well-known regulatory molecule of muscle development and fat deposition. GHR-AS-EST is predominantly expressed in the chicken liver and muscle tissues. GHR-AS-EST sequence conservation among vertebrates is weak. GHR-AS-EST forms an RNA-RNA duplex with GHBP to increase its stability, and regulates the expression of GHR sense transcripts at both the mRNA and protein levels. Further, GHR-AS-EST promotes cell proliferation by stimulating the expression of signaling factors in the JAK2/STAT pathway, and contributes to fat deposition via downregulating the expression of signaling factors in the JAK2/SOCS pathway in LMH hepatocellular carcinoma cells. We expect that the discovery of a NAT for a regulatory gene associated with cell proliferation and lipolysis will further our understanding of the molecular regulation of both muscle development and fat deposition.

Lung adenocarcinomas (LUAD) arise from precancerous lesions such as atypical adenomatous hyperplasia, which progress into adenocarcinoma in situ and minimally invasive adenocarcinoma, then finally into invasive adenocarcinoma. The cellular heterogeneity and molecular events underlying this stepwise progression remain unclear. In this study, we perform single-cell RNA sequencing of 268,471 cells collected from 25 patients in four histologic stages of LUAD and compare them to normal cell types. We detect a group of cells closely resembling alveolar type 2 cells (AT2) that emerged during atypical adenomatous hyperplasia and whose transcriptional profile began to diverge from that of AT2 cells as LUAD progressed, taking on feature characteristic of stem-like cells. We identify genes related to energy metabolism and ribosome synthesis that are upregulated in early stages of LUAD and may promote progression. MDK and TIMP1 could be potential biomarkers for understanding LUAD pathogenesis. Our work shed light on the underlying transcriptional signatures of distinct histologic stages of LUAD progression and our findings may facilitate early diagnosis.

Ticks act as important vectors of infectious agents, and several emerging tick-borne viruses have recently been identified to be associated with human diseases in northeastern China. However, little is known about the tick virome in northeastern China.

Liver X receptors (LXRs) are transcriptional regulators of cellular and systemic cholesterol homeostasis. Under conditions of excess cholesterol, LXR activation induces the expression of several genes involved in cholesterol efflux, facilitates cholesterol esterification by promoting fatty acid synthesis, and inhibits cholesterol uptake by the low-density lipoprotein receptor. The fact that sterol content is maintained in a narrow range in most cell types and in the organism as a whole suggests that extensive crosstalk between regulatory pathways must exist. However, the molecular mechanisms that integrate LXRs with other lipid metabolic pathways are incompletely understood. Here we show that ligand activation of LXRs in mouse liver not only promotes cholesterol efflux, but also simultaneously inhibits cholesterol biosynthesis. We further identify the long non-coding RNA LeXis as a mediator of this effect. Hepatic LeXis expression is robustly induced in response to a Western diet (high in fat and cholesterol) or to pharmacological LXR activation. Raising or lowering LeXis levels in the liver affects the expression of genes involved in cholesterol biosynthesis and alters the cholesterol levels in the liver and plasma. LeXis interacts with and affects the DNA interactions of RALY, a heterogeneous ribonucleoprotein that acts as a transcriptional cofactor for cholesterol biosynthetic genes in the mouse liver. These findings outline a regulatory role for a non-coding RNA in lipid metabolism and advance our understanding of the mechanisms that coordinate sterol homeostasis.

A long and ever-expanding roster of small (∼20-30 nucleotides) RNAs has emerged during the last decade, and most can be subsumed under the three main headings of microRNAs (miRNAs), Piwi-interacting RNAs (piRNAs), and short interfering RNAs (siRNAs). Among the three categories, miRNAs is the most quickly expanded group. The most recent number of identified miRNAs is 16,772 (Sanger miRbase, April 2011). However, there are insufficient publications on their primary forms, and no tissue-specific small RNAs precursors have been reported in the epididymis. Here, we report the identification in rats of an epididymis-specific, chimeric, noncoding RNA that is spliced from two different chromosomes (chromosomes 5 and 19), which we named HongrES2. HongrES2 is a 1.6 kb mRNA-like precursor that gives rise to a new microRNA-like small RNA (mil-HongrES2) in rat epididymis. The generation of mil-HongrES2 is stimulated during epididymitis. An epididymis-specific carboxylesterase named CES7 had 100% cDNA sequence homology at the 3'end with HongrES2 and its protein product could be downregulated by HongrES2 via mil-HongrES2. This was confirmed in vivo by initiating mil-HongrES2 over-expression in rats and observing an effect on sperm capacitation.

Animal growth and development are regulated by neural and endocrine growth axes, in which cell proliferation plays key roles. Recently, many research showed that circular RNAs were involved in hepatocyte and myoblast proliferation. Previously, we identified a circular RNA derived from the chicken GHR gene, named circGHR. However, the function of circGHR is unclear. The objective of this study was to investigate circGHR expression pattern and its roles in cell proliferation. Results indicated that circGHR was a closed-loop structure molecule, and it was richer in the nucleus of hepatocytes and myoblast. Real-time PCR showed that circGHR was increased from E13 to the 7th week in the liver but decreased in the thigh and breast muscle. The CCK-8 assay displayed that circGHR promoted cell proliferation. Simultaneously, the biomarker genes PCNA, CCND1, and CDK2 and the linear transcripts GHR and GHBP were upregulated when circGHR was overexpressed. Altogether, these data exhibited that circGHR could promote cell proliferation possibly by regulating GHR mRNA and GHBP expression.

Non-covalent complexes of glycolytic enzymes, called metabolons, were postulated in the 1970s, but the concept has been controversial. Here we show that a c-Myc-responsive long noncoding RNA (lncRNA) that we call glycoLINC (gLINC) acts as a backbone for metabolon formation between all four glycolytic payoff phase enzymes (PGK1, PGAM1, ENO1, and PKM2) along with lactate dehydrogenase A (LDHA). The gLINC metabolon enhances glycolytic flux, increases ATP production, and enables cell survival under serine deprivation. Furthermore, gLINC overexpression in cancer cells promotes xenograft growth in mice fed a diet deprived of serine, suggesting that cancer cells employ gLINC during metabolic reprogramming. We propose that gLINC makes a functional contribution to cancer cell adaptation and provide the first example of a lncRNA-facilitated metabolon.

RBM15, an RNA binding protein, determines cell-fate specification of many tissues including blood. We demonstrate that RBM15 is methylated by protein arginine methyltransferase 1 (PRMT1) at residue R578, leading to its degradation via ubiquitylation by an E3 ligase (CNOT4). Overexpression of PRMT1 in acute megakaryocytic leukemia cell lines blocks megakaryocyte terminal differentiation by downregulation of RBM15 protein level. Restoring RBM15 protein level rescues megakaryocyte terminal differentiation blocked by PRMT1 overexpression. At the molecular level, RBM15 binds to pre-messenger RNA intronic regions of genes important for megakaryopoiesis such as GATA1, RUNX1, TAL1 and c-MPL. Furthermore, preferential binding of RBM15 to specific intronic regions recruits the splicing factor SF3B1 to the same sites for alternative splicing. Therefore, PRMT1 regulates alternative RNA splicing via reducing RBM15 protein concentration. Targeting PRMT1 may be a curative therapy to restore megakaryocyte differentiation for acute megakaryocytic leukemia.

Penile cancer (PeCa) is a relatively rare tumor entity but possesses higher morbidity and mortality rates especially in developing countries. To date, the concrete pathogenic signaling pathways and core machineries involved in tumorigenesis and progression of PeCa remain to be elucidated. Several studies suggested miRNAs, which modulate gene expression at posttranscriptional level, were frequently mis-regulated and aberrantly expressed in human cancers. However, the miRNA profile in human PeCa has not been reported before. In this present study, the miRNA profile was obtained from 10 fresh penile cancerous tissues and matched adjacent non-cancerous tissues via next-generation sequencing. As a result, a total of 751 and 806 annotated miRNAs were identified in normal and cancerous penile tissues, respectively. Among which, 56 miRNAs with significantly different expression levels between paired tissues were identified. Subsequently, several annotated miRNAs were selected randomly and validated using quantitative real-time PCR. Compared with the previous publications regarding to the altered miRNAs expression in various cancers and especially genitourinary (prostate, bladder, kidney, testis) cancers, the most majority of deregulated miRNAs showed the similar expression pattern in penile cancer. Moreover, the bioinformatics analyses suggested that the putative target genes of differentially expressed miRNAs between cancerous and matched normal penile tissues were tightly associated with cell junction, proliferation, growth as well as genomic instability and so on, by modulating Wnt, MAPK, p53, PI3K-Akt, Notch and TGF-β signaling pathways, which were all well-established to participate in cancer initiation and progression. Our work presents a global view of the differentially expressed miRNAs and potentially regulatory networks of their target genes for clarifying the pathogenic transformation of normal penis to PeCa, which research resource also provides new insights into future investigations aimed to explore the in-depth mechanisms of miRNAs and other small RNAs including piRNAs in penile carcinogenesis regulation and effective target-specific theragnosis.

MicroRNAs (miRNAs) are considered to be involved in the pathogenic initiation and progression of chronic nonbacterial prostatitis (CNP); however, the comprehensive expression profile of dysregulated miRNAs, relevant signaling pathways, and core machineries in CNP have not been fully elucidated. In the current research, CNP rat models were established through the intraprostatic injection of carrageenan into the prostate. Then, next-generation sequencing was performed to explore the miRNA expression profile in CNP. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) bioinformatical analyses were conducted to reveal the enriched biological processes, molecular functions, and cellular components and signaling pathways. As a result, 1224, 1039, and 1029 known miRNAs were annotated in prostate tissues from the blank control (BC), normal saline injection (NS), and carrageenan injection (CAR) groups (n = 3 for each group), respectively. Among them, 84 miRNAs (CAR vs BC) and 70 miRNAs (CAR vs NS) with significantly different expression levels were identified. Compared with previously reported miRNAs with altered expression in various inflammatory diseases, the majority of deregulated miRNAs in CNP, such as miR-146b-5p, miR-155-5p, miR-150-5p, and miR-139-5p, showed similar expression patterns. Moreover, bioinformatics analyses have enriched mitogen-activated protein kinase (MAPK), cyclic adenosine monophosphate (cAMP), endocytosis, mammalian target of rapamycin (mTOR), and forkhead box O (FoxO) signaling pathways. These pathways were all involved in immune response, which indicates the critical regulatory role of the immune system in CNP initiation and progression. Our investigation has presented a global view of the differentially expressed miRNAs and potential regulatory networks containing their target genes, which may be helpful for identifying the novel mechanisms of miRNAs in immune regulation and effective target-specific theragnosis for CNP.

Prostate cancer (PCa) is a common lethal malignancy in men. RNA binding proteins (RBPs) have been proven to regulate the biological processes of various tumors, but their roles in PCa remain less defined. In the present study, we used bioinformatics analysis to identify RBP genes with prognostic and diagnostic values. A total of 59 differentially expressed RBPs in PCa were obtained, comprising 28 upregulated and 31 downregulated RBP genes, which may play important roles in PCa. Functional enrichment analyses showed that these RBPs were mainly involved in mRNA processing, RNA splicing, and regulation of RNA splicing. Additionally, we identified nine RBP genes (EXO1, PABPC1L, REXO2, MBNL2, MSI1, CTU1, MAEL, YBX2, and ESRP2) and their prognostic values by a protein-protein interaction network and Cox regression analyses. The expression of these nine RBPs was validated using immunohistochemical staining between the tumor and normal samples. Further, the associations between the expression of these nine RBPs and pathological T staging, Gleason score, and lymph node metastasis were evaluated. Moreover, these nine RBP genes showed good diagnostic values and could categorize the PCa patients into two clusters with different malignant phenotypes. Finally, we constructed a prognostic model based on these nine RBP genes and validated them using three external datasets. The model showed good efficiency in predicting patient survival and was independent of other clinical factors. Therefore, our model could be used as a supplement for clinical factors to predict patient prognosis and thereby improve patient survival.

The high-throughput sequencing technology, RNA-Seq, has been widely used to quantify gene and isoform expression in the study of transcriptome in recent years. Accurate expression measurement from the millions or billions of short generated reads is obstructed by difficulties. One is ambiguous mapping of reads to reference transcriptome caused by alternative splicing. This increases the uncertainty in estimating isoform expression. The other is non-uniformity of read distribution along the reference transcriptome due to positional, sequencing, mappability and other undiscovered sources of biases. This violates the uniform assumption of read distribution for many expression calculation approaches, such as the direct RPKM calculation and Poisson-based models. Many methods have been proposed to address these difficulties. Some approaches employ latent variable models to discover the underlying pattern of read sequencing. However, most of these methods make bias correction based on surrounding sequence contents and share the bias models by all genes. They therefore cannot estimate gene- and isoform-specific biases as revealed by recent studies.

Tumor metastasis is the predominant cause of lethality in cancer. We found that Neobractatin (NBT), a natural compound isolated from Garcinia bracteata, could efficiently inhibit breast and lung cancer cells metastasis. However, the mechanisms of NBT inhibiting cancer metastasis remain unclear. Based on the RNA-sequencing result and transcriptome analysis, Muscleblind-like 2 (MBNL2) was found to be significantly upregulated in the cells treated with NBT. The Cancer Genome Atlas (TCGA) database analysis indicated that the expression of MBNL2 in breast and lung carcinoma tumor tissues was significantly lower compared to normal tissues. We thus conducted to investigate the antimetastatic role of MBNL2. MBNL2 overexpression mimicked the effect of NBT on breast cancer and lung cancer cell motility and metastasis, in addition significantly enhanced the inhibition effect of NBT. MBNL2 knockdown furthermore partially eliminated the inhibitory effect of NBT on metastasis. Mechanistically, we demonstrated that NBT- and MBNL2-mediated antimetastasis regulation significantly correlated with the pAKT/epithelial-mesenchymal transition (EMT) pathway. Subsequent in vivo study showed the same metastasis inhibition effect in NBT and MBNL2 in MDA-MB-231 xenografts mouse model. This study suggest that NBT possesses significant antitumor activity in breast and lung cancer cells that is partly mediated through the MBNL2 expression and enhancement in metastasis via the pAKT/EMT signaling pathway.

Glaucoma and traumatic optic nerve crush (ONC) injury result in progressive loss of retinal ganglion cells (RGCs) and defects in visual function. In clinical trials of Traditional Chinese Medicine, acupuncture has been widely used for the treatment of ocular diseases. However, the molecular mechanisms of acupuncture treatment are still unclear. In this study, we used technique of RNA sequencing (RNA-seq) to study the effects of acupuncture treatment on retinal transcriptome after axotomy injury. RNA-seq results revealed that 436 genes including 31 transcription factors (TFs) were changed after injury, among them were many well-known neural degeneration related TFs such as Jun, Ddit3, Atf3, and Atf4. Interestingly, acupuncture treatment at acupoint GB20 (Fengchi) significantly reversed a series of differential expressed genes (DEGs) induced by optic nerve injury. While treatments at BL1 (Jingming) or GB20 sham control acupoint-GV16 (Fengfu), led to limited DEG reversal. In contrast, treatments at these two sites further enhanced the trend of DEG expression induced by axotomy injury. At last, retina immunostaining results revealed that only GB20 acupoint treatment increased RGC survival, in consistent with RNA-seq results. Therefore, our study first reported that acupuncture treatment regulated retinal transcriptome and reversed the gene expression induced by axotomy injury, and GB20 acupoint treatment increased RGC survival, which will provide novel therapeutic targets for treatment of ocular diseases.

Long non-coding RNA (lncRNA) has been associated with human diseases. To study the role of lncRNA in the pathogenic mechanism of acute fulminant myocarditis (AFM), we used a microarray to analyze lncRNA and messenger RNA (mRNA) expression in leukocyte samples from AFM patients and normal children. In total, using a 2/0.5-fold change and P < 0.05 as the cutoff criteria, we found that 3,101 lncRNAs and 2,170 mRNAs were differentially expressed in AFM patients. Quantitative real-time polymerase chain reaction (RT-qPCR) analysis was used to verify the microarray data. Eight differentially expressed molecules were randomly selected, including 3 upregulated lncRNAs, 3 downregulated lncRNAs, and 2 upregulated mRNAs. Among them, 7 expression profiles were consistent with the microarray results. Gene Ontology enrichment and Kyoto encyclopedia of genes and genomes (KEGG) pathway analysis were used to investigate the biological functions of these genes. Establishment of a lncRNA-mRNA co-expression network and lncRNA target predication were performed to study the molecular interactions of these molecules. Our study is the first to use microarrays to examine the lncRNA and mRNA expression profiles associated with AFM, and the results indicate that the immune system plays an important role in AFM. These findings may provide a new perspective for the pathogenesis, diagnosis, and therapy of AFM.