key: cord-315072-b28yikvj
authors: Giotis, Efstathios S.; Robey, Rebecca C.; Skinner, Natalie G.; Tomlinson, Christopher D.; Goodbourn, Stephen; Skinner, Michael A.
title: Chicken interferome: avian interferon-stimulated genes identified by microarray and RNA-seq of primary chick embryo fibroblasts treated with a chicken type I interferon (IFN-α)
date: 2016-08-05
journal: Vet Res
DOI: 10.1186/s13567-016-0363-8
sha: 
doc_id: 315072
cord_uid: b28yikvj

Viruses that infect birds pose major threats—to the global supply of chicken, the major, universally-acceptable meat, and as zoonotic agents (e.g. avian influenza viruses H5N1 and H7N9). Controlling these viruses in birds as well as understanding their emergence into, and transmission amongst, humans will require considerable ingenuity and understanding of how different species defend themselves. The type I interferon-coordinated response constitutes the major antiviral innate defence. Although interferon was discovered in chicken cells, details of the response, particularly the identity of hundreds of stimulated genes, are far better described in mammals. Viruses induce interferon-stimulated genes but they also regulate the expression of many hundreds of cellular metabolic and structural genes to facilitate their replication. This study focusses on the potentially anti-viral genes by identifying those induced just by interferon in primary chick embryo fibroblasts. Three transcriptomic technologies were exploited: RNA-seq, a classical 3′-biased chicken microarray and a high density, “sense target”, whole transcriptome chicken microarray, with each recognising 120–150 regulated genes (curated for duplication and incorrect assignment of some microarray probesets). Overall, the results are considered robust because 128 of the compiled, curated list of 193 regulated genes were detected by two, or more, of the technologies. ELECTRONIC SUPPLEMENTARY MATERIAL: The online version of this article (doi:10.1186/s13567-016-0363-8) contains supplementary material, which is available to authorized users.

The interferon (IFN) response is one of the most important arms of host innate immunity against virus infection [1, 2] . Infected cells are able to recognise foreign nucleic acids and induce the synthesis and secretion of type I IFN (IFN-α and IFN-β) and type III IFN (IFN-λ), which bind to receptors on the surface of neighbouring cells and trigger the transcriptional regulation of genes involved in the antiviral state. Studies in mammals have demonstrated that there are several hundred such IFN-regulated genes (IRGs). Because the vast majority are up-regulated they are overwhelmingly referred to as IFN-stimulated genes (ISGs) so, hereafter, they will be referred to generically as ISGs (or specifically as chicken ISGs, ChISGs), except where the more generic term avoids confusion. Induction of ISGs involves the JAK/STAT signalling pathway: STAT1 is either recruited directly to target promoters for a relatively weak activation or, more commonly, is recruited in a complex called ISGF3 in association with STAT2 and IRF9 [1, 3] .

ISGs are the focus of considerable current attention with regard to: (i) their antiviral activity, (ii) an increasing appreciation of the complexity of their regulation and (iii) their targeting by virus-encoded modulators of IFN-induced responses [1, 3, 4] . These studies require comprehensive catalogues of the ISGs, especially where system-wide approaches are undertaken. Even though many key mammalian ISGs have been known for some time, it is with the relatively recent advent of transcriptomic technologies that the full complement has been catalogued (mainly using microarrays [5] ; see also Schoggins et al. [6] ).

In contrast to the mammalian IFN system our equivalent knowledge of the avian system has lagged behind. Although IFN was discovered in chickens in 1957 [7] the first chicken IFN gene was characterised in 1994 [8] and the key chicken ISG, PKR, was identified in 2004 [9] . The derivation of the chicken genome sequence, first drafted in 2004 [10] , did not greatly advance our understanding of chicken ISGs because of the incomplete nature of the Gallus gallus genome assembly, even at v4 (Galgal4), which might be partly due to the fact that the chicken karyotype has six pairs of macrochromosomes (but 33 pairs of microchromosomes), and the difficulties in annotating immunity genes, which are some of the most divergent between mammals and birds [11] . However, it has become apparent that key genes of the innate immune system, such as the transcription factors IRF9 and one member of the IRF3/IRF7 dyad [12, 13; unpublished] , are absent from avian species, indicative of significant functional differences between them and mammals. Moreover, for reasons that are not understood, the cytosolic pattern recognition receptor, RIG-I, appears to have been lost from chicken as well as other galliformes [13, 14] .

To generate a chicken ISG database we have compared data from three transcriptomic technology platforms: (i) the classical 3′-biased GeneChip Chicken Genome Array (32K; Affymetrix, High Wycombe, UK), (ii) the Chicken Gene 1.0 Sense Target (ST) whole transcriptome Array (Affymetrix) and (iii) Illumina (Little Chesterford, UK) RNA-seq. This three-way comparison allowed a high level of cross-validation of data from each technology, beyond what would normally be achieved by qRT-PCR. It also allows subsequent studies, constrained to use any particular technology, to be more broadly compared. We monitored IRG expression in chicken embryo fibroblast (CEF) induced for 6 h with 1000 units recombinant chicken IFN-α (rChIFN1; hereafter routinely referred to as IFN), a time chosen to reflect predominantly primary signalling targets. The expression data for selected genes were also validated by PCR and qRT-PCR. Overlapping data show generally high degrees of concordance in the identity of the IRGs and their relative levels of regulation by IFN, with disparity mainly where multiple microarray probes exist for single genes. The study was presented in a preliminary form as a poster at the International Cytokine and Interferon Society (ICIS) meeting ("Cytokines 2015"; October 11-14, 2015) in Bamberg, Germany [15] .

Freshly isolated CEF were provided by the former Institute for Animal Health (Compton, UK, now The Pirbright Institute, Pirbright, UK). Cells were seeded in T25 flasks (Greiner Bio One, Kremsmünster, Austria; 5.6 × 10 6 cells/flask) and cultured overnight in 5.5 mL 199 media (Gibco Thermo Fisher Scientific, Paisley, UK) supplemented with 8% heat-inactivated newborn bovine serum (NBCS; Gibco), 10% tryptose phosphate broth (TPB; Sigma-Aldrich, Gillingham, UK), 2% nystatin (Sigma-Aldrich) and 0.1% penicillin streptomycin (Gibco).

Recombinant chicken IFN-α (rChIFN1) was prepared as previously reported [16] and was added in culture media to a final concentration of 1000 units/mL. Confluent cells were treated with IFN or mock-treated and incubated for six hours before harvesting. Cells were stored at −80 °C in RNAlater (Sigma-Aldrich) until RNA extraction. The experiment was repeated in triplicate with three different batches of CEF.

Total RNA was extracted from cells using an RNeasy kit (Qiagen, Crawley, UK) according to the manufacturer's instructions. On-column DNA digestion was performed using RNase-free DNase (Qiagen) to remove contaminating genomic DNA. RNA samples were quantified using a Nanodrop Spectrophotometer (Thermo Fisher Scientific, Paisley, UK) and checked for quality using a 2100 Bioanalyzer (Agilent Technologies, Wokingham, UK). All RNA samples had an RNA integrity number (RIN) ≥9.6.

RNA samples were processed for microarray with the GeneChip ® Chicken Genome Array (Affymetrix) using the GeneChip ® 3′ IVT Express Kit (Affymetrix) or for microarray with the Chicken Gene 1.0 ST Array (Affymetrix) using the Ambion (Paisley, UK) WT Expression Kit for Affymetrix GeneChip ® Whole Transcript (WT) Expression Arrays (Ambion) and the GeneChip WT Terminal Labelling and Controls Kit (Ambion), following the manufacturers' instructions, as described previously [17] .

Total RNA (100 ng) was used as input and quality checks were performed using the 2100 Bioanalyzer at all stages suggested by the manufacturer. RNA samples were processed in two batches of 18 but batch mixing was used at every stage to avoid creating experimental bias. Hybridisation of RNA to chips and scanning of arrays was performed by the Medical Research Council's Clinical Sciences Centre (CSC) Genomics Laboratory (Hammersmith Hospital, London, UK). RNA was hybridised to GeneChip Chicken Genome Array chips (Affymetrix) in a GeneChip Hybridization Oven (Affymetrix), the chips were stained and washed on a GeneChip Fluidics Station 450 (Affymetrix), and the arrays were scanned in a Gene-Chip Scanner 3000 7G with autoloader (Affymetrix).

cDNA was synthesised from RNA samples from untreated and IFN-treated CEF using the QuantiTect ® Reverse Transcriptase system (Qiagen) according to the manufacturer's instructions. The Glyceraldehyde 3-phosphate dehydrogenase (GAPDH) was used as a reference gene. All target gene expression levels were calculated relative to GAPDH expression levels and the target gene expression level in −2 h uninfected CEF using the comparative C T method (also referred to as the 2 −ΔΔCT method).

Triplicate untreated (control) and IFN-treated CEF were processed for transcriptome analysis by RNA-seq. The cell samples used were identical to those used for the microarray analyses. Total RNA was extracted as for microarrays (above) and RNA libraries were prepared for deep sequencing using the TruSeq RNA Sample Preparation Kit (Illumina) according to the manufacturer's instructions. Total RNA (2.5 μg) was used as an input for each library. A total of six RNA adapter indices were randomly assigned to the 12 samples to allow multiplexing of libraries. At the end of the protocol, libraries were quantified using a Nanodrop Spectrophotometer and checked for quality using a 2100 Bioanalyzer High Sensitivity DNA chip (Agilent Technologies). RNA library qPCR quantification, multiplexing and sequencing was performed by the Medical Research Council's Clinical Sciences Centre (CSC) Genomics Laboratory, Hammersmith Hospital, London, UK. Libraries were quantified using the KAPA Biosystems (London, UK) library quantification kit (KK4824) on an ABI 7500 FAST qPCR machine (Applied Biosystems). Libraries were then diluted to a 2 nM stock solution, pooled for multiplexing, denatured and diluted to a final molarity of 20 pM. Libraries were loaded on to the flow cell (8-16 pM per lane) for clustering and cluster generation was performed by the Illumina cBot using version 3 kits. Sequencing of the flow cell was then carried out on the Illumina HiSeq 2000 using the version 3 kits. Data were processed using 

Microarray data were processed using workflows in GENESPRING ™ (Agilent) and PARTEK ™ (Partek Inc., St Louis, MO, USA) commercial software suites. Data (.CEL files) were analysed and statistically filtered using either Partek Genomic Suite 6.6 (Partek GS) or Genespring version 7.2 (Agilent Technologies) software. Input files were normalized with either GCRMA or Genespring algorithms for gene array on core metaprobesets. A one-way ANOVA was performed using either software across all samples. Statistically significant genes were identified using mixed model analysis of variance with a false discovery rate (Benjamini-Hochberg test) of P < 0.05. Fold-change values <±3.0 were removed.

RNA-seq data were imported into CLC bio's Genomics Workbench (CLC Bio, Aarhus, Denmark; now Qiagen), quality-controlled and thereafter processed using that package (versions 6 and 7).

After quality control, the reads were subjected to quality trimming then mapped against ENSEMBL Galgal4 annotated genes (release 75 [18] ) for quantitative analysis of expression. Fold change and False Discovery Rates (FDR) were calculated using Kal's Z test [19] , with pooled data, or Baggerly's test [20] , using separate triplicates.

Initially, we used the 32K GeneChip ® Chicken Genome Array (Affymetrix) because, as well as displaying probes for 32 773 chicken transcripts, it displays probes for 684 transcripts from 17 different viral pathogens of chickens, which offers advantages to those studying virus infections in a chicken background. Subsequently, we used the more refined Chicken Gene 1.0 ST Array (Affymetrix) because it offers a higher probe density against 18 214 chicken genes and should allow detection of transcript isoforms, including non-polyadenylated and alternatively polyadenylated, though it does not include probes for viral genes.

Separate weekly batches of CEF, produced from pools of eggs from the same flock (Rhode Island Red) held in SPF-like conditions at the former Compton Laboratory of the Institute for Animal Health (now The Pirbright Laboratory) served as biological replicates. Principal component analysis of the microarray data (data not shown) indicated limited variation between batches so, thereafter, biological triplicates were used routinely.

IRGs were identified from expression analysis data determined using the 32K GeneChip following IFN treatment (1000 units, 6 h) of CEF. After quantile normalization, significant hits were identified with GENESPRING using an unpaired T test with asymptotic p-value computation and Benjamini-Hochberg multiple testing correction to generate false discovery rates (FDR). A matrix of FDR (from <0.001 to 1) plotted against fold change (FC; from 1.0 to >3) is shown in Table 2 . A relatively conservative FDR of <0.01 returned 250 differentially expressed probesets. Overlaying this with a value for FC for which changes in expression might reasonably be expected to be readily and reliably assayed using other technologies, namely >3, reduced the number of selected, significant probesets to a manageable 181 (180 up, 1 down). These settings were therefore chosen for further analysis. For 23 of these probe sets, no currently recognised genes were automatically assigned. Of the remaining 158 probe sets, 29 were assigned to genes recognised in duplicate by other probe sets. Consequently 129 recognised genes were identified as differentially expressed (the down-regulated transcript was not, at that time, assigned to a recognised gene).

With the Chicken Gene 1.0 ST Array, 157 probe sets demonstrated differential expression (156 up, 1 down) at the same settings (FC > 3, FDR < 0.01). Amongst these, there were five duplicated probe sets and 27 that were not automatically assigned to recognised genes therefore 125 recognised genes were uniquely identified as differentially regulated.

Illumina RNA-seq yielded a total of 170 million reads (100 bases; paired) for the mock-treated CEF triplicate samples and 167 million for the IFN-treated samples. Upon quality trimming and mapping to ENSEMBL Galgal4 annotated genes (release 75), using CLC Bio's Genomic Workbench, 138 recognised genes were identified as differentially regulated (137 up, 1 down) using Kal's proportion-based Z test [19; as implemented in the CLC Bio package] at the same settings (FC > 3, FDR < 0.01). Kal's is performed on the pooled reads from IFN-treated and untreated samples. It is perhaps, therefore, more widely applicable; it also returned a number of IRGs comparable to those returned by the microarrays. Triplicate-based analysis using Baggerly's proportionbased Beta-binomial test [20; as implemented in the CLC Bio package] at the same settings (FC > 3, FDR < 0.01) returned an additional 37 up-regulated genes.

Comparison of the complete raw gene lists from the three technologies using the most compatible identifier (essentially the Gene Symbol) with an online Venn Diagram tool (Venn Diagram Generator; [21] ) demonstrated that 233 recognised genes were identified as differentially regulated. Of these, 51 were identified in common by all three technologies and a further 57 were identified by two out of three technologies, meaning that 108 were identified by at least two technologies. A total of 125 were therefore each identified only by individual technologies ( Figure 1A) .

As well as comparing the identities of the differentially regulated genes, the correlation of expression of the genes identified by the different platforms was examined in terms of both level and rank of FC (Figures 2A and  B) . For instance, comparing RNA-seq data with the 32K GeneChip data, Spearman correlation values were 0.93 for FC level and rank. Considering the current state of assembly and annotation of the chicken genome, the correlation of ISGs in terms of gene identity as well as the level and rank of induction as indicated by all three technology platforms is reassuring. Nevertheless the platform transcriptomic data were validated for selected genes by RT-PCR (data not shown) and by qRT-PCR ( Figure 3A) . A 6 h time point was chosen for microarray and RNAseq analysis of IFN treatment as it has been widely used and is known to result in significant levels of a broad range of ISGs in mammals, making it suitable for defining the chicken interferome. Use of this single time point does not, however, provide unequivocal insight into mechanistic interpretation of ISG induction; for instance, it does not discriminate between strictly ISRE-dependent induction of ISGs and ISRE-independent induction of ISGs by mechanisms that might include immediate high-level induction of IRF1, which has been observed in mammalian systems [22] [23] [24] . Kinetic analysis of the induction of expression of a subset of ISGs was therefore conducted at 45, 90, 180 and 360 min post application of IFN (see Figure 3B ). Even among highly-induced ISGs, different temporal profiles were observed, from the rapid accumulation of IFIT5 (1000-fold by 90 min) and RSAD2 (which remain at steady levels to 360 min) to the steadier, sustained accumulation of Mx and the more modestly induced STAT1; with LGP2 and TRIM25 peaking at 180 min. Although differences in mRNA stability and turnover will influence the profiles, this identification of the ISGs will allow detailed analysis of their promoters to investigate elements (and the factors that bind them) that contribute to the complexity of the observed induction patterns.

Of the 51 IRGs initially identified by all three technologies, 47 ). This suggests either that the mammalian equivalents are ISGs but that they are not included as such in Interferome or that they are not ISGs in mammals. The raw lists were refined by manual "curation", allowing for synonyms of recognised genes (for instance ISG12-2 versus ISG12(2)) and, after bioinformatic analysis using BLAST, etc., assigning recognised gene identifiers to probe sets that previously lacked them. At the end of this process ( Figure 1B ; Additional files 1, 2), it was apparent that some (n = 12) differentially regulated genes identified by the microarrays were also identified as differentially regulated by RNA-seq but that they fell outside of the strict FC > 3 and FDR < 0.01 parameters, reflecting unsurprising disparity in the sensitivity of the three technologies. Those genes that were expressed down to FC > 2.5 or with an FDR up to < 0.05 were, therefore, also incorporated to produce a final list ( Figure 1C ; Additional files 1, 2).

It is obvious that this manual curation of the data, to allow for alternative Gene ID nomenclature used by the three technologies and for differences in sensitivity, introduced minor changes to the figures from the automatic comparisons cited above (Figure 1 ; Additional files 1, 2). Curation, therefore, reduced the number of IRGs from 233 to 193. It also increased the number of differentially expressed genes detected by two out of three technologies from 108 to 118 (compare Figures 1A and B) . Relaxing the criteria for detection of differentially regulated genes by RNA-seq (to FC > 2.5 and/or FDR < 0.05) further increased the number of genes detected by all three technologies from 70 to 72 (representing 37%) or by at least two of the technologies from 118 to 128 (66%), leaving 65 genes detected by single technologies (compare Figures 1B and C) , with 29 of those detected by RNA-seq alone (using the Kal's test, at FC > 3.0 and FDR < 0.01; Additional files 1, 2).

Of the 37 additional ISGs identified by RNA-seq as significant (FC > 3 FDR < 0.01) by the more sensitive Baggerly's test but not by Kal's (Table 3) , two were also identified as significant by Kal's using the relaxed criteria (FDR < 0.05). Baggerly's, therefore, identified 35 ISGs additional to those described in the above analyses using RNA-seq (Kal's analysis) and the microarrays (Table 3) 

Analysis of RNA-seq data depends directly on the extant annotated genome sequence. Perhaps not surprisingly therefore, RNA-seq identified the largest proportion of genes amongst the set of 193 unique IRGs that we compiled (150; 78%). Nevertheless, the microarrays each identified 63% of the genes (122 and 121) . Congruence was highest, and almost identical, between RNA-seq and each microarray (98 and 95; 51 ± 1%; all percentages referring to the total of 193 unique IRGs). Between microarrays it fell to 41% (79). For two-way-only comparisons, the distribution of unique genes between the microarrays was symmetrical (42 and 43; 22%). Between RNA-seq and each microarray, unique genes were biased >2-fold towards RNA-seq: 52 (27%) versus 24 (12%) against the Genechip and 55 (28%) versus 26 (13%) against the ST Array.

Clearly in simple terms of numbers of IRGs identified, RNA-seq outperforms the microarrays. This is probably attributable to the historic nature of the array design based on earlier genome assemblies and annotations, with consequent effects on overall coverage (which might disproportionately affect conditionally expressed genes such as those of the innate immune responses). Nevertheless, the ability of microarrays to quantify expression of 50% (about 100) of such a large pool of important genes will often prove sufficient for the experimental objectives where other considerations might affect the choice of technology (see below).

Moving away from actual numbers of genes, it is worth noting that deeper analysis (in the form of validation by alternative approaches) will, by definition, be required to determine which of the genes identified uniquely as IRGs by individual technologies are actually IRGs.

Genomic loci for each of the predicted ISGs were visually inspected using Genomic Workbench's genome browser, displaying tracks showing: gene, transcript, exon and ORF annotations for the current chicken genome build as well as read-mapping for control and IFN-treated reads [27] . On occasions, such inspection revealed the presence of non-annotated, inducibly-transcribed regions, representing exons, whole genes or even gene families. Examples include those previously described at the chicken IFITM locus [28; data not shown], at the HERC locus (described below) or downstream of CCL19 (LOC100857191; "C-C motif chemokine 26-like"; Figure 4 ). Systematic analysis of these ISGs is outside the scope of this manuscript but the data deposited from this study (European Nucleotide Archive (ENA) study number PRJEB7620 [44] ) will facilitate ongoing study and improved annotation. In some cases, although not currently annotated on the ENSEMBL chicken genome, the genes have IDs in NCBI and were identified as ISGs by one of the microarrays. Examples of these include LOC415756, LOC415922 ("guanylate-binding protein 4-like") and LOC422513 ("hect domain and RLD 4-like", a member of the HERC family, discussed below). 

About 10% of the reads from CEFS did not map to the current chicken genome. The unmapped reads combined from the control and IFN-treated samples were assembled into contigs using the de novo assembly function of Genomic Workbench. The RNA-seq function of Genomic Workbench was then used to quantitate expression of the contigs in control and IFN-treated samples. One of the most highly-expressed contigs was one which, when analysed by BLAST, proved to represent a homologue of STAT2, which is missing from the current ENSEMBL annotated reference chicken genome In (B) PYURF shows 24-fold suppression by IFN but the sequence surrounding PYURF shows 87-fold induction from the right-hand end of the unannotated, antisense LOC422513 and considerably higher upregulation from the left-hand end (due to its lower uninduced levels), consistent with these representing homologues of IFN-inducible human genes HERC6 and HERC5. assembly (Galgal4; release 84), though at NCBI it has recently been placed as a Refseq gene on chromosome 33 in the new assembly Galgal5 (an annotated form of which has not yet been released and is currently not scheduled for release). The de novo assembled contig sequence was used to derive primers for RT-PCR; characterisation of chicken STAT2 will be reported elsewhere.

The data on differential expression showed an overwhelming over-representation of genes up-regulated by IFN. For each of the technologies, only one gene was detected as down-regulated. Corresponding GeneIDs were PYURF (PIGY upstream reading frame; ENSGALG00000026229) for RNA-seq and PIGY (phosphatidylinositol glycan anchor biosynthesis, class Y; NCBI GeneID: 101748971) for the ST array. The down-regulated 32K GeneChip probe (Gga.8802.1.S1_at), though not mapped to a known gene at the time of initial processing, according to the Affymetrix NetAffx ™ Analysis Center [29] is now also assigned as PYURF. In humans, PIGY and PYURF represent different open reading frames on the same spliced transcript of a gene on Hs chromosome 4 located downstream of HERC6 then HERC5. The PYURF/PIGY gene is overlapped on the opposite strand by HERC3, which extends downstream to be followed by FAM13A. Similarly, the chicken PIGY (NCBI) and PYURF (Ensembl) genes map to a locus lying upstream of HERC3 then FAM13A on Gg chromosome 4 (see Figure 4) , with HERC-like LOC422513 ("hect domain and RLD 4-like") starting upstream but spanning and extending downstream of the chicken PYURF. Our RNA-seq data ( Figure 4) indicate that this locus is poorly annotated and demonstrates complex regulation of the component genes by IFN. Thus, although the PIGY/ PYURF transcript is down-regulated by IFN, as recorded by all three technologies, it appears to be closely flanked upstream and downstream by still unannotated multiple exons that are clearly strongly induced by IFN (Figure 4 ). Sequences within these upstream and downstream regions (which are represented by the single NCBI Refseq (Gal-gal5) gene, LOC422513, but appear as though they may represent two separate genes, Figure 4 ) bear homology with genes of the HERC family, consistent with the fact that HERC5 neighbours the human PIGY/PYURF gene and that HERC3 neighbours the chicken PIGY/PYURF gene. The chicken HERC3 gene shows no evidence of induction by IFN.

Description of the interferon-inducibility of the ChISGs serves as the first step in understanding the regulation of their expression and their role in anti-viral (and potentially broader anti-microbial) activities. There is considerable current interest in the antiviral responses of particular cell types, particularly those of the lymphoid, myeloid and dendritic lineages. However, the definition of a wide variety of these cell types is not so advanced in avian species so we felt it best to produce baseline data for readily available, primary cells, namely chick embryo fibroblasts (CEF) as they are highly responsive to IFN. They also remain important for commercial production of vaccine viruses (including human vaccines) as well as for the routine isolation and diagnosis of avian pathogens.

Given the currently incomplete nature of the chicken genome assembly (even at Galgal5) and of its annotation (as currently available for Galgal4 and even as awaited for Galgal5) it is inevitable that updates will continue to be released but the primary data reported here, and publicly-available, for microarrays and RNA-seq, can always be applied to updated microarray assignments as well as to subsequent genome assemblies and annotations.

All things being equal, RNA-seq would seem to be the method of choice for transcriptomic analysis of chicken IFN responses, particularly given its ability to produce high-resolution quantitative and qualitative data. Moreover the data are readily portable and can be easily mined by others with different research focus. They can also be applied immediately to newly released genome assemblies and annotations (whether global or local), whereas microarray analysis must await the generation of annotation updates for each technology.

However, although the cost of sequencing has fallen, and will probably continue to do so, there remain considerable overheads to handling large data sets from extensive, complicated experiments, especially in terms of computing and data storage capacity, as well as speed of processing and archiving. For such experiments, microarrays continue to offer a tractable approach, capable of quickly quantifying and comparing the expression of the central core of IRGs producing relatively compact data for rapid analysis and easy archiving.

Induction of innate responses with PAMPS will trigger different or broader ranges of responses by virtue of the fact that they will trigger other or more pathways than just the IFN-pathway. For instance we (Giotis et al. unpublished) and others [12] have begun to analyse the responses induced by the dsRNA analogue poly[I:C]. Regulation of ISG expression might affect the innate responses observed in different cell lines or tissues so it will be important to understand the mechanisms involved. Additionally, we have observed suppression of ISG induction in the spontaneously immortalized chicken fibroblast cell line, DF-1 [30] , due to their enhanced basal expression of the regulatory ISG, SOCS1 (Giotis et al., unpublished) . Identification of the ISGs means that their promoters, enhancers and other regulatory elements can be systematically analysed to help understand the complex kinetics of expression of their expression (Figure 4 ).

Several studies have investigated changes in host gene expression in response to infection in vivo or in culture with particular avian viruses [31] [32] [33] [34] [35] [36] [37] [38] [39] . Although many of these genes will represent innate (and potentially antiviral) host responses, the majority will be involved in the metabolic, cell cycle and ultrastructural changes that the virus has to induce to facilitate replication. Furthermore, it is not unusual for viruses to modulate the expression of signalling molecules key to the antiviral responses or of antiviral effectors themselves. For instance, we have shown that even an attenuated strain of fowlpox virus blocks induction of IFN-β (ChIFN2) and is highly resistant to the antiviral activity induced by IFN [16, 40] .

The results of existing and future studies of infection in vivo or in culture with particular avian viruses can now be compared with data presented here for ISG induction by IFN to look for evidence of modulation of ISG expression by viruses, whether that be modulation of individual ISGs, subsets [4] or the complete set. For instance, fowlpox virus blocks essentially all ISG expression but a mutant defective in the fpv012 ankyrin repeat/F-box protein identified by Laidlaw et al. [40] induces modest levels of a subset of the ISGs (Giotis et al., unpublished) . Such analyses can be extended to important avian zoonotic viruses and pathogens with huge impact on the global poultry industry. Although this study relates to type I IFN, extensive comparison with the effects of type III IFN could now be conducted, extending on the qRT-PCR comparison made by Masuda et al., who looked at induction of Mx and OAS by IFN-β, IFN-γ and IFN-λ [41] .

Interferons and viruses: an interplay between induction, signalling, antiviral responses and virus countermeasures

Inborn errors of anti-viral interferon immunity in humans

Interferon-stimulated genes: a complex web of host defenses

Pathogenic influenza viruses and coronaviruses utilize similar and contrasting approaches to control interferon-stimulated gene responses

Functional classification of interferon-stimulated genes identified using microarrays

A diverse range of gene products are effectors of the type I interferon antiviral response

Virus interference. I. The interferon

Chicken interferon gene: cloning, expression, and analysis

Characterization of the chicken PKR: polymorphism of the gene and antiviral activity against vesicular stomatitis virus

International Chicken Genome Sequencing C (2004) Sequence and comparative analysis of the chicken genome provide unique perspectives on vertebrate evolution

Evidence of the adaptive evolution of immune genes in chicken

Functional analysis of chicken IRF7 in response to dsRNA analog Poly(I:C) by integrating overexpression and knockdown

Defense genes missing from the flight division

Innate sensing of viruses by pattern recognition receptors in birds

ID: 217: transcriptomic analysis of the chicken interferome

Genetic screen of a library of chimeric poxviruses identifies an ankyrin repeat protein involved in resistance to the avian type I interferon response

Species difference in ANP32A underlies influenza A virus polymerase host restriction

Dynamics of gene expression revealed by comparison of serial analysis of gene expression transcript profiles from yeast grown on two different carbon sources

Differential expression in SAGE: accounting for normal between-library variation

Involvement of the IRF-1 transcription factor in antiviral responses to interferons

Constitutive expression of an ISGF2/IRF1 transgene leads to interferon-independent activation of interferon-inducible genes and resistance to virus infection

IFN regulatory factor-1 bypasses IFN-mediated antiviral effects through viperin gene induction

Interferome v2.0: an updated database of annotated interferon-regulated genes

Chicken interferon-inducible transmembrane protein 3 restricts influenza viruses and lyssaviruses in vitro

The DF-1 chicken fibroblast cell line: transformation induced by diverse oncogenes and cell death resulting from infection by avian leukosis viruses

Transcriptomic profiling of virus-host cell interactions following chicken anaemia virus (CAV) infection in an in vivo model

Molecular responses to the influenza A virus in chicken trachea-derived cells

Early host responses to avian influenza A virus are prolonged and enhanced at transcriptional level depending on maturation of the immune system

Transcriptional analysis of host responses to Marek's disease viral infection

Analysis of the early immune response to infection by infectious bursal disease virus in chickens differing in their resistance to the disease

A comparative analysis of host responses to avian influenza infection in ducks and chickens highlights a role for the interferon-induced transmembrane proteins in viral resistance

Analysis of the crow lung transcriptome in response to infection with highly pathogenic H5N1 avian influenza virus

Integrated analysis of microRNA expression and mRNA transcriptome in lungs of avian influenza virus infected broilers

Differential expression of microRNAs in Marek's disease virus-transformed T-lymphoma cell lines

Genetic screen of a mutant poxvirus library identifies an ankyrin repeat protein involved in blocking induction of avian type I interferon

Biological effects of chicken type III interferon on expression of interferon-stimulated genes in chickens: comparison with type I and type II interferons

We are grateful for the skilled support of Laurence Game, Nathalie Lambie and Adam Giess of the Medical Research Council's (MRC) Clinical Sciences Centre's (CSC) Genomics Facility in conducting microarray analysis and Illumina sequencing. We gratefully acknowledge Sarah Butcher and Geraint Barton of the Bioinformatics Support Service at Imperial College London for their advice.

The datasets supporting the conclusions of this article are available from the following repositories: European Bioinformatics Institute (EBI) ArrayExpress accession numbers E-MTAB-3711 (for the 32K GeneChip; [42] ) and E-MTAB-3712 (for the ST array; [43] ). European Nucleotide Archive (ENA) study number PRJEB7620 (for Illumina RNA-seq; [44] ).

Additional file 1. Table of Additional file 2. Detailed information on ChISGs identified by RNA-seq, and microarray technologies (1). Technologies identifying significant IRGs are listed as "1" RNA-seq (using Kal's Z test); "2" Affymetrix 32K GeneChip Chicken Genome Array and "3" Chicken Gene 1.0 ST Array' . ChISGs significant by one or both microarrays and RNA-seq using Kal's Z test under relaxed criteria (FC > 2.5 or FDR < 0.05) are indicated by "(1)". "+" after the technology identifier indicates that IFN-induced RNA-seq read density was observed at the location of the unannotated gene. (2) Interferome status [45] . (3) Human homologue data (HUGO) [46] . (4) Mouse orthologue data (MGI) [47] .

IFN: interferon; IRGs: IFN-regulated genes; ISGs: IFN-stimulated genes; CEF: chicken embryo fibroblasts; rChIFN1: recombinant chicken IFN-α; RIN: RNA integrity number; qRT-PCR: quantitative real-time PCR; GAPDH: glyceraldehyde 3-phosphate dehydrogenase; FC: fold change; FDR: false discovery rate.

The authors declare that they have no competing interests.

ESG and RCC design of the study, data acquisition and analysis, drafting the manuscript. NGS data compilation and analysis, drafting the manuscript. CDT design, production, curation and maintenance of ChISG Browser website. SG design of the study, critically reviewing the manuscript. MAS design of the study, data analysis, finalizing manuscript. All authors read and approved the final manuscript.