Gene expression profile and long non-coding RNA analysis, using RNA-Seq, in chicken embryonic fibroblast cells infected by avian leukosis virus J

Transcriptomic and Translatomic Analyses Reveal Insights into the Signaling Pathways of the Innate Immune Response in the Spleens of SPF Chickens Infected with Avian Reovirus

Avian reovirus (ARV) infection is prevalent in farmed poultry and causes viral arthritis and severe immunosuppression. The spleen plays a very important part in protecting hosts against infectious pathogens. In this research, transcriptome and translatome sequencing technology were combined to investigate the mechanisms of transcriptional and translational regulation in the spleen after ARV infection. On a genome-wide scale, ARV infection can significantly reduce the translation efficiency (TE) of splenic genes. Differentially expressed translational efficiency genes (DTEGs) were identified, including 15 upregulated DTEGs and 396 downregulated DTEGs. These DTEGs were mainly enriched in immune regulation signaling pathways, which indicates that ARV infection reduces the innate immune response in the spleen. In addition, combined analyses revealed that the innate immune response involves the effects of transcriptional and translational regulation. Moreover, we discovered the key gene IL4I1, the most significantly upregulated gene at both the transcriptional and translational levels. Further studies in DF1 cells showed that overexpression of IL4I1 could inhibit the replication of ARV, while inhibiting the expression of endogenous IL4I1 with siRNA promoted the replication of ARV. Overexpression of IL4I1 significantly downregulated the mRNA expression of IFN-β, LGP2, TBK1 and NF-κB; however, the expression of these genes was significantly upregulated after inhibition of IL4I1, suggesting that IL4I1 may be a negative feedback effect of innate immune signaling pathways. In addition, there may be an interaction between IL4I1 and ARV σA protein, and we speculate that the IL4I1 protein plays a regulatory role by interacting with the σA protein. This study not only provides a new perspective on the regulatory mechanisms of the innate immune response after ARV infection but also enriches the knowledge of the host defense mechanisms against ARV invasion and the outcome of ARV evasion of the host's innate immune response.

Long Noncoding RNA AROD Inhibits Host Antiviral Innate Immunity via the miR-324-5p-CUEDC2 Axis

Long noncoding RNAs (lncRNAs) are a class of noncoding RNAs that are involved in multiple biological processes. Here, we report a mechanism through which the lnc-AROD-miR-324-5p-CUEDC2 axis regulates the host innate immune response, using influenza A virus (IAV) as a model. We identified that host lnc-AROD without protein-coding capability is composed of 975 nucleotides. Moreover, lnc-AROD inhibited interferon-β expression, as well as interferon-stimulated genes ISG15 and MxA. Furthermore, in vivo assays confirmed that lnc-AROD overexpression increased flu virus pathogenicity and mortality in mice. Mechanistically, lnc-AROD interacted with miR-324-5p, leading to decreased binding of miR-324-5p to CUEDC2. Collectively, our findings demonstrated that lnc-AROD is a critical regulator of the host antiviral response via the miR-324-5p-CUEDC2 axis, and lnc-AROD functions as competing endogenous RNA. Our results also provided evidence that lnc-AROD serves as an inhibitor of the antiviral immune response and may represent a potential drug target. IMPORTANCE lnc-AROD is a potential diagnostic and discriminative biomarker for different cancers. However, so far the mechanisms of lnc-AROD regulating virus replication are not well understood. In this study, we identified that lnc-AROD is downregulated during RNA virus infection. We demonstrated that lnc-AROD enhanced CUEDC2 expression, which in turn inhibited innate immunity and favored IAV replication. Our studies indicated that lnc-AROD functions as a competing endogenous RNA that binds miR-324-5p and reduces its inhibitory effect on CUEDC2. Taken together, our findings reveal that lnc-AROD plays an important role during the host antiviral immune response.

Long noncoding RNAs expression profile of RIP2 knockdown in chicken HD11 macrophages associated with avian pathogenic E. coli (APEC) infection

Avian pathogenic E. coli (APEC) has been detected to cause many acute and chronic diseases, resulting in huge economic losses to the poultry industry. Previous experiments have identified the effect of receptor interacting serine/threonine kinase 2 (RIP2) gene in APEC infection. Moreover, increasing evidence indicates that long noncoding RNAs (lncRNAs) play important roles in the anti-bacteria responses. However, little is known about the functions of lncRNAs, especially related to RIP2, in response to APEC. Therefore, we tried to reveal lncRNAs potentially involved in the immune and inflammatory response against APEC infection, with a particular focus on those possibly correlated with RIP2. A total of 1856 and 1373 differentially expressed (DE) lncRNAs were identified in knockdown of RIP2 cells following APEC infection (shRIP2+APEC) vs. APEC and shRIP2 vs. wild type cells (WT), respectively, which were mainly enriched in lysosome, phagosome, NOD-like receptor signaling pathway, TGF-beta signaling pathway. Significantly, TCONS_00009695 regulated by RIP2 could directly alter the expression of target BIRC3 to modulate cytokines and to participate in immune and inflammatory response against APEC infection. Our findings aid to a better understanding of host responses to APEC infection and provide new directions for understanding the potential association between lncRNAs and APEC pathogenesis.

Long Non-Coding RNA BNIP3 Inhibited the Proliferation of Bovine Intramuscular Preadipocytes via Cell Cycle

The intramuscular fat (or marbling fat) content is an essential economic trait of beef cattle and improves the flavor and palatability of meat. Several studies have highlighted the correlation between long non-coding RNAs (lncRNAs) and intramuscular fat development; however, the precise molecular mechanism remains unknown. Previously, through a high-throughput sequencing analysis, we found a lncRNA and named it a long non-coding RNA BNIP3 (lncBNIP3). The 5' RACE and 3' RACE explored 1945 bp total length of lncBNIP3, including 1621 bp of 5'RACE, and 464 bp of 3'RACE. The nucleoplasmic separation and FISH results explored the nuclear localization of lncBNIP3. Moreover, the tissue expression of lncBNIP3 was higher in the longissimus dorsi muscle, followed by intramuscular fat. Furthermore, down-regulation of lncBNIP3 increased the 5-Ethynyl-2'- deoxyuridine (EdU)-EdU-positive cells. The flow cytometry results showed that the number of cells in the S phase was significantly higher in preadipocytes transfected with si-lncBNIP3 than in the control group (si-NC). Similarly, CCK8 results showed that the number of cells after transfection of si-lncBNIP3 was significantly higher than in the control group. In addition, the mRNA expressions of proliferative marker genes CyclinB1 (CCNB1) and Proliferating Cell Nuclear Antigen (PCNA) in the si-lncBNIP3 group were significantly higher than in the control group. The Western Blot (WB) results also showed that the protein expression level of PCNA transfection of si-lncBNIP3 was significantly higher than in the control group. Similarly, the enrichment of lncBNIP3 significantly decreased the EdU-positive cells in the bovine preadipocytes. The results of flow cytometry and CCK8 assay also showed that overexpression of lncBNIP3 inhibited the proliferation of bovine preadipocytes. In addition, the overexpression of lncBNIP3 significantly inhibited the mRNA expressions of CCNB1 and PCNA. The WB results showed that the overexpression of lncBNIP3 significantly inhibited the expression of the CCNB1 protein level. To further explore the mechanism of lncBNIP3 on the proliferation of intramuscular preadipocytes, RNA-seq was performed after interference with si-lncBNIP3, and 660 differentially expressed genes (DEGs) were found, including 417 up-regulated DEGs and 243 down-regulated DEGs. The KEGG pathway analysis showed that the cell cycle was the most significant pathway for the functional enrichment of DEGs, followed by the DNA replication pathway. The RT-qPCR quantified the expression of twenty DEGs in the cell cycle. Therefore, we speculated that lncBNIP3 regulated intramuscular preadipocyte proliferation through the cell cycle and DNA replication pathways. To further confirm this hypothesis, the cell cycle inhibitor Ara-C was used to inhibit DNA replication of the S phase in intramuscular preadipocytes. Herein, Ara-C and si-lncBNIP3 were simultaneously added to the preadipocytes, and the CCK8, flow cytometry, and EdU assays were performed. The results showed that the si-lncBNIP3 could rescue the inhibitory effect of Ara-C in the bovine preadipocyte proliferation. In addition, lncBNIP3 could bind to the promoter of cell division control protein 6 (CDC6), and down-regulation of lncBNIP3 promoted the transcription activity and the expression of CDC6. Therefore, the inhibitory effect of lncBNIP3 on cell proliferation might be understood through the cell cycle pathway and CDC6 expression. This study provided a valuable lncRNA with functional roles in intramuscular fat accumulation and revealed new strategies for improving beef quality.

The origin and loss of interferon regulatory factor 10 (IRF10) in different lineages of vertebrates

The vertebrate IFN regulatory factor (IRF) family consists of 11 members that exert distinct roles in a variety of biological processes, including antiviral defense, regulation of cell proliferation, differentiation and apoptosis. Of these, IRF10 is widely present in different vertebrate lineages, but appears to have been lost in primates and rodents. To understand the evolutionary occurrence of IRF10, we performed comparative analyses of currently available genomic data in a taxonomically diverse set of vertebrates, and found that IRF10 originated after the divergence of chondrichthyans from gnathostomes. Phylogenetically, vertebrate IRF10 is much more closely related to IRF4 than to IRF8 or IRF9, although these four IRFs may have a common ancestor. In addition, the loss of IRF10 in Euarchontoglires might be resulted from mutation accumulation, and the rate of mutations in rodents appears to be higher than in the primate lineage. In primates, the gene-disruptive mutations may have occurred at least prior to the separation of new world monkey and old world primates, roughly 40 million years ago. Overall, we propose a detailed evolutionary scenario for IRF10, which may help us understand the evolutionary mechanisms in the expansion and contraction of the IRF family.

Global Maps of Avian Leukosis Viruses: Research Trends and Themes Based on Networking

Avian leukosis virus (ALV) has a tremendous adverse impact on the poultry industry. Since its discovery, research on different aspects of ALV have been published. Due to the vast academic emphasis and economic importance of the ALV infection in poultry worldwide, this bibliometric analysis explored the scientific output associated with ALV utilizing the Web of Science (Core Collection) database. The relevant data were collected using the search query "AVIAN LEUKOSIS VIRUS", further refined by document types (article, book chapter, and proceedings paper). Finally, 1060 items with full records were imported in Plaintext and tab-delimited formats. The data analysis was carried out using MS Excel, VOS viewer, and R (Biblioshiny) software. Chinese and American research institutions produced the majority of papers during study time period. The Journal of Virology and Avian Diseases appeared as the favorite journal/source for publications. Apart from the avian leukosis virus and the ALV-J, the important keywords mentioned included avian leukosis virus subgroup j, chicken, and retrovirus. The analysis revealed substantial findings on ALV research, with a strong research response from the USA and China.

Transcription analysis of chicken embryo fibroblast cells infected with the recombinant avian leukosis virus isolate GX14FF03

Infection with recombinant avian leukosis virus (ALV) has previously been linked to malignancies and immunosuppression. However, the processes behind the unique pathophysiology of recombinant ALV are poorly understood. In this study, we analyzed gene expression patterns in chicken fibroblast cells (CEFs) infected with the recombinant ALV isolate GX14FF03 and used the RNA-seq technique to perform a complete analysis of the transcribed mRNAs. A total of 907 significant differentially expressed genes (SDEGs) were identified. Among these SDEGs, the most significantly upregulated gene was interleukin 8-like 1 (IL8L1), while the most significantly downregulated gene was fibroblast growth factor 16 (FGF16). The 907 SDGEs were highly enriched (p < 0.05) for 252 Gene Ontology (GO) terms, including 197 BP, 3 CC, and 52 MF. According to KEGG data analysis, SDEGs are implicated in eight significant pathways (p < 0.05). Furthermore, protein-protein interaction (PPI) network analysis revealed that IL8L1 interacts with 17 genes. These findings shed light on the molecular mechanisms involved in recombinant ALV infection by showing the mRNA expression profile in CEFs infected with GX14FF03 virus.

An Endogenous Retroviral LTR-Derived Long Noncoding RNA lnc-LTR5B Interacts With BiP to Modulate ALV-J Replication in Chicken Cells

Infection with the avian leukosis virus subgroup J (ALV-J) impairs host genes and facilitates the establishment of chronic infection and the viral life cycle. However, the involvement of long noncoding RNAs (lncRNAs) in ALV-J infection remains largely unknown. In this study, we identified a novel chicken lncRNA derived from LTR5B of the ERV-L family (namely lnc-LTR5B), which is significantly downregulated in ALV-J infected cells. lnc-LTR5B was localized in the cytoplasm and was relatively high expressed in the chicken lung and liver. Notably, the replication of ALV-J was inhibited by the overexpression of lnc-LTR5B but enhanced when lnc-LTR5B expression was knocked down. We further confirmed that lnc-LTR5B could bind to the binding immunoglobulin protein (BiP), a master regulator of endoplasmic reticulum (ER) function. Mechanistically, lnc-LTR5B serves as a competing endogenous RNA for BiP, restricting its physical availability. Upon ALV-J infection, the reduction of lnc-LTR5B released BiP, which facilitated its translocation to the cell surface. This is crucial for ALV-J entry as well as pro-survival signaling. In conclusion, we identified an endogenous retroviral LTR-activated lnc-LTR5B that is involved in regulating the cell surface translocation of BiP, and such regulatory machinery can be exploited by ALV-J to complete its life cycle and propagate.

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.

Tryptophan Metabolites and Aryl Hydrocarbon Receptor in Severe Acute Respiratory Syndrome, Coronavirus-2 (SARS-CoV-2) Pathophysiology

The metabolism of tryptophan is intimately associated with the differential regulation of diverse physiological processes, including in the regulation of responses to severe acute respiratory syndrome, coronavirus-2 (SARS-CoV-2) infection that underpins the COVID-19 pandemic. Two important products of tryptophan metabolism, viz kynurenine and interleukin (IL)4-inducible1 (IL41)-driven indole 3 pyruvate (I3P), activate the aryl hydrocarbon receptor (AhR), thereby altering the nature of immune responses to SARS-CoV-2 infection. AhR activation dysregulates the initial pro-inflammatory cytokines production driven by neutrophils, macrophages, and mast cells, whilst AhR activation suppresses the endogenous antiviral responses of natural killer cells and CD8+ T cells. Such immune responses become further dysregulated by the increased and prolonged pro-inflammatory cytokine suppression of pineal melatonin production coupled to increased gut dysbiosis and gut permeability. The suppression of pineal melatonin and gut microbiome-derived butyrate, coupled to an increase in circulating lipopolysaccharide (LPS) further dysregulates the immune response. The AhR mediates its effects via alterations in the regulation of mitochondrial function in immune cells. The increased risk of severe/fatal SARS-CoV-2 infection by high risk conditions, such as elderly age, obesity, and diabetes are mediated by these conditions having expression levels of melatonin, AhR, butyrate, and LPS that are closer to those driven by SARS-CoV-2 infection. This has a number of future research and treatment implications, including the utilization of melatonin and nutraceuticals that inhibit the AhR, including the polyphenols, epigallocatechin gallate (EGCG), and resveratrol.

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.

Interleukin 4 inducible 1 gene (IL4I1) is induced in chicken phagocytes by Salmonella Enteritidis infection

In attempt to identify genes that are induced in chickens by Salmonella Enteritidis we identified a new highly inducible gene, interleukin 4 induced 1 gene (IL4I1). IL4I1 reached its peak expression (458× induction) in the cecum of newly hatched chickens 4 days post-infection and remained upregulated for an additional 10 days. IL4I1 was expressed and induced in macrophages and granulocytes, both at the mRNA and protein level. IL4I1 was expressed and induced also in CD4 and γδ T-lymphocytes though at a 50-fold lower level than in phagocytes. Expression of IL4I1 was not detected in CD8 T lymphocytes or B lymphocytes. Mutation of IL4I1 in chicken HD11 macrophages did not affect their bactericidal capacity against S. Enteritidis but negatively affected their oxidative burst after PMA stimulation. We therefore propose that IL4I1 is not directly involved in bactericidal activity of phagocytes and, instead, it is likely involved in the control of inflammatory response and signaling to T and B lymphocytes.

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.

LncRNA expression profile during autophagy and Malat1 function in macrophages

Long noncoding RNAs (lncRNAs) are a class of functional non-coding transcripts that are longer than 200 nt and regulate gene expression via diverse mechanisms in eukaryotes. In fact, they have emerged as critical epigenetic and transcriptional regulators of autophagy in mammals in response to various stressors. Autophagy not only plays a crucial role in maintaining cellular homeostasis, but it is also essential to immunity, targets intracellular pathogens for degradation, modulates inflammation, and participates in adaptive immune responses. However, the expression profile of lncRNA and its role in regulating autophagy in macrophages have been poorly defined. Here, we used transcriptomic and bioinformatics to analysis LncRNA expression profile during autophagy and functional studies to evaluate the function of the metastasis-associated lung adenocarcinoma transcript-1 (Malat1) lncRNA in macrophages. A total of 1112 putative lncRNAs (240 novel lncRNAs) were identified, including 831 large intergenic, 129 intronic, and 152 anti-sense lncRNA, of which 59 differentially expressed transcripts exhibited a greater than 1.5-fold change under different conditions. The interaction of Malat1 lncRNA with microRNA (mir)-23-3p and lysosomal-associated membrane protein 1 (Lamp1) was found, Malat1 releases inhibition of Lamp1 expression in macrophages through competitive adsorption of mir-23-3p. The results of this study provide a better understanding of lncRNA function in macrophages and a basis for further investigation into the roles and mechanisms of ncRNA in immunology, particularly the functions of Malat1 and mir-23-3p in the pathogenesis of macrophages.