Swine Genome Sequencing Consortium (SGSC): a strategic roadmap for sequencing the pig genome

Design of a high density SNP genotyping assay in the pig using SNPs identified and characterized by next generation sequencing technology

BACKGROUND

The dissection of complex traits of economic importance to the pig industry requires the availability of a significant number of genetic markers, such as single nucleotide polymorphisms (SNPs). This study was conducted to discover several hundreds of thousands of porcine SNPs using next generation sequencing technologies and use these SNPs, as well as others from different public sources, to design a high-density SNP genotyping assay.

METHODOLOGY/PRINCIPAL FINDINGS

A total of 19 reduced representation libraries derived from four swine breeds (Duroc, Landrace, Large White, Pietrain) and a Wild Boar population and three restriction enzymes (AluI, HaeIII and MspI) were sequenced using Illumina's Genome Analyzer (GA). The SNP discovery effort resulted in the de novo identification of over 372K SNPs. More than 549K SNPs were used to design the Illumina Porcine 60K+SNP iSelect Beadchip, now commercially available as the PorcineSNP60. A total of 64,232 SNPs were included on the Beadchip. Results from genotyping the 158 individuals used for sequencing showed a high overall SNP call rate (97.5%). Of the 62,621 loci that could be reliably scored, 58,994 were polymorphic yielding a SNP conversion success rate of 94%. The average minor allele frequency (MAF) for all scorable SNPs was 0.274.

CONCLUSIONS/SIGNIFICANCE

Overall, the results of this study indicate the utility of using next generation sequencing technologies to identify large numbers of reliable SNPs. In addition, the validation of the PorcineSNP60 Beadchip demonstrated that the assay is an excellent tool that will likely be used in a variety of future studies in pigs.

Comparative genomics of Toll-like receptor signalling in five species

BACKGROUND

Over the last decade, several studies have identified quantitative trait loci (QTL) affecting variation of immune related traits in mammals. Recent studies in humans and mice suggest that part of this variation may be caused by polymorphisms in genes involved in Toll-like receptor (TLR) signalling. In this project, we used a comparative approach to investigate the importance of TLR-related genes in comparison with other immunologically relevant genes for resistance traits in five species by associating their genomic location with previously published immune-related QTL regions.

RESULTS

We report the genomic localisation of TLR1-10 and ten associated signalling molecules in sheep and pig using in-silico and/or radiation hybrid (RH) mapping techniques and compare their positions with their annotated homologues in the human, cattle and mouse whole genome sequences. We also report medium-density RH maps for porcine chromosomes 8 and 13. A comparative analysis of the positions of previously published relevant QTLs allowed the identification of homologous regions that are associated with similar health traits in several species and which contain TLR related and other immunologically relevant genes. Additional evidence was gathered by examining relevant gene expression and association studies.

CONCLUSION

This comparative genomic approach identified eight genes as potentially causative genes for variations of health related traits. These include susceptibility to clinical mastitis in dairy cattle, general disease resistance in sheep, cattle, humans and mice, and tolerance to protozoan infection in cattle and mice. Four TLR-related genes (TLR1, 6, MyD88, IRF3) appear to be the most likely candidate genes underlying QTL regions which control the resistance to the same or similar pathogens in several species. Further studies are required to investigate the potential role of polymorphisms within these genes.

Refined localization of the Escherichia coli F4ab/F4ac receptor locus on pig chromosome 13

Diarrhoea in newborn and weaned pigs caused by enterotoxigenic Escherichia coli (ETEC) expressing F4 fimbriae leads to considerable losses in pig production. In this study, we refined the mapping of the receptor locus for ETEC F4ab/F4ac adhesion (F4bcR) by joint analysis of Nordic and Swiss data. A total of 236 pigs from a Nordic experimental herd, 331 pigs from a Swiss experimental herd and 143 pigs from the Swiss performing station were used for linkage analysis. Genotyping data of six known microsatellite markers, two newly developed markers (MUC4gt and HSA125gt) and an intronic SNP in MUC4 (MUC4-8227) were used to create the linkage map. The region for F4bcR was refined to the interval SW207-S0075 on pig chromosome 13. The most probable position of F4bcR was in the SW207-MUC4 region. The order of six markers was supported by physical mapping on the BAC fingerprint contig from the Wellcome Trust Sanger Institute. Thus, the region for F4bcR could be reduced from 26 to 14 Mb.

Porcine Toll-like receptors: the front line of pathogen monitoring and possible implications for disease resistance

Toll-like receptors (TLRs) are the most famous pattern-recognition receptors (PRRs); they monitor pathogen-associated molecular patterns and play a critical role in activation of the immune system against infection. TLR gene mutations may affect the gene products in terms of their ligand-binding ability or their signal transduction ability after ligand binding; such changes have a great influence on pathogen monitoring and disease resistance. Thirteen mammalian TLRs have been identified, and genes corresponding to all 10 TLR genes identified in humans have been fully cloned in pigs. Porcine TLR gene coding sequences possess a large number of nonsynonymous single nucleotide polymorphisms (SNPs). They are concentrated in ectodomains, and may increase the variability of pathogen recognition in pig populations. We summarize the current knowledge of TLR molecules in mammals and livestock (particularly pigs) and speculate on the relationship between SNPs in porcine TLRs and their application to vaccine design and disease-resistance breeding.

Genomic sequence encoding diversity segments of the pig TCR delta chain gene demonstrates productivity of highly diversified repertoire

To better understand the function and diversity of gammadelta T cells, we determined the genomic sequence encoding diversity (D) segments of the porcine TCR delta chain and its upstream regions, because pigs and other artiodactyls have relatively high proportions of gammadelta T cells. The revealed sequence contained 28 variable (V) alpha/delta segments, including 4 TRDV1 and at least 6 Ddelta segments, a much higher number than in humans and mice. All 6 of the Ddelta segments that had canonical recombination signal sequences were functionally utilized in expressed TCR delta chain genes. The multiplicity of Ddelta segments enabled the use of more than 3 Ddelta segments in a single functional TCR delta chain. The increased number of TCR delta segments was acquired by the duplication of the germline sequence, which occurred after the divergence of artiodactyls from primates and rodents. These data demonstrate that the pig is able to generate a highly diversified repertoire of TCR delta chain molecules.

Mining for single nucleotide polymorphisms in pig genome sequence data

Background

Single nucleotide polymorphisms (SNPs) are ideal genetic markers due to their high abundance and the highly automated way in which SNPs are detected and SNP assays are performed. The number of SNPs identified in the pig thus far is still limited.

Results

A total of 4.8 million whole genome shotgun sequences obtained from the NCBI trace-repository with center name "SDJVP", and project name "Sino-Danish Pig Genome Project" were analysed for the presence of SNPs. Available BAC and BAC-end sequences and their naming and mapping information, all obtained from SangerInstitute FTP site, served as a rough assembly of a reference genome. In 1.2 Gb of pig genome sequence, we identified 98,151 SNPs in which one of the sequences in the alignment represented the polymorphism and 6,374 SNPs in which two sequences represent an identical polymorphism. To benchmark the SNP identification method, 163 SNPs, in which the polymorphism was represented twice in the sequence alignment, were selected and tested on a panel of three purebred boar lines and wild boar. Of these 163 in silico identified SNPs, 134 were shown to be polymorphic in our animal panel.

Conclusion

This SNP identification method, which mines for SNPs in publicly available porcine shotgun sequences repositories, provides thousands of high quality SNPs. Benchmarking in an animal panel showed that more than 80% of the predicted SNPs represented true genetic variation.

An animal model to evaluate the function and regulation of the adaptively evolving stress protein SEP53 in oesophageal bile damage responses

Squamous epithelium in mammals has evolved an atypical stress response involving down-regulation of the classic HSP70 protein and induction of sets of proteins including one named SEP53. This atypical stress response might be due to the unusual environmental pressures placed on squamous tissue. In fact, SEP53 plays a role as an anti-apoptotic factor in response to DNA damage induced by deoxycholic acid stresses implicated in oesophageal reflux disease. SEP53 also has a genetic signature characteristic of an adaptively and rapidly evolving gene, and this observation has been used to imply a role for SEP53 in immunity. Physiological models of squamous tissue are required to further define the regulation and function of SEP53. We examined whether porcine squamous epithelium would be a good model to study SEP53, since this animal suffers from a bile-reflux disease in squamous oesophageal tissue. We have (1) cloned and sequenced the porcine SEP53 locus from porcine bacterial artificial chromosome genomic DNA, (2) confirmed the strikingly divergent nature of the C-terminal portion of the SEP53 gene amongst mammals, (3) discovered that a function of the conserved N-terminal domain of the gene is to maintain cytoplasmic localisation, and (4) examined SEP53 expression in normal and diseased porcine pars oesophagea. SEP53 expression in porcine tissue was relatively confined to gastric squamous epithelium, consistent with its expression in normal human squamous epithelium. Immunohistochemical staining for SEP53 protein in normal and damaged pars oesophagea demonstrated significant stabilisation of SEP53 protein in the injured tissue. These results suggest that porcine squamous epithelium would be a robust physiological model to examine the evolution and function of the SEP53 stress pathway in modulating stress-induced responses in squamous tissue.

A high utility integrated map of the pig genome

BACKGROUND

The domestic pig is being increasingly exploited as a system for modeling human disease. It also has substantial economic importance for meat-based protein production. Physical clone maps have underpinned large-scale genomic sequencing and enabled focused cloning efforts for many genomes. Comparative genetic maps indicate that there is more structural similarity between pig and human than, for example, mouse and human, and we have used this close relationship between human and pig as a way of facilitating map construction.

RESULTS

Here we report the construction of the most highly continuous bacterial artificial chromosome (BAC) map of any mammalian genome, for the pig (Sus scrofa domestica) genome. The map provides a template for the generation and assembly of high-quality anchored sequence across the genome. The physical map integrates previous landmark maps with restriction fingerprints and BAC end sequences from over 260,000 BACs derived from 4 BAC libraries and takes advantage of alignments to the human genome to improve the continuity and local ordering of the clone contigs. We estimate that over 98% of the euchromatin of the 18 pig autosomes and the X chromosome along with localized coverage on Y is represented in 172 contigs, with chromosome 13 (218 Mb) represented by a single contig. The map is accessible through pre-Ensembl, where links to marker and sequence data can be found.

CONCLUSION

The map will enable immediate electronic positional cloning of genes, benefiting the pig research community and further facilitating use of the pig as an alternative animal model for human disease. The clone map and BAC end sequence data can also help to support the assembly of maps and genome sequences of other artiodactyls.

Impact of genomics on animal agriculture and opportunities for animal health

Livestock production is an important source of animal protein worldwide. In the developed world meat consumption will remain steady but demand is forecast to grow enormously in developing countries. The use of genomics will speed genetic improvement and increase levels of production quickly in the developed world but might face problems in the developing world, including scientific, economic and political challenges. Considerable increases in public and private research funding will be required to develop and utilize novel tools and collections of detailed trait information on appropriate animals. The development of policies protecting the environment and managing all genetic resources will also be needed. Advances in livestock genomics have major implications for increasing food output as well as improving human health.

Genetic resources, genome mapping and evolutionary genomics of the pig (Sus scrofa)

The pig, a representative of the artiodactyla clade, is one of the first animals domesticated, and has become an important agriculture animal as one of the major human nutritional sources of animal based protein. The pig is also a valuable biomedical model organism for human health. The pig's importance to human health and nutrition is reflected in the decision to sequence its genome (3X). As an animal species with its wild ancestors present in the world, the pig provides a unique opportunity for tracing mammalian evolutionary history and defining signatures of selection resulting from both domestication and natural selection. Completion of the pig genome sequencing project will have significant impacts on both agriculture and human health. Following the pig whole genome sequence drafts, along with large-scale polymorphism data, it will be possible to conduct genome sweeps using association mapping, and identify signatures of selection. Here, we provide a description of the pig genome sequencing project and perspectives on utilizing genomic technologies to exploit pig genome evolution and the molecular basis for phenotypic traits for improving pig production and health.

Detection of novel quantitative trait loci for cutaneous melanoma by genome-wide scan in the MeLiM swine model

Human cutaneous melanoma is a complex trait inherited in about 10% of cases. Although 2 high-risk genes, CDKN2A and CDK4, and 1 low risk gene, MC1R, have been identified, susceptibility genes remain to be discovered. Here, we attempted to determine new genomic regions linked to melanoma using the pig MeLiM strain, which develops hereditary cutaneous melanomas. We applied quantitative trait loci (QTL) mapping method to a significant genome-wide scan performed on 331 backcross pigs derived from this strain. QTLs were detected at chromosome-wide level for a melanoma synthetic trait corresponding to the development of melanoma. The peak positions on Sus scrofa chromosomes (SSC) were at 49.4 and 88.0 cM (SSC1), 56.0 cM (SSC13), 86.5 cM (SSC15) and 39.8 cM (SSC17), and, on SSC2, at 16.9 cM, in families derived from F1 males only (p < 0.05, except for SSC13, p < 0.01). Analysis of 7 precise specific traits revealed highly significant QTLs on SSC10 (ulceration), on SSC12 (presence of melanoma at birth), on SSC13 (lesion type), and on SSC16 and SSC17 (number of aggressive melanomas) at the respective positions 42.0, 95.6, 81.0, 45.3 and 44.8 cM (p < 0.001 and p < 0.05 respectively at the chromosome- and genome-wide levels). We also showed that MeLiM MC1R*2 allele, which determines black coat colour in pigs, predisposes significantly to melanoma. Interactions were observed between MC1R and markers located on SSC1 (p < 0.05). Taken together, these results indicate that MeLiM swine is a model for human multigenic diseases. Comparative mapping revealed human regions of interest to search for new melanoma susceptibility candidates.

PEDE (Pig EST Data Explorer) has been expanded into Pig Expression Data Explorer, including 10 147 porcine full-length cDNA sequences

We formerly released the porcine expressed sequence tag (EST) database Pig EST Data Explorer (PEDE; ), which comprised 68 076 high-quality ESTs obtained by using full-length-enriched cDNA libraries derived from seven tissues. We have added eight tissues and cell types to the EST analysis and have integrated 94 555 additional high-quality ESTs into the database. We also fully sequenced the inserts of 10 147 of the cDNA clones that had undergone EST analysis; the sequences and annotation of the cDNA clones were stored in the database. Further, we constructed an interface that can be used to perform various searches in the database. The PEDE database is the primary resource of expressed pig genes that are supported by full-length cDNA sequences. This resource not only enables us to pick cDNA clones of interest for a particular analysis, but it also confirms and thus contributes to the sequencing integrity of the pig genome, which is now being compiled by an international consortium (). PEDE has therefore evolved into what we now call ‘Pig Expression Data Explorer’.

Evolutionary breakpoints through a high-resolution comparative map between porcine chromosomes 2 and 16 and human chromosomes

This study reports a high-resolution comparative map between human chromosomes and porcine chromosomes 2 (SSC2) and 16 (SSC16), pointing out new homologies and evolutionary breakpoints. SSC2 is of particular interest because of the presence of several important QTLs. Among 226 porcine ESTs selected according to their expected localization, 151 were RH mapped and ordered on SSC2. This study confirmed the extensive conservation between SSC2 and HSA11 and HSA19 and refined the homology with HSA5 (three blocks defined). Furthermore the SSC2q pericentromeric region was shown to be homologous to another human chromosome (HSA1). A complex organization of these syntenies was demonstrated on SSC2q. Our strategy led us to improve also the SSC16 RH map by adding 45 markers. Two-color fluorescence in situ hybridization of markers representative of each synteny confirmed block order. Finally, 29 breakpoints were identified in both species, and porcine BACs containing two breakpoints were isolated.

The genome and the nucleus: a marriage made by evolution. Genome organisation and nuclear architecture

Genomes are housed within cell nuclei as individual chromosome territories. Nuclei contain several architectural structures that interact and influence the genome. In this review, we discuss how the genome may be organised within its nuclear environment with the position of chromosomes inside nuclei being either influenced by gene density or by chromosomes size. We compare interphase genome organisation in diverse species and reveal similarities and differences between evolutionary divergent organisms. Genome organisation is also discussed with relevance to regulation of gene expression, development and differentiation and asks whether large movements of whole chromosomes are really observed during differentiation. Literature and data describing alterations to genome organisation in disease are also discussed. Further, the nuclear structures that are involved in genome function are described, with reference to what happens to the genome when these structures contain protein from mutant genes as in the laminopathies.