Polymorphisms in the Perilipin gene may affect carcass traits of Chinese meat-type chickens.
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
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 cmx60 cmx30 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 cmx55 cmx36 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).
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[degrees]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[degrees]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 [micro]L containing 200 ng genomic DNA, 1 [micro]L each of forward and reverse primers (10 pmol/[micro]L), and 25 [micro]L 2xTaq PCR Master Mix (Beijing Tianwei Biology Technique Corporation, Beijing, China). The reaction mixture was predenatured for 4 min at 94[degrees]C, followed by 35 cycles at 94[degrees]C for 30 s, 57[degrees]C for 30 s, 72[degrees]C for 1 min, and ended with a full extension cycle at 72[degrees]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).
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.sup.2]) 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.sup.2]-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:
[Y.sub.ijk] = [mu] + [P.sub.i] + [S.sub.j] + [G.sub.k] + [(PG).sub.ik] + [(PS).sub.ij] + [(GS).sub.ij] + [(PSG).sub.ijk] + [e.sub.ijk]
Where [Y.sub.ijk] = the dependent variable; [mu] = the population mean; [P.sub.i] = fixed effect of population; [S.sub.j] = fixed effect of sex; [G.sub.k] = fixed effect of genotype/diplotype. [(PG).sub.ik] = interaction effect between population and genotype/diplotype; [(PS).sub.ij] = interaction effect between population and sex; [(SG).sub.jk] = interaction effect between genotype/diplotype and sex; [(PSG).sub.ijk] = interaction effect between population and genotype/diplotype and sex; [e.sub.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.
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).
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.
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.
The [x.sup.2]-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.
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).
The linkage disequilibrium analysis between the three single nucleotide polymorphisms
The pattern of pairwise LD structure as measured by D' and [r.sup.2] 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.sup.2] for them were 0.601, 0.539, and 0.309, respectively.
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.sup.2] 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.sup.2] 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.
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).
Barrett, J. C. 2009. Haploview: Visualization and analysis of snp genotype data. Cold Spring Harb Protoc. http://dx.doi.org/10.1101/pdb.ip71
Benjamini, Y. and Y. Hochberg. 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, P. E., J. T. Tansey, and M. A. Welte. 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, A. G. 2004. The role of haplotypes in candidate gene studies. Genet. Epidemiol. 27:321-333.
Fan, H. J., X. L. Liu, J. Wang, and S. S. Hou. 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., L. Zan, H. Wang, and Y. Yang. 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., M. Mazzoni, P. Zambonelli, G. Lalatta-Costerbosa, A. Tronca, V. Russo, and R. Davoli. 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, Z. Y, P. P. Lin, Y. N. Yuan, S. S. Zhou, B. F. Liu, J. H. Liu, C. Linag, L. Y. Qiao, and W. Z. Liu. 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, H. B. and Y. L. Dai. 2011. China "new generalized meaning" agricultural test lead taiwan hongbiao native chicken. Guide to Chinese Poultry 10:41-41.
Kern, P. A., G. Di Gregorio, T. Lu, N. Rassouli, and G Ranganathan. 2004. Perilipin expression in human adipose tissue is elevated with obesity. J. Clin. Endocrinol. Metab. 89: 1352-1358.
Kimchi-Sarfaty, C., J. M. Oh, I.-W. Kim, Z. E. Sauna, A. M. Calcagno, S. V. Ambudkar, and M. M. Gottesman. 2007. A" silent" polymorphism in the MDR1 gene changes substrate specificity. Science 315(5811):525-528.
Kimmel, A. R., D. L. Brasaemle, M. McAndrews-Hill, C. Sztalryd, and C. Londos. 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. and R. Thompson. 1990. Efficiency of marker-assisted selection in the improvement of quantitative traits. Genetics 124:743-756.
Lei, Q. X., Y. Zhou, H. X. Han, G. M. Li, D. G. Cao, F. W. Li, and Y Lu. 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., C. Sztalryd, J. T. Tansey, and A. R. Kimmel. 2005. Role of pat proteins in lipid metabolism. Biochimie 87:45-49.
Marcinkiewicz, A., D. Gauthier, A. Garcia, and D. L. Brasaemle. 2006. The phosphorylation of serine 492 of perilipin a directs lipid droplet fragmentation and dispersion. J. Biol. Chem. 281: 11901-11909.
Miura, S., J.-W. Gan, J. Brzostowski, M. J. Parisi, C. J. Schultz, C. Londos, B. Oliver, and A. R. Kimmel. 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, R. W. and N. L. Kaplan. 2002. On the advantage of haplotype analysis in the presence of multiple disease susceptibility alleles. Genet. Epidemiol. 23:221-233.
Nei, M. and A. Roychoudhury. 1974. Sampling variances of heterozygosity and genetic distance. Genetics 76:379-390.
Petracci, M. and C. Cavani. 2012. Muscle growth and poultry meat quality issues. Nutrients 4:1-12.
Rosa, M. D., A. Gambacorta, B. Nicolaus, B. Chappe, and P. Albrecht. 1983. Isoprenoid ethers; backbone of complex lipids of the archaebacterium sulfolobus solfataricus. Biochim. Biophy. Acta (BBA)-Lipids and Lipid Metabolism 753:249-256.
Rubin, C. J., M. C. Zody, J. Eriksson, J. R. S. Meadows, E. Sherwood, M. T. Webster, L. Jiang, M. Ingman, T. Sharpe, S. Ka, F. Hallbook, F. Besnier, O. Carlborg, B. Bed'hom, M. Tixier-Boichard, P. Jensen, P. Siegel, K. Lindblad-Toh, and L. Andersson. 2010. Whole-genome resequencing reveals loci under selection during chicken domestication. Nature 464:587-591.
Ruiz, J. R., E. Larrarte, J. Margareto, R. Ares, P. Alkorta, and I. Labayen. 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., H. Yu, Y. Lin, K. Sun, Y. Zhang, Y. Song, R. Hui, and J. Chen. 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., N. J. Smith, and P. Donnelly. 2001. A new statistical method for haplotype reconstruction from population data. Am. J. Hum. Genet. 68:978-989.
Sun, H. J. and G. L. Wang, 2006. Study progress on the evaluation system of quality chicken meat. China Poult. 28:38-42.
Vaiman, D., D. Mercier, K. Moazami-Goudarzi, A. Eggen, R. Ciampolini, A. Lepingle, R. Velmala, J. Kaukinen, S.-L. Varvio, P. Martin, H. Leveziel, and G. Guerin. 1994. A set of 99 cattle microsatellites: Characterization, synteny mapping, and polymorphism. Mamm. Genome 5:288-297.
Wang, Y, D. Shu, L. Li, H. Qu, C. Yang, and Q. Zhu. 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, L. H. Xiao, X. L. Zhao, Y. P. Liu, and Q. Zhu. 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, Q. Zhu, L. Yang, and Y. P. Liu. 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, H. L., H. J. Fan, X. L. Liu, Y. Wu, and S. S. Hou. 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, X. L., Y. P. Liu, X. S. Jiang, H. R. Du, and Q. Zhu. 2009. Association of polymorphisms of chicken adipose differentiation-related protein gene with carcass traits. J. Poult. Sci. 46:87-94.
Zhao, X. L., Y. P. Liu, Y. Luo, Y. Zhou, and Q. Zhu. 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, X. L., P. B. Siegel, Y. P. Liu, Y. Wang, E. R. Gilbert, Q. Zhu, and L. Zhang. 2012. Housing system affects broiler characteristics of local chinese breed reciprocal crosses. Poult. Sci. 91:2405-2410.
Zhou, Y, Q. X. Lei, F. W. Li, J. B. Gao, W. Liu, Y. Lu, and D. G Chao. 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, Y. P. Liu, X. S. Jiang, H. R. Du, and Q. Zhu. 2009. Study on association of single nucleotide polymorphism of mef2a gene with carcass traits in chicken. Chinese J. Anim. Vet. Sci. 8: 010.
Submitted Jul. 23, 2014; Revised Aug. 20, 2014; Accepted Jan. 12, 2015
Lu Zhang (1), Qing Zhu (1), Yiping Liu (1), Elizabeth R. Gilbert (2), Diyan Li (1), Huadong Yin (1), Yan Wang (1), Zhiqin Yang (1), Zhen Wang (1), Yuncong Yuan (1), and Xiaoling Zhao *
Department of Animal Science, Sichuan Agricultural University, Ya'an Sichuan 625014, China
* Corresponding Author: Xiaoling Zhao. Tel: +86-835-2886080, Fax: +86-835-2886080, E-mail: firstname.lastname@example.org
(1) Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, Sichuan Agricultural University, Chengdu, Sichuan 611130, China.
(2) Department of Animal and Poultry Sciences, Virginia Polytechnic Institute and State University, Blacksburg, VA24061, USA.
Table 1. Primers for screening SNPs in the chicken PLIN1 gene Primer Forward primer Reverse primer number (5'-3') (5'-3') 1 TCCAGCAGAGGAGGAGTGTGT TTCTGCAAAGGCTGATTCTT 2 CAACCCCTGCACATCGT CAGCAC GGTATGAGC TATTCC 3 GTCCCCATCTGATGCAC CCACCACATCCCTCCTGATTG 4 GGCACAGCTCCCAGCCTTAAT GGCCCAGAGGAGAAATAC 5 GCTTGCCTCCAAACTGA AACAGGTACTCCATCAGCTTC 6 TGAAGCCCAAGCAGACGTGT GGGGATCCAGACGACCAGTTC 7 CCAAGCACCGGTTTCTGTT TGACCCCATCTCACCATTATC 8 GGGTGCCGTAAATGACGCT GTGCTTAGGGGCACTGATGAG 9 CAGGCAGTCCAAGGACGA GGTACAGATT GCTGCGTATGA 10 CTCTTCATACGCAGCAATCTG CAATTTGGGACCGAGGTT Primer Annealing Product Coverage number temperature length ([degrees]C) (bp) 1 56.9 239 5'flanking region 2 57.5 200 Exon 1 3 57.5 245 Exon 2 4 56.1 221 Exon 3 5 56.6 231 Exon 4 6 60.8 221 Exon 5 7 58.1 209 Exon 6 8 59.8 238 Exon 7 9 59.0 203 Exon 8 10 58.7 214 Exon 8 SNPs, single nucleotide polymorphisms. Table 2. Population effects on carcass traits and IMF in chickens Line Size BW (g) PCW (%) DH 40 2,160.88 [+ or -] 25.09 (b) 89.68 [+ or -] 0.48# (a) HB 40 2,822.73 [+ or -] 25.09# (a) 86.71 [+ or -] 0.48 (b) MB 40 1,501.20 [+ or -] 25.09@ (c) 85.78 [+ or -] 0.48@ (b) EM 40 2,676.18 [+ or -] 25.09 (a) 85.88 [+ or -] 0.48 (b) BH1 40 2,257.53 [+ or -] 25.09 (b) 85.98 [+ or -] 0.48 (b) BH2 40 2,321.50 [+ or -] 25.09 (b) 86.75 [+ or -] 0.48 (b) p-value [less than or equal to] [less than or equal to] (3) 0.0001 0.0001 Line Traits PSEW (%) PEW (%) DH 83.66 [+ or -] 0.56# (a) 70.01 [+ or -] 0.50 (b) HB 79.83 [+ or -] 0.56 (b) 62.73 [+ or -] 0.50 (c) MB 78.95 [+ or -] 0.56@ (b) 74.05 [+ or -] 0.50# (a) EM 79.85 [+ or -] 0.56 (b) 61.81 [+ or -] 0.50 (c) BH1 79.29 [+ or -] 0.56 (b) 61.23 [+ or -] 0.50 (c) BH2 79.02 [+ or -] 0.56 (b) 60.90 [+ or -] 0.50@ (c) p-value [less than or equal to] [less than or equal to] (3) 0.0001 0.0001 Line PBM (%) PLM (%) DH 8.31 [+ or -] 0.21@ (d) 11.85 [+ or -] 0.24 (c) HB 15.74 [+ or -] 0.21# (a) 22.45 [+ or -] 0.24# (a) MB 12.92 [+ or -] 0.21 (b) 19.42 [+ or -] 0.24 (b) EM 9.33 [+ or -] 0.21 (c) 11.28 [+ or -] 0.24 (c) BH1 8.64 [+ or -] 0.21 (c)d 11.25 [+ or -] 0.24@ (c) BH2 8.75 [+ or -] 0.21 (c)d 11.54 [+ or -] 0.24 (c) p-value [less than or equal to] [less than or equal to] (3) 0.0001 0.0001 Line PAF (%) IMF (%) DH 2.79 [+ or -] 0.25 (b) 1.88 [+ or -] 0.19 (c) HB 2.21 [+ or -] 0.25 (bc) 3.76 [+ or -] 0.19# (a) MB 1.09 [+ or -] 0.25@ (c) 1.18 [+ or -] 0.19@ (c) EM 4.69 [+ or -] 0.25 (a) 2.99 [+ or -] 0.19 (b) BH1 4.81 [+ or -] 0.25 (a) 2.93 [+ or -] 0.19 (b) BH2 5.30 [+ or -] 0.25# (a) 2.93 [+ or -] 0.19 (b) p-value [less than or equal to] [less than or equal to] (3) 0.0001 0.0001 IMF, intramuscular fat content; BW, body weight; PCW, percentage of carcass weight; PSEW, percentage of semi-eviscerated weight; PEW, percentage of eviscerated weight; PBM, percentage of breast muscle; PLM, percentage of leg muscle; PAF, percentage of abdominal fat; DH, Daheng chicken; HB, Hongbiao chicken; MB, Mountainous Black-bone chicken; EM, Erlang Mountainous chickens; BH1, line of BH1 from Banghe company; BH2, line of BH2 from Banghe company. (1)Probability of the F-test for Type III SS. Values represent least squares means plus or minus standard error. Values in a row without same superscript letter differ (p < 0.05, Tukey's test). The bold numbers represent the greatest value; underlined numbers represent the smallest value. Note: Numbers represent the greatest value are indicated with #. Note: Numbers represent the smallest value are indicated with @. Table 3. Genotype and allele frequencies and genetic diversity parameters for the PLIN1 gene in six breeds Mutation SNP Line Size Genotype frequencies (2) BB AB AA C/T (1) rs315831750 DH 40 0.3000 0.6000 0.1000 HB 40 0.1250 0.8000 0.0750 MB 40 0.2500 0.6250 0.1250 EM 40 0.6000 0.2750 0.1250 BH1 40 0.3000 0.5750 0.1250 BH2 40 0.4250 0.5500 0.0250 G/A (1) rs313726543 DH 40 0.3750 0.4500 0.1750 HB 40 0.2000 0.6750 0.1250 MB 40 0.3500 0.5500 0.1000 EM 40 0.6250 0.2750 0.1000 BH1 40 0.1750 0.6000 0.2250 BH2 40 0.3000 0.6000 0.1000 T/C (1) rs80724063 DH 40 0.5750 0.4250 0.0000 HB 40 0.5000 0.5000 0.0000 MB 40 0.4000 0.4750 0.1250 EM 40 0.6750 0.3000 0.0250 BH1 40 0.4250 0.3750 0.2000 BH2 40 0.4500 0.5000 0.0500 Mutation SNP Line Allele [x.sup.2] frequencies B A C/T (1) rs315831750 DH 0.6000 0.4000 2.500 HB 0.5250 0.4750 14.590 ** MB 0.5625 0.4375 2.913 EM 0.7375 0.2625 3.358 BH1 0.5875 0.4125 1.389 BH2 0.7000 0.3000 3.832 G/A (1) rs313726543 DH 0.6000 0.4000 0.156 HB 0.5375 0.4625 5.116 * MB 0.6250 0.3750 1.202 EM 0.7625 0.2375 2.318 BH1 0.4750 0.5250 1.648 BH2 0.6000 0.4000 2.500 T/C (1) rs80724063 DH 0.7875 0.2125 2.913 HB 0.7500 0.2500 4.444 * MB 0.6375 0.3625 0.031 EM 0.8250 0.1750 0.061 BH1 0.6125 0.3875 1.764 BH2 0.7000 0.3000 1.451 Mutation SNP Line Ho He Ne PIC C/T (1) rs315831750 DH 0.5200 0.4800 1.9231 0.3648 HB 0.5013 0.4988 1.9950 0.3744 MB 0.5078 0.4922 1.9692 0.3711 EM 0.6128 0.3872 1.6318 0.3122 BH1 0.5153 0.4847 1.9406 0.3672 BH2 0.5800 0.4200 1.7241 0.3318 G/A (1) rs313726543 DH 0.5200 0.4800 1.9231 0.3648 HB 0.5028 0.4972 1.9888 0.3736 MB 0.5313 0.4688 1.8824 0.3589 EM 0.6378 0.3622 1.5679 0.2966 BH1 0.5013 0.4988 1.9950 0.3744 BH2 0.5200 0.4800 1.9231 0.3648 T/C (1) rs80724063 DH 0.6653 0.3347 1.5031 0.2787 HB 0.6250 0.3750 1.6000 0.3047 MB 0.5378 0.4622 1.8594 0.3554 EM 0.7113 0.2888 1.4060 0.2471 BH1 0.5253 0.4747 1.9036 0.3620 BH2 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. (1) Minor allele is in bold. (2) 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 ([chi square]-test). Table 4. Effect of single SNP genotype on carcass traits Location Trait Genotype (1) BB AB AA rs315831750 PBM 9.53 [+ or -] 11.13 [+ or -] 11.23 [+ or -] (%) 0.34 (b) 0.26 (a) 0.64 (a) PLM 13.04 [+ or -] 15.48 [+ or -] 15.04 [+ or -] (%) 0.53 (b) 0.40 (a) 1.01 (ab) PAF 4.20 [+ or -] 3.10 [+ or -] 3.27 [+ or -] (%) 0.27 (a) 0.21 (b) 0.52 (ab) Location Trait p-value (2) pFDR rs315831750 PBM 0.0007 0.0113 (%) PLM 0.0013 0.0113 (%) PAF 0.0059 0.0342 (%) SNP, single nucleotide polymorphism; pFDR, the positive false discovery rate; PBM, percentage of breast muscle; PLM, percentage of leg muscle; PAF, percentage of abdominal fat. (1) AA minor allele homozygote, AB heterozygote, BB major allele homozygote. (2) Probability of the F-test for Type III SS. The least square means within a row with different superscripts differ significantly (Tukey's test, p < 0.05). Table 5. Haplotypes constructed with single nucleotide polymorphisms (SNP) and their frequencies Haplotype SNP SNP1 SNP2 SNP3 Frequency (1) 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 (1) Haplotype frequency estimates with PHASE 2.0.2. Table 6. Diplotype effects on PEW and PBM in chickens Trait H1H7 (28) (1) H2H2 (10) PEW (%) 65.21 [+ or -] 1.19 (b) 64.53 [+ or -] 1.98 (b) PBM (%) 9.35 [+ or -] 0.39 (ab) 11.27 [+ or -] 0.65# (a) Trait H2H7 (53) H3H7 (12) PEW (%) 66.70 [+ or -] 0.86 (ab) 71.81 [+ or -] 1.81# (a) PBM (%) 9.79 [+ or -] 0.28 (ab) 10.22 [+ or -] 0.59 (ab) Trait H6H7 (13) H7H7 (59) PEW (%) 64.38 [+ or -] 1.74 (b) 64.29 [+ or -] 0.82@ (b) PBM (%) 9.77 [+ or -] 0.57 (ab) 9.16 [+ or -] 0.27@ (b) Trait p-value (2) pFDR PEW (%) 0.007 0.011 PBM (%) 0.048 0.034 PEW, percentage of eviscerated weight; PBM, percentage of breast muscle; pFDR, the positive false discovery rate. (1) The numbers in parentheses represent the number of diplotypes. (2) Probability of the F-test for Type III SS. Values represent least squares means plus or minus standard error. Superscripts within a trait represent a difference among breeds at p < 0.05 (Tukey's test). Bolded values represent the diplotypes with the greatest values. Underlined values represent the diplotypes with the smallest values. Note: Values represent the diplotypes with the greatest values are indicated with #. Note: Underlined values represent the diplotypes with the smallest values are indicated with #.
|Printer friendly Cite/link Email Feedback|
|Author:||Zhang, Lu; Zhu, Qing; Liu, Yiping; Gilbert, Elizabeth R.; Li, Diyan; Yin, Huadong; Wang, Yan; Yang,|
|Publication:||Asian - Australasian Journal of Animal Sciences|
|Date:||Jun 1, 2015|
|Previous Article:||Selenium in pig nutrition and reproduction: boars and semen quality--a review.|
|Next Article:||Multiple genes related to muscle identified through a joint analysis of a two-stage genome-wide association study for racing performance of 1,156...|