Expansion of the pig comparative map by expressed sequence tags (EST) mapping

Epistatic QTL pairs associated with meat quality and carcass composition traits in a porcine Duroc × Pietrain population

Background

Quantitative trait loci (QTL) analyses in pig have revealed numerous individual QTL affecting growth, carcass composition, reproduction and meat quality, indicating a complex genetic architecture. In general, statistical QTL models consider only additive and dominance effects and identification of epistatic effects in livestock is not yet widespread. The aim of this study was to identify and characterize epistatic effects between common and novel QTL regions for carcass composition and meat quality traits in pig.

Methods

Five hundred and eighty five F₂ pigs from a Duroc × Pietrain resource population were genotyped using 131 genetic markers (microsatellites and SNP) spread over the 18 pig autosomes. Phenotypic information for 26 carcass composition and meat quality traits was available for all F₂ animals. Linkage analysis was performed in a two-step procedure using a maximum likelihood approach implemented in the QxPak program.

Results

A number of interacting QTL was observed for different traits, leading to the identification of a variety of networks among chromosomal regions throughout the porcine genome. We distinguished 17 epistatic QTL pairs for carcass composition and 39 for meat quality traits. These interacting QTL pairs explained up to 8% of the phenotypic variance.

Conclusions

Our findings demonstrate the significance of epistasis in pigs. We have revealed evidence for epistatic relationships between different chromosomal regions, confirmed known QTL loci and connected regions reported in other studies. Considering interactions between loci allowed us to identify several novel QTL and trait-specific relationships of loci within and across chromosomes.

A robust linkage map of the porcine autosomes based on gene-associated SNPs

BACKGROUND

Genetic linkage maps are necessary for mapping of mendelian traits and quantitative trait loci (QTLs). To identify the actual genes, which control these traits, a map based on gene-associated single nucleotide polymorphism (SNP) markers is highly valuable. In this study, the SNPs were genotyped in a large family material comprising more than 5,000 piglets derived from 12 Duroc boars crossed with 236 Danish Landrace/Danish Large White sows. The SNPs were identified in sequence alignments of 4,600 different amplicons obtained from the 12 boars and containing coding regions of genes derived from expressed sequence tags (ESTs) and genomic shotgun sequences.

RESULTS

Linkage maps of all 18 porcine autosomes were constructed based on 456 gene-associated and six porcine EST-based SNPs. The total length of the averaged-sex whole porcine autosome was estimated to 1,711.8 cM resulting in an average SNP spacing of 3.94 cM. The female and male maps were estimated to 2,336.1 and 1,441.5 cM, respectively. The gene order was validated through comparisons to the cytogenetic and/or physical location of 203 genes, linkage to evenly spaced microsatellite markers as well as previously reported conserved synteny. A total of 330 previously unmapped genes and ESTs were mapped to the porcine autosome while ten genes were mapped to unexpected locations.

CONCLUSION

The linkage map presented here shows high accuracy in gene order. The pedigree family network as well as the large amount of meiotic events provide good reliability and make this map suitable for QTL and association studies. In addition, the linkage to the RH-map of microsatellites makes it suitable for comparison to other QTL studies.

A muscle transcriptome analysis identifies positional candidate genes for a complex trait in pig

Muscle tenderness is an important complex trait for meat quality and thus for genetic improvement through animal breeding. However, the physiological or genetic control of tenderness development in muscle is still poorly understood. In this work, using transcriptome analysis, we found a relationship between gene expression variability and tenderness. Muscle (longissimus dorsi) samples from 30 F(2) pigs were characterized by Warner-Bratzler Shear Force (WBSF) on cooked meat as a measurement of muscle tenderness. Gene expression levels were measured using microarrays for 17 muscle samples selected to represent a range of WBSF values. Using a linear regression model, we determined that samples with WBSF values above 30 N could be effectively analysed for genes exhibiting a significant association of their expression level on shear force (false discovery rate <0.05). These genes were shown to be involved in three functional networks: cell cycle, energy metabolism and muscle development. Twenty-two genes were mapped on the pig genome and 12 were found to be located in regions previously reported to contain quantitative trait loci (QTL) affecting pig meat tenderness (chromosomes 2, 6 and 13). Some genes appear therefore as positional candidate genes for QTL.

High-resolution comprehensive radiation hybrid maps of the porcine chromosomes 2p and 9p compared with the human chromosome 11

We are constructing high-resolution, chromosomal 'test' maps for the entire pig genome using a 12,000-rad WG-RH panel (IMNpRH2(12,000-rad))to provide a scaffold for the rapid assembly of the porcine genome sequence. Here we present an initial, comparative map of human chromosome (HSA) 11 with pig chromosomes (SSC) 2p and 9p. Two sets of RH mapping vectors were used to construct the RH framework (FW) maps for SSC2p and SSC9p. One set of 590 markers, including 131 microsatellites (MSs), 364 genes/ESTs, and 95 BAC end sequences (BESs) was typed on the IMNpRH2(12,000-rad) panel. A second set of 271 markers (28 MSs, 138 genes/ESTs, and 105 BESs) was typed on the IMpRH(7,000-rad) panel. The two data sets were merged into a single data-set of 655 markers of which 206 markers were typed on both panels. Two large linkage groups of 72 and 194 markers were assigned to SSC2p, and two linkage groups of 84 and 168 markers to SSC9p at a two-point LOD score of 10. A total of 126 and 114 FW markers were ordered with a likelihood ratio of 1000:1 to the SSC2p and SSC9p RH(12,000-rad) FW maps, respectively, with an accumulated map distance of 4046.5 cR(12,000 )and 1355.2 cR(7,000 )for SSC2p, and 4244.1 cR(12,000) and 1802.5 cR(7,000) for SSC9p. The kb/cR ratio in the IMNpRH2(12,000-rad) FW maps was 15.8 for SSC2p, and 15.4 for SSC9p, while the ratio in the IMpRH(7,000-rad) FW maps was 47.1 and 36.3, respectively, or an approximately 3.0-fold increase in map resolution in the IMNpRH(12,000-rad) panel over the IMpRH(7,000-rad) panel. The integrated IMNpRH(12,000-rad) andIMpRH(7,000-rad) maps as well as the genetic and BAC FPC maps provide an inclusive comparative map between SSC2p, SSC9p and HSA11 to close potential gaps between contigs prior to sequencing, and to identify regions where potential problems may arise in sequence assembly.

Genetic analysis of the cytoplasmic dynein subunit families

Cytoplasmic dyneins, the principal microtubule minus-end-directed motor proteins of the cell, are involved in many essential cellular processes. The major form of this enzyme is a complex of at least six protein subunits, and in mammals all but one of the subunits are encoded by at least two genes. Here we review current knowledge concerning the subunits, their interactions, and their functional roles as derived from biochemical and genetic analyses. We also carried out extensive database searches to look for new genes and to clarify anomalies in the databases. Our analysis documents evolutionary relationships among the dynein subunits of mammals and other model organisms, and sheds new light on the role of this diverse group of proteins, highlighting the existence of two cytoplasmic dynein complexes with distinct cellular roles.

A high-resolution comparative map between pig chromosome 17 and human chromosomes 4, 8, and 20: identification of synteny breakpoints

We report on the construction of a high-resolution comparative map of porcine chromosome 17 (SSC17) focusing on evolutionary breakpoints with human chromosomes. The comparative map shows high homology with human chromosome 20 but suggests more limited homologies with other human chromosomes. SSC17 is of particular interest in studies of chromosomal organization due to the presence of QTLs that affect meat quality and carcass composition. A total of 158 pig ESTs available in databases or developed by the Sino-Danish Pig Genome Sequencing Consortium were mapped using the INRA-University of Minnesota porcine radiation hybrid panel. The high-resolution map was further anchored by fluorescence in situ hybridization. This study confirmed the extensive conservation between SSC17 and HSA20 and enabled the gene order to be determined. The homology of the SSC17 pericentromeric region was extended to other human chromosomes (HSA4, HSA8) and the chromosomal breakpoint boundaries were accurately defined. In total 15 breakpoints were identified.

Piggy-BACing the human genome

Using the INRA-Minnesota porcine radiation hybrid panel, we have constructed a human-pig comparative map composed of 2274 loci, including 206 ESTs and 2068 BAC-end sequences, assigned to 34 linkage groups. The average spacing between comparative anchor loci is 1.15 Mb based on human genome sequence coordinates. A total of 51 conserved synteny groups that include 173 conserved segments were identified. This radiation hybrid map has the highest resolution of any porcine map to date and its integration with the porcine linkage map (reported here) will greatly facilitate the positional cloning of genes influencing complex traits of both agricultural and biomedical interest. Additionally, this map will provide a framework for anchoring contigs generated through BAC fingerprinting efforts and assist in the selection of a BAC minimal tiling path and assembly of the first sequence-ready map of the porcine genome.

Elucidation of correspondence between swine chromosome 4 and human chromosome 1 by assigning 27 genes to the ImpRH map, and development of microsatellites in the proximity of 14 genes

Loci affecting swine intramuscular fat content, backfat thickness, carcass weight, and daily weight gain were assigned to regions of swine chromosome (SSC) 4, which were shown to correspond to human chromosome (HSA) 1p22--> q25 by ZOO-FISH, bidirectional chromosome painting, as well as by the linkage map of genes. In order to select candidate genes responsible for the above traits from the human genome database, precise correspondence between SSC4 and HSA1 is a prerequisite. In the present study, 27 genes, PTGFR, GBP1, GBP2, GFI1, GCLM, ABCD3, EXTL2, KCNA3, ADORA3, KCND3, WNT2B, NRAS, SYCP1, PTGFRN, IGSF2, NOTCH2, S100A10, SHC1, SSR2, LMNA, CCT3, CD5L, PEA15, FCER1G, EAT2, DDR2, and LAMB3, located in the HSA1 region corresponding to SSC4 or possibly SSC4, were assigned to the IMpRH map. The alignment of genes from centromere to telomere in the SSC4 q arm is basically conserved in HSA1p22-->q25 with the direction from the q arm to the p arm, which is in good agreement with results from linkage mapping. In addition, the present study first demonstrated that WNT2B residing in the middle of the HSA1 region was assigned to SSC18 with a high lod score (> 5), and that at least three intrachromosomal rearrangements occurred in the region in the process of swine and human evolution. PTGFR, and LAMB3 localized at both ends of the HSA1 region were assigned to SSC6 and SSC9, respectively, which is consistent with regional correspondence reported earlier. In the course of the above analysis, microsatellite markers were developed in the proximity of eleven genes localized on SSC4, and three genes on other swine chromosomes.

Study of candidate genes for glycolytic potential of porcine skeletal muscle: identification and analysis of mutations, linkage and physical mapping and association with meat quality traits in pigs

Several genes (PRKAA2, PRKAB1, PRKAB2, PRKAG3, GAA, GYS1, PYGM, ALDOA, GPI, LDHA, PGAM2 and PKM2), chosen according to their role in the regulation of the energy balance and in the glycogen metabolism and glycolysis of the skeletal muscle, were studied. Eleven single nucleotide polymorphisms (SNPs) were identified in six of these genes (PRKAB1, GAA, PYGM, LDHA, PGAM2 and PKM2). Allele frequencies were analyzed in seven different pig breeds for these loci and for a polymorphism already described for GPI and for three polymorphic sites already reported at the PRKAG3 locus (T30N, G52S and I199V). Linkage mapping assigned PYGM and LDHA to porcine chromosome (SSC) 2, PKM2 to SSC7, GAA to SSC12, PRKAB1 to SSC14 and PGAM2 to SSC18. Physical mapping, obtained by somatic cell hybrid panel analysis, confirmed the linkage assignments of PRKAB1 and GAA and localized ALDOA, PRKAB2 and GYS1 to SSC3, SSC4 and SSC6, respectively. Pigs selected for the association study, for which several meat quality traits were measured, were first genotyped at the PRKAG3 R200Q polymorphic site (RN locus), in order to exclude carriers of the 200Q allele, and then were genotyped for all the mutations considered in this work. Significant associations (P < or = 0.001) were observed for the PRKAG3 T30N and G52S polymorphic sites with meat colour (L* at 24 h post mortem). PGAM2 and PKM2 were significantly associated (P = 0.01) with drip loss percentage and glycogen content at one hour post mortem, respectively.

cDNA cloning, genomic structure, chromosomal mapping, and expression analysis of parotid secretory protein in pig

A novel cDNA has been isolated from pig parotid glands by 3' and 5' rapid amplification of cDNA ends and designated parotid secretory protein (PSP). The open reading frame of this cDNA covers 714 bases, encoding 238 amino acids, which show 56% identity with human PSP at the level of the primary protein structure. The PSP genomic sequence comprises eight exons and seven introns, is approximately 22 kb in size, determined by sequencing, and maps to pig chromosome 17q21-q23. RT-PCR, dot blot, and Northern blot analyses demonstrated that PSP is strongly expressed in parotid glands, but is not present in heart, liver, lung, kidney, muscle, or stomach. A search for functionally significant protein motifs revealed consensus sequences for casein kinase II phosphorylation and N-myristoylation. We observed a unique amino acid sequence pattern consisting of the residues Leu-X(6)-Leu-X(6)-Leu-X(7)-Leu-X(6)-Leu-X(6)-Leu near the amino-terminal portion of the protein, which is similar to the leucine zipper.

Physical and linkage mapping of mammary-derived expressed sequence tags in cattle

This study describes the physical and linkage mapping of 42 gene-associated markers developed from mammary gland-derived expressed sequence tags to the cattle genome. Of the markers, 25 were placed on the USDA reference linkage map and 37 were positioned on the Roslin 3000-rad radiation hybrid (RH) map, with 20 assignments shared between the maps. Although no novel regions of conserved synteny between the cattle and the human genomes were identified, the coverage was extended for bovine chromosomes 3, 7, 15, and 29 compared with previously published comparative maps between human and bovine genomes. Overall, these data improve the resolution of the human-bovine comparative maps and will assist future efforts to integrate bovine RH and linkage map data.

Tropomyosin 4 expression is enhanced in dedifferentiating smooth muscle cells in vitro and during atherogenesis

Dedifferentiation of smooth muscle cells (SMC) from the contractile to the synthetic phenotype is a key event in atherosclerosis. A comparable phenotypic change from the contractile to the synthetic state is rapidly incurred when SMC are grown in culture. To identify genes that characterize the contractile and synthetic phenotypes, we performed differential display reverse transcription polymerase chain reactions on RNA from porcine arterial contractile SMC obtained directly from medial tissues and from SMC made synthetic by cell culturing. One of the differentially expressed cDNAs we identified encoded tropomyosin 4 (TM4). Whereas basal levels of TM4 existed in contractile SMC, the amount of TM4 transcripts strongly increased in synthetic SMC (33% vs. 86-106%; p < 0.005). Induction of foam cell formation had no additional enhancing effect on the expression of TM4 in cultivated SMC. We also tested whether TM4 expression was correspondingly enhanced during atherogenesis. The number of TM4-expressing SMC increased with plaque development as demonstrated by simultaneous in situ hybridization and immunohistochemistry. We compared the localization patterns of myosin heavy chain isoforms in normal arteries and lesions of increasing severity and determined that TM4 expression was relegated mainly to SMC of the synthetic phenotype in the media and intima during atherogenesis. The present study demonstrates that upregulation of TM4 mRNA is a relevant marker of dedifferentiation in vascular SMC.

Comparative mapping in the pig: localization of 214 expressed sequence tags

In total, 214 ESTs (Expressed Sequence Tags) were assigned to the porcine gene map by using somatic cell hybrid mapping, radiation hybrid mapping, and FISH. The ESTs were isolated from a porcine small intestine cDNA library on the basis of significant sequence identity with human annotated genes. In total, 390 primer pairs were designed primarily in the 3' UTR of the sequences. Overall, 58.6% of the ESTs were successfully mapped by this approach. In total, 191 of the localizations are in agreement with the human comparative map, strongly indicating that these represent true orthologous genes. The remaining 23 ESTs provide new comparative mapping data, which should be considered as preliminary until confirmed by other studies. Our mapping efforts provide a significant contribution to the porcine map as well as to the comparative map for human and pig.

Chromosomal assignments for porcine genes encoding enzymes in hepatic metabolic pathways

Increasing the number of mapped genes will facilitate (1) the identification of potential candidate genes for a trait of interest within quantitative trait loci regions and (2) comparative mapping. The metabolic activities of the liver are essential for providing fuel to peripheral organs, for regulation of amino acid, carbohydrate and lipid metabolism and for homoeostasis of vitamins, minerals and electrolytes. We aimed to identify and map genes coding for enzymes active in the liver by somatic cell genetics in order to contribute to the improvement of the porcine gene map. We mapped 28 genes of hepatic metabolic pathways including six genes whose locations could be confirmed and 22 new assignments. Localization information in human was available for all but one gene. In total 24 genes were assigned to in the expected chromosomal regions on the basis of the currently available information on the comparative human and pig map while for four genes our results suggest a new correspondence or extended regions of conservation between porcine and human chromosomes.

Measuring conservation of contiguous sets of autosomal markers on bovine and porcine genomes in relation to the map of the human genome

Based on published information, we have identified 991 genes and gene-family clusters for cattle and 764 for pigs that have orthologues in the human genome. The relative linear locations of these genes on human sequence maps were used as "rulers" to annotate bovine and porcine genomes based on a CSAM (contiguous sets of autosomal markers) approach. A CSAM is an uninterrupted set of markers in one genome (primary genome; the human genome in this study) that is syntenic in the other genome (secondary genome; the bovine and porcine genomes in this study). The analysis revealed 81 conserved syntenies and 161 CSAMs between human and bovine autosomes and 50 conserved syntenies and 95 CSAMs between human and porcine autosomes. Using the human sequence map as a reference, these 991 and 764 markers could correlate 72 and 74% of the human genome with the bovine and porcine genomes, respectively. Based on the number of contiguous markers in each CSAM, we classified these CSAMs into five size groups as follows: singletons (one marker only), small (2-4 markers), medium (5-10 markers), large (11-20 markers), and very large (> 20 markers). Several bovine and porcine chromosomes appear to be represented as di-CSAM repeats in a tandem or dispersed way on human chromosomes. The number of potential CSAMs for which no markers are currently available were estimated to be 63 between human and bovine genomes and 18 between human and porcine genomes. These results provide basic guidelines for further gene and QTL mapping of the bovine and porcine genomes, as well as insight into the evolution of mammalian genomes.

An updated linkage and comparative map of porcine chromosome 18

Swine chromosome 18 (SSC18) has the poorest marker density in the USDA-MARC porcine linkage map. In order to increase the marker density, seven genes from human chromosome 7 (HSA7) expected to map to SSC18 were selected for marker development. The genes selected were: growth hormone releasing hormone receptor (GHRHR), GLI-Kruppel family member (GLI3), leptin (LEP), capping protein muscle Z-line alpha 2 subunit (CAPZA2), beta A inhibin (INHBA), T-cell receptor beta (TCRB) and T-cell receptor gamma (TCRG). Large-insert clones (YACs, BACs and cosmids) that contained these genes, as well as two previously mapped microsatellite markers (SW1808 and SW1984), were identified and screened for microsatellites. New microsatellite markers were developed from these clones and mapped. Selected clones were also physically assigned by fluorescence in situ hybridization (FISH). Fifteen new microsatellite markers were added to the SSC18 linkage map resulting in a map of 28 markers. Six genes have been included into the genetic map improving the resolution of the SSC18 and HSA7 comparative map. Assignment of TCRG to SSC9 has identified a break in conserved synteny between SSC18 and HSA7.

Construction of a whole-genome radiation hybrid panel for high-resolution gene mapping in pigs

We have developed a panel of 152 whole-genome radiation hybrids by fusing irradiated diploid pig lymphocytes or fibroblasts with recipient hamster permanent cells. The number and size of the porcine chromosome fragments retained in each hybrid clone were checked by fluorescence in situ hybridization with a SINE probe or by primed in situ labeling (PRINS) with SINE-specific primers. A strategy based on the interspersed repetitive sequence polymerase chain reaction (IRS-PCR) was developed for selected clones to determine if the large fragments painted by the SINE probe corresponded to one pig chromosome or to different fragments of several chromosomes. This strategy was buttressed by a double PRINS approach using primers specific for alpha-satellite sequences of two different groups of swine chromosomes. Genome retention frequency was estimated for each clone by PCR with 32 markers localized on different porcine chromosomes. Of the 152 hybrids produced, 126 were selected on the basis of cytogenetic content and chromosome retention frequency to construct a radiation hybrid map of swine chromosome 8. Our initial results for this chromosome indicate that the resolution of the radiation hybrid map is 18 times higher than that obtained by linkage analysis.

Linkage mapping and comparative analysis of bovine expressed sequence tags (ESTs)

Bovine expressed sequence tags (ESTs) containing microsatellites are suitable markers for both linkage and comparative maps. We isolated clones from a bovine fetal thigh skeletal muscle cDNA library that were positive for a (CA)10 probe. Thirty individual clones were isolated and characterised by sequencing. Sequences from the 5' and 3' ends of a clone were considered as separate ESTs until a contiguous sequence was identified. A total of 47 ESTs were sequenced from the 5' and/or 3' ends and full sequence was obtained for the 30 clones. BLAST nucleotide analysis identified significant homology to known mammalian coding regions for 31 of the bovine ESTs, 30 of which also matched human ESTs or sequence-tagged sites (STS). The remaining 16 bovine ESTs represented novel transcripts. Microsatellites were isolated in 27 of the ESTs, 11 of which were developed into markers and placed on the MARC bovine linkage map. Human cytogenetic map positions were available for 20 of the 30 human EST orthologs, and a putative bovine map position for 17 of the sequences could be inferred using comparative mapping data. These results demonstrated that mapping bovine ESTs containing microsatellites is a plausible strategy to increase the density of gene markers on the bovine linkage and comparative maps.

A primary male autosomal linkage map of the horse genome

A primary male autosomal linkage map of the domestic horse (Equus caballus) has been developed by segregation analysis of 140 genetic markers within eight half-sib families. The family material comprised four Standardbred trotters and four Icelandic horses, with a total of 263 offspring. The marker set included 121 microsatellite markers, eight protein polymorphisms, five RFLPs, three blood group polymorphisms, two PCR-RFLPs, and one single strand conformation polymorphism (SSCP). One hundred markers were arranged into 25 linkage groups, 22 of which could be assigned physically to 18 different chromosomes (ECA1, ECA2, ECA3, ECA4, ECA5, ECA6, ECA7, ECA9, ECA10, ECA11, ECA13, ECA15, ECA16, ECA18, ECA19, ECA21, ECA22, and ECA30). The average distance between linked markers was 12.6 cM and the longest linkage group measured 103 cM. The total map distance contained within linkage groups was 679 cM. If the distances covered outside the ends of linkage groups and by unlinked markers were included, it was estimated that the marker set covered at least 1500 cM, that is, at least 50% of the genome. A comparison of the relationship between genetic and physical distances in anchored linkage groups gave ratios of 0.5-0.8 cM per Mb of DNA. This would suggest that the total male recombinational distance in the horse is 2000 cM; this value is lower than that suggested by chiasma counts. The present map should provide an important framework for future genome mapping in the horse.

Improvement of the porcine transcription map: localization of 33 genes, of which 24 are orthologous

From a resource of porcine ESTs, 52 transcripts were selected for regional chromosomal assignments in a somatic cell hybrid panel. Except for six ESTs, the chosen transcripts represented genes where the BLASTX database searches showed high similarity scores (>90%) with a part of the single pass 5' sequence to human, bovine, mouse, or pig entries. PCR primers for hybrid cell analysis of the ESTs were positioned in the 3'UTR of the sequences. Confident regional assignments to pig chromosomes were obtained for 33 of the 52 porcine ESTs. Comparative human mapping data were available for 27 of these. Twenty-four proved to be orthologous genes now placed on the porcine transcription map. The data presented provide further comparative data for 13 autosomes and the X chromosome.