Transcriptome Analysis of Genes Responding to Infection of Leghorn Male Hepatocellular Cells With Fowl Adenovirus Serotype 4

Introduction

Fowl adenovirus serotype 4 (FAdV-4) is a member of the Aviadenovirus genus, Adenoviridae family. FAdV-4 is a non-enveloped virus with a double-stranded DNA genome. It is ~43.7 kb in length and encodes 11 and 32 structural and non-structural proteins, respectively (1). FAdV-4 causes hydropericardium hepatitis syndrome (HHS), a severe disease of poultry worldwide, which is characterized by pericardial effusion and associated with high mortality rates (30 to 80%), resulting in major financial losses all over the world. The first outbreak of HHS was recorded in 1987 in Angara Goth, Pakistan. Since then, FAdV-4 has become prevalent in several other parts of the world (2–4). Prior to 2014, HHS was only sporadically reported in poultry in China. However, since 2015, a significant increase in the frequency of HHS outbreaks has been noted in various production areas of broilers in China, causing major economic burden (5–7).

Chicken has huge commercial value globally in the food industry (8). However, FAdV-4 infection is a significant threat to the chicken industry (9). Therefore, understanding the molecular response of FAdV-4 disease-resistance is vital for prevention of FAdV-4 infection. Important innate signaling pathways activated during the molecular response to infection may assist development of vaccines and anti-viral drugs (10–12). Leghorn male hepatocellular (LMH) cells belong to an immortalized chicken hepatic cell line and are the primary target cells of FAdV-4; LMH cells therefore play an important role in in vitro studies of host–FAdV-4 interactions (1, 13, 14). Previous research has attempted to elucidate the mechanism of FAdV-4–host interactions. These studies provide valuable and important mechanistic clues regarding virus–host interactions; however, our understanding of the FAdV-4–host cell interaction mechanism remains limited.

The aim of this study was to gain more comprehensive insight into the infection of host cells by FAdV-4. We used transcriptomics for the first time to our knowledge to investigate gene expression in LMH cells at five characteristic stages (12, 24, 36, 48, and 60 h) post infection (hpi) with FAdV-4. We then elucidated the molecular mechanisms (e.g., viral invasion, immunosuppression, and metabolic pathways) affected by viral infection. Our findings provide novel information regarding the mechanisms underlying virus–host interactions and the pathogenesis of FAdV-4.

Materials and Methods

Virus and Cells

FAdV-4 isolate SX17 (GenBank: MF592716.1) and LMH cells were kindly provided by Dr. Jing-Yu Wang (Northwest A&F University, China). LMH cells (4 × 10⁵ cells/well) were cultured in six-well plates maintained at 37°C in a 5% CO₂ atmosphere in Dulbecco's modified Eagle's medium (DMEM, Gibco) with 10% fetal bovine serum (FBS). Cells were cultured for ~24 h to reach 80% confluence and subsequently infected with FAdV-4 (MOI = 0.1); uninfected LMH cells were used as a mock control. FAdV-4 infected with LMH Cells caused cytopathic effect (CPE), viral titers and viral loads as previously described (15).

RNA Isolation, Library Construction, and RNA Sequencing (RNA-Seq)

TRIzol reagent (Invitrogen; USA) was used to isolate total RNA from 18 samples (five groups) according to the manufacturer's instructions. RNA purity and concentration were determined by ultraviolet absorbance using a Nanodrop Nano Photometer spectrophotometer (IMPLEN, USA), while RNA integrity was evaluated using an Agilent 2100 Bioanalyzer (Agilent Technologies, USA). Total RNA (3 μg) was used to construct a cDNA library using an RNA Library Prep Kit for Illumina (New England Biolabs, USA), according to the instructions provided by the manufacturer. In brief, total RNA was treated with DNaseI (Invitrogen, Life Technologies, USA) to remove genomic DNA contamination. This was followed by enrichment of poly(A) mRNA using oligo d(T) magnetic beads (Invitrogen). mRNA was next treated with fragmentation buffer, and cleaved RNA was copied into first-strand cDNA fragments using reverse transcriptase (Invitrogen) and random hexamer primers. Second-strand cDNA fragments were synthesized using Second Strand Synthesis Enzyme Mix, including DNA polymerase I, buffer, deoxyribonucleoside triphosphate, and RNase H. After purification and paired-end repair, cDNA fragments were ligated to sequencing adapters and amplified by PCR to obtain the final paired-end library. Libraries were subsequently sequenced using the Illumina HiSeq 4000 sequencing platform (Novogene, Shanghai, China).

Transcriptome Data Analysis

To ensure high accuracy and reliability of the results of subsequent analyses, the original ordinal number was determined by filtering to obtain high-quality sequencing data (clean data). The original sequencing data were filtered to remove low-quality reads as follows: (1) removal of adapter sequences with a quality score < 20 and sequences with an N base rate of raw reads > 10%; and (2) removal of adapter and quality-trimmed short segments <25 bp in length. The remaining reads were termed “clean reads” and stored in FASTQ format. For mapping clean reads, Bowtie 2 was used to build reference genes and HISAT (hierarchical indexing for spliced alignment of transcripts) was used to reference the genome.

For gene expression analysis, gene expression levels were calculated and RSEM (RNA-Seq by Expectation-Maximization) was used to estimate the unigene FPKM value. Default criteria of | log2 fold-change | > 1 and false discovery rate < 0.05 denoted significant changes in differentially expressed genes (DEGs). Gene Ontology (GO), and pathway annotation and enrichment analyses were performed using the National Center for Biotechnology Information GO (http://www.geneontology.org) and Kyoto Encyclopedia of Genes and Genomes (KEGG) (http://www.genome.jp/kegg/) pathway databases. Cluster and Java Tree view software were used to conduct hierarchical cluster analysis of gene expression patterns. All gene expression data presented were obtained from three independent replicates. Original sequence data were submitted to the NCBI database (accession number: SUB9560975).

Quantitative Reverse-Transcription PCR (qRT-PCR) to Validate Gene Expression

Nine genes were randomly selected for qRT-PCR analysis to validate the reproducibility and repeatability of the DEGs identified from transcriptome sequencing. Primer sequences used for qRT-PCR are listed in Supplementary Table 1. All primers were synthesized by Sangon Biotech (Shanghai, China). Total RNA was extracted from cell-line pellets using TRIzol reagent (Invitrogen, USA), according to the instructions provided by the manufacturer. RNA samples were reverse transcribed into cDNAs using a PrimeScript™ RT reagent Kit (TaKaRa, China). Relative gene expression values were calculated using the 2^−ΔΔCT method as previously described (16). mRNA levels were calculated with reference to those of the housekeeping gene β-actin. Each reaction was performed in triplicate.

Results

Growth Kinetics of FAdV-4-Infected LMH Cells

We determined the growth kinetics of FAdV-4 replication in primary LMH cells by monitoring the cytopathic effect and viral titers at 0, 12, 24, 36, 48, and 60 hpi. Viral load was observed at 12 hpi and increased slightly between 12 and 24 hpi. Subsequently, viral load increased significantly from 24 and 48 hpi; viral replication plateaued at 48 hpi (Figure 1). The wide range of infection time points used in this study allowed us to identify fluctuations in gene expression. This approach may enhance our understanding of the molecular response to infection with FAdV-4.

FIGURE 1

Figure 1. Characteristics of FAdV-4-infected LMH cells assessed via viral replication analysis. One-step growth curve of FAdV-4-infected LMH cells at a multiplicity of infection of 0.1, obtained using quantitative real-time PCR. Each time point represents the mean ± SEM of the Hexon gene from the FAdV-4 genome. SEM, standard error of the mean.

RNA-Seq and Read Assembly

We used LMH cells at five time points (12, 24, 36, 48, and 60 hpi) after infection with FAdV-4 (each condition with three biological replicates); the data obtained at these time points were compared with those obtained at 0 hpi (three biological replicates of mock control). A total of 330.568 GB of clean data (sequencing data after quality control) were obtained. The average amount of clean data for each sample was 15.741 GB. The percentage of Q30 bases was > 92.28%, and the GC content ranged between 46.09 and 55.24%. Ribosomal RNA sequences were removed, and low-quality reads were discarded. We then mapped 18 high-quality reads to the chicken (Gallus domesticus) genome using HISAT (Table 1). The mean percentage of high-quality reads aligned to the reference genome was 84.96%. Qualifying sample data were evaluated according to multiple metrics. RNA sequencing raw data were deposited in the National Center for Biotechnology Information Gene Expression Omnibus (GEO) under accession number SRP149860.

TABLE 1

Table 1. Statistics of the RNA-seq dataset.

Identification of Differentially Expressed Transcripts

A total of 8,242 DEGs were identified at five important time points (12, 24, 36, 48, and 60 hpi) (Figure 2A). At 0 and 12, 12 and 24, 24 and 36, 36 and 48, and 48 and 60 hpi, we identified 589, 956, 724, 1,301, and 518 upregulated DEGs and 535, 736, 1,045, 450, and 1,442 downregulated DEGs (Supplementary Files 1–5). We identified 37 common upregulated and downregulated DEGs at all time points (Supplementary Files 11). Moreover, 10, 410, 336, 239, and 534 DEGs were unique to comparisons of 0 and 12, 12 and 24, 24 and 36, 36 and 48, and 48 and 60 hpi, respectively (Figure 2B).

FIGURE 2

Figure 2. Identification and clustering of DEGs. (A) Number of differentially expressed genes at five time points after infection with FAdV-4. (B) Up Set plot showing unique and common DEGs in response to infection with FAdV-4 at five time points.

Hierarchical clustering analysis revealed gradual changes in gene expression across at 0 and 12, 12 and 24, 24 and 36, 36 and 48, and 48 and 60 hpi; the left side of the hierarchical clustering gene tree shows further classification of DEGs (Supplementary Files 10).

GO Enrichment Analyses of DEGs

We conducted GO analysis to identify the biological processes involved in FAdV-4 infection. Among DEGs between 0 and 12 hpi, 8,348 unigenes were assigned GO terms in the cellular component (676), molecular function (1,347), and biological process (6,325) categories, respectively. For the 12 and 24 hpi comparison, 9,742 unigenes were assigned to 834, 1,651, and 7,257 GO terms in the aforementioned categories, respectively. For the 24 and 36 hpi comparison, 9,695 unigenes were assigned to 844, 1,681, and 7,170 GO terms in the aforementioned categories, respectively. For the 36 and 48 hpi comparison, 9,856 unigenes were assigned to 831, 1,636, and 7,389 GO terms in the aforementioned categories, respectively. Finally, for the 48 and 60 hpi comparison, 10,482 unigenes were assigned to 906, 1,707, and 7,869 GO terms in the aforementioned categories, respectively. We next selected the 10 most significant molecular functions and biological processes with altered gene expression among viral infection stages (Figure 3).

FIGURE 3

Figure 3. Analysis of DEGs using Gene Ontology (GO) terms. Top 10 GO categories significantly (P-value < 0.05) enriched in biological processes at five time points. The x-axis represents the number of genes, and the y-axis represents functional groups.

KEGG Pathway Enrichment Analysis of DEGs

Enrichment analysis of KEGG pathways can further identify signal transduction pathways in which DEGs participate. We mapped 8,242 DEGs to referential canonical pathways via KEGG database analysis. DEGs identified at five time points, 0 and 12 hpi (620 DEGs), 12 and 24 hpi (1,226 DEGs), 24 and 36 hpi (724 DEGs), 36 and 48 hpi (1,301 DEGs), and 48 and 60 (1,960 DEGs), were enriched to 138, 324, 321, 329, and 326 KEGG pathways, respectively. The 10 most enriched pathways were determined according to P-value < 0.05 (Figure 4). We identified many DEGs playing important roles related to viral invasion, immunosuppression, and metabolic pathways (Supplementary Files 6–8).

FIGURE 4

Figure 4. Analysis of DEGs using Kyoto Encyclopedia of Genes and Genomes (KEGG) terms. Top 10 KEGG pathways significantly (P-value < 0.05) enriched in biological processes at five time points. The x-axis represents the number of genes, and the y-axis represents functional groups. Red and blue portions of the bar represent the number of upregulated and downregulated genes, respectively.

In order to further analyze the differential expression host immunosuppression genes response during FAdV-4 infection. KEGG enrichment analyze, the DEGs were significantly enriched in immune-related signaling pathways, such as, TNF signaling pathway, IL-17 signaling pathway, Apoptosis, P13K-Akt signaling pathway, cytokine-cytokine receptor interaction, Th17 cell differentiation, Th1 and Th2 cell differentiation, Toll-like receptor signaling path way (Supplementary Files 9A–E). Meanwhile, Apoptosis and immune-related signaling pathway the changes in important gene expression shown by the heatmap (Figures 5A,B).

FIGURE 5

Figure 5. Heat maps shows the hierarchical clustering of genes involved in apoptosis genes (A) and Immune-related genes (B) at five time points. The color in the heat map represents the normalized isotig number. The color in the heat map represents gene expression changes (green, high gene expression; pink, low gene expression; white, no gene expression). DEGs are labeled on the left of the heat-map line.

Validation of RNA-Seq Data Using qRT-PCR

To further evaluate our library of differentially expressed unigenes, we used qRT-PCR to detect the expression of nine pro-inflammatory response genes (early growth response 1 [EGR1], suppressor of cytokine signaling 3 [SOCS3], thrombospondin 1 [THBS1], interleukin-8 [IL8], suppressor of cytokine signaling 1 [SOCS1], TNF receptor superfamily member 8 [CD30], C-C motif chemokine ligand 4 [CCL4], phosphoenolpyruvate carboxykinase1 [PCK1], Apoptosis regulator proteins Bcl-2 [Bcl-2] and TNF receptor superfamily member 21 [TNFRSF21] (Figure 6). Comparison of the fold-change in expression of the selected genes observed by RNA-seq and qRT-PCR showed that the two data sets were largely consistent.

FIGURE 6

Figure 6. Validation of RNA-seq results using quantitative reverse-transcription PCR (qRT-PCR). Expression profiles of nine DEGs playing important roles in viral defense. qRT-PCR verified that mRNA levels of these genes were increased after viral infection, consistent with the altered expression profiles observed in transcriptome data. β-actin was used for normalization. The x-axis represents the different time points, and the y-axis represents the normalized fold-change of transcripts.

Discussion

First emerging in 2015 in China, cases of fowl adenovirus serotype 4 (FAdV-4) infection in broilers, laying hens, ducks, and geese have caused economic and social losses globally, including in India, South Korea, Hungary, Japan, Europe, and China (17, 18). As a recently discovered adenovirus, research on FAdV-4 remains insufficient. In this study, we employed transcriptome sequencing technology to explore crucial factors affecting host–virus interaction at five characteristic stages (12, 24, 36, 48, and 60 hpi) of FAdV-4 infection in LMH cells. Many important genes associated with viral invasion, host immunosuppression, and host metabolic pathways were differentially expressed.

When FAdV-4 infects a host, a complex battle between host cells and the virus takes place. Continuous replication of the virus damages host cells to different degrees and alters their physiological functions, leading to differential gene expression. DEGs identified following FAdV-4 infection of LMH cells can be divided into five infection periods based on their differential expression patterns. At the first stage of infection (0 to 12 hpi), we detected 1,124 DEGs between 0 and 12 hpi. At 12 hpi, entry of FAdV-4 into LMH cells induces a strong immune response, and host cells enter a static phase of balanced offensive and defensive attacks (19, 20). The second stage of infection between 12 and 24 hpi is a crucial time point for FAdV-4 and host cell interaction; beyond 24 hpi, the pathogen prevails, resulting in marked viral replication. We identified 1,692 DEGs between 12 and 24 hpi (20). These results were similar to those observed in previous studies (19). In the third stage of infection (24 to 36 hpi), we detected 1,769 DEGs between 24 and 36 hpi; during this time, the virus controls the metabolic machinery of host cells to produce more progeny virus. In the fourth stage of infection (36 and 48 hpi), we detected 1,751 DEGs between 48 and 36 hpi; among these different genes, expression of 1,301 DEGs increased dramatically while 450 DEGs were downregulated. The peak of FAdV-4 DNA replication occurs at 48 hpi and affects the host cell cycle, antiviral response, and cell proliferation as the virus changes its strategy from reproduction to efficient release and spread of progeny through destruction of the cell. In the late stage of virus infection (48 to 60 hpi), we detected 1,906 DEGs between 48 and 60 hpi; during this time, the number of DEGs was significantly increased compared with earlier stages of infection, with the majority being downregulated at 60 hpi compared with 48 hpi. This downregulation was non-specific and a consequence of imminent cell death. At this stage, FAdV-4 changes its strategy from reproduction to release of viral progeny and spread to adjacent cells through destruction of the cell, leading to gradual destruction of infected tissue.

From these five characteristic stages of FAdV-4 infection of host cells, we determined that regulation of host cell gene expression by FAdV-4 is specific, with FAdV-4 invading host cells to replicate its genome, produce new infectious virions, and spread virions to adjacent cells. These results are consistent with findings from previous studies on human adenovirus (19, 20). In addition, at these five important time points, we found 37 common differentially upregulated or downregulated genes, Among these common differential gene, the IL-8,SOCS3, NFKBIA, CCL4, NFAIP3, and TNIP2, which are associated with immune signaling pathways. Especially, TNFAIP3 belong to tumor necrosis factor were influenced during FAdV-4 infection process and play vital roles in the inflammatory response (21, 22). However, IL-8 uncontrolled expression can cause several inflammatory responses infected tissues and ultimately lead to death (23, 24). These results indicated that these common upregulated or downregulated genes play a common role for these genes in host response to viral infection. Importantly, our study covered a wide range of time points of LMH cell infection with FAdV-4, allowing us to identify a large number of DEGs relevant to viral invasion, immunity response, apoptosis, and metabolic pathways. In-depth analysis of these pathways will help us understand the mechanisms involved in FAdV-4–host interaction.

Viral Invasion

Viral entry into cells is an extremely complex process involving multiple signaling pathways and receptor proteins. It is the premise for virus replication to cause pathological changes and consequently induce disease. During adenovirus infection, viral invasion occurs through adhesion to the cell surface and binding to host receptors, followed by membrane fusion and viral entry into the cell (20). In the early stages of FAdV-4 infection, we found that several important DEGs (e.g., THBS1, EGR1, and SOCS3) were upregulated. As previously reported, THBS1 plays a key role in pathophysiology of hantavirus infection (25, 26). EGR1 facilitates the replication of enterovirus 71 by directly targeting viral genomic RNA (27). The mitogen-activated protein kinase (MEK) signaling pathway can activate EGR1 during Foot-and-mouth disease infection (28). An important finding of the current and previous studies is that SOCS proteins are upregulated following viral infection; expression of SOCS1 or SOCS3 following viral infection promotes the replication of multiple viruses (29, 30). Further studies are warranted to clarify the main role and mechanism of SOCS protein following FADV-4 infection. Based on these observations, THBS1, EGR1, and SOCS3 are a good candidate gene for further functional studies.

KEGG analysis revealed that DEGs at 12 and 24 hpi were mainly enriched in the regulatory signaling pathways of receptor binding, including cell adhesion molecules, phosphatidylinositol 3-kinase-AKT (PI3K-AKT) pathway, toll-like receptor pathway, and endocytosis pathway, which are closely related to adenovirus infection (8, 19, 20). Previous studies have shown that cell adhesion molecules, as important cellular receptors, help numerous other viruses (e.g., human adenoviruses) to attach to and enter into target cells (31). The PI3K-AKT signaling pathway (PIK3AP1, PIK3CB, and PIK3R5) is closely related to adenovirus infection and specifically involves mediation of PI3K activation (32). The toll-like receptor signaling pathway can recognize viral nucleic acids participating in FAdV-4 infection (33, 34). The endocytosis pathway plays a critical role in mediation of different receptors during adenoviruses entry into host cells (35, 36).

Therefore, further investigation the role of candidate genes and signal pathways may play essential roles in the regulation of FAdV-4 infection, will not only provide new molecular targets for the development of antiviral drugs but also contribute to research vaccines in poultry.

Metabolic Pathways

As strictly intracellular parasites, viruses depend on manipulation of the host metabolism to obtain the necessary energy and nutrients for their own survival and proliferation. Considerable evidence has shown that viral infection (e.g., adenovirus 5, influenza A virus, Newcastle disease virus, and herpes simplex virus-1) alters the metabolism of host cells (37–40). Different viruses require unique metabolic changes to successfully spread their progeny (41). In this study, we identified numerous significantly altered DEGs, which were enriched in important metabolism pathways, including carbohydrate metabolism, amino acids, lipid metabolism, and biosynthesis of secondary metabolites (42) (Supplementary File 6).

Carbohydrate metabolism is the primary carbon and energy source and plays a vital role in mammalian cells during homeostasis. Numerous metabolomic studies have identified an increase in carbohydrate metabolism caused by viral infection (39). For example, adenovirus can promote anabolic glucose metabolism in host cells and viral replication through activation of early region 4 open reading frame 1 (E4ORF1)-induced MYC. Studies of duck hepatitis A virus genotype 3 also demonstrated that glucose metabolism is a crucial hub during virus infection (43). In our study, genes associated with the glycolysis metabolic pathway, including PCK1, phosphofructokinase, muscle (PFKM), and alcohol dehydrogenase 1C (class I), gamma polypeptide (ADH1C), were upregulated at 12 and 48 hpi. Subsequently, during the rapid replication period of FAdV-4, expression levels of genes encoding various carbohydrate branch metabolites involved in fructose and mannose metabolism, glycosaminoglycan biosynthesis, tricarboxylic acid cycle, and the pentose phosphate pathway were significantly increased (44).

Amino acids are the most basic building blocks of living organisms. They serve as the essential substrate for the biosynthesis of many important biomolecules. In addition, amino acids can provide intermediate metabolites for glycolysis and the tricarboxylic acid cycle (45). Previous studies suggested that infection of human breast and bronchial epithelial cells with wild-type adenovirus 5 can increase glucose consumption. Additionally, adenovirus infection can selectively increase the expression of asparagine synthetase, which regulates the replication of adenovirus (46). Recent studies have found that FAdV-4 influences the metabolism of arginine to promote its replication (47). In the present study, KEGG enrichment analysis showed that amino acid metabolism pathways were upregulated during FAdV-4 infection; moreover, the production of arginine, glycine, serine, threonine, tyrosine, valine, leucine, isoleucine, alanine, aspartate, glutamate, cysteine, and methionine was increased during infection with FAdV-4. Collectively, these finding suggest that amino acid metabolism may contribute to FAdV-4 replication and deserves attention in future research studies.

Lipid is an important component of human cell tissue and a key structural component of the cell membrane, which plays a crucial role in a variety of biological processes. Numerous studies have demonstrated that viruses can take advantage of several core lipid metabolic pathways to successfully adhere to and invade host cells (38, 48). In the present study, genes involved in glycerol lipid metabolism, ether lipid metabolism, and sphingolipid metabolism were upregulated by FAdV-4 infection (Supplementary File 6).

In addition to the three major metabolic processes of host cells described above, we found that steroid biosynthesis and nitrogen metabolism were also altered after infection with FAdV-4. These are important host metabolic pathways that control the process and outcome of virus infections (49, 50). Further studies on metabolic genes and signaling pathways associated with FAdV-4 infection may provide some reference value for studying the pathogenesis and new diagnostic methods of FAdV-4.

Immunosuppression

FAdV-4 is the primary pathogen of HHS, which has an immunosuppressive potential (51, 52). Previous studies have shown that apoptosis and inflammatory signaling activation are potent inducers of immunosuppressive mechanisms (4, 53). Interestingly, in this study, we observed that a substantial number of DEGs were enriched in the apoptosis-related and inflammatory pathways.

Apoptosis is a physiological and pathological phenomenon. Excessive apoptosis can lead to physical dysfunction, which results in a bystander effect of induced immunosuppression (4, 54). In the present study, we observed that infection of LMH cells with FAdV-4 activated a number of DEGs that regulate the apoptosis, apoptosis-fly, and apoptosis-multiple species signaling pathways, mainly containing pro-survival genes e.g., TMBIM6 (Bax inhibitor), BCL2L1 (Bcl-xL), BCL2 apoptosis regulator (BCL2); pro-apoptotic genes e.g., caspase 8 (CASP8), caspase 9 (CASP9), transcription factor subunit (FOS), Jun proto-oncogene, AP-1 transcription factor subunit (JUN), tubulin alpha 1a (TUBA1A), TUBA1B, TUBA1C, and TUBA3E (Figure 5A). CASP8 and CASP9 are the indicator of the mitochondrion-mediated apoptosis pathway. Once activated, this gene triggers characteristic apoptotic changes in cells (4, 55). Interestingly, we observed that BCL2 and BCL2 associated X, apoptosis regulator (BAX) were released. The BCL2 dimer regulates apoptosis by combining with BAX to form a heterodimer. This process leads to loss of the inhibitory effect of BCL2 on apoptosis, maintaining the stability of the host internal environment (56, 57). Furthermore, JUN and FOS regulate the gene encoding the Fas cell surface death receptor (FAS) ligand and play an important role as potent inducers of apoptosis (58–60).

A moderate inflammatory response can maintain immune system homeostasis (61, 62). However, persistent infection with a pathogenic microorganism may induce severe inflammation. In this study, infection of LMH cells with FAdV-4 activated numerous immune signaling pathways (Supplementary File 7). Our results indicated that these immune pathways were associated with several common DEGs, mainly those encoding pro-inflammatory genes e.g., IL1, IL-6, IL-8, IL-15, C-X3-C motif chemokine ligand 1 (CX3CL1), C-X-C motif chemokine receptor 4 (CXCR4), C-X-C motif chemokine ligand 12 (CXCL12), TNF superfamily member 10 (TNFSF10), TNFSF11, TNFSF13B, TNFSF18, TNFSF21 and anti-inflammatory genes e.g., SOCS3, CDKN1A, IL10, CDKN2A (Figure 5B). IL-6 and IL-8 are the most important pro-inflammatory cytokines in FAdV-4 infection (25). IL21R belongs to the Th17 cell differentiation pathway; this is a newly discovered subset of effector T helper cells, which mainly produce cytokines. Previous studies have demonstrated that IL21R is crucial for the differentiation of native T cells into Th17 cells and also plays a key role in the development of autoimmune disease (63). The Th1 and Th2 cell differentiation pathway was also activated, as evidenced by the upregulation of genes encoding cytokine-secreting Th1 (pro-inflammatory) cytokines (IL-17, IL-2). Th1 participates in cellular immunity, delays hypersensitivity inflammatory response, and plays a key role in host resistance to intracellular bacteria and other pathogens (64). TNFRSF13B and TNFSF21, members of the tumor necrosis factor receptor family, affect multiple cytokine pathways. They play an important role in initiating signal transduction and participate in complex immune responses, pathological injury, and the mediation of immunosuppression (61, 65, 66). We observed that, in LMH cells infected with FAdV-4, TNFRSF13B was downregulated, whereas TNFSF21 was upregulated. These findings are consistent with data obtained from studies of the above immune-related diseases (34, 61).

Altogether, the balance between these pro-survival genes, pro-apoptotic genes, pro-inflammatory and anti-inflammatory factors may be play important role for FAdV-4 immunosuppression. Further explore those of genes were enriched in the apoptosis-related and inflammatory pathways is essential to clarify the pathogenesis mechanisms of avian adenovirus and identify potential therapeutic strategies.

Conclusion

This transcriptomic study offers novel insights for better understanding the response of LMH cells to infection with FAdV-4. Our data highlight the mechanisms involved in FAdV-4 invasion and FAdV-4-induced immunosuppression, as well as the effects of metabolic pathways on FAdV-4 infection. The basic data obtained in this study indicate the complexity of host–FAdV-4 interaction and provide new ideas for the development of effective vaccines against FAdV-4.

Data Availability Statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.

Author Contributions

XW conceived and designed the experiments and wrote sections of the manuscript. XW and BW performed the experiments. XZ, LM, and XL analyzed the data. MZ contributed reagents, materials, and analysis tools. All authors contributed to manuscript revision, read, and approved the submitted version.

Funding

This work was supported by the National Natural Science Foundation of China [Grant No. NSFC32002313] and the Key Scientific Research Project of Henan Higher Education Institutions of China [Grant No. 20B230001].

Conflict of Interest

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.

Publisher's Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fvets.2022.871038/full#supplementary-material

References

1. Zhao M, Duan X, Wang Y, Gao L, Cao H, Li X, et al. A novel role for PX, a structural protein of fowl adenovirus serotype 4 (FAdV4), as an apoptosis-inducer in leghorn male hepatocellular cell. Viruses. (2020) 12:228. doi: 10.3390/v12020228

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Pan Q, Yang Y, Gao Y, Qi X, Liu C, Zhang Y, et al. An inactivated novel genotype fowl adenovirus 4 protects chickens against the hydropericardium syndrome that recently emerged in China. Viruses. (2017) 9:216. doi: 10.3390/v9080216

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Yu G, Wang Y, Zhang M, Lin Y, Diao Y. Pathogenic, phylogenetic, and serological analysis of group I fowl adenovirus serotype 4 SDSX isolated from Shandong, China. Front Microbiol. (2018) 9:2772. doi: 10.3389/fmicb.2018.02772

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Niu Y, Sun Q, Shi Y, Ding Y, Li Z, Sun Y, et al. Immunosuppressive potential of fowl adenovirus serotype 4. Poult Sci. (2019) 98:3514–22. doi: 10.3382/ps/pez179

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Chen Z, Shi S, Qi B, Lin S, Chen C, Zhu C, et al. Hydropericardium syndrome caused by fowl adenovirus serotype 4 in replacement pullets. J Vet Med Sci. (2018) 81:245–51. doi: 10.1292/jvms.18-0168

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Wang Z, Zhao JJV. Pathogenesis of hypervirulent fowl adenovirus serotype 4: the contributions of viral and host factors. Virus. (2019) 11:741. doi: 10.3390/v11080741

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Zsabc D, Bxabc D, Slabc D, Xrabc D, Yhabc D, Jlabc D, et al. Genetic characterization of fowl adenovirus serotype 4 isolates in Southern China reveals potential cross-species transmission. Infect Genet Evol. (2019) 75:103928. doi: 10.1016/j.meegid.2019.103928

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Wang Y, Wang Y, Wang B, Mei X, Jiang S, Li JPS, et al. Protocatechuic acid improved growth performance, meat quality, and intestinal health of Chinese yellow-feathered broilers. Poult Sci. (2019) 98:3138–49. doi: 10.3382/ps/pez124

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Shao H, Lu Y, Wang W, Li T, Zhang J, Wan Z, et al. Two novel monoclonal antibodies against fiber-1 protein of FAdV-4 and their application in detection of FAdV-4/10. BMC Vet Res. (2019) 15:232. doi: 10.1186/s12917-019-1987-5

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Quigley M, Martinez J, Huang X, Yang Y. A critical role for direct TLR2-MyD88 signaling in CD8 T-cell clonal expansion and memory formation following vaccinia viral infection. Blood. (2009) 113:2256–64. doi: 10.1182/blood-2008-03-148809

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Li C, Wang T, Zhang Y, Wei F. Evasion mechanisms of the type I interferons responses by influenza A virus. Crit Rev Microbiol. (2020) 46:420–32. doi: 10.1080/1040841X.2020.1794791

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Yu J, Sun X, Goie JYG, Zhang Y. Regulation of host immune responses against influenza A virus infection by mitogen-activated protein kinases (MAPKs). Microorganisms. (2020) 8:1067. doi: 10.3390/microorganisms8071067

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Xie Q, Wang W, Li L, Kan Q, Fu H, Geng T, et al. Domain in Fiber-2 interacted with KPNA3/4 significantly affects the replication and pathogenicity of the highly pathogenic FAdV-4. Virulence. (2021) 12:754–65. doi: 10.1080/21505594.2021.1888458

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Zhang W, Guo X, Yin F, Zou X, Hou W, Lu Z. Fiber modifications enable fowl adenovirus 4 vectors to transduce human cells. J Gene Med. (2021) 23:e3368. doi: 10.1002/jgm.3368

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Wen B, Wang X, Yang L, Wang T, Hou X, Qi X, et al. Transcriptome analysis reveals the potential role of long noncoding RNAs in regulating fowl adenovirus serotype 4-induced apoptosis in leghorn male hepatocellular cells. Viruses. (2021) 13:1623. doi: 10.3390/v13081623

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Yue H, Lei XW, Yang FL, Li MY, Tang C. Reference gene selection for normalization of PCR analysis in chicken embryo fibroblast infected with H5N1 AIV. Virol Sin. (2010) 25:425–31. doi: 10.1007/s12250-010-3114-4

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Hou L, Chen X, Wang J, Li J, Yang H. A tandem mass tag-based quantitative proteomic analysis of fowl adenovirus serotype 4-infected LMH cells. Vet Microbiol. (2021) 255:109026. doi: 10.1016/j.vetmic.2021.109026

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Xie S, Zhang J, Chen H, Wang W, Wang P, Xie Q, et al. Development of colloidal gold-based test strip for rapid detection of serotype 4 fowl adenovirus. J Virol Methods. (2021) 296:114231. doi: 10.1016/j.jviromet.2021.114231

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Zhao H, Granberg F, Pettersson UJV. How adenovirus strives to control cellular gene expression. Virology. (2007) 363:357–75. doi: 10.1016/j.virol.2007.02.013

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Wang XP, Qi XF, Yang B, Chen SY, Wang JY. RNA-Seq analysis of duck embryo fibroblast cell gene expression during the early stage of egg drop syndrome virus infection. Poult Sci. (2018) 98:404–12. doi: 10.3382/ps/pey318

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Alsamman AM, Zayed H. The transcriptomic profiling of SARS-CoV-2 compared to SARS, MERS, EBOV, and H1N1. PLoS ONE. (2020) 15:e0243270. doi: 10.1371/journal.pone.0243270

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Mo Y, Liu Y, Lu A, Zhang H, Tang L. Role of circRNAs in viral infection and their significance for diagnosis and treatment (Review). Int J Mol Med. (2021) 47:88. doi: 10.3892/ijmm.2021.4921

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Yiemwattana I, Chaisomboon N, Jamdee K. Antibacterial and anti-inflammatory potential of Morus alba stem extract. Open Dent J. (2018) 12:265–74. doi: 10.2174/1874210601812010265

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Meng K, Yuan X, Yu J, Zhang Y, Ai W, Wang Y. Identification, pathogenicity of novel fowl adenovirus serotype 4 SDJN0105 in Shandong, China and immunoprotective evaluation of the newly developed inactivated oil-emulsion FAdV-4 vaccine. Viruses. (2019) 11:627. doi: 10.3390/v11070627

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Hamzić E, Buitenhuis B, Hérault F, Hawken R, Abrahamsen MS, Servin B. Genome-wide association study and biological pathway analysis of the Eimeria maxima response in broilers. Genet Sel Evol. (2015) 47:91. doi: 10.1186/s12711-015-0170-0

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Khaiboullina SF, Morzunov SP, St Jeor SC, Rizvanov AA, Lombardi VC. Hantavirus infection suppresses thrombospondin-1 expression in cultured endothelial cells in a strain-specific manner. Front Microbiol. (2016) 7:1077. doi: 10.3389/fmicb.2016.01077

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Song Y, Cheng X, Yang X, Zhao R, Wang P, Han Y, et al. Early growth response-1 facilitates enterovirus 71 replication by direct binding to the viral genome RNA. Int J Biochem Cell Biol. (2015) 62:36–46. doi: 10.1016/j.biocel.2015.02.012

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Zhu Z, Du X, Li P, Zhang X, Yang F, Cao W, et al. Early growth response gene-1 suppresses foot-and-mouth disease virus replication by enhancing type i interferon pathway signal transduction. Front Microbiol. (2018) 9:2326. doi: 10.3389/fmicb.2018.02326

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Oh SJ, Gim JA, Lee JK, Park H, Shin OS. Coxsackievirus B3 infection of human neural progenitor cells results in distinct expression patterns of innate immune genes. Viruses. (2020) 12:325. doi: 10.3390/v12030325

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Seong RK, Lee JK, Shin OS. Zika Virus-induction of the suppressor of cytokine signaling 1/3 contributes to the modulation of viral replication. Pathogens. (2020) 9:163. doi: 10.3390/pathogens9030163

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Persson BD, Lenman A, Frängsmyr L, Schmid M, Ahlm C, Plückthun A, et al. Lactoferrin-hexon interactions mediate CAR-independent adenovirus infection of human respiratory cells. J Virol. (2020) 94:e00542-20. doi: 10.1128/JVI.00542-20

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Rajala MS, Rajala RV, Astley RA, Butt AL, Chodosh J. Corneal cell survival in adenovirus type 19 infection requires phosphoinositide 3-kinase/Akt activation. J Virol. (2005) 79:12332–41. doi: 10.1128/JVI.79.19.12332-12341.2005

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Bolduc JF, Ouellet M, Hany L, Tremblay MJ. Toll-like receptor 2 ligation enhances HIV-1 replication in activated CCR6+ CD4+ T cells by increasing virus entry and establishing a more permissive environment to infection. J Virol. (2017) 91:e01402-16. doi: 10.1128/JVI.01402-16

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Zhang J, Zou Z, Huang K, Lin X, Chen H, Jin JVM, et al. Insights into leghorn male hepatocellular cells response to fowl adenovirus serotype 4 infection by transcriptome analysis. Vet Microbiol. (2018) 214:65–74. doi: 10.1016/j.vetmic.2017.12.007

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Kasman LM, Barua S, Lu P, Rege K, Voelkel-Johnson C. Polymer-enhanced adenoviral transduction of CAR-negative bladder cancer cells. Mol Pharm. (2009) 6:1612–9. doi: 10.1021/mp9000958

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Kim SY, Lee SJ, Lim SJ. Formulation and in vitro and in vivo evaluation of a cationic emulsion as a vehicle for improving adenoviral gene transfer. Int J Pharm. (2014) 475:49–59. doi: 10.1016/j.ijpharm.2014.08.024

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Sen S, Kaminiski R, Deshmane S, Langford D, Khalili K, Amini S, et al. Role of hexokinase-1 in the survival of HIV-1-infected macrophages. Cell Cycle. (2015) 14:980–9. doi: 10.1080/15384101.2015.1006971

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Mesquita I, Estaquier J. Viral manipulation of the host metabolic network. Exp Suppl. (2018) 109:377–401. doi: 10.1007/978-3-319-74932-7_10

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Liu P, Yin Y, Gong Y, Qiu X, Sun Y, Tan L, et al. In vitro and in vivo metabolomic profiling after infection with virulent newcastle disease virus. Viruses. (2019) 11:962. doi: 10.3390/v11100962

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Tian X, Zhang K, Min J, Chen C, Cao Y, Ding C, et al. Metabolomic analysis of influenza A virus A/WSN/1933 (H1N1) infected A549 cells during first cycle of viral replication. Viruses. (2019) 11:1007. doi: 10.3390/v11111007

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Sanchez EL, Lagunoff M. Viral activation of cellular metabolism. Virology. (2015) 479–80:609–18. doi: 10.1016/j.virol.2015.02.038

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Raniga K, Liang C. Interferons: reprogramming the metabolic network against viral infection. Viruses. (2018) 10:36. doi: 10.3390/v10010036

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Liang S, Tang J, Wang X, Wen Z, Xie M, Wu Y, et al. Effect of duck hepatitis A virus genotype 3 infection on glucose metabolism of Pekin ducklings and underlying mechanism. Gene. (2020) 748:144710. doi: 10.1016/j.gene.2020.144710

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Thai M, Graham NA, Braas D, Nehil M, Komisopoulou E, Kurdistani SK, et al. Adenovirus E4ORF1-induced MYC activation promotes host cell anabolic glucose metabolism and virus replication. Cell Metab. (2014) 19:694–701. doi: 10.1016/j.cmet.2014.03.009

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Zeng D, Cui J, Yin Y, Xiong Y, Liu M, Guan S, et al. Metabolomics analysis in different development stages on SP0 generation of rice seeds after spaceflight. Front Plant Sci. (2021) 12:700267. doi: 10.3389/fpls.2021.700267

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Krall AS, Xu S, Graeber TG, Braas D, Christofk HR. Asparagine promotes cancer cell proliferation through use as an amino acid exchange factor. Nat Commun. (2016) 7:11457. doi: 10.1038/ncomms11457

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Lin Z, Huang R, Zhou J, Chen Y, Xu L, Gao Y, et al. Fowl adenovirus serotype 4 influences arginine metabolism to benefit replication. Avian Dis. (2020) 64:16–22. doi: 10.1637/0005-2086-64.1.16

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Zhang J, Lan Y, Sanyal S. Modulation of lipid droplet metabolism-a potential target for therapeutic intervention in flaviviridae infections. Front Microbiol. (2017) 8:2286. doi: 10.3389/fmicb.2017.02286

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Kumar N, Mishra B, Mehmood A, Mohammad A, Mukhtar MS. Integrative network biology framework elucidates molecular mechanisms of SARS-CoV-2 pathogenesis. iScience. (2020) 23:101526. doi: 10.1016/j.isci.2020.101526

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Zhang X, Zhao Q, Wu C, Xie Z, Ci X, Li H, et al. Nitrate is crucial for the proliferation of gut Escherichia coli caused by H9N2 AIV infection and effective regulation by Chinese herbal medicine ageratum-liquid. Front Microbiol. (2020) 11:555739. doi: 10.3389/fmicb.2020.555739

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Pei Y, Corredor JC, Griffin BD, Krell PJ, Nagy É. Fowl adenovirus 4 (FAdV-4)-based infectious clone for vaccine vector development and viral gene function studies. Viruses. (2018) 10:97. doi: 10.3390/v10020097

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Mohamed Sohaimi N, Bejo MH, Omar AR, Ideris A, Mat Isa N. Molecular characterization of fowl adenovirus isolate of Malaysia attenuated in chicken embryo liver cells and its pathogenicity and immunogenicity in chickens. PLoS ONE. (2019) 14:e0225863. doi: 10.1371/journal.pone.0225863

PubMed Abstract | CrossRef Full Text | Google Scholar

53. Shivachandra SB, Sah RL, Singh SD, Kataria JM, Manimaran K. Immunosuppression in broiler chicks fed aflatoxin and inoculated with fowl adenovirus serotype-4 (FAV-4) associated with hydropericardium syndrome. Vet Res Commun. (2003) 27:39–51. doi: 10.1023/A:1022058623634

PubMed Abstract | CrossRef Full Text | Google Scholar

54. Brune W. Inhibition of programmed cell death by cytomegaloviruses. Virus Res. (2011) 157:144–50. doi: 10.1016/j.virusres.2010.10.012

PubMed Abstract | CrossRef Full Text | Google Scholar

55. Qi X, Xu J, Wang Z, Wang X, Wang J. The role of hexon in egg drop syndrome virus (EDSV) inducing apoptosis in duck embryo fibroblast cells. Res Vet Sci. (2017) 114:395–400. doi: 10.1016/j.rvsc.2017.07.015

PubMed Abstract | CrossRef Full Text | Google Scholar

56. Fan X, Mu L. The role of heme oxygenase-1 (HO-1) in the regulation of inflammatory reaction, neuronal cell proliferation and apoptosis in rats after intracerebral hemorrhage (ICH). Neuropsychiatr Dis Treat. (2017) 13:77–85. doi: 10.2147/NDT.S120496

PubMed Abstract | CrossRef Full Text | Google Scholar

57. Qian J, Cao X, Shen Q, Cai YF, Lu W, Yin H, et al. Thrombopoietin promotes cell proliferation and attenuates apoptosis of aplastic anemia serum-treated 32D cells via activating STAT3/STAT5 signaling pathway and modulating apoptosis-related mediators. Cell Transplant. (2021) 30:963689720980367. doi: 10.1177/0963689720980367

PubMed Abstract | CrossRef Full Text | Google Scholar

58. Zhang X, Zhang L, Yang H, Huang X, Otu H, Libermann TA, et al. c-Fos as a proapoptotic agent in TRAIL-induced apoptosis in prostate cancer cells. Cancer Res. (2007) 67:9425–34. doi: 10.1158/0008-5472.CAN-07-1310

PubMed Abstract | CrossRef Full Text | Google Scholar

59. Asim M, Chaturvedi R, Hoge S, Lewis ND, Singh K, Barry DP, et al. Helicobacter pylori induces ERK-dependent formation of a phospho-c-Fos c-Jun activator protein-1 complex that causes apoptosis in macrophages. J Biol Chem. (2010) 285:20343–57. doi: 10.1074/jbc.M110.116988

PubMed Abstract | CrossRef Full Text | Google Scholar

60. Ge R, Wang Z, Zeng Q, Xu X, Olumi AF. F-box protein 10, an NF-kappaB-dependent anti-apoptotic protein, regulates TRAIL-induced apoptosis through modulating c-Fos/c-FLIP pathway. Cell Death Differ. (2011) 18:1184–95. doi: 10.1038/cdd.2010.185

PubMed Abstract | CrossRef Full Text | Google Scholar

61. Guo Y, Su A, Tian H, Zhai M, Li W, Tian Y, et al. Transcriptomic analysis of spleen revealed mechanism of dexamethasone-induced immune suppression in chicks. Genes. (2020) 11:513. doi: 10.3390/genes11050513

PubMed Abstract | CrossRef Full Text | Google Scholar

62. Zhang K, Huang Q, Deng S, Yang Y, Li J, Wang S. Mechanisms of TLR4-mediated autophagy and nitroxidative stress. Front Cell Infect Microbiol. (2021) 11:766590. doi: 10.3389/fcimb.2021.766590

PubMed Abstract | CrossRef Full Text | Google Scholar

63. Sonderegger I, Kisielow J, Meier R, King C, Kopf M. IL-21 and IL-21R are not required for development of Th17 cells and autoimmunity in vivo. Eur J Immunol. (2008) 38:1833–8. doi: 10.1002/eji.200838511

PubMed Abstract | CrossRef Full Text | Google Scholar

64. Wu X, Shou Q, Chen C, Cai H, Zhang J, Tang S, et al. An herbal formula attenuates collagen-induced arthritis via inhibition of JAK2-STAT3 signaling and regulation of Th17 cells in mice. Oncotarget. (2017) 8:44242–54. doi: 10.18632/oncotarget.17797

PubMed Abstract | CrossRef Full Text | Google Scholar

65. Schmidt CS, Liu J, Zhang T, Song HY, Sandusky G, Mintze K, et al. Enhanced B cell expansion, survival, and humoral responses by targeting death receptor 6. J Exp Med. (2003) 197:51–62. doi: 10.1084/jem.20020617

PubMed Abstract | CrossRef Full Text | Google Scholar

66. Dulwich KL, Giotis ES, Gray A, Nair V, Skinner MA, Broadbent AJ. Differential gene expression in chicken primary B cells infected ex vivo with attenuated and very virulent strains of infectious bursal disease virus (IBDV). J Gen Virol. (2017) 98:2918–30. doi: 10.1099/jgv.0.000979

PubMed Abstract | CrossRef Full Text | Google Scholar