Identification of Gene Expression Signatures in the Chicken Intestinal Intraepithelial Lymphocytes in Response to Herb Additive Supplementations
Article information
Abstract
Anethole and garlic have an immune modulatory effects on avian coccidiosis, and these effects are correlated with gene expression changes in intestinal epithelial lymphocytes (IELs). In this study, we integrated gene expression datasets from two independent experiments and investigated gene expression profile changes by anethole and garlic respectively, and identified gene expression signatures, which are common targets of these herbs as they might be used for the evaluation of the effect of plant herbs on immunity toward avian coccidiosis. We identified 4,382 and 371 genes, which were differentially expressed in IELs of chickens supplemented with garlic and anethole respectively. The gene ontology (GO) term of differentially expressed genes (DEGs) from garlic treatment resulted in the biological processes (BPs) related to proteolysis, e.g., “modification-dependent protein catabolic process”, “proteolysis involved in cellular protein catabolic process”, “cellular protein catabolic process”, “protein catabolic process”, and “ubiquitin-dependent protein catabolic process”. In GO analysis, one BP term, “Proteolysis”, was obtained. Among DEGs, 300 genes were differentially regulated in response to both garlic and anethole, and 234 and 59 genes were either up- or down-regulated in supplementation with both herbs. Pathway analysis resulted in enrichment of the pathways related to digestion such as “Starch and sucrose metabolism” and “Insulin signaling pathway”. Taken together, the results obtained in the present study could contribute to the effective development of evaluation system of plant herbs based on molecular signatures related with their immunological functions in chicken IELs.
INTRODUCTION
Coccidia of the genus Eimeria infects the intestinal epithelia of chickens, turkeys, and some mammalian hosts. And the infection causes substantial economic losses to the poultry industry such as decreased nutrient absorption, retarded growth of broilers, and decreased egg production of layers (Lillehoj and Li, 2004; Lillehoj et al., 2007). In-feed anticoccidial drugs and antibiotics growth promoters showed effectiveness in loss of weight gain by avian coccidiosis. However, there is an increasing governmental restriction in the use of antibiotics in livestock production worldwide (Lillehoj and Lillehoj, 2000; Lillehoj et al., 2011). Recently, the alternative to antibiotics has been explored to enhance innate immunity of poultry (Lee et al., 2010; Lillehoj and Lee, 2012). Especially, plant extracts are considered as feed additives to reduce the risk of coccidiosis owing to the beneficial effects of phytonutrients on infectious diseases and cancer (Duke et al., 2003; Veldhuizen et al., 2006; Burt et al., 2007; Chang et al., 2008; Tsubura et al., 2011).
Garlic (Allium sativum) is a medicinal plant that has been used for the treatment and prevention of some diseases and has been shown to have antimicrobial, antitumor, antioxidant, and immunostimulatory properties (Imai et al., 1994; Tsao and Yin, 2001). The pharmacologic effects of garlic on cardiovascular diseases are mediated through alteration of cholesterol and triglyceride levels in the bloodstream (Kyo et al., 2001). The medicinal effects of garlic are derived from its flavonoid and organo-sulfur components. Studies in mammalian systems have indicated that these active pharmacologic components act on multiple intracellular signaling pathways. In chickens, garlic compounds have also been shown to lower cholesterol and exert immunomodulating and antimicrobial activities (Chowdhury et al., 2002; Kim and Kim, 2007; Choi et al., 2010; Hanieh et al., 2010). However, compared with human studies, there is relatively limited information on the use of garlic in poultry veterinary medicine.
In addition, garlic has been suggested for the treatment of parasites and other intestinal diseases (Iciek et al., 2009). For example, crude extracts of garlic reduced or eliminated Hymenolepis, Aspiculuris, Histomonas, and Eimeria parasites in animal models of infection (Hafez and Hauck, 2006; Toulah and Al-Rawi, 2007; Ayaz et al., 2008; Abdel-Ghaffar et al., 2011; Dkhil et al., 2011). Of these parasites, coccidia protozoa comprise a subclass of single-celled eukaryotic microorganisms belonging to the phylum Apicomplexa and class Conoidasida. Coccidia of the genus Eimeria infect the intestinal epithelia of chickens, turkeys, and some mammalian hosts. Garlic has been shown to exert anticoccidial effects in mouse and rabbit infection models (Dkhil et al., 2011).
Anethole, 1-methoxy-4-(1-propenyl) benzene, is an aromatic, unsaturated ether that occurs as a major component of the essential oils of anise (Pimpinella anisum), star anise (Illicium verum), fennel (Foeniculum vulgare), and liquorice (Glycyrrhiza glabra), and is widely used as a flavoring agent (Windholz et al., 1983). In mammalian systems, anethole had demonstrable anti-carcinogenic, anti-oxidant, and anti-inflammatory activities (Al-Harbi et al., 1995; Chainy et al., 2000; Freire et al., 2005), as well as antimicrobial properties against bacteria, fungi, and nematodes (Kubo and Fujita, 2001; Lee et al., 2002). Broiler chickens fed with an anethole-supplemented diet showed resistance to coccidiosis (Kim et al., 2013a,b). These broad-spectrum pharmacologic activities suggest that anethole may also protect avian hosts against infectious diseases. However, only limited studies have been reported describing the use of anethole in veterinary medicine (Cardozo et al., 2005).
In this study, by integrating the public microarray data of the studies that investigated the effect of garlic and anethole on transcriptomes of the intestinal intraepithelial lymphocytes (IEL) of chicken, we aimed to identify the common signatures, which are targeted by both garlic and anethole to provide comprehensive strategy to evaluate plant herbs.
MATERIALS AND METHODS
Microarray data acquisition
To investigate common effects of two herbs, garlic and anethole, we used the Agilent Chicken Gene Expression Microarray dataset from two previous studies on their effects as the feed additives that were obtained from Gene Expression Omnibus. In one of the studies, a compound including garlic-derived secondary metabolites, propyl thiosulphinate and propyl thiosulphinate oxide, was added to the feed to enhance intestinal immunity of chicken (accession number GSE36302) (Kim et al., 2013a). The other study used anethole as the dietary supplement (accession number GSE41250) (Kim et al., 2013b). The design of microarray image analysis of both studies was a comparison of treatment (garlic and anethole) with non-supplemented control, identifying immunological effects of herbs.
The dataset used in this study was in reference to these previous studies, which contains details on feeding procedures of each herb (Kim et al., 2013a, b). The feeding procedures of each herb have been reported in Kim et al. studies (2013a, b). Briefly, chickens were fed from hatch with a diet supplemented with each herb or a standard diet alone (control). At 14 days post hatch, IEL were isolated from chickens given the herb-supplemented or non-supplemented diets. Total RNAs were isolated from the cells using Trizol (Invitrogen, Carlsbad, CA, USA) and amplified with cyanine 3 (Cy3)- or Cy5-labeled CTP. Then, the labeled RNAs were hybridized to a Chicken Gene Expession Microarray (Agilent Technologies, Santa Clara, CA, USA). In this study, microarray image analysis was performed to identify common and specific effects of two different herbs on transcriptional changes in chicken gut.
Identification and functional analysis of differentially expressed genes
R package ‘limma’ was used to normalize and qualify microarray images. Median signal intensities were corrected by adaptive background correction (Ritchie et al., 2007) and normalized by locally-weighted scatterplot smoothing method. The log2-transformed fold changes and standard errors were estimated by fitting a linear model and empirical Bayes statistics was applied for smoothing standard errors. Differentially expressed genes (DEG) were filtered by cutoff 0.1 of false discovery rate, with adjusted p-value of two-sample t-test. Annotation of DEGs and biological function analysis were performed using DAVID Bioinformatics Resources (http://david.abcc.ncifcrf.gov/).
RESULTS
Identification of differentially expressed genes
In the present study, microarray datasets from two independent studies were combined to analyze the discriminative effects of two different herbs, garlic, and anethole. Both of herbs were fed to the chickens as feed additives. In the results, garlic-fed chickens had greater number of modified genes than those fed with anethole. The numbers of DEGs were 4,682 (2,613 genes upregulated and 2,069 genes downregulated) and 371 (271 genes upregulated and 100 genes downregulated) in the supplementation of garlic-derived compound and anethole, respectively. The heatmap and hierarchical clustering for each experiment was displayed in supplemental Figure 1. The number of common DEGs found by both garlic and anethole additives was 300 (Figure 1). Detailed information for common DEGs was described in supplemental Table 1. The DEGs were used to identify their biological functions and the pathways in which the genes are involved.
Biological function analysis
In the gene ontology (GO) analysis, 38 and 1 terms of biological process (BPs) were significantly identified from the DEGs by the treatment of garlic and anethole, respectively (Table 1). The only BP that was classified from the DEGs by anethole treatment was the term of “Proteolysis”. Among the BPs from the DEGs by garlic treatment, there were several terms related with proteolysis such as “Modification-dependent protein catabolic process”, “Proteolysis involved in cellular protein catabolic process”, “Cellular protein catabolic process”, “Protein catabolic process”, and “Ubiquitin-dependent protein catabolic process”.
Figure 2 showed the correlated expression of the common DEGs changed by feeding of both herbs. Among the 300 common DEGs, 234 and 59 genes were up- and down-regulated in the treatment of both herbs, respectively. The remaining 7 genes did not show any correlated expression for both herbs. The terms of BPs, cellular components and molecular functions (MFs) and by GO analysis for the DEGs with same directional expression changes by garlic and anethole were listed in Table 2, 3, and 4, respectively. KEGG pathway terms enriched by common DEGs from the treatment of two different herbs were listed in Table 5. Both herbs showed common effects on 10 pathways, including adipocytokine signaling, vascular smooth muscle contraction, focal adhesion, starch and sucrose metabolism, drug metabolism, metabolism of xenobiotics by cytochrome P450, tight junction, insulin signaling pathway, glutathione metabolism and lysine degradation (Table 5).
DISCUSSIONS
For the comparative analysis of the effects of two different herbs, we collected microarray datasets from the independent trials for each herb, garlic and anethole, and identified GO terms and pathways enriched by the DEGs. In terms of number of DEGs, garlic showed more gene expression than anethole. In the present study, the most enriched category of BP affected by anethole was “Proteolysis” and those of BPs by garlic contained 8 processes related with “Catabolic processes” (Table 1). Proteolysis is the breakdown of proteins into smaller polypeptides or amino acids. The catabolic process refers to the chemical reactions and pathways resulting in the breakdown of substances, including the breakdown of carbon compounds with the liberation of energy for use by the cell or organism (De Bolster, 1997). These two GO terms are related to feed intake. The improvement of feed intake and daily live weight gain was observed in 1% anise supplement in Broiler chicken (Al-Kassie, 2008). This improvement may be due to the appetizing effect of active ingredient (such as anethole) in anise (Cabuk et al., 2003). In previous studies, garlic powder supplement and oil improved feed digestibility in cattle and sheep (Klevenhusen et al., 2011; Manasri et al., 2012). Therefore, the results suggest that anethole and garlic may function as important herbs for feed digestibility in chicken.
Interestingly, most of common DEGs by the treatment of anethole and garlic showed same directional changes in gene expression (Figure 2). Furthermore, in GO analysis in Table 2, 3, and 4, for the DEGs with same direction changes, the BPs identified from the down-regulated genes by anethole and garlic were “Proteolysis” and “Digestion”. The MFs from the down-regulated genes by both herbs were “Glutathione transferase activity”, “Peptidase activity, acting on L-amino acid peptides” and “Peptidase activity”. Together, the two herbs analyzed in this study could have important roles in the function of digestion in the chicken gut and it may be regulated by the down-regulated genes.
In the pathway analysis (Table 5), the DEGs, which were altered by both herbs constitute the pathways related to digestion such as “Starch and sucrose metabolism” and “Insulin signaling pathway”. “Glutathione metabolism” is critical for the development of effective strategies to improve health and many diseases since glutathione plays important roles in antioxidant defense, nutrient metabolism, and regulation of cellular events (including gene expression, DNA and protein synthesis, cell proliferation and apoptosis, signal transduction, cytokine production and immune response, and protein glutathionylation) (Wu et al., 2004). Glutathione deficiency contributes to oxidative stress, which plays a key role in aging and the pathogenesis of many diseases (including kwashiorkor, seizure, Alzheimer’s disease, Parkinson’s disease, liver disease, cystic fibrosis, sickle cell anemia, HIV/AIDS, cancer, heart attack, stroke, and diabetes) (Townsend et al., 2003). Therefore, the common digestive properties of anethole and garlic may be utilized as feed supplement for chicken to improve the animal’s health.
Supplementary Information
ACKNOWLEDGMENTS
This work was supported by a grant from the Next-Generation BioGreen 21 Program (No. PJ01104401), Rural Development Administration, Republic of Korea.
Notes
CONFLICT OF INTEREST
We certify that there is no conflict of interest with any financial organization regarding the material discussed in the manuscript.