The conservation and signatures of lincRNAs in Marek's disease of chicken

Novel lncRNA 803 related to Marek's disease inhibits apoptosis of DF-1 cells

Marek's disease (MD) is a neoplastic disease that significantly affects the poultry industry. Long non-coding RNAs (lncRNAs) are crucial regulatory factors in various biological processes, including tumourigenesis. However, the involvement of novel lncRNAs in the course of MD virus (MDV) infection is still underexplored. Here, we present the first comprehensive characterization of differentially expressed lncRNAs in chicken spleen at different stages of MDV infection. A series of differentially expressed lncRNAs was identified at each stage of MDV infection through screening. Notably, our investigation revealed a novel lncRNA, lncRNA 803, which exhibited significant differential expression at different stages of MDV infection and was likely to be associated with the p53 pathway. Further analyses demonstrated that the overexpression of lncRNA 803 positively regulated the expression of p53 and TP53BP1 in DF-1 cells, leading to the inhibition of apoptosis. This is the first study to focus on the lncRNA expression profiles in chicken spleens during MDV pathogenesis. Our findings highlight the potential role of the p53-related novel lncRNA 803 in MD pathogenesis and provide valuable insights for decoding the molecular mechanism of MD pathogenesis involving non-coding RNA.RESEARCH HIGHLIGHTS Differentially expressed lncRNAs in spleens of chickens infected with Marek's disease virus at different stages were identified for the first time.The effects of novel lncRNA 803 on p53 pathway and apoptosis of DF-1 cells were reported for the first time.

The Identification and Function of Linc01615 on Influenza Virus Infection and Antiviral Response

Influenza virus infection poses a great threat to human health globally each year. Non-coding RNAs (ncRNAs) in the human genome have been reported to participate in the replication process of the influenza virus, among which there are still many unknowns about Long Intergenic Non-Coding RNAs (LincRNAs) in the cell cycle of viral infections. Here, we observed an increased expression of Linc01615 in A549 cells upon influenza virus PR8 infection, accompanied by the successful activation of the intracellular immune system. The knockdown of Linc01615 using the shRNAs promoted the proliferation of the influenza A virus, and the intracellular immune system was inhibited, in which the expressions of IFN-β, IL-28A, IL-29, ISG-15, MX1, and MX2 were decreased. Predictions from the catRAPID website suggested a potential interaction between Linc01615 and DHX9. Also, knocking down Linc01615 promoted influenza virus proliferation. The subsequent transcriptome sequencing results indicated a decrease in Linc01615 expression after influenza virus infection when DHX9 was knocked down. Further analysis through cross-linking immunoprecipitation and high-throughput sequencing (CLIP-seq) in HEK293 cells stably expressing DHX9 confirmed the interaction between DHX9 and Linc01615. We speculate that DHX9 may interact with Linc01615 to partake in influenza virus replication and that Linc01615 helps to activate the intracellular immune system. These findings suggest a deeper connection between DHX9 and Linc01615, which highlights the significant role of Linc01615 in the influenza virus replication process. This research provides valuable insights into understanding influenza virus replication and offers new targets for preventing influenza virus infections.

Transgenerational epigenetic inheritance and immunity in chickens that vary in Marek's disease resistance

Marek's disease virus (MDV), a naturally oncogenic, highly contagious alpha herpesvirus, induces a T cell lymphoma in chickens that causes severe economic loss. Marek's disease (MD) outcome in an individual is attributed to genetic and environmental factors. Further investigation of the host-virus interaction mechanisms that impact MD resistance is needed to achieve greater MD control. This study analyzed genome-wide DNA methylation patterns in 2 highly inbred parental lines 63 and 72 and 5 recombinant congenic strains (RCS) C, L, M, N, and X strains from those parents. Lines 63 and 72, are MD resistant and susceptible, respectively, whereas the RCS have different combinations of 87.5% Line 63 and 12.5% Line 72. Our DNA methylation cluster showed a strong association with MD incidence. Differentially methylated regions (DMRs) between the parental lines and the 5 RCS were captured. MD-resistant and MD-susceptible markers of DNA methylation were identified as transgenerational epigenetic inheritable. In addition, the growth of v-src DNA tumors and antibody response against sheep red blood cells differed among the 2 parental lines and the RCS. Overall, our results provide very solid evidence that DNA methylation patterns are transgenerational epigenetic inheritance (TEI) in chickens and also play a vital role in MD tumorigenesis and other immune responses; the specific methylated regions may be important modulators of general immunity.

Investigation of lncRNA in Bos taurus Mammary Tissue during Dry and Lactation Periods

This study aims to collect RNA-Seq data from Bos taurus samples representing dry and lactating mammary tissue, identify lncRNA transcripts, and analyze findings for their features and functional annotation. This allows for connections to be drawn between lncRNA and the lactation process. RNA-Seq data from 103 samples of Bos taurus mammary tissue were gathered from publicly available databases (60 dry, 43 lactating). The samples were filtered to reveal 214 dry mammary lncRNA transcripts and 517 lactating mammary lncRNA transcripts. The lncRNAs met common lncRNA characteristics such as shorter length, fewer exons, lower expression levels, and less sequence conservation when compared to the genome. Interestingly, several lncRNAs showed sequence similarity to genes associated with strong hair keratin intermediate filaments. Human breast cancer research has associated strong hair keratin filaments with mammary tissue cellular resilience. The lncRNAs were also associated with several genes/proteins that linked to pregnancy using expression correlation and gene ontology. Such findings indicate that there are crucial relationships between the lncRNAs found in mammary tissue and the development of the tissue, to meet both the animal's needs and our own production needs; these relationships should be further investigated to ensure that we continue to breed the most resilient, efficient dairy cattle.

Genome-wide analysis of long noncoding RNA profiles in pseudorabies-virus-infected PK15 cells

Recently, an increasing number of studies have shown that long noncoding RNAs (lncRNAs) are involved in host metabolism after infection with pseudorabies virus (PRV). In our study, via RNA sequencing analysis, a total of 418 mRNAs, 137 annotated lncRNAs, and 312 new lncRNAs were found to be differentially expressed. These lncRNAs were closely associated with metabolic regulation and immunity-related signalling pathways, including the T-cell receptor signalling pathway, chemokine signalling pathway, mitogen-activated protein kinase (MAPK) signalling pathway, TNF signalling pathway, Ras signalling pathway, calcium signalling pathway, and phosphatidylinositol signalling system. Real-time PCR indicated that several mRNAs and lncRNAs involved in the regulation of the immune effector process, T-cell receptor signalling pathway, TNF signalling pathway, MAPK signalling pathway, and chemokine signalling pathways were significantly expressed. These mRNAs and lncRNAs might play a role in PRV infection.

High-density single nucleotide polymorphism markers reveal the population structure of 2 local chicken genetic resources

Italy counts a large number of local chicken populations, some without a recognized genetic structure, such as Val Platani (VPL) and Cornuta (COS), which represent noteworthy local genetic resources. In this study, the genotype data of 34 COS and 42 VPL, obtained with the Affymetrix Axiom600KChicken Genotyping Array, were used with the aim to investigate the genetic diversity, the runs of homozygosity (ROH) pattern, as well as the population structure and relationship within the framework of other local Italian and commercial chickens. The genetic diversity indices, estimated using different approaches, displayed moderate levels of genetic diversity in both populations. The identified ROH hotspots harbored genes related to immune response and adaptation to local hot temperatures. The results on genetic relationship and population structure reported a clear clustering of the populations according to their geographic origin. The COS formed a nonoverlapping genomic cluster and clearly separated from the other populations, but showed evident proximity to the Siciliana breed (SIC). The VPL highlighted intermediate relationships between the COS-SIC group and the rest of the sample, but closer to the other Italian local chickens. Moreover, VPL showed a complex genomic structure, highlighting the presence of 2 subpopulations that match with the different source of the samples. The results obtained from the survey on genetic differentiation underline the hypothesis that Cornuta is a population with a defined genetic structure. The substructure that characterizes the Val Platani chicken is probably the consequence of the combined effects of genetic drift, small population size, reproductive isolation, and inbreeding. These findings contribute to the understanding of genetic diversity and population structure, and represent a starting point for designing programs to monitor and safeguard these local genetic resources, in order to define a possible official recognition program as breeds.

Critical roles of non-coding RNAs in lifecycle and biology of Marek's disease herpesvirus

Over the past two decades, numerous non-coding RNAs (ncRNAs) have been identified in different biological systems including virology, especially in large DNA viruses such as herpesviruses. As a representative oncogenic alphaherpesvirus, Marek's disease virus (MDV) causes an important immunosuppressive and rapid-onset neoplastic disease of poultry, namely Marek's disease (MD). Vaccinations can efficiently prevent the onset of MD lymphomas and other clinical disease, often heralded as the first successful example of vaccination-based control of cancer. MDV infection is also an excellent model for research into virally-induced tumorigenesis. Recently, great progress has been made in understanding the functions of ncRNAs in MD biology. Herein, we give a review of the discovery and identification of MDV-encoded viral miRNAs, focusing on the genomics, expression profiles, and emerging critical roles of MDV-1 miRNAs as oncogenic miRNAs (oncomiRs) or tumor suppressor genes involved in the induction of MD lymphomas. We also described the involvements of host cellular miRNAs, lincRNAs, and circRNAs participating in MDV life cycle, pathogenesis, and/or tumorigenesis. The prospects, strategies, and new techniques such as the CRISPR/Cas9-based gene editing applicable for further investigation into the ncRNA-mediated regulatory mechanisms in MDV pathogenesis/oncogenesis were also discussed, together with the possibilities of future studies on antiviral therapy and the development of new efficient MD vaccines.

Genome-Wide Analysis of Long Noncoding RNA Profiles in Seneca Valley Virus–Infected PK15 Cells

Long noncoding RNAs (lncRNAs) have been demonstrated to play key roles in various biological processes. However, the contributions of lncRNAs to Seneca Valley virus (SVV) infection and host defense remain largely unknown. In this study, differentially expressed lncRNAs and mRNAs in SVV-infected PK15 cells were detected by genome-wide analysis. A total of 14,127 lncRNAs and 63,562 mRNAs were identified, and 1,780 lncRNAs were differentially expressed. The functional prediction of SVV-induced lncRNAs showed high associations with biological regulation and many immunity-related signaling pathways, including the B-cell receptor pathway, RIG-I-like receptor signaling pathway, and NF-kappa B (NF-κB) signaling pathway. We next screened lncRNAs and target genes related to immune response pathways and further demonstrated their differential expression in SVV-infected PK15 cells. Our study investigated the function of lncRNAs involved in SVV infection and provided new insight into the pathogenic mechanisms of SVV.

Whole transcriptome analysis of RNA expression profiles reveals the potential regulating action of long noncoding RNA in lactating cows fed a high concentrate diet

Subacute ruminal acidosis (SARA) is a common metabolic disease in the dairy farming industry which is usually caused by an excessive amount of high concentrate diet. SARA not only threatens animal welfare but also leads to economic losses in the farming industry. The liver plays an important role in the distribution of nutritional substances and metabolism; however, a high concentrate diet can cause hepatic metabolic disorders and liver injury. Recently, noncoding RNA has been considered as a critical regulator of hepatic disease, however, its role in the bovine liver is limited. In this study, 12 mid-lactating dairy cows were randomly assigned to a control (CON) group (40% concentrate of dry matter, n = 6) and a SARA group (60% concentrate of dry matter, n = 6). After 21 d of treatment, all cows were sacrificed, and liver tissue samples were collected. Three dairy cows were randomly selected from the CON and SARA groups respectively to perform whole transcriptome analysis. More than 20,000 messenger RNA (mRNA), 10,000 long noncoding RNA (lncRNA), 3,500 circular RNA (circRNA) and 1,000 micro RNA (miRNA) were identified. Furthermore, 43 mRNA, 121 lncRNA and 3 miRNA were differentially expressed, whereas no obvious differentially expressed circRNA were detected between the 2 groups. Gene Ontology (GO) annotation revealed that the differentially expressed genes were mainly enriched in oxidoreductase activity, stress, metabolism, the immune response, cell apoptosis, and cell proliferation. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis showed that the deferentially expressed genes were highly enriched in the phosphatidylinositol 3 kinase (PI3K)-serine/threonine kinase (AKT) signaling pathway (P < 0.05). According to KEGG pathway analysis, the differentially expressed lncRNA (DElncRNA) target genes were mainly related to proteasomes, peroxisomes, and the hypoxia-inducible factor-1 signaling pathway (P < 0.005). Further bioinformatics and integrative analyses revealed that the lncRNA were strongly correlated with mRNA; therefore, it is reasonable to speculate that lncRNA potentially play important roles in the liver dysfunction induced by SARA. Our study provides a valuable resource for future investigations on the mechanisms of SARA to facilitate an understanding of the importance of lncRNA, and offer functional RNA information.

LncRNAs in domesticated animals: from dog to livestock species

Animal genomes are pervasively transcribed into multiple RNA molecules, of which many will not be translated into proteins. One major component of this transcribed non-coding genome is the long non-coding RNAs (lncRNAs), which are defined as transcripts longer than 200 nucleotides with low coding-potential capabilities. Domestic animals constitute a unique resource for studying the genetic and epigenetic basis of phenotypic variations involving protein-coding and non-coding RNAs, such as lncRNAs. This review presents the current knowledge regarding transcriptome-based catalogues of lncRNAs in major domesticated animals (pets and livestock species), covering a broad phylogenetic scale (from dogs to chicken), and in comparison with human and mouse lncRNA catalogues. Furthermore, we describe different methods to extract known or discover novel lncRNAs and explore comparative genomics approaches to strengthen the annotation of lncRNAs. We then detail different strategies contributing to a better understanding of lncRNA functions, from genetic studies such as GWAS to molecular biology experiments and give some case examples in domestic animals. Finally, we discuss the limitations of current lncRNA annotations and suggest research directions to improve them and their functional characterisation.

Insights Into mRNA and Long Non-coding RNA Profiling RNA Sequencing in Uterus of Chickens With Pink and Blue Eggshell Colors

The blue egg is both of biological interest and economic importance for consumers, egg retailers, and scientists. To date, the genetic mechanisms underlying pigment have mainly focused on protein-coding genes. However, the underpinning mechanism of non-coding RNAs on the pigment deposition among different eggshell colors remains unknown. In this study, RNA sequencing was employed to profile the uterine gland transcriptome (lncRNA and mRNA) of 15 Changshun blue eggshell layers, to better understand the genetic mechanisms of deposition of blue eggshell color. Results showed that differentially expressed mRNAs, GO terms, and KEGG pathways among pink-eggshell and blue-eggshell chickens were mainly targeting immune- and transporter-related terms with the SLC family, IgJ, CD family, and MTMR genes. Furthermore, the progesterone-mediated oocyte maturation and cortisol synthesis and secretion pathway with targeted gene PGR and Pbx1 were significantly enriched between blue- and pink-eggshell chickens. Integrating analysis of lncRNA and mRNA profiles predicted 4 and 25 lncRNA–gene pairs by antisense and cis analysis. They were relative to immune, nerve, and lipids and amino acid metabolisms, porphyrin, and chlorophyll metabolism with targeted gene FECH and oxidative phosphorylation and cardiac muscle contraction pathways with targeted gene COX6A1. Within blue-eggshell chickens, the GO terms hindbrain tangential cell migration and phosphatidylinositol monophosphate phosphatase activity with targeted gene Plxna2 and MTRM1 were identified. Integrating analysis of lncRNA and mRNA profiles predicted 8 and 22 lncRNA–gene pairs. Most pathways were mainly enriched on lipid-related metabolisms as found in mRNA sequencing. The lncRNAs did exert similar functions in color formation by modulating pigment disposition and immune- and lipid-related metabolisms. Our results provide a catalog of chicken uterine lncRNAs and genes worthy of further studies to understand their roles in the selection for blue eggshell color layers.

Expression profile of long non-coding RNAs in porcine lymphnode response to porcine circovirus type 2 infection

Porcine circovirus type 2 (PCV2) can cause various clinical diseases in pigs, resulting in huge losses for the pig farms all over the world. In order to develop a new strategy to control PCV2, it is essential to understand its mechanisms firstly, especially PCV2 interferes with the host's innate immunity. In the present study, lncRNA and mRNA expression profiles in porcine lymphnode response to PCV2 infection were deeply sequenced and analyzed. 3271 novel lncRNAs were identified in all. 1898 mRNAs and 282 lncRNAs showed differential expression between control and PCV2-infected groups. The bioinformatics analysis including lncRNA-mRNA co-expression network construction, as well as GO and KEGG pathway analysis focused on the DEGs was carried out. The results indicated that lncRNAs might participate in PCV2 infection-induced the pathogenesis of immunosuppression through regulating the host's immune responses, biological regulation, response to stimulus, cellular component organization or biogenesis and metabolism. And these differentially expressed lncRNAs might play important roles in response to PCV2 infection in the host's innate immune system. These findings provided a large-scale survey of dysregulated lncRNAs after PCV2 infection, especially the lncRNAs responded to host's innate immune within the lymphnode. This study will provide a novel insight into the lncRNAs' functions and the possible immunosuppressive mechanism induced by PCV2 infection. However, further research will be required to verify the characteristic function of the dysregulated lncRNAs.

Transcriptome-wide N6-methyladenosine modification profiling of long non-coding RNAs during replication of Marek's disease virus in vitro

Background

The newly discovered reversible N6-methyladenosine (m⁶A) modification plays an important regulatory role in gene expression. Long non-coding RNAs (lncRNAs) participate in Marek’s disease virus (MDV) replication but how m⁶A modifications in lncRNAs are affected during MDV infection is currently unknown. Herein, we profiled the transcriptome-wide m⁶A modification in lncRNAs in MDV-infected chicken embryo fibroblast (CEF) cells.

Results

Methylated RNA immunoprecipitation sequencing results revealed that the lncRNA m⁶A modification is highly conserved with MDV infection increasing the expression of lncRNA m⁶A modified sites compared to uninfected cell controls. Gene Ontology and the Kyoto Encyclopedia of Genes and Genomes pathway analysis revealed that lncRNA m⁶A modifications were highly associated with signaling pathways associated with MDV infection.

Conclusions

In this study, the alterations seen in transcriptome-wide m⁶A occurring in lncRNAs following MDV-infection suggest this process plays important regulatory roles during MDV replication. We report for the first time profiling of the alterations in transcriptome-wide m⁶A modification in lncRNAs of MDV-infected CEF cells.

Role of microRNA and long non-coding RNA in Marek's disease tumorigenesis in chicken

Marek's disease virus (MDV), the causative agent of Marek's disease (MD), results in highly infectious phymatosis, lymphatic tissue hyperplasia, and neoplasia. MD is associated with high morbidity and mortality rate. Non-coding RNAs (ncRNAs) entails long non-coding RNA (lncRNA) and microRNA (miRNA). Numerous studies have reported that specific miRNAs and lncRNAs participate in multiple cellular processes, such as proliferation, migration, and tumor cell invasion. Specialized miRNAs and lncRNAs militate a similar role in MD tumor oncogenesis. Despite its growing popularity, only a few reviews are available on ncRNA in MDV tumor oncogenes. Herein, we summarized the role of the miRNAs and lncRNAs in MD tumorigenesis. Altogether, we brought forth the research issues, such as MD prevention, screening, regulatory network formation, novel miRNAs, and lncRNAs analysis in MD that needs to be explored further. This review provides a theoretical platform for the further analysis of miRNAs and lncRNAs functions and the prevention and control of MD and malignancies in domestic animals.

Epigenetic Regulation by Non-Coding RNAs in the Avian Immune System

The identified non-coding RNAs (ncRNAs) include circular RNAs, long non-coding RNAs, microRNAs, ribosomal RNAs, small interfering RNAs, small nuclear RNAs, piwi-interacting RNAs, and transfer RNAs, etc. Among them, long non-coding RNAs, circular RNAs, and microRNAs are regulatory RNAs that have different functional mechanisms and were extensively participated in various biological processes. Numerous research studies have found that circular RNAs, long non-coding RNAs, and microRNAs played their important roles in avian immune system during the infection of parasites, virus, or bacterium. Here, we specifically review and expand this knowledge with current advances of circular RNAs, long non-coding RNAs, and microRNAs in the regulation of different avian diseases and discuss their functional mechanisms in response to avian diseases.

Circular RNA expression profiling reveals that circ-PLXNA1 functions in duck adipocyte differentiation

Adipocytes are derived from pluripotent mesenchymal stem cells through adipogenesis. Pre-adipocyte differentiation in poultry greatly influences fat deposition and meat quality. Circular RNAs (circRNAs) have an important function in cancer and some differentiation processes. Herein, high-throughput transcriptome sequencing was used to detect circRNAs present in cherry valley duck pre-adipocyte and adipocyte differentiation over 3 days. We identified 9,311 circRNAs and 141 differentially expressed circRNAs. Sequencing results were verified through qRT-PCR using seven randomly selected circRNAs, and competing endogenous RNA (ceRNA) networks were exhibited by ten important circRNAs in duck adipocyte differentiation. circRNA plexin A1 (circ-PLXNA1) was detected in duck adipocytes and mainly expressed in adipose, leg muscle and liver. Inhibition of circ-PLXNA1 limited the differentiation of duck adipocyte. There were four corresponding miRNAs for circ-PLXNA1 and 313 target genes for those miRNAs. CeRNA“circ-PLXNA1/miR-214/CTNNB1 axis” was focused and verified by a dual-luciferase reporter experiment. After co-transfection of cells with si-circ-PLXNA1 and miR-214 mimics, the expression level of CTNNB1 was down-regulated, triglyceride content and the adipogenic capacity of preadipocytes decreased. While there were no significant change after si-CTNNB1 transfection. All these results provide further insight into the circRNAs, especially for circ-PLXNA1 in duck adipocyte differentiation.

Genome-wide analysis of long noncoding RNA profiles in Vero cells infected with porcine epidemic diarrhea virus

Porcine epidemic diarrhea (PED) is an acute enteric disease caused by porcine epidemic diarrhea virus (PEDV). In China, variant PEDV causes severe watery diarrhea, vomiting, and dehydration in piglets, leading to very high morbidity and mortality. However, the pathogenesis of PEDV is still not fully understood. In our study, we analyzed the long noncoding RNA (lncRNA) and mRNA expression profiles of PEDV GDgh16 in infected Vero cells at 60 h postinfection. A total of 61,790 annotated mRNAs, 14,247 annotated lncRNAs and 1290 novel lncRNAs were identified. A total of 227 annotated lncRNAs and 13 novel lncRNAs were significantly and differentially expressed after viral infection. The Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) databases were used to identify genes adjacent to the lncRNAs, and it was found that these lncRNAs were enriched in pathways related to immune and antiviral responses. Next, we selected candidate lncRNAs and their predicted target genes for study. RT-qPCR demonstrated that these lncRNAs and genes were differentially expressed after PEDV infection. Our study investigated the function of lncRNAs involved in PEDV infection, providing new insight into the pathogenic mechanisms of PEDV.

Long noncoding RNA and mRNA profiling of hypothalamic-pituitary-mammary gland axis in lactating sows under heat stress

Heat stress (HS) seriously affects sow lactation performance and Long non-coding RNAs (lncRNAs) play vital roles in the regulation of transcription and post transcription. However, the mechanism of lncRNAs expression affecting lactation performance on the hypothalamus-pituitary-mammary axis of sows is still unclear. In this study, we performed RNA sequencing and bioinformatics analysis of the hypothalamus, pituitary, and mammary gland tissues of lactating sows under HS and thermal comfort. In total, the analysis identified 658, 6021, and 6745 differently expressed (DE) mRNAs, 26, 126, and 169 DE lncRNAs between comparison groups in the hypothalamus, pituitary, and mammary glands, respectively. The hormone genes and most DE mRNAs encoding heat shock protein were differently expressed in the HS group. In addition, 2, 60, and 86 pairs of DE lncRNAs and mRNAs correlation were observed in those tissues, respectively. Some lncRNAs may be involved in the regulation of lactation performance in the HS sows.

Analysis of Long Non-Coding RNAs and mRNAs Associated with Lactation in the Crop of Pigeons (Columba livia)

Pigeons have the ability to produce milk and feed their squabs. The genetic mechanisms underlying milk production in the crops of 'lactating' pigeons are not fully understood. In this study, RNA sequencing was employed to profile the transcriptome of lncRNA and mRNA in lactating and non-'lactating' pigeon crops. We identified 7066 known and 17,085 novel lncRNAs. Of these lncRNAs, 6166 were differentially expressed. Among the 15,138 mRNAs detected, 6483 were differentially expressed, including many predominant genes with known functions in the milk production of mammals. A GO annotation analysis revealed that these genes were significantly enriched in 55, 65, and 30 pathways of biological processes, cellular components, and molecular functions, respectively. A KEGG pathway enrichment analysis revealed that 12 pathways (involving 544 genes), including the biosynthesis of amino acids, the propanoate metabolism, the carbon metabolism and the cell cycle, were significantly enriched. The results provide fundamental evidence for the better understanding of lncRNAs' and differentially expressed genes' (DEGs) regulatory role in the molecular pathways governing milk production in pigeon crops. To our knowledge, this is the first genome-wide investigation of the lncRNAs in pigeon crop associated with milk production. This study provided valuable resources for differentially expressed lncRNAs and mRNAs, improving our understanding of the molecular mechanism of pigeon milk production.

Linc-GALMD1 Regulates Viral Gene Expression in the Chicken

A rapidly increasing number of reports on dysregulated long intergenic non-coding RNA (lincRNA) expression across numerous types of cancers indicates that aberrant lincRNA expression may be a major contributor to tumorigenesis. Marek’s disease (MD) is a T cell lymphoma of chickens induced by Marek’s disease virus (MDV). Although we have investigated the roles of lincRNAs in bursa tissue of MDV-infected chickens in previous studies, the molecular mechanisms of lincRNA functions in T cells remain poorly understood. In the present study, Linc-GALMD1 was identified from CD4+ T cells and MSB1 cells, and its expression was significantly downregulated in MD-resistant line of birds in response to MDV challenge. Furthermore, loss-of-function experiments indicated that linc-GALMD1 significantly affected the expression of 290 genes in trans. Through integrated analysis of differentially expressed genes (DEGs) induced by MDV and linc-GALMD1, we found that IGLL1 gene expression levels had a positive correlation with the degree of MD infection and could potentially serve as an indicator for clinical diagnosis of MD. Moreover, an interaction between MDV and linc-GALMD1 was also observed. Accordingly, chicken embryonic fibroblast cells were inoculated with MDV with and without the linc-GALMD1 knockdown, and the data showed that linc-GALMD1 could repress MDV gene expression during the course of MDV infection. These findings uncovered a role of linc-GALMD1 as a viral gene regulator and suggested a function of linc-GALMD1 contributing to tumor suppression by coordinating expression of MDV genes and tumor-related genes and regulating immune responses to MDV infection.

Identification and analysis of long non-coding RNAs that are involved in inflammatory process in response to transmissible gastroenteritis virus infection

Background

Transmissible gastroenteritis virus (TGEV) infection can cause acute inflammation. Long noncoding RNAs (lncRNAs) play important roles in a number of biological process including inflammation response. However, whether lncRNAs participate in TGEV-induced inflammation in porcine intestinal epithelial cells (IPECs) is largely unknown.

Results

In this study, the next-generation sequencing (NGS) technology was used to analyze the profiles of lncRNAs in Mock and TGEV-infected porcine intestinal epithelial cell-jejunum 2 (IPEC-J2) cell line. A total of 106 lncRNAs were differentially expressed. Many differentially expressed lncRNAs act as elements to competitively attach microRNAs (miRNAs) which target to messenger RNA (mRNAs) to mediate expression of genes that related to toll-like receptors (TLRs), NOD-like receptors (NLRs), tumor necrosis factor (TNF), and RIG-I-like receptors (RLRs) pathways. Functional analysis of the binding proteins and the up/down-stream genes of the differentially expressed lncRNAs revealed that lncRNAs were principally related to inflammatory response. Meanwhile, we found that the differentially expressed lncRNA TCONS_00058367 might lead to a reduction of phosphorylation of transcription factor p65 (p-p65) in TGEV-infected IPEC-J2 cells by negatively regulating its antisense gene promyelocytic leukemia (PML).

Conclusions

The data showed that differentially expressed lncRNAs might be involved in inflammatory response induced by TGEV through acting as miRNA sponges, regulating their up/down-stream genes, or directly binding proteins.

Genome-Wide Identification of Long Non-Coding RNAs and Their Regulatory Networks Involved in Apis mellifera ligustica Response to Nosema ceranae Infection

Long non-coding RNAs (lncRNAs) are a diverse class of transcripts that structurally resemble mRNAs but do not encode proteins, and lncRNAs have been proven to play pivotal roles in a wide range of biological processes in animals and plants. However, knowledge of expression patterns and potential roles of honeybee lncRNA response to Nosema ceranae infection is completely unknown. Here, we performed whole transcriptome strand-specific RNA sequencing of normal midguts of Apis mellifera ligustica workers (Am7CK, Am10CK) and N. ceranae-inoculated midguts (Am7T, Am10T), followed by comprehensive analyses using bioinformatic and molecular approaches. A total of 6353 A. m. ligustica lncRNAs were identified, including 4749 conserved lncRNAs and 1604 novel lncRNAs. These lncRNAs had minimal sequence similarities with other known lncRNAs in other species; however, their structural features were similar to counterparts in mammals and plants, including shorter exon and intron length, lower exon number, and lower expression level, compared with protein-coding transcripts. Further, 111 and 146 N. ceranae-responsive lncRNAs were identified from midguts at 7-days post-inoculation (dpi) and 10 dpi compared with control midguts. Twelve differentially expressed lncRNAs (DElncRNAs) were shared by Am7CK vs. Am7T and Am10CK vs. Am10T comparison groups, while the numbers of unique DElncRNAs were 99 and 134, respectively. Functional annotation and pathway analysis showed that the DElncRNAs may regulate the expression of neighboring genes by acting in cis and trans fashion. Moreover, we discovered 27 lncRNAs harboring eight known miRNA precursors and 513 lncRNAs harboring 2257 novel miRNA precursors. Additionally, hundreds of DElncRNAs and their target miRNAs were found to form complex competitive endogenous RNA (ceRNA) networks, suggesting that these DElncRNAs may act as miRNA sponges. Furthermore, DElncRNA-miRNA-mRNA networks were constructed and investigated, the results demonstrated that a portion of the DElncRNAs were likely to participate in regulating the host material and energy metabolism as well as cellular and humoral immune host responses to N. ceranae invasion. Our findings revealed here offer not only a rich genetic resource for further investigation of the functional roles of lncRNAs involved in the A. m. ligustica response to N. ceranae infection, but also a novel insight into understanding the host-pathogen interaction during honeybee microsporidiosis.

Conserved regions in long non-coding RNAs contain abundant translation and protein-RNA interaction signatures

The mammalian transcriptome includes thousands of transcripts that do not correspond to annotated protein-coding genes and that are known as long non-coding RNAs (lncRNAs). A handful of lncRNAs have well-characterized regulatory functions but the biological significance of the majority of them is not well understood. LncRNAs that are conserved between mice and humans are likely to be enriched in functional sequences. Here, we investigate the presence of different types of ribosome profiling signatures in lncRNAs and how they relate to sequence conservation. We find that lncRNA-conserved regions contain three times more ORFs with translation evidence than non-conserved ones, and identify nine cases that display significant sequence constraints at the amino acid sequence level. The study also reveals that conserved regions in intergenic lncRNAs are significantly enriched in protein–RNA interaction signatures when compared to non-conserved ones; this includes sites in well-characterized lncRNAs, such as Cyrano, Malat1, Neat1 and Meg3, as well as in tens of lncRNAs of unknown function. This work illustrates how the analysis of ribosome profiling data coupled with evolutionary analysis provides new opportunities to explore the lncRNA functional landscape.

Marek's Disease Virus Infection Induced Mitochondria Changes in Chickens

Mitochondria are crucial cellular organelles in eukaryotes and participate in many cell processes including immune response, growth development, and tumorigenesis. Marek’s disease (MD), caused by an avian alpha-herpesvirus Marek’s disease virus (MDV), is characterized with lymphomas and immunosuppression. In this research, we hypothesize that mitochondria may play roles in response to MDV infection. To test it, mitochondrial DNA (mtDNA) abundance and gene expression in immune organs were examined in two well-defined and highly inbred lines of chickens, the MD-susceptible line 7₂ and the MD-resistant line 6₃. We found that mitochondrial DNA contents decreased significantly at the transformation phase in spleen of the MD-susceptible line 7₂ birds in contrast to the MD-resistant line 6₃. The mtDNA-genes and the nucleus-genes relevant to mtDNA maintenance and transcription, however, were significantly up-regulated. Interestingly, we found that POLG2 might play a potential role that led to the imbalance of mtDNA copy number and gene expression alteration. MDV infection induced imbalance of mitochondrial contents and gene expression, demonstrating the indispensability of mitochondria in virus-induced cell transformation and subsequent lymphoma formation, such as MD development in chicken. This is the first report on relationship between virus infection and mitochondria in chicken, which provides important insights into the understanding on pathogenesis and tumorigenesis due to viral infection.

Integrated analysis of lncRNA and mRNA repertoires in Marek's disease infected spleens identifies genes relevant to resistance

BACKGROUND

Marek's disease virus (MDV) is an oncogenic herpesvirus that can cause T-cell lymphomas in chicken. Long noncoding RNA (lncRNA) is strongly associated with various cancers and many other diseases. In chickens, lncRNAs have not been comprehensively identified. Here, we profiled mRNA and lncRNA repertoires in three groups of spleens from MDV-infected and non-infected chickens, including seven tumorous spleens (TS) from MDV-infected chickens, five spleens from the survivors (SS) without lesions after MDV infection, and five spleens from noninfected chickens (NS), to explore the underlying mechanism of host resistance in Marek's disease (MD).

RESULTS

By using a precise lncRNA identification pipeline, we identified 1315 putative lncRNAs and 1166 known lncRNAs in spleen tissue. Genomic features of putative lncRNAs were characterized. Differentially expressed (DE) mRNAs, putative lncRNAs, and known lncRNAs were profiled among three groups. We found that several specific intergroup differentially expressed genes were involved in important biological processes and pathways, including B cell activation and the Wnt signaling pathway; some of these genes were also found to be the hub genes in the co-expression network analyzed by WGCNA. Network analysis depicted both intergenic correlation and correlation between genes and MD traits. Five DE lncRNAs including MSTRG.360.1, MSTRG.6725.1, MSTRG.6754.1, MSTRG.15539.1, and MSTRG.7747.5 strongly correlated with MD-resistant candidate genes, such as IGF-I, CTLA4, HDAC9, SWAP70, CD72, JCHAIN, CXCL12, and CD8B, suggesting that lncRNAs may affect MD resistance and tumorigenesis in chicken spleens through their target genes.

CONCLUSIONS

Our results provide both transcriptomic and epigenetic insights on MD resistance and its pathological mechanism. The comprehensive lncRNA and mRNA transcriptomes in MDV-infected chicken spleens were profiled. Co-expression analysis identified integrated lncRNA-mRNA and gene-gene interaction networks, implying that hub genes or lncRNAs exert critical influence on MD resistance and tumorigenesis.

Expression patterns of novel circular RNAs in chicken cells after avian leukosis virus subgroup J infection

Avian leukosis virus subgroup J (ALV-J) is an oncogenic retrovirus that causes severe economic losses to the poultry industry worldwide. Circular RNAs (circRNAs) are a class of non-coding RNAs that has been described in various biological systems and pathogenic processes. However, the immune mechanisms in response to circRNAs remain unknown. In this study, high-throughput transcriptome sequencing was used to detect circRNAs present in chicken macrophage (HD11) and chick embryo fibroblast (CEF) cells infected with ALV-J. We identified 7684 circRNAs from diverse genomic locations in CEF and HD11 after ALV-J infection, these RNAs showed complex expression patterns that differed based on the cells type and infection time. In total, 302 differentially expressed (DE) circRNAs and 164 DE circRNAs were identified in CEF and HD11 after ALV-J infection, respectively. CircRNA7419-associated with KDM4C- and circRNA6679 and circRNA6680-associated with TNFAIP6- were involved in the immune response upon ALV-J infection in CEF. Host genes were analyzed through further bioinformatics analysis. The result confirmed that a large number of DE circRNAs corresponded to several immune-associated or tumor-associated terms and pathways, such as Mucin type O-Glycan biosynthesis, MAPK signaling pathway, B cell receptor signaling, and Wnt signaling pathway in CEF, as well as Jak-STAT signaling pathway, apoptosis, and MAPK signaling pathway in HD11. CircRNAs related to the B cell receptor signaling pathway in CEF, and the Jak-STAT signaling pathway in HD11, were selected for circRNA-miRNA interaction network analyses. Our study indicates that circRNAs expression was altered by ALV-J infection in both CEF and HD11, and may play a key role in the progression of ALV-J infection.

Porcine endemic diarrhea virus infection regulates long noncoding RNA expression

Long noncoding RNAs (lncRNAs) have been implicated in various life processes. However, the lncRNA expression and potential functions in porcine endemic diarrhea virus (PEDV) infection and host defense are still poorly understood. In this study, we investigated the lncRNA expression profiles during PEDV infection in intestinal porcine epithelial cell-jejunum 2 (IPEC-J2) cell lines by next-generation sequencing and identified 6188 novel lncRNAs. The functional annotation analysis revealed that these lncRNAs might be associated with many immunity-related genes. We next selected candidate lncRNAs related to immune response pathways and further identified their differential expression in PEDV-infected IPEC-J2 cells and newborn piglets. Our results demonstrated that PEDV infection regulated lncRNA expression patterns in both the IPEC-J2 cell line and piglet ileum. These findings provide the first large-scale survey of lncRNAs associated with PEDV infection, specifically the lncRNAs responsible for the activation of the immune system within the ileum.

Analysis of expression profiles of long noncoding RNAs and mRNAs in brains of mice infected by rabies virus by RNA sequencing

Rabies, caused by rabies virus (RABV), is still the deadliest infectious disease. Mechanism of host immune response upon RABV infection is not yet fully understood. Accumulating evidences suggest that long noncoding RNAs (lncRNAs) plays key roles in host antiviral responses. However, expression profile and function of lncRNAs in RABV infection remain unclear. In the present study, expression profile of lncRNAs and mRNAs profiles were investigated in RABV-infected brain tissues of mice by RNA sequencing. A total of 140 lncRNAs and 3,807 mRNAs were differentially expressed in RABV-infected animals. The functional annotation and enrichment analysis using Gene Oncology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) revealed that differentially expressed transcripts were predominantly involved in signaling pathways related to host immune response. The expression profiles of the selected lncRNAs in brains of mice during RABV infections were verified by quantitative real time polymerase chain reaction (qRT-PCR). To our knowledge, this is the first report to profile the lncRNA expression in RABV infected mice. Our findings provide insights into understanding the role of lncRNAs in host immune response against RABV infection.

Genome Wide Identification of Novel Long Non-coding RNAs and Their Potential Associations With Milk Proteins in Chinese Holstein Cows

Long non-coding RNAs (lncRNAs) have emerged as a novel class of regulatory molecules involved in various biological processes. However, their role in milk performance is unknown. Here, whole transcriptome RNA sequencing was used to generate the lncRNA transcriptome profiles in mammary tissue samples from 6 Chinese Holstein cows with 3 extremely high and 3 low milk protein percentage phenotypes. In this study, 6,450 lncRNA transcripts were identified through 5 stringent steps and filtration by coding potential. In total, 31 lncRNAs and 18 novel genes were identified to be differentially expressed in high milk protein samples (HP) relative to low milk protein samples (LP), respectively. Differentially expressed lncRNAs were selected to predict target genes through bioinformatics analysis, followed by the integration of differentially expressed mRNA data, gene function, gene ontology (GO) and pathway, genome wide association study (GWAS) and quantitative trait locus (QTL) information, as well as network analysis to further characterize potential interactions. Several lncRNAs were found (such as XLOC_059976) that could be used as candidate markers for milk protein content prediction. This is the first study to perform global expression profiling of lncRNAs and mRNAs related to milk protein traits in dairy cows. These results provide important information and insights into the synthesis of milk proteins, and potential targets for the future improvement of milk quality.

Molecular mechanisms of long noncoding RNAs and their role in disease pathogenesis

LncRNAs are long non-coding regulatory RNAs that are longer than 200 nucleotides. One of the major functions of lncRNAs is the regulation of specific gene expression at multiple steps including, recruitment and expression of basal transcription machinery, post-transcriptional modifications and epigenetics [1]. Emerging evidence suggests that lncRNAs also play a critical role in maintaining tissue homeostasis during physiological and pathological conditions, lipid homeostasis, as well as epithelial and smooth muscle cell homeostasis, a topic that has been elegantly reviewed [2-5]. While aberrant expression of lncRNAs has been implicated in several disease conditions, there is paucity of information about their contribution to the etiology of diseases [6]. Several studies have compared the expression of lncRNAs under normal and cancerous conditions and found differential expression of several lncRNAs, suggesting thereby an involvement of lncRNAs in disease processes [7, 8]. Furthermore, the ability of lncRNAs to influence epigenetic changes also underlies their role in disease pathogenesis since epigenetic regulation is known to play a critical role in many human diseases [1]. LncRNAs thus are not only involved in homeostatic functioning but also play a vital role in the progression of many diseases, thereby underscoring their potential as novel therapeutic targets for the alleviation of a variety of human disease conditions.

Dynamic Changes in the Splenic Transcriptome of Chickens during the Early Infection and Progress of Marek's Disease

Gallid alphaherpesvirus 2 (GaHV2) is an oncogenic avian herpesvirus inducing Marek’s disease (MD) and rapid-onset T-cell lymphomas. To reveal molecular events in MD pathogenesis and tumorigenesis, the dynamic splenic transcriptome of GaHV2-infected chickens during early infection and pathogenic phases has been determined utilizing RNA-seq. Based on the significant differentially expressed genes (DEGs), analysis of gene ontology, KEGG pathway and protein-protein interaction network has demonstrated that the molecular events happening during GaHV2 infection are highly relevant to the disease course. In the ‘Cornell Model’ description of MD, innate immune responses and inflammatory responses were established at early cytolytic phase but persisted until lymphoma formation. Humoral immunity in contrast began to play a role firstly in the intestinal system and started at late cytolytic phase. Neurological damage caused by GaHV2 is first seen in early cytolytic phase and is then sustained throughout the following phases over a long time period. During the proliferative phase many pathways associated with transcription and/or translation were significantly enriched, reflecting the cell transformation and lymphoma formation. Our work provides an overall view of host responses to GaHV2 infection and offers a meaningful basis for further studies of MD biology.

Long intergenic non-coding RNA GALMD3 in chicken Marek's disease

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.

A subset of conserved mammalian long non-coding RNAs are fossils of ancestral protein-coding genes

Background

Only a small portion of human long non-coding RNAs (lncRNAs) appear to be conserved outside of mammals, but the events underlying the birth of new lncRNAs in mammals remain largely unknown. One potential source is remnants of protein-coding genes that transitioned into lncRNAs.

Results

We systematically compare lncRNA and protein-coding loci across vertebrates, and estimate that up to 5% of conserved mammalian lncRNAs are derived from lost protein-coding genes. These lncRNAs have specific characteristics, such as broader expression domains, that set them apart from other lncRNAs. Fourteen lncRNAs have sequence similarity with the loci of the contemporary homologs of the lost protein-coding genes. We propose that selection acting on enhancer sequences is mostly responsible for retention of these regions. As an example of an RNA element from a protein-coding ancestor that was retained in the lncRNA, we describe in detail a short translated ORF in the JPX lncRNA that was derived from an upstream ORF in a protein-coding gene and retains some of its functionality.

Conclusions

We estimate that ~ 55 annotated conserved human lncRNAs are derived from parts of ancestral protein-coding genes, and loss of coding potential is thus a non-negligible source of new lncRNAs. Some lncRNAs inherited regulatory elements influencing transcription and translation from their protein-coding ancestors and those elements can influence the expression breadth and functionality of these lncRNAs.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-017-1293-0) contains supplementary material, which is available to authorized users.

Analyses of Long Non-Coding RNA and mRNA profiling using RNA sequencing in chicken testis with extreme sperm motility

Sperm motility is the most important indicator in evaluating roosters' fecundity. However, the genetic mechanisms underlying chicken sperm motility is not yet clear. Long non-coding RNA (lncRNA) play epigenetic roles in reproduction. In this study, RNA sequencing was employed to profile the testis transcriptome (lncRNA and mRNA) of six Beijing-you cocks divergent in sperm motility. In total, 2,597 lncRNAs were identified in the chicken testis, including 1,267 lincRNAs, 975 anti-sense lncRNAs, and 355 intronic lncRNAs. They shared similar features with previous studies. Of these lncRNAs, 124 were differentially expressed. Among 17,690 mRNAs detected in this study, 544 were differentially expressed, including a bunch of genes with known functions on sperm motility. GO annotation analysis revealed these genes were involved in ATP binding, cilium assembly, and oxidation-reduction processes. Integrating analysis of lncRNA and mRNA profiles predicted 10 lncRNA-gene pairs, including 8 co-regulated and 2 inversely-regulated pairs. To the best of our knowledge, this is the first genome-wide investigation of the lncRNAs in the chicken testis associated with sperm motility. Our results provided a catalog of chicken testis lncRNAs and genes worthy of further studies to understand their roles in cocks' reproductive performance regulation.

Annotating long intergenic non-coding RNAs under artificial selection during chicken domestication

Background

Numerous biological functions of long intergenic non-coding RNAs (lincRNAs) have been identified. However, the contribution of lincRNAs to the domestication process has remained elusive. Following domestication from their wild ancestors, animals display substantial changes in many phenotypic traits. Therefore, it is possible that diverse molecular drivers play important roles in this process.

Results

We analyzed 821 transcriptomes in this study and annotated 4754 lincRNA genes in the chicken genome. Our population genomic analysis indicates that 419 lincRNAs potentially evolved during artificial selection related to the domestication of chicken, while a comparative transcriptomic analysis identified 68 lincRNAs that were differentially expressed under different conditions. We also found 47 lincRNAs linked to special phenotypes.

Conclusions

Our study provides a comprehensive view of the genome-wide landscape of lincRNAs in chicken. This will promote a better understanding of the roles of lincRNAs in domestication, and the genetic mechanisms associated with the artificial selection of domestic animals.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-1036-6) contains supplementary material, which is available to authorized users.

Downregulated long noncoding RNA ALDBGALG0000005049 induces inflammation in chicken muscle suffered from selenium deficiency by regulating stearoyl-CoA desaturase

Long non-coding RNAs (lncRNAs) have been demonstrated to play a pivotal role in proliferation and differentiation of muscles. However, the study on the roles of lncRNAs in Selenium (Se) deficiency induced muscle injury is still unclear. In this study, deep sequencing was performed to profile lncRNAs and mRNAs of the muscles from the Se deficiency (-Se group) and control (C group) chickens. The results revealed that 38 lncRNAs (23 up-regulated and 15 down-regulated) and 687 mRNAs (285 up-regulated and 402 down-regulated) were significantly dysregulated expressed, and the significantly dysregulated pathway including Phagosome, Cardiac muscle contraction, and Peroxisome Proliferator-Activated Receptor (PPAR) in -Se group. The regulatory relationship between ALDBGALG0000005049 and stearoyl-CoA desaturase (SCD), which involved in PPAR pathway was verified. The results also showed that the decreased expressions of SCD, PPARα, PPARβ and PPARγ, and the increased expressions of interleukin (IL)-1β, IL-6, IL-8, and chemokine (C-C motif) ligand 4 (CCL4) along with silencing of ALDBGALG0000005049 in chicken myoblasts. Moreover, increased expressions of IL-1β, IL-6, IL-8, and CCL4 and inflammatory cell infiltration in microstructure of chicken muscles treated with Se deficiency were observed. This study displayed an overview of aberrantly expressed lncRNAs and mRNAs profiles and PPAR pathway, and revealed that downregulation of ALDBGALG0000005049 caused inflammation by regulating SCD in chicken muscle resulted from Se deficiency.

Using RNA sequencing to identify putative competing endogenous RNAs (ceRNAs) potentially regulating fat metabolism in bovine liver

RNA sequencing has been extensively used to study specific gene expression patterns to discover potential key genes related to complex traits of interest in animals. Of note, a new regulatory mechanism builds a large-scale regulatory network among transcriptome, where lncRNAs act as competing endogenous RNAs (ceRNAs) to sponge miRNAs to regulate the expression of miRNA target genes post-transcriptionally. In this study, we sequenced the cDNA and sRNA libraries of nine liver samples from three Holstein cows during dry period, early lactation, and peak of lactation with HiSeq platform. As a result, we identified 665 genes, 57 miRNAs and 33 lncRNAs that displayed differential expression patterns across periods. Subsequently, a total of 41ceRNA pairs (lncRNA-mRNA) sharing 11 miRNAs were constructed including 30 differentially expressed genes. Importantly, 12 among them were presented in our large metabolic networks, and predicted to influence the lipid metabolism through insulin, PI3K-Akt, MAPK, AMPK, mTOR, and PPAR signaling pathways, thus, these genes were considered as the most promising candidates for milk fat formation. To our knowledge, this is first investigation to profile the ceRNA regulatory networks of liver transcriptome that could affect milk fat synthesis in bovine, providing a new view of the regulatory mechanism of RNAs.

Regulation of crucial lncRNAs in differentiation of chicken embryonic stem cells to spermatogonia stem cells

Regulation of crucial lncRNAs involved in differentiation of chicken embryonic stem cells (ESCs) to spermatogonia stem cells (SSCs) was explored by sequencing the transcriptome of ESCs, primordial germ cells (PGCs) and SSCs with RNA-Seq; analytical bioinformatic methods were used to excavate candidate lncRNAs. We detected expression of candidate lncRNAs in ESCs, PGCs and SSCs and forecasted related target genes. Utilizing wego, david and string, function and protein-protein interactions of target genes were analyzed. Finally, based on string analysis, interaction diagrams and relevant signaling pathways were established. Our results indicate a total of 9657 lncRNAs in ESCs, PGCs and SSCs, with 3549 defined as significantly different. We screened 20 candidate lncRNAs, each demonstrating a greater than eight-fold difference in |logFC| value between groups (ESCs vs. PGCs, ESCs vs. SSCs and PGCs vs. SSCs) or specifically expressed in an individual cell type. qRT-PCR results indicated that expression tendencies of candidate lncRNAs were consistent with RNA-Seq. Fifteen cis and four trans target genes were forecasted. Based on wego and string analyses, we found lnc-SSC1, lnc-SSC5, lnc-SSC2 and lnc-ESC2 negatively regulated target genes SUFU, EPHA3, KLF3, ARL3 and TRIM8, whereas SHH, NOTCH, TGF-β, cAMP/cGMP and JAK/STAT signaling pathways were promoted, causing differentiation of ESCs into SSCs. Our findings represent a preliminary unveiling of lncRNA-associated regulatory mechanisms during differentiation of chicken ESCs into SSCs, filling a research void in male germ cell differentiation related to lncRNA. Our results also provide basic information for improving in vitro induction systems for differentiation of chicken ESCs into SSCs.