The liver transcriptome of two full-sibling Songliao black pigs with extreme differences in backfat thickness

Kai Xing，F(xiàn)eng Zhu，Liwei Zhai，Huijie Liu，Zhijun Wang，Zhuocheng Hou*and Chuduan Wang*

The liver transcriptome of two full-sibling Songliao black pigs with extreme differences in backfat thickness

Kai Xing1，F(xiàn)eng Zhu1，Liwei Zhai1，Huijie Liu1，Zhijun Wang2，Zhuocheng Hou1*and Chuduan Wang1*

Background:Fatness traits in animals are important for their growth，meat quality，reproductive performance，and immunity.The liver is the principal organ of the regulation of lipid metabolism，and this study used massive parallelized high-throughput sequencing technologies to determine the porcine liver tissue transcriptome architecture of two full-sibling Songliao black pigs harboring extremely different phenotypes of backfat thickness.

Backfat thickness，Liver，Pig，RNA-Seq

Background

innate tendency to over consume food，means that pigs are a good animal model for studying obesity［4］.

The extent of fat deposition can be determined by triacylglycerol synthesis and storage，lipid mobilization，and fatty acid oxidation［5］.The liver is one of the most important organs to regulate appetite and body weight in pigs，as well as playing a key role in regulating several metabolic processes［6］.In pigs，de novo cholesterol synthesis and fatty acid oxidation mainly take place in the liver［7］.Lipid hydrolysis in adipose tissue results in free fatty acid，which is combined with plasma albumin then transported to the liver for use as an energy source through oxidation［8］.

RNA-Seq technology for transcriptome profiling has previously been used to explore the transcriptome of pig liver tissue.Such transcriptomes were recently compared between a full-sibling（full-sib）pair of F2 females from a White Duroc×Erhualian resource population with extreme phenotypes in growth and fat deposition［9］.Other comparisons include a backcross of two female groups（H and L）with extreme intramuscular fatty acid composition（25%Iberian×75%Landrace）［10］，and Duroc×F2（Leicoma×German Landrace）cross

Consumer choice of pig meat is an important factor that dictates the principles of swine breeding worldwide. One aspect of this is the deposition of fat in the muscle and backfat，which is associated with growth rate，meat quality，and reproductive performance［1］.Backfat thickness is highly correlated with body fat radio，carcass crosssectional fat area ratios，and intramuscular fat，so is a good indicator of fat deposition in pigs［2］.Therefore，the selection of backfat thickness using B-mode real-time ultrasound is a practical and economical method in pig breeding for increasing feeding efficiency，carcass value，and consumer acceptance of pork［3］.Because of the increasing problem of human obesity worldwide，it is beneficial to regulate fat deposition in pig breeding through molecular markers. Additionally，their similarity to humans in body size and other physiological/anatomical features，including theirpigs with divergent skatole levels in backfat［11］.Although several studies have previously attempted to identify the genes and pathways involved in fatty traits in the liver，to our knowledge，they either lacked sufficient RNA-Seq liver samples or did not take into account the effect of different genetic background noise when analyzing liver regulation.

The Songliao black pig is a Chinese domestic breed with powerful stress resistance and good reproductive and fat deposition capabilities.It is therefore a good model for studying fatty deposition.In the present study，we used RNA-Seq to obtain the liver transcriptomes of two full-sib Songliao black pigs with a high variation in backfat thickness.The main aim of this study was to elucidate the genes and pathways involved in lipid metabolism in liver tissue using RNA deep sequencing technology.

Methods

Experimental design，animals，and phenotypes

The Songliao black female pig population（average age，217（range，216-218）days；average live weight，100 kg（range，92.5-116.4 kg）was housed in consistent and standard environmental conditions with natural，uncontrolled room temperature and light.Animals were fed three times a day and had access to water ad libitum.Pedigree information was available for all animals.Live backfat thickness was measured on the last 3/4 rib using B-mode real-time ultrasound（HS1500，Honda，Japan）.We analyzed a total of 53 individuals with full/half-sibs for backfat thickness to identify pairs with two divergent phenotypes.To minimize the noise of different genetic back grounds，full-sibs were selected as a priority.

We set out to compare transcriptome changes between two groups with a high variation in backfat thickness:pigs with higher backfat thickness（BH）and those with lower backfat thickness（BL）which had a backfat thickness 2-3 times lower than that of BH pigs.The chosen animals also had to have a similar backfat thickness within the same group（BH/BL）after adjustment for live body weight. Based on our criteria，experimental samples were made up of two pairs of pigs with extreme backfat thickness differences，both of which were full-sibs.

The chosen pigs were slaughtered according to guidelines for the ethical use and treatment of animals in experiments in China.Liver tissue was separated and stored in liquid nitrogen until analyzed.Total liver RNA was extracted using the total RNA extraction Kit（Bioteke，China）according to the manufacturer’s recommendations.The quality of total RNA was assessed by the 2100 Bioanalyzer（Agilent，USA）.

mRNA library construction and sequencing

mRNA was isolated from total RNA samples using oligo（DT）magnetic beads（Invitrogen，USA）.Purified mRNA was first fragmented by the RNA fragmentation kit（Ambion，USA），then a one paired-end library was prepared for each sample according to the manufacturer’s instructions.mRNA libraries were individually sequenced for foursamples（two from each of the BH and BL group）using the Illumina High-seq 2000 sequencing system.The libraries were sequenced using a multiplexed paired-ends protocol with 180 bp of data collected per run.The average insert size for the paired-end libraries was 180 bp.A total of four paired-end mRNA libraries were constructed individually for four liver samples.

Mapping and counting reads

Quality control and reads statistics were determined using FASTQC（http://www.bioinformatics.babraham.ac.uk/projects/fastqc/）.All reads were trimmed 20-bp from the 5′end according to the reads quality distributions.After removal of the sequencing adapt and low-complexity reads，all RNA-Seq reads were mapped on the reference pig genome（Sscrofa10.2）using TopHat v2.0.1 software［12］with default parameters.The annotation database Ensembl Genes v67 was used as a reference.Additionally，the intersect from BED Tools was used to count the number of reads mapping to exons，introns，and intergenic positions in the genome［13］.The reads count was measured using easy RNASeq software［14］to quantify the raw reads mapped on each gene.

Differential expression and novel transcript analysis

The trimmed mean of M-values（TMM）was used to normalize gene expression levels［15］.After normalization，the NOISeq package implanted in the R computation environment was used to detect differentially expressed genes（DEGs）between two groups［16］.This method infers the noise distribution from the data and performs pair wise comparisons of the samples to identify DEGs.To measure expression level changes between two conditions，NOISeq takes into consideration two statistics:M（the log2-ratio of the two conditions）and D（the absolute value of the difference between conditions）.The probability thresholds were P≥0.8 and the TMM value in the lower expressed sample was≥1.The higher the probability，the greater the change in expression between the two groups.Using a probability threshold of 0.8 means that the gene is 4 times more likely to be differentially expressed than non-differentially expressed［17］.Novel transcripts，differentially expressed spliced transcripts，and primary transcripts were also detected using the Cufflinks suite of software for RNA-Seq［12］.

Table 1 The traits of backfat thickness and related fat deposition

Table 2 The number of reads obtained and percentages of mapped reads per sample

Functional enrichment analysis of differentially

expressed genes

Because the pig genome is poorly annotated，pig gene IDs were converted to human gene IDs using BioMart.DEG lists were submitted to the Database for Annotation，Visualization and Integrated Discovery（DAVID）bioinformatics resource for enrichment analysis of the significant overrepresentation of GO biological processes（GO-BP），molecular function（GO-MF），cellular component（GOCC），and KEGG-pathway category［18］.In all tests，P-values were calculated using Benjamini-corrected modified Fisher’s exact test and≤0.05 was taken as a threshold of significance. To further identify the DEG interaction network in the liver，the SearchTool for the Retrieval of Interacting Genes（STRING）was used，which is based on a known proteinprotein interaction database program.It generates a network of interactions from a variety of sources，including different interaction databases，text mining，genetic interactions，and shared pathway interactions［19］.

Quantitative PCR and data analysis

Total RNA was extracted from the liver and converted into cDNA using the Revert Aid?First Strand cDNA Synthesis Kit（Thermo Fisher Scientific Inc，USA）following the manufacturer’s protocol.cDNA samples were analyzed with real-time reverse transcriptase（RT）-PCR using the Light Cycler?480 Real-Time PCR System（Roche，USA）. RT-PCR reactions were performed in a final volume of 20 μl with the Roche SYBR Green PCR Kit（Roche）according to the manufacturer’s instructions.Pig GAPDH was used as an internal standard to correct the cDNA input.Triplicate RT-qPCRs were performed for each cDNA and the average Ct was used for further analysis.Relative quantification values were calculated using the 2-ΔΔCtmethod.

Data availability

Complete data sets have been submitted to NCBI Sequence Read Archive（SRA）under Accession no.SRP035376，Bioproject:PRJNA234465.

Results

Analysis of RNA deep sequencing data

The backfat thickness of the pig carcass and weight of the kidney were shown to differ greatly between groups（Table 1）.In general，individuals in the BH group had twice the backfat thickness compared with those in the BL group.Our experimental population is a conserved breed that has not undergone extensive selection nor hybridized with other breeds.Individuals therefore have very similar genetic backgrounds.

We obtained approximately 53 million paired-end clean reads of 90 bp for each sample，and high percentages of mapped reads ranging from 84.80 to 92.90%. The number of reads and percentages of mapped reads were similar between the two groups（Table 2）.Most mapped reads were located within an exon，with percent ages ranging from 72.85 to 76.21%.Other reads mapped within the untranslated region，introns，and intergenic regions.The percentages of reads in each region are shown in Table 2.

The total number of genes expressed in the liver in the four samples ranged from 16，815 to 17，025（Additional file 1: Table S1），with numbers of expressed genes being similar between the two groups.Correlations between biological replicate samples showed that the expressed genes were very highly reproducible，suggesting that a major fraction of the liver transcriptome is conserved between groups.

To confirm changes in transcript levels between BH and BL groups，six genes related to fatty acid synthesis or lipid metabolism were selected for RT-PCR analysis: ACACA，LDHA，ELOVL6，CYP1A2，PDK1，and SCD. We designed RT-PCR primers（Additional file 2:Table S2）for these genes，using GAPDH as a reference.When gene expression levels were compared，a strong correlation between RT-qPCR and RNA-Seq platforms was observed（0.72），confirming the high reproducibility of the data.For all six genes，the fold-change ratios between H and L groups were consistent with the RNASeq data（Figure 1）.

Differentially expressed genes between H and L groups

We quantified transcript expression levels in TMM to normalize gene expression data across different samples.Differential gene expression in liver tissue was calculated from the raw reads using the NOISeq package.We treated high/low backfat thickness samples as biological replicates because they showed similar phenotypes（Table 1）.This process identified 92 liver DEGs between pigs with extreme high and low backfat thickness levels（Additional file 3:Table S3）We also identified 587 and 690 genes that were only expressed in the BH or BL group，respectively（Figure 2）.Comparing shared DEGs between the two different pairs of pigs，we found that DEGs of biological replicates were more homogeneous with fewer false positives. Because of the sample limitations for each pair，we onlypresent functional analysis of DEGs obtained by treating the two pig pairs as two biological replicates.

Figure 1 Comparison of qPCR and RNA-Seq expression ratios（BH and BL groups）for selected genes.

Figure 2 Gene expression in BH and BL groups.A.Red points represent genes that are significantly differentially expressed.B.Venn diagram showing genes only expressed in the BH group（blue circle），only expressed in the BL group（yellow circle），and common to both groups（intersection）.

Table 3 Pathways enriched in up-regulated genes in the liver（BH vs.BL）

Novel transcripts

We used the reference annotation based transcript algorithm implemented in Cufflinks to identify novel transcripts that were not annotated in the current Ensembl pig gene annotation database.A total of 8，226 novel transcripts were detected for FPKM≥1 in the Songliao black pigs.Thirtyone spliced transcripts and 33 primary transcripts were significantly differentially expressed between the two groups（Additional file 4:Table S4）.

Functional enrichment analysis of the DEGs

Of the 92 DEGs，39 were up-regulated and 53 were down-regulated in the BH group compared with the BL group.Pig gene IDs were converted to human gene IDs，but five genes did not match with their human homologs（ENSSSCG00000008012，ENSSSCG00000016695，ENSSSCG00000030368，ENSSSCG00000010427，and ENSSSCG00000012881）.

To gain an insight into the liver tissue processes that differ between the BL and BH groups of pigs，differentially up-regulated and down-regulated genes underwent separate pathway analysis and gene ontology analysis using DAVID.Human homologs were recognized for 32 of the 39 up-regulated genes，and up-regulated genes were found to be involved in the metabolism of xenobiotics，drugs，retinol，and tryptophan（Table 3）.The expression of DEGs in the enriched pathways is shown in Table 4.Following GO analysis，the DEGs were shown to be related to biological processes such as amino acid biosynthesis andmetabolism，small molecule metabolism，and oxidation reduction.Forty-five of the 53 down-regulated genes had annotations in DAVID.No significant pathway was found to be associated with down-regulated genes.GO analysis showed that the biological processes enriched by DEGs are complex and relate to protein transport and enzyme activity（Figure 3）.

Table 4 Expression levels of genes in up-regulated pathways in the liver（BH vs.BL）

Protein-protein interaction analysis

To gain a better understanding of the biological relationships between genes，the integral DEG list was inputted into the STRING database.Most proteins encoded by DEGs were dissociative，and two protein-protein interaction networks were identified（Figure 4）.In one network，genes were associated with metabolism，detoxification，and superoxides，while genes in the other network were related to the heat stress response.

Discussion

In the present study，the percentages of mapped reads obtained per individual animal（84.80-92.90%）were higher than those in previous porcine liver transcriptome studies: 61.4-65.6% ［9］，71.42-77.75% ［10］，and 43-84% ［11］. Additionally，the percentage of reads within exons was higher than in previous studies，indicating that our results are more effective and credible.

Few studies have investigated how the liver transcriptome affects fat deposition.Our current gene expression analyses show that DEGs and pathways may play very important roles in this，and that genes involved in material metabolism and fatty acid transport are more active in the fatter groups of animals.These findings are consistent with those of previous studies.Previously，the expression levels of CA3，SERPINA6，GATM，GSTA2，and ALAS1 were also found to be significantly different between the two groups［9］.CA3 is associated with the internal fat rate and backfat thickness ofthe pig［20］，while SERPINA6 within the quantitative trait loci is associated with cortisol levels，fat，and muscle content，so can be considered a key regulator of obesity susceptibility［21］.These findings indicate that CA3 and SERPINA6 are strongly associated with fatty traits.

It is noteworthy that several DEGs identified in the present study（FABP-1，LCN2，PLIN2，CYP1A1，CYP1A2，CYP2A6，and CYP26A1）are involved in lipid metabolism. As a member of the family of fatty acid-binding proteins，F(xiàn)ABP-1 was expressed in all 12 tissues studied in a previous investigation，but transcript levels were more abundant in the liver and small intestine［22］.FABP-1 is involved in the transport of fatty acids to the site of β-oxidation，as well as the synthesis of triacylglycerol and phospholipids.It is also shown to interact with peroxisome proliferator-activated receptor（PPAR）α，which regulates fatty acid catabolism［23］. Compared with leaner pigs，F(xiàn)ABP mRNA and protein levels were expressed at higher levels in the muscle of fatter animals［24］.In the present study，F(xiàn)ABP-1 was up-regulated in the liver samples of BH pigs compared with BL.Consistent with our results，F(xiàn)ABP-1 was also up-regulated in obese ob/ob mice compared with controls［25］，suggesting that it’s increased expression may enhance fatty acid transport.

Figure 3 DEG GO analysis in liver tissue.Top three charts show gene ontology annotation processes（biological process（BP），cellular component（CC），and molecular functions（MF））of up-regulated genes.Lower three charts show processes（BP，CC，and MF）of down-regulated genes.

Figure 4 STRING analysis shows that DEGs are involved in known and predicted protein-protein interactions.STRING analysis of DEGs in the livers of BH and BL pigs.Network nodes represent genes shown in Additional file 2∶Table S2.Lines of different color represent seven types of evidence used in predicting associations.Red line∶fusion evidence；green line∶neighborhood evidence；blue line∶co-occurrence evidence；purple line∶experimental evidence；yellow line∶text mining evidence；light blue line∶database evidence；black line∶co-expression evidence.

The adipokine lipocalin 2（LCN2）is a member of the lipocalin family that transports small lipophilic ligands and is highly expressed by fat cells both in vitro and in vivo［26］.Levels of LCN2 mRNA are dramatically increased in the adipose tissue and liver of ob/ob mice，while reducing LCN2 expression leads to a decrease in PPARγ expression［27］.However，we found that LCN2 was down-regulated in the BH group compared with the BL group of pigs. Perilipin 2（PLIN2）is a cytosolic protein that promotes the formation and stabilization of intracellular lipid droplets，which are organelles involved in the storage of lipiddepots.PLIN2 polymorphisms have been associated with carcass traits including backfat thickness in pigs［28］，while increased PLIN2 mRNA expression was detected in the skeletal muscle of pigs with higher intermuscular fat［28，29］.Our results contradict this by revealing a lower PLIN2 mRNA expression level in the livers of pigs with higher backfat thickness.Finally，CYP1A1，CYP1A2，CYP2A6，and CYP26A1 belong to the highly diverse CYP450 super family，and also showed a differential expression pattern in the current study.CYPs have been shown to play critical roles in catalyzing metabolism reactions and in the oxidation of unsaturated fatty acids［30］. Because several differences appear to exist in observed expression levels between studies，further investigation into expression differences of these genes is warranted to elucidate their role in fat deposition.

The present study identified more pathways that were up-regulated than down-regulated in the fatter group of pigs.Most of these pathways are related to the metabolism of substances such as xenobiotics，drugs，retinol，and tryptophan.The retinoid metabolism pathway has previously been shown to be up-regulated in the liver. Retinoid is crucial for most forms of life，and many studies have identified an association between retinoid and lipid metabolism［31］.Retinoids regulate metabolism by activating specific nuclear receptors，including the retinoic acid receptor and the retinoid X receptor，an obligate heterodimeric partner for other nuclear receptors including PPARs.This helps to coordinate energy balance［32］and suggests that the liver has a larger burden in fatter compared with leaner pigs.

STRING analysis in the present study indicated that genes regulating metabolism influence those related to heat stress.Heat stress was previously shown to affect lipid peroxidation，causing serious damage to membrane lipids，lipoprotein，and other lipid-containing structures［33］.In pigs，body weight is positively correlated with heat stress［34］，but our current results show that mRNA levels of genes related to heat stress are down-regulated in the livers of fatter compared with leaner pigs.Future investigations into the interaction between body weight and heat stress may therefore identify novel methods to show how body weight affects healthy individuals.

Conclusions

This study undertook transcriptome analysis between two groups of Songliao black pigs with different backfat thicknesses.A total of 92 DEGs were identified between BH and BL groups.In concordance with the phenotypic differences，these genes belonged to pathways and gene networks related to lipid metabolism，regulation，and transport.Additionally，the identified DEGs related to heat stress could provide a new method of understanding and combating obesity.Our findings will be of use in understanding liver lipid regulation and in the design of new selection strategies to improve pig production.

Additional files

Additional file 1:Table S1.Gene expression count in four samples.

Additional file 2:Table S2.Primer sequences of six genes related to lipid metabolism for qRT-PCR.

Additional file 3:Table S3.Liver DEGs between BH and BL pigs with biological replicates.

Additional file 4:Table S4.Differentially expressed spliced transcripts and primary transcripts.

Abbreviations

BH∶Group with higher backfat thickness；BL∶Group with lower backfat thickness；DAVID∶Database for annotation，visualization and integrated discovery；DEG∶Differentially expressed gene；Full-sib∶Full-sibling；GO-BP∶Gene ontology biological processes；GO-CC∶GO cellular component；GO-MF∶GO molecular function；LCN2∶Lipocalin 2；PLIN2∶Perilipin 2；PPAR∶Proliferator-activated receptor；TMM∶Trimmed mean of M-values；RT∶Reverse transcriptase；STRING∶Search tool for the retrieval of interacting genes.

Competing interests

The authors declare that they have no competing interests.

Authors’contributions

KX carried out the experiment and drafted the manuscript.ZCH and CDW conceived the study，participated in its design and coordination，and helped draft the manuscript.FZ，HJL and ZJW helped sample，experiment and analysis of data.All authors read and approved the final manuscript.

Acknowledgment

This research was financially supported by the innovation research team for modern agricultural industry and technology in Beijing.

Author details

1National Engineering Laboratory for Animal Breeding and MOA Key Laboratory of Animal Genetics and Breeding，Department of Animal Genetics and Breeding，China Agricultural University，100193 Beijing，China.2Tianjin Ninghe Primary Pig Breeding Farm，Ninghe 301500，Tianjin，China.

Received：30 March 2014 Accepted：27 May 2014

Published：6 June 2014

1. Organ R∶Key in regulating appetite and body weight.；2012.

2. Suzuki K，Inomata K，Katoh K，Kadowaki H，Shibata T∶Genetic correlations among carcass cross-sectional fat area ratios，production traits，intramuscular fat，and serum leptin concentration in Duroc pigs.J Anim Sci 2009，87：2209-2215.

3. Fontanesi L，Schiavo G，Galimberti G，Calò DG，Scotti E，Martelli PL，Buttazzoni L，Casadio R，Russo V∶A genome wide association study for backfat thickness in Italian Large White pigs highlights new regions affecting fat deposition including neuronal genes.BMC Genomics 2012，13：583.

4. Houpt KA，Houpt TR，Pond WG∶The pig as a model for the study of obesity and of control of food intake：a review.Yale J Biol Med 1979，52：307.

5. Reiter SS，Halsey CH，Stronach BM，Bartosh JL，Owsley WF，Bergen WG∶Lipid metabolism related gene-expression profiling in liver，skeletal muscle and adipose tissue in crossbred Duroc and Pietrain pigs.Comp Biochem Physiol Part D Genomics Proteomics 2007，2：200-206.

6. Fam BC，Joannides CN，Andrikopoulos S∶The liver：Key in regulating appetite and body weight.Adipocyte 2012，1：259-264.

7. Mu?oz R，Estany J，Tor M，Doran O∶Hepatic lipogenic enzyme expression in pigs is affected by selection for decreased backfat thickness at constant intramuscular fat content.Meat Sci 2013，93：746-751.

8. Nguyen P，Leray V，Diez M，Serisier S，Bloc’h JL，Siliart B，Dumon H∶Liver lipid metabolism.J Anim Physiol Anim Nutr 2008，92：272-283.

9. Chen C，Ai H，Ren J，Li W，Li P，Qiao R，Ouyang J，Yang M，Ma J，Huang L∶A global view of porcine transcriptome in three tissues from a full-sib pair with extreme phenotypes in growth and fat deposition by paired-end RNA sequencing.BMC Genomics 2011，12：448.

10.Ramayo-Caldas Y，Mach N，Esteve-Codina A，Corominas J，Castelló A，Ballester M，Estellé J，Ibá?ez-Escriche N，F(xiàn)ernández AI，Pérez-Enciso M∶Liver transcriptome profile in pigs with extreme phenotypes of intramuscular fatty acid composition.BMC Genomics 2012，13：547.

11.Gunawan A，Sahadevan S，Cinar MU，Neuhoff C，Gro?e-Brinkhaus C，F(xiàn)rieden L，Tesfaye D，Tholen E，Looft C，Wondim DS∶Identification of the Novel Candidate Genes and Variants in Boar Liver Tissues with Divergent Skatole Levels Using RNA Deep Sequencing.PLoS One 2013，8：e72298.

12.Trapnell C，Roberts A，Goff L，Pertea G，Kim D，Kelley DR，Pimentel H，Salzberg SL，Rinn JL，Pachter L∶Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks.Nat Protoc 2012，7：562-578.

13.Quinlan AR，Hall IM∶BEDTools：a flexible suite of utilities for comparing genomic features.Bioinformatics 2010，26：841-842.

14.Delhomme N，Padioleau I，F(xiàn)urlong EE，Steinmetz LM∶easyRNASeq：a bioconductor package for processing RNA-Seq data.Bioinformatics 2012，28：2532-2533.

15.Dillies MA，Rau A，Aubert J，Hennequet-Antier C，Jeanmougin M，Servant N，Keime C，Marot G，Castel D，Estelle J，Guernec G，Jagla B，Jouneau L，Lalo? D，Le Gall C，Brigitte Scha?ffer B，Le Crom S，Guedj MM，Jaffrézic F∶A comprehensive evaluation of normalization methods for Illumina high-throughput RNA sequencing data analysis.Brief Bioinform 2013，14：671-683.

16.Tarazona S，García-Alcalde F，Dopazo J，F(xiàn)errer A，Conesa A∶Differential expression in RNA-seq：a matter of depth.Genome Res 2011，21：2213-2223.

17.Tarazona S，F(xiàn)urió-Tar? P，F(xiàn)errer A，Conesa A∶NOISeq:Differential Expression in RNA-seq；2013.

18.Huang DW，Sherman BT，Lempicki RA∶Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc 2009，4：44-57.

19.Franceschini A，Szklarczyk D，F(xiàn)rankild S，Kuhn M，Simonovic M，Roth A，Lin J，Minguez P，Bork P，von Mering C∶STRING v9.1：protein-protein interaction networks，with increased coverage and integration.Nucleic Acids Res 2013，41：D808-D815.

20.Wu J，Zhou D，Deng C，Xiong Y，Lei M，Li F，Jiang S，Zuo B，Zheng R∶Expression pattern and polymorphism of three microsatellite markers in the porcine CA3 gene.Genet Sel Evol 2008，40：227-239.

21.Ousova O，Guyonnet-Duperat V，Iannuccelli N，Bidanel JP，Milan D，Genet C，Llamas B，Yerle M，Gellin J，Chardon P，Emptoz-Bonneton A，Pugeat M，Mormède P，Moisan M∶Corticosteroid binding globulin：a new target for cortisol-driven obesity.Mol Endocrinol 2004，18：1687-1696.

22.Jiang YZ，Li XW，Yang GX∶Sequence characterization，tissue-specific expression and polymorphism of the porcine（Sus scrofa）liver-type fatty acid binding protein gene.Yi Chuan Xue Bao 2006，33：598-606.

23.Wolfrum C，Borrmann CM，Borchers T，Spener F∶Fatty acids and hypolipidemic drugs regulate peroxisome proliferator-activated receptors alpha-and gamma-mediated gene expression via liver fatty acid binding protein：a signaling path to the nucleus.Proc Natl Acad Sci U S A 2001，98：2323-2328.

24.Zhao S，Ren L，Chen L，Zhang X，Cheng M，Li W，Zhang Y，Gao S∶Differential expression of lipid metabolism related genes in porcine muscle tissue leading to different intramuscular fat deposition.Lipids 2009，44：1029-1037.

25.Chen R，Cao Y，Ma X，Wang Y，Zhang H，Zhang W，Zhou D∶Characterization of lipid metabolism-related genes in mouse fatty liver；2012.

26.Yan Q-W，Yang Q，Mody N，Graham TE，Hsu C-H，Xu Z，Houstis NE，Kahn BB，Rosen ED∶The adipokine lipocalin 2 is regulated by obesity and promotes insulin resistance.Diabetes 2007，56：2533-2540.

27.Zhang J，Wu Y，Zhang Y，LeRoith D，Bernlohr DA，Chen X∶The role of lipocalin 2 in the regulation of inflammation in adipocytes and macrophages.Mol Endocrinol 2008，22：1416-1426.

28.Davoli R，Gandolfi G，Braglia S，Comella M，Zambonelli P，Buttazzoni L，Russo V∶New SNP of the porcine Perilipin 2（PLIN2）gene，association with carcass traits and expression analysis in skeletal muscle.Mol Biol Rep 2011，38：1575-1583.

29.Gandolfi G，Mazzoni M，Zambonelli P，Lalatta-Costerbosa G，Tronca A，Russo V，Davoli R∶Perilipin 1 and perilipin 2 protein localization and gene expression study in skeletal muscles of European cross-breed pigs with different intramuscular fat contents.Meat Sci 2011，88：631-637.

30.Lewis DF∶57 varieties：the human cytochromes P450.Pharmacogenomics 2004，5：305-318.

31.Keller H，Dreyer C，Medin J，Mahfoudi A，Ozato K，Wahli W∶Fatty acids and retinoids control lipid metabolism through activation of peroxisome proliferator-activated receptor-retinoid X receptor heterodimers. Proc Natl Acad Sci U S A 1993，90：2160-2164.

32.Ziouzenkova O，Plutzky J∶Retinoid metabolism and nuclear receptor responses：New insights into coordinated regulation of the PPAR-RXR complex.FEBS Lett 2008，582：32-38.

33.Altan ?，Pabu?cuo?lu A，Altan A，Konyalio?lu S，Bayraktar H∶Effect of heat stress on oxidative stress，lipid peroxidation and some stress parameters in broilers.Br Poultry Sci 2003，44：545-550.

34.Ingram D，Legge K∶Effects of environmental temperature on food intake in growing pigs.Comp Biochem Physiol A Physiol 1974，48：573-581.

Cite this article as:Xing et al.∶The liver transcriptome of two full-sibling Songliao black pigs with extreme differences in backfat thickness.Journal of Animal Science and Biotechnology 2014 5∶32.

Submit your next manuscript to BioMed Central and take full advantage of:

? Convenient online submission

? Thorough peer review

? No space constraints or color figure charges

? Immediate publication on acceptance

? Inclusion in PubMed， CAS， Scopus and Google Scholar

? Research which is freely available for redistribution

Submit your manuscript at

www.biomedcentral.com/submit

10.1186/2049-1891-5-32

*

Correspondence∶zchou@cau.edu.cn；cdwang@cau.edu.cn

1National Engineering Laboratory for Animal Breeding and MOA Key Laboratory of Animal Genetics and Breeding，Department of Animal Genetics and Breeding，China Agricultural University，100193 Beijing，China Full list of author information is available at the end of the article

?2014 Xing et al.；licensee BioMed Central Ltd.This is an Open Access article distributed under the terms of the Creative

Commons Attribution License（http∶//creativecommons.org/licenses/by/4.0），which permits unrestricted use，distribution，and reproduction in any medium，provided the original work is properly credited.The Creative Commons Public Domain

Dedication waiver（http∶//creativecommons.org/publicdomain/zero/1.0/）applies to the data made available in this article，unless otherwise stated.

Results:The total number of reads produced for each sample was in the region of 53 million，and 8，226 novel transcripts were detected.Approximately 92 genes were differentially regulated in the liver tissue，while 31 spliced transcripts and 33 primary transcripts showed significantly differential expression between pigs with higher and lower backfat thickness.Genes that were differentially expressed were involved in the metabolism of various substances，small molecule biochemistry，and molecular transport.

Conclusions:Genes involved in the regulation of lipids could play an important role in lipid and fatty acid metabolism in the liver.These results could help us understand how liver metabolism affects the backfat thickness of pigs.