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

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.

Epigenetic Factor MicroRNAs Likely Mediate Vaccine Protection Efficacy against Lymphomas in Response to Tumor Virus Infection in Chickens through Target Gene Involved Signaling Pathways

Epigenetic factors, including microRNAs (miRNAs), play an important role in affecting gene expression and, therefore, are involved in various biological processes including immunity protection against tumors. Marek's disease (MD) is a highly contagious disease of chickens caused by the MD virus (MDV). MD has been primarily controlled by vaccinations. MD vaccine efficacy might, in part, be dependent on modulations of a complex set of factors including host epigenetic factors. This study was designed to identify differentially expressed miRNAs in the primary lymphoid organ, bursae of Fabricius, in response to MD vaccination followed by MDV challenge in two genetically divergent inbred lines of White Leghorns. Small RNA sequencing and bioinformatic analyses of the small RNA sequence reads identified hundreds of miRNAs among all the treatment groups. A small portion of the identified miRNAs was differentially expressed within each of the four treatment groups, which were HVT or CVI988/Rispens vaccinated line 63-resistant birds and line 72-susceptible birds. A direct comparison between the resistant line 63 and susceptible line 72 groups vaccinated with HVT followed by MDV challenge identified five differentially expressed miRNAs. Gene Ontology analysis of the target genes of those five miRNAs revealed that those target genes, in addition to various GO terms, are involved in multiple signaling pathways including MAPK, TGF-β, ErbB, and EGFR1 signaling pathways. The general functions of those pathways reportedly play important roles in oncogenesis, anti-cancer immunity, cancer cell migration, and metastatic progression. Therefore, it is highly likely that those miRNAs may, in part, influence vaccine protection through the pathways.

Genomic analysis of Nigerian indigenous chickens reveals their genetic diversity and adaptation to heat-stress

Indigenous poultry breeds from Africa can survive in harsh tropical environments (such as long arid seasons, excessive rain and humidity, and extreme heat) and are resilient to disease challenges, but they are not productive compared to their commercial counterparts. Their adaptive characteristics are in response to natural selection or to artificial selection for production traits that have left selection signatures in the genome. Identifying these signatures of positive selection can provide insight into the genetic bases of tropical adaptations observed in indigenous poultry and thereby help to develop robust and high-performing breeds for extreme tropical climates. Here, we present the first large-scale whole-genome sequencing analysis of Nigerian indigenous chickens from different agro-climatic conditions, investigating their genetic diversity and adaptation to tropical hot climates (extreme arid and extreme humid conditions). The study shows a large extant genetic diversity but low level of population differentiation. Using different selection signature analyses, several candidate genes for adaptation were detected, especially in relation to thermotolerance and immune response (e.g., cytochrome P450 2B4-like, TSHR, HSF1, CDC37, SFTPB, HIF3A, SLC44A2, and ILF3 genes). These results have important implications for conserving valuable genetic resources and breeding improvement of chickens for thermotolerance.

A chromosome-level genome assembly for the Silkie chicken resolves complete sequences for key chicken metabolic, reproductive, and immunity genes

A set of high-quality pan-genomes would help identify important genes that are still hidden/incomplete in bird reference genomes. In an attempt to address these issues, we have assembled a de novo chromosome-level reference genome of the Silkie (Gallus gallus domesticus), which is an important avian model for unique traits, like fibromelanosis, with unclear genetic foundation. This Silkie genome includes the complete genomic sequences of well-known, but unresolved, evolutionarily, endocrinologically, and immunologically important genes, including leptin, ovocleidin-17, and tumor-necrosis factor-α. The gap-less and manually annotated MHC (major histocompatibility complex) region possesses 38 recently identified genes, with differentially regulated genes recovered in response to pathogen challenges. We also provide whole-genome methylation and genetic variation maps, and resolve a complex genetic region that may contribute to fibromelanosis in these animals. Finally, we experimentally show leptin binding to the identified leptin receptor in chicken, confirming an active leptin ligand-receptor system. The Silkie genome assembly not only provides a rich data resource for avian genome studies, but also lays a foundation for further functional validation of resolved genes.

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.

Integrated analysis of methylation profiles and transcriptome of MDV-infected chicken spleens reveal hypomethylation of CD4 and HMGB1 genes might promote MD tumorigenesis

Marek's disease (MD) is a lymphoproliferative neoplastic disease caused by Marek's disease virus (MDV). Previous studies have showed that DNA methylation was involved in MD development, but systematic studies are still lacking. Herein, we performed whole genome bisulfite sequencing (WGBS) and RNA-seq in MDV-infected tumorous spleens (IN), noninfected spleens (NoIN), and survivor (SUR) spleens of chickens to identify the genes playing important roles in MD tumor transformation. We generated the first genome-wide DNA methylation profile of MDV-infected, noninfected, and survivor chickens. Combined the WGBS and RNA-Seq, we found that the expression of 25% differential expression genes (DEGs) were significantly correlated with methylation of CpG sites in their gene bodies or promoters. Further, we focused on the DEGs with differentially methylated regions (DMRs) on genes' body and promoter, and it showed the expression of 60% DEGs were significantly correlated with methylation of CpG sites in DMRs. Finally, we identified 8 genes, including CD4, CTLA4, DTL, HMGB1, LGMN, NUP210, RAD52, and ZAP70, and their expression was negatively correlated with methylation of DMRs in their promoters in both IN vs. NoIN and IN vs. SUR. These 8 genes showed specifically high expression in IN groups and clustered in module turquoise analyzed by WGCNA. Out of 8 genes, CD4 and HMGB1 were drop in QTLs associated with MD resistance. Thus, we overexpressed the 2 genes to simulate their high expression in the IN group and found they significantly promoted MDCC-MSB-1 cell proliferation, which revealed they might play promoting roles in MD tumorigenesis in IN due to their high expression induced by hypomethylation.

A comprehensive analysis of avian lymphoid leukosis-like lymphoma transcriptomes including identification of LncRNAs and the expression profiles

Avian lymphoid leukosis-like (LL-like) lymphoma has been observed in some experimental and commercial lines of chickens that are free of exogenous avian leukosis virus. Reported cases of avian lymphoid leukosis-like lymphoma incidences in the susceptible chickens are relatively low, but the apathogenic subgroup E avian leukosis virus (ALV-E) and the Marek’s disease vaccine, SB-1, significantly escalate the disease incidence in the susceptible chickens. However, the underlying mechanism of tumorigenesis is poorly understood. In this study, we bioinformatically analyzed the deep RNA sequences of 6 lymphoid leukosis-like lymphoma samples, collected from susceptible chickens post both ALV-E and SB-1 inoculation, and identified a total of 1,692 novel long non-coding RNAs (lncRNAs). Thirty-nine of those novel lncRNAs were detected with altered expression in the LL-like tumors. In addition, 13 lncRNAs whose neighboring genes also showed differentially expression and 2 conserved novel lncRNAs, XLOC_001407 and XLOC_022595, may have previously un-appreciated roles in tumor development in human. Furthermore, 14 lncRNAs, especially XLOC_004542, exhibited strong potential as competing endogenous RNAs via sponging miRNAs. The analysis also showed that ALV subgroup E viral gene Gag/Gag-pol and the MD vaccine SB-1 viral gene R-LORF1 and ORF413 were particularly detectable in the LL-like tumor samples. In addition, we discovered 982 novel lncRNAs that were absent in the current annotation of chicken genome and 39 of them were aberrantly expressed in the tumors. This is the first time that lncRNA signature is identified in avian lymphoid leukosis-like lymphoma and suggests the epigenetic factor, lncRNA, is involved with the avian lymphoid leukosis-like lymphoma formation and development in susceptible chickens. Further studies to elucidate the genetic and epigenetic mechanisms underlying the avian lymphoid leukosis-like lymphoma is indeed warranted.

Identification of Long Noncoding RNAs Involved in Eyelid Pigmentation of Hereford Cattle

Several ocular pathologies in cattle, such as ocular squamous cell carcinoma and infectious keratoconjunctivitis, have been associated with low pigmentation of the eyelids. The main objective of this study was to analyze the transcriptome of eyelid skin in Hereford cattle using strand-specific RNA sequencing technology to characterize and identify long noncoding RNAs (lncRNAs). We compared the expression of lncRNAs between pigmented and unpigmented eyelids and analyzed the interaction of lncRNAs and putative target genes to reveal the genetic basis underlying eyelid pigmentation in cattle. We predicted 4,937 putative lncRNAs mapped to the bovine reference genome, enriching the catalog of lncRNAs in Bos taurus. We found 27 differentially expressed lncRNAs between pigmented and unpigmented eyelids, suggesting their involvement in eyelid pigmentation. In addition, we revealed potential links between some significant differentially expressed lncRNAs and target mRNAs involved in the immune response and pigmentation. Overall, this study expands the catalog of lncRNAs in cattle and contributes to a better understanding of the biology of eyelid pigmentation.

Long Noncoding RNA Expression Rofiles Elucidate the Potential Roles of lncRNA- XR_003496198 in Duck Hepatitis A Virus Type 1 Infection

Duck hepatitis A virus type 1 (DHAV-1) is a highly lethal virus that severely affects the duck industry worldwide. Long noncoding RNAs (lncRNAs) exert crucial roles in pathogen attacks. Here, we conducted deep transcriptome analysis to investigate the dynamic changes of host lncRNAs profiles in DHAV-1-infected duck embryo fibroblasts. We identified 16,589 lncRNAs in total and characterized their genomic features. Moreover, 772 and 616 differentially expressed lncRNAs (DELs) were screened at 12 and 24 h post-infection. Additionally, we predicted the DELs’ cis- and trans-target genes and constructed lncRNA-target genes regulatory networks. Functional annotation analyses indicated that the putative target genes of DELs participated in diverse vital biological processed, including immune responses, cellular metabolism, and autophagy. For example, we confirmed the dysregulation of pattern recognition receptors (TLR3, RIG-I, MDA5, LGP2, cGAS), signal transducers (STAT1), transcription factors (IRF7), immune response mediators (IL6, IL10, TRIM25, TRIM35, TRIM60, IFITM1, IFITM3, IFITM5), and autophagy-related genes (ULK1, ULK2, EIF4EBP2) using RT-qPCR. Finally, we confirmed that one DHAV-1 induced lncRNA-XR_003496198 is likely to inhibit DHAV-1 replication in DEFs. Our study comprehensively analyzed the lncRNA profiles upon DHAV-1 infection and screened the target genes involved in the innate immune response and autophagy signaling pathway, thereby revealing the essential roles of duck lncRNAs and broadening our understanding of host-virus interactions.

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.

Transcriptome Analysis Revealed Potential Mechanisms of Resistance to Trichomoniasis gallinae Infection in Pigeon (Columba livia)

Trichomoniasis gallinae (T. gallinae) is one of the most pathogenic parasites in pigeon, particularly in squabs. Oral cavity is the main site for the host-parasite interaction. Herein, we used RNA-sequencing technology to characterize lncRNA and mRNA profiles and compared transcriptomic dynamics of squabs, including four susceptible birds (S) from infected group, four tolerant birds (T) without parasites after T. gallinae infection, and three birds from uninfected group (N), to understand molecular mechanisms underlying host resistance to this parasite. We identified 29,809 putative lncRNAs and characterized their genomic features subsequently. Differentially expressed (DE) genes, DE-lncRNAs and cis/trans target genes of DE-lncRNAs were further compared among the three groups. The KEGG analysis indicated that specific intergroup DEGs were involved in carbon metabolism (S vs. T), metabolic pathways (N vs. T) and focal adhesion pathway (N vs. S), respectively. Whereas, the cis/trans genes of DE-lncRNAs were enriched in cytokine-cytokine receptor interaction, toll-like receptor signaling pathway, p53 signaling pathway and insulin signaling pathway, which play crucial roles in immune system of the host animal. This suggests T. gallinae invasion in pigeon mouth may modulate lncRNAs expression and their target genes. Moreover, co-expression analysis identified crucial lncRNA-mRNA interaction networks. Several DE-lncRNAs including MSTRG.82272.3, MSTRG.114849.42, MSTRG.39405.36, MSTRG.3338.5, and MSTRG.105872.2 targeted methylation and immune-related genes, such as JCHAIN, IL18BP, ANGPT1, TMRT10C, SAMD9L, and SOCS3. This implied that DE-lncRNAs exert critical influence on T. gallinae infections. The quantitative exploration of host transcriptome changes induced by T. gallinae infection broaden both transcriptomic and epigenetic insights into T. gallinae resistance and its pathological mechanism.

Gga-miR-181a modulates ANP32A expression and inhibits MDCC-MSB-1 cell

Marek's disease (MD), a highly contagious T cell lymphoid neoplasia disease of chickens, causes huge economic losses to the poultry industry. It is the only one tumor disease which can be prevented by vaccine in chickens; therefore, MD is considered to be an excellent model to study the pathogenesis of virus-induced cancer. Recently, abundant evidences have verified that miRNAs are regulators in the process of neoplastic transformation. In our previous study on miRNome analysis of MDV-induced lymphoma in chicken, we found that gga-miR-181a was downregulated drastically in MDV-infected spleens. To further investigate the role of gga-miR-181a in MDV-induced lymphomagenesis, we performed cell migration assay, and the results suggested that gga-miR-181a suppressed the migration of MDV-transformed lymphoid cell (MSB-1). Subsequently, luciferase reporter gene assay revealed that acidic nuclear phosphoprotein 32A (ANP32A) was a functional target gene of gga-miR181a. Real-time PCR and western blot assay showed that the mRNA and protein levels of ANP32A were downregulated in gga-miR-181a mimic group at 48-h and 96-h post-transfection, respectively, indicating that ANP32A was modulated by gga-miR-181a. All the results suggested that gga-miR-181a was an inhibitor in MSB-1 cell migration. ANP32A was a direct target gene of gga-miR-181a and they were implicated in MD lymphoma tumorigenesis.

A signature of estimate-stromal-immune score-based genes associated with the prognosis of lung adenocarcinoma

Background

Immune and stromal component evaluation is necessary to establish accurate prognostic markers for the prediction of clinical outcomes in lung adenocarcinoma (LUAD). We aimed to develop a gene signature based on the Estimation of STromal and Immune cells in MAlignant Tumor tissues using Expression data (ESTIMATE)-stromal-immune score in LUAD.

Methods

The transcriptomic profiles of patients with LUAD were obtained from The Cancer Genome Atlas (TCGA), and the immune and stromal scores were derived using the ESTIMATE algorithm. The prognostic signature genes were selected from the differentially expressed genes (DEGs) using the robust partial likelihood-based cox proportional hazards regression method. The negative log-likelihood and the Akaike Information Criterion (AIC) were used to identify the optimal gene signature. The validation was carried out in 2 independent datasets from the Gene Expression Omnibus (GSE68571 and GSE72094).

Results

Patients with high ESTIMATE, stromal, and immune scores had better overall survivals (P=0.0035, P=0.066, and P=0.0077). The expression of thirty-seven genes was related to ESTIMATE-stromal-immune score. A risk stratification model was developed based on a gene signature containing CD74, JCHAIN, and PTGDS. The ESTIMATE-stromal-immune risk score was revealed to be a prognostic factor (P=0.009) after multivariate analysis. Four groups were classified based on this risk stratification model, yielding increasing survival outcomes (log-rank test, P=0.0051). This risk stratification model and other clinicopathological factors were combined to generate a nomogram. The calibration curves showed perfect agreement between the nomogram-predicted outcomes and those actually observed. Similar observations were made in 2 independent cohorts.

Conclusions

The gene signature based on the ESTIMATE-stromal-immune score could predict the prognosis of patients with LUAD.

Identification and analysis of long non-coding RNAs and mRNAs in chicken macrophages infected with avian infectious bronchitis coronavirus

BACKGROUND

Avian infectious bronchitis virus (IBV) is a gamma coronavirus that severely affects the poultry industry worldwide. Long non-coding RNAs (lncRNAs), a subset of non-coding RNAs with a length of more than 200 nucleotides, have been recently recognized as pivotal factors in the pathogenesis of viral infections. However, little is known about the function of lncRNAs in host cultured cells in response to IBV infection.

RESULTS

We used next-generation high throughput sequencing to reveal the expression profiles of mRNAs and lncRNAs in IBV-infected HD11 cells. Compared with the uninfected cells, we identified 153 differentially expressed (DE) mRNAs (106 up-regulated mRNAs, 47 down-regulated mRNAs) and 181 DE lncRNAs (59 up-regulated lncRNAs, 122 down-regulated lncRNAs) in IBV-infected HD11 cells. Moreover, gene ontology (GO) and pathway enrichment analyses indicated that DE mRNAs and lncRNAs were mainly involved in cellular innate immunity, amino acid metabolism, and nucleic acid metabolism. In addition, 2640 novel chicken lncRNAs were identified, and a competing endogenous RNA (ceRNAs) network centered on gga-miR-30d and miR-146a-5p was established.

CONCLUSIONS

We identified expression profiles of mRNAs and lncRNAs during IBV infection that provided new insights into the pathogenesis of IBV.

Identification and Expression Analysis of Long Non-coding RNA in Large Yellow Croaker (Larimichthys crocea) in Response to Cryptocaryon irritans Infection

Large-scale transcription studies have revealed numerous lncRNAs (long non-coding RNAs). lncRNAs have been proposed to participate in the regulation of a diverse range of biological processes, including transcriptional regulation. Although lncRNAs have attracted increasing attention, the studies in large yellow croaker (Larimichthys crocea) are still rare, and they lack systematic analysis. In this study, 101 RNA-seq datasets varied in ages, sexes, and tissues were retrieved from the NCBI database to generate a comprehensive catalog of large yellow croaker transcriptome database. A set of 14,599 high-confidence lncRNAs from 13,673 loci were identified and characterized. Furthermore, RNA-seq datasets obtained from the infection of C. irritans were employed to investigate the differential expression pattern of lncRNAs and analyze potential biological functions. A total of 77 differentially expressed lncRNAs targeting to 567 protein-coding genes were identified by using expression analysis. Several immune genes, including TLR5, CD2AP, and MMP9, were highlighted. With GO enrichment and KEGG pathway analysis, the immune-related terms or pathways were enriched. This study created a comprehensive dataset of lncRNAs for large yellow croaker, which would be helpful for the researches of functional roles of lncRNAs in large yellow croaker.

Gene regulation network analysis reveals core genes associated with survival in glioblastoma multiforme

Glioblastoma multiforme (GBM) is a very serious mortality of central nervous system cancer. The microarray data from GSE2223, GSE4058, GSE4290, GSE13276, GSE68848 and GSE70231 (389 GBM tumour and 67 normal tissues) and the RNA-seq data from TCGA-GBM dataset (169 GBM and five normal samples) were chosen to find differentially expressed genes (DEGs). RRA (Robust rank aggregation) method was used to integrate seven datasets and calculate 133 DEGs (82 up-regulated and 51 down-regulated genes). Subsequently, through the PPI (protein-protein interaction) network and MCODE/ cytoHubba methods, we finally filtered out ten hub genes, including FOXM1, CDK4, TOP2A, RRM2, MYBL2, MCM2, CDC20, CCNB2, MYC and EZH2, from the whole network. Functional enrichment analyses of DEGs were conducted to show that these hub genes were enriched in various cancer-related functions and pathways significantly. We also selected CCNB2, CDC20 and MYBL2 as core biomarkers, and further validated them in CGGA, HPA and CCLE database, suggesting that these three core hub genes may be involved in the origin of GBM. All these potential biomarkers for GBM might be helpful for illustrating the important role of molecular mechanisms of tumorigenesis in the diagnosis, prognosis and targeted therapy of GBM cancer.

Mapping QTL Associated with Resistance to Avian Oncogenic Marek's Disease Virus (MDV) Reveals Major Candidate Genes and Variants

Marek’s disease (MD) represents a significant global economic and animal welfare issue. Marek’s disease virus (MDV) is a highly contagious oncogenic and highly immune-suppressive α-herpes virus, which infects chickens, causing neurological effects and tumour formation. Though partially controlled by vaccination, MD continues to have a profound impact on animal health and on the poultry industry. Genetic selection provides an alternative and complementary method to vaccination. However, even after years of study, the genetic mechanisms underlying resistance to MDV remain poorly understood. The Major Histocompatability Complex (MHC) is known to play a role in disease resistance, along with a handful of other non-MHC genes. In this study, one of the largest to date, we used a multi-facetted approach to identify quantitative trait locus regions (QTLR) influencing resistance to MDV, including an F₆ population from a full-sib advanced intercross line (FSIL) between two elite commercial layer lines differing in resistance to MDV, RNA-seq information from virus challenged chicks, and genome wide association study (GWAS) from multiple commercial lines. Candidate genomic elements residing in the QTLR were further tested for association with offspring mortality in the face of MDV challenge in eight pure lines of elite egg-layer birds. Thirty-eight QTLR were found on 19 chicken chromosomes. Candidate genes, microRNAs, long non-coding RNAs and potentially functional mutations were identified in these regions. Association tests were carried out in 26 of the QTLR, using eight pure lines of elite egg-layer birds. Numerous candidate genomic elements were strongly associated with MD resistance. Genomic regions significantly associated with resistance to MDV were mapped and candidate genes identified. Various QTLR elements were shown to have a strong genetic association with resistance. These results provide a large number of significant targets for mitigating the effects of MDV infection on both poultry health and the economy, whether by means of selective breeding, improved vaccine design, or gene-editing technologies.

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.

Review: Long non-coding RNA in livestock

Less than 2% of mammalian genomes code for proteins, but 'the majority of its bases can be found in primary transcripts' - a phenomenon termed the pervasive transcription, which was first reported in 2007. Even though most of the transcripts do not code for proteins, they play a variety of biological functions, with regulation of gene expression appearing as the most common one. Those transcripts are divided into two groups based on their length: small non-coding RNAs, which are maximally 200 bp long, and long non-coding RNAs (lncRNAs), which are longer than 200 nucleotides. The advances in next-generation sequencing methods provided a new possibility of investigating the full set of RNA molecules in the cell. In this review, we summarized the current state of knowledge on lncRNAs in three major livestock species - Sus scrofa, Bos taurus and Gallus gallus, based on the literature and the content of biological databases. In the NONCODE database, the largest number of identified lncRNA transcripts is available for pigs, but cattle have the largest number of lncRNA genes. Poultry is represented by less than a half of records. Genomic annotation of lncRNAs showed that the majority of them are assigned to introns (pig, poultry) or intergenic (cattle). The comparison with well-annotated human and mouse genomes indicates that such annotation is a result of lack of proper lncRNA annotation data. Since lncRNAs play an important role in genomic studies, their characterization in farm animals' genomes is critical in bridging the gap between genotype and phenotype.

Genome-wide analysis of differentially expressed profiles of mRNAs, lncRNAs and circRNAs in chickens during Eimeria necatrix infection

BACKGROUND

Eimeria necatrix, the most highly pathogenic coccidian in chicken small intestines, can cause high morbidity and mortality in susceptible birds and devastating economic losses in poultry production, but the underlying molecular mechanisms in interaction between chicken and E. necatrix are not entirely revealed. Accumulating evidence shows that the long-non-coding RNAs (lncRNAs) and circular RNAs (circRNAs) are key regulators in various infectious diseases. However, the expression profiles and roles of these two non-coding RNAs (ncRNAs) during E. necatrix infection are still unclear.

METHODS

The expression profiles of mRNAs, lncRNAs and circRNAs in mid-segments of chicken small intestines at 108 h post-infection (pi) with E. necatrix were analyzed by using the RNA-seq technique.

RESULTS

After strict filtering of raw data, we putatively identified 49,183 mRNAs, 818 lncRNAs and 4153 circRNAs. The obtained lncRNAs were classified into four types, including 228 (27.87%) intergenic, 67 (8.19%) intronic, 166 (20.29%) anti-sense and 357 (43.64%) sense-overlapping lncRNAs; of these, 571 were found to be novel. Five types were also predicted for putative circRNAs, including 180 exonic, 54 intronic, 113 antisense, 109 intergenic and 3697 sense-overlapping circRNAs. Eimeria necatrix infection significantly altered the expression of 1543 mRNAs (707 upregulated and 836 downregulated), 95 lncRNAs (49 upregulated and 46 downregulated) and 13 circRNAs (9 upregulated and 4 downregulated). Target predictions revealed that 38 aberrantly expressed lncRNAs would cis-regulate 73 mRNAs, and 1453 mRNAs could be trans-regulated by 87 differentially regulated lncRNAs. Additionally, 109 potential sponging miRNAs were also identified for 9 circRNAs. GO and KEGG enrichment analysis of target mRNAs for lncRNAs, and sponging miRNA targets and source genes for circRNAs identified associations of both lncRNAs and circRNAs with host immune defense and pathogenesis during E. necatrix infection.

CONCLUSIONS

To the best of our knowledge, the present study provides the first genome-wide analysis of mRNAs, lncRNAs and circRNAs in chicken small intestines infected with E. necatrix. The obtained data will offer novel clues for exploring the interaction mechanisms between chickens and Eimeria spp.

Genome-wide analysis of circular RNAs involved in Marek's disease tumourigenesis in chickens

Marek's disease (MD), induced by Marek's disease virus (MDV), is a lymphotropic neoplastic disease and causes huge economic losses to the poultry industry. Non-coding RNAs (ncRNAs) play important regulatory roles in disease pathogenesis. To investigate host circular RNA (circRNA) and microRNA (miRNA) expression profile, RNA sequencing was performed in tumourous spleens (TS), spleens from the survivors (SS) without any lesion after MDV infection, and non-infected chicken spleens (NS). A total of 2,169 circRNAs were identified and more than 80% of circRNAs were derived from exon. The flanking introns of 1,744 exonic circRNAs possessed 579 reverse complementary matches (RCMs), which mainly overlapped with chicken repeat 1 family (CR1F). It suggested that CR1F mediated the cyclization of exons by intron pairing. Out of 2,169 circRNAs, 113 were differentially expressed circRNAs (DECs). The Q-PCR and Rnase R digestion experiments showed circRNA possessed high stability compared with their linear RNAs. Integrated with previous transcriptome data, we profiled regulatory networks of circRNA/long non-coding RNA (lncRNA)-miRNA-mRNA. Extensive competing endogenous RNA (ceRNA) networks were predicted to be involved in MD tumourigenesis. Interestingly, circZMYM3, an intronic circRNA, interacted with seven miRNAs which targeted some immune genes, such as SWAP70 and CCL4. Gga-miR-155 not only interacted with circGTDC1 and circMYO1B, but also targeted immune-related genes, such as GATA4, which indicated the roles of non-coding RNAs played to mediate immune responsive genes. Collectively, this is the first study that integrated RNA expression profiles in MD model. Our results provided comprehensive interactions of ncRNAs and mRNA in MD tumourigenesis.

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.