- Research article
- Open Access
Transcriptomic analysis of immune response to bacterial lipopolysaccharide in zebra finch (Taeniopygia guttata)
BMC Genomicsvolume 20, Article number: 647 (2019)
Despite the convergence of rapid technological advances in genomics and the maturing field of ecoimmunology, our understanding of the genes that regulate immunity in wild populations is still nascent. Previous work to assess immune function has relied upon relatively crude measures of immunocompetence. However, with next-generation RNA-sequencing, it is now possible to create a profile of gene expression in response to an immune challenge. In this study, captive zebra finch (Taeniopygia guttata; adult males) were challenged with bacterial lipopolysaccharide (LPS) or vehicle to stimulate the innate immune system. 2 hours after injection, birds were euthanized and hypothalami, spleen, and red blood cells (RBCs) were collected. Taking advantage of the fully sequenced genome of zebra finch, total RNA was isolated, sequenced, and partially annotated in these tissue/cells.
In hypothalamus, there were 707 significantly upregulated transcripts, as well as 564 and 144 in the spleen and RBCs, respectively, relative to controls. Also, 155 transcripts in the hypothalamus, 606 in the spleen, and 61 in the RBCs were significantly downregulated. More specifically, a number of immunity-related transcripts (e.g., IL-1β, RSAD2, SOCS3) were upregulated among tissues/cells. Additionally, transcripts involved in metabolic processes (APOD, LRAT, RBP4) were downregulated.
These results suggest a potential trade-off in expression of genes that regulate immunity and metabolism in birds challenged with LPS. This finding is consistent with a hypothermic response to LPS treatment in small birds. Unlike mammals, birds have nucleated RBCs, and these results support a novel transcriptomic response of avian RBCs to immune challenge.
Resources are rarely ubiquitous within an environment. This is especially true for organisms that migrate and/or experience seasonal alterations in life-histories. The energetically costly activities of reproduction, molt, and migration are typically cyclic . In most birds, reproduction occurs on a seasonal basis when food is most abundant and generally does not conflict with other life-history events, such as molting and migration [2, 3]. Resources are also needed to combat pathogens, and immune defense is allocated strategically among life-history activities (e.g., sexual signaling, reproduction, growth; [4, 5]). These trade-offs of immune defense with other behavioral and physiological activities can serve as proximate underpinnings that shape life-history decisions .
The innate immune system involves rapid, indiscriminate responses , initiated by pattern recognition receptors that are highly conserved among vertebrates [8, 9]. The acute phase response (APR) is the first level of defense against infection as part of the innate immune system . This rapid response, within hours of infection, is characterized by behavioral and physiological alterations. Stereotypical behavioral changes are known as “sickness behaviors” and include lethargy, somnolence, reduced food and water intake, and decreased activities such as grooming [10, 11]. Physiological changes occur through upregulation of genes that lead to translation of proteins that regulate fever and inflammation (interleukin-1β (IL-1β), interleukin-6 (IL-6), tumor necrosis factor-alpha (TNF-α)), and downregulation of genes associated with reproduction and growth ( [12, 13]. The APR can also be triggered by exposure to bacterial lipopolysaccharide (LPS) derived from gram-negative bacteria in a dose-dependent manner . Importantly, use of LPS avoids the confounding effect of pathogen manipulation upon the immune system, such that measurement of host immune responses can be accurately quantified .
Most studies of the immune system have focused on adaptive immunity which takes longer to develop than innate immunity [15, 16]. Of the studies that have examined the innate immune system, and more specifically the APR, most have used mammals and domesticated birds [8, 17, 18]. Studies on free-living organisms are less common. The interdisciplinary field of ecoimmunology aims to understand how immune responses are related to host fitness, along with environmental and genetic variability, in non-model organisms . However, ecoimmunology studies have typically used few markers to measure immunocompetence . With newer technology, analysis of gene expression in these non-model organisms can give further understanding into the APR and trade-offs that occur with other life-history functions.
High-throughput RNA-sequencing (RNA-seq) presents advantages in non-model systems in conjunction with those model systems already in use. This method allows for the mapping and quantifying of whole transcriptomes and does not depend on an existing genomic sequence . RNA-seq technology can compare expression levels in different tissues as well as facilitate identification of genes that regulate the response to infection. The relatively low cost and increased sensitivity of RNA-seq compared to other sequencing methods  also makes it a viable option for non-model system studies.
Previous studies using RNA-seq on zebra finch (Taeniopygia guttata) documented many immune-related genes, such as the major histocompatibility complex, to be constitutively expressed in a tissue-specific manner . The zebra finch is one of the first birds with a fully sequenced genome . One study has assessed the effect of West Nile virus (WNV) infection, in which transcriptomes were analyzed 2 days post-inoculation . However, no study to date has measured rapid transcriptomic responses of birds to an immune challenge that is isolated from pathogen manipulation. The previous work in zebra finch makes it an ideal candidate organism to pinpoint functionally important immune system genes.
This study examines the transcriptomic response of hypothalamus, spleen, and red blood cells (RBCs) in zebra finch following acute challenge with LPS. We predicted an upregulation in expression of immune-regulated genes and a corresponding down-regulation of genes associated with growth and/or reproduction, which would suggest a molecular mechanism that could mediate a trade-off between immune defense and other life-history activities.
Transcriptome and gene ontology (GO) analysis
Sequencing generated 17.5 to 63.6 million reads per sample with a mean of 32 million and standard deviation of 14.8 million. The alignment rate ranged from 69.8 to 84.4% with a mean of 80.9 across the 24 samples (Table 1). In hypothalamus, 707 genes were significantly upregulated (of which 628 had associated gene symbols) among LPS-treated birds compared with control samples. In spleen and RBCs, 564 genes (439 with associated gene symbols) and 144 genes (121 with associated gene symbols) were significantly upregulated, respectively, in LPS birds relative to controls (Additional files 1, 2 and 3; Fig. 1a, b). Additionally, 155 genes (134 with associated gene symbols) in the hypothalamus, 606 genes (517 with associated gene symbols) in the spleen, and 61 genes (all associated with gene symbols) in RBCs were shown to be significantly downregulated in LPS birds relative to control birds (Additional files 4, 5 and 6; Fig. 1b). These numbers include genes receiving overlap in differentially expressed genes (DEGs; Fig. 1a, b).
To examine variation among replicates and tissues, principal component analysis resulted in PC1 and PC2 loadings that explained 17.6 and 10.9% of the variance in individual gene expression (measured in Fragments per Kilobase of transcript per Million mapped reads (FPKM)), respectively. Clusters separated by PC1 vs. PC2 demonstrate a separation of samples according to tissue type and treatment (Fig. 2). More specifically, PC1 delineates large variation between RBC samples and other tissues, while PC2 reflects variation among all tissue types (Fig. 2).
GO analysis of 628 upregulated DEGs in the hypothalamus showed 92 to be involved in the stress response (Additional file 7). Several of the most upregulated genes appeared in this category were radical S-adenosyl methionine domain containing 2 (RSAD2), suppressor of cytokine signaling 3 (SOCS3), connective tissue growth factor (CTGF), G protein-coupled receptor 75 (GPR75), interferon regulatory factor 1 (IRF1), and eukaryotic translation initiation factor 4E binding protein 1 (EIF4EBP1; Additional file 7, Fig. 3). In addition, macromolecule, cellular, and protein localization and protein catabolic processes were other GO terms identified (Additional file 7, Fig. 4).
Differential expression analysis from the transcriptomic response resulted in 134 downregulated genes that were successfully associated with gene symbols. GO analysis of the total downregulated DEGs indicated processes associated with metabolism: small molecule metabolic process, regulation of hormone levels, feeding behavior, and hormone metabolic process (Additional file 1, Fig. 4) that included calcitonin related polypeptide alpha (CALCA), retinol binding protein 4 (RBP4), carboxyl ester lipase (CEL), acetylserotonin O-methyltransferase (ASMT), phosphoenolpyruvate carboxykinase 1 (PCK1), apolipoprotein D (APOD), and tyrosine hydroxylase (TYH; Additional file 1, Fig. 3). Pro-melanin-concentrating hormone (PMCH), in addition to CALCA and TYH, were some of the genes present in the overrepresented function of feeding behavior (Additional file 7, Fig. 3). The defense gene, avian beta defensin 9 (AvBD9), was also downregulated.
The serine peptidase inhibitor, kazal type 4 (SPINK4) gene was the most upregulated of all genes in our global results (14.71-fold change). From GO analysis, functions associated with the 439 upregulated genes in spleen include immune system process, response to external biotic stimulus, response to other organism, response to biotic stimulus, immune response, stress response, and defense response (Additional file 7, Fig. 4). Some of the top 20 DEGs, such as ankyrin repeat and BTB domain containing 2 (ABTB2), activating transcription factor 3 (ATF3), and aldehyde dehydrogenase 1 family, member A3 (ALDH1A3), were involved in the above stress response categories. Several other DEGs were associated with more than one function, demonstrating the multi-functionality of these genes. Interleukin-1 beta (IL-1β), interleukin-8 (IL-8), C-C motif chemokine 4 homolog (CCL4) (Additional file 7, Fig. 3) are involved in the stress response, immune function, and cytokine production, while SOCS3 is involved in the stress response and cytokine production, and 2′-5′-oligoadenylate synthase-like protein 1-like (OASL) is involved in the stress response and immune function (Fig. 3).
Functions resulting from GO analysis are almost exclusively attributed to metabolic activity: small molecule, organic acid, carboxylic acid, and oxoacid metabolic processes, as well as oxidation-reduction process (Additional file 7, Fig. 4). Lecithin retinal acyltransferase (LRAT), hydroxyacid oxidase 2 (HAO2), tyrosine aminotransferase (TAT), dihydropyrimidinase (DPYS), fatty-acid amide hydrolase 1-like (FAAH), serpin family C member 1 (SERPINC1), and RBP4 are DEGs associated with metabolic process (Additional file 7, Fig. 3). Additionally, the LRAT gene experienced the largest fold-change in downregulation (− 9.51-fold change) of all DEGs across all tissues analyzed.
Red blood cells
From the GO analysis on the 121 upregulated DEGs of the RBCs, immune and cytokine function (e.g. innate immune response, leukocyte activation, cytokine production, type I interferon production), immune system development, hematopoesis, hematopoetic or lymphoid organ development, as well as stress response, was observed (Additional file 7, Fig. 4). C-C motif chemokine ligand 5 (CCL5), toll-like receptor 3 (TLR3), leucine-rich repeat kinase 1 (LRRK1), serum amyloid A like 1 (SAAL1), mitogen-activated protein kinase kinase kinase 8 (MAP3K8), and interferon induced with helicase C domain 1 (IFIH1) were some of the genes involved with immune and cytokine functions (Additional file 7, Fig. 3).
There were 61 significantly downregulated genes in the RBCs. GO analysis indicated the following functional terms: regulation of cellular component organization, receptor metabolic process, regeneration, regulation of vesicle function, and negative regulation of interleukin 8 (IL-8; Additional file 7, Fig. 4). Annexin A1 isoform p37-like (ANXA1) and arrestin beta 1 (ARRB1) were among the most downregulated genes that contribute to the negative regulation of IL-8 (Additional file 7, Fig. 3).
Five genes were significantly upregulated in all three tissues tested (Table 2). Two of the five common genes— CCL5 and RSAD2—were reported among the top 30 overrepresented functional categories. Additionally, a single gene, glutathione S-transferase class-alpha variant 2 (GSTA2), was downregulated in all three tissues. The GSTA2 gene is involved in the metabolic process according to its GO annotation (Table 2).
The aims of this study were to determine the transcriptomic response of zebra finch to LPS challenge using RNA-seq and to assess whether trade-offs between immune function and other life-history activities (growth/reproduction) could be recapitulated at the molecular level. However, gene transcripts that directly influenced growth or reproduction were not transcriptionally down-regulated in response to LPS. Rather, transcripts responsible for regulating metabolism exhibited down-regulation. Because we examined the genomic response to LPS on an acute time scale (2 h post-challenge), it is possible that genes responsible for controlling growth and/or reproduction may occur at a later time frame following immune challenge. For example, in white-crowned sparrows (Zonotrichia leucophyrs), the hypothalamic-pituitary-gonadal (HPG) axis is suppressed 24 h after LPS injection as measured by a reduction in plasma luteinizing hormone concentration .
Genes involved with the stress response were upregulated within the hypothalamus in response to LPS challenge. This response is consistent with rapid activation of the hypothalamic-pituitary-adrenal (HPA) axis that occurs during the APR, including in birds [25, 26]. Given the well-known function of the spleen in regulating immunity, it was not surprising that immune function and cytokine production genes were activated in the spleen. However, the increased expression of immune-related and cytokine production-related genes in avian RBCs indicates a novel function of nucleated RBCs. Most studies to date have assumed that RBCs play little or no role in regulating immunity, but only a few have examined non-mammalian RBCs [27, 28].
The downregulated transcripts of the hypothalamus and spleen imply a reduction in metabolism based on the Gene Ontology (GO) analysis (Additional file 7; Fig. 4). Although many endotherms exhibit fever when challenged with LPS, in small songbirds, such as the zebra finch, hypothermia results with a concurrent decrease in metabolic rate , which likely explains why genes associated with metabolism were down-regulated rather than up-regulated. It is hypothesized that in small endotherms, the heat loss from a fever response would be more detrimental than the purported advantages . Hyperthermia, experienced by most of the larger endotherms (pigs, goats, and chickens), and hypothermia in small birds, is thought to be regulated by the same inflammatory and stress signals that trigger fever . In the hypothalamus, feeding behavior and hormone-related processes were also down-regulated. A decrease in feeding behavior (anorexia) is a common occurrence observed in animals injected with LPS [8, 25].
Of the total number of transcripts differentially expressed in each tissue, only six were shared among all three, and included five upregulated genes and one downregulated gene (Fig. 1a, b; Table 2). The majority of differentially expressed genes (DEGs) were tissue specific, with the hypothalamus experiencing the largest shift in unique upregulated DEGs and spleen having the most downregulated. These data suggest that different tissues have a specialized role when responding to an immune challenge, although the processes they are involved with regulating could be similar. The most DEGs within overrepresented GO functional categories are the most informative for the study and were chosen for further examination within those functions.
The hypothalamus plays an integral role in homeostasis by regulating body temperature, appetite, reproductive behaviors, and circadian rhythms. The APR, triggered by LPS treatment, activates the HPA axis and suppresses the HPG axis via stimulation by inflammatory cytokines such as IL-1β [8, 31].
Some of the notable genes involved with stress response were RSAD2, SOCS3, and EIF4EBP1. Also known as viperin, RSAD2 is upregulated in response to H5N1 avian flu  and dengue infections  through involvement in interferon and cytokine signaling pathways. RSAD2 is upregulated in chicken kidney samples following avian infectious bronchitis virus exposure , suggesting that RSAD2 might be induced during early infection. SOCS3 and EIF4EBP1 are associated with the inhibition of interferon production via a negative feedback loop that prevents interferon and other inflammatory cytokines from inducing immunopathology . SOCS3 was upregulated by LPS challenge with a 3.462-fold change increase compared to control samples. Ranaware et al.  reported a 4.2-fold change increase following H5N1infection in chicken lung tissue. It is also possible that a virus, such as H5N1 or influenza A, may evade the host’s immune response by inducing cytokine suppressors (SOCS1 and SOCS3), that in turn inhibit interferon and toll-like receptor signaling [36, 37]. EIF4EBP1 augments the innate immune system as a translational repressor of interferon production . EIF4EBP1 acts by repressing expression of interferon regulatory factor 7 (IRF7), thereby repressing interferon production . Without EIF4EBP1 repression, IRF7 triggers interferon production, creating a positive feedback loop for further IRF7 and interferon production. IRF7 was not significantly upregulated in these samples (0.964-fold change, q-value 0.191), and interferon was not present. However, there was upregulation in IRF1 (2.975-fold change; Additional file 2). The products of IRF1 also trigger an interferon response during early infection. A moderate upregulation in IRF1 transcription was seen in chicken lung following H9N2 infection  and IRF1 knock-out mice cannot survive West Nile virus infection . In hypothalamus, IL-1β was not significantly upregulated (q-value of 0.053), but this does not discount a possible biological significance, especially when examining the fold change (5.132) and results found in the other tissues examined (fold change value 7.184 in spleen and 2.129 in RBC).
Genes associated with metabolism and cellular and protein localization were down-regulated in LPS-treated birds. The CALCA gene encodes several peptide hormones, including calcitonin and calcitonin-gene-related peptide (CGRP), that are linked to glucose and lipid metabolism  . CGRP is also involved in controlling inflammation through anti-inflammatory actions [42, 43]. In chronic inflammation disorders, such as inflammatory response syndrome in dogs, CALCA is upregulated , where it inhibits tumor necrosis factor (TNF) and interleukin 12 (IL-12) cytokines . Tsujikawa et al.  documented this anti-inflammatory effect 24 h after LPS treatment in mice, where cytokine levels were highest 3–6 h post-treatment.
Another downregulated gene that influences both metabolism and inflammation is RBP4, a negative acute phase inflammatory reactant. Decreased transcription in favor of the amino acids being used for producing positive acute phase reactants, such as IL-1β and TNF, that also inhibit RBP4 likely potentiates the APR [46,47,48]. Metabolic GO terms associated with this gene include ‘gluconeogenesis’ and ‘positive regulation of insulin secretion,’ and it has been linked with human obesity . The RBP4 protein is required for the transport and utilization of retinol (vitamin A) from the liver to peripheral tissues  where retinol derivatives are involved in lipid metabolism . However, retinol has been shown to enhance recovery from infection via actions in lymphoid tissues and its role in immune cell development (neutrophils, macrophages, natural killer cells; ). During the APR, already circulating retinol can be utilized, but less will be transported from the liver due to the downregulation of RBP4 transcription . Consequently, the immune system of malnourished individuals, with less circulating retinol before infection, may have a decreased ability to combat infection.
There were also downregulated genes involved with mediating feeding behavior. PMCH encodes melanin-concentrating hormone (MCH), with confirmed activity in appetite control by the use of MCH agonists and antagonists . Saito et al.  and Adams et al.  demonstrated the role of this hormone in olfaction, where the downregulation or deletion of PMCH impaired food-seeking behaviors. CALCA and TYH are also in this category for their metabolic cascades driving energy homeostasis by CALCA  and brain function products by TYH .
The differing roles of genes that were up- or downregulated in the hypothalamus following LPS treatment reflects a potential trade-off between functions. Several of the most upregulated genes code for proteins that influence cytokine pathways. Control of cytokine activity is evident in our results through the upregulation of the cytokine inhibitors SOCS3 and IEF4EBP1. All of the downregulated genes discussed here participate in metabolism and/or feeding behaviors. CALCA has anti-inflammatory actions [42, 43] which are depressed during the APR in this study. The reduction in feeding behaviors by the downregulation of CALCA, PMCH, and TYH may lower the risk of being a predatory target when searching for food when the animal is not in optimal health. The TYH enzyme product facilitates catecholamine synthesis (dopamine, epinephrine, norepinephrine), further affecting brain functions such as attention, memory, and cognition . The above results suggest a physiological shift towards combating illness, and a reduction in unnecessary energy expenditure and behaviors that could increase risk of predation.
One of the largest organs of the lymphatic system, the spleen filters fluids in the body, including blood. During filtration, excess blood can be stored, old RBCs removed, and iron recycled . The spleen also plays a prominent role in the peripheral immune system. When macrophages in the spleen interact with bacterial or viral components, the innate immune system is triggered. Additionally, the spleen is involved in the adaptive immune system through B cell production of antibodies and T cell activation .
In response to immune challenge by LPS, the spleen experienced an upregulation in genes involved in immune function, cytokine production, and stress response (Additional file 7, Fig. 4). This finding is consistent with the role that the spleen plays in responding to bacterial infections. The potent cytokine gene, IL-1β, is involved in the stress response, immune function, and cytokine production GO function categories. IL-1β is one of the most studied pro-inflammatory cytokines for its crucial role in initiating and regulating the immune system, and has been well documented in chickens by LPS injection , infectious bursal disease , and different bacterial infections [61, 62]. In wild birds that act as natural reservoirs for disease and encounter domesticated animals, IL-1β is less understood. Park et al.  showed that house finch splenocytes had several cytokines significantly induced by IL-1β, including interleukin 2 (IL-2), interleukin 10 (IL-10), and chemokine C-X-C motif ligand 1 (CXCL1).
Two chemokines, IL-8 and CCL4, were also upregulated. IL-8 products induce lysosomal enzyme release from neutrophils and chemotactic activity for basophils and adaptive T cells . IL-8 gene is greatly stimulated by IL-1β, but not by interferon  and is elevated in respiratory distress syndrome in humans . CCL4 also encodes products with chemotactic activity. CCL4 attracts natural killer cells, monocytes, and in later adaptive responses, T cells . Of interest is the change in expression levels of CCL4 with H5N1 infection. Compared to seasonal influenza viruses, such as H1N1 and H3N2, CCL4 is more strongly induced by H5N1 . Cheung et al.  even reported the downregulation of CCL4 in human samples with H1N1 infection. This relationship of H5N1 and CCL4 may have additional implications for wild and domestic avian species.
A gene involved with the stress response in spleen in response to LPS was ALDH1A3. This gene produces the rate limiting dehydrogenase enzyme for converting retinol (vitamin A) to retinoic acid, a hormone-like metabolite with the ability to modulate immune responses by promoting inflammation . Retinoic acid activates and increases cytokine production in the innate immune system and is used by dendritic, B, and T cells in the adaptive immune system . Increased conversion of retinol by the upregulation of ALDH1A3 also decreases the limited serum retinol resulting from the downregulation in transcription of its transporter, RBP4. Without this transporter, retinol remains stored in the liver and unusable in other parts of the body . The limit on retinol is further compounded by the fact RBP4 is downregulated in both the hypothalamus and spleen.
SOCS3, involved in both the stress response and cytokine production, was upregulated in the hypothalamus and spleen (Additional file 7, Fig. 3). SOCS3 helps regulate the immune system response by preventing the induction of immunopathology from cytokines . Cytokine genes were among the most upregulated in spleen compared to the hypothalamus, suggesting that increased expression of SOCS3 facilitates the dampening of splenic inflammatory responses from cytokines.
LPS induced a down-regulation of splenic genes that regulate metabolism. One of these, LRAT, is involved in retinol metabolism as discussed above. Specifically, it encodes the predominant enzyme for esterification of retinol [72, 73]. LRAT negatively regulates retinoic acid biosynthesis . The downregulation of LRAT, along with the upregulation of ALDH1A3 (see above), can result in increased biosynthesis of retinoic acid. The increased presence of retinoic acid further suggests the importance of retinoic acid in inflammation. As in the hypothalamus, RBP4 was downregulated in the spleen.
Among other downregulated genes involved in metabolism, HAO2 encodes an enzyme that mediates the oxidation of fatty acids [75, 76]. Likewise, the hydrolase encoding gene, FAAH, is involved in the hydrolysis of 2-arachidonoylglycerol, a lipid that stimulates phagocytic activities in macrophages during a defense response . Also, by the downregulation of TAT, the amino acid tyrosine is not catabolized [78, 79], possibly making the amino acid further available for cytokine production. Another downregulated gene, DPYS, encodes a rate-limiting enzyme that plays a role in pyrimidine catabolism . Inhibiting these genes reduces degradation of materials needed by the immune response, such as amino acids and immune system mediators.
In spleen, there was polarized expression of immune system genes and metabolic genes as in hypothalamus. The spleen shared the upregulation in the cytokine suppressor, SOCS3, with the hypothalamus. RBP4, the negative acute phase reactant, was a downregulated gene shared with the hypothalamus, in the metabolic functional category. From specific functions of genes in the metabolic category, there was a reduction in catabolic processes of fatty acids, amino acids, and pyrimidines. The most differentially expressed gene in either direction was also found within the spleen: SPINK4 (14.710-fold change) and LRAT (− 9.510-fold change). The downregulation of LRAT, which normally inhibits retinoic acid biosynthesis, along with the upregulation of ALDH1A3, further suggests the importance of retinoic acid to the immune system. The known functions of the spleen, compared to the hypothalamus, and the absence of behavioral-based genes is expected (Fig. 4).
Red blood cells
In contrast to mammalian RBCs, avian, reptilian, and amphibian RBCs contain a nucleus . The best reported defense that anucleated RBCs have against pathogens is the release of reactive oxygen species from hemoglobin when the cell is lysed, breaking down lipids, proteins, and DNA of nearby pathogens . Little attention has been given to the function of nucleated RBCs. It has been assumed these nucleated RBCs merely participate in gas exchange, with little or no protein synthesis, like mammalian anucleated RBCs . Studies using avian whole blood show an immune response [84, 85] but the results cannot differentiate between RBC and leukocyte activity.
More recently, a few studies examined the activity of nucleated RBCs in immune function. There is now evidence of active transcriptional machinery within these cells that do react to various stimuli. A differential response of genes involved in the immune system and metabolic processes was shown in trout and chicken cultured RBCs . To our knowledge this is the first study to use RNA-seq for quantitative analysis of nucleated RBCs to immune challenge in a wild avian species.
The reaction of the nucleated RBCs in the zebra finch to LPS involve genes associated with the immune response. Genes that regulate inflammation, cytokine production, and the stress response were upregulated within top DEGs (Additional file 7, Fig. 4). One of these, the chemokine CCL5, is expressed early in an immune response, activated through signaling by the presence of pro-inflammatory cytokines and by contact with pathogens [86, 87]. Elevated expression of CCL5 has been noted in hantavirus, reovirus, adenovirus, and influenza virus (H5N1) infections, but severe acute respiratory syndrome (SARS) coronavirus is an especially strong inducer in human lung tissue .
Three upregulated genes are involved in a cascade of signaling: MAP3K8, TLR3, and IFIH1. MAP3K8, activated by IL-1β, is critically involved in TLR signaling . TLR3 propagates signals after being induced by interferons and cytokines and is involved in RIG-I and MDA5 pathways . Upregulation of TLR and IFIH1 (MDA5) is documented in chickens infected with H5N1, low pathogenic avian influenza (LPAIV) H7N1, and avian Tembusu virus (ATMUV) that cause great economic loss and present a danger to humans . Once activated, the gene encoding MDA5 initiates further signal transduction for cytokine secretion. MDA5 is the primary influenza A virus sensor for chickens . The upregulation of TLR3 is one gene that can be used to assess RIG-I or MDA5 expression, depending on the avian species.
Few functional categories were produced by GO analysis of the downregulated RBC genes, and few genes represented these functions. Out of 61 downregulated genes, the functional category of regulation of cellular component organization only contained 12 genes from the DEG list. Given that many of the genes in zebra finch are not annotated thoroughly, this could lend to a lack in functional categorization.
Two genes were involved in the negative regulation of IL-8. IL-8 has chemotactic activities for basophils, as well as stimulating phospholipase enzymes that interact with neutrophils [64, 91]. Both downregulated genes in this category have anti-inflammatory activity. ANXA1 inhibits phospholipase, interfering with IL-8 downstream mediation of inflammation. ARRB1 inhibits upstream processes in TLR signaling that lead to IL-8 expression [92, 93].
These novel findings show that there are transcriptomic adjustments during the APR of the immune response. In fact, the six most upregulated DEGs are involved in the immune system and/or cytokine production. The upregulation of CCL5 in this these cells may be advantageous because of the gene product’s role in leukocyte recruitment. Despite there being few downregulated genes represented in the GO functional categories, ANXA1 and ARRB1 are involved in the negative regulation of IL-8. With dampened expression of ANXA1 and ARRB1, IL-8 involved signaling is not inhibited. More can be elucidated about the complete transcriptome changes from immune challenge by examining the other 115 upregulated and 58 downregulated genes.
Of the five genes shared between all three tissues, only two contributed to the functions produced by GO analysis: CCL5 and RSAD. While not within the top most upregulated in all tissues, they were both statistically significantly expressed. The first of these genes, CCL5, as discussed, is a chemokine that regulates leukocyte trafficking . It has also been highly conserved during evolution, with sequence homology in mammals, fish, and birds . The second, RSAD2, regulates interferon and cytokine signaling pathways. The exclusion of the other genes was based on their associated GO terms.
The single downregulated gene that was shared, GSTA2, has the biological process GO term: ‘metabolic process.’ The encoded enzyme is reported to be involved in cellular detoxification and excretion of several xenobiotic substances, with strong implications for ingestion of mycotoxin aflatoxin B1 in turkeys . Turkeys are especially susceptible to the toxin because GSTA2 and the enzyme product are unresponsive to aflatoxin B1. IL-6 was significantly upregulated in the spleen compared to controls (1.834-fold change, q-value 0.030). This is further explained by GSTA2’s anti-inflammatory properties and activation by glucocorticoids. Interestingly, GSTA2 is associated with plumage dichromatism; GSTA2 and APOD are involved with carotenoid uptake, binding, and deposition in pheasants . In this study, APOD was downregulated in the hypothalamus (− 3.050-fold change, q-value 0.002). The role of GSTA2 in carotenoid processes suggests a downregulation much like with RBP4 and LRAT, by reducing interference in the use of retinoic acid.
Given these results we have a better understanding of the global transcriptomic changes in the hypothalamus, spleen, and nucleated RBCs of zebra finch following activation of the APR by bacterial LPS. This study provides insight into physiological changes during the APR, with many genes involved in the immune system being upregulated and genes in metabolic pathways being downregulated. Because the environment of free-living organisms is widely variable in resources and conditions, it is important to understand these potential trade-offs and implications on seasonally important functions, such as migration, reproduction, and growth. For example, parental care has been shown to decrease in house sparrows (Passer domesticus) during illness , as well as LPS injection reducing territorial aggression in white-crowned sparrows .
White blood cells
Comparing transcriptome changes of RBCs and WBCs will give further insight into the APR of zebra finch. However, complications with obtaining an appropriate amount of RNA for sequencing prevented that analysis in this study. Although RNA-seq requires much less RNA than other methods , the concentration of WBCs in whole blood is low. Combined with the size of the zebra finch, other methods or techniques will be needed to accomplish this in the future.
Further validation of these results can be done through RT-PCR. This would be best accomplished with multiple gene primers based on the sheer number of DEGs in our data. Custom primer design is recommended.
This study provides insight into the rapid molecular responses that occur during immune system activation in a passerine bird species. The APR is an evolutionarily conserved rapid response to eliminate and control infection, and within 2 hours of LPS administration there is upregulation of immune genes and a downregulation in genes putatively involved with metabolic pathways in the hypothalamus and spleen. The upregulation in immune related genes in the hypothalamus, spleen, and red blood cells supports part of the hypothesized transcriptomic response to immune challenge in zebra finch. However, the current results can only suggest a trade-off in immune function and other important aspects of life. Further research is needed to examine this potential trade-off using experimental methods.
By using LPS, a component of gram-negative bacteria, pathogen manipulation is avoided. Tissue specificity is also shown by the DEGs not shared between and among tissues (Figs. 1 and 2). Of the total DEGs, only five upregulated and one downregulated gene were shared. We additionally present novel insight into transcriptional changes in nucleated RBCs of a non-model avian organism. These findings on the immunological response of avian RBCs will lay a framework for further investigations into the function of nucleated RBCs, a topic that to date has been poorly studied.
Male zebra finch (Taeniopygia guttata) taken from our breeding colony were housed in individual cages (16.5 × 11.8 × 22, Petsmart Co.) in the animal facility at Western Kentucky University for 7 days. Food (Finch and Canary Breeding and Molting Seed Blend, Lady Gouldian Finch, Irvine, CA; lettuce/cabbage leaves) and water were provided ad libitum. Animals were exposed to a 12 h light: 12 h dark photoperiod (lights on at 0700) and an ambient temperature of 23 ± 1 °C. Between 4 and 5 h after lights on, birds were quickly captured (under 30 s) and given an i.p. injection of either bacterial lipopolysaccharide (LPS; 2 mg/kg of BW; n = 8 Escherichia coli; serotype 026:B6; Sigma-L8274) dissolved in 0.9% saline (vehicle) or vehicle alone (n = 8). Two hours later, animals were euthanized using isoflurane until unconscious (< 20 s) and then rapidly decapitated. Brain and spleen were collected from half of the treatment group (n = 4) and half of the control group (n = 4) and immediately placed in RNAlater (AM7021; Thermo Fischer Scientific). Tissues were stored in RNAlater at − 20 °C until processed. The hypothalamus was dissected out of the brain before RNA isolation (see below). Whole trunk blood was collected from the remaining four treatment and four control birds and placed into tubes containing EDTA and kept on ice until processed (see “Blood” below).
EDTA-treated blood was layered on top of a single step density gradient medium (PolymorphPrep; Axis-Shield) consisting of 1-part blood: 2-parts gradient medium and then further processed following the manufacturer’s protocol. Layers of peripheral blood mononuclear cells (PBMC), polymorphonuclear leukocytes (PMN), and white blood cells (WBCs) were extracted and pooled per individual animal. Slides were prepared from erythrocyte (red blood cells; RBCs) and PBMC/PMN samples and stained using Hemacolor Rapid (Millipore Sigma) following the manufacturer’s protocol. Each slide was surveyed under 400 x magnification for 10 min to assess for contamination before RNA isolation.
RNA was isolated from PBMC/PMNs and RBCs using Trizol (Ambion) lysing and a RNeasy Mini Kit (Qiagen). The hypothalamus and spleen were homogenized and RNA was extracted using a RNeasy Mini Kit (Qiagen) following the manufacturer’s instructions. RNA yield and purity were assessed using a NanoDrop N-D 2000 spectrophotometer (Thermo Fischer Scientific) and an Agilent 2100 Bioanalyzer (Agilent), respectively. Total RNA concentrations ranged from 10.0 to 53.6 ng/μl for hypothalami, 17.2 to 77.2 ng/μl for spleen, and 4.6 to 20.8 ng/μl for RBCs. RNA Integrity Numbers (RIN) were > 7.5. Unfortunately, PBMC/PMNs yielded RNA concentrations that were too low for sequencing purposes. Low RNA yield can be attributed to the low concentration of WBCs in whole blood and the small volume of whole blood we were able to obtain from zebra finch (typically < 150 μl).
Library preparation and sequencing
Samples were analyzed at the University of Louisville Genomics Core (Louisville, KY) for library preparation and sequencing. SMART-Seq v4 Ultra Low Input RNA kits (Takara Bio USA, Inc.) were used to prepare cDNA libraries. Samples were barcoded with Illumina TruSeq Adapters. After library clean-up using Agencourt AMPure XP Beads, quality was assessed on an Agilent Bioanalyzer using the Agilent DNA 1000 Kit. This confirmed the final fragment size for all samples to be approximately 400 bp, as expected from the protocol. Samples were sequenced twice on an Illumina NextSeq 500 sequencer, with four biological replicates and four lanes per replicate.
Bioinformatics analysis was performed at the Kentucky Biomedical Research Infrastructure Network (KBRIN) Bioinformatics Core at the University of Louisville Genomics Core. Raw sequencing files were downloaded from Illumina’s BaseSpace. Quality scores for raw sequences were sufficiently high, above the recommended score of 20 (base call accuracy of 99%), to proceed with analysis [97, 98]. The Tuxedo Suite pipeline was used for analysis of data . High quality single-end reads (152 million for hypothalamus, 140 million for spleen, and 326 million for RBCs) were successfully mapped to the zebra finch reference genome (taeGut18.104.22.168) from Ensembl (http://www.ensembl.org; ) using the TopHat2 (v. 2.0.13) tool. Differential expression between LPS and saline conditions was determined from fragments per kilobase million (FPKM) normalized read counts using Cuffdiff2 (v. 2.2.1) of the Tuxedo Suite programs and log base 2 transformed. The replicates were handled as a group (LPS-injected birds were treated as one group and saline-injected birds in another group). The Cuffdiff2 program calculates the mean value and variance as it assigns p-values and adjusted p-values. All fold change expression values presented from this study are provided as log base 2. Genes with Benjamini-Hochberg false discovery rate (FDR) (q-value) ≤ 0.05 and a log base 2 transformed fold change ≥1 were considered significantly differentially expressed.
To assess variation within groups and among tissues, we conducted a principal component analysis using normalized read counts (FPKM) for 18,168 tested genes locations as input using the R statistical function ‘prcomp.’
Gene ontology analysis
Gene annotations came from Ensembl BioMarts as part of the KBRIN Bioinformatics Core analysis [101,102,103]. Transcripts with differential expression from sequence code but no known gene ID were excluded from further analysis. Gene ontology (GO) enrichment analysis was performed on differentially expressed genes (DEGs) using the ShinyGO web-based tool (http://bioinformatics.sdstate.edu/go/) (v. 0.6). GO functional categories falling under biological processes, with a q-value ≤0.05 following hypergeometric testing, were considered significantly overrepresented. In the hierarchy of GO, a gene can be represented in more than one category because of the functional versatility of genes, but only once within each category. Likewise, genes appear in the categories represented by both the specific (child term) and broad (parent term) categories under these biological processes.
Availability of data and materials
Raw sequencing data and analyses for this study are available at GEO (GSE121348): https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE121348.
Ankyrin repeat and BTB domain containing 2
Aldehyde dehydrogenase 1 family, member A
Annexin A1 isoform p37-like
Acute phase response
Arrestin beta 1
Activating transcription factor 3
Avian Tembusu virus
Avian beta defensin 9
Calcitonin related polypeptide alpha
C-C motif chemokine 4 homolog
Carboxyl ester lipase
Connective tissue growth factor
Differentially expressed genes
DEXH (Asp-Glu-X-His) box polypeptide 58
Eukaryotic translation initiation factor 4E binding protein 1
Fatty-acid amide hydrolase 1-like
False discovery rate
Fragments per Kilobase of transcript per Million mapped reads
G protein-coupled receptor 75
Glutathione S-transferase class-alpha variant 2
Influenza A virus subtype H1N1
Influenza A virus subtype H3N2
Influenza A virus subtype H5N1
Influenza A virus subtype H7N1
Hydroxyacid oxidase 2
Interferon induced with helicase C domain 1
Interferon regulatory factor 1
Transcription factor AP-1
Laboratory of genetics and physiology 2
Low pathogenic avian influenza
Lecithin retinal acyltransferase
Leucine-rich repeat kinase 1
Lymphocyte antigen 86 gene
mitogen-activated protein kinase kinase kinase 8
Melanin concentrating hormone
Melanoma differentiation-associated gene 5
NEWCASTLE disease virus
2′-5′-oligoadenylate synthase-like protein 1-like
Peripheral blood mononuclear cells
Principal component analysis
Phosphoenolpyruvate carboxykinase 1
Red blood cells
Retinol binding protein 4
Retinoic acid-inducible gene I
Radical S-adenosyl methionine domain containing 2
Reverse transcription polymerase chain reaction
Serum amyloid A like 1
Severe acute respiratory syndrome
Serpin family C member 1
Suppressor of cytokine signaling 3
Synoviocyte proliferation-associated in collagen-induced arthritis 1
Serine peptidase inhibitor, kazal type 4
Toll-like receptor 3
Tumor necrosis factor-alpha
White blood cells
West Nile virus
Ebling FJ. On the value of seasonal mammals for identifying mechanisms underlying the control of food intake and body weight. Horm Behav. 2014;66(1):56–65.
Dawson A, King VM, Bentley GE, Ball GF. Photoperiodic control of seasonality in birds. J Biol Rhythm. 2001;16:365–80.
Wingfield JC. Organization of vertebrate annual cycles: implications for control mechanisms. Philos Trans R Soc B. 2008;363(1490):425–41.
Norris K, Evans MR. Ecological immunology: life history trade-offs and immune defense in birds. Behav Ecol. 2000;11(1):19–26.
Sheldon BC, Verhulst S. Ecological immunology: costly parasite defences and trade-offs in evolutionary ecology. Trends Ecol Evol. 1996;11(8):317–21.
Ricklefs RE, Wikelski M. The physiology/life-history nexus. Trends Ecol Evol. 2002;17:462–8.
Demas GE, Zysling DA, Beechler BR, Muehlenbein MP, French SS. Beyond phytohaemagglutinin: assessing vertebrate immune function across ecological contexts. J Anim Ecol. 2011;80(4):710–30.
Owen-Ashley NT, Wingfield JC. Acute phase responses of passerines: characterization and seasonal variation. J Ornithol. 2007;148(Suppl 2):S583–S91.
Suffredini AF, Fantuzzi G, Badolato R, Oppenheim JJ, O'Grady NP. New insights into the biology of the acute phase response. J Clin Immunol. 1999;19(4):203–14.
Ashley NT, Wingfield JC. Sickness behavior in vertebrates: Allostasis, life-history modulation, and hormonal regulation. In: Demas GE, Nelson RJ, editors. EcoImmunology. Oxford: Oxford University Press; 2012. p. 45–91.
Hart BL. Biological basis of the behavior of sick animals. Neurosci Biobehav Rev. 1988;12:123–37.
Ashley NT, Weil ZM, Nelson RJ. Inflammation: mechanisms, costs, and natural variation. Annu Rev Ecol Evol Syst. 2012;42:385–406.
Medzhitov R. Origin and physiological roles of inflammation. Nature. 2008;454(7203):428–35.
Skold-Chiriac S, Nord A, Nilsson JA, Hasselquist D. Physiological and behavioral responses to an acute-phase response in zebra finches: immediate and short-term effects. Physiol Biochem Zool. 2014;87(2):288–98.
Hedrick SM. The acquired immune system: a vantage from beneath. Immunity. 2004;21(5):607–15.
Janeway CA Jr, Medzhitov R. Innate immune recognition. Annu Rev Immunol. 2002;20:197–216.
Bayne CJ, Gerwick L. The acute phase response and innate immunity of fish. Dev Comp Immunol. 2001;25(8–9):725–43.
Cray C, Zaias J, Altman NH. Acute phase response in animals: a review. Comp Med. 2009;59(6):517–26.
Demas GE, Nelson RJ. Ecoimmunology. Oxford: Oxford University Press; 2012.
Demas GE. The energetics of immunity: a neuroendocrine link between energy balance and immune function. Horm Behav. 2004;45:173–80.
Wang Z, Gerstein M, Snyder M. RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev Genet. 2009;10(1):57–63.
Ekblom R, French L, Slate J, Burke T. Evolutionary analysis and expression profiling of zebra finch immune genes. Genome Biol Evol. 2010;2:781–90.
Warren WC, Clayton DF, Ellegren H, Arnold AP, Hillier LW, Kunstner A, et al. The genome of a songbird. Nature. 2010;464(7289):757–62.
Newhouse DJ, Hofmeister EK, Balakrishnan CN. Transcriptional response to West Nile virus infection in the zebra finch (Taeniopygia guttata). R Soc Open Sci. 2017;4(6):170296.
Owen-Ashley NT, Turner M, Hahn TP, Wingfield JC. Hormonal, behavioral, and thermoregulatory responses to bacterial lipopolysaccharide in captive and free-living white-crowned sparrows (Zonotrichia leucophrys gambelii). Horm Behav. 2006;49:15–29.
Ashley NT, Demas GE. Neuroendocrine-immune circuits, phenotypes, and interactions. Horm Behav. 2017;87:25–34.
Morera D, MacKenzie SA. Is there a direct role for erythrocytes in the immune response? Vet Res. 2011;42:89.
Morera D, Roher N, Ribas L, Balasch JC, Donate C, Callol A, et al. RNA-Seq reveals an integrated immune response in nucleated erythrocytes. PLoS One. 2011;6(10):e26998.
Jones CA, Edens FW, Denbow DM. Influence of age on the temperature response of chickens to Escherichia coli and Salmonella typhimurium endotoxins. Poult Sci. 1983;62:1553–8.
Romanovsky AA, Almeida MC, Aronoff DM, Ivanov AI, Konsman JP, Steiner AA, et al. Fever and hypothermia in systemic inflammation: recent discoveries and revisions. Front Biosci. 2005;10:2193–216.
Beishuizen A, Thijs LG. Endotoxin and the hypothalamo-pituitary-adrenal (HPA) axis. J Endotoxin Res. 2003;9(1):3–24.
Ranaware PB, Mishra A, Vijayakumar P, Gandhale PN, Kumar H, Kulkarni DD, et al. Genome wide host gene expression analysis in chicken lungs infected with avian influenza viruses. PLoS One. 2016;11(4):e0153671.
Schoggins JW. Interferon-stimulated genes: roles in viral pathogenesis. Curr Opin Virol. 2014;6:40–6.
Cong F, Liu X, Han Z, Shao Y, Kong X, Liu S. Transcriptome analysis of chicken kidney tissues following coronavirus avian infectious bronchitis virus infection. BMC Genomics. 2013;14:743.
Chaves de Souza JA, Nogueira AV, Chaves de Souza PP, Kim YJ, Silva Lobo C, Pimentel Lopes de Oliveira GJ, et al. SOCS3 expression correlates with severity of inflammation, expression of proinflammatory cytokines, and activation of STAT3 and p38 MAPK in LPS-induced inflammation in vivo. Mediat Inflamm. 2013;2013:650812.
Pothlichet J, Chignard M, Si-Tahar M. Cutting edge: innate immune response triggered by influenza A virus is negatively regulated by SOCS1 and SOCS3 through a RIG-I/IFNAR1-dependent pathway. J Immunol. 2008;180(4):2034–8.
Smith J, Speed D, Law AS, Glass EJ, Burt DW. In-silico identification of chicken immune-related genes. Immunogenetics. 2004;56(2):122–33.
Wick G, Backovic A, Rabensteiner E, Plank N, Schwentner C, Sgonc R. The immunology of fibrosis: innate and adaptive responses. Trends Immunol. 2010;31(3):110–9.
Colina R, Costa-Mattioli M, Dowling RJ, Jaramillo M, Tai LH, Breitbach CJ, et al. Translational control of the innate immune response through IRF-7. Nature. 2008;452(7185):323–8.
Brien JD, Daffis S, Lazear HM, Cho H, Suthar MS, Gale M Jr, et al. Interferon regulatory factor-1 (IRF-1) shapes both innate and CD8(+) T cell immune responses against West Nile virus infection. PLoS Pathog. 2011;7(9):e1002230.
Bartelt A, Jeschke A, Muller B, Gaziano I, Morales M, Yorgan T, et al. Differential effects of Calca-derived peptides in male mice with diet-induced obesity. PLoS One. 2017;12(6):e0180547.
Assas BM, Pennock JI, Miyan JA. Calcitonin gene-related peptide is a key neurotransmitter in the neuro-immune axis. Front Neurosci. 2014;8:23.
Russell FA, King R, Smillie SJ, Kodji X, Brain SD. Calcitonin gene-related peptide: physiology and pathophysiology. Physiol Rev. 2014;94(4):1099–142.
Giunti M, Peli A, Battilani M, Zacchini S, Militerno G, Otto CM. Evaluation of CALC-I gene (CALCA) expression in tissues of dogs with signs of the systemic inflammatory response syndrome. J Vet Emerg Crit Care (San Antonio, Tex : 2001). 2010;20(5):523–7.
Tsujikawa K, Yayama K, Hayashi T, Matsushita H, Yamaguchi T, Shigeno T, et al. Hypertension and dysregulated proinflammatory cytokine production in receptor activity-modifying protein 1-deficient mice. Proc Natl Acad Sci U S A. 2007;104(42):16702–7.
Moraes-Vieira PM, Yore MM, Dwyer PM, Syed I, Aryal P, Kahn BB. RBP4 activates antigen-presenting cells, leading to adipose tissue inflammation and systemic insulin resistance. Cell Metab. 2014;19(3):512–26.
Rosales FJ, Ritter SJ, Zolfaghari R, Smith JE, Ross AC. Effects of acute inflammation on plasma retinol, retinol-binding protein, and its mRNA in the liver and kidneys of vitamin A-sufficient rats. J Lipid Res. 1996;37(5):962–71.
Zabetian-Targhi F, Mahmoudi MJ, Rezaei N, Mahmoudi M. Retinol binding protein 4 in relation to diet, inflammation, immunity, and cardiovascular diseases. Adv Nutr (Bethesda, Md). 2015;6(6):748–62.
Kotnik P, Keuper M, Wabitsch M, Fischer-Posovszky P. Interleukin-1beta downregulates RBP4 secretion in human adipocytes. PLoS One. 2013;8(2):e57796.
Racke MK, Burnett D, Pak SH, Albert PS, Cannella B, Raine CS, et al. Retinoid treatment of experimental allergic encephalomyelitis. IL-4 production correlates with improved disease course. J Immunol (Baltimore, Md : 1950). 1995;154(1):450–8.
Klor HU, Weizel A, Augustin M, Diepgen TL, Elsner P, Homey B, et al. The impact of oral vitamin A derivatives on lipid metabolism - What recommendations can be derived for dealing with this issue in the daily dermatological practice? J Dtsch Dermatol Ges Journal of the German Society of Dermatology : JDDG. 2011;9(8):600–6.
Stephensen CB. Vitamin a, infection, and immune function. Annu Rev Nutr. 2001;21:167–92.
Naufahu J, Cunliffe AD, Murray JF. The roles of melanin-concentrating hormone in energy balance and reproductive function: Are they connected? Reproduction (Cambridge, England). 2013;146(5):R141–50.
Saito Y, Cheng M, Leslie FM, Civelli O. Expression of the melanin-concentrating hormone (MCH) receptor mRNA in the rat brain. J Comp Neurol. 2001;435(1):26–40.
Adams AC, Domouzoglou EM, Chee MJ, Segal-Lieberman G, Pissios P, Maratos-Flier E. Ablation of the hypothalamic neuropeptide melanin concentrating hormone is associated with behavioral abnormalities that reflect impaired olfactory integration. Behav Brain Res. 2011;224(1):195–200.
Daubner SC, Le T, Wang S. Tyrosine hydroxylase and regulation of dopamine synthesis. Arch Biochem Biophys. 2011;508(1):1–12.
Mebius RE, Kraal G. Structure and function of the spleen. Nat Rev Immunol. 2005;5(8):606–16.
Balazs M, Martin F, Zhou T, Kearney J. Blood dendritic cells interact with splenic marginal zone B cells to initiate T-independent immune responses. Immunity. 2002;17(3):341–52.
Weining KC, Sick C, Kaspers B, Staeheli P. A chicken homolog of mammalian interleukin-1 beta: cDNA cloning and purification of active recombinant protein. Eur J Biochem. 1998;258(3):994–1000.
Heggen CL, Qureshi MA, Edens FW, Barnes HJ. Alterations in macrophage-produced cytokines and nitrite associated with poult enteritis and mortality syndrome. Avian Dis. 2000;44(1):59–65.
Kogut MH, He H, Kaiser P. Lipopolysaccharide binding protein/CD14/ TLR4-dependent recognition of salmonella LPS induces the functional activation of chicken heterophils and up-regulation of pro-inflammatory cytokine and chemokine gene expression in these cells. Anim Biotechnol. 2005;16(2):165–81.
Lavric M, Maughan MN, Bliss TW, Dohms JE, Bencina D, Keeler CL Jr, et al. Gene expression modulation in chicken macrophages exposed to mycoplasma synoviae or Escherichia coli. Vet Microbiol. 2008;126(1–3):111–21.
Park M, Kim S, Adelman JS, Leon AE, Hawley DM, Dalloul RA. Identification and functional characterization of the house finch interleukin-1beta. Dev Comp Immunol. 2017;69:41–50.
Mukaida N, Harada A, Matsushima K. Interleukin-8 (IL-8) and monocyte chemotactic and activating factor (MCAF/MCP-1), chemokines essentially involved in inflammatory and immune reactions. Cytokine Growth Factor Rev. 1998;9(1):9–23.
John M, Au BT, Jose PJ, Lim S, Saunders M, Barnes PJ, et al. Expression and release of interleukin-8 by human airway smooth muscle cells: inhibition by Th-2 cytokines and corticosteroids. Am J Respir Cell Mol Biol. 1998;18(1):84–90.
Harada A, Sekido N, Akahoshi T, Wada T, Mukaida N, Matsushima K. Essential involvement of interleukin-8 (IL-8) in acute inflammation. J Leukoc Biol. 1994;56(5):559–64.
Cheung R, Malik M, Ravyn V, Tomkowicz B, Ptasznik A, Collman RG. An arrestin-dependent multi-kinase signaling complex mediates MIP-1beta/CCL4 signaling and chemotaxis of primary human macrophages. J Leukoc Biol. 2009;86(4):833–45.
Peiris JS, Cheung CY, Leung CY, Nicholls JM. Innate immune responses to influenza a H5N1: friend or foe? Trends Immunol. 2009;30(12):574–84.
Cheung CY, Poon LL, Lau AS, Luk W, Lau YL, Shortridge KF, et al. Induction of proinflammatory cytokines in human macrophages by influenza A (H5N1) viruses: a mechanism for the unusual severity of human disease? Lancet (London, England). 2002;360(9348):1831–7.
Broadhurst MJ, Leung JM, Lim KC, Girgis NM, Gundra UM, Fallon PG, et al. Upregulation of retinal dehydrogenase 2 in alternatively activated macrophages during retinoid-dependent type-2 immunity to helminth infection in mice. PLoS Pathog. 2012;8(8):e1002883.
Pino-Lagos K, Benson MJ, Noelle RJ. Retinoic acid in the immune system. Ann N Y Acad Sci. 2008;1143:170–87.
O'Byrne SM, Wongsiriroj N, Libien J, Vogel S, Goldberg IJ, Baehr W, et al. Retinoid absorption and storage is impaired in mice lacking lecithin:retinol acyltransferase (LRAT). J Biol Chem. 2005;280(42):35647–57.
Zhao W, Wang D, Zhao J, Zhao W. Bioinformatic analysis of retinal gene function and expression in diabetic rats. Exp Ther Med. 2017;14(3):2485–92.
Ghosh R. Lecithin:retinol acyltransferase and retinyl esters: is balance the essence in carcinogenesis? Cancer Biol Ther. 2009;8(13):1226–7.
Barawkar DA, Meru A, Bandyopadhyay A, Banerjee A, Deshpande AM, Athare C, et al. Potent and selective inhibitors of long chain l-2-Hydroxy acid oxidase reduced blood pressure in DOCA salt-treated rats. ACS Med Chem Lett. 2011;2(12):919–23.
Jones JM, Morrell JC, Gould SJ. Identification and characterization of HAOX1, HAOX2, and HAOX3, three human peroxisomal 2-hydroxy acid oxidases. J Biol Chem. 2000;275(17):12590–7.
Lee JH, Hou X, Kummari E, Borazjani A, Edelmann MJ, Ross MK. Endocannabinoid hydrolases in avian HD11 macrophages identified by chemoproteomics: inactivation by small-molecule inhibitors and pathogen-induced downregulation of their activity. Mol Cell Biochem. 2018;444(1–2):125–41.
Grossman A, Mavrides C. Studies on the regulation of tyrosine aminotransferase in rats. J Biol Chem. 1967;242(7):1398–405.
Iynedjian PB, Auberger P, Guigoz Y, Le Cam A. Pretranslational regulation of tyrosine aminotransferase and phosphoenolpyruvate carboxykinase (GTP) synthesis by glucagon and dexamethasone in adult rat hepatocytes. Biochem J. 1985;225(1):77–84.
Kikugawa M, Kaneko M, Fujimoto-Sakata S, Maeda M, Kawasaki K, Takagi T, et al. Purification, characterization and inhibition of dihydropyrimidinase from rat liver. Eur J Biochem. 1994;219(1–2):393–9.
Zhang ZW, Cheng J, Xu F, Chen YE, Du JB, Yuan M, et al. Red blood cell extrudes nucleus and mitochondria against oxidative stress. IUBMB Life. 2011;63(7):560–5.
Jiang N, Tan NS, Ho B, Ding JL. Respiratory protein-generated reactive oxygen species as an antimicrobial strategy. Nat Immunol. 2007;8(10):1114–22.
Heegaard ED, Brown KE. Human parvovirus B19. Clin Microbiol Rev. 2002;15(3):485–505.
Videvall E, Cornwallis CK, Palinauskas V, Valkiunas G, Hellgren O. The avian transcriptome response to malaria infection. Mol Biol Evol. 2015;32(5):1255–67.
Watson H, Videvall E, Andersson MN, Isaksson C. Transcriptome analysis of a wild bird reveals physiological responses to the urban environment. Sci Rep. 2017;7:44180.
Cartier L, Hartley O, Dubois-Dauphin M, Krause KH. Chemokine receptors in the central nervous system: role in brain inflammation and neurodegenerative diseases. Brain Res Brain Res Rev. 2005;48(1):16–42.
Marques RE, Guabiraba R, Russo RC, Teixeira MM. Targeting CCL5 in inflammation. Expert Opin Ther Targets. 2013;17(12):1439–60.
Mielke LA, Elkins KL, Wei L, Starr R, Tsichlis PN, O'Shea JJ, et al. Tumor progression locus 2 (Map 3k8) is critical for host defense against Listeria monocytogenes and IL-1 beta production. J Immunol (Baltimore, Md : 1950). 2009;183(12):7984–93.
Chen S, Cheng A, Wang M. Innate sensing of viruses by pattern recognition receptors in birds. Vet Res. 2013;44:82.
Chen S, Luo G, Yang Z, Lin S, Chen S, Wang S, et al. Avian Tembusu virus infection effectively triggers host innate immune response through MDA5 and TLR3-dependent signaling pathways. Vet Res. 2016;47(1):74.
L'Heureux GP, Bourgoin S, Jean N, McColl SR, Naccache PH. Diverging signal transduction pathways activated by interleukin-8 and related chemokines in human neutrophils: interleukin-8, but not NAP-2 or GRO alpha, stimulates phospholipase D activity. Blood. 1995;85(2):522–31.
Anjum SG, Xu W, Nikkholgh N, Basu S, Nie Y, Thomas M, et al. Regulation of toll signaling and inflammation by beta-arrestin and the SUMO protease Ulp1. Genetics. 2013;195(4):1307–17.
Drewniak A, Tool AT, Geissler J, van Bruggen R, van den Berg TK, Kuijpers TW. Toll-like receptor-induced reactivity and strongly potentiated IL-8 production in granulocytes mobilized for transfusion purposes. Blood. 2010;115(22):4588–96.
Kim JE, Bauer MM, Mendoza KM, Reed KM, Coulombe RA Jr. Comparative genomics identifies new alpha class genes within the avian glutathione S-transferase gene cluster. Gene. 2010;452(2):45–53.
Gao GQ, Song LS, Tong B, Li GP. Expression levels of GSTA2 and APOD genes might be associated with carotenoid coloration in golden pheasant (Chrysolophus pictus) plumage. Dongwuxue Yan jiu = Zoological research. 2016;37(3):144–50.
Bonneaud C, Mazuc J, Gonzalez G, Haussy C, Chastel O, Faivre B, et al. Assessing the cost of mounting an immune response. Am Nat. 2003;161:367–79.
Ewing B, Green P. Base-calling of automated sequencer traces using phred. II. Error probabilities. Genome Res. 1998;8(3):186–94.
Ewing B, Hillier L, Wendl MC, Green P. Base-calling of automated sequencer traces using phred. I. Accuracy assessment. Genome Res. 1998;8(3):175–85.
Trapnell C, Roberts A, Goff L, Pertea G, Kim D, Kelley DR, et al. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and cufflinks. Nat Protoc. 2012;7(3):562–78.
Kersey PJ, Allen JE, Allot A, Barba M, Boddu S, Bolt BJ, et al. Ensembl genomes 2018: an integrated omics infrastructure for non-vertebrate species. Nucleic Acids Res. 2018;46(D1):D802–d8.
Consortium TGO. Expansion of the gene ontology knowledgebase and resources. Nucleic Acids Res. 2017;45(D1):D331–d8.
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, et al. Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000;25(1):25–9.
Kinsella RJ, Kahari A, Haider S, Zamora J, Proctor G, Spudich G, et al. Ensembl BioMarts: a hub for data retrieval across taxonomic space. Database : the journal of biological databases and curation. 2011;2011:bar030.
We thank Sabine Waigel and the University of Louisville Genomics Core for running samples and assistance with analyzing RNA-seq data. We thank Rachel Hoffman for husbandry and care of our zebra finch colony.
Part of this work was performed with assistance of the UofL Genomics Facility, which is supported by NIH P20GM103436 (KY IDeA Networks of Biomedical Research Excellence), NIH P30GM106396 (UofL J. G. Brown Cancer Center Phase III CoBRE), the J. G. Brown Foundation, and user fees. The contents of this work are solely the responsibility of the authors and do not represent the official views of the NIH or the National Institute for General Medical Sciences (NIGMS). Support for JHC and ECR was provided by P20 GM103436 (Nigel Cooper, PI). NTA was supported by NIH R15GM117534 and NSF IOS-1557882. None of the funding agencies listed above were involved in the design of the study or collection, analysis, and/or interpretation of data presented herein.
Ethics approval and consent to participate
All experimental and animal handling procedures were approved by the Institutional Animal Care and Use Committee at Western Kentucky University (#16–01), and followed the National Institutes of Health “Guide for the Use and Care of Laboratory Animals” and international ethical standards.
Consent for publication
The authors declare they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Upregulated DEGs in hypothalamus. Raw DEG data for upregulated genes in hypothalamus. (XLSX 344 kb)
Upregulated DEGs in spleen. Raw DEG data for upregulated genes in spleen. (XLSX 267 kb)
Upregulated DEGs in red blood cells. Raw DEG data for upregulated genes in red blood cells. (XLSX 42 kb)
Downregulated DEGs in hypothalamus. Raw DEG data for downregulated genes in hypothalamus. (XLSX 141 kb)
Downregulated DEGs in spleen. Raw DEG data for downregulated genes in spleen. (XLSX 253 kb)
Downregulated DEGs in red blood cells. Raw DEG data for downregulated genes in red blood cells. (XLSX 26 kb)
Gene ontology analyses for each tissue. Lists of 30 most significantly overrepresented GO terms for upregulated and downregulated genes of the hypothalamus, spleen, and RBCs. (DOCX 32 kb)