Transcriptional profiling of liver tissues in chicken embryo at day 16 and 20 using RNA sequencing reveals differential antioxidant enzyme activity

PLOS ONE, Feb 2018

Considering the high proportion of polyunsaturated fatty acids, the antioxidant defense of chick embryo tissues is vital during the oxidative stress experienced at hatching. In order to better understand the mechanisms of the defense system during chicken embryo development, we detected the activity of antioxidant enzymes during the incubation of chicken embryo. Results showed that the activity of superoxide dismutase (SOD) and (GSH-PX) in livers were higher than those in hearts. Based on these results, liver tissues were used as the follow-up study materials, which were obtained from chicken embryo at day 16 and day 20. Thus, we used RNA sequencing (RNA-Seq) analysis to identify the transcriptome from 6 liver tissues. In total, we obtained 45,552,777–45,462,856 uniquely mapped reads and 18,837 mRNA transcripts, across the 6 liver samples. Among these, 1,154 differentially expressed genes (p<0.05, foldchange≥1) were identified between the high and low groups, and 1,069 GO terms were significantly enriched (p<0.05). Of these, 10 GO terms were related to active oxygen defense and antioxidant enzyme activity. GO enrichment and KEGG pathway analysis indicated that GSTA2, GSTA4, MGST1, GPX3, and HAO2 participated in glutathione metabolism, and were considered as the most promising candidate genes affecting the antioxidant enzyme activity of chicken embryo at day 16 and day 20. Using RNA-Seq and differential gene expression, our study here investigated the complexity of the liver transcriptome in chick embryos and analyzed the key genes associated with the antioxidant enzyme.

A PDF file should load here. If you do not see its contents the file may be temporarily unavailable at the journal website or you do not have a PDF plug-in installed and enabled in your browser.

Alternatively, you can download the file locally and open with any standalone PDF reader:

http://journals.plos.org/plosone/article/file?id=10.1371/journal.pone.0192253&type=printable

Transcriptional profiling of liver tissues in chicken embryo at day 16 and 20 using RNA sequencing reveals differential antioxidant enzyme activity

February Transcriptional profiling of liver tissues in chicken embryo at day 16 and 20 using RNA sequencing reveals differential antioxidant enzyme activity Shaohua Yang 1 2 3 Lu Lu Wang 1 2 3 Zhaoyuan Shi 1 2 3 Xiaoqian Ou 1 2 3 Wei Wang 0 1 3 Xue Chen 1 2 3 Guoqing Liu 1 2 3 0 Agricultural Products Quality and Safety Supervision and Management Bureau , Xuancheng, Anhui , P. R. China 1 Fundamental Research Funds for Industry- University-Research in Hefei University of Technology (XC2016JZBZ03), the Fundamental Research Funds for the Central Universities 2 College of Food Science and Engineering, Hefei University of Technology , Hefei, Anhui , P. R. China 3 Editor: Salah A. Sheweita, Alexandria University , EGYPT Considering the high proportion of polyunsaturated fatty acids, the antioxidant defense of chick embryo tissues is vital during the oxidative stress experienced at hatching. In order to better understand the mechanisms of the defense system during chicken embryo development, we detected the activity of antioxidant enzymes during the incubation of chicken embryo. Results showed that the activity of superoxide dismutase (SOD) and (GSH-PX) in livers were higher than those in hearts. Based on these results, liver tissues were used as the follow-up study materials, which were obtained from chicken embryo at day 16 and day 20. Thus, we used RNA sequencing (RNA-Seq) analysis to identify the transcriptome from 6 liver tissues. In total, we obtained 45,552,777±45,462,856 uniquely mapped reads and 18,837 mRNA transcripts, across the 6 liver samples. Among these, 1,154 differentially expressed genes (p<0.05, foldchange 1) were identified between the high and low groups, and 1,069 GO terms were significantly enriched (p<0.05). Of these, 10 GO terms were related to active oxygen defense and antioxidant enzyme activity. GO enrichment and KEGG pathway analysis indicated that GSTA2, GSTA4, MGST1, GPX3, and HAO2 participated in glutathione metabolism, and were considered as the most promising candidate genes affecting the antioxidant enzyme activity of chicken embryo at day 16 and day 20. Using RNA-Seq and differential gene expression, our study here investigated the complexity of the liver transcriptome in chick embryos and analyzed the key genes associated with the antioxidant enzyme. Introduction Oxidative stress is always common in poultry production. During embryo growth, more oxygen is required to provide energy. However, elevated oxygen concentrations lead to high levels of reactive oxygen species (ROS) [ 1 ], which may which may cause protein and lipid oxidation [ 2 ]. Therefore, ROS has been proposed to cause many diseases and pathological changes (JZ2017HGTA0228), the innovation Project for College Students (2017CXCYS235), Wanjiang Institute of Poultry Technology (XC2015AKKG002), and the Science and Technology Bureau of AnHui (XC2015XKKJ02; XC2015AKKG003). during chick embryo development [ 3 ], especially in the cardiovascular system. Wells showed that excess ROS could have a teratogenic effect on developing embryos[ 4 ] as well as induce neural tube defects[ 5 ]. Meanwhile, it was found that ROS can also cause myocardial hypertrophy in the developing chick embryo [ 6 ]. The damaging effects of ROS can be exerted on the developing embryos in a directly or indirect style. Thus, antioxidant defences play a very important role in chick embryonic development. In fact, the integrated antioxidant systems within the egg and embryonic tissues are crucial for the protection of chick embryo in its development. Of these, the main antioxidant enzymes, superoxide dismutases (SOD), glutathione peroxidase (GSH-PX), glutathione S transferases (GSTs), Peroxidase (POD) and catalase(CAT), can clearly serve as a major defense line against ROS during the oxidative stress experienced at hatching in chick embryo[ 7 ]. During the incubation period, SOD converts highly reactive superoxide anions into H2O2 and O2, Catalase (CAT) catalyzes the dismutation of H2O2 to form the neutral products O2 and H2O, and GSH-PX catalyzes the reductive destruction of hydrogen and lipid hydroperoxides with glutathione as an electron donor [ 8 ]. Meanwhile, the end product of SOD is decomposed then to scavenge ROS [ 9 ]. Additionally, GSH-PX degrades hydrogen peroxide and other peroxides. The activities of these antioxidant enzymes are close relevant to higher levels of environmental oxygen. Using the chick embryo, van Golde et al[ 10 ] investigated the relations between hyperoxia and antioxidant enzyme activity, and found that SOD activity had a 2- to 10-foldincrease and Catalase and GPx enzyme activities remain almost the same in heart, liver, intestine and lungs during incubation at different time points. Starrs et al [ 11 ] examined the activities of catalase, SOD and GPx in the developing lungs of the chicken and showed that SOD activity decreased, whereas catalase and GPx activities were significantly increased in late incubation. Nevertheless, Dhage et al [ 12 ] revealed there was a significant increase of the SOD activity in the chick embryo from Day 4±11 (units/mg protein), respectively. Considering these inconsistent results, we still do not know the mechanism of genes that regulate antioxidant activity (SOD, CAT, and GSH-PX) during embryo development. Therefore, clarification of the differential gene expression underlying antioxidant enzyme activity in embryo development will have both biological and economic significance. As a research method, RNA-Seq provides a comprehensive and accurate tool for gene expression pattern analysis [ 13 ]. Among these, RNA-Seq results show gene structure, gene biological function, and new transcripts[14±15]. There have been several studies on various types of transcriptomes using RNA-Seq techniques, such as on fish [16±17], mice[ 18 ], cows[ 19 ], and pigs [ 20 ]. However, no studies on the chicken embryo transcriptome by RNA-Seq have been published. Thus, for a better understanding of the adaptive mechanism of antioxidant enzyme activity during the oxidative stress experienced at hatching in embryos, we examined liver transcriptome data from different incubation days in order to determine the key genes that were associated with antioxidant enzyme activity. Materials and methods Ethics statement All procedures for animal handling were reviewed and approved by the Institutional Animal Care and Use Committee (IACUC) of Hefei University of Technology (Permit Number: DK838). Materials Fertilized eggs were purchased from a hatchery (Changlv Native Products, Nanjing, China), and were incubated at 37.8ÊC and 60% relative humidity in an incubator (Photoshop Solar 2 / 14 Energy Co., Zibo, China). Fertilized eggs were studied at day 14, 15, 16, 17, 18, 19, and 20 of incubation. SOD, GSH-PX, Peroxidase (POD) and total protein quantitative assay kits were provided by Jiancheng Bioengineering Institute (Nanjing, Jiangsu, China). Preparation of tissues Fertilized eggs at day 14, 15, 16, 17, 18, 19, and 20 of incubation were obtained, peeled, and the heart and liver tissues were removed, respectively. The removed tissues were rinsed with cold saline (0.86%) to remove the blood, immediately frozen in liquid nitrogen, and placed in a centrifuge stored in an ultra-low temperature freezer. Determination of activity of antioxidant enzymes The collected tissue samples were weighed and made into a 20% tissue homogenate by adding the appropriate amount of saline (0.86%), according to the weight of the volume, then centrifuged at 4,000 rpm/min for 10 min at 4ÊC. The supernatant was obtained for further analysis. The activities of SOD, POD and GSH-PX were measured using commercial assay kits purchased from Nanjing Jiancheng Bioengineering Institute (Nanjing, Jiangsu, China). The antioxidant experiment was conducted in triplicate, and results were analyzed using SPSS 16.0 software (SPSS, Inc, Chicago, IL, USA). Comparison between groups was analyzed by Oneway analysis of variance (ANOVA) followed by Duncan's multiple range tests and the results were considered statistically significant at P < 0.05. RNA isolation and quality assessment Liver tissues at day16 and day 20 were chosen for transcriptome study. The total RNA was extracted from the embryo tissues using the Trizol method (Invitrogen, Carlsbad, CA) according to the manufacturer's instructions. RNA degradation and contamination was monitored on 1% agarose gels. Furthermore, the Nanophotometer (IMPLEN, CA, USA) was used to check RNA purity, and RNA concentration was measured using Qubit RNA Assay Kit in Qubit 2.0 Fluorometer (Life Technologies, CA, USA). The RNA integrity was assessed with the RNA Nano 6000 Assay Kit of the Bioanalyzer 2100 system (Agilent Technologies, CA, USA). Library construction and RNA sequencing A total amount of 3 μg RNA per sample was used as input material for the RNA sample preparations. Sequencing libraries were processed using NEBNext1 Ultra™ RNA Library Prep Kit for Illumina1 (NEB, USA), following the manufacturer's recommendations. Cluster generation of the index-coded samples was conducted on a cBot Cluster Generation System using Tru Seq PE Cluster Kit v3-c Bot-HS (Illumia) according to the manufacturer's instructions. After cluster generation, the library preparations were sequenced on an Illumina Hiseq platform, and 125-bp/150-bp paired-end reads were generated. Quality control Raw data (raw reads) of Fastq format were first processed through in-house perl scripts. In this step, clean data (clean reads) were acquired by removing reads containing adapter, reads containing ploy-N, and low quality reads from raw data. At the same time, Q20 (the proportion of bases with a phredbase quality score greater than 20; i.e., the proportion of read bases whose error rate is less than 1%), Q30 (the proportion of bases with a phredbase quality score greater than 30; i.e., the proportion of read bases whose error rate is less than 1%), and GC content of the clean data were calculated. All the downstream analyses were based on the clean data. 3 / 14 Reads mapping Reference genome and gene model annotation files were downloaded from a genome website directly (ftp://ftp.ensembl.org/pub/release-83/gtf/gallus_gallus/). Index of the reference genome was established using Bowtie v2.2.3. Because TopHat can generate a database of splice junctions based on the gene model annotation file (which possesses better mapping results than other non-splice mapping tools), paired-end clean reads were aligned to the reference genome using TopHat v2.0.12. Quantification of gene expression level The number of reads mapped to each gene was calculated by HTSeq v0.6.1. FPKM, expected number of Fragments Per Kilobase of transcript sequence per Millions base pairs sequenced, is currently the most commonly used method for estimating gene expression levels. FPKM of each gene was calculated based on the length of the gene and reads count mapped to this gene. Differential expression analysis Differential expression analysis was conducted using the DESeq R package (1.18.0). DESeq facilitates accurate comparisons between antioxidant enzyme activity of liver tissues by normalizing the number of reads, and provides statistical routines for determining differential expression in digital gene expression date using a model based on the negative binomial distribution. The resulting P-values were adjusted using the Benjamini and Hochberg's approach for controlling the false discovery rate. Genes with an adjusted P-value <0.05 and log2 (Fold change) 1 found by DESeq were assigned as differentially expressed. GO and KEGG enrichment analysis of differentially expressed genes Gene Ontology (GO) enrichment analysis of differentially expressed genes was performed by the GOseq R package, in which gene length bias was corrected. GO terms with P< 0.05 were considered significantly enriched by differentially expressed genes. KEGG provides comprehensive database resources for research of high-level functions and utilities of biological systems (http://www.genome.jp/kegg/). Statistical enrichment of differential expression genes in KEGG pathways was evaluated by KOBAS software. Real-time quantitative reverse-transcription-PCR (qRT-PCR) To verify the accuracy and repeatability of the transcription sequencing results, 10 differentially expressed genes were randomly selected to be detected using qRT-PCR. Designed by Primer3 (http://fokker.wi.mit.edu/primer3/input.htm), the primer sequences were shown in S1 File. The housekeeping gene GAPDH was used to correct the mRNA levels of differentially expressed genes. qRT-PCR was carried out in triplicate with the LightCycler 1 480 SYBR Green I Master Kit (Roche Applied Science, Penzberg, Germany), in a 15 μ L reaction on a LightCycler480 (Roche), using the following program: 95ÊC for 8 min, 45 cycles of 95ÊC for 10 s, 60ÊC for 15 s, 72ÊC for 10 s, and 72ÊC for 10 min. Results Antioxidant enzyme activities of embryos Antioxidant enzyme activities of SOD, GSH-Px and POD were detected in the incubation period. As shown in Fig 1A, the activity of SOD in the liver was higher than that in the heart tissue during incubation at different time points. In addition, there was a slight increase in 4 / 14 Fig 1. Effect of incubation day on SOD and GSH-PX activity in embryonic liver and heart. Note: The SOD activity (A) and GSH-PX (B) of heart and liver tissues were determined at days 14 to 20. The value of each fraction was the mean± standard deviation (n = 3), different letters (a, b) above columns indicate significant differences (p<0.05) in liver tissues, different letters (A, B) above columns indicate significant differences (p<0.05) in heart tissues. SOD activity in the liver, while it decreased in the heart with the increase in the incubation period. Fig 1A showed that SOD activity in the liver had a significant change from days 14 to 20 (25.86% increase, P<0.05) while Fig 1B showed that GSH-PX activity in the heart had an extremely significant change from days 14 to 20 (43.75% increase, P<0.01). Overall, the activity of GSH-PX in the liver tissue was significantly higher than in that in the heart. Meanwhile, the activity of POD in the heart and liver is shown in Fig 2. The results showed that POD activity in the liver increased initially and then decreased, and its activity had a significant change from days 14 to 20 (28.69% decrease, P<0.05). In the heart, the activity of POD increased steadily. According to the antioxidant activities, liver tissues at day 16 and 20 were used for further study. Sequencing and mapping of the liver transcriptome In total, we obtained 51,868,410±58,937,096 paired-end reads per sample. After removing the sequencing adaptors and poor quality reads, we acquired a total of 336,724,356 clean reads, and the total read length was 50.5 gigabases (GB) for six samples (Table 1). The data sets analyzed are available in the NCBI GenBank (https://www.ncbi.nlm.nih.gov/genbank) and the BioProject ID is PRJNA416967 (SRP123539). In this study, the Pearson correlation method was used for the calculation of the correlation coefficient (R2). The result showed that R2of biologically repeated samples was higher than 0.92, indicating that the similarity of the three biological replicates within each group was sufficiently high. Different gene expression between high and low groups for SOD activity Using the DEseq R package, the differential gene expression profile between liver tissues in chicken embryo at day 16 and day 20 was examined. In total, 18,837 genes were plotted. Meanwhile, 1,154 (571 down and 583 up) differentially expressed genes were identified (DEGs) at an FDR (false discovery rate) adjusted p-value <0.05, and absolute value of fold change 1 in 5 / 14 Fig 2. Effect of incubation day on POD activity in embryonic liver and heart. Note: The POD activity of heart and liver tissues was determined at days 14 to 20. The value of each fraction was the mean± standard deviation (n = 3), different letters (a, b) above columns indicate significant differences (p<0.05) in liver tissues, different letters (A, B) above columns indicate significant differences (p<0.05) in heart tissues. two samples. Volcano plots of the two comparison groups that are differentially expressed illustrate the distinct transcriptional profiles, displayed in Fig 3. The details of all DEGs are shown in S2 File. GO and pathway analysis of the DGEs To further evaluate the function of differentially expressed genes, GO classification analysis was used to annotate all genes identified from liver tissues, and formed into three categories: cellular components, biological processes, and molecular functions. These 1154 differential expressed genes were analyzed by GO enrichment, and 7731 GO terms were obtained. Of these, 1069 GO terms (13.83%) were significantly enriched (p<0.05) (S3 File). GO terms with P-value less than 0.05 were considered significantly enriched by DEGs. Among these GO terms, there are some GO terms related with cell cycle process, regulation of cell cycle and response to stress etc. In addition, there are 10 GO terms significantly related with antioxidant enzyme activity (P<0.01), such as glutathione transferase activity, a response to reactive oxygen species (Table 2). Meanwhile, we conducted metabolic pathway analysis using KOBAS software; the details of the significant pathway in the two-comparison group are shown in S4 File. The results of the KEGG analysis showed that several important pathways, such as ªglyoxylate and dicarboxylate metabolismº, ªcarbon metabolismº, and the ªp53 signaling pathwayº were significant. Real-time quantitative PCR Ten genes (PAPSS1, CCNB3, DYNLL1, GGT5, CLGN, ULK2, SPP1, VAV2, CEP170B, and SARS) were analyzed to confirm expression profiles and validate the transcriptome analysis 6 / 14 Note: Day 20 means the liver tissues in chicken embryo at day 20 and Day 20±1, Day 20±2, Day 20±3 mean the three biological replicate liver tissues in Day20 group, the rest of groups share the same name rules. Q20, the proportion of bases with a phred base quality score greater than 20; i.e., the proportion of read bases whose error rate is less than 1%.Q30, the proportion of bases with a phred base quality score greater than 30; i.e., the proportion of read bases whose error rate is less than 0.1%. results (Fig 4). The results showed that the gene expression levels were all consistent with mRNA-Seq results, which confirmed that the results obtained from the transcriptome sequencing platform were accurate. Candidate genes According to the GO database, within significantly enriched biological processes, 8 biological processes were related to reactive oxygen species and antioxidant enzyme activities. Among these processes, the most significant one was ªresponse to reactive oxygen speciesº. Integrated analysis of DEGs, GO and pathway results, and gene function allow us to suggest GSTA2, Fig 3. Volcano plot displaying DEGs within two different comparison groups. Note: the y-axis shows the mean expression value of log10(q-value), and the x-axis displays the log2fold change value. The blue dots represent the transcripts that did not reach statistical significance (q > 0.05); the red (up-regulated) and green dots (down-regulated) represent those whose expression levels were significantly different (q < 0.05); the blue dots represent the transcripts did not reach statistical significance (q > 0.05). 7 / 14 Note: GO ID indicates the label information in the Gene Ontology database, GO term refers to the description information of Gene Ontology function, P<0.05 means the function is an enriched item MGST1, GSTA4, GPX3 and HAO2 (Table 3) as the 5 promising candidate genes for affecting the activity of antioxidant enzymes during chicken embryo development. Discussion RNA-Seq has been applicated in many fields in chickens. Monson [ 21 ] established the liver library of domesticated turkey and wild turkey with AFB1 treatment by RNA-Seq and Fig 4. Validation of the gene expression profile by real-time PCR. Note: The x-axis represents the gene name, the yaxis represents the log2Ratio (Day20/Day16), different color columns represent data from RT-PCR or RNA-Seq. 8 / 14 obtained 89.2Gb of sequence. Sporer [ 22 ] utilized the microarray to detect the key genes in turkey skeletal muscle development and identified DEGs between two genetic lines of turkeys. Zhuo [ 23 ] revealed the differences of gene expression in the abdominal fat with high and low feed efficiency commercial broiler chickens. RNA-Seq was first used for analysis of early embryonic development in cattle and provided a method for further study mammalian embryonic development [ 24 ]. Based on those studies, we evaluated the whole genome transcriptome profile of chicken embryo liver tissues on different incubation days using RNA-Seq, with the aim of determining key genes that regulate the activity of antioxidant enzymes in this study. Compared with mammals, chicken embryonic development is carried out in a semi-closed system, and these natural antioxidants (SOD, GSH-PX and CAT) have been suggested to play a central role at hatching [ 25 ]. Therefore, the aim of this study is to detect the differentially expressed genes associated with the activities of antioxidant enzymes during chick embryonic development. In the present study, the activity of SOD and GSH-PX was significantly higher in the liver than in the heart. In the heart, GSH-PX activity value was half of that in the liver, which was consistent with results shown by Surai [ 26 ]. Additionally, CAT in the liver might be less important in ROS defense[27]. Thus, the results of antioxidant enzyme activity revealed that the liver was an important source of antioxidants. The different activities of antioxidant enzymes in the liver may be caused by different regulation of gene expression. According to the antioxidant enzyme activity detected in this study, the cDNA library was established by the liver. The transcript was deeply sequenced and the package of DESeq and Cuffdiff was used for analyses. We obtained 1,154 differentially expressed genes, while some of the genes have a known function, e.g., GSTA, GGT, GPX3, studies have also reported that these genes are associated with antioxidant enzyme activity [28±30]. Of these, DEGs, APOA4, LTC4S, GSTA2, MGST1, GSTA4, GPX3, and PEX11A were suggested to be the promising candidate genes for affecting the mechanisms of the defense system in response to ROS during chicken embryo development. Apolipoprotein (apo) A-IV (APOA4) was detected in the processes of ªresponse to reactive oxygen speciesº, ªresponse to oxidative stressº, ªresponse to lipid hydroperoxideº, ªcellular response to oxidative stressº, and ªresponse to hydroperoxideº. APOA4 is a 46 kDa glycoprotein and encodes a protein consisting of 396 amino acid residues. Kumar et al. [ 31 ] reported that APOA4 is an important mediator of lipid metabolism and has antioxidant activity. Studies have shown that the antioxidant effect of APOA4 is mediated by its ability to bind to the surface of the abundant lipoprotein lipase within atherosclerotic plaques [ 32 ]. Ostoa et al [ 33 ] 9 / 14 suggested that APOA4 and APOA1 accumulate in diseased artery walls to help protect against lipid peroxidation. As a member of the MAPEG (Membrane Associated Proteins in Eicosanoid and Glutathione metabolism) family of transmembrane proteins, leukotrieneC4 synthase (LTC4S) was down-regulated in this study. LTC4S encodes an enzyme that catalyzes the first step in biosynthesis of cysteinyl leukotrienes (LT), which possess important functions in inflammation [ 34 ]. Zhang et al. [ 35 ] showed that LTC4S and Orai3 can promote vascular smooth muscle cell migration and neointima formation by changing Akt signaling. Glutathione S-transferase alpha 2 (GSTA2) is a member of a family of glutathione S-transferases (GSTs), located in a cluster of similar genes and pseudogenes on chromosome 6. GSTA2 plays a role in detoxification by adding glutathione to target electrophilic compounds. Located in a cluster mapped to chromosome 6, GSTA2 is the most abundantly expressed glutathione S-transferase in the liver [36±37]. Additionally, GSTA2 has glutathione peroxidase activity, thereby protecting the cells from reactive oxygen species and the products of peroxidation. Tetlow [ 38 ] suggested that GSTA2 is a majorline of defense against oxidative stress. Glutathione S-transferase alpha4 (GSTA4) encodes a glutathione S-transferase belonging to the alpha class. The alpha class genes are highly related and encoded enzymes with glutathione peroxidase activity, that have a function in the detoxification of lipid peroxidation products [ 39 ]. Shearn [ 40 ] suggested that GSTA4 is a phase 2 detoxifying enzymes, and that its expression increases in response to oxidative stress. The microsomal glutathione S-transferase 1 (MGST1) gene encodes a protein that catalyzes the combination of glutathione to electrophiles and the reduction of lipid hydroperoxides. In addition, it scavenges reactive intermediates through its glutathione dependent transferase and peroxidase activities [ 41 ]. Glutathione peroxidase 3 (GPX3) belongs to the glutathione peroxidase family, which catalyzes the reduction of organic hydroperoxides and hydrogen peroxide by glutathione, and protects cells against oxidative damage [ 42 ]. Olson et al. [ 43 ] suggested that GPX3 was the only known selenocysteine-containing extracellular form of glutathione peroxidase. Barrett et al. [ 44 ] reported that knockdown of GPX3 in the human colon cancer cell line Caco2 caused an increase in ROS production. In addition, Xu et al.[ 45 ] suggested that vitamin E improves the antioxidant defense mechanisms, and enriches the GPX3 mRNA and protein expression levels, thereby enhancing the testicular antioxidant capacity. Peroxisomal biogenesis factor 11 alpha (PEX11A) is the richest ingredient of the peroxisomal membrane, and essential for proliferation of peroxisomes [ 46 ]. RodrÂõguez-Serrano et al [ 47 ] reported that in Arabidopsis, PEX11A lines exhibited higher levels of lipid peroxidation content and lower expression of genes involved in antioxidative defense and signaling. Weng et al. [ 48 ] found that the deficiency of PEX11A is related to peroxisome abundance reduction. KEGG metabolic pathway was used to analyze the function of differentially expressed genes, and 16 out of 139pathway terms was significantly enriched. Of these, ªGlutathione metabolismº and ªGlyoxylate and dicarboxylate metabolismº were pathway terms related to ROS defense. As a low molecular weight tripeptide, glutathione (GSH) is only present in a small quantity in the oxidation process[ 49 ]. Reduced GSH reduces the peroxide toH2O, and free radical reactions in vivo can be maintained. ªGlutathione metabolismº pathway involved 9 differentially expressed genes, which were all down-regulated. Presumably, the mechanism is that the up-regulation of GPX3 (1.11.1.9) promotes the conversion of GSH to GSSG, and glutathione reductase (GSR) (1.8.1.7) catalyzes the reduction of GSSG to GSH with NADPH as an H donor. Meanwhile, changes in the differentially expressed genes indicate that the glutathione metabolism in the liver is accelerated during the later stage of hatching. Additionally, hydroxyacid oxidase 2 (HAO2) was up-regulated in the conversion process of Glycolate to H2O2. The up-regulation of HAO2 indicated that it can catalyze more glycolate and generate 10 / 14 more H2O2. Mattu et al. [ 50 ] also revealed that the increased expression of HAO2 caused increased ROS production and lipid peroxidation. SOD, GSH-PX, and POD together play a role in the fight against ROS together in the embryo; the SOD activity is decreased with the generation of H2O2, when scavenging O2-. While adding CAT and GSH, H2O2can be decomposed. In this study, most genes in the ªglutathione metabolismº pathway were down-regulated, indicating more GSH involved in scavenging H2O2, thereby protecting SOD activity. PEX11a down-regulated revealed the decrease in POD activity, and more POD participating in ROS defense with GSH. In this comprehensive analysis of GO enrichment and the KEGG pathway, we discovered that some different gene expression not only exists in significantly enriched GO terms, but also participates significantly in the KEGG pathway; these genes were MGST1, GSTA2, GSTA4, GPX3, HAO2. This result revealed that these 5 genes may be key genes affecting ROS defense and antioxidant enzymes. These genes were associated with glutathione metabolism, and thus affected enzyme activity. This discovery also confirmed that GSH-PX played a primary role in antioxidant defense in the liver. Further research is required to understand the molecular mechanisms of these candidate genes on ROS defense and antioxidant enzyme in chickens. Conclusions Chicken embryo liver tissues at day 16 and day 20 were used as the follow-up study materials. In this study, we provided a comprehensive analysis of the complexity of the liver tissue transcriptome, and identified 1,154 differentially expressed genes between different incubation periods, with high and low antioxidant enzyme activity. GO enrichment and pathway analysis revealed 5 key genes affecting ROS defense and antioxidant enzymes, including MGST1, GSTA2, GSTA4, GPX3 and HAO2. Supporting information S1 File. PCR primers for qRT-PCR validation of 10 DEGs between the two different comparison groups. (DOCX) S2 File. List of DEGs in liver tissues within two different comparison groups. (XLSX) S3 File. GO enrichment list of DEGs between high and low antioxidant enzyme activity in biological processes. (XLSX) S4 File. List of significant KEGG pathway categories for DEGs. (XLSX) Author Contributions Conceptualization: Shaohua Yang, Guoqing Liu. Data curation: Shaohua Yang, Lu Lu Wang, Wei Wang. Formal analysis: Lu Lu Wang, Wei Wang, Xue Chen. Funding acquisition: Shaohua Yang, Wei Wang, Guoqing Liu. Investigation: Xue Chen, Guoqing Liu. Methodology: Xiaoqian Ou, Xue Chen. 11 / 14 Project administration: Shaohua Yang, Guoqing Liu. Resources: Lu Lu Wang, Zhaoyuan Shi. Software: Lu Lu Wang, Zhaoyuan Shi, Xiaoqian Ou. Supervision: Shaohua Yang, Guoqing Liu. Validation: Zhaoyuan Shi. Visualization: Zhaoyuan Shi, Xiaoqian Ou. Writing ± original draft: Lu Lu Wang. Writing ± review & editing: Shaohua Yang, Guoqing Liu. 12 / 14 27. Wilson JX, Lui EM, Del Maestro RF. Developmental profiles of antioxidant enzymes and trace metals in chick embryo. Mech Ageing Dev. 1992; 65: 51±64. PMID: 1405790 13 / 14 1. Turrens JF . Mitochondrial formation of reactive oxygen species . J Physiol . 2003 ; 552 : 335 ± 344 . https:// doi.org/10.1113/jphysiol. 2003 .049478 PMID: 14561818 2. He RR , Li Y , Li XD , Yi RN , Wang XY , Tsoi B , et al. A New Oxidative Stress Model , 2 , 2 -Azobis(2 -Amidinopropane) Dihydrochloride Induces Cardiovascular Damages in Chicken Embryo . PLOS one . 2013 ; 8: e57732 . https://doi.org/10.1371/journal.pone. 0057732 PMID: 23469224 3. Covarrubias L , HernaÂndez-GarcõÂa D , Schnabel D , Salas-Vidal E , Castro-ObregoÂn S . Function of reactive oxygen species during animal development: Passive or active? . Dev Biol . 2008 ; 320 : 1± 11 . https://doi.org/10.1016/j.ydbio. 2008 . 04 .041 PMID: 18555213 4. Wells PG , McCallum GP , Chen CS , Henderson JT , Lee CJ , Perstin J , et al. Oxidative stress in developmental origins of disease: teratogenesis, neuro developmental deficits, and cancer . Toxicol Sci . 2009 ; 108 : 4± 18 . https://doi.org/10.1093/toxsci/kfn263 PMID: 19126598 5. Han ZJ , Song G , Cui Y , Xia HF , Ma X . Oxidative stress is implicated in arsenic-induced neural tube defects in chick embryos . Int. J. Devl Neuroscience . 2011 ; 29 : 673 ± 680 . 6. Li Y , Wang XY , Zhang Z, Cheng X , Li XD , Chuai M , et al. Excess ROS induced by AAPH causes myocardial hypertrophy in the developing chick embryo . International Journal of Cardiology . 2014 ; 176 : 62 ± 73 . https://doi.org/10.1016/j.ijcard. 2014 . 06 .044 PMID: 25037699 7. Yigit AA , Panda AK , Cherian G. The avian embryo and its antioxidant defense system . Worlds Poult Sci J , 2014 ; 70 : 563 ± 574 . 8. Gill SS , Tuteja N. Reactive oxygen species and antioxidant machinery in abiotic stress tolerance in crop plants . Plant Physiol Biochem . 2010 ; 48 : 909 ± 930 . https://doi.org/10.1016/j.plaphy. 2010 . 08 .016 PMID: 20870416 9. Agarwal A , Gupta S , Sharma RK . Role of oxidative stress in female reproduction . Reprod Biol Endocrinol . 2005 ; 3 : 28 . https://doi.org/10.1186/ 1477 -7827-3-28 PMID: 16018814 10. van Golde JC , Borm PJ , Wolfs MC , Rhijnsburger EH , Blanco CE . Induction of antioxidant enzyme activity by hyperoxia (60% O2) in the developing chick embryo . Journal of Physiology . 1998 ; 509 : 289 ± 296 . https://doi.org/10.1111/j.1469- 7793 . 1998 .289bo. x PMID: 9547401 11. Starrs AP , Orgeig S , Daniels CB , Davies M , Lopatko OV . Antioxidant enzymes in the developing lungs of egg-laying and metamorphosing vertebrates . J Exp Biol . 2001 ; 204 : 3973 ± 3981 . PMID: 11807115 12. Dhage PA , Kamble LK , Bhargava SY . Localization and distribution of superoxide dismutase-1 in the neural tube morphogenesis of chick embryo . Int J Dev Neuro sci . 2017 ; 56 : 1 ± 9 . 13. Zhang WY , Dong R , Diao S , Du J , Fan ZP , Wang F . Differential long noncoding RNA/mRNA expression profiling and functional network analysis during osteogenic differentiation of human bone marrow mesenchymal stem cells . Stem Cell Reasearch and Therapy . 2017 ; 8 : 30 . 14. Yang D , Liu Q , Yang M , Wu H , Wang Q , Xiao J , et al. RNA-seq liver transcriptome analysis reveals an activated MHC-I pathway and an inhibited MHC-II pathway at the early stage of vaccine immunization in zebrafish . BMC Genomics . 2012 ; 13 : 319 . https://doi.org/10.1186/ 1471 -2164-13-319 PMID: 22805612 15. Wang ZH , Gerstein M , Snyder M. RNA-Seq: a revolutionary tool for transcriptomics . Nat Rev Genet . 2009 ; 10 : 57 ± 63 . https://doi.org/10.1038/nrg2484 PMID: 19015660 16. Cheng X , Zhao SG , Yue Y , Li HW , Wu JP . Comparative analysis of the liver tissue transcriptomes of Mongolian and Lanzhou fat-tailed sheep . GMR . 2016 ; 15 . 17. Magnanou E , Klopp C , Noirot C , Besseau L , Falcon J . Generation and characterization of the sea bass Dicentrarchus labrax brain and liver transcriptomes . Gene . 2014 ; 544 : 56 ± 66 . https://doi.org/10.1016/j. gene. 2014 . 04 .032 PMID: 24768179 18. Tran NT , Huang CH . Gene expression and gene ontology enrichment analysis for H3K4me3 and H3K4me1 in mouse liver and mouse embryonic stem cell using ChIP-Seq and RNA-Seq . Gene RegulSyst Bio . 2014 ; 8 : 33 ± 43 . 19. Cui X , Hou Y , Yang S , Xie Y , Zhang S , Sun D , et al. Transcriptional profiling of mammary gland in Holstein cows with extremely different milk protein and fat percentage using RNA sequencing . BMC Genomics . 2014 ; 15 : 226 . https://doi.org/10.1186/ 1471 -2164-15-226 PMID: 24655368 20. Sodhi SS , Song KD , Ghosh M , Sharma N , Lee SJ , Kim JH , et al. Comparative transcriptomic analysis by RNA-seq to discern differential expression of genes in liver and muscle tissues of adult Berkshire and Jeju Native Pig . Gene . 2014 ; 546 : 233 ± 242 . https://doi.org/10.1016/j.gene. 2014 . 06 .005 PMID: 24910116 21. Monson MS , Cardona CJ , Coulombe RA , Reed KM . Hepatic Transcriptome Responses of Domesticated and Wild Turkey Embryos to Aflatoxin B . Toxins , 2016 ; 8 : 16 . 22. Sporer KRB , Tempelman RJ , Ernst CW , Reed KM , Velleman SG , Strasburg GM . Transcriptional profiling identifies differentially expressed genes in developing turkey skeletal muscle . BMC Genomics . 2011 ; 12 : 143 . https://doi.org/10.1186/ 1471 -2164-12-143 PMID: 21385442 23. Zhuo Z , Lamont SJ , Lee WR , Abasht B . RNA-Seq Analysis of Abdominal Fat Reveals Differences between Modern Commercial Broiler Chickens with High and Low Feed Efficiencies . Plos One . 2015 ; 10 . 24. Huang W , Khatib H . Comparison of transcriptomic landscapes of bovine embryos using RNA-Seq . BMC Genomics . 2010 ; 11 : 711 . https://doi.org/10.1186/ 1471 -2164-11-711 PMID: 21167046 25. Niu ZY , Min YN , Liu FZ . Dietary vitamin E improves meat quality and antioxidant capacity in broilers by upregulating the expression of antioxidant enzyme genes . Journal of Applied Animal Research . 2017 ; https://doi.org/10.1080/09712119. 2017 .1309321 26. Surai PF . Vitamin E in avian reproduction . Poul Avian Biol Rev . 1999 ; 10 : 1± 60 . 28. Bawzduarte BG , Mendozacarrera F , Zamoraginez I , Torresrasgado E. Glutathione Peroxidase 3 serum levels and GPX3 gene polymorphismsin subjects with metabolic syndrome . Arch Med Res . 2014 ; 45 : 375 ± 382 . https://doi.org/10.1016/j.arcmed. 2014 . 05 .001 PMID: 24819036 29. Xu C , Zuo Z , Liu K , Jia H , Zhang Y , Luo H . Transcriptome analysis of the Tan sheep testes: differential expression of antioxidant enzyme-related genes and proteins in response to dietary vitamin E supplementation . Gene . 2016 ; 579 : 47 . https://doi.org/10.1016/j.gene. 2015 . 12 .045 PMID: 26723511 30. Pan L , Zhang XQ , Wang J.P , Ma X , Zhou M , Huang L , et al. Transcriptional profiles of drought-related genes in modulating metabolic processes and antioxidant defenses in Lolium multiflorum . FrontPlant Sci . 2016 ; 7 . 31. Kumar P , Nandi S , Tan TZ , Ler SG , Chia KS , Lim WY , et al. Highly sensitive and specific novel biomarkers for the diagnosis of transitional bladder carcinoma . Oncotarget . 2015 ; 6 : 13539 ± 13549 . https://doi. org/10.18632/oncotarget.3841 PMID: 25915536 32. Wong WM , Gerry AB , Putt W , Roberts JL , Weinberg RB , Humphries SE , et al. Common variants of apolipoprotein A-IV differ in their ability to inhibit low density lipoprotein oxidation . Atherosclerosis . 2007 192: 266 ± 274 . https://doi.org/10.1016/j.atherosclerosis. 2006 . 07 .017 PMID: 16945374 33. Ostos MA , Conconi M , Vergnes L , Baroukh N , Ribalta J , Girona J , et al. Antioxidative and antiatherosclerotic effects of human apolipoprotein A-IV in apolipoprotein E-deficient mice . Arterioscler Thromb Vasc Biol . 2001 21: 1023 ± 1028 . PMID: 11397714 34. Ahmad A , Ytterberg AJ , Thulasingam M , Tholander F , Bergman T , Zheng L. Structural insight into substrate selection and catalysis of lipid phosphate phosphatase PgpB in the cell membrane . J Bio Chem . 2016 ; 291:jbc . M116 .737874. 35. Zhang W , Zhang XX , GonzaÂlez-Cobos JC , Stolwijk JA , Matrougui K , Trebak M. LeukotrieneC4synthase , a critical enzyme in the activation of store-independent Orai1/Orai3channels, is required for neointimal hyperplasia . J Bio Chem . 2015 ; 290 : 5015 . 36. Coles BF , Kadlubar FF . Human alpha class glutathione S-transferases: genetic polymorphism, expression , and susceptibility to disease. Methods Enzymol . 2005 ; 401 : 9± 42 . https://doi.org/10.1016/S0076- 6879 ( 05 ) 01002 - 5 PMID: 16399377 37. Ahn J , Gammon MD , Santella RM , Gaudet MM , Britton JA , Teitelbaum SL , et al. Effects of glutathione S-transferase A1 (GSTA1) genotype and potential modifiers on breast cancer risk . Carcinogenesis . 2006 ; 27 : 1876 ± 1882 . https://doi.org/10.1093/carcin/bgl038 PMID: 16624829 38. Tetlow N , Coggan M , Casarotto MG , Board PG . Functional polymorphism of human glutathione transferase A3: effects on xenobiotic metabolism and steroid biosynthesis . Pharmacogenetics . 2004 ; 14 : 111 ± 116 . 39. Yang Y , Sharma R , Sharma A , Awasthi S , Awasthi YC . Lipid peroxidation and cell cycle signaling: 4- hydroxynonenal, a key molecule in stress mediated signaling . Acta BiochimicaPolonica . 2003 ; 50 : 319 ± 336 . 40. Shearn CT , Fritz KS , Shearn AH , Saba LM , Mercer KE , Engi B , et al. Deletion of GSTA4-4 results in increased mitochondrial post-translational modification of proteins by reactive aldehydes following chronic ethanol consumption in mice . Redox Biol . 2016 ; 7 : 68 ± 77 . https://doi.org/10.1016/j.redox. 2015 . 11 .013 PMID: 26654979 41. Schaffert CS . Role of MGST1 in reactive intermediate-induced injury . World JGastroenterol . 2011 ; 17 : 2552 ± 2557 . 42. Wang Y , Zhan X , Zhang X , Wu R , Yuan D . Comparison of different forms of dietary selenium supplementation on growth performance, meat quality, selenium deposition, and antioxidant property in broilers . Biol Trace Elem Res . 2011 ; 143 : 261 ± 273 . https://doi.org/10.1007/s12011-010 -8839-2 PMID: 20838922 43. Olson GE , Whitin JC , Hill KE , Winfrey VP , Motley AK , Austin LM , et al. Extracellular glutathione peroxidase (Gpx3) binds specifically to basement membranes of mouse renal cortex tubule cells . Aip Renal Physiol . 2010 ; 298 : 1244 ± 1253 . 44. Barrett CW , Ning W , Chen X , Smith JJ , Washington MK , Hill KE , et al. Tumor suppressor function of the plasma glutathione peroxidase Gpx3 in colitis-associated carcinoma . Cancer Res . 2013 ; 73 : 1245 ± 1255 . https://doi.org/10.1158/ 0008 - 5472 .CAN- 12 -3150 PMID: 23221387 45. Xu CC , Zuo ZY , Liu K , Jia HN , Zhang YW , Luo HL . Transcriptome analysis of the Tan sheep testes: Differential expression of antioxidant enzyme-related genes and proteins in response to dietary vitamin E supplementation . Gene . 2015 ; 579 : 47 ± 51 . https://doi.org/10.1016/j.gene. 2015 . 12 .045 PMID: 26723511 46. Opaliński è , Veenhuis M , Van IJ . Peroxisomes: membrane events accompanying peroxisome proliferation . Int J Biochem Cell Biol . 2011 ; 43 : 847 ± 851 . https://doi.org/10.1016/j.biocel. 2011 . 03 .006 PMID: 21419861 47. RodrÂõguez-Serrano M , Romero-Puertas MC , Sanz-FernaÂndez M , Hu J , Sandalio LM . Peroxisomes extend peroxules in a fast response to stress via a reactive oxygen species-mediated induction of the peroxin PEX11a . Plant Physiol . 2016 ; 171 : 1665 ± 1674 . https://doi.org/10.1104/pp. 16 .00648 PMID: 27208303 48. Weng H , Ji X , Naito Y , Endo K , Ma X , Takahashi R , et al. Pex11α deficiency impairs peroxisome elongation and division and contributes to nonalcoholic fatty liver in mice . Am J Physiol Endocrinol Metab . 2013 ; 304 : 187 ± 196 . 49. Reed MC , Thomas RL , Pavisic J , James SJ , Ulrich CM , Nijhout HF . A mathematical model of glutathione metabolism . TheorBiol Med Model. 2008 ; 5 : 8 . 50. Mattu S , Fornari F , Quagliata L , Perra A , Angioni MM , Petrelli A , et al. The metabolic gene HAO2 is down regulated in hepatocellular carcinoma and predicts metastasis and poor survival . J Hepatol . 2016 ; 64 : 891 ± 898 . https://doi.org/10.1016/j.jhep. 2015 . 11 .029 PMID: 26658681


This is a preview of a remote PDF: http://journals.plos.org/plosone/article/file?id=10.1371/journal.pone.0192253&type=printable

Shaohua Yang, Lu Lu Wang, Zhaoyuan Shi, Xiaoqian Ou, Wei Wang, Xue Chen, Guoqing Liu. Transcriptional profiling of liver tissues in chicken embryo at day 16 and 20 using RNA sequencing reveals differential antioxidant enzyme activity, PLOS ONE, 2018, DOI: 10.1371/journal.pone.0192253