Gene expression analysis of host spleen responses to Marek's disease virus infection at late tumor transformation phase

Serum Insulin-Like Growth Factor-Binding Protein 7 Deriving from Spleen and Lung Could Be Used for Early Recognition of Cardiac Surgery-Associated Acute Kidney Injury

INTRODUCTION

The utility of arithmetic product of urinary tissue metalloproteinase inhibitor 2 (TIMP2) and insulin-like growth factor-binding protein 7 (IGFBP7) concentrations has been widely accepted on early diagnosis of acute kidney injury (AKI). However, which organ is the main source of those two factors and how the concentration of IGFBP7 and TIMP2 changed in serum during AKI still remain to be defined.

METHODS

In mice, gene transcription and protein levels of IGFBP7/TIMP2 in the heart, liver, spleen, lung, and kidney were measured in both ischemia-reperfusion injury (IRI)- and cisplatin-induced AKI models. Serum IGFBP7 and TIMP2 levels were measured and compared in patients before cardiac surgery and at inclusion (0 h), 2 h, 6 h, and 12 h after intensive care unit (ICU) admission, and compared with serum creatinine (SCr), blood urea nitrogen (BUN), estimated glomerular filtration rate (eGFR), and serum uric acid (UA).

RESULTS

In mouse IRI-AKI model, compared with the sham group, the expression levels of IGFBP7 and TIMP2 did not change in the kidney, but significantly upregulated in the spleen and lung. Compared with patients who did not develop AKI, the concentration of serum IGFBP7 at as early as 2 h after ICU admission (sIGFBP7-2 h) was significantly higher in patients who developed AKI. The relationships between sIGFBP7-2 h in AKI patients and log2 (SCr), log2 (BUN), log2 (eGFR), and log2 (UA) were statistically significant. The diagnostic performance of sIGFBP7-2 h measured by the macro-averaged area under the receiver operating characteristic curve was 0.948 (95% CI, 0.853-1.000; p < 0.001).

CONCLUSION

The spleen and lung might be the main source of serum IGFBP7 and TIMP2 during AKI. The serum IGFBP7 value demonstrated good predictive accuracy for AKI following cardiac surgery within 2 h after ICU admission.

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.

Exploration of Alternative Splicing (AS) Events in MDV-Infected Chicken Spleens

Marek’s disease (MD) was an immunosuppression disease induced by Marek’s disease virus (MDV). MD caused huge economic loss to the global poultry industry, but it also provided an ideal model for studying diseases induced by the oncogenic virus. Alternative splicing (AS) simultaneously produced different isoform transcripts, which are involved in various diseases and individual development. To investigate AS events in MD, RNA-Seq was performed in tumorous spleens (TS), spleens from the survivors (SS) without any lesion after MDV infection, and non-infected chicken spleens (NS). In this study, 32,703 and 25,217 AS events were identified in TS and SS groups with NS group as the control group, and 1198, 1204, and 348 differently expressed (DE) AS events (p-value < 0.05 and FDR < 0.05) were identified in TS vs. NS, TS vs. SS, SS vs. NS, respectively. Additionally, Function enrichment analysis showed that ubiquitin-mediated proteolysis, p53 signaling pathway, and phosphatidylinositol signaling system were significantly enriched (p-value < 0.05). Small structural variations including SNP and indel were analyzed based on RNA-Seq data, and it showed that the TS group possessed more variants on the splice site region than those in SS and NS groups, which might cause more AS events in the TS group. Combined with previous circRNA data, we found that 287 genes could produce both circular and linear RNAs, which suggested these genes were more active in MD lymphoma transformation. This study has expanded the understanding of the MDV infection process and provided new insights for further analysis of resistance/susceptibility mechanisms.

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.

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.

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.

Imaging Mass Spectrometry and Proteome Analysis of Marek's Disease Virus-Induced Tumors

The highly oncogenic alphaherpesvirus Marek's disease virus (MDV) causes immense economic losses in the poultry industry. MDV induces a variety of symptoms in infected chickens, including neurological disorders and immunosuppression. Most notably, MDV induces transformation of lymphocytes, leading to T cell lymphomas in visceral organs with a mortality of up to 100%. While several factors involved in MDV tumorigenesis have been identified, the transformation process and tumor composition remain poorly understood. Here we developed an imaging mass spectrometry (IMS) approach that allows sensitive visualization of MDV-induced lymphoma with a specific mass profile and precise differentiation from the surrounding tissue. To identify potential tumor markers in tumors derived from a very virulent wild-type virus and a telomerase RNA-deficient mutant, we performed laser capture microdissection (LCM) and thereby obtained tumor samples with no or minimal contamination from surrounding nontumor tissue. The proteomes of the LCM samples were subsequently analyzed by quantitative mass spectrometry based on stable isotope labeling. Several proteins, like interferon gamma-inducible protein 30 and a 70-kDa heat shock protein, were identified that are differentially expressed in tumor tissue compared to surrounding tissue and naive T cells. Taken together, our results demonstrate for the first time that MDV-induced tumors can be visualized using IMS, and we identified potential MDV tumor markers by analyzing the proteomes of virus-induced tumors.IMPORTANCE Marek's disease virus (MDV) is an oncogenic alphaherpesvirus that infects chickens and causes the most frequent clinically diagnosed cancer in the animal kingdom. Not only is MDV an important pathogen that threatens the poultry industry but it is also used as a natural virus-host model for herpesvirus-induced tumor formation. In order to visualize MDV-induced lymphoma and to identify potential biomarkers in an unbiased approach, we performed imaging mass spectrometry (IMS) and noncontact laser capture microdissection. This study provides a first description of the visualization of MDV-induced tumors by IMS that could be applied also for diagnostic purposes. In addition, we identified and validated potential biomarkers for MDV-induced tumors that could provide the basis for future research on pathogenesis and tumorigenesis of this malignancy.

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.

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.

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.

The inhibitory effects of gga-miR-199-3p, gga-miR-140-3p, and gga-miR-221-5p in Marek's disease tumorigenesis

Marek's disease (MD) is lymphoproliferative neoplastic disease in chickens, which is caused by Marek's disease virus (MDV). Our previous study profiled microRNA (miRNA) transcriptome in MD lymphoma, and found that gga-miR-199-3p, gga-miR-140-3p, and gga-miR-221-5p were down-regulated in MD lymphoma. In this study, we further investigated their differential expression between MDV-infected spleens and noninfected spleens at 4, 7, 14, 21, and 28 d postinfection (dpi) to elucidate whether deregulation of them was specific to late tumor transformation phase or not. The results showed that gga-miR-199-3p was down-regulated at 14 and 28 dpi, and the expression of gga-miR-140-3p was decreased at 14 dpi, which indicated that deregulation of these miRNAs appeared since early stage of MD tumor transformation. Additionally, the inhibitory effects of gga-miR-199-3p, gga-miR-140-3p, and gga-miR-221-5p on MDV-transformed lymphoid cell line (MSB1) cells proliferation were observed, which suggested that these miRNAs acted as MD tumor suppressors. Their aberrant expression at early tumor transformation phase and suppressive role in cell proliferation indicated that they were involved in MD lymphoma transformation, and might play crucial roles in MD tumorigenesis.

Breed-specific transcriptome response of spleen from six to eight week old piglet after infection with Streptococcus suis type 2

Different pig breeds have shown differential susceptibility to the pathogen infection; however, molecular mechanisms of the infection susceptibility are not fully understood. Streptococcus suis type 2 (SS2) is an important zoonotic pathogen. To identify the genes responsible for infection susceptibility, pigs from two different breeds (Enshi black and Landrace) were inoculated with SS2 and their spleen transcriptome profiles were investigated in the present study. The differentially expressed genes (DEGs) were analyzed from infected versus control pigs in each breed, and then compared between both pig breeds. Enshi black pig showed more DEGs than Landrace (830 vs. 611) and most of these were due to down-regulated genes (543 vs. 387). However some DEGs were uniquely expressed in one breed, some were expressed in opposite direction in both breeds. A number of candidate genes and pathways are identified which might be involved in susceptibility to SS2, for example, MMP9 and Resistin were only significantly expressed in Landrace. NPG3 and PMAP23 were up-regulated in Landrace whereas down-regulated in Enshi black. LENG8 in control Landrace have inherently higher expression than control Enshi black. IGKV6 is down-regulated in Landrace but up-regulated in Enshi black. Overall, the transcriptome profiles are consistent with the clinical signs, i.e. the Enshi black is more susceptible to SS2 infection than Landrace. This is the first study to identify differential gene expression between indigenous and modern commercial pigs after in vivo SS2 infection using RNA-seq. The significant DEGs in splenic profiles between two pig breeds suggested considerable involvement of genetic background in susceptibility to the SS2 infection in pigs.

gga-miR-26a targets NEK6 and suppresses Marek's disease lymphoma cell proliferation

MicroRNA (miRNA) are a class of highly conserved, small noncoding RNA that emerge as key posttranscriptional regulators in various neoplastic transformations. Our previous study profiling the miRNA transcriptome in Marek's disease virus (MDV)-induced lymphoma revealed many novel and differentially expressed miRNA, including gga-miR-26a, which was downregulated in MDV-infected spleens of chickens. In this study, differential expression of gga-miR-26a between MDV-infected and noninfected spleens at 4, 7, 14, 21, and 28 d postinfection was analyzed by real-time PCR. The results showed gga-miR-26a were downregulated in MDV-infected spleens at cytolytic infection, latency, and tumor transformation phases. Subsequent cell proliferation assay revealed cell viability was lower in gga-miR-26a mimic transfection group than that in negative controls. Target genes of gga-miR-26a were identified by luciferase reporter gene assay. The results showed significant interaction between gga-miR-26a and Never In Mitosis Gene A (NIMA)-related kinase 6 (NEK6) gene. Subsequent gain of function experiment and Western blot assay showed that mRNA and protein levels of NEK6 were downregulated after gga-miR-26 mimic was transfected into MDV-transformed lymphoid cell line (MSB-1), indicating that NEK6 was modulated by gga-miR-26a. The expression of NEK6 showed a higher trend in MDV-infected samples including tumorous spleen and MD lymphoma from liver than that in noninfected controls. The results suggested that gga-miR-26a inhibited MSB-1 cell proliferation. Gga-miR-26a and its direct target, NEK6, might play important roles in MDV infection.

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

MicroRNAs (miRNAs) are small RNA molecules that regulate gene expression. Emerging evidence suggests that differential miRNA expression is associated with viral infection and tumorigenesis. Recently discovered microRNAs in the Marek's disease virus (MDV) genome have been suggested to have regulatory roles during MDV oncogenesis. To gain more insight into the molecular mechanisms of the tumorigenesis of MDV, we used microarrays to screen host and viral miRNAs that were sensitive to infection by MDV. Microarray analysis showed significant differential expression of 79 miRNAs, which was confirmed by qRT-PCR analysis. These data suggest that differentially expressed miRNAs may have major roles in MDV-induced tumorigenesis. In addition, we found two clades of chicken miRNAs had increased expression in splenic tumors and non-tumorous spleen tissues from GA-infected chickens. Thus, the expression of these miRNAs can be considered signatures for MDV infection and tumorigenesis.

Differential expression of Toll-like receptor genes in lymphoid tissues between Marek's disease virus-infected and noninfected chickens

Toll-like receptors (TLR) are trans-membrane sensors recognizing invading microbes. Toll-like receptors play a central role in initiating immune responses against several pathogens. In this study, we investigated the response of TLR and downstream genes to Marek's disease virus (MDV) infection. Forty 1-d-old chicks were randomly divided into 2 groups, with 20 chicks infected with MDV and 20 chicks mock-infected. Four chickens were euthanized respectively from infected and age-matched noninfected groups at 4, 7, 14, 21, and 28 d postinfection (dpi). Bursas, spleens, and thymuses were removed. The differential expression of TLR genes, including TLR3, TLR5, TLR7, TLR15, and TLR21, and downstream genes of TLR7, including MyD88, TRAF3, TRAF6, IFNA, IFNB, and IL6, in lymphoid tissues of MDV-infected and noninfected chickens was determined by real-time PCR. The results showed that the change of TLR genes was different in 3 lymphoid tissues. Expression of TLR7 and MyD88 was upregulated at 14 dpi and downregulated at 28 dpi in MDV-infected compared with noninfected spleens. The TRAF6 and IFNB were upregulated, and TRAF3, IFNA, and IL6 genes showed increasing trends in MDV-infected compared with noninfected spleens at 14 dpi. The expression of TLR3 and TLR15 genes was downregulated in MDV-infected compared with noninfected spleens at 28 dpi. The results indicated that TLR7 and its downstream genes were a response to MDV infection at 14 dpi. However, the function of TLR was impaired when the infection entered the tumor transformation phase. In bursas, TLR3 and TLR15 genes were upregulated at 7 and 4 dpi, respectively. It indicated that TLR3 and TLR15 might be involved in response to MDV infection in bursa at early phases. However, no differential expression of TLR genes was observed between MDV-infected and noninfected thymuses, which indicated that the thymus had little response to MDV infection mediated by TLR.

A systematic analysis of miRNA transcriptome in Marek's disease virus-induced lymphoma reveals novel and differentially expressed miRNAs

Marek's disease is a lymphoproliferative neoplastic disease of the chicken, which poses a serious threat to poultry health. Marek's disease virus (MDV)-induced T-cell lymphoma is also an excellent biomedical model for neoplasia research. Recently, miRNAs have been demonstrated to play crucial roles in mediating neoplastic transformation. To investigate host miRNA expression profiles in the tumor transformation phase of MDV infection, we performed deep sequencing in two MDV-infected samples (tumorous spleen and MD lymphoma from liver), and two non-infected controls (non-infected spleen and lymphocytes). In total, 187 and 16 known miRNAs were identified in chicken and MDV, respectively, and 17 novel chicken miRNAs were further confirmed by qPCR. We identified 28 down-regulated miRNAs and 11 up-regulated miRNAs in MDV-infected samples by bioinformatic analysis. Of nine further tested by qPCR, seven were verified. The gga-miR-181a, gga-miR-26a, gga-miR-221, gga-miR-222, gga-miR-199*, and gga-miR-140* were down-regulated, and gga-miR-146c was up-regulated in MDV-infected tumorous spleens and MD lymphomas. In addition, 189 putative target genes for seven differentially expressed miRNAs were predicted. The luciferase reporter gene assay showed interactions of gga-miR-181a with MYBL1, gga-miR-181a with IGF2BP3, and gga-miR-26a with EIF3A. Differential expression of miRNAs and the predicted targets strongly suggest that they contribute to MDV-induced lymphomagenesis.