Porcine and Bovine genome annotation and QTL dissection

Grant: AG 5438-7-3171
Award Amount: $2,010,767
PI: L.B. Schook, G.A. Rohrer
Start Date: June 1, 2007
End Date: September 30, 2014
 

Completion of the human genome sequence provided the starting point for understanding the genetic complexity of man and how genetic variation contributes to diverse phenotype and disease.  Model organisms have played an invaluable role in capturing this information.  Additional species must be sequenced to effectively extrapolate genetic information from comparative (veterinary) medicine to human medicine.

During domestication, cattle and swine underwent intense selection pressures for various phenotypes.  These selective pressures have differentiated subpopulations and produced phenotypes extremely relevant to human health research. The USDA Agricultural Research Service, the USDA's chief scientific research agency, has provided significant funding in support of sequencing efforts.

During the past five years considerable progress has been made on the research objectives of this project. The work in bovine genomics identified the causative variants for two defects (tibial hemimelia and pulmonary hypoplasia with anasarca) as well as provided a whole-genome, high resolution identity-by-descent of a father-son pair of influential dairy bulls. Porcine sequencing efforts successfully described the transcriptome of multiple tissues, characterized a family of genes important in maintaining pig health (Toll-like receptors, TLR), and co-developed the Illumina PorcineSNP60 BeadChip, which was then used to characterize pig evolution.

GOALS & OBJECTIVES

The objective of this project was to develop linkage disequilibrium (LD) haplotype maps of genes responsible for immune function and disease susceptibility in the exotic pig breeds. We will begin with Toll-like receptor genes. Toll-like receptors (TLR) 3, 7 and 8 recognize economically detrimental viruses including rotavirus and African swine fever virus. Gene specific resequencing of domestic as well as wild suids will be used to determine the influence of selective pressures on immune function during the domestication and selection of commercial swine populations. Use RNAseq to examine differential expression (DE) of genes within the muscle of animals that differ in genotype at the paternally expressed IGF2 locus. Results from previous transcriptome sequencing will be used in conjunction with deep-sequencing of the muscle transcriptome from pigs of differing genotypes and developmental time points to examine how genotype influences the transcriptome and to determine molecular phenotypes resulting from the presence or absence of the IGF2 intron 3-g30nA mutation.

Genomic DNA samples of various exotic sus species including S.barbatus, S.verrucosus, S. celebenesis and Phachoerus african have been genotyped using the Porcine SNP60Beadchip. Genotype data will be used to construct haplotype maps of regions encompassing immune regulation genes. Resequencing of specific immune genes in both commercial and exotic swine populations will identify variations responsible for altering innate immunity and host susceptibility to infectious agents. Sequence analysis of the toll-like receptor genes has begun. We will follow those with major histocompatibility complex, interleukins and immunoglobulins and will be included to provide insight to phylogenetic divergence of gene families. Though several groups have established the presence and increased prevalence of the IGF2 intron 3-G3072A mutation in the US commercial swine population, it is unclear whether increased muscling results from greater hyperplasia during fetal development and/or greater postnatal hypertrophy. Furthermore, the molecular impact of increased IGF2 expression levels remains unclear. For this objective, IGF2 mutant and wild-type animals will be produced with a common genetic background. First, a complete growth curve of body weight and length will be established from fetal development through market weight. Muscle hyperplasia, hypertrophy and fiber type will be investigated. Finally, global gene expression, throughout the life span of market animals will be established. Boars having the genotype G/A at the IGF2 locus will be mated to commercial sows which are A/A at the IGF2 locus. These matings will result in offspring with a paternal G for the IGF2 locus (lighter muscled) or a paternal A at the IGF2 locus (heavier muscled). IGF2 genotype will be determined by a real-time PCR allele discrimination assay. Both genotypes will be produced in each litter with identical sow genotypes. Key events in fetal and postnatal development will be targeted to determine changes in gene expression, muscle cell size and number and muscle fiber type, if appropriate. Muscle fibers form in two waves of fusion with primary fiber formation occurring from 35 to 50 days of gestation and secondary hyperplasia beginning between day 50 and 60 of gestation. Fiber formation is thought to be complete by day 90 of gestation. All further growth of muscle occurs through hypertrophy as muscle fiber number is set. To target developmental events, Longissimus dorsi muscle will be collected from fetuses at gestational days 50 and 90. Samples from pigs will also be collected at birth, weaning (21 days) and when harvested at market weight (25 weeks). These samples will be used to determine gene expression levels, muscle fiber hyperplasia and hypertrophy and fiber type in IGF2 mutant and wild type pigs. The W. M. Keck Center for Comparative and Functional Genomics at the University of Illinois will perform deep-sequencing to determine differential expression of mRNA between IGF2 Gpat and IGF2 Apat animals. Based on body weight and loin eye area, four animals from each genotype/time point combination which are closest to the mean will be selected. When possible, littermates will be used.

PUBLICATIONS

Snelling, W.M., R. Chiu, J.E. Schein, M. Hobbs, C.A. Abbey, D.L. Adelson, J. Aerts, G.L. Bennett, I.E. Bosdet, M. Boussaha, R. Brauning, A.R. Caetano, M.M. Costa, A.M. Crawford, B.P. Dalrymple, A. Eggen, A. Everts-Van Der Wind, S. Floriot, M. Gautier, C.A. Gill, R.D. Green, R. Holt, O. Jann, S.J. Jones, S.M. Kappes, J.W. Keele, P.J. De Jong, D.M. Larkin, H.A. Lewin, J.C. Mcewan, S. Mckay, M.A. Marra, C.A. Mathewson, L.K. Matukumalli, S.S. Moore, B. Murdoch, F.W. Nicholas, K. Osoegawa, A. Roy, H. Salih, L. Schibler, R.D. Schnabel, L. Silveri, L.C. Skow, T.P. Smith, T.S. Sonstegard, J.F. Taylor, R. Tellam, C.P. Van Tassell, J.L. Williams, J.E. Womack, N.H. Wye, G. Yang, S. Zhao. (2007). A physical map of the bovine genome. Genome Biol. 8, R165.