Genotype imputation accuracy in a F2 pig population using high density and low density SNP panels

Background: F2 resource populations have been used extensively to map QTL segregating between pig breeds. A limitation associated with the use of these resource populations for fine mapping of QTL is the reduced number of founding individuals and recombinations of founding haplotypes occurring in the population. These limitations, however, become advantageous when attempting to impute unobserved genotypes using within family segregation information. A trade-off would be to re-type F2 populations using high density SNP panels for founding individuals and low density panels (tagSNP) in F2 individuals followed by imputation. Subsequently a combined meta-analysis of several populations would provide adequate power and resolution for QTL mapping, and could be achieved at relatively low cost. Such a strategy allows the wealth of phenotypic information that has previously been obtained on experimental resource populations to be further mined for QTL identification. In this study we used experimental and simulated high density genotypes (HD-60K) from an F2 cross to estimate imputation accuracy under several genotyping scenarios. Results: Selection of tagSNP using physical distance or linkage disequilibrium information produced similar imputation accuracies. In particular, tagSNP sets averaging 1 SNP every 2.1 Mb (1,200 SNP genome-wide) yielded imputation accuracies (IA) close to 0.97. If instead of using custom panels, the commercially available 9K chip is used in the F2, IA reaches 0.99. In order to attain such high imputation accuracy the F0 and F1 generations should be genotyped at high density. Alternatively, when only the F0 is genotyped at HD, while F1 and F2 are genotyped with a 9K panel, IA drops to 0.90. Conclusions: Combining 60K and 9K panels with imputation in F2 populations is an appealing strategy to re-genotype existing populations at a fraction of the cost.Fil: Gualdron Duarte, Jose Luis. Michigan State University; Estados Unidos. Universidad de Buenos Aires. Facultad de Agronomia. Departamento de Producción Animal; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas; ArgentinaFil: Bates, Ronald O.. Michigan State University; Estados UnidosFil: Ernst, Catherine W.. Michigan State University; Estados UnidosFil: Raney, Nancy E.. Michigan State University; Estados UnidosFil: Cantet, Rodolfo Juan Carlos. Universidad de Buenos Aires. Facultad de Agronomia. Departamento de Producción Animal; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas; ArgentinaFil: Steibel, Juan P.. Michigan State University; Estados Unido

Estimation of linkage disequilibrium in four US pig breeds

<p>Abstract</p> <p>Background</p> <p>The success of marker assisted selection depends on the amount of linkage disequilibrium (LD) across the genome. To implement marker assisted selection in the swine breeding industry, information about extent and degree of LD is essential. The objective of this study is to estimate LD in four US breeds of pigs (Duroc, Hampshire, Landrace, and Yorkshire) and subsequently calculate persistence of phase among them using a 60 k SNP panel. In addition, we report LD when using only a fraction of the available markers, to estimate persistence of LD over distance.</p> <p>Results</p> <p>Average <it>r²</it>between adjacent SNP across all chromosomes was 0.36 for Landrace, 0.39 for Yorkshire, 0.44 for Hampshire and 0.46 for Duroc. For markers 1 Mb apart, <it>r²</it>ranged from 0.15 for Landrace to 0.20 for Hampshire. Reducing the marker panel to 10% of its original density, average <it>r²</it>ranged between 0.20 for Landrace to 0.25 for Duroc. We also estimated persistence of phase as a measure of prediction reliability of markers in one breed by those in another and found that markers less than 10 kb apart could be predicted with a maximal accuracy of 0.92 for Landrace with Yorkshire.</p> <p>Conclusions</p> <p>Our estimates of LD, although in good agreement with previous reports, are more comprehensive and based on a larger panel of markers. Our estimates also confirmed earlier findings reporting higher LD in pigs than in American Holstein cattle, especially at increasing marker distances (> 1 Mb). High average LD (<it>r²</it>> 0.4) between adjacent SNP found in this study is an important precursor for the implementation of marker assisted selection within a livestock species.</p> <p>Results of this study are relevant to the US purebred pig industry and critical for the design of programs of whole genome marker assisted evaluation and selection. In addition, results indicate that a more cost efficient implementation of marker assisted selection using low density panels with genotype imputation, would be feasible for these breeds.</p

Identification of Carcass and Meat Quality QTL in an F2 Duroc × Pietrain Pig Resource Population Using Different Least-Squares Analysis Models

A three-generation resource population was constructed by crossing pigs from the Duroc and Pietrain breeds. In this study, 954 F2 animals were used to identify quantitative trait loci (QTL) affecting carcass and meat quality traits. Based on results of the first scan analyzed with a line-cross (LC) model using 124 microsatellite markers and 510 F2 animals, 9 chromosomes were selected for genotyping of additional markers. Twenty additional markers were genotyped for 954 F2 animals and 20 markers used in the first scan were genotyped for 444 additional F2 animals. Three different Mendelian models using least-squares for QTL analysis were applied for the second scan: a LC model, a half-sib (HS) model, and a combined LC and HS model. Significance thresholds were determined by false discovery rate (FDR). In total, 50 QTL using the LC model, 38 QTL using the HS model, and 3 additional QTL using the combined LC and HS model were identified (q < 0.05). The LC and HS models revealed strong evidence for QTL regions on SSC6 for carcass traits (e.g., 10th-rib backfat; q < 0.0001) and on SSC15 for meat quality traits (e.g., tenderness, color, pH; q < 0.01), respectively. QTL for pH (SSC3), dressing percent (SSC7), marbling score and moisture percent (SSC12), CIE a* (SSC16), and carcass length and spareribs weight (SSC18) were also significant (q < 0.01). Additional marker and animal genotypes increased the statistical power for QTL detection, and applying different analysis models allowed confirmation of QTL and detection of new QTL

Application of alternative models to identify QTL for growth traits in an F2 Duroc x Pietrain pig resource population

<p>Abstract</p> <p>Background</p> <p>A variety of analysis approaches have been applied to detect quantitative trait loci (QTL) in experimental populations. The initial genome scan of our Duroc x Pietrain F₂resource population included 510 F₂animals genotyped with 124 microsatellite markers and analyzed using a line-cross model. For the second scan, 20 additional markers on 9 chromosomes were genotyped for 954 F₂animals and 20 markers used in the first scan were genotyped for 444 additional F₂animals. Three least-squares Mendelian models for QTL analysis were applied for the second scan: a line-cross model, a half-sib model, and a combined line-cross and half-sib model.</p> <p>Results</p> <p>In total, 26 QTL using the line-cross model, 12 QTL using the half-sib model and 3 additional QTL using the combined line-cross and half-sib model were detected for growth traits with a 5% false discovery rate (FDR) significance level. In the line-cross analysis, highly significant QTL for fat deposition at 10-, 13-, 16-, 19-, and 22-wk of age were detected on SSC6. In the half-sib analysis, a QTL for loin muscle area at 19-wk of age was detected on SSC7 and QTL for 10th-rib backfat at 19- and 22-wk of age were detected on SSC15.</p> <p>Conclusions</p> <p>Additional markers and animals contributed to reduce the confidence intervals and increase the test statistics for QTL detection. Different models allowed detection of new QTL which indicated differing frequencies for alternative alleles in parental breeds.</p

Transcriptional profiling identifies differentially expressed genes in developing turkey skeletal muscle

<p>Abstract</p> <p>Background</p> <p>Skeletal muscle growth and development from embryo to adult consists of a series of carefully regulated changes in gene expression. Understanding these developmental changes in agriculturally important species is essential to the production of high quality meat products. For example, consumer demand for lean, inexpensive meat products has driven the turkey industry to unprecedented production through intensive genetic selection. However, achievements of increased body weight and muscle mass have been countered by an increased incidence of myopathies and meat quality defects. In a previous study, we developed and validated a turkey skeletal muscle-specific microarray as a tool for functional genomics studies. The goals of the current study were to utilize this microarray to elucidate functional pathways of genes responsible for key events in turkey skeletal muscle development and to compare differences in gene expression between two genetic lines of turkeys. To achieve these goals, skeletal muscle samples were collected at three critical stages in muscle development: 18d embryo (hyperplasia), 1d post-hatch (shift from myoblast-mediated growth to satellite cell-modulated growth by hypertrophy), and 16wk (market age) from two genetic lines: a randombred control line (RBC2) maintained without selection pressure, and a line (F) selected from the RBC2 line for increased 16wk body weight. Array hybridizations were performed in two experiments: Experiment 1 directly compared the developmental stages within genetic line, while Experiment 2 directly compared the two lines within each developmental stage.</p> <p>Results</p> <p>A total of 3474 genes were differentially expressed (false discovery rate; FDR < 0.001) by overall effect of development, while 16 genes were differentially expressed (FDR < 0.10) by overall effect of genetic line. Ingenuity Pathways Analysis was used to group annotated genes into networks, functions, and canonical pathways. The expression of 28 genes involved in extracellular matrix regulation, cell death/apoptosis, and calcium signaling/muscle function, as well as genes with miscellaneous function was confirmed by qPCR.</p> <p>Conclusions</p> <p>The current study identified gene pathways and uncovered novel genes important in turkey muscle growth and development. Future experiments will focus further on several of these candidate genes and the expression and mechanism of action of their protein products.</p

Rapid screening for phenotype-genotype associations by linear transformations of genomic evaluations

Background: Currently, association studies are analysed using statistical mixed models, with marker effects estimated by a linear transformation of genomic breeding values. The variances of marker effects are needed when performing the tests of association. However, approaches used to estimate the parameters rely on a prior variance or on a constant estimate of the additive variance. Alternatively, we propose a standardized test of association using the variance of each marker effect, which generally differ among each other. Random breeding values from a mixed model including fixed effects and a genomic covariance matrix are linearly transformed to estimate the marker effects. Results: The standardized test was neither conservative nor liberal with respect to type I error rate (false-positives), compared to a similar test using Predictor Error Variance, a method that was too conservative. Furthermore, genomic predictions are solved efficiently by the procedure, and the p-values are virtually identical to those calculated from tests for one marker effect at a time. Moreover, the standardized test reduces computing time and memory requirements. The following steps are used to locate genome segments displaying strong association. The marker with the highest − log(p-value) in each chromosome is selected, and the segment is expanded one Mb upstream and one Mb downstream of the marker. A genomic matrix is calculated using the information from those markers only, which is used as the variance-covariance of the segment effects in a model that also includes fixed effects and random genomic breeding values. The likelihood ratio is then calculated to test for the effect in every chromosome against a reduced model with fixed effects and genomic breeding values. In a case study with pigs, a significant segment from chromosome 6 explained 11% of total genetic variance. Conclusions: The standardized test of marker effects using their own variance helps in detecting specific genomic regions involved in the additive variance, and in reducing false positives. Moreover, genome scanning of candidate segments can be used in meta-analyses of genome-wide association studies, as it enables the detection of specific genome regions that affect an economically relevant trait when using multiple populations.Fil: Gualdron Duarte, Jose Luis. Universidad de Buenos Aires. Facultad de Agronomia. Departamento de Producción Animal; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas; ArgentinaFil: Cantet, Rodolfo Juan Carlos. Universidad de Buenos Aires. Facultad de Agronomia. Departamento de Producción Animal; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas; ArgentinaFil: Bates, Ronald O.. Michigan State University; Estados UnidosFil: Ernst, Catherine W.. Michigan State University; Estados UnidosFil: Raney, Nancy E.. Michigan State University; Estados UnidosFil: Steibel, Juan P.. Michigan State University; Estados Unido

Gene expression profiling in hepatic tissue of newly weaned pigs fed pharmacological zinc and phytase supplemented diets

<p>Abstract</p> <p>Background</p> <p>Zinc (Zn) is an essential trace element. However, Zn bioavailability from commonly consumed plants may be reduced due to phytic acid. Zn supplementation has been used to treat diarrheal disease in children, and in the U.S. swine industry at pharmacological levels to promote growth and fecal consistency, but underlying mechanisms explaining these beneficial effects remain unknown. Moreover, adding supplemental phytase improves Zn bioavailability. Thus, we hypothesized that benefits of pharmacological Zn supplementation result from changes in gene expression that could be further affected by supplemental phytase. The goal of this study was to investigate the effects of feeding newly weaned pigs dietary Zn (150, 1,000, or 2,000 mg Zn/kg) as Zn oxide with or without phytase [500 phytase units (FTU)/kg] for 14 d on hepatic gene expression. Liver RNA from pigs fed 150, 1,000, or 2,000 mg Zn/kg, or 1,000 mg Zn/kg with phytase (n = 4 per treatment) was reverse transcribed and examined using the differential display reverse transcription polymerase chain reaction technique. Liver RNA from pigs fed 150 or 2,000 mg Zn/kg (n = 4 per treatment) was also evaluated using a 70-mer oligonucleotide microarray.</p> <p>Results</p> <p>Expressed sequence tags for 61 putatively differentially expressed transcripts were cloned and sequenced. In addition, interrogation of a 13,297 element oligonucleotide microarray revealed 650 annotated transcripts (FDR ≤ 0.05) affected by pharmacological Zn supplementation. Seven transcripts exhibiting differential expression in pigs fed pharmacological Zn with sequence similarities to genes encoding <it>GLO1</it>, <it>PRDX4</it>, <it>ACY1</it>, <it>ORM1</it>, <it>CPB2</it>, <it>GSTM4</it>, and <it>HSP70.2 </it>were selected for confirmation. Relative hepatic <it>GLO1 </it>(<it>P </it>< 0.0007), <it>PRDX4 </it>(<it>P </it>< 0.009) and <it>ACY1 </it>(<it>P </it>< 0.01) mRNA abundances were confirmed to be greater in pigs fed 1,000 (n = 8) and 2,000 (n = 8) mg Zn/kg than in pigs fed 150 (n = 7) mg Zn/kg. Relative hepatic <it>HSP70.2 </it>(P < 0.002) mRNA abundance was confirmed to be lower in pigs fed 2,000 mg Zn/kg than in pigs fed 150 or 1,000 mg Zn/kg.</p> <p>Conclusion</p> <p>Results suggest that feeding pharmacological Zn (1,000 or 2,000 mg Zn/kg) affects genes involved in reducing oxidative stress and in amino acid metabolism, which are essential for cell detoxification and proper cell function.</p

Refining genomewide association for growth and fat deposition traits in an F₂ pig population

The identification of genomic regions that affect additive genetic variation and contain genes involved in controlling growth and fat deposition has enormous impact in the farm animal industry (e.g., carcass merit and meat quality). Therefore, a genomewide association study was implemented in an F₂ pig population using a 60,000 SNP marker panel for traits related to growth and fat deposition. Estimated genomic EBV were linearly transformed to calculate SNP effects and to identify genomic positions possibly associated with the genetic variability of each trait. Genomic segments were then defined considering the markers included in a region 1 Mb up- and downstream from the SNP with the smallest -value and a false discovery rate < 0.05 for each trait. The significance for each 2-Mb segment was tested using the Bonferroni correction. Significant SNP were detected on SSC2, SSC3, SSC5, and SSC6, but 2-Mb segment significant effects were observed on SSC3 for weight at birth (wt_birth) and on SSC6 for 10th-rib backfat and last-rib backfat measured by ultrasound at different ages. Furthermore, a 6-Mb segment on SSC6 was also considered because the 2-Mb segments for 10 different fat deposition traits were overlapped. Although the segment effects for each trait remain significant, the proportion of additive variance explained by this larger segment was slightly smaller in some traits. In general, the results confirm the presence of genetic variability for wt_birth on SSC3 (18.0-20.2 Mb) and for fat deposition traits on SSC6 (133.8-136.0 Mb). Within these regions, fibrosin (FBRS) and myosin light chain, phosphorylatable, fast skeletal muscle (MYLPF) genes could be considered as candidates for the wt_birth signal on SSC3, and the SERPINE1 mRNAbinding protein 1 gene (SERBP1) may be a candidate for the fat deposition trait signals on SSC6.Facultad de Ciencias VeterinariasInstituto de Genética Veterinari