Transcriptional and Bioinformatic Analysis Provide a Relationship between Host Response Changes to Marek's Disease Viruses Infection and an Integrated Long Terminal Repeat
Ning Cui
Xianyao Li1 
Shuai Su- 1Department of Preventive Veterinary Medicine, College of Veterinary Medicine, Shandong Agricultural University, Tai'an, China
- 2Shandong Provincial Key Laboratory of Animal Biotechnology and Disease Control and Prevention, Shandong Agricultural University, Tai'an, China
- 3Department of Animal Nutrition and Feed Science, College of Animal Science, South China Agricultural University, Guangzhou, China
- 4Qingdao Animal Husbandry and Veterinary Research Institute, Qingdao, China
GX0101, Marek's disease virus (MDV) strain with a long terminal repeat (LTR) insert of reticuloendotheliosis virus (REV), was isolated from CVI988/Rispens vaccinated birds showing tumors. We have constructed a LTR deleted strain GX0101ΔLTR in our previous study. To compare the host responses to GX0101 and GX0101ΔLTR, chicken embryo fibroblasts (CEF) cells were infected with two MDV strains and a gene-chip containing chicken genome was employed to examine gene transcription changes in host cells in the present study. Of the 42,368 chicken transcripts on the chip, there were 2199 genes that differentially expressed in CEF infected with GX0101 compared to GX0101ΔLTR significantly. Differentially expressed genes were distributed to 25 possible gene networks according to their intermolecular connections and were annotated to 56 pathways. The insertion of REV LTR showed the greatest influence on cancer formation and metastasis, followed with immune changes, atherosclerosis, and nervous system disorders in MDV-infected CEF cells. Based on these bio functions, GX0101 infection was predicated with a greater growth and survival inhibition but lower oncogenicity in chickens than GX0101ΔLTR, at least in the acute phase of infection. In summary, the insertion of REV LTR altered the expression of host genes in response to MDV infection, possibly resulting in novel phenotypic properties in chickens. Our study has provided the evidence of retroviral insertional changes of host responses to herpesvirus infection for the first time, which will promote to elucidation of the possible relationship between the LTR insertion and the observed phenotypes.
Introduction
Marek's disease (MD) is caused by MDV and characterized by bursal-thymic atrophy with several distinct pathologic syndromes that include early mortality syndrome, cytolytic infection, immunodepression, transient paralysis, persistent neurological diseases, atherosclerosis, local lesions, and transplants in chickens (Calnek, 2001). It has long been recognized as a model for the evolution of pathogen virulence derived by vaccine (Gandon et al., 2001; Witter, 2001).
Vaccination for MD is an outstanding example of successful disease control in the area of veterinary medicine. However, frequent outbreaks of MD occur on individual farms or regions in recent years. Several reports lend credence to the implication of exceptionally virulent MDV isolates in vaccine failures (Witter et al., 1980; Powell and Lombardini, 1986; Suresh et al., 2015; Cui et al., 2016). The MDV field strain GX0101 with a long terminal repeat (LTR) region of reticuloendotheliosis virus (REV) insert was isolated from CVI988/Rispens vaccinated birds showing tumors in 2005 (Zhang and Cui, 2005). The pathogenesis of GX0101 was higher than virulent MDV (vMDV) GA strain but lower than very virulent MDV (vvMDV) strain Md5. The outbreak of MD caused by GX0101 was not attributed to a greater virulence but might associated with higher transmission capacity (Sun A. J. et al., 2010). It has also been demonstrated that the REV LTR could be integrated into the MDV genome following both long- and short-term co-infections of chicken or duck embryo fibroblasts (DEF) (Isfort et al., 1992). RM1 with a REV LTR insert was attenuated for oncogenicity but not for its immunosuppressive effect or in vivo replication (Witter et al., 1997). Several possible consequences of the REV LTR insertion into herpesviruses, including the transmission of retroviral information by herpesviruses, the activation, or inactivation of herpesvirus genes, the alteration of herpesvirus biological properties, etc., have also been discussed (Isfort et al., 1994; Jones et al., 1996; Witter et al., 1997). But the mechanism for the differential pathogenesis caused by the herpesviruses with/without a REV LTR is complex and not fully understood.
With the completion of the chicken genome sequencing and annotation endeavor and the availability of upgrade chicken-specific microarrays, gene expression profiling has become a valuable tool for evaluating host–pathogen interactions and was extensively applied to study the host response to MDV. Global gene expression profiling has been used to identify genes that are differentially expressed in response to MDV infection (Chen et al., 2011), vaccination with MD vaccines (Morgan et al., 2001), and between MD-resistant and -susceptible chicken lines (Sarson et al., 2008). Numerous genes involved in host response to MDV infection were defined in the microarray studies, which take the place to explore the mechanism involving MDV pathogenesis from the aspect of tumorigenesis, immunity, and host susceptibility. While this article has focused on the influence of insertion of REV LTR on global gene expression profiling of MDV-infected cells. The gene transcription profile of CEF cells infected with different MDV strains (with/without a REV LTR) reported here provides a unique opportunity to identify possible genes involved in the mechanism of the differential virulence of MDVs. Most significantly, this study provides an initial lead to explore the possible mechanism of emerging MDVs from vaccinated chicken flocks.
Materials and Methods
Cell Culture and MDV Infection
Ten-day-old specific-pathogen free white leghorn chicken embryos for preparation of chicken embryo fibroblast (CEF) cultures were from SPAFAS Co. (Jinan, China). Primary CEF cells collected from one embryo were seeded onto three individual flasks with the density of 5 × 106 cells/flask. Cells in one of the flasks were infected by GX0101 (1.5 × 105 PFU/flask), the other flask infected by GX0101 ΔLTR (1.5 × 105 PFU/flask), and the last one was mock infected as control. In total, four chicken embryos as experimental replicates were used in the current study.
Total RNA Extraction
Total RNA was extracted from GX0101-infected and GX0101ΔLTR-infected CEF cells at 56 h post infection by using TRIzol Reagent (Life Technologies, Carlsbad, CA, US) following the manufacturer's instructions. RNA integrity was checked on an Agilent Bioanalyzer 2100 (Agilent Technologies, Santa Clara, CA, USA). The total RNA was then purified using an RNeasy Mini Kit (QIAGEN, GmBH, Germany) and an RNase-Free DNase Set (QIAGEN, GmBH, Germany). The RNA quality was assessed by formaldehyde agarose gel electrophoresis and was quantified spectrophotometrically. High-quality RNA with an A260: A280 ratio between 1.8 and 2.0, and intact ribosomal 28S and 18S bands were used for microarray and real-time-qPCR (RT-qPCR) analysis.
Probe Labeling and Microarray Hybridization
A 4 × 44 K Agilent custom chicken oligo microarray (array ID: 042688) was used to compare transcription profile of GX0101-infected and GX0101ΔLTR-infected CEF cells. Four biological replicates were used in each group with dye balance. Fluorescently labeled complementary RNA (cRNA) probes were generated by using the Two Color Microarray Quick Labeling kit (Agilent Technologies, Santa Clara, CA, USA) and following the manufacturer's instructions. Labeled cRNA were purified by RNeasy mini kit (QIAGEN, GmBH, Germany). Each slide on the 4 × 44 K microarray was hybridized with 825 ng Cy3-labeled cRNA and 825 ng Cy5-labeled cRNA using Gene Expression Hybridization Kit (Agilent Technologies, Santa Clara, CA, USA) in Hybridization Oven (Agilent Technologies, Santa Clara, CA, USA), according to the manufacturer's instructions. After 17 h hybridization, slides were washed in staining dishes (Thermo Shandon, Waltham, MA, USA) with Gene Expression wash Buffer Kit (Agilent Technologies, Santa Clara, CA, USA). The labeling, hybridization and washing procedures were followed according to Agilent's recommendation and described in detail previously (Chiang et al., 2008). Slides were scanned by Agilent Microarry Scanner (Cat#G2565CA, Agilent technologies, Santa Clara, CA, US) and Feature Extraction software 10.7 (Agilent technologies, Santa Clara, CA, US).
Microarray Data Normalization and Analysis
Signal intensity of each probe was filtered against negative controls before normalization. Comparisons were made between CEF cells infected by MDV GX0101 and CEF cells infected by GX0101ΔLTR. Data normalization was performed using locally weighted scatter plot smoothing (LOWESS) by R project (http://www.r-project.org). The normalized natural log intensities were analyzed by SAS using a mixed model (SAS, Cary, NC) with fixed effect of infection with MDV (GX0101 or GX0101ΔLTR) and dye (Cy5 or Cy3), interaction between CEF cells and treatment, and random effect of slide and array. P-value and relative fold changes between each comparison for each gene were calculated. After analysis using the method, a gene was considered to be significantly differentially expressed only if the log2 median of the ratios of the Cy5: Cy3 signal was greater than 1.00-fold or lower than −1.00-fold with P ≤ 0.05.
Ingenuity Pathway Analysis
Functional interpretation of significantly differentially expressed genes were analyzed in the context of gene ontology and molecular networks by using Ingenuity Pathway Analysis (IPA) software (Ingenuity Systems; www.ingenuity.com). Those genes with known gene-probe ID numbers and corresponding expression fold-changes were uploaded into the software. In IPA, the analysis was done with P ≤ 0.05 as the cut-off point. The IPA analysis assigned those genes into subcategories within each category, and supplied an appropriate P-value and the number of genes identified. The differentially expressed genes were also mapped to genetic networks in the IPA database and ranked by scores that denoted the probability that a collection of genes equal to or greater than the number in a network could be achieved by chance alone. IPA uses a right-tailed Fisher's exact test to calculate the P-value for functional categories, networks, and canonical pathway analysis.
Confirmation of Microarray Results by RT-qPCR
Validation of differential gene expression was performed for a number of genes that were found to be differentially expressed in the microarray analysis by RT-qPCR. The primers (Table 1) were designed by using the PRIMER3 program and were based on published target sequences. RT-qPCR was performed using a 7500 System (ABI, Singapore) with SYBR®; Premix Ex Taq™II (TaKaRa, China). The temperature program used for the amplification was 95°C for 30 s and 40 cycles of 95°C for 5 s, and 60°C for 34 s. Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) was used as the endogenous reference gene. The relative fold change of the differentially expressed genes was calculated through the 2−ΔΔCq method (Livak and Schmittgen, 2001). Triplicate RT-qPCRs were performed on each cDNA to guarantee the reproducibility of the amplification.
Table 1. Validation of microarray data by real-time RT-qPCR.
Results
Gene Transcription Changes in Host Cells Infected with Different MDVs
In order to characterize the relationship between host response changes to MDV infection and an integrated LTR further, microarray analysis with GX0101-infected, and GX0101ΔLTR-infected CEF cells were performed. The infected cells for microarray analysis were collected at 56 h post-infection, based on the transform status (defined MDV specific plaques). Transcription level of 2199 normalized genes were significantly altered in the GX0101-infected cells during the transformation stage, as compared to the GX0101ΔLTR-infected cells (shown in Supplementary Table). The insertion of LTR caused 720 up-regulated genes and 1479 down-regulated genes in MDV-infected cells.
Validation of Differential Transcripts
RT-qPCR was used to verify a subset of 15 differential expressed genes between GX0101-infected and GX0101ΔLTR-infected cells during the transform stage of MDV infection. We found that the RT-qPCR results of the 15 selective genes were consistent with the microarray analysis (Table 1).
Annotation of the Differentially Expressed Genes
A subset of 460 differentially expressed genes have been characterized with specific bio functions in IPA. The bio function of other genes remain to be determined. The 460 differentially expressed genes were distributed to three categories including molecular and cellular functions, disorders and diseases, physiological responses, and physiological system development. It can be seen that the top molecular and cellular functions were cellular movement, cellular development, cellular growth and proliferation, cellular assembly and organization, cellular function and maintenance, cell death and survival, cell morphology ranked according to –log (P-value) (Figure 1). Based on the activation z-score, top bio functions with a predicated increased activation state were organismal death, Hypoplasia, morphology of embryonic tissue, growth failure, while top bio functions with a decreased activation state were kidney development, quantity of gonadal cells, size of body, cell survival (Figure 2).
Figure 1. Functional gene ontology for differentially expressed genes. The top 16 functional groups, with larger numbers of genes, were categorized by using the IPA program. *indicated the top molecular and cellular functions.
Figure 2. Predicted activation state of 460 differentially expressed genes analyzed with the IPA program. Important focus molecules belonging to different category were displayed according to their function annotation. Predicted activation state (increased/decreased) are presented according to p-value.
Gene Network Analysis
To further refine the analysis of the 2199 genes, we investigated their biological interactions by using the IPA program, and found 460 differential expressed genes with functional relationships were mapped to interactive gene networks. Insertion of REV LTR caused 25 possible gene networks in MDV-infected CEF cells (P ≤ 0.05). The top 10 gene networks were listed in Table 2. Gene networks with high-scoring are showed in Figure 3 and are associated with hematological disease, metabolic disease, connective tissue disorders, cell-to-cell signaling and interaction, nervous system development and function, tissue development, cardiovascular disease, embryonic development, organismal development, cellular movement, cellular function and maintenance, cellular growth, and proliferation. The remain top gene networks are involved in bio functions related to humoral immune response, protein synthesis, amino acid metabolism, ophthalmic disease, renal and urological system development and function, dermatological diseases and conditions, gastrointestinal disease, molecular transport, small molecule biochemistry, connective tissue development and function, digestive system development and function (Table 2).
Table 2. The highest scoring networks (P ≤ 0.05).
Figure 3. Significant four-gene networks of 460 differentially expressed genes analyzed with the IPA program. Molecular interactions among important focus molecules are displayed. Node (gene) and edge (gene relationship) symbols are described in the key. The intensity of the node color indicates the degree of up-(red) or down- (green) regulation; The color-change scale bar is a log2 scale. Symbols for each molecule are presented according to their molecular function and the type of interactions they participate in.
Canonical Pathway Analysis
As analyzed by the IPA canonical pathway analysis program, 460 differential expressed genes were annotated to 56 metabolic and signaling pathways in the IPA canonical pathway library (Table 3, P ≤ 0.05). More than 20 pathways directly related to tumors were modulated in GX0101-infected CEF cells compared to GX0101ΔLTR-infected ones. Four pathways directly related to atherosclerosis caused by MDV infection were regulated by the insertion of REV LTR: atherosclerosis signaling; cysteine biosynthesis III; extrinsic prothrombin activation pathway; cysteine biosynthesis/homocysteine degradation. A series of complex immune changes that related to crosstalk between dendritic cells and natural killer cells, communication between innate and adaptive immune cells and T helper cell differentiation were also found.
Table 3. List of the genes in significant canonical pathways (P ≤ 0.05).
Discussion
The insertion of REV LTR into MDV genome is associated with alteration of biological properties of MDVs, and may further provides a selective force for virus evolution (Isfort et al., 1994; Jones et al., 1996; Witter et al., 1997). Herein, we studied the host response changes to MDV infection by an integrated LTR in vitro using high-throughout microarray analysis. We found that some differentially expressed genes that have already been reported as such for the MD system were also significantly altered in our study (reviewed by Haq et al., 2010), but majority differentially expressed genes were reported for the first time.
The insertion of REV LTR changed the expressed transcription profiles of host cells after MDV infection that associated with cell-to-cell signaling and interaction appear to be a common strategy for the host responses to stress caused by virus infection. The oxidants generated in MD causes significant damage to DNA (Keles et al., 2010). In response to stress, HSP90 promotes the maturation, structural maintenance and proper regulation of specific target proteins involved for instance in cell cycle control and signal transduction as a molecular chaperone (Chen et al., 2015). While Hsp90 is an vital protein that mediates the infectivity and replication of some herpes-viruses, such as herpes simplex virus type 1 (Burch and Weller, 2005), Epstein-Barr virus (Sun X. et al., 2010) and Kaposi's sarcoma-associated herpesvirus (Wen and Damania, 2010), the transcription level is reduced in MDV-infected chickens (Hu et al., 2012). The inserted REV LTR activated transcription of HSP90 and might play a negative role in regulating cell cycles after MDV infection. A set of proteins that interact with HSP90 and involve in promoter responses to many activators and repressors were also modulated. CDC5L, DNA binding and poly(A) RNA binding related proteins, that take part in signal transduction involved in DNA damage checkpoint, were down-regulated; on the contrary, RFWD3 has a positive regulation role in defense response to virus by host regulation of DNA damage checkpoint, was up-regulated (Gong and Chen, 2011). Thus, the insertion of REV LTR regulated host cell cycles and DNA damage checkpoint in response to stress by MDV infection, both of which were closely associated with tumor formation. In fact, we found that more than 20 pathways directly related to tumors were modulated in GX0101-infected CEF cells compared to GX0101ΔLTR-infected CEF cells. From various cancer formation and metastasis signals, it can be seen that a set of matrix metalloproteases (MMPs) (MMP27, MMP23B, THBS2, MMP13, MMP11) as well as the MMPs receptor ITGB3 were inhibited in GX0101-infected CEFs. MMPs are members of an enzyme family that are critical for maintaining tissue allostasis and are required for tumor invasion and metastasis by destroy extracellular matrix (Clark et al., 2008). In line with this, oncostatin M signaling that inhibits the proliferation of some tumor cells (Krona et al., 2007) was activated. RAR signal that plays a role in anti-cancer aspect by suppress the proliferation of tumor cells and promote apoptosis was also activated (Brigger et al., 2015). Therefore, there is a lower chances for GX0101 with a REV LTR to induce tumors in vivo. This provides some clues and support the report that the insertion of REV LTR into MDV genome is closely related to the tumor formation in vivo (Jones et al., 1996; Witter et al., 1997). The REV LTR shows strong promoter or enhancer activity, and various genes are transactivated depending on the location of the insertion (Jones et al., 1993, 1996). The insertion of REV LTR into RM1 attenuate the oncogenicity of the virus in vivo replication (Witter et al., 1997). The field MDV strain GX0101 with a REV LTR, shows a lower but not significant oncogenicity in vivo (Sun A. J. et al., 2010). Our study have provided some evidence that the GX0101 with a REV LTR indeed showed a lower oncogenicity from the aspect of host gene expression profiles.
The growth inhibition and early mortality syndrome caused by MDV infection might be altered by the inserted REV LTR in MDV genome. Our results showed that REV LTR caused predicated increased activation state of organismal death and growth failure, while decreased activation state of size of body, cell survival. In line with this, Network 3 clustered a set of genes which had top functions associated embryonic development and organismal development. WNT signal was activated as indicated by the up-regulation of WBP2, WNT6, and WNT16 genes. WNT signal integrates signal (like retinoic acid, FGF, TGF-β, and BMP) delivered by other pathways and regulates the embryonic and organismal development. In addition, protein kinase A (Pka) signal pathway was supposed to be significantly up-regulated in the current gene expression mode. Pka system is one of the G protein-coupled system signal transduction pathways. GPCRs is activated by various external stimuli, which further activates the downstream MAPK cascade signaling as revealed by network 4. Activation of MAPK cascade is an essential event in various cellular activity including cellular growth and proliferation, and cellular movement and death, those are the top functions of network 4. In our study, we found that LTR up-regulated key activators in the pathway like IGFBP3, INHBA, and MAP2K3, which initiate downstream MAPK cascade signaling. IGFBP3 can positive regulate the insulin-like growth factor receptor and growth factor; INHBA takes part in activating the Erk5/BMK1 MAPK cascade signal leading to cell growth, differentiation and development; MAP2K3 can activate TAK1 that mediates p38 MAPK activation and plays a role in inflammation, apoptosis, growth and differentiation. Other differential expressed genes identified (Table 2) that associated to the cell differentiation and proliferation, and cellular function and maintenance were also up-regulated, for example, EGLN3, that responses to hypoxia and regulate cell proliferation; PHLDA2, that plays a key role in organ morphogenesis, regulation of gene expression, cell migration, and embryonic development; TP53I11, that responses to stress and regulate cell transcriptional activity and apoptosis. The combined results from the above studies indicated that the insertion of REV LTR into MDV GX0101 genome had a greater growth and survival inhibition on chickens caused by GX0101 infection. Previous in vivo study demonstrated that GX0101 with a REV LTR insert in the genome was attenuated for pathogenic effects with lower growth and mortality rates of the infected birds during 100 days post infection (Sun A. J. et al., 2010). The in vivo and in vitro experiments seem contradictory and this is probably because that the MDV infection cellular model may not represent the whole infection phase, while indicates the productive infection phase of MDV infection in vivo.
It has long been recognized that atherosclerosis and nervous system disorders are clinically distinct disease syndromes caused by MDV infection (Fabricant and Fabricant, 1999). Four pathways directly related to atherosclerosis caused by MDV infection were regulated by the insertion of REV LTR. PLA2G10, an indicator of atherosclerosis inflammatory (Piñón and Kaski, 2006), was significantly up-regulated in GX0101-infected cells. On the contrary, most of the genes in these pathways were down-regulated, such as CBS, CTH, F7, MGMT, MAT1A, COL1A2, ALB, CXCR4, APOA1, MMP13. It was reported that MAT1A and APOA1 were required to maintain normal cardiovascular system (Andrikoula and McDowell, 2008). Myosins, that plays a vital role in muscle contraction and involves in a wide range of other motility activities in eukaryotes, were predicted to be down-regulated by a set of calcium ion binding proteins and associated proteins including MYL1, TNNC1, TNNT2, MYH7B, TEAD1, and WNT11 (Leslie, 2015). A group of genes involved in F-actin protein were highly down-regulated (NEXN, TNNT3, TRIM55, TTN, MYOT). Actin family proteins have a close relationship with atherosclerosis pathogenesis (Berceli et al., 1991). Therefore, The insertion of REV LTR can alter the atherosclerosis status caused by MDV infection. Besides, four genes were clustered to play a protective role in neurodegenerative diseases, including two up-regulated genes, PRKACB and NTS, and two down-regulated genes, TAC1 and ACE. ID1, CHRNG, and S100B were the three genes that directly related to a predicted down-regulated nicotinic acetylcholine receptor that play important roles in the peripheral nervous system (Pillai et al., 2011). These genes might related to the modified nervous system disorders in MDV-infected chickens caused by the insertion of REV LTR into MDV genome; the functions of all remain to be investigated in the future.
The insertion of REV LTR is also associated with a series of complex immune changes that related to crosstalk between dendritic cells and natural killer cells, communication between innate and adaptive immune cells and T helper cell differentiation in MDV infections. The gene array showed increased transcription of IL-6, IL-8, IRF1, CD83, HLA-DRB5, IFNGR2, CCL20 while decreased transcription of IL-4, IL-15, TLR7, TNFSF10, IL13RA1. IL-6 is involved in the acute phase immune response and involved in lymphocyte and monocyte differentiation (Mayer et al., 1991). It plays an essential role when B-cells are finally differentiated into Ig-secreting cells. IL-4 induces the expression of MHC class II molecules on resting B-cells and activates the process of several immune cell activation including B-cells (Ryzhov et al., 2004). CD83 may be involved in the antigen presentation or the cellular interactions after the activation of lymphocyte and negative regulates IL-4 production (Roy et al., 2004). IL13RA1 is an alternate accessory protein to the gamma chain of the common cytokine receptor in IL-4 signal transduction pathway (Miloux et al., 1997). Taken together, the insertion of REV LTR into MDV genome might stimulate the differentiation of B-cells to Ig-secreting cells, while inhibited the B-cell activation process. Therefore, the host humoral immunity was dysregulated by the insertion of REV LTR into MDV genome, even though that no obvious changes of antibody titer was observed in vivo study (Sun A. J. et al., 2010). IL-15 that stimulates the proliferation of T-lymphocytes and mediate the crosstalk between dendritic cells and natural killer cells was down-regulated (Boudreau et al., 2011). TRL7 immediate innate immune response that plays a role in defense response to virus was also down regulated. This indicated that the REV LTR may help GX0101 to introduce a immune evasion mechanism to escape the host innate and adaptive immune responses. IL-8 is secreted by several immune cells in response to an inflammatory stimulus and serves as a chemotactic factor that attracts a set of cell types except monocytes. Neutrophils can be activated by IL-8. CCL20 is hemotactic factor that mainly attracts lymphocytes. The increased transcription of IL-8 and CCL20 indicated a up-regulated inflammatory status in GX0101-infected cells compared to GX0101ΔLTR. In addition, retinoate biosynthesis I was also influenced in MDV-infected cells by the insertion of REV LTR with the activation of retinoic acids receptors (RAR). Retinoic acids (RA) participates in apoptosis process of T lymphocytes through the interaction of RARs and RXRs (Szondy et al., 1998) and mediates the Th1 to Th2 polarization via RARs (Iwata et al., 2003). RA is closely related to thymocytes apoptosis. Thus, the insertion of REV LTR might relate to thymus atrophy in MDV infections by influencing of the RA signal. In addition, the microarray analysis of host cells in responses to different MDVs infection also detected a various of other pathways that take part in virus infection. For instance, clathrin-mediated endocytosis signaling that mediated virus endocytosis by the inward budding of plasma membrane vesicles was modified (Reis et al., 2015). Wnt/GSK-3β signaling that plays a role in the pathogenesis of influenza was also modulated (Hiyoshi et al., 2015). The precise role of the indicated pathways in the differential host responses after the insertion of REV LTR into MDV genome have not been verified, further investigations are being performed to identify unique and more deeply involved interactions between host cells and different MDVs.
Taken together, our study using microarray analysis of CEF cells infected with MDVs with/without REV LTR, provides novel insights on the mechanistic basis of how REV LTR alters the pathogenesis of the productive infection phase of MDV infection in vitro. In the study, a number of differential transcript genes were clustered to a large number of biological pathways according to their biological functions. The insertion of REV LTR showed the greatest influence on cancer formation and metastasis, followed with immune changes, atherosclerosis and nervous system disorders. Based on the bio functions, GX0101 infection was predicated with a greater growth and survival inhibition but lower oncogenicity on chickens than GX0101ΔLTR, at least in the acute phase of infection. More detailed investigation should be carried out to eliminate the possible relationship between these differential genes with the selective advantages for GX0101 to become the predominant strain that could be isolated at higher frequency.
Author Contributions
NC and HH collection and assembly of the data, manuscript writing, and data analysis; XL, NC, and SS discussion, manuscript revision; SS and ZC concept and design, data analysis, manuscript revision, and final approval of the manuscript.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
The work was supported by the National Natural Science Foundation of China under Grant (31402235).
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fcimb.2016.00046
Abbreviations
MD, Marek's disease; MDV, Marek's disease virus; REV, reticuloendotheliosis virus; LTR, long terminal repeat; CEF, chicken embryo fibroblasts; DEF, duck embryo fibroblasts; GAPDH, Glyceraldehyde-3-phosphate dehydrogenase.
References
Andrikoula, M., and McDowell, I. F. (2008). The contribution of ApoB and ApoA1 measurements to cardiovascular risk assessment. Diabetes Obes. Metab. 10, 271–278. doi: 10.1111/j.1463-1326.2007.00714.x
Berceli, S. A., Borovetz, H. S., Sheppeck, R. A., Moosa, H. H., Warty, V. S., Armany, M. A., et al. (1991). Mechanisms of vein graft atherosclerosis: LDL metabolism and endothelial actin reorganization. J. Vasc. Surg. 13, 336–347. doi: 10.1016/0741-5214(91)90227-L
Boudreau, J. E., Stephenson, K. B., Wang, F., Ashkar, A. A., Mossman, K. L., Lenz, L. L., et al. (2011). IL-15 and type I interferon are required for activation of tumoricidal NK cells by virus-infected dendritic cells. Cancer Res. 71, 2497–2506. doi: 10.1158/0008-5472.CAN-10-3025
Brigger, D., Schläfli, A. M., Garattini, E., and Tschan, M. P. (2015). Activation of RARα induces autophagy in SKBR3 breast cancer cells and depletion of key autophagy genes enhances ATRA toxicity. Cell Death Dis. 6:e1861. doi: 10.1038/cddis.2015.236
Burch, A. D., and Weller, S. K. (2005). Herpes simplex virus type 1 DNA polymerase requires the mammalian chaperone hsp90 for proper localization to the nucleus. J. Virol. 79, 10740–10749. doi: 10.1128/JVI.79.16.10740-10749.2005
Calnek, B. W. (2001). Pathogenesis of Marek's disease virus infection. Curr. Top. Microbiol. Immunol. 255, 25–55. doi: 10.1007/978-3-642-56863-3_2
Chiang, H. I., Swaggerty, C. L., Kogut, M. H., Dowd, S. E., Li, X., Pevzner, I. Y., et al. (2008). Gene expression profiling in chicken heterophils with Salmonella enteritidis stimulation using a chicken 44 K Agilent microarray. BMC Genomics 9:526. doi: 10.1186/1471-2164-9-526
Chen, C., Li, H., Xie, Q., Shang, H., Ji, J., Bai, S., et al. (2011). Transcriptional profiling of host gene expression in chicken liver tissues infected with oncogenic Marek's disease virus. J. Gen. Virol. 92, 2724–2733. doi: 10.1099/vir.0.034066-0
Chen, X., Zhao, C., Li, X., Wang, T., Li, Y., Cao, C., et al. (2015). Terazosin activates Pgk1 and Hsp90 to promote stress resistance. Nat. Chem. Biol. 11, 19–25. doi: 10.1038/nchembio.1657
Clark, I. M., Swingler, T. E., Sampieri, C. L., and Edwards, D. R. (2008). The regulation of matrix metalloproteinases and their inhibitors. Int. J. Biochem. Cell Biol. 40, 1362–1378. doi: 10.1016/j.biocel.2007.12.006
Cui, N., Su, S., Sun, P., Zhang, Y., Han, N., and Cui, Z. (2016). Isolation and pathogenic analysis of virulent Marek's disease virus field strain in China. Poult. Sci. pii:pew073. doi: 10.3382/ps/pew073
Fabricant, C. G., and Fabricant, J. (1999). Atherosclerosis induced by infection with Marek's disease herpesvirus in chickens. Am. Heart J. 138, 465–468. doi: 10.1016/S0002-8703(99)70276-0
Gandon, S., Mackinnon, M. J., Nee, S., and Read, A. F. (2001). Imperfect vaccines and the evolution of pathogen virulence. Nature 414, 751–756. doi: 10.1038/414751a
Gong, Z., and Chen, J. (2011). E3 ligase RFWD3 participates in replication checkpoint control. J. Biol. Chem. 286, 22308–22313. doi: 10.1074/jbc.M111.222869
Haq, K., Brisbin, J. T., Thanthrige-Don, N., Heidari, M., and Sharif, S. (2010). Transcriptome and proteome profiling of host responses to Marek's disease virus in chickens. Vet. Immunol. Immunopathol. 138, 292–302. doi: 10.1016/j.vetimm.2010.10.007
Hiyoshi, M., Indalao, I. L., Yano, M., Yamane, K., Takahashi, E., and Kido, H. (2015). Influenza A virus infection of vascular endothelial cells induces GSK-3β-mediated β-catenin degradation in adherens junctions, with a resultant increase in membrane permeability. Arch. Virol. 160, 225–234. doi: 10.1007/s00705-014-2270-5
Hu, X., Qin, A., Qian, K., Shao, H., Yu, C., Xu, W., et al. (2012). Analysis of protein expression profiles in the thymus of chickens infected with Marek's disease virus. Virol. J. 9:256. doi: 10.1186/1743-422X-9-256
Isfort, R., Jones, D., Kost, R., Witter, R., and Kung, H. J. (1992). Retrovirus insertion into herpesvirus in vitro and in vivo. Proc. Natl. Acad. Sci. U.S.A. 89, 991–995. doi: 10.1073/pnas.89.3.991
Isfort, R., Witter, R., and Kung, H. J. (1994). Retrovirus insertion into herpesvirus. Trends Microbiol. 2, 174–177. doi: 10.1016/0966-842X(94)90668-8
Iwata, M., Eshima, Y., and Kagechika, H. (2003). Retinoic acids exert direct effects on T cells to suppress Th1 development and enhance Th2 development via retinoic acid receptors. Int. Immunol. 15, 1017–1025. doi: 10.1093/intimm/dxg101
Jones, D., Brunovskis, P., Witter, R., and Kung, H. J. (1996). Retroviral insertional activation in a herpesvirus: transcriptional activation of US genes by an integrated long terminal repeat in a Marek's disease virus clone. J. Virol. 70, 2460–2467.
Jones, D., Isfort, R., Witter, R., Kost, R., and Kung, H. J. (1993). Retroviral insertions into a herpesvirus are clustered at the junctions of the short repeat and short unique sequences. Proc. Natl. Acad. Sci. U.S.A. 90, 3855–3859. doi: 10.1073/pnas.90.9.3855
Keles, H., Fidan, A. F., Cigerci, I. H., Kucukkurt, I., Karadas, E., and Dundar, Y. (2010). Increased DNA damage and oxidative stress in chickens with natural Marek's disease. Vet. Immunol. Immunopathol. 133, 51–58. doi: 10.1016/j.vetimm.2009.07.003
Krona, A., Åman, P., Örndal, C., and Josefsson, A. (2007). Oncostatin M-induced genes in human astrocytomas. Int. J. Oncol. 31, 1457–1463. doi: 10.3892/ijo.31.6.1457
Livak, K. J., and Schmittgen, T. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 25, 402–408. doi: 10.1006/meth.2001.1262
Mayer, P., Geissler, K., Valent, P., Ceska, M., Bettelheim, P., and Liehl, E. (1991). Recombinant human interleukin 6 is a potent inducer of the acute phase response and elevates the blood platelets in nonhuman primates. Exp. Hematol. 19, 688–696.
Miloux, B., Laurent, P., Bonnin, O., Lupker, J., Caput, D., Vita, N., et al. (1997). Cloning of the human IL-13R alpha1 chain and reconstitution with the IL4R alpha of a functional IL-4/IL-13 receptor complex. FEBS Lett. 401, 163–166. doi: 10.1016/S0014-5793(96)01462-7
Morgan, R., Sofer, L., Anderson, A., Bernberg, E., Cui, J., and Burnside, J. (2001). Induction of host gene expression following infection of chicken embryo fibroblasts with oncogenic Marek's disease virus. J. Virol. 75, 533–539. doi: 10.1128/JVI.75.1.533-539.2001
Piñón, P., and Kaski, J. C. (2006). Inflammation, atherosclerosis and cardiovascular disease risk: PAPP-A, Lp-PLA2 and cystatin C. New insights or redundant information? Rev. Esp. Cardiol. 59, 247–258. doi: 10.1016/S1885-5857(06)70028-7
Pillai, S., Rizwani, W., Li, X., Rawal, B., Nair, S., Schell, M. J., et al. (2011). ID1 facilitates the growth and metastasis of non-small cell lung cancer in response to nicotinic acetylcholine receptor and epidermal growth factor receptor signaling. Mol. Cell. Biol. 31, 3052–3067. doi: 10.1128/MCB.01311-10
Powell, P. C., and Lombardini, F. (1986). Isolation of very virulent pathotypes of Marek's disease virus from vaccinated chickens in Europe. Vet. Rec. 118, 688–691. doi: 10.1136/vr.118.25.688
Reis, C. R., Chen, P. H., Srinivasan, S., Aguet, F., Mettlen, M., and Schmid, S. L. (2015). Crosstalk between Akt/GSK3β signaling and dynamin-1 regulates clathrin-mediated endocytosis. EMBO J. 34, 2132–2146. doi: 10.15252/embj.201591518
Roy, K. C., Bandyopadhyay, G., Rakshit, S., Ray, M., and Bandyopadhyay, S. (2004). IL-4 alone without the involvement of GM-CSF transforms human peripheral blood monocytes to a CD1a(dim), CD83(+) myeloid dendritic cell subset. J. Cell Sci. 117, 3435–3445. doi: 10.1242/jcs.01162
Ryzhov, S., Goldstein, A. E., Matafonov, A., Zeng, D., Biaggioni, I., and Feoktistov, I. (2004). Adenosine-activated mast cells induce IgE synthesis by B lymphocytes: an A2B-mediated process involving Th2 cytokines IL-4 and IL-13 with implications for asthma. J. Immunol. 172, 7726–7733. doi: 10.4049/jimmunol.172.12.7726
Sarson, A. J., Parvizi, P., Lepp, D., Quinton, M., and Sharif, S. (2008). Transcriptional analysis of host responses to Marek's disease virus infection in genetically resistant and susceptible chickens. Anim. Genet. 39, 232–240. doi: 10.1111/j.1365-2052.2008.01710.x
Sun, A. J., Xu, X. Y., Petherbridge, L., Zhao, Y. G., Nair, V., and Cui, Z. Z. (2010). Functional evaluation of the role of reticuloendotheliosis virus long terminal repeat (LTR) integrated into the genome of a field strain of Marek's disease virus. Virology 397, 270–276. doi: 10.1016/j.virol.2009.11.017
Sun, X., Barlow, E. A., Ma, S., Hagemeier, S. R., Duellman, S. J., Burgess, R. R., et al. (2010). Hsp90 inhibitors block outgrowth of EBV-infected malignant cells in vitro and in vivo through an EBNA1-dependent mechanism. Proc. Natl. Acad. Sci. U.S.A. 107, 3146–3151. doi: 10.1073/pnas.0910717107
Suresh, P., Johnson Rajeswar, J., Sukumar, K., Harikrishnan, T. J., and Srinivasan, P. (2015). Pathotyping of recent Indian field isolates of Marek's disease virus serotype 1. Acta Virol. 59, 156–165. doi: 10.4149/av_2015_02_156
Szondy, Z., Reichert, U., and Fésüs, L. (1998). Retinoic acids regulate apoptosis of T lymphocytes through an interplay between RAR and RXR receptors. Cell Death Differ. 5, 4–10. doi: 10.1038/sj.cdd.4400313
Wen, K. W., and Damania, B. (2010). Hsp90 and Hsp40/Erdj3 are required for the expression and anti-apoptotic function of KSHV K1. Oncogene 29, 3532–3544. doi: 10.1038/onc.2010.124
Witter, R. L. (2001). Protective efficacy of Marek's disease vaccines. Curr. Top. Microbiol. Immunol. 255, 57–90. doi: 10.1007/978-3-642-56863-3_3
Witter, R. L., Li, D., Jones, D., Lee, L. F., and Kung, H. J. (1997). Retroviral insertional mutagenesis of a herpesvirus: a Marek's disease virus mutant attenuated for oncogenicity but not for immunosuppression or in vivo replication. Avian Dis. 41, 407–421. doi: 10.2307/1592197
Witter, R. L., Sharma, J. M., and Fadly, A. M. (1980). Pathogenicity of variant Marek's disease virus isolants in vaccinated and unvaccinated chickens. Avian Dis. 24, 210–232. doi: 10.2307/1589781
Keywords: Marek's disease virus, microarrays, bioinformatic analysis, long terminal repeat, virus evolution
Citation: Cui N, Li X, Chen C, Hao H, Su S and Cui Z (2016) Transcriptional and Bioinformatic Analysis Provide a Relationship between Host Response Changes to Marek's Disease Viruses Infection and an Integrated Long Terminal Repeat. Front. Cell. Infect. Microbiol. 6:46. doi: 10.3389/fcimb.2016.00046
Received: 14 March 2016; Accepted: 07 April 2016;
Published: 26 April 2016.
Edited by:
Yongqun He, University of Michigan, USAReviewed by:
Guanghua Fu, Fujian Academy of Agricultural Sciences, ChinaJiabo Ding, China Institute of Veterinary Drugs Control, China
Copyright © 2016 Cui, Li, Chen, Hao, Su and Cui. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Shuai Su, ssu6307@163.com;
Zhizhong Cui, zzcuisdau@126.com