Distinct expression pattern of miRNAs in Marek's disease virus infected-chicken splenic tumors and non-tumorous spleen tissues

Comparative transcriptomics analysis identifies crucial genes and pathways during goose spleen development

As the largest peripheral lymphoid organ in poultry, the spleen plays an essential role in regulating the body's immune capacity. However, compared with chickens and ducks, information about the age- and breed-related changes in the goose spleen remains scarce. In this study, we systematically analyzed and compared the age-dependent changes in the morphological, histological, and transcriptomic characteristics between Landes goose (LG; Anser anser) and Sichuan White goose (SWG; Anser cygnoides). The results showed a gradual increase in the splenic weights for both LG and SWG until week 10, while their splenic organ indexes reached the peak at week 6. Meanwhile, the splenic histological indexes of both goose breeds continuously increased with age, reaching the highest levels at week 30. The red pulp (RP) area was significantly higher in SWG than in LG at week 0, while the splenic corpuscle (AL) diameter was significantly larger in LG than in SWG at week 30. At the transcriptomic level, a total of 1710 and 1266 differentially expressed genes (DEGs) between week 0 and week 30 were identified in spleens of LG and SWG, respectively. Meanwhile, a total of 911 and 808 DEGs in spleens between LG and SWG were identified at weeks 0 and 30, respectively. Both GO and KEGG enrichment analysis showed that the age-related DEGs of LG or SWG were dominantly enriched in the Cell cycle, TGF-beta signaling, and Wnt signaling pathways, while most of the breed-related DEGs were enriched in the Neuroactive ligand-receptor interaction, Cytokine-cytokine receptor interaction, ECM-receptor interaction, and metabolic pathways. Furthermore, through construction of protein-protein interaction networks using significant DEGs, it was inferred that three hub genes including BUB1, BUB1B, and TTK could play crucial roles in regulating age-dependent goose spleen development while GRIA2, GRIA4, and RYR2 could be crucial for the breed-specific goose spleen development. These data provide novel insights into the splenic developmental differences between Chinese and European domestic geese, and the identified crucial pathways and genes are helpful for a better understanding of the mechanisms regulating goose immune functions.

Characterization of miRNAs in Embryonic, Larval, and Adult Lumpfish Provides a Reference miRNAome for Cyclopterus lumpus

Simple Summary

Lumpfish (Cyclopterus lumpus) is an emergent aquaculture species, and its miRNA repertoire is still unknown. miRNAs are critical post-transcriptional modulators of teleost gene expression. Therefore, a lumpfish reference miRNAome was characterized by small RNA sequencing and miRDeep analysis of samples from different organs and developmental stages. The resulting miRNAome, an essential reference for future expression analyses, consists of 443 unique mature miRNAs from 391 conserved and eight novel miRNA genes. Enrichment of specific miRNAs in particular organs and developmental stages indicates that some conserved lumpfish miRNAs regulate organ and developmental stage-specific functions reported in other teleosts.

Abstract

MicroRNAs (miRNAs) are endogenous small RNA molecules involved in the post-transcriptional regulation of protein expression by binding to the mRNA of target genes. They are key regulators in teleost development, maintenance of tissue-specific functions, and immune responses. Lumpfish (Cyclopterus lumpus) is becoming an emergent aquaculture species as it has been utilized as a cleaner fish to biocontrol sea lice (e.g., Lepeophtheirus salmonis) infestation in the Atlantic Salmon (Salmo salar) aquaculture. The lumpfish miRNAs repertoire is unknown. This study identified and characterized miRNA encoding genes in lumpfish from three developmental stages (adult, embryos, and larvae). A total of 16 samples from six different adult lumpfish organs (spleen, liver, head kidney, brain, muscle, and gill), embryos, and larvae were individually small RNA sequenced. Altogether, 391 conserved miRNA precursor sequences (discovered in the majority of teleost fish species reported in miRbase), eight novel miRNA precursor sequences (so far only discovered in lumpfish), and 443 unique mature miRNAs were identified. Transcriptomics analysis suggested organ-specific and age-specific expression of miRNAs (e.g., miR-122-1-5p specific of the liver). Most of the miRNAs found in lumpfish are conserved in teleost and higher vertebrates, suggesting an essential and common role across teleost and higher vertebrates. This study is the first miRNA characterization of lumpfish that provides the reference miRNAome for future functional studies.

Exosomal miRNA profiling from H5N1 avian influenza virus-infected chickens

Exosomes are membrane vesicles containing proteins, lipids, DNA, mRNA, and micro RNA (miRNA). Exosomal miRNA from donor cells can regulate the gene expression of recipient cells. Here, Ri chickens were divided into resistant (Mx/A; BF2/B21) and susceptible (Mx/G; BF2/B13) trait by genotyping of Mx and BF2 genes. Then, Ri chickens were infected with H5N1, a highly pathogenic avian influenza virus (HPAIV). Exosomes were purified from blood serum of resistant chickens for small RNA sequencing. Sequencing data were analysed using FastQCv0.11.7, Cutadapt 1.16, miRBase v21, non-coding RNA database, RNAcentral 10.0, and miRDeep2. Differentially expressed miRNAs were determined using statistical methods, including fold-change, exactTest using edgeR, and hierarchical clustering. Target genes were predicted using miRDB. Gene ontology analysis was performed using gProfiler. Twenty miRNAs showed significantly different expression patterns between resistant control and infected chickens. Nine miRNAs were up-regulated and 11 miRNAs were down-regulated in the infected chickens compared with that in the control chickens. In target gene analysis, various immune-related genes, such as cytokines, chemokines, and signalling molecules, were detected. In particular, mitogen-activated protein kinase (MAPK) pathway molecules were highly controlled by differentially expressed miRNAs. The result of qRT-PCR for miRNAs was identical with sequencing data and miRNA expression level was higher in resistant than susceptible chickens. This study will help to better understand the host immune response, particularly exosomal miRNA expression against HPAIV H5N1 and could help to determine biomarkers for disease resistance.

Analysis of the microRNA expression profiles of chicken dendritic cells in response to H9N2 avian influenza virus infection

MicroRNA (miRNA) plays a key role in virus-host interactions. Here, we employed deep sequencing technology to determine cellular miRNA expression profiles in chicken dendritic cells infected with H9N2 avian influenza virus (AIV). A total of 66 known and 36 novel miRNAs were differently expressed upon H9N2 infection, including 72 up-regulated and 30 down-regulated miRNAs. Functional analysis showed that the predicted targets of these miRNAs were significantly enriched in several pathways including endocytosis, notch, lysosome, p53, RIG-I-like and NOD-like receptor signaling pathways. These data provide valuable information for further investigating the roles of miRNA in AIV pathogenesis and host defense response.

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.

Understanding the cross-talk between host and virus in poultry from the perspectives of microRNA

In poultry, viral infections (e.g., Marek's disease virus, avian leukosis virus, influenza A virus, and so on) can cause devastating mortality and economic losses. Because viruses are solely dependent on host cells to propagate, they alter the host intracellular microenvironment. Thus, understanding the virus-host interaction is important for antiviral immunity and drug development in the poultry industry. MicroRNAs are crucial posttranscriptional regulators of gene expression in a wide spectrum of biological processes, including viral infection. Recently, microRNAs have been identified as key players in virus-host interactions. In this review, we will discuss the intricacies involved in the virus-host cross-talk mediated by host and viral microRNAs in poultry (i.e., chicken and ducks), as well as recent trends and challenges in this field. These findings may provide some insights into the rapidly developing area of research regarding viral pathogenesis and antiviral immunity in poultry production.

Research Progress on MicroRNAs Involved in the Regulation of Chicken Diseases

MicroRNAs (miRNAs) are small, non-coding RNA molecules that inhibit protein translation from target mRNAs. Accumulating evidence suggests that miRNAs can regulate a broad range of biological pathways, including cell differentiation, apoptosis, and carcinogenesis. With the development of miRNAs, the investigation of miRNA functions has emerged as a hot research field. Due to the intensive farming in recent decades, chickens are easily influenced by various pathogen transmissions, and this has resulted in large economic losses. Recent reports have shown that miRNAs can play critical roles in the regulation of chicken diseases. Therefore, the aim of this review is to briefly discuss the current knowledge regarding the effects of miRNAs on chickens suffering from common viral diseases, mycoplasmosis, necrotic enteritis, and ovarian tumors. Additionally, the detailed targets of miRNAs and their possible functions are also summarized. This review intends to highlight the key role of miRNAs in regard to chickens and presents the possibility of improving chicken disease resistance through the regulation of miRNAs.

Integrated Analysis of miRNA and mRNA Expression Profiles in Spleen of Specific Pathogen-Free Chicken Infected with Avian Reticuloendotheliosis Virus Strain SNV

The Reticuloendotheliosis virus (REV) primarily causes avian severe immunosuppression, in addition to other symptoms, which include avian dwarfing syndrome and chronic tumors in lymphoid and other tissue. To date, REV’s molecular mechanisms leading to immunosuppression have not been fully elucidated. In the current study, we aimed to elucidate the role of microRNAs (miRNA) in regulating gene expression during REV infections. Therefore, we used a high-dose spleen necrosis virus (SNV) model of REV to inoculate one-day-old specific pathogen-free (SPF) chickens, thereby inducing congenital infections. We analyzed miRNA and mRNA expression profiles using Next Generation Sequencing (NGS) in a total of 19 spleen samples that were collected at 7, 14, and 21 days post infection (dpi). The results showed that 63 differentially expressed miRNAs (DEmiRNAs) (30 known miRNAs and 33 novel miRNAs) and 482 differentially expressed target genes (DETGs) were identified. Integration analysis identified 886 known miRNA–mRNA and 580 novel miRNA–mRNA interaction pairs, which involved miRNAs that were inversely correlated with the above DETGs. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis revealed that the DETGs were considerably enriched in the immune-relevant pathways category, such as immune system, cell growth and death, signaling molecules and interaction, signal transduction, etc. We further verified selected immune-relevant miRNA and their DETGs while using quantitative RT-PCR (qRT-PCR). Overall, our data revealed valuable immune-related miRNA–mRNA interaction information that occurred during REV infections, thereby broadening our understanding of the REV-induced immunosuppression.

Gga-miR-219b targeting BCL11B suppresses proliferation, migration and invasion of Marek's disease tumor cell MSB1

Marek’s disease (MD), caused by Marek’s disease virus (MDV), is a lymphotropic neoplastic disease. Previous miRNAome analysis showed gga-miR-219b was significantly downregulated in MDV-induced lymphoma, and one of its potential target genes, B-cell chronic lymphocytic /lymphoma 11B (BCL11B) was predicted. In this study, we further investigated the function of gga-miR-219b, and the gain/loss of function assay showed gga-miR-219b inhibited cell migration and reduced cell proliferation by promoting apoptosis not by cell cycle arrest. Gga-miR-219b also suppressed expression of two cell invasion-related genes MMP2 and MMP9. The results indicated suppressive effect of gga-miR-219b on MD tumorigenesis. The gene BCL11B was verified as a direct target gene of gga-miR-219b. RNA interference was performed to block BCL11B. As expected, the effects triggered by BCL11B downregulation were in accordance with that triggered by gga-miR-219b overexpression, suggesting that BCL11B was a stimulative regulator of MD transformation. Moreover, both gga-miR-219b and BCL11B influenced the expression of Meq gene, the most important oncogene in MDV. Additionally, gene expression level of anti-apoptotic genes BCL2 and BCL2L1 was downregulated and pro-apoptotic gene TNFSF10 was upregulated in MSB1 cells with gga-miR-219b overexpression or BCL11B knockdown, which suggested gga-miR-219b promoted cell apoptosis via regulating gene expression in the apoptosis pathways.

gga-miR-2127 downregulates the translation of chicken p53 and attenuates chp53-mediated innate immune response against IBDV infection

Infectious bursal disease (IBD) is characterized by the immune suppression of infected birds. The molecular mechanism by which IBD virus (IBDV) suppresses the host immune system remains to be elucidated. The tumor suppressor protein p53 can inhibit the replication of various viruses, but its effect on IBDV remains unknown. This study established an in vitro infection model based on DF-1 cells (chicken embryo fibroblast cell line) to investigate the antiviral effects of chicken p53 (chp53) on IBDV infection. The expression level and activity of chp53 remarkably increased in IBDV-infected DF-1 cells. The overexpression of chp53 inhibited IBDV replication and upregulated the expression of multiple chicken antiviral innate immunity genes (IPS-1, IRF3, PKR, OAS, and Mx), whereas the suppression of chp53 led to the opposite effect. This result indicates that chp53 activates the antiviral innate immune response of chickens to IBDV infection. Bioinformatics analysis and dual-luciferase reporter assay showed that gga-miR-2127 targeted the 3'UTR of chp53. qRT-PCR and western blot revealed that gga-miR-2127 overexpression in DF-1 cells not only downregulated the expression levels of chp53 and of the antiviral innate immunity genes in chickens but also promoted IBDV replication. Our results suggest that gga-miR-2127 downregulates chp53 mRNA translation by targeting its 3'UTR and attenuates chp53-mediated antiviral innate immune response against IBDV.

Enrichment and verification of differentially expressed miRNAs in bursa of Fabricius in two breeds of duck

Objective

The bursa of Fabricius (BF) is a central humoral immune organ belonging specifically to avians. Recent studies had suggested that miRNAs were active regulators involved in the immune processes. This study was to investigate the possible differences of the BF at miRNA level between two genetically disparate duck breeds.

Methods

Using Illumina next-generation sequencing, the miRNAs libraries of ducks were established.

Results

The results showed that there were 66 differentially expressed miRNAs and 28 novel miRNAs in bursa. A set of abundant miRNAs (i.e., let-7, miR-146a-5p, miR-21-5p, miR-17~92) which are involved in immunity and disease were detected and the predicted target genes of the novel miRNAs were associated with duck high anti-adversity ability. By gene ontology analysis and enriching KEGG pathway, the targets of differential expressed miRNAs were mainly involved in immunity and disease, supporting that there were differences in the BF immune functions between the two duck breeds. In addition, the metabolic pathway had the maximum enriched target genes and some enriched pathways that were related to cell cycle, protein synthesis, cell proliferation and apoptosis. It indicted that the difference of metabolism may be one of the reasons leading the immune difference between the BF of two duck breeds.

Conclusion

This data lists the main differences in the BF at miRNAs level between two genetically disparate duck breeds and lays a foundation to carry out molecular assisted breeding of poultry in the future.

Chicken gga-miR-130a targets HOXA3 and MDFIC and inhibits Marek's disease lymphoma cell proliferation and migration

Marek's disease (MD) is an infectious disease of chickens caused by MD virus (MDV), which is a herpesvirus that initiates tumor formation. Studies have indicated that microRNAs (miRNAs) are linked with the development of cancers or tumors. Previously, gga-miR-130a was discovered downregulated in MDV-infected tissues. Here, we aimed to explore the further function of gga-miR-130a in MD. The expression of gga-miR-130a in MDV-infected and uninfected spleens was detected by quantitative real-time PCR (qRT-PCR). Subsequently, proliferation and migration assays of MDV-transformed lymphoid cells (MSB1) were carried out by transfecting gga-miR-130a. The target genes of gga-miR-130a were predicted using TargetScan and miRDB and clustered through Gene Ontology analysis. The target genes were validated by western blot, qRT-PCR, and a dual luciferase reporter assay. Our results show that the expression of gga-miR-130a was reduced in MDV-infected spleens. Gga-miR-130a showed an inhibitory effect on MSB1 cell proliferation and migration. Two target genes, homeobox A3 (HOXA3) and MyoD family inhibitor domain containing (MDFIC), were predicted and clustered to cell proliferation. Results indicate that gga-miR-130a regulates HOXA3 and MDFIC at the protein level but not at the mRNA level. Moreover, the gga-miR-130a binding sites of two target genes have been confirmed. We conclude that gga-miR-130a can arrest MSB1 cell proliferation and migration, and target HOXA3 and MDFIC, which are both involved in the regulation of cell proliferation. Collectively, gga-miR-130a plays a critical role in the tumorigenesis associated with chicken MD.

Chicken gga-miR-103-3p Targets CCNE1 and TFDP2 and Inhibits MDCC-MSB1 Cell Migration

Marek’s disease (MD) is a highly contagious viral neoplastic disease caused by Marek’s disease virus (MDV), which can lead to huge economic losses in the poultry industry. Recently, microRNAs (miRNAs) have been found in various cancers and tumors. In recent years, 994 mature miRNAs have been identified through deep sequencing in chickens, but only a few miRNAs have been investigated further in terms of their function. Previously, gga-miR-103-3p was found downregulated in MDV-infected samples by using Solexa deep sequencing. In this study, we further verified the expression of gga-miR-103-3p among MDV-infected spleen, MD lymphoma from liver, noninfected spleen, and noninfected liver, by qPCR. The results showed that the expression of gga-miR-103-3p was decreased in MDV-infected tissues, which was consistent with our previous study. Furthermore, two target genes of gga-miR-103-3p, cyclin E1 (CCNE1) and transcription factor Dp-2 (E2F dimerization partner 2) (TFDP2), were predicted and validated by luciferase reporter assay, qPCR, and western blot analysis. The results suggested that CCNE1 and TFDP2 are direct targets of gga-miR-103-3p in chickens. Subsequent cell proliferation and migration assay showed that gga-miR-103-3p suppressed MDCC-MSB1 migration, but did not obviously modulate MDCC-MSB1 cell proliferation. In conclusion, gga-miR-103-3p targets the CCNE1 and TFDP2 genes, and suppresses cell migration, which indicates that it might play an important role in MD tumor transformation.

Chicken gga-miR-181a targets MYBL1 and shows an inhibitory effect on proliferation of Marek's disease virus-transformed lymphoid cell line

Marek's disease (MD), caused by Marek's disease virus (MDV), is a lymphoproliferative neoplastic disease of chickens and is characterized by MD lymphoma in multiple visceral organs of chicken. It causes great damage to poultry health. Recently, miRNA has been reported to be involved in Marek's disease lymphomagenesis. Our previous study showed that gga-miR-181a was downregulated in MDV-induced lymphoma, and its target gene, v-myb myeloblastosis viral oncogene homolog-like 1 (MYBL1), was predicted. In this study, the interaction between gga-miR-181a and MYBL1 was further verified by detecting protein expression levels of MYBL1 after transfecting miR-181a mimic into MD lymphoma cell line, MSB1. The result showed that protein level of MYBL1 was lower in gga-miR-181a mimic transfecting group than that in the negative control group at 96 h post transfection, which indicated that MYBL1 was a target gene of gga-miR-181a. Additionally, we found that the expression of MYBL1 was higher in MDV-infected samples than that in non-infected controls, which agreed with the proposition that miRNA showed a negatively correlated expression pattern with its target gene. We observed the inhibitory effect of gga-miR-181a on MSB1 cell proliferation. Collectively, the aberrant expression of gga-miR-181a and MYBL1 in MD lymphoma suggested that they might be involved in MD tumor transformation and played important roles.

Identification and Characterization of MicroRNAs in Snakehead Fish Cell Line upon Snakehead Fish Vesiculovirus Infection

MicroRNAs (miRNAs) play important roles in mediating multiple biological processes in eukaryotes and are being increasingly studied to evaluate their roles associated with cellular changes following viral infection. Snakehead fish Vesiculovirus (SHVV) has caused mass mortality in snakehead fish during the past few years. To identify specific miRNAs involved in SHVV infection, we performed microRNA deep sequencing on a snakehead fish cell line (SSN-1) with or without SHVV infection. A total of 205 known miRNAs were identified when they were aligned with the known zebrafish miRNAs, and nine novel miRNAs were identified using MiRDeep2 software. Eighteen and 143 of the 205 known miRNAs were differentially expressed at three and 24 h post-infection (poi), respectively. From the differentially-expressed miRNAs, five were randomly selected to validate their expression profiles using quantitative reverse transcription polymerase chain reaction (qRT-PCR), and their expression profiles were consistent with the microRNA sequencing results. In addition, the target gene prediction of the SHVV genome was performed for the differentially-expressed host miRNAs, and a total of 10 and 58 differentially-expressed miRNAs were predicted to bind to the SHVV genome at three and 24 h poi, respectively. The effects of three selected miRNAs (miR-130-5p, miR-214 and miR-216b) on SHVV multiplication were evaluated using their mimics and inhibitors via qRT-PCR and Western blotting. The results showed that all three miRNAs were able to inhibit the multiplication of SHVV; whereas the mechanisms underlying the SHVV multiplication inhibited by the specific miRNAs need to be further characterized in the future.

MicroRNA profile analysis of Epithelioma papulosum cyprini cell line before and after SVCV infection

MicroRNAs (miRNAs) play significant roles in regulating almost all of the biological processes in eukaryotes. An accumulating body of evidence shows that miRNAs are associated with cellular changes following viral infection. Spring viremia of carp virus (SVCV) is the pathogen of Spring viremia of carp (SVC), which results in heavy losses in the cultured common carp (Cyprinus carpio) industry in many countries. To study the involvement of miRNAs during SVCV infection, we adopted the Solexa sequencing technology to sequence small RNA libraries from the Epithelioma papulosum cyprini (EPC) cell line before and after infection with SVCV. In this study, a total of 161 conserved and 26 novel miRNAs were identified. Subsequently, the expression patterns of these miRNAs were compared between the uninfected (control library, M) and SVCV-infected (infection library, E) libraries. In addition, to verify the Solexa sequencing results, the expression patterns of 14 randomly selected miRNAs were validated by qRT-PCR. The targets of the significantly differentially expressed miRNAs were then predicted, and the miRNAs that could directly target the SVCV genome were also predicted. No miRNA encoded by SVCV itself was detected. To the best of our knowledge, this study presents the first miRNA profiling assessment in association with fish rhabdovirus infection, and the data presented lay a foundation for further investigations to determine the roles of miRNAs in regulating the molecular mechanism during SVCV infection.