Next Article in Journal
Meiotic Chromosome Contacts as a Plausible Prelude for Robertsonian Translocations
Next Article in Special Issue
Transcriptome Analysis Reveals Inhibitory Effects of Lentogenic Newcastle Disease Virus on Cell Survival and Immune Function in Spleen of Commercial Layer Chicks
Previous Article in Journal
Altered Expression of GABAergic Markers in the Forebrain of Young and Adult Engrailed-2 Knockout Mice
Previous Article in Special Issue
Transcriptomic Analysis of the Chicken MDA5 Response Genes
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 11
Issue 4
10.3390/genes11040385
genes-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Knockout of IRF7 Highlights its Modulator Function of Host Response Against Avian Influenza Virus and the Involvement of MAPK and TOR Signaling Pathways in Chicken

by
Tae Hyun Kim
1,2,
Colin Kern
1 and
Huaijun Zhou
1,2,*
1
Department of Animal Science, University of California, Davis, CA 95616, USA
2
Integrative Genetics and Genomics Graduate Group, University of California, Davis, CA 95616, USA
*
Author to whom correspondence should be addressed.
Genes 2020, 11(4), 385; https://doi.org/10.3390/genes11040385
Submission received: 22 March 2020 / Accepted: 31 March 2020 / Published: 2 April 2020
(This article belongs to the Special Issue Avian Viruses and Host Responses)
Browse Figures
Review Reports Versions Notes

Abstract

:
Interferon regulatory factor 7 (IRF7) is known as the master transcription factor of the type I interferon response in mammalian species along with IRF3. Yet birds only have IRF7, while they are missing IRF3, with a smaller repertoire of immune-related genes, which leads to a distinctive immune response in chickens compared to in mammals. In order to understand the functional role of IRF7 in the regulation of the antiviral response against avian influenza virus in chickens, we generated IRF7-/- chicken embryonic fibroblast (DF-1) cell lines and respective controls (IRF7wt) by utilizing the CRISPR/Cas9 (clustered regularly interspaced short palindromic repeats/CRISPR-associated protein 9) system. IRF7 knockout resulted in increased viral titers of low pathogenic avian influenza viruses. Further RNA-sequencing performed on H6N2-infected IRF7-/- and IRF7wt cell lines revealed that the deletion of IRF7 resulted in the significant down-regulation of antiviral effectors and the differential expression of genes in the MAPK (mitogen-activated protein kinase) and mTOR (mechanistic target of rapamycin) signaling pathways. Dynamic gene expression profiling of the host response between the wildtype and IRF7 knockout revealed potential signaling pathways involving AP1 (activator protein 1), NF-κB (nuclear factor kappa B) and inflammatory cytokines that may complement chicken IRF7. Our findings in this study provide novel insights that have not been reported previously, and lay a solid foundation for enhancing our understanding of the host antiviral response against the avian influenza virus in chickens.

1. Introduction

The avian influenza virus (AIV) has significant impact on the poultry industry worldwide [1,2]. For instance, the highly pathogenic avian influenza (HPAI) outbreaks in the U.S. between 2014 and 2015 resulted in almost twelve and eight percent loss of layer chickens and meat turkeys, respectively (more than 50 million birds in total) [3,4]. In addition, AIVs are zoonotic and have posed great challenges to global public health [5,6]. However, current passive prevention practices such as quarantining and culling have been mainly applied to control AIV, partially due to our limited understanding of the antiviral response in birds compared to in mammals [7]. Improving our understanding of the avian host antiviral response could help to develop better strategies for controlling AIV in chickens and limiting zoonotic transmission, possibly by identifying potential biomarkers and therapeutic targets against AIV infection [8,9].
The coordination of Interferon Regulatory Factor 7 (IRF7) and IRF3 on the initiation of the type I interferon (IFN) response to establish the antiviral state is well characterized in mammals; however, chickens may utilize a different modulation system due to a smaller repertoire of immune-related genes in avian species [9,10]. Particularly, chickens are missing IRF3, and only have IRF7, which has less than 40% of amino acid sequence identity to its mammalian orthologues [11,12,13]. Therefore, studying the chicken IRF7 could be a great starting point to elucidate the underlying cellular mechanisms of host immune regulation against AIV in chickens.
In our previous study using stable IRF7 overexpression and knockdown in chicken DF-1 cell lines with double-stranded RNA (dsRNA) analog poly(I:C) (polyinosinic-polycytidylic acid) inductions, we demonstrated the conserved function of IRF7 as a type I IFN modulator [14]. Another study also suggested a similar role of IRF7 by the siRNA knockdown of IRF7 in chicken embryonic fibroblasts (CEFs), which limited IFNA, IFNB and STAT1 (signal transducer and activator of transcription 1) expression and increased Newcastle disease virus replication [15]. While these knockdown approaches such as RNA interference (RNAi) have successfully generated partial IRF7 loss-of-function phenotypes in chickens, and generated some novel insights into the role of IRF7 in the regulation of the host response to virus infection in poultry, RNAi has its own technical limitation of incompleteness of knockouts [16]. In this study, we decided to apply a CRISPR/Cas9 (clustered regularly interspaced short palindromic repeats/CRISPR-associated protein 9) approach to completely knockout (KO) the chicken IRF7 in the DF-1 cell line, in order to elucidate, in more depth, the underlying molecular mechanism of the modulation of the host response by IRF7 in chickens.
To understand the functional role of chicken IRF7 on the regulation of the host response to AIV at the genome-wide level, we applied a KO assay followed by a functional genomics approach (RNA-seq). First, the complete KO of IRF7 was achieved by the mutation chain reaction (MCR) which utilizes the CRISPR/Cas9 system to induce spontaneous homozygous mutation [17]. The viral replication phenotype was measured using two LPAIV (low pathogenic AIV) strains to examine the consequence of IRF7 deletion in vitro. Next, we applied RNA-seq to identify transcriptomic landscape differences behind the phenotype resulting from IRF7 KO and further reveal genes and signaling pathways associated with the regulatory function of chicken IRF7 in the context of AIV infection.

2. Materials and Methods

2.1. Construction of IRF7 Knockout MCR Plasmids

We constructed MCR plasmids targeting chicken IRF7 by modifying the pSpCas9(BB)-2A-Puro (px459, a gift from Feng Zhang-Addgene plasmid 62988) [18]. Three guide RNAs (gRNAs) targeting a region in exon 2 of IRF7 were designed, and homology arms (~0.5 kb each), in which the two genomic sequences are immediately adjacent to either side of the target double strand break site, were amplified using specific primers: HA1 forward 5′-cca gag agc agg gag ggg aa-3′, HA1 reverse 5′-gtc cgg gtc gat cca gca ga-3′, HA2 forward 5′-cgc cgt atc ttc cgc atc cc-3′, HA2 reverse 5′-agc ctc cct tga agc ccc tg-3′—and assembled with PciI/NotI double digested px459 by Gibson assembly [19]. Each gRNA was cloned into the HAs assembled px459 and three different MCR-IRF7-Exon2 vectors were generated (Figure 1A).

2.2. IRF7 Knockout Cell Clone Generation

Immortalized CEF, DF-1 cells (ATCC, Manassas, VA, USA) were cultured in Dulbecco’s Modified Eagle’s medium (Thermo Fisher Scientific, Waltham, MA, USA) supplemented with 10% fetal bovine serum (Thermo Fisher Scientific) and 1x Antibiotic-antimycotic (Thermo Fisher Scientific), and incubated at 37 °C in a humidified atmosphere containing 5% CO2. Each MCR-IRF7-Exon2 vector was transfected into DF-1 cells using the Lipofectamine 3000 reagent (Thermo Fisher Scientific) according to the manufacturer’s protocol. Approximately 5 × 103 cells were seeded into a 150 mm culture dish (Thermo Fisher Scientific) 48 h after transfection with puromycin (3 μg/mL) added to the 20 mL culture media. We selected 12 cell clones for each gRNA target by Pyrex® cloning cylinders (Corning, Corning, NY, USA) according to the manufacturer’s protocol after two weeks of transfection. Each clone was sub-cultured into a 48-well plate where 10% of cells were used for genotyping, and 90% of cells were cryopreserved, when each well reached confluency for subsequent analyses.

2.3. IRF7 Knockout Clone Genotyping

For genomic DNA extraction, we used 50 µL of the QuickExtract DNA extraction solution (Epicentre, San Diego, CA, USA) for each clone according to the manufacturer’s protocol. PCR analysis was performed to detect junctional bands corresponding to MCR insertion into the IRF7 exon2 locus, as well as the absence of a PCR band derived from the wild-type locus (Figure 1A, Table S1).

2.4. Cas9-Puro Cell Line Generation

To generate the piggyBac-Cas9-Puro plasmid (Figure 1A), the SpCas9-Puro cassette was cloned in from px459 into the modified PB513B plasmid backbone (System Biosciences, Mountain View, CA, USA). The PiggyBac transposase plasmid was co-transfected with the piggyBac-Cas9-Puro into DF-1 cells using the Lipofectamine 3000 reagent (Thermo Fisher Scientific) to induce efficient and stable integration. Puromycin (3 μg/mL) was added to the culture media 48 h after transfection, and stably integrated cell lines were selected for two weeks. Cas9 activity was examined by multiplex gRNA transfection followed by genomic DNA PCR. Two gRNAs targeting the exon 2 and exon 10 loci of IRF7 were co-transfected into the Cas9-Puro cell line, and the deletions of the genomic DNA region between gRNA targets were examined by genomic DNA PCR (Figure S1, Table S1).

2.5. Quantitative Reverse Transcriptase PCR

Total RNA was isolated from approximately 1 million cells using the Direct-zol RNA MiniPrep Kit (Zymo research, Irvine, CA), and complement DNA (cDNA) was synthesized from total RNA (500 ng) using Verso cDNA Synthesis Kit (Thermo Fisher Scientific). Quantitative reverse transcriptase PCR (qRT-PCR) was performed using the QuantStudio 3 Real-Time PCR system (Thermo Fisher Scientific) with the SYBR Select Master Mix (Life Technologies). IRF7, IFNA and IFNB expressions were normalized to that of the GAPDH (glyceraldehyde 3-phosphaste dehydrogenase) gene using the ΔΔCT method (Table S1) [14,20].

2.6. AIV and In Vitro Infection

A/Chicken/California/2000 (H6N2) and A/Chicken/California/1999 (H10N7) low pathogenic avian influenza virus (LPAIV) strains were kindly provided by Dr. Rodrigo Gallardo (University of California, Davis) and Dr. Peter Woolcock (University of California, Davis, California Animal Health & Food Safety (CAHFS)), respectively. Each LPAIV was propagated in Madin-Darby Canine Kidney (MDCK) cells as described previously [21]. All in vitro AIV infections were performed using CellBIND 12 well tissue culture plates (Corning, Corning, NY, USA) with 1 × 106 cells per well at the seeding. For the viral replication kinetics phenotype, established DF-1 cell lines were infected with either H6N2 or H10N7 at a multiplicity of infection (MOI) of 0.01 with 0.05 μg/mL TPCK-trypsin in DMEM. Culture supernatants were collected at 0, 12, 24 and 36 hours post-infection (hpi) and the viral titer of each sample was measured by the plaque assay using MDCK cells [22]. For the transcriptome profiling study, cell lines were infected with either mock or H6N2 at an MOI of 1 with 0.05 μg/mL TPCK-trypsin in DMEM. Trizol reagent (Thermo Fisher Scientific) was directly added to the cell monolayer, after it was washed twice with PBS, to extract the total RNA at 0, 2, 6 and 12 hpi.

2.7. RNA Sequencing and Data Analysis

A total of 24 cDNA libraries were prepared from three biological replicates of each time point (0 h (mock), 2-, 6- and 12 hpi IRF7wt or IRF7-/-). Directional RNA sequencing libraries were prepared from poly-adenylated RNA and sequenced with the Illumina HiSeq4000, which generated over 20 million 100 bp paired-end reads per sample, except for 2 hpi samples where 150 bp paired-end reads were generated. The read files from RNA-seq analysis have been deposited in NCBI’s Gene Expression Omnibus with accession number GSE119590. We checked the quality of each library using fastQC (Babraham Institute, Cambridge, UK, version 0.11.6) and trimmed the adaptor sequence with TrimGalore (Babraham Institute, Cambridge, UK, version 0.4.5). We aligned the trimmed fastq files to the galGal5 chicken genome using the STAR aligner with NCBI annotation release 103 [23]. Unmapped reads were aligned against the H6N2 (A/chicken/CA/6643/2001) genome [24]. Raw read counts were extracted by HTSeq from each aligned bam file and used to identify differentially expressed genes (DEGs) [25]. Both DESeq2 and EdgeR R packages were used to identify DEGs and the DEG sets from both packages were combined (false discovery rate (FDR) < 5% or 1% in any one of the packages) [26,27,28]. The combined DEG lists were further filtered by removing the low expression genes if both samples of a given comparison had fragments per kilobase per million reads (FPKM) value less than 1. Functional annotations for significantly differentially expressed genes were performed using DAVID 6.8 [29,30]. The enriched gene ontology (GO) terms on biological processes and the pathways obtained from DAVID functional analysis were filtered for significance by gene count ≥5 and p-value < 0.05. Additional pathway analyses were performed using the Ingenuity Pathway Analysis (IPA) software (QIAGEN, Germantown, MD, USA). RNA expression of the MCR target region was visualized by the Integrative Genomics Viewer (IGV, version 2.4.13) [31]. Heat maps and k-means clustering analysis diagrams were analyzed using in-house R and Python scripts. We applied the gap statistics and sum of squared error methods to evaluate the optimum number of clusters [32].

3. Results

3.1. Establishment of IRF7 Knockout DF-1 Cell Clones

We applied the MCR method, which uses homology directed repair (HDR), to disrupt chicken IRF7 in a highly efficient manner. Three gRNAs targeting within a 6bp region of IRF7 exon 2 were designed to address the potential off-target effects induced by the CRISPR/Cas9 system, since the MCR insert constitutively expresses both SpCas9 and gRNA (Figure 1A) [33]. We were able to use an identical set of homology arms (HAs) flanking the gRNA-SpCas9-Puro transgenes for all the gRNAs, as the target sites of the gRNAs were adjacent to each other (Figure 1A). PCR analysis of the IRF7 exon 2 locus in individual cell clones confirmed the precise gRNA- and HDR-directed genomic insertion of the MCR construct in the IRF7 exon 2 locus (Figure 1B). By the MCR and puromycin selection, we were able to achieve 50% of bi-allele KO cell clones, which only required genotyping less than 10 clones per gRNA. The genotypes we observed were either bi-allele genomic insertions or single-allele insertions. We randomly selected one clone for each gRNA target, which served as biological replicates for IRF7-/- (Figure 1B, marked asterisk). We also confirmed disrupted IRF7 transcription and interrupted mRNA splicing in the IRF7-/- cells at the MCR region from RNA-seq analysis (Figure 1C).
Since the IRF7-/- cells have a targeted insertion of SpCas9 and puromycin resistance genes in their genome, we generated corresponding control IRF7wt cells that have SpCas9 and puromycin resistant genes randomly integrated by piggyBac transposition (Figure 1A, right panel) [34]. We assessed the activity of SpCas9 and puromycin resistance genes in the IRF7wt cells by co-transfecting two gRNAs and selecting with puromycin, respectively. SpCas9 activity was confirmed by successfully detecting the 2.7kb-deletion alleles as a result of gRNA multiplexing (Figure S1).

3.2. LPAIV In Vitro Infection

We next sought to determine whether IRF7 deficiency affects the viral replication phenotype in vitro. IRF7-/- and IRF7wt cells were infected with H6N2 at an MOI of 0.01 and then supernatant was collected at 0, 12, 24 and 36 hpi to measure viral titers (Figure 1D). We observed significantly an increased H6N2 viral replication level in the IRF7-/- cells compared with in the IRF7wt cells at 12 hpi. Another strain of LPAIV H10N7 showed a higher viral replication level as a result of IRF7 knockout (Figure S2A). Transcriptome analysis corresponds to this pattern where we observed a relatively lower fraction of host mRNA and a higher fraction of viral RNA in IRF7-/- cells than in IRF7wt cells at 6 and 12 hpi, yet not at a significant level (Figure S2B).

3.3. The Effect of IRF7 Knockout on the Host Response Against H6N2 Infection

To elucidate the molecular mechanisms of chicken IRF7 in the context of AIV infection in depth, we directly contrasted the transcriptomes between the IRF7wt and IRF7-/- cells at each time point. Differential gene-expression analysis indicated that 180 genes were significantly upregulated, and 164 genes were significantly downregulated in IRF7-/- cells compared to their expressions in IRF7wt cells without infection (FDR < 5%; Figure 2A and Table S2). Upon infection, there were 138 down- and 100 up-regulated DEGs at 2 hpi due to IRF7 deletion, and the number of DEGs increased to 441 down- and 483 up-regulated genes at 6 hpi. At 12 hpi, fewer genes were differentially expressed (366 up- and 300 down-regulated) compared to at 6 hpi (Figure S2). Of note, the numbers of DEGs observed comparing IRF7wt with IRF7-/- cells were similar between up-regulated genes and down-regulated genes at all time points. We used these DEGs to perform functional enrichment analyses by DAVID to identify enriched biological process gene ontology (GO) functions and pathways that were regulated by IRF7 (Figure 2B). Upon infection at 2 hpi, we observed enriched immune response related terms such as antigen processing and presentation, as well as the regulation of the inflammatory response. The mTOR (mechanistic target of rapamycin) signaling GO term and its downstream process actin polymerization were also enriched at 6 hpi. Protein processing associated terms were significantly enriched across all time points (Figure 2B). Heatmaps with fold changes of genes in the influenza A pathway and significantly enriched pathways across all time points using all individual DEGs are presented in Figure 2C–E. IL18 (interleukin 18), MAPK13 (mitogen-activated protein kinase 13; p38) and CHUK (conserved helix-loop-helix ubiquitous kinase; IKKA) were up-regulated in influenza A and the MAPK signaling pathways, whereas MHC (major histocompatibility complex) protein family genes (BF1, BF2, DMA, DMB2 and BLB2) displayed down-regulated expression across all time points once IRF7 was deleted.

3.4. Comparative Analysis of DEGs Between IRF7 Knockout and Overexpression

Our previous work has identified candidate genes that were regulated by IRF7 in the context of H6N2 infection by applying an inducible overexpression-based gain-of-function approach [35]. We infected IRF7 overexpression cells (CuO-IRF7: 2-fold IRF7 overexpression) and their respective controls with H6N2 (1 MOI) and performed RNA-seq at 0 h (mock) and 6 hpi. To identify candidate genes regulated by IRF7, we directly compared the datasets of IRF7 KO and overexpression that were contrasted to their respective controls at 0h and 6 hpi (IRF7-/-/IRF7wt vs. CuO-IRF7/Control; FDR < 5%, Figure 3). A total of 245 common DEGs between the two datasets were identified, which were grouped into four distinct clusters: same direction of regulations (upregulated or downregulated) or opposite directions of the regulation between knockout and overexpression. For example, cell adhesion molecules (CAMs) such as laminins, integrins and collagens had up-regulated expression if IRF7 expression was altered in any direction. On the other hand, genes enriched for the cellular response to laminar fluid shear stress or small GTPase-mediated signal transduction were down-regulated for both IRF7 KO and overexpression. Of particular interest, opposite directions of regulation were observed. TOR-signaling and I-kappaB phosphorylation-associated DEGs were up-regulated due to IRF7 deletion and down-regulated once IRF7 was overexpressed. In addition, the expressions of antigen processing and presentation-related genes and immune response genes were positively correlated with IRF7 expression. Interestingly, the MHC I family gene BF2 was down-regulated in both IRF7 KO and overexpression, whereas the expression of MHC II family genes (DMA, DMB2, BLB2) corresponded with the IRF7 expression level.

3.5. Differential Expression Analysis of the Host Response upon H6N2 Infection

Next, we compared the host responses between the cell lines upon H6N2 infection, by comparing infected cells to the uninfected cells at each time point of infection, and identified which genes and pathways were potentially modulated by IRF7 in the host response (Figure 4A, Table S3). To further analyze which biological process was possibly suppressed due to IRF7 deletion, or alternatively activated in the antiviral response to compensate for the absence of IRF7, we first compared the list of DEGs between cell lines at each time point (Figure 4B). Next, we took the DEGs that were specific to each cell line and performed GO enrichment analysis at each time point (Figure 4C). Enriched terms demonstrated IRF7 to be involved in a wide range of biological processes including mitosis, cell cycle, protein modification, transcription and translation. At 6 hpi, genes involved in NF-κB (nuclear factor-κB) signaling regulation were enriched only in the IRF7 KO cells. However, we were not able to specify which functions are mainly regulated by IRF7 with this analysis, mainly because the expression levels of individual genes did not show significant differences between the cell lines, despite the differences in the DEG lists and functional annotations.

3.6. Dynamic Gene Expression Profile Analysis Revealed an Alternative Pathway to IRF7 Regulation

To evaluate the global gene expression profile alteration due to IRF7 KO in the host response to H6N2 infection, we performed k-means clustering among the 5091 DEGs that were obtained from previous analyses in Figure 2 and Figure 4. A total of six clusters were identified based on their expression patterns (Figure 5A, Table S4). A predominant expression pattern shift was observed from Cluster V to Cluster II due to the deletion of IRF7, which was almost 60% (630 out of 1067, Figure 5B, Table S4) of the genes from Cluster V. Gene ontology based on the IRF7wt DEGs of each cluster revealed that the genes in Cluster V have enriched functions associated with translation, inflammatory response and immune response (Figure 5C). These 630 shifted genes were similarly enriched with GO terms such as translation, immune response and defense response to pathogens (Table S4). Moreover, the expression dynamics of the genes changed from the stable-up-down pattern to constant up pattern once IRF7 was deleted. Based on this observation, we hypothesized that the deletion of IRF7 resulted in an insufficient host response against AIV, hence the expression of alternative inflammatory response and immune response genes might be needed to compensate for the effect of deleted IRF7. This hypothesis was demonstrated in the influenza A pathway (Figure 5D, Figure S3) in which seven out of twelve genes from Cluster V in the IRF7wt cells were shifted to Cluster II in IRF7-/- cells. These genes include key immune and inflammatory response genes, such as inflammatory cytokines (IL12 (interleukin 12) and transcription factors (FOS (Fos proto-oncogene; AP1 transcription factor subunit) and RELA (RELA proto-oncogene; NF-κB subunit, p65)), as well as their upstream molecules (IKBKE (inhibitor of kappa light polypeptide gene enhancer in B-cells, kinase epsilon) and NFKBIB (NF-κB inhibitor beta)).

4. Discussion

IRF7 was identified as the regulator of the type I IFN response in an Irf7-/- mouse model, in which the Irf7-/- mouse was significantly more vulnerable to viral infections, with markedly inhibited type I IFN expression, than the wild type [10,36]. The significant role of IRF7 in the regulation of the host anti-viral response was further elucidated in humans, in which the IRF7-deficient patient had impaired type I and III IFN responses against primary infection and experienced a life-threatening influenza infection [37].
Interaction of highly homologous IRF3 and IRF7 (oligomerization) is recognized as a key regulatory mechanism of the transcriptional activation of type I IFNs in initiating the host antiviral state in mammals [38]. Chicken IRF7 was first identified and named as IRF3 due to its high degree of similarity to mammalian IRF3, but was renamed as IRF7 when subsequent studies revealed that IRF3 is absent in avian species [39,40]. We took a de novo transcript assembly approach to identify potential chicken IRF3 at the transcript level, in order to eliminate potential annotation issues in the chicken reference assembly, yet we did not detect the chicken IRF3 [14]. Based on the function of IRF7 homologs in mammals and missing a member of the IRF3/IRF7 dyad, we hypothesized that chicken IRF7 is a central transcription factor in the antiviral response against AIV infection in chickens.
In this study, we employed a functional genomic approach to elucidate the functional role of chicken IRF7 in the host response against AIV infection in chickens. Our work has first revealed that one of type I IFN, IFNB, but not IFNA, was selectively responsive to IRF7 KO in chicken DF-1 cells. The type I IFN induction models from in vitro studies using mouse embryonic fibroblast or primary tracheal epithelial cells of Irf7-/- mice showed an abolished production of type I IFNs upon viral infection [10,41]. In our previous studies, we consistently observed a significant upregulation of IFNB expression as a result of IRF7 overexpression in the DF-1 cell lines upon dsRNA induction or H6N2 infection [14,35]. On the other hand, IFNA mRNA levels were not significantly affected by the level of IRF7 expression [14,35], and a similar pattern was observed in this study (Figure 1E–G). Although DF-1 cells are capable of mounting a type I IFN response, the level of response is possibly lower than that of primary CEFs due to a higher suppressor of cytokine signaling 1 (SOCS1) expression in DF-1 [42]. Therefore, additional studies in different cell types such as lymphocytes, macrophages or dendritic cells could further address the context and cell type specific type I IFN regulation by IRF7.
We also demonstrated that the deletion of IRF7 resulted in the higher replication of LPAIVs in vitro, which further suggests an anti-viral role of IRF7 in chickens. This phenotype was in agreement with the susceptible phenotype observed in Irf7-/- mice against viral infections (herpes simplex virus, encephalomyocarditis virus or influenza virus) [10,43]. Functional enrichment of DEGs contrasting between IRF7wt and IRF7-/- cells identified the MAPK pathway as one of the significant pathways altered due to IRF7 KO (Figure 2D). The MAPK pathway is an important signaling cascade that controls a wide range of cellular responses such as proliferation, differentiation, and the immune response [44]. Among the known four MAPK family members, p38 MAPK (MAPK13) was upregulated as a result of IRF7 KO. The inflammatory response due to activation of p38 by influenza virus infection has been reported, and the pro-viral function of p38 MAPK was confirmed in a mouse model where the inhibition of p38 MAPK significantly protected mice from lethal influenza infection (H5N1) [45]. In this perspective, the upregulation of the MAPK pathway in the IRF7-/- cells could possibly explain the increased viral replication phenotype of the IRF7 KO cells. Our observations were in line with the study comparing transcriptomes between the wildtype and Irf3-/-Irf7-/- double KO (DKO) mice infected with PR8 (H1N1), in which the upregulation of MAPK signaling pathway was observed to correlate with susceptible phenotypes of mutants such as significant weight loss, a higher mortality and a higher viral load with infection [43]. A validation of the chicken IRF7 KO susceptibility phenotype against AIV at a systemic level is warranted to elucidate the underlying molecular mechanism.
The integration of both gain- and loss-of-function approaches employed in this study provided a unique opportunity to distinguish between the genes that were either modulated by IRF7 or altered as a result of IRF7 perturbation. For instance, we observed differential expression of genes that are involved in the membrane trafficking pathways, such as endocytosis, lysosomes and phagosomes, as a result of IRF7 KO after infection, as well as at baseline (Figure 2). Of particular note, both MHC I and MHC II family genes were downregulated due to IRF7 deletion. Further analysis comparing the DEGs of IRF7 KO and overexpression revealed that only MHC II family genes were positively correlated to IRF7 expression, while the expressions of MHC I genes were downregulated in both IRF7 KO and overexpression (Figure 3). This suggest an important role of chicken IRF7 in the regulation of antigen presentation against viral infection through the MHC II family of genes. Cell adhesion molecules, particularly integrins, play important roles as mediators in surveillance and the regulation of immune cell functions [46], and we observed the up-regulation of integrins as a result of altered IRF7 expressions in any direction, even in the absence of infection. This could imply that the IRF7-dependent host response in chicken is a finely-tuned cell process, and IRF7 may be an important regulator of cell state determination between the homeostatic state and antiviral microenvironments in the cells throughout the host-pathogen interaction.
In addition, a negative correlation between mTOR and IRF7 was corroborated: IRF7 KO resulted in the upregulation of MTOR, while the overexpression of IRF7 led to the downregulation of MTOR (Figure 3, Figure S4). The mTOR signaling pathway is highly conserved in eukaryotes from yeast to humans and plays a central regulator role in cellular activities such as cell proliferation, metabolism, survival and differentiation [47,48]. This may suggest chicken IRF7 directly regulates mTOR signaling to modulate cellular processes in the host response against AIV. This is consistent with the previous report that indicated the repressive translational regulation of Irf7 mRNA by the mTOR signaling pathway [49]. The mTOR signaling pathway is triggered during the host innate immune response and reconfigures cellular metabolism to provide the cellular environment for a subsequent immune response in a cell-type-specific manner [47,50,51]. The mTOR signaling pathway coordinates metabolic processes through forming two distinct multi-protein complexes, mTOR complex 1 (mTORC1) and mTORC2 [47]. Our functional enrichment analyses of IRF7 KO revealed alterations of pathways usually mediated by mTORC1 such as the metabolic process, protein synthesis (translation and ribosome) and protein turnover (protein ubiquitination) within the host response (Figure 3 and Figure 5). Moreover, IRF7 deletion resulted in the upregulation of RICTOR (rapamycin insensitive companion of mTOR), a component of mTORC2, which functions as an effector of PI3K/Akt signaling and cytoskeleton rearrangement [47]. Accordingly, we observed an expression pattern shift of PI3K and upregulation of AKT3 in the IRF7-/- cells compared to in the IRF7wt (Figure S3). Our previous study provided another line of evidence supporting the mTOR-mediated IRF7 modulation of the host response, in which GO terms such as arginine and proline metabolism, pyrimidine metabolism and the regulation of the actin cytoskeleton were significantly enriched from the DEGs of IRF7 overexpression vs. wildtype [35].
Comparing the host responses between the IRF7wt and IRF7-/- cells enabled the uncovering of the pathways that work together with IRF7 in the host antiviral response against H6N2 infection. Significantly more genes (1200+) differentially expressed between the IRF7-/- cells and IRF7wt cells at 12 hpi (Figure 4A) highlighted the potential pathways that might compensate for the host response, due to the lack of IRF7, in response to AIV. Dynamic gene expression profiling applying k-means clustering further revealed a set of genes that had their expression pattern shifted due to IRF7 deletion with AIV infection (Figure 5). The shifts of immune response- and inflammatory response-related genes again emphasized a significant modulator role of IRF7. Particularly, the constant augmentation of transcription factor NF-κB (nuclear factor kappa B; subunit RELA, inhibitor NFKBIB) and AP1 (activator protein 1; subunit FOS) expression during infection due to IRF7 deletion suggests that these factors complement IRF7 in the host antiviral response by modulating inflammatory cytokines and type I IFNs in chickens, as in mammals (Figure 5D, Figure S3) [52]. This is in line with the observation from a Irf3-/-Irf7-/- (DKO) mice study in which the DKO mice displayed impaired interferon activation and an increased inflammatory response upon influenza A virus infection [43]. In addition, Irf3-/-Irf7-/- DKO mice as well as Irf3-/-Irf5-/-Irf7-/- triple KO mice infected with dengue virus had demonstrated enriched inflammatory process-related genes at 24 hpi compared to the wildtype animals [53].
A schematic model of early host response associated with chicken IRF7 is proposed in Figure 6, based on the current study and our previous in vitro studies [14,35]. The entry of AIV is recognized by pattern recognition receptors (PRRs) sensing viral moieties. TLR3 (toll-like receptor 3) and MDA5 (melanoma differentiation-associated protein 5, IFIH1) respond to dsRNA, and TLR7 responds to single-stranded RNA (ssRNA) during AIV infection in chickens [54,55,56]. This leads to a subsequent signaling cascade via adaptor molecules and protein kinases that activate the key transcription factors (AP1, NF-κB and IRF7). The MAPK and PI3K/Akt pathways are potentially engaged in this signaling module in addition to the classical type I IFN induction pathway activating IRF7 by IkB kinases (IKK) TBK1 (TANK binding kinase 1), IKKε (IKBKE) or IKKα (CHUK) [13,57,58]. Overall, our data argue that IRF7 is a major modulator of the host antiviral response, and that the transcription factors NF-κB and AP1 are also capable of regulating both IFNA and INFB in addition to the inflammatory cytokines. Among the wide range of biological processes in the antiviral responses that are modulated by IRF7, two modes of regulatory mechanism were speculated: a cascade of downstream type I IFN responses, and the direct modulation of genes by IRF7 binding to their regulatory regions (Figure 6). Genome-wide chromatin immuno-precipiation (ChIP)-seq approach using CRISPR-Cas9-mediated knock-in could further dissect the molecular basis of the IRF7-mediated antiviral response.
In summary, we integrated genome editing and RNA sequencing to reveal some novel findings on the role of chicken IRF7 in the context of AIV infection at the whole transcriptome level. This study suggests the conserved function of chicken IRF7 as a major modulator of the antiviral response, despite the considerable differences between the chicken and mammalian innate immune systems. Moreover, our results suggest a novel insight on the chicken type I IFN response, in which MAPK and mTOR signaling may play important roles in modulating the antiviral response against AIV. Lastly, we found evidence that implicates the involvement of other important mediators working together with IRF7 in the host antiviral response in chickens. Although further systemic-level studies using different cell types are warranted, our results have provided the first line of evidence that has enhanced our understanding of the role of IRF7 on the host antiviral response to AIV in chickens.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4425/11/4/385/s1, Figure S1 (A) piggyBac SpCas9-puro plasmid construct. (B) Locations of gRNAs and PCR primers used for analysis of the genomic deletion are indicated on the IRF7 locus. (C) Genomic DNA PCR analysis showing the SpCas9 mediated large deletion of IRF7 allele (DEL) in IRF7wt cells transfected with two gRNAs targeting each end of IRF7 CDS (lane 2: 3.7 kb between DEL-F and DEL-R–2.7kb deletion between two gRNAs = 1 kb amplicon). Short (30 sec) extension time limited the amplification of 3.7 kb un-deleted allele. Un-deleted alleles are detected in all groups as analyzed in the population (INT). Figure S2 (A) Culture supernatants were collected at each time point and viral titer was determined by plaque assay as (plaque forming unit (PFU)+1)/mL using Madin-Darby kidney (MDCK) cells. (B) Alignment rates of RNA-seq libraries against reference chicken genome (galGal5) and H6N2 genome (C) Venn diagram of DEGs between the comparisons of IRF7wt and IRF7-/- at all time points post infection (0,2,6,12 hpi). (D) Heatmap showing relative gene expression (log2(FPKM+1)) of 63 common DEGs reflected in (C) with all samples. Each row represents a DEG and each column represents a sample. Figure S3 KEGG pathway view for differentially expressed genes (DEGs) in Influenza A pathway. DEGs of IRF7 knockout (vs. IRF7wt), H6N2 infection (vs. 0 h) and expression cluster shift from IRF7wt to IRF7-/- are indicated in slightly modified ‘Influenza A pathway: gga05164’. IRF7 knockout DEGs are indicated as blue (down-) or red (upregulation) if the gene is called as DEG at any time point. H6N2 infection DEGs are indicated with asterisk (*) regardless of cell line. Expression cluster shift is indicated with green line. Figure S4 DEGs of (A) IRF7 overexpression or (B) knockout at 6 hpi in mTOR signaling pathway using Ingenuity Pathway Analysis (IPA) software. In silico molecular activity predictor function of IPA was applied based on the expression dataset to visualize the overall effect on the pathway Table S1. Primers used in this study Table S2. Differentially expressed genes between of IRF7-/- and control (IRF7wt) DF-1 cell lines at each time point Table S3. Differentially expressed genes upon H6N2 infection in each cell line Table S4. 1. Cluster shift of Differentially expressed genes due to IRF7 deletion, 2. Enriched gene ontology terms (Biological Process) of cluster shifted genes from 5 to 2.

Author Contributions

Conceptualization, T.H.K. and H.Z.; Formal analysis, T.H.K. and C.K.; Funding acquisition, T.H.K. and H.Z.; Investigation, T.H.K.; Project administration, H.Z.; Resources, H.Z.; Writing–original draft, T.H.K., C.K. and H.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by United States Department of Agriculture (USDA), National Institute of Food and Agriculture (NIFA), Multistate Research Project NRSP8 and NC1170 (HZ), and the California Agricultural Experimental Station (HZ). TK was partially supported by Fellowships Grant Program Award #2017-67011-26762 from the USDA NIFA.

Acknowledgments

Two LPAIV strains used in this study were gifts from the labs of Rodrigo Gallardo and Peter Woolcock.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Olsen, B.; Munster, V.J.; Wallensten, A.; Waldenstrom, J.; Osterhaus, A.D.; Fouchier, R.A. Global patterns of influenza a virus in wild birds. Science 2006, 312, 384–388. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Swayne, D.E. Impact of vaccines and vaccination on global control of avian influenza. Avian Dis. 2012, 56, 818–828. [Google Scholar] [CrossRef] [PubMed]
  3. Jhung, M.A.; Nelson, D.I.; Centers for Disease Control and Prevention. Outbreaks of avian influenza A (H5N2), (H5N8), and (H5N1) among birds—United States, December 2014–January 2015. MMWR Morb. Mortal. Wkly. Rep. 2015, 64, 111. [Google Scholar] [PubMed]
  4. Ramos, S.; MacLachlan, M.; Melton, A. Impacts of the 2014-2015 Highly Pathogenic Avian Influenza Outbreak on the U.S. Poultry Sector; US Department of Agriculture World Agriculture Outlook Board: Washington, DC, USA, 2017.
  5. Wang, X.; Jiang, H.; Wu, P.; Uyeki, T.M.; Feng, L.; Lai, S.; Wang, L.; Huo, X.; Xu, K.; Chen, E.; et al. Epidemiology of avian influenza A H7N9 virus in human beings across five epidemics in mainland China, 2013–2017: An epidemiological study of laboratory-confirmed case series. Lancet Infect. Dis. 2017, 17, 822–832. [Google Scholar] [CrossRef]
  6. Lai, S.; Qin, Y.; Cowling, B.J.; Ren, X.; Wardrop, N.A.; Gilbert, M.; Tsang, T.K.; Wu, P.; Feng, L.; Jiang, H.; et al. Global epidemiology of avian influenza A H5N1 virus infection in humans, 1997–2015: A systematic review of individual case data. Lancet Infect. Dis. 2016, 16, e108–e118. [Google Scholar] [CrossRef] [Green Version]
  7. Goossens, K.E.; Ward, A.C.; Lowenthal, J.W.; Bean, A.G. Chicken interferons, their receptors and interferon-stimulated genes. Dev. Comp. Immunol. 2013, 41, 370–376. [Google Scholar] [CrossRef] [PubMed]
  8. Downing, T.; Cormican, P.; O‘Farrelly, C.; Bradley, D.G.; Lloyd, A.T. Evidence of the adaptive evolution of immune genes in chicken. BMC Res. Notes 2009, 2, 254. [Google Scholar] [CrossRef] [Green Version]
  9. Magor, K.E.; Miranzo Navarro, D.; Barber, M.R.; Petkau, K.; Fleming-Canepa, X.; Blyth, G.A.; Blaine, A.H. Defense genes missing from the flight division. Dev. Comp. Immunol. 2013, 41, 377–388. [Google Scholar] [CrossRef]
  10. Honda, K.; Yanai, H.; Negishi, H.; Asagiri, M.; Sato, M.; Mizutani, T.; Shimada, N.; Ohba, Y.; Takaoka, A.; Yoshida, N.; et al. IRF-7 is the master regulator of type-I interferon-dependent immune responses. Nature 2005, 434, 772–777. [Google Scholar] [CrossRef]
  11. Huang, B.; Qi, Z.T.; Xu, Z.; Nie, P. Global characterization of interferon regulatory factor (IRF) genes in vertebrates: Glimpse of the diversification in evolution. BMC Immunol. 2010, 11, 22. [Google Scholar] [CrossRef] [Green Version]
  12. Santhakumar, D.; Rubbenstroth, D.; Martinez-Sobrido, L.; Munir, M. Avian Interferons and Their Antiviral Effectors. Front. Immunol. 2017, 8, 49. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Cheng, Y.; Zhu, W.; Ding, C.; Niu, Q.; Wang, H.; Yan, Y.; Sun, J. IRF7 Is Involved in Both STING and MAVS Mediating IFN-beta Signaling in IRF3-Lacking Chickens. J. Immunol. 2019, 203, 1930–1942. [Google Scholar] [CrossRef] [PubMed]
  14. Kim, T.H.; Zhou, H. Functional Analysis of Chicken IRF7 in Response to dsRNA Analog Poly(I:C) by Integrating Overexpression and Knockdown. PLoS ONE 2015, 10, e0133450. [Google Scholar] [CrossRef]
  15. Wang, Z.; Ning, Z.; Sun, M.; Gao, S.; Kang, Y.; Xie, P.; Ren, T. Interferon regulatory factor 7- (IRF7-) mediated immune response affects Newcastle disease virus replication in chicken embryo fibroblasts. Acta Vet. Hung. 2014, 62, 500–511. [Google Scholar] [CrossRef] [Green Version]
  16. Boutros, M.; Ahringer, J. The art and design of genetic screens: RNA interference. Nat. Rev. Genet. 2008, 9, 554–566. [Google Scholar] [CrossRef]
  17. Gantz, V.M.; Bier, E. Genome editing. The mutagenic chain reaction: A method for converting heterozygous to homozygous mutations. Science 2015, 348, 442–444. [Google Scholar] [CrossRef] [Green Version]
  18. Ran, F.A.; Hsu, P.D.; Wright, J.; Agarwala, V.; Scott, D.A.; Zhang, F. Genome engineering using the CRISPR-Cas9 system. Nat. Protoc. 2013, 8, 2281–2308. [Google Scholar] [CrossRef] [Green Version]
  19. Gibson, D.G.; Young, L.; Chuang, R.Y.; Venter, J.C.; Hutchison, C.A., 3rd; Smith, H.O. Enzymatic assembly of DNA molecules up to several hundred kilobases. Nat. Methods 2009, 6, 343–345. [Google Scholar] [CrossRef]
  20. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2(∆∆C(T)) Method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef]
  21. Eisfeld, A.J.; Neumann, G.; Kawaoka, Y. Influenza A virus isolation, culture and identification. Nat. Protoc. 2014, 9, 2663–2681. [Google Scholar] [CrossRef] [Green Version]
  22. Huprikar, J.; Rabinowitz, S. A simplified plaque assay for influenza viruses in Madin-Darby kidney (MDCK) cells. J. Virol. Methods 1980, 1, 117–120. [Google Scholar] [CrossRef]
  23. Dobin, A.; Davis, C.A.; Schlesinger, F.; Drenkow, J.; Zaleski, C.; Jha, S.; Batut, P.; Chaisson, M.; Gingeras, T.R. STAR: Ultrafast universal RNA-seq aligner. Bioinformatics 2013, 29, 15–21. [Google Scholar] [CrossRef] [PubMed]
  24. Webby, R.J.; Woolcock, P.R.; Krauss, S.L.; Webster, R.G. Reassortment and interspecies transmission of North American H6N2 influenza viruses. Virology 2002, 295, 44–53. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Anders, S.; Pyl, P.T.; Huber, W. HTSeq—A Python framework to work with high-throughput sequencing data. Bioinformatics 2015, 31, 166–169. [Google Scholar] [CrossRef] [PubMed]
  26. Love, M.I.; Huber, W.; Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014, 15, 550. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Robinson, M.D.; McCarthy, D.J.; Smyth, G.K. edgeR: A Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 2010, 26, 139–140. [Google Scholar] [CrossRef] [Green Version]
  28. McCarthy, D.J.; Chen, Y.; Smyth, G.K. Differential expression analysis of multifactor RNA-Seq experiments with respect to biological variation. Nucleic Acids Res. 2012, 40, 4288–4297. [Google Scholar] [CrossRef] [Green Version]
  29. Dennis, G., Jr.; Sherman, B.T.; Hosack, D.A.; Yang, J.; Gao, W.; Lane, H.C.; Lempicki, R.A. DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biol. 2003, 4, P3. [Google Scholar] [CrossRef]
  30. Huang, D.W.; Sherman, B.T.; Lempicki, R.A. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat. Protoc. 2009, 4, 44–57. [Google Scholar] [CrossRef]
  31. Robinson, J.T.; Thorvaldsdottir, H.; Winckler, W.; Guttman, M.; Lander, E.S.; Getz, G.; Mesirov, J.P. Integrative genomics viewer. Nat. Biotechnol. 2011, 29, 24–26. [Google Scholar] [CrossRef] [Green Version]
  32. Tibshirani, R.; Walther, G.; Hastie, T. Estimating the number of clusters in a data set via the gap statistic. J. Royal Stat. Soc. Ser. B Stat. Methodol. 2001, 63, 411–423. [Google Scholar] [CrossRef]
  33. Fu, Y.; Foden, J.A.; Khayter, C.; Maeder, M.L.; Reyon, D.; Joung, J.K.; Sander, J.D. High-frequency off-target mutagenesis induced by CRISPR-Cas nucleases in human cells. Nat. Biotechnol. 2013, 31, 822–826. [Google Scholar] [CrossRef] [Green Version]
  34. Wu, S.C.; Meir, Y.J.; Coates, C.J.; Handler, A.M.; Pelczar, P.; Moisyadi, S.; Kaminski, J.M. piggyBac is a flexible and highly active transposon as compared to sleeping beauty, Tol2, and Mos1 in mammalian cells. Proc. Natl. Acad. Sci. USA 2006, 103, 15008–15013. [Google Scholar] [CrossRef] [Green Version]
  35. Kim, T.H.; Zhou, H. Overexpression of Chicken IRF7 Increased Viral Replication and Programmed Cell Death to the Avian Influenza Virus Infection Through TGF-Beta/FoxO Signaling Axis in DF-1. Front. Genet. 2018, 9. [Google Scholar] [CrossRef]
  36. Wilk, E.; Pandey, A.K.; Leist, S.R.; Hatesuer, B.; Preusse, M.; Pommerenke, C.; Wang, J.; Schughart, K. RNAseq expression analysis of resistant and susceptible mice after influenza A virus infection identifies novel genes associated with virus replication and important for host resistance to infection. BMC Genom. 2015, 16, 655. [Google Scholar] [CrossRef]
  37. Ciancanelli, M.J.; Huang, S.X.; Luthra, P.; Garner, H.; Itan, Y.; Volpi, S.; Lafaille, F.G.; Trouillet, C.; Schmolke, M.; Albrecht, R.A.; et al. Infectious disease. Life-threatening influenza and impaired interferon amplification in human IRF7 deficiency. Science 2015, 348, 448–453. [Google Scholar] [CrossRef] [Green Version]
  38. Honda, K.; Taniguchi, T. IRFs: Master regulators of signalling by Toll-like receptors and cytosolic pattern-recognition receptors. Nat. Rev. Immunol. 2006, 6, 644–658. [Google Scholar] [CrossRef]
  39. Grant, C.E.; Vasa, M.Z.; Deeley, R.G. cIRF-3, a new member of the interferon regulatory factor (IRF) family that is rapidly and transiently induced by dsRNA. Nucleic Acids Res. 1995, 23, 2137–2146. [Google Scholar] [CrossRef] [Green Version]
  40. Cormican, P.; Lloyd, A.T.; Downing, T.; Connell, S.J.; Bradley, D.; O’Farrelly, C. The avian Toll-Like receptor pathway—Subtle differences amidst general conformity. Dev. Comp. Immunol. 2009, 33, 967–973. [Google Scholar] [CrossRef]
  41. Crotta, S.; Davidson, S.; Mahlakoiv, T.; Desmet, C.J.; Buckwalter, M.R.; Albert, M.L.; Staeheli, P.; Wack, A. Type I and type III interferons drive redundant amplification loops to induce a transcriptional signature in influenza-infected airway epithelia. PLoS Pathog. 2013, 9, e1003773. [Google Scholar] [CrossRef] [Green Version]
  42. Giotis, E.S.; Ross, C.S.; Robey, R.C.; Nohturfft, A.; Goodbourn, S.; Skinner, M.A. Constitutively elevated levels of SOCS1 suppress innate responses in DF-1 immortalised chicken fibroblast cells. Sci. Rep. 2017, 7, 17485. [Google Scholar] [CrossRef] [Green Version]
  43. Hatesuer, B.; Hoang, H.T.; Riese, P.; Trittel, S.; Gerhauser, I.; Elbahesh, H.; Geffers, R.; Wilk, E.; Schughart, K. Deletion of Irf3 and Irf7 Genes in Mice Results in Altered Interferon Pathway Activation and Granulocyte-Dominated Inflammatory Responses to Influenza A Infection. J. Innate Immun. 2017, 9, 145–161. [Google Scholar] [CrossRef]
  44. Pearson, G.; Robinson, F.; Beers Gibson, T.; Xu, B.E.; Karandikar, M.; Berman, K.; Cobb, M.H. Mitogen-activated protein (MAP) kinase pathways: Regulation and physiological functions. Endocr. Rev. 2001, 22, 153–183. [Google Scholar] [CrossRef] [Green Version]
  45. Borgeling, Y.; Schmolke, M.; Viemann, D.; Nordhoff, C.; Roth, J.; Ludwig, S. Inhibition of p38 mitogen-activated protein kinase impairs influenza virus-induced primary and secondary host gene responses and protects mice from lethal H5N1 infection. J. Biol. Chem. 2014, 289, 13–27. [Google Scholar] [CrossRef] [Green Version]
  46. Harjunpaa, H.; Llort Asens, M.; Guenther, C.; Fagerholm, S.C. Cell Adhesion Molecules and Their Roles and Regulation in the Immune and Tumor Microenvironment. Front. Immunol. 2019, 10, 1078. [Google Scholar] [CrossRef] [Green Version]
  47. Saxton, R.A.; Sabatini, D.M. mTOR Signaling in Growth, Metabolism, and Disease. Cell 2017, 169, 361–371. [Google Scholar] [CrossRef]
  48. Laplante, M.; Sabatini, D.M. mTOR signaling at a glance. J. Cell Sci. 2009, 122, 3589–3594. [Google Scholar] [CrossRef] [Green Version]
  49. Colina, R.; Costa-Mattioli, M.; Dowling, R.J.; Jaramillo, M.; Tai, L.H.; Breitbach, C.J.; Martineau, Y.; Larsson, O.; Rong, L.; Svitkin, Y.V.; et al. Translational control of the innate immune response through IRF-7. Nature 2008, 452, 323–328. [Google Scholar] [CrossRef]
  50. Weichhart, T.; Hengstschlager, M.; Linke, M. Regulation of innate immune cell function by mTOR. Nat. Rev. Immunol. 2015, 15, 599–614. [Google Scholar] [CrossRef]
  51. Jones, R.G.; Pearce, E.J. MenTORing Immunity: mTOR Signaling in the Development and Function of Tissue-Resident Immune Cells. Immunity 2017, 46, 730–742. [Google Scholar] [CrossRef] [Green Version]
  52. McNab, F.; Mayer-Barber, K.; Sher, A.; Wack, A.; O’Garra, A. Type I interferons in infectious disease. Nat. Rev. Immunol. 2015, 15, 87–103. [Google Scholar] [CrossRef]
  53. Carlin, A.F.; Plummer, E.M.; Vizcarra, E.A.; Sheets, N.; Joo, Y.; Tang, W.; Day, J.; Greenbaum, J.; Glass, C.K.; Diamond, M.S.; et al. An IRF-3-, IRF-5-, and IRF-7-Independent Pathway of Dengue Viral Resistance Utilizes IRF-1 to Stimulate Type I and II Interferon Responses. Cell Rep. 2017, 21, 1600–1612. [Google Scholar] [CrossRef] [Green Version]
  54. Karpala, A.J.; Lowenthal, J.W.; Bean, A.G. Activation of the TLR3 pathway regulates IFNβ production in chickens. Dev. Comp. Immunol. 2008, 32, 435–444. [Google Scholar] [CrossRef]
  55. Karpala, A.J.; Stewart, C.; McKay, J.; Lowenthal, J.W.; Bean, A.G. Characterization of chicken Mda5 activity: Regulation of IFN-β in the absence of RIG-I functionality. J. Immunol. 2011, 186, 5397–5405. [Google Scholar] [CrossRef]
  56. Cornelissen, J.B.; Post, J.; Peeters, B.; Vervelde, L.; Rebel, J.M. Differential innate responses of chickens and ducks to low-pathogenic avian influenza. Avian Pathol. 2012, 41, 519–529. [Google Scholar] [CrossRef]
  57. Wang, Y.; Yin, Y.; Lan, X.; Ye, F.; Tian, K.; Zhao, X.; Yin, H.; Li, D.; Xu, H.; Liu, Y.; et al. Molecular characterization, expression of chicken TBK1 gene and its effect on IRF3 signaling pathway. PLoS ONE 2017, 12, e0177608. [Google Scholar] [CrossRef]
  58. Hiscott, J.; Lin, R.; Nakhaei, P.; Paz, S. MasterCARD: A priceless link to innate immunity. Trends Mol. Med. 2006, 12, 53–56. [Google Scholar] [CrossRef]
Genes 11 00385 g001 550
Figure 1. The establishment of IRF7-/- and control (IRF7wt) DF-1 cell lines. (A) Left panel: a scheme outlining the IRF7 knockout by the mutagenic chain reaction (MCR). An MCR plasmid targeting the second exon of chicken IRF7 was designed. Three adjacent guide RNAs (gRNAs) targeting the IRF7 exon 2 locus were selected, and each gRNA was cloned into the designed MCR plasmid. The locations of the PCR primers used for the analysis of the genomic insertion site are indicated. Right panel: the SpCas9-puromycin resistant gene construct from the MCR plasmid was cloned into the piggyBac transposon expression vector to generate control (IRF7wt) DF-1 cells. (B) PCR analysis of bi-allelic MCR insertion cell clones for each gRNA is shown. The absence of a PCR band derived from the IRF7 exon2 locus (Exon2) and the presence of junctional bands corresponding to MCR insertion into the chromosomal IRF7 exon 2 locus (HA1, HA2) are shown for MCR cell clones. (C) RNA-seq alignment visualization of IRF7wt (upper panel) and IRF7-/- (lower panel) in non-infected cells (0 h). The genomic region of exons 1 and 2 of IRF7 is shown. The gray vertical dashed line indicates a putative double strand break and homology directed repair site. (D) Culture supernatants were collected at each time point and the viral titer was determined by the plaque assay as (plaque forming unit (PFU)+1)/ml using Madin-Darby kidney (MDCK) cells. (E–G) The qRT-PCR analysis of IRF7, IFNA and IFNB transcripts in IRF7-/- and IRF7wt cells upon H6N2 infection (1 multiplicity of infection (MOI)). All data are shown as mean ± SEM from six biological replicates (* p < 0.05, ** p < 0.01, *** p < 0.001, NS: not significant; two tailed T-test).
Figure 1. The establishment of IRF7-/- and control (IRF7wt) DF-1 cell lines. (A) Left panel: a scheme outlining the IRF7 knockout by the mutagenic chain reaction (MCR). An MCR plasmid targeting the second exon of chicken IRF7 was designed. Three adjacent guide RNAs (gRNAs) targeting the IRF7 exon 2 locus were selected, and each gRNA was cloned into the designed MCR plasmid. The locations of the PCR primers used for the analysis of the genomic insertion site are indicated. Right panel: the SpCas9-puromycin resistant gene construct from the MCR plasmid was cloned into the piggyBac transposon expression vector to generate control (IRF7wt) DF-1 cells. (B) PCR analysis of bi-allelic MCR insertion cell clones for each gRNA is shown. The absence of a PCR band derived from the IRF7 exon2 locus (Exon2) and the presence of junctional bands corresponding to MCR insertion into the chromosomal IRF7 exon 2 locus (HA1, HA2) are shown for MCR cell clones. (C) RNA-seq alignment visualization of IRF7wt (upper panel) and IRF7-/- (lower panel) in non-infected cells (0 h). The genomic region of exons 1 and 2 of IRF7 is shown. The gray vertical dashed line indicates a putative double strand break and homology directed repair site. (D) Culture supernatants were collected at each time point and the viral titer was determined by the plaque assay as (plaque forming unit (PFU)+1)/ml using Madin-Darby kidney (MDCK) cells. (E–G) The qRT-PCR analysis of IRF7, IFNA and IFNB transcripts in IRF7-/- and IRF7wt cells upon H6N2 infection (1 multiplicity of infection (MOI)). All data are shown as mean ± SEM from six biological replicates (* p < 0.05, ** p < 0.01, *** p < 0.001, NS: not significant; two tailed T-test).
Genes 11 00385 g001
Genes 11 00385 g002 550
Figure 2. IRF7 deletion altered transcriptome regulation in DF-1 cells (A) Volcano plots of transcriptomic differences between IRF7-/- and IRF7wt cells at each time point; the colored dots correspond to significant differentially expressed genes (DEGs; False discovery rate (FDR) < 5 %). The numbers of DEGs are shown at the bottom of each plot, and the numbers of genes that have more than 2-fold changes are in the parentheses. (B) Gene ontology (GO) and pathway analysis by DAVID 6.8 using DEGs from each time point. The numbers of genes enriched in each biological process are in parentheses. (C–E) Heatmaps showing the expression fold change of DEGs in influenza A, MAPK (mitogen-activated protein kinase) signaling and the phagosome pathways at all time points (IRF7-/-/IRF7wt). Lists of DEGs from all time points in each pathway were combined to generate each heatmap.
Figure 2. IRF7 deletion altered transcriptome regulation in DF-1 cells (A) Volcano plots of transcriptomic differences between IRF7-/- and IRF7wt cells at each time point; the colored dots correspond to significant differentially expressed genes (DEGs; False discovery rate (FDR) < 5 %). The numbers of DEGs are shown at the bottom of each plot, and the numbers of genes that have more than 2-fold changes are in the parentheses. (B) Gene ontology (GO) and pathway analysis by DAVID 6.8 using DEGs from each time point. The numbers of genes enriched in each biological process are in parentheses. (C–E) Heatmaps showing the expression fold change of DEGs in influenza A, MAPK (mitogen-activated protein kinase) signaling and the phagosome pathways at all time points (IRF7-/-/IRF7wt). Lists of DEGs from all time points in each pathway were combined to generate each heatmap.
Genes 11 00385 g002
Genes 11 00385 g003 550
Figure 3. Comparative analysis of IRF7 knockout and overexpression. The expression profile of 245 common DEGs between IRF7 knockout (IRF7-/-/IRF7wt) and IRF7 overexpression (CuO-IRF7/Control) [35] at 0 h (mock) and 6 hpi with H6N2 infection (heatmap). Each cluster’s representative genes are indicated on the left of the heatmap, and biological process gene ontology terms and their p-values are indicated on the right.
Figure 3. Comparative analysis of IRF7 knockout and overexpression. The expression profile of 245 common DEGs between IRF7 knockout (IRF7-/-/IRF7wt) and IRF7 overexpression (CuO-IRF7/Control) [35] at 0 h (mock) and 6 hpi with H6N2 infection (heatmap). Each cluster’s representative genes are indicated on the left of the heatmap, and biological process gene ontology terms and their p-values are indicated on the right.
Genes 11 00385 g003
Genes 11 00385 g004 550
Figure 4. IRF7 deletion altered the host response against H6N2 infection (A) The number of differentially expressed genes (DEGs; FDR < 1%) for IRF7wt and IRF7-/- cell lines upon H6N2 infection at 2, 6 and 12 h post infections (hpi). Each horizontal comparison indicates the number of DEGs upon H6N2 infection at each time point compared to mock infection (0 h). The numbers of genes that have 2-fold expression changes are in the parentheses. (B) A Venn diagram of DEGs between IRF7wt and IRF7-/- cells upon H6N2 infection at each time point. (C) Gene ontology functional analysis was performed using each cell-specific DEG. The numbers of genes enriched in each biological process are in parentheses.
Figure 4. IRF7 deletion altered the host response against H6N2 infection (A) The number of differentially expressed genes (DEGs; FDR < 1%) for IRF7wt and IRF7-/- cell lines upon H6N2 infection at 2, 6 and 12 h post infections (hpi). Each horizontal comparison indicates the number of DEGs upon H6N2 infection at each time point compared to mock infection (0 h). The numbers of genes that have 2-fold expression changes are in the parentheses. (B) A Venn diagram of DEGs between IRF7wt and IRF7-/- cells upon H6N2 infection at each time point. (C) Gene ontology functional analysis was performed using each cell-specific DEG. The numbers of genes enriched in each biological process are in parentheses.
Genes 11 00385 g004
Genes 11 00385 g005 550
Figure 5. Gene expression profiles shifted as a result of IRF7 deletion in the host response. (A) The k-means clustering of differentially expressed genes (DEGs; n = 5091) and expression cluster shift of genes due to IRF7 knockout. Respective gene expression is shown in gray, and the representative expression pattern is highlighted in different colors. A chord diagram indicates the number of DEGs shifted clusters from IRF7wt to IRF7-/-. (B) The number of genes in each cluster. A complete list of genes in each cluster are in Table S4. (C) The biological process gene ontology terms of IRF7wt DEGs in each cluster. (D) DEGs in the influenza A pathway and their expression clusters between IRF7wt and IRF7-/- are shown on the left. The colors of the clusters correspond to the colors used in (A). A heatmap of gene expression of the pathway is in the right panel. See Figure S3 for each gene in the influenza A pathway.
Figure 5. Gene expression profiles shifted as a result of IRF7 deletion in the host response. (A) The k-means clustering of differentially expressed genes (DEGs; n = 5091) and expression cluster shift of genes due to IRF7 knockout. Respective gene expression is shown in gray, and the representative expression pattern is highlighted in different colors. A chord diagram indicates the number of DEGs shifted clusters from IRF7wt to IRF7-/-. (B) The number of genes in each cluster. A complete list of genes in each cluster are in Table S4. (C) The biological process gene ontology terms of IRF7wt DEGs in each cluster. (D) DEGs in the influenza A pathway and their expression clusters between IRF7wt and IRF7-/- are shown on the left. The colors of the clusters correspond to the colors used in (A). A heatmap of gene expression of the pathway is in the right panel. See Figure S3 for each gene in the influenza A pathway.
Genes 11 00385 g005
Genes 11 00385 g006 550
Figure 6. A proposed model of the molecular mechanism of chicken IRF7 modulation in the host response to avian influenza virus (AIV) infection A schematic model of the early antiviral response in chickens regarding IRF7 is proposed based on the DEGs in this study. AIV infection is sensed by pattern recognition receptors. TLR3 and TLR7 recognize endosomal/phagosomal double-stranded RNA (dsRNA) and single-stranded RNA (ssRNA) respectively. MDA5 and PKR respond to cytosolic dsRNA. This leads to a subsequent signaling cascade via adaptor molecules and protein kinases that activate the key transcription factors AP1, NF-κB and IRF7. The MAPK and PI3K/Akt pathways are potentially involved in the activation process. Activated IRF7 translocates to the nucleus and activates the transcription of IFNβ. AP1 and NF-κB also translocate to the nucleus and induce type I IFNs and inflammatory cytokines. Type I IFNs activate ISGs via the Jak–STAT pathway in an autocrine/paracrine manner through their receptors (IFNAR). IFIT5, interferon induced protein with tetratricopeptide repeats 5; JAK1/2, Janus kinase 1/2; LY6E, lymphocyte antigen 6 family member E; MKK3/6, MAPK kinase 3/6; MYD88, myeloid differentiation primary response 88; PKB, protein kinase B; STAT1/2, signal transducer and activator of transcription 1/2; TRAF6, TNF receptor associated factor 6; TRIF (TICAM1), toll like receptor adaptor molecule 1; TYK2, tyrosine kinase 2. Red and blue arrows indicate the up- and down-regulation of each gene respectively as a result of IRF7 knockout in the dataset. Further investigation is required for the items indicated by question marks.
Figure 6. A proposed model of the molecular mechanism of chicken IRF7 modulation in the host response to avian influenza virus (AIV) infection A schematic model of the early antiviral response in chickens regarding IRF7 is proposed based on the DEGs in this study. AIV infection is sensed by pattern recognition receptors. TLR3 and TLR7 recognize endosomal/phagosomal double-stranded RNA (dsRNA) and single-stranded RNA (ssRNA) respectively. MDA5 and PKR respond to cytosolic dsRNA. This leads to a subsequent signaling cascade via adaptor molecules and protein kinases that activate the key transcription factors AP1, NF-κB and IRF7. The MAPK and PI3K/Akt pathways are potentially involved in the activation process. Activated IRF7 translocates to the nucleus and activates the transcription of IFNβ. AP1 and NF-κB also translocate to the nucleus and induce type I IFNs and inflammatory cytokines. Type I IFNs activate ISGs via the Jak–STAT pathway in an autocrine/paracrine manner through their receptors (IFNAR). IFIT5, interferon induced protein with tetratricopeptide repeats 5; JAK1/2, Janus kinase 1/2; LY6E, lymphocyte antigen 6 family member E; MKK3/6, MAPK kinase 3/6; MYD88, myeloid differentiation primary response 88; PKB, protein kinase B; STAT1/2, signal transducer and activator of transcription 1/2; TRAF6, TNF receptor associated factor 6; TRIF (TICAM1), toll like receptor adaptor molecule 1; TYK2, tyrosine kinase 2. Red and blue arrows indicate the up- and down-regulation of each gene respectively as a result of IRF7 knockout in the dataset. Further investigation is required for the items indicated by question marks.
Genes 11 00385 g006

Share and Cite

MDPI and ACS Style

Kim, T.H.; Kern, C.; Zhou, H. Knockout of IRF7 Highlights its Modulator Function of Host Response Against Avian Influenza Virus and the Involvement of MAPK and TOR Signaling Pathways in Chicken. Genes 2020, 11, 385. https://doi.org/10.3390/genes11040385

AMA Style

Kim TH, Kern C, Zhou H. Knockout of IRF7 Highlights its Modulator Function of Host Response Against Avian Influenza Virus and the Involvement of MAPK and TOR Signaling Pathways in Chicken. Genes. 2020; 11(4):385. https://doi.org/10.3390/genes11040385

Chicago/Turabian Style

Kim, Tae Hyun, Colin Kern, and Huaijun Zhou. 2020. "Knockout of IRF7 Highlights its Modulator Function of Host Response Against Avian Influenza Virus and the Involvement of MAPK and TOR Signaling Pathways in Chicken" Genes 11, no. 4: 385. https://doi.org/10.3390/genes11040385

APA Style

Kim, T. H., Kern, C., & Zhou, H. (2020). Knockout of IRF7 Highlights its Modulator Function of Host Response Against Avian Influenza Virus and the Involvement of MAPK and TOR Signaling Pathways in Chicken. Genes, 11(4), 385. https://doi.org/10.3390/genes11040385

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop