Polymorphisms in the Perilipin Gene May Affect Carcass Traits of Chinese Meat-type Chickens

Article information

Asian-Australas J Anim Sci. 2015;28(6):763-770

doi : https://doi.org/10.5713/ajas.14.0554

Lu Zhang ¹, Qing Zhu ¹, Yiping Liu ¹, Elizabeth R. Gilbert ², Diyan Li ¹, Huadong Yin ¹, Yan Wang ¹, Zhiqin Yang ¹, Zhen Wang ¹, Yuncong Yuan ¹, Xiaoling Zhao

Department of Animal Science, Sichuan Agricultural University, Ya’an Sichuan 625014, China.

¹Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, Sichuan Agricultural University, Chengdu, Sichuan 611130, China.

²Department of Animal and Poultry Sciences, Virginia Polytechnic Institute and State University, Blacksburg, VA24061, USA.

^*Corresponding Author: Xiaoling Zhao. Tel: +86-835-2886080, Fax: +86-835-2886080, E-mail: zhaoxiaoling04@163.com

Received 2014 July 23; Revised 2014 August 20; Accepted 2015 January 12.

Abstract

Improved meat quality and greater muscle yield are highly sought after in high-quality chicken breeding programs. Past studies indicated that polymorphisms of the Perilipin gene (PLIN1) are highly associated with adiposity in mammals and are potential molecular markers for improving meat quality and carcass traits in chickens. In the present study, we screened single nucleotide polymorphisms (SNPs) in all exons of the PLIN1 gene with a direct sequencing method in six populations with different genetic backgrounds (total 240 individuals). We evaluated the association between the polymorphisms and carcass and meat quality traits. We identified three SNPs, located on the 5′ flanking region and exon 1 of PLIN1 on chromosome 10 (rs315831750, rs313726543, and rs80724063, respectively). Eight main haplotypes were constructed based on these SNPs. We calculated the allelic and genotypic frequencies, and genetic diversity parameters of the three SNPs. The polymorphism information content (PIC) ranged from 0.2768 to 0.3750, which reflected an intermediate genetic diversity for all chickens. The CC, CT, and TT genotypes influenced the percentage of breast muscle (PBM), percentage of leg muscle (PLM) and percentage of abdominal fat at rs315831750 (p<0.05). Diplotypes (haplotype pairs) affected the percentage of eviscerated weight (PEW) and PBM (p<0.05). Compared with chickens carrying other diplotypes, H3H7 had the greatest PEW and H2H2 had the greatest PBM, and those with diplotype H7H7 had the smallest PEW and PBM. We conclude that PLIN1 gene polymorphisms may affect broiler carcass and breast muscle yields, and diplotypes H3H7 and H2H2 could be positive molecular markers to enhance PEW and PBM in chickens.

Keywords: Chickens; Perilipin Gene; Polymorphism; Carcass Traits; Association Analysis

INTRODUCTION

Consumption of poultry meat has grown rapidly over the last 50 years, as poultry meat is a high-quality dietary protein source for humans. The worldwide demands for meat will continue; meanwhile, the cost of poultry production, especially the feed and labor, continues to also rise. These factors have urged the breeders, nutritionists, and farmers to develop strategies to improve broiler feed efficiency and muscle yield (Sun and Wang, 2006; Petracci and Cavani, 2012). Molecular marker assisted selection is an effective way to improve these traits within a short breeding period (Lande and Thompson, 1990).

Studies indicate that the Perilipin (PLIN1) gene affects body weight (BW) and fat deposition in animals (Londos et al., 2005; Bickel et al., 2009). Perilipin, the lipid droplet coat protein from the lipid droplets related proteins family, exists in a variety of eukaryotic species ranging from amoebas to vertebrates (Miura et al., 2002). Perilipins share a highly conserved sequence and avidity for the surfaces of intracellular, neutral lipid storage droplets (Kimmel et al., 2010). Genetic variation in PLIN1 was associated with carcass traits and adiposity in humans (Ruiz et al., 2011), pigs (Gandolfi et al., 2011), cattle (Fan et al., 2010), sheep (Gao et al., 2012), ducks (Zhang et al., 2013), and chickens (Zhou et al., 2014). In humans, the rs6496589 G allele frequencies of PLIN1 were significantly associated with a lower risk of central obesity (Song et al., 2014). The various PLIN1 genotypes were associated with significant differences in carcass weight (CW), carcass net weight, and percentage of abdominal fat (PAF) in ducks (Zhang et al., 2013).

Research on PLIN1 gene polymorphisms in poultry has been reported, and developmental variants of PLIN1 expression have a relationship with fat and lipid parameters (Zhao et al., 2009). Certainly not all genetic variation is correlated with carcass and adipose traits; at least one study indicated that there was no significant effect on carcass and adipose traits when the PLIN1 gene harbored a G-T mutation at 2,224 bp (Lei et al., 2011). A previous study showed that the diplotypes of three single nucleotide polymorphisms (SNPs) in PLIN1, g.2272C>T, g.2319C>T and g.2467G>A, were significantly associated with BW, CW and leg muscle weight (LMW), and significantly associated with abdominal fat weight (AW) and PAF in Luqin chickens (Zhou et al., 2014).

Previous studies utilized polymerase chain reaction (PCR)-restriction fragment length polymorphism (RFLP), which has a lower resolution ratio and accuracy than direct sequencing. Thus, we used direct sequencing to analyze polymorphisms in the PLIN1 gene, which is located on chromosome 10 and contains 8 exons and 7 introns (http://asia.ensembl.org/Gallus_gallus/Transcript/Exons?db=core;g=ENSGALG00000023395;r=10:12315435-12319628;t=ENSGALT00000038622). We screened the genetic polymorphisms in all exons and partial introns of the chicken PLIN1 gene in six broiler populations to discern potential associations between the nucleotide polymorphisms and carcass traits, as well as intramuscular fat content (IMF).

MATERIALS AND METHODS

Chicken populations

In total, 240 chickens (male:female = 1:1) from six populations were used in the present study. The native Chinese breeds, Mountainous Black-bone (MB; n = 40) and Erlang Mountainous chickens (EM; n = 40) were indigenous breeds from the Sichuan province, China, with spotty feathers and black skin (Wang et al., 2012; Wang et al., 2014). The Hongbiao chickens (HB; n = 40) were indigenous breeds from Taiwan with red feathers and yellow skin (Gu and Dai, 2011). The commercial breeds, Daheng (DH; n = 40), BH1 (n = 40) and BH2 chickens (n = 40) were characterized by spotty feathers and yellow skin. The Daheng and BH1 and BH2 chickens were obtained from Sichuan Daheng Poultry Breeding Company and Banghe Agricultural Science and Technology Company, China (Zhou et al., 2009).

Ten thousand fertilized eggs from each population were hatched on the same day, and six thousand chicks were housed in batteries (battery size: 200 cm×60 cm×30 cm). There were 60 birds per pen with 100 pens of each population. At day 49, all individuals were transferred to single cages (size: 28 cm×55 cm×36 cm). Throughout the trial, birds had free access to feed and water. Diets fed were 19% crude protein (CP) and 2,897 kcal of ME/kg from 1 to 28 d, 17% CP and 2,998 kcal of ME/kg from 28 to 43 d, and 15% CP and 3,819 kcal of ME/kg from 43 to 91 d, which satisfied the dietary requirement of yellow-feather chickens (Zhao et al., 2012).

Data collection

We randomly sampled 240 individuals at 90 d of age for performance data collection and genotyping. Prior to slaughter, whole blood was collected from each chicken and stored at −20°C. Carcass traits measured included live BW after 12 h without access to feed, CW, semi-eviscerated weight, eviscerated weight (EW), breast muscle weight, LMW, and AW, as described by Zhao (Zhao et al., 2012). Intramuscular fat was extracted from muscle fibers by using the Soxhlet extraction technique (Rosa et al., 1983) and was defined as the content of IMF per 100 g of dried breast muscle. The dried breast muscle was prepared by baking the fresh muscle at 105°C in an oven for 120 min and then cooling to room temperature in desiccators.

Procedure for polymerase chain reaction and genotyping

We extracted genomic DNA from whole blood using the standard phenol/chloroform method (Wang et al., 2007). The Gallus gallus PLIN1 gene sequence (accession number NC 006097.3) was the basis for designing primer pairs via the OLIGO 6.0 program (Molecular Biology Insights, Inc., Cascade, CO, USA). Primers covered all 8 exons of the PLIN1 gene (Table 1). The PCR was performed in a final volume of 50 μL containing 200 ng genomic DNA, 1 μL each of forward and reverse primers (10 pmol/μL), and 25 μL 2×Taq PCR Master Mix (Beijing Tianwei Biology Technique Corporation, Beijing, China). The reaction mixture was predenatured for 4 min at 94°C, followed by 35 cycles at 94°C for 30 s, 57°C for 30 s, 72°C for 1 min, and ended with a full extension cycle at 72°C for 8 min. The PCR products were electrophoresed on a 1% agarose gel and were visualized on a Gel DocTMEQ170-8060 after staining with ethidium bromide. Samples showing different banding patterns on gels were reamplified, purified, and sequenced by a commercial sequencing company (Beijing Liuhe Huada Gene Technology Company, Beijing, China).

Table 1

Primers for screening SNPs in the chicken PLIN1 gene

Haplotype construction

The PHASE 2.0.2 was employed to build the most likely pair of haplotypes (diplotypes), using all three SNPs. This software implements a Bayesian statistical method for reconstructing haplotypes from population genotype data (Stephens et al., 2001). This method is applicable to genotypic data at linked loci from a population sample and may reduce error for haplotype construction by >5.0%. The parameters for PHASE input contained the number of individuals and SNPs, position of SNPs in the nucleotide sequence, symbol “sss” and the genotypes of all SNPs (Stephens et al., 2001). Diplotypes with a frequency less than 5.0% were excluded in the association analysis of diplotypes and carcass traits.

Linkage disequilibrium and statistical analysis

The pattern of pairwise linkage disequilibrium (LD) coefficient (D′) and correlation coefficient (r²) between the SNPs was calculated using the program Haploview 4.1 (Barrett, 2009). Allelic frequencies were determined by direct counting from observed genotypes. The x²-test was used to determine Hardy-Weinberg equilibrium (HWE) of the mutation. Population genetic indices, including gene homozygosity (Ho), gene heterozygosity (He), effective allele numbers (Ne), and polymorphism information content (PIC) were calculated according to Nei and Roychoudhury (Nei and Roychoudhury, 1974).

Data were analyzed by the GLM procedure of SAS 8.0 (SAS Inst. Inc., Cary, NC, USA). The genetic effects of diplotypes on carcass traits and IMF were analyzed by the following model:

Yijk=μ+Pi+Sj+Gk+(PG)ik+(PS)ij+(GS)ij+(PSG)ijk+eijk

Where Y_ijk = the dependent variable; μ = the population mean; P_i = fixed effect of population; S_j = fixed effect of sex; G_k = fixed effect of genotype/diplotype. (PG)_ik = interaction effect between population and genotype/diplotype; (PS)_ij = interaction effect between population and sex; (SG)_jk = interaction effect between genotype/diplotype and sex; (PSG)_ijk = interaction effect between population and genotype/diplotype and sex; e_ijk = random error. Statistical significance threshold was determined as p<0.05, unless otherwise specified. The false discovery rate (FDR) was calculated for corrections of multiple testing for each p-value (Benjamini and Hochberg, 1995). A significant association was defined as a value of 0.05 for the adjusted p-values.

RESULTS

Population effect on carcass and intramuscular fat traits

The population effects on chicken carcass traits and IMF are summarized in Table 2. There was a population effect on carcass traits and IMF (p<0.001). The IMF, percentage of leg muscle (PLM), percentage of breast muscle (PBM), and BW of HB chickens were greater than other chickens (p<0.001), whereas the BW, percentage of carcass weight (PCW), PSEW, PAF, and IMF were lower in MB than other chickens (p<0.001).

Table 2

Population effects on carcass traits and IMF in chickens

Variations in the chicken PLIN1 gene

A direct sequencing method was developed for screening the SNPs in the chicken PLIN1 gene. Of the 10 fragments that covered the 5′ flanking regions and entire exons of the PLIN1 gene, only one was polymorphic and three SNPs, rs315831750, rs313726543, and rs80724063 were identified (Figure 1). The rs315831750 is located in the 5′ flanking regions and the rs313726543 and rs80724063 are located in exon 1 of PLIN1.

Figure 1

Sequencing analysis of the Perilipin (PLIN1) gene. Arrows indicate mutation sites.

Genotype and allele frequencies and genetic diversity parameters for single nucleotide polymorphisms in the PLIN1 gene

The allelic and genotypic frequencies, and genetic diversity parameters (Ho, He, Ne, and PIC) of the three SNPs are shown in Table 3. In locus rs315831750, the allele T predominated due to its higher frequency than allele C in all populations. In locus rs313726543, the allele A predominated due to its higher frequency than allele G in most of the populations. In locus rs80724063, the frequency of allele C was higher than allele T in all populations, and the frequency of genotype TT was the smallest in all six populations and not detected in DH and MB chickens.

Table 3

Genotype and allele frequencies and genetic diversity parameters for the PLIN1 gene in six breeds

The x²-test performed to examine the HWE, revealed that the three loci in DH, MB, EM, BH1, and BH2 chickens, but not HB, were in HWE (p>0.05), and the PIC ranged from 0.2471 to 0.3736 (Table 3).

Effect of single nucleotide polymorphism genotype on carcass traits

Effects of single SNP genotype on carcass traits are summarized in Table 4. The single SNP effect was significant for PBM, PLM, and PAF. Chickens with the homozygous genotype CC had the greatest PBM, the heterozygous genotype CT the greatest PLM and smallest PAF, and the homozygous genotype TT the smallest PBM and PLM, but the greatest PAF at rs315831750. The single SNP genotype had no significant effect on carcass traits at the other two SNPs.

Table 4

Effect of single SNP genotype on carcass traits

Haplotype construction and association with carcass and intramuscular fat traits

Eight haplotypes were reconstructed based on the three SNPs. Among them, four haplotypes H1 (CGC), H2 (CGT), H6 (TGT), and H7 (TAC) were prevalent and accounted for 94.44% of the estimates (Table 5). Two main haplotypes, H2 (CGT) and H7 (TAC), accounted for 76.35% of the estimates. The haplotype H7 (TAC) had a higher frequency than other haplotypes, accounting for 55.79% of the estimates. Four haplotypes, H3 (CAC), H4 (CAT), H5 (TGC), and H8 (TAT), which had a frequency lower than 5%, were not included in Table 6. Diplotypes affected percentage of eviscerated weight (PEW) and PBM (p<0.05). The HB chickens were not analyzed in diplotypes effects on carcass traits and IMF because HB chickens were not in HWE. Compared with chickens carrying other diplotypes, H3H7 chickens had greater PEW and H2H2 greater PLM than other diplotypes, and those with diplotype H7H7 had less PEW and PBM than other diplotypes (Table 6).

Table 5

Haplotypes constructed with single nucleotide polymorphisms (SNP) and their frequencies

Table 6

Diplotype effects on PEW and PBM in chickens

The linkage disequilibrium analysis between the three single nucleotide polymorphisms

The pattern of pairwise LD structure as measured by D′ and r² between the SNPs was calculated using the program Haploview 4.1, and the haploblock structure presented in Figure 2. We performed LD and haplotype analysis for the three SNPs in all chickens, except for HB chickens. According to the data, the D′ value between SNP1 and SNP2, SNP2 and SNP3, and SNP1 and SNP3 were 0.81, 0.92, and 0.66, respectively. Moreover, r² for them were 0.601, 0.539, and 0.309, respectively.

Figure 2

Haploblock structures of the 5′flanking and exon 1 region as revealed by Haploview program.

DISCUSSION

There are many reports of polymorphisms in the PLIN1 gene in mammals, and dozens of SNP loci for the PLIN1 gene were associated with lipid metabolism and carcass traits (Marcinkiewicz et al., 2006; Fan et al., 2010). In humans, elevation of PLIN1 expression was correlated with increased adiposity (Kern et al., 2004), and PLIN1 polymorphisms affect diet-induced changes in body fat and energy metabolism (Ruiz et al., 2011). Polymorphisms in PLIN1 gene exon 2 and intron 2 have an influence on carcass and adipose traits, and PLIN1 was identified as a candidate gene affecting carcass and meat quality traits in ducks (Fan et al., 2011).

Research on polymorphisms of the PLIN1 gene in chickens has been reported, but the PCR-RFLP used in variation scaning had lower resolution ratio, compared with the direct sequencing used in the present study (Lei et al., 2011; Zhou et al., 2014). Moreover, we screened all exons of this gene for genetic variations in quite a few populations. Variations observed in PLIN1 mRNA sequences were rs315831750, rs313726543, and rs80724063 (Rubin et al., 2010). The rs313726543 and rs80724063 are synonymous mutations, and the rs315831750 leads to a non-coding triplet just upstream of the initiation ATG in the reference sequence becoming a coding initiation ATG, thus the variant will have two adjacent methionines instead of just one in the reference sequence. The role of these SNPs in affecting the efficiency of transcription or translation of the PLIN1 gene remains unknown. A seemingly silent SNP may not necessarily be neutral and may enact a function (Kimchi-Sarfaty et al., 2007).

We calculated the allelic and genotypic frequencies and genetic diversity parameters of the three SNPs in the six chicken populations. The six populations in this study included both indigenous and commercial breeds, showing that these SNPs exist in different genetic background populations, not only in a discrete population. The HB chicken breed was not in dynamic equilibrium for the three loci, which perhaps was a consequence of long-term artificial selection. Since 1996, HB chickens have been subjected to 16 generations of breeding and is famous for its production performance (Gu and Dai, 2011). Therefore, in this study all the associations with carcass traits and IMF in chickens did not include the HB chickens. The IMF, PLM, PBM, and BW in HB chickens were greater than in other chickens, and the BW, PCW, PSEW, PAF, and IMF were lower in MB than in other chickens (Table 2).

Polymorphic information content >0.5 indicates high polymorphism, 0.25<PIC<0.5 reflects medium polymorphism, and PIC<0.25 represents low polymorphism (Vaiman et al., 1994). In this study, the PIC ranged from 0.2471 to 0.3736 (Table 3), which reflected an intermediate genetic diversity of the chicken PLIN1 gene in all chicken populations studied. Thus, these loci may serve as potential molecular markers, but whether they are appropriate as molecular markers for marker-assisted selection should be validated in more groups.

We constructed haplotypes for all six chicken populations. The current study used PHASE software to infer the most likely haplotype pairs for each individual, with probabilities ranging from 5% to 100%. For association analysis between a candidate gene and the studied phenotypes, there are differing views (Zhao et al., 2009). Some researchers propose that haplotype analysis may be more powerful than single SNP analysis due to accumulated allelic effects at different loci and/or interactions between these loci, whereas others suggest that if a connection between SNPs and a certain phenotype is driven by only a single SNP, then a single SNP approach will be superior to the haplotype-based approach (Morris and Kaplan, 2002; Clark, 2004).

The single SNP genotype had a large effect on PBM, PLM and PAF (p<0.05), and the diplotype had a large effect on PEW and PBM (p<0.05). All of these results indicate that additive effects of a single SNP might not explain our results. More likely, there are interactions among these loci or even linkage disequilibrium. We performed linkage disequilibrium analysis for the three SNPs in all chickens, but did not include HB chickens. The pattern of pairwise LD structure as measured by D′ and r² between the SNPs was calculated using the program Haploview 4.1, The D′ values between SNP1 and SNP2, SNP2 and SNP3, and SNP1 and SNP3 were 0.81, 0.92, and 0.66, respectively. Moreover, r² for them were 0.601, 0.539, and 0.309, respectively. This showed that SNP1 and SNP2, SNP2, and SNP3 had high LD, the SNP1 and SNP3 had medium LD.

We observed significant relationships between single SNP genotype and traits PBM, PLM, and PAF, and between haplotype pairs and traits PEW and PBM. The homozygous genotype CC chickens had the greatest PBM, the heterozygous genotype CT had the greatest PLM and smallest PAF at rs315831750. Birds with the H3H7 had greater PEW and H2H2 had greater PLM than other diplotypes. Thus, we infer that diplotype H3H7 and H2H2 are ideal genetic markers for selecting greater PEW and PLM in chickens, respectively. Meanwhile, as there were significant population effects for PEW and PBM, it is necessary to analyze the effects of PLIN1 polymorphisms using populations with different genetic backgrounds.

ACKNOWLEDGMENTS

We thank Xiao-Song Jiang and Hua-Rui Du for their help in managing the birds and collecting the data. This work was financially supported by the China Agricultural Research System (CARS-41) and the National Natural Science Foundation of China (Grant No: 31402070).

References

Barrett JC. 2009. Haploview: Visualization and analysis of snp genotype data. Cold Spring Harb Protoc http://dx.doi.org/10.1101/pdb.ip71.

Benjamini Y, Hochberg Y. 1995;Controlling the false discovery rate: A practical and powerful approach to multiple testing. J R Stat Soc Series B (Methodological) 57:289–300.

Bickel PE, Tansey JT, Welte MA. 2009;Pat proteins, an ancient family of lipid droplet proteins that regulate cellular lipid stores. Biochim Biophys Acta Cell Biol Lipids 1791:419–440.

Clark AG. 2004;The role of haplotypes in candidate gene studies. Genet Epidemiol 27:321–333.

Fan HJ, Liu XL, Wang J, Hou SS. 2011;Analysis of plin gene pcr-rflp polymorphism and its association with carcass and adipose traits in meat ducks. Acta Agriculturae Boreali-Occidentalis Sinica 6:004.

Fan Y, Zan L, Wang H, Yang Y. 2010;Study on the relationship between polymorphism of plin gene and carcass and meat quality traits in qinchuan cattle. Chin J Anim Vet Sci 41:268–273.

Gandolfi G, Mazzoni M, Zambonelli P, Lalatta-Costerbosa G, Tronca A, Russo V, Davoli R. 2011;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 88:631–637.

Gao ZY, Lin PP, Yuan YN, Zhou SS, Liu BF, Liu JH, Linag C, Qiao LY, Liu WZ. 2012;Study on the polymorphism of plin gene and its association with tail and slaughter traits in sheep. J Shanxi Agric Univ (Natural Science Edition) 32:158–164.

Gu HB, Dai YL. 2011;China “new generalized meaning” agricultural test lead taiwan hongbiao native chicken. Guide to Chinese Poultry 10:41–41.

Kern PA, Di Gregorio G, Lu T, Rassouli N, Ranganathan G. 2004;Perilipin expression in human adipose tissue is elevated with obesity. J Clin Endocrinol Metab 89:1352–1358.

Kimchi-Sarfaty C, Oh JM, Kim I-W, Sauna ZE, Calcagno AM, Ambudkar SV, Gottesman MM. 2007;A “silent” polymorphism in the MDR1 gene changes substrate specificity. Science 315(5811):525–528.

Kimmel AR, Brasaemle DL, McAndrews-Hill M, Sztalryd C, Londos C. 2010;Adoption of PERILIPIN as a unifying nomenclature for the mammalian PAT-family of intracellular lipid storage droplet proteins. J Lipid Res 51:468–471.

Lande R, Thompson R. 1990;Efficiency of marker-assisted selection in the improvement of quantitative traits. Genetics 124:743–756.

Lei QX, Zhou Y, Han HX, Li GM, Cao DG, Li FW, Lu Y. 2011;Polymorphism in plin gene intron 5 and its association with carcass and fatness traits in chicken. China Anim Husb Vet Med 12:036.

Londos C, Sztalryd C, Tansey JT, Kimmel AR. 2005;Role of pat proteins in lipid metabolism. Biochimie 87:45–49.

Marcinkiewicz A, Gauthier D, Garcia A, Brasaemle DL. 2006;The phosphorylation of serine 492 of perilipin a directs lipid droplet fragmentation and dispersion. J Biol Chem 281:11901–11909.

Miura S, Gan J-W, Brzostowski J, Parisi MJ, Schultz CJ, Londos C, Oliver B, Kimmel AR. 2002;Functional conservation for lipid storage droplet association among perilipin, ADRP, and TIP47 (PAT)-related proteins in mammals, drosophila, and dictyostelium. J Biol Chem 277:32253–32257.

Morris RW, Kaplan NL. 2002;On the advantage of haplotype analysis in the presence of multiple disease susceptibility alleles. Genet Epidemiol 23:221–233.

Nei M, Roychoudhury A. 1974;Sampling variances of heterozygosity and genetic distance. Genetics 76:379–390.

Petracci M, Cavani C. 2012;Muscle growth and poultry meat quality issues. Nutrients 4:1–12.

Rosa MD, Gambacorta A, Nicolaus B, Chappe B, Albrecht P. 1983;Isoprenoid ethers; backbone of complex lipids of the archaebacterium sulfolobus solfataricus. Biochim Biophy Acta (BBA)-Lipids and Lipid Metabolism 753:249–256.

Rubin CJ, Zody MC, Eriksson J, Meadows JRS, Sherwood E, Webster MT, Jiang L, Ingman M, Sharpe T, Ka S, Hallbook F, Besnier F, Carlborg O, Bed’hom B, Tixier-Boichard M, Jensen P, Siegel P, Lindblad-Toh K, Andersson L. 2010;Whole-genome resequencing reveals loci under selection during chicken domestication. Nature 464:587–591.

Ruiz JR, Larrarte E, Margareto J, Ares R, Alkorta P, Labayen I. 2011;Preliminary findings on the role of PLIN1 polymorphisms on body composition and energy metabolism response to energy restriction in obese women. Br J Nutr 106:486–490.

Song W, Yu H, Lin Y, Sun K, Zhang Y, Song Y, Hui R, Chen J. 2015;A functional variant in the exon 5 of PLIN1 reduces risk of central obesity by possible regulation of lipid storage. Biochem Biophys Res Commun 456:896–900.

Stephens M, Smith NJ, Donnelly P. 2001;A new statistical method for haplotype reconstruction from population data. Am J Hum Genet 68:978–989.

Sun HJ, Wang GL. 2006;Study progress on the evaluation system of quality chicken meat. China Poult 28:38–42.

Vaiman D, Mercier D, Moazami-Goudarzi K, Eggen A, Ciampolini R, Lépingle A, Velmala R, Kaukinen J, Varvio S-L, Martin P, Levéziel H, Guérin G. 1994;A set of 99 cattle microsatellites: Characterization, synteny mapping, and polymorphism. Mamm Genome 5:288–297.

Wang Y, Shu D, Li L, Qu H, Yang C, Zhu Q. 2007;Identification of single nucleotide polymorphism of h-fabp gene and its association with fatness traits in chickens. Asian Australas J Anim Sci 20:1812–1819.

Wang Y, Xiao LH, Zhao XL, Liu YP, Zhu Q. 2014;Identification of snps in cellular retinol binding protein 1 and cellular retinol binding protein 3 genes and their associations with laying performance traits in erlang mountainous chicken. Asian Australas J Anim Sci 27:1075–1081.

Wang Y, Zhu Q, Yang L, Liu YP. 2012;Ontogenic expression pattern and genetic polymorphisms of the fatty acid transport protein 4 (FATP4) gene in Chinese chicken populations. Int J Mol Sci 13:6820–6835.

Zhang HL, Fan HJ, Liu XL, Wu Y, Hou SS. 2013;Molecular cloning of the perilipin gene and its association with carcass and fat traits in chinese ducks. Genet Mol Res 12:1582–1592.

Zhao XL, Liu YP, Jiang XS, Du HR, Zhu Q. 2009;Association of polymorphisms of chicken adipose differentiation-related protein gene with carcass traits. J Poult Sci 46:87–94.

Zhao XL, Liu YP, Luo Y, Zhou Y, Zhu Q. 2009;Study on the relationship between developmental variants of plin gene expression and fatness traits in chickens. Chinese J Anim Vet Sci 40:149–154.

Zhao XL, Siegel PB, Liu YP, Wang Y, Gilbert ER, Zhu Q, Zhang L. 2012;Housing system affects broiler characteristics of local chinese breed reciprocal crosses. Poult Sci 91:2405–2410.

Zhou Y, Lei QX, Li FW, Gao JB, Liu W, Lu Y, Chao DG. 2014;Association on single nucleotide polymorphism of perilipin gene (PLIN) with carcass and fatness traits in luqin chicken (Gallus gallus). J Agric Biotechnol 22:1002–1008.

Zhou Y, Liu YP, Jiang XS, Du HR, Zhu Q. 2009;Study on association of single nucleotide polymorphism of mef2a gene with carcass traits in chicken. Chinese J Anim Vet Sci 8:010.

Article information Continued

This is an open-access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/3.0/), wwhich permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.

Primer number	Forward primer (5′- 3′)	Reverse primer (5′- 3′)	Annealing temperature (°C)	Product length (bp)	Coverage
1	TCCAGCAGAGGAGGAGTGTGT	TTCTGCAAAGGCTGATTCTT	56.9	239	5′flanking region
2	CAACCCCTGCACATCGT	CAGCACGGTATGAGCTATTCC	57.5	200	Exon 1
3	GTCCCCATCTGATGCAC	CCACCACATCCCTCCTGATTG	57.5	245	Exon 2
4	GGCACAGCTCCCAGCCTTAAT	GGCCCAGAGGAGAAATAC	56.1	221	Exon 3
5	GCTTGCCTCCAAACTGA	AACAGGTACTCCATCAGCTTC	56.6	231	Exon 4
6	TGAAGCCCAAGCAGACGTGT	GGGGATCCAGACGACCAGTTC	60.8	221	Exon 5
7	CCAAGCACCGGTTTCTGTT	TGACCCCATCTCACCATTATC	58.1	209	Exon 6
8	GGGTGCCGTAAATGACGCT	GTGCTTAGGGGCACTGATGAG	59.8	238	Exon 7
9	CAGGCAGTCCAAGGACGA	GGTACAGATTGCTGCGTATGA	59.0	203	Exon 8
10	CTCTTCATACGCAGCAATCTG	CAATTTGGGACCGAGGTT	58.7	214	Exon 8

Line	Size	Traits
DH	40	2,160.88±25.09^b	89.68±0.48^a	83.66±0.56^a	70.01±0.50^b	8.31±0.21^d	11.85±0.24^c	2.79±0.25^b	1.88±0.19^c
HB	40	2,822.73±25.09^a	86.71±0.48^b	79.83±0.56^b	62.73±0.50^c	15.74±0.21^a	22.45±0.24^a	2.21±0.25^bc	3.76±0.19^a
MB	40	1,501.20±25.09^c	85.78±0.48^b	78.95±0.56^b	74.05±0.50^a	12.92±0.21^b	19.42±0.24^b	1.09±0.25^c	1.18±0.19^c
EM	40	2,676.18±25.09^a	85.88±0.48^b	79.85±0.56^b	61.81±0.50^c	9.33±0.21^c	11.28±0.24^c	4.69±0.25^a	2.99±0.19^b
BH1	40	2,257.53±25.09^b	85.98±0.48^b	79.29±0.56^b	61.23±0.50^c	8.64±0.21^cd	11.25±0.24^c	4.81±0.25^a	2.93±0.19^b
BH2	40	2,321.50±25.09^b	86.75±0.48^b	79.02±0.56^b	60.90±0.50^c	8.75±0.21^cd	11.54±0.24^c	5.30±0.25^a	2.93±0.19^b
p-value³		≤0.0001	≤0.0001	≤0.0001	≤0.0001	≤0.0001	≤0.0001	≤0.0001	≤0.0001

Table 3

Genotype and allele frequencies and genetic diversity parameters for the PLIN1 gene in six breeds

Mutation	SNP	Line	Size	Genotype frequencies2			Allele frequencies		x²	Ho	He	Ne	PIC

				BB	AB	AA	B	A
C/T1	rs315831750	DH	40	0.3000	0.6000	0.1000	0.6000	0.4000	2.500	0.5200	0.4800	1.9231	0.3648
		HB	40	0.1250	0.8000	0.0750	0.5250	0.4750	14.590**	0.5013	0.4988	1.9950	0.3744
		MB	40	0.2500	0.6250	0.1250	0.5625	0.4375	2.913	0.5078	0.4922	1.9692	0.3711
		EM	40	0.6000	0.2750	0.1250	0.7375	0.2625	3.358	0.6128	0.3872	1.6318	0.3122
		BH1	40	0.3000	0.5750	0.1250	0.5875	0.4125	1.389	0.5153	0.4847	1.9406	0.3672
		BH2	40	0.4250	0.5500	0.0250	0.7000	0.3000	3.832	0.5800	0.4200	1.7241	0.3318
G/A1	rs313726543	DH	40	0.3750	0.4500	0.1750	0.6000	0.4000	0.156	0.5200	0.4800	1.9231	0.3648
		HB	40	0.2000	0.6750	0.1250	0.5375	0.4625	5.116*	0.5028	0.4972	1.9888	0.3736
		MB	40	0.3500	0.5500	0.1000	0.6250	0.3750	1.202	0.5313	0.4688	1.8824	0.3589
		EM	40	0.6250	0.2750	0.1000	0.7625	0.2375	2.318	0.6378	0.3622	1.5679	0.2966
		BH1	40	0.1750	0.6000	0.2250	0.4750	0.5250	1.648	0.5013	0.4988	1.9950	0.3744
		BH2	40	0.3000	0.6000	0.1000	0.6000	0.4000	2.500	0.5200	0.4800	1.9231	0.3648
T/C1	rs80724063	DH	40	0.5750	0.4250	0.0000	0.7875	0.2125	2.913	0.6653	0.3347	1.5031	0.2787
		HB	40	0.5000	0.5000	0.0000	0.7500	0.2500	4.444*	0.6250	0.3750	1.6000	0.3047
		MB	40	0.4000	0.4750	0.1250	0.6375	0.3625	0.031	0.5378	0.4622	1.8594	0.3554
		EM	40	0.6750	0.3000	0.0250	0.8250	0.1750	0.061	0.7113	0.2888	1.4060	0.2471
		BH1	40	0.4250	0.3750	0.2000	0.6125	0.3875	1.764	0.5253	0.4747	1.9036	0.3620
		BH2	40	0.4500	0.5000	0.0500	0.7000	0.3000	1.451	0.5800	0.4200	1.7241	0.3318

SNP, single nucleotide polymorphism; Ho, observed homozygosity; He, expected heterozygosity; Ne, effective allele numbers; PIC, polymorphism information content.

Minor allele is in bold.

AA minor allele homozygote, AB heterozygote, BB major allele homozygote.

Asterisks represent loci that deviate from the Hardy-Weinberg equilibrium.

^**

p<0.01;

p<0.05 (x²-test).

Location	Trait	Genotype1			p-value2	pFDR

		BB	AB	AA
rs315831750	PBM (%)	9.53±0.34^b	11.13±0.26^a	11.23±0.64^a	0.0007	0.0113
	PLM (%)	13.04±0.53^b	15.48±0.40^a	15.04±1.01^ab	0.0013	0.0113
	PAF (%)	4.20±0.27^a	3.10±0.21^b	3.27±0.52^ab	0.0059	0.0342

Trait	H1H7 (28)1	H2H2 (10)	H2H7 (53)	H3H7 (12)	H6H7 (13)	H7H7 (59)	p-value2	pFDR
PEW (%)	65.21±1.19^b	64.53±1.98^b	66.70±0.86^ab	71.81±1.81^a	64.38±1.74^b	64.29±0.82^b	0.007	0.011
PBM (%)	9.35±0.39^ab	11.27±0.65^a	9.79±0.28^ab	10.22±0.59^ab	9.77±0.57^ab	9.16±0.27^b	0.048	0.034

Line	Size	Traits

		BW (g)	PCW (%)	PSEW (%)	PEW (%)	PBM (%)	PLM (%)	PAF (%)	IMF (%)
DH	40	2,160.88±25.09^b	89.68±0.48^a	83.66±0.56^a	70.01±0.50^b	8.31±0.21^d	11.85±0.24^c	2.79±0.25^b	1.88±0.19^c
HB	40	2,822.73±25.09^a	86.71±0.48^b	79.83±0.56^b	62.73±0.50^c	15.74±0.21^a	22.45±0.24^a	2.21±0.25^bc	3.76±0.19^a
MB	40	1,501.20±25.09^c	85.78±0.48^b	78.95±0.56^b	74.05±0.50^a	12.92±0.21^b	19.42±0.24^b	1.09±0.25^c	1.18±0.19^c
EM	40	2,676.18±25.09^a	85.88±0.48^b	79.85±0.56^b	61.81±0.50^c	9.33±0.21^c	11.28±0.24^c	4.69±0.25^a	2.99±0.19^b
BH1	40	2,257.53±25.09^b	85.98±0.48^b	79.29±0.56^b	61.23±0.50^c	8.64±0.21^cd	11.25±0.24^c	4.81±0.25^a	2.93±0.19^b
BH2	40	2,321.50±25.09^b	86.75±0.48^b	79.02±0.56^b	60.90±0.50^c	8.75±0.21^cd	11.54±0.24^c	5.30±0.25^a	2.93±0.19^b
p-value³		≤0.0001	≤0.0001	≤0.0001	≤0.0001	≤0.0001	≤0.0001	≤0.0001	≤0.0001

Haplotype	SNP			Frequency1

	SNP1	SNP2	SNP3
H1	C	G	C	0.1136
H2	C	G	T	0.2056
H3	C	A	C	0.0400
H4	C	A	T	0.0008
H5	T	G	C	0.0010
H6	T	G	T	0.0673
H7	T	A	C	0.5579
H8	T	A	T	0.0137