Investigation of a QTL Region for Loin Eye Area and Fatness on Pig Chromosome 1
Total Page:16
File Type:pdf, Size:1020Kb
Investigation of a QTL region for loin eye area and fatness on pig Chromosome 1 Laura Grapes, Max F. Rothschild Department of Animal Science and Center for Integrated Animal Genomics, 2255 Kildee Hall, Iowa State University, Ames, Iowa 50011, USA Received: 22 December 2005 / Accepted: 14 March 2006 quantitative trait loci (QTL) that impact traits of Abstract economic importance. However, very few of the Previously, quantitative trait loci (QTL) for tenth-rib genes and/or mutations underlying the observed backfat (TENTHRIB) and loin eye area (LEA) were QTL have been identified (Kim et al. 2000; Van Laere identified on pig Chromosome 1 (SSC 1) near micro- et al. 2003). In fact, few genes have been identified as satellite S0008 from a three-generation Berkshire · significant positional candidates for observed QTL Yorkshire cross (BY). This work attempted to refine effects. Ciobanu et al. (2004) showed that mutations these QTL positions and identify genes associated in the calpastatin gene, which is located in a region with these QTL. Genotypes of BY (n = 555) were suggestive for average instron force, tenderness, determined by PCR-RFLP or PCR tests for 13 poly- juiciness, and chewiness QTL, were significantly morphisms identified in BY F0 individuals for candi- associated with these traits. Clearly, further work is date genes, BAC end sequences, and genomic clones. necessary to identify the genes underlying the Using least-squares regression interval mapping, the results from the numerous porcine genome scans. LEA QTL was estimated at S0008; the TENTHRIB Malek et al. (2001a) reported genome-wise sig- QTL position was shifted approximately 1 cM nificant QTL for loin eye area (LEA) and tenth-rib downstream from S0008. Of the genes/sequences backfat (TENTHRIB) on pig Chromosome 1 (SSC 1) mapped in the QTL region, CL349415 was signifi- near microsatellite S0008. These QTL were identi- cantly associated with TENTHRIB (p = 0.02) and fied in a three-generation Berkshire · Yorkshire solute carrier family 2, member 12 (SLC2A12) was breed cross (BY). Alleles originating from the Berk- significantly associated with LEA (p = 0.02). These shire breed interestingly were shown to be cryptic results suggest that the gene(s) responsible for the and increase LEA and decrease TENTHRIB; these LEA and TENTHRIB QTL effects are tightly linked to QTL explained 4.21% and 4.78% of the phenotypic S0008 or that the high informativeness of S0008 rel- variance for LEA and TENTHRIB, respectively, in ative to surrounding markers is influencing the QTL BY F2 individuals. In the same chromosomal region, position estimates. In addition, janus kinase 2 (JAK2) a chromosome-wise significant QTL for average was mapped to a suggestive LEA QTL region and backfat (AVBFAT) and suggestive QTL for marbling showed association with LEA (p = 0.009), fatness, score (MARB), total lipid percentage (TOTLIPPR), color, and pH traits in BY. lumbar fat (LUMBAR), and cholesterol (CHOLES) were identified in BY (Malek et al. 2001a,b). Thom- sen et al. (2004) performed further QTL analyses in the same BY population after adding 33 microsatel- lite markers to the overall genome map, 3 of which Introduction were located on SSC 1. The LEA and TENTHRIB Several whole-genome scans have been conducted in QTL identified by Malek et al. (2001a) were con- a variety of pig breed crosses (Malek et al. 2001a,b; firmed, but the AVBFAT QTL dropped just below Milan et al. 2002; Ovilo et al. 2002; Varona et al. significance (Thomsen et al. 2004). These results 2002) to identify chromosomal regions containing indicate that there are gene(s) located near S0008 that have a large impact on LEA and fatness traits in Correspondence to: Max F. Rothschild; E-mail: mfrothsc@ Berkshire and/or Yorkshire pigs. Hence, the objec- iastate.edu tive of this work was to fine map these QTL regions DOI: 10.1007/s00335-005-0188-7 Volume 17, 657 668 (2006) Ó Springer Science+Business Media, Inc. 2006 657 À 658 L. Grapes and M.F. Rothschild: FINE MAPPING QTL FOR LOIN EYE AREA AND FATNESS ON SSC 1 as an initial step toward identifying the underlying 2001) and the gene encoding it (RRAGD) maps to HSA genes. 6q15. The nuclear receptor 2 E1 gene (NR2E1) maps to HSA 6q21 and was selected as a representative gene from this human region to determine the comparative Materials and methods map position in pig. Connective tissue growth factor Resource population. A three-generation resource functions in insulin-like growth factor binding population was created and managed as described by and stimulation of cell proliferation (reviewed in Malek et al. (2001a). Briefly, two purebred Berkshire Brigstock 1999) and the gene (CTGF) maps to HSA 6q grandsires were crossed with nine purebred Yorkshire 23.2. Solute carrier family 2, member 12 gene granddams, producing nine litters from which 8 F1 (SLC2A12) maps to HSA 6q 23.2, is expressed in adi- sires and 26 F1 dams were selected. These F1 individ- pose tissue, and has glucose transport properties (re- uals were crossed to produce 525 F2 individuals. viewed in Joost and Thorens 2001). Ectonucleotide pyrophosphatase 1 gene (ENPP1) also maps to HSA Phenotypic data collection. Phenotypic values 6q23.2, about 140 kb upstream of CTGF, and muta- for loin eye area, tenth-rib backfat, average backfat, tions in the gene have been significantly associated lumbar backfat, last rib backfat, marbling score, total with insulin resistance and obesity (Pizzuti et al. lipid percent, and cholesterol were collected on pigs 1999; Meyre et al. 2005). Aldehyde dehydrogenase 8 as described in Malek et al. (2001a,b). family, member A1 gene (ALDH8A1) maps to HSA 6q23.3 and plays a role in the pathway of 9-cis-retinoic Genetic marker development. Microsatellite acid biosynthesis (Lin and Napoli 2000). In addition, marker genotypes for SSC 1 described in Malek et al. the janus kinase 2 gene (JAK2) was selected for study (2001a) and Thomsen et al. (2004) were available. To because of its role in muscle cell differentiation and increase marker density in the region significant for proliferation and its known location on SSC 1 (Wang LEA and TENTHRIB QTL, several candidate genes, et al. 2004). BAC end sequences, and genomic clones were se- After candidate genes CTGF and ENPP1 were lected for polymorphism discovery and genotyping. mapped in the BY flanking S0008, pig BAC end se- quence bE1I10T7 was identified as homologous to Candidate gene, BAC end sequence, and clone human sequence between CTGF and ENPP1 using the selection. Candidate genes were selected based upon Porcine BES Search tool available from The Wellcome comparative human pig chromosomal information Trust Sanger Institute (http://www.sanger.ac.uk/ and known functionalÀ information. The microsatel- cgi-bin/Projects/S_scrofa/BESsearch.cgi) and selected lite marker S0008, nearest to the most significant for marker development. Two porcine genomic positions for LEA and TENTHRIB QTL, had been clones, CL349415 and CL353852, were also selected physically mapped to SSC 1p22-23 (Robic et al. for marker development based on their comparative 1996). This region is syntenic to the q arm of human human map positions (Rogatcheva et al. 2006). For Chromosome 6 (HSA 6; Fronicke et al. 1996). How- NR2E1, no annotated pig sequence was available. In- ever, the synteny between the p arm of SSC 1 and the stead, one BAC end sequence, 199F3SP6, identified as q arm of HSA 6 is not well defined, and gene order in homologous to the human gene position for NR2E1 this region of SSC 1p has been inverted compared to using the Porcine BES Search tool, was selected. HSA 6q. Therefore, candidate genes from HSA 6q12 q23 were selected for marker development and Polymorphism identification. For all genes and subsequentÀ mapping in the pig to determine overall sequences, PCR primers were designed from avail- gene order on SSC 1p and increase marker density in able human or pig sequence using the default set- the QTL regions. tings of Primer3 web-based software (http:// The Ph.D. finger protein 3 gene (PHF3) maps to frodo.wi.mit.edu/cgi-bin/primer3/primer3_www.cgi) HSA 6q12 and was chosen to determine whether to amplify genomic DNA from the 11 F0 animals of genes from this human region may map to the pig the BY population. Accession numbers of human QTL region of interest. The potassium channel gene and pig sequences from which primers were de- (KCNQ5) maps to HSA 6q13 and is involved in volt- signed, the sequences of those primers, and the age-gated potassium transfer (Lerche et al. 2000). The resulting fragment sizes are given in Table 1. Se- gene encoding phosphoglucomutase 3 (PGM3) maps quences from the F0 individuals were compared for to HSA 6q14.2 and is involved in carbohydrate each gene/clone/BAC end sequence using Sequen- metabolism (Hopkinson and Harris 1968). Ras-related cher v3.0 software (Gene Codes Corporation, Ann GTP binding D is a G protein that acts in numerous Arbor, MI) to identify polymorphisms. Most poly- cell processes and signaling pathways (Sekiguchi et al. morphisms discovered were single nucleotide poly- L. Grapes and M.F. Rothschild: FINE MAPPING QTL FOR LOIN EYE AREA AND FATNESS ON SSC 1 659 Table 1. Pig/human sequences used to design initial PCR primers, primer sequences, and resulting fragment sizes for genes and sequences examined on SSC 1 Gene/Sequence (Accession No. / TIGR gene index) Primer sequences Fragment size (bp) PHF3 exon 16 (NM_015153) F: 5¢ GCCTCCGAATGTCTTTAACC 3¢ 237 R: 5¢ GATTTCCATTTCCAGGACCTC 3¢ KCNQ5 exon 14 (TC145264) F: 5¢ AGCTCCGCACAGAACAGC