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

The Innovative Informatics Approaches of High-Throughput Technologies in Livestock: Spearheading the Sustainability and Resiliency of Agrigenomics Research

For more than a decade, next-generation sequencing (NGS) has been emerging as the mainstay of agrigenomics research. High-throughput technologies have made it feasible to facilitate research at the scale and cost required for using this data in livestock research. Scale frameworks of sequencing for agricultural and livestock improvement, management, and conservation are partly attributable to innovative informatics methodologies and advancements in sequencing practices. Genome-wide sequence-based investigations are often conducted worldwide, and several databases have been created to discover the connections between worldwide scientific accomplishments. Such studies are beginning to provide revolutionary insights into a new era of genomic prediction and selection capabilities of various domesticated livestock species. In this concise review, we provide selected examples of the current state of sequencing methods, many of which are already being used in animal genomic studies, and summarize the state of the positive attributes of genome-based research for cattle (Bos taurus), sheep (Ovis aries), pigs (Sus scrofa domesticus), horses (Equus caballus), chickens (Gallus gallus domesticus), and ducks (Anas platyrhyncos). This review also emphasizes the advantageous features of sequencing technologies in monitoring and detecting infectious zoonotic diseases. In the coming years, the continued advancement of sequencing technologies in livestock agrigenomics will significantly influence the sustained momentum toward regulatory approaches that encourage innovation to ensure continued access to a safe, abundant, and affordable food supplies for future generations.

The Oncopig as an Emerging Model to Investigate Copper Regulation in Cancer

Emerging evidence points to several fundamental contributions that copper (Cu) has to promote the development of human pathologies such as cancer. These recent and increasing identification of the roles of Cu in cancer biology highlights a promising field in the development of novel strategies against cancer. Cu and its network of regulatory proteins are involved in many different contextual aspects of cancer from driving cell signaling, modulating cell cycle progression, establishing the epithelial-mesenchymal transition, and promoting tumor growth and metastasis. Human cancer research in general requires refined models to bridge the gap between basic science research and meaningful clinical trials. Classic studies in cultured cancer cell lines and animal models such as mice and rats often present caveats when extended to humans due to inherent genetic and physiological differences. However, larger animal models such as pigs are emerging as more appropriate tools for translational research as they present more similarities with humans in terms of genetics, anatomical structures, organ sizes, and pathological manifestations of diseases like cancer. These similarities make porcine models well-suited for addressing long standing questions in cancer biology as well as in the arena of novel drug and therapeutic development against human cancers. With the emergent roles of Cu in human health and pathology, the pig presents an emerging and valuable model to further investigate the contributions of this metal to human cancers. The Oncopig Cancer Model is a transgenic swine model that recapitulates human cancer through development of site and cell specific tumors. In this review, we briefly outline the relationship between Cu and cancer, and how the novel Oncopig Cancer Model may be used to provide a better understanding of the mechanisms and causal relationships between Cu and molecular targets involved in cancer.

Non-responsiveness to cardioprotection by ischaemic preconditioning in Ossabaw minipigs with genetic predisposition to, but without the phenotype of the metabolic syndrome

The translation of successful preclinical and clinical proof-of-concept studies on cardioprotection to the benefit of patients with reperfused acute myocardial infarction has been difficult so far. This difficulty has been attributed to confounders which patients with myocardial infarction typically have but experimental animals usually not have. The metabolic syndrome is a typical confounder. We hypothesised that there may also be a genuine non-responsiveness to cardioprotection and used Ossabaw minipigs which have the genetic predisposition to develop a diet-induced metabolic syndrome, but before they had developed the diseased phenotype. Using a prospective study design, a reperfused acute myocardial infarction was induced in 62 lean Ossabaw minipigs by 60 min coronary occlusion and 180 min reperfusion. Ischaemic preconditioning by 3 cycles of 5 min coronary occlusion and 10 min reperfusion was used as cardioprotective intervention. Ossabaw minipigs were stratified for their single nucleotide polymorphism as homozygous for valine (V/V) or isoleucine (I/I)) in the γ-subunit of adenosine monophosphate-activated protein kinase. Endpoints were infarct size and area of no-reflow. Infarct size (V/V: 54 ± 8, I/I: 54 ± 13% of area at risk, respectively) was not reduced by ischaemic preconditioning (V/V: 55 ± 11, I/I: 46 ± 11%) nor was the area of no-reflow (V/V: 57 ± 18, I/I: 49 ± 21 vs. V/V: 57 ± 21, I/I: 47 ± 21% of infarct size). Bioinformatic comparison of the Ossabaw genome to that of Sus scrofa and Göttingen minipigs identified differences in clusters of genes encoding mitochondrial and inflammatory proteins, including the janus kinase (JAK)-signal transducer and activator of transcription (STAT) pathway. The phosphorylation of STAT3 at early reperfusion was not increased by ischaemic preconditioning, different from the established STAT3 activation by cardioprotective interventions in other pig strains. Ossabaw pigs have not only the genetic predisposition to develop a metabolic syndrome but also are not amenable to cardioprotection by ischaemic preconditioning.

Trajectory of livestock genomics in South Asia: A comprehensive review

Livestock plays a central role in sustaining human livelihood in South Asia. There are numerous and distinct livestock species in South Asian countries. Several of them have experienced genetic development in recent years due to the application of genomic technologies and effective breeding programs. This review discusses genomic studies on cattle, buffalo, sheep, goat, pig, horse, camel, yak, mithun, and poultry. The frontiers covered in this review are genetic diversity, admixture studies, selection signature research, QTL discovery, genome-wide association studies (GWAS), and genomic selection. The review concludes with recommendations for South Asian livestock systems to increasingly leverage genomic technologies, based on the lessons learned from the numerous case studies. This paper aims to present a comprehensive analysis of the dichotomy in the South Asian livestock sector and argues that a realistic approach to genomics in livestock can ensure long-term genetic advancements.

Temporal and sex-dependent gene expression patterns in a renal ischemia-reperfusion injury and recovery pig model

Men are more prone to acute kidney injury (AKI) and chronic kidney disease (CKD), progressing to end-stage renal disease (ESRD) than women. Severity and capacity to regenerate after AKI are important determinants of CKD progression, and of patient morbidity and mortality in the hospital setting. To determine sex differences during injury and recovery we have generated a female and male renal ischemia/reperfusion injury (IRI) pig model, which represents a major cause of AKI. Although no differences were found in blood urea nitrogen (BUN) and serum creatinine (SCr) levels between both sexes, females exhibited higher mononuclear infiltrates at basal and recovery, while males showed more tubular damage at injury. Global transcriptomic analyses of kidney biopsies from our IRI pig model revealed a sexual dimorphism in the temporal regulation of genes and pathways relevant for kidney injury and repair, which was also detected in human samples. Enrichment analysis of gene sets revealed five temporal and four sexual patterns governing renal IRI and recovery. Overall, this study constitutes an extensive characterization of the time and sex differences occurring during renal IRI and recovery at gene expression level and offers a template of translational value for further study of sexual dimorphism in kidney diseases.

Articular cartilage and synovium may be important sources of post-surgical synovial fluid inflammatory mediators

The primary source of synovial fluid inflammatory mediators is currently unknown and may include different tissues comprising the joint, including the synovium and articular cartilage. Prior work in a porcine model has demonstrated that anterior cruciate ligament (ACL) surgery leads to significant changes in early gene expression in the synovium and articular cartilage, which are the same whether concomitant ligament restoration is performed or not. In this study, 36 Yucatan minipigs underwent ACL surgery, and a custom multiplex assay was used to measure synovial fluid protein levels of MMP-1, MMP-2, MMP-3, MMP-7, MMP-9, MMP-12, MMP-13, IL-1α, IL-2, IL-4, IL-6, IL-8, IL-10, IL-12, IL-18, GM-CSF, and TNFα in 18 animals at 1 and 4 weeks after surgery. Linear regressions were used to evaluate the relationships between synovial fluid protein levels and the previously reported gene expression levels in the articular cartilage and synovium from the same animal cohort. Synovial fluid levels of MMP-13 and IL-6 were significantly correlated with synovial gene expression (P=.003 and P<.001 respectively), while IL-1α levels were significantly correlated with articular cartilage gene expression (P=.037). The synovium may be an important source of MMP-13 and IL-6, and the articular cartilage may be an important source of IL-1α in post-surgical inflammation. In developing treatments for post-surgical inflammation, the synovium may therefore be a promising target for modulating inflammatory mediators such as MMP-13 and IL-6 in the synovial fluid.

[Opportunities and risks of the use of genetic resistances to infectious diseases in pigs - an overview]

Demands for health, performance and welfare in pigs, as well as the desire for consumer protection and reduced antibiotic use, require optimal measures in advance of disease development. This includes, in principle, the use of genetically more resistant lines and breeding animals, whose existence has been proven for a wide range of pathogen-host interactions. In addition, attempts are being made to identify the gene variants responsible for disease resistance in order to force the selection of suitable populations, also using modern biotechnical technics. The present work is intended to provide an overview of the research status achieved in this context and to highlight opportunities and risks for the future.The evaluation of the international literature shows that genetic disease resistance exist in many areas of swine diseases. However, polygenic inheritance, lack of animal models and the influence of environmental factors during evaluation render their implementation in practical breeding programs demanding. This is where modern molecular genetic methods, such as Gene Editing, come into play. Both approaches possess their pros and cons, which are discussed in this paper. The most important infectious diseases in pigs, including general diseases and epizootics, diseases of the respiratory and digestive tract and diseases of the immune system are taken into account.

Remnant of Unrelated Amniote Sex Chromosomal Linkage Sharing on the Same Chromosome in House Gecko Lizards, Providing a Better Understanding of the Ancestral Super-Sex Chromosome

Comparative chromosome maps investigating sex chromosomal linkage groups in amniotes and microsatellite repeat motifs of a male house gecko lizard (Hemidactylus frenatus, HFR) and a flat-tailed house gecko lizard (H. platyurus, HPL) of unknown sex were examined using 75 bacterial artificial chromosomes (BACs) from chicken and zebra finch genomes. No massive accumulations of microsatellite repeat motifs were found in either of the gecko lizards, but 10 out of 13 BACs mapped on HPL chromosomes were associated with other amniote sex chromosomes. Hybridization of the same BACs onto multiple different chromosome pairs suggested transitions to sex chromosomes across amniotes. No BAC hybridization signals were found on HFR chromosomes. However, HFR diverged from HPL about 30 million years ago, possibly due to intrachromosomal rearrangements occurring in the HFR lineage. By contrast, heterochromatin likely reshuffled patterns between HPL and HFR, as observed from C-positive heterochromatin distribution. Six out of ten BACs showed partial homology with squamate reptile chromosome 2 (SR2) and snake Z and/or W sex chromosomes. The gecko lizard showed shared unrelated sex chromosomal linkages-the remnants of a super-sex chromosome. A large ancestral super-sex chromosome showed a correlation between SR2 and snake W sex chromosomes.

Investigating the impact of reference assembly choice on genomic analyses in a cattle breed

BACKGROUND

Reference-guided read alignment and variant genotyping are prone to reference allele bias, particularly for samples that are greatly divergent from the reference genome. A Hereford-based assembly is the widely accepted bovine reference genome. Haplotype-resolved genomes that exceed the current bovine reference genome in quality and continuity have been assembled for different breeds of cattle. Using whole genome sequencing data of 161 Brown Swiss cattle, we compared the accuracy of read mapping and sequence variant genotyping as well as downstream genomic analyses between the bovine reference genome (ARS-UCD1.2) and a highly continuous Angus-based assembly (UOA_Angus_1).

RESULTS

Read mapping accuracy did not differ notably between the ARS-UCD1.2 and UOA_Angus_1 assemblies. We discovered 22,744,517 and 22,559,675 high-quality variants from ARS-UCD1.2 and UOA_Angus_1, respectively. The concordance between sequence- and array-called genotypes was high and the number of variants deviating from Hardy-Weinberg proportions was low at segregating sites for both assemblies. More artefactual INDELs were genotyped from UOA_Angus_1 than ARS-UCD1.2 alignments. Using the composite likelihood ratio test, we detected 40 and 33 signatures of selection from ARS-UCD1.2 and UOA_Angus_1, respectively, but the overlap between both assemblies was low. Using the 161 sequenced Brown Swiss cattle as a reference panel, we imputed sequence variant genotypes into a mapping cohort of 30,499 cattle that had microarray-derived genotypes using a two-step imputation approach. The accuracy of imputation (Beagle R2) was very high (0.87) for both assemblies. Genome-wide association studies between imputed sequence variant genotypes and six dairy traits as well as stature produced almost identical results from both assemblies.

CONCLUSIONS

The ARS-UCD1.2 and UOA_Angus_1 assemblies are suitable for reference-guided genome analyses in Brown Swiss cattle. Although differences in read mapping and genotyping accuracy between both assemblies are negligible, the choice of the reference genome has a large impact on detecting signatures of selection that already reached fixation using the composite likelihood ratio test. We developed a workflow that can be adapted and reused to compare the impact of reference genomes on genome analyses in various breeds, populations and species.

LDscaff: LD-based scaffolding of de novo genome assemblies

Background

Genome assembly is fundamental for de novo genome analysis. Hybrid assembly, utilizing various sequencing technologies increases both contiguity and accuracy. While such approaches require extra costly sequencing efforts, the information provided millions of existed whole-genome sequencing data have not been fully utilized to resolve the task of scaffolding. Genetic recombination patterns in population data indicate non-random association among alleles at different loci, can provide physical distance signals to guide scaffolding.

Results

In this paper, we propose LDscaff for draft genome assembly incorporating linkage disequilibrium information in population data. We evaluated the performance of our method with both simulated data and real data. We simulated scaffolds by splitting the pig reference genome and reassembled them. Gaps between scaffolds were introduced ranging from 0 to 100 KB. The genome misassembly rate is 2.43% when there is no gap. Then we implemented our method to refine the Giant Panda genome and the donkey genome, which are purely assembled by NGS data. After LDscaff treatment, the resulting Panda assembly has scaffold N50 of 3.6 MB, 2.5 times larger than the original N50 (1.3 MB). The re-assembled donkey assembly has an improved N50 length of 32.1 MB from 23.8 MB.

Conclusions

Our method effectively improves the assemblies with existed re-sequencing data, and is an potential alternative to the existing assemblers required for the collection of new data.

Optimization of sperm RNA processing for developmental research

Recent studies have demonstrated the significance of sperm RNA function as a transporter of important information directing the course of life. To determine the message contained in sperm RNA, it is necessary to optimize transcriptomic research tools. The current study was performed to optimize the processing of sperm RNA from sample storage to quantitative real-time PCR and assess the corresponding method with to evaluate male fertility and its representative markers, equatorin (EQTN) and peroxiredoxin (PRDX). Following successive steps of the Minimum Information for Publication of Quantitative Real-Time PCR Experiments guidelines, several options were compared using boar spermatozoa. To evaluate the optimized procedures, the relationship between mRNA expression of EQTN and PRDX in spermatozoa and the fertility (litter size) of 20 individual boars was assessed. Unexpectedly, DNase treatment during RNA isolation had the deleterious effect by decreasing the RNA concentration by 56% and eliminating the correlation between EQTN and PRDX4 mRNA expression and male fertility. Moreover, when sperm RNA was processed using the corresponding method, the results showed the highest exon sequence expression, male fertility prediction power, and consistency. This optimized protocol for predicting male fertility can be used to study the transport of messages directing the life course from spermatozoon to offspring.

An improved pig reference genome sequence to enable pig genetics and genomics research

Background

The domestic pig (Sus scrofa) is important both as a food source and as a biomedical model given its similarity in size, anatomy, physiology, metabolism, pathology, and pharmacology to humans. The draft reference genome (Sscrofa10.2) of a purebred Duroc female pig established using older clone-based sequencing methods was incomplete, and unresolved redundancies, short-range order and orientation errors, and associated misassembled genes limited its utility.

Results

We present 2 annotated highly contiguous chromosome-level genome assemblies created with more recent long-read technologies and a whole-genome shotgun strategy, 1 for the same Duroc female (Sscrofa11.1) and 1 for an outbred, composite-breed male (USMARCv1.0). Both assemblies are of substantially higher (>90-fold) continuity and accuracy than Sscrofa10.2.

Conclusions

These highly contiguous assemblies plus annotation of a further 11 short-read assemblies provide an unprecedented view of the genetic make-up of this important agricultural and biomedical model species. We propose that the improved Duroc assembly (Sscrofa11.1) become the reference genome for genomic research in pigs.

Of Mice, Dogs, Pigs, and Men: Choosing the Appropriate Model for Immuno-Oncology Research

The immune system plays dual roles in response to cancer. The host immune system protects against tumor formation via immunosurveillance; however, recognition of the tumor by immune cells also induces sculpting mechanisms leading to a Darwinian selection of tumor cell variants with reduced immunogenicity. Cancer immunoediting is the concept used to describe the complex interplay between tumor cells and the immune system. This concept, commonly referred to as the three E's, is encompassed by 3 distinct phases of elimination, equilibrium, and escape. Despite impressive results in the clinic, cancer immunotherapy still has room for improvement as many patients remain unresponsive to therapy. Moreover, many of the preclinical results obtained in the widely used mouse models of cancer are lost in translation to human patients. To improve the success rate of immuno-oncology research and preclinical testing of immune-based anticancer therapies, using alternative animal models more closely related to humans is a promising approach. Here, we describe 2 of the major alternative model systems: canine (spontaneous) and porcine (experimental) cancer models. Although dogs display a high rate of spontaneous tumor formation, an increased number of genetically modified porcine models exist. We suggest that the optimal immuno-oncology model may depend on the stage of cancer immunoediting in question. In particular, the spontaneous canine tumor models provide a unique platform for evaluating therapies aimed at the escape phase of cancer, while genetically engineered swine allow for elucidation of tumor-immune cell interactions especially during the phases of elimination and equilibrium.

A Nutrigenomics Approach Using RNA Sequencing Technology to Study Nutrient-Gene Interactions in Agricultural Animals

Thorough understanding of animal gene expression driven by dietary nutrients can be regarded as a bottom line of advanced animal nutrition research. Nutrigenomics (including transcriptomics) studies the effects of dietary nutrients on cellular gene expression and, ultimately, phenotypic changes in living organisms. Transcriptomics can be applied to investigate animal tissue transcriptomes at a defined nutritional state, which can provide a holistic view of intracellular RNA expression. As a novel transcriptomics approach, RNA sequencing (RNA-Seq) technology can monitor all gene expressions simultaneously in response to dietary intervention. The principle and history of RNA-Seq are briefly reviewed, and its 3 principal steps are described in this article. Application of RNA-Seq in different areas of animal nutrition research is summarized. Lastly, the application of RNA-Seq in swine science and nutrition is also reviewed. In short, RNA-Seq holds significant potential to be employed for better understanding the nutrient-gene interactions in agricultural animals.

Translating Human Cancer Sequences Into Personalized Porcine Cancer Models

The global incidence of cancer is rapidly rising, and despite an improved understanding of cancer molecular biology, immune landscapes, and advancements in cytotoxic, biologic, and immunologic anti-cancer therapeutics, cancer remains a leading cause of death worldwide. Cancer is caused by the accumulation of a series of gene mutations called driver mutations that confer selective growth advantages to tumor cells. As cancer therapies move toward personalized medicine, predictive modeling of the role driver mutations play in tumorigenesis and therapeutic susceptibility will become essential. The development of next-generation sequencing technology has made the evaluation of mutated genes possible in clinical practice, allowing for identification of driver mutations underlying cancer development in individual patients. This, combined with recent advances in gene editing technologies such as CRISPR-Cas9 enables development of personalized tumor models for prediction of treatment responses for mutational profiles observed clinically. Pigs represent an ideal animal model for development of personalized tumor models due to their similar size, anatomy, physiology, metabolism, immunity, and genetics compared to humans. Such models would support new initiatives in precision medicine, provide approaches to create disease site tumor models with designated spatial and temporal clinical outcomes, and create standardized tumor models analogous to human tumors to enable therapeutic studies. In this review, we discuss the process of utilizing genomic sequencing approaches, gene editing technologies, and transgenic porcine cancer models to develop clinically relevant, personalized large animal cancer models for use in co-clinical trials, ultimately improving treatment stratification and translation of novel therapeutic approaches to clinical practice.

Sequencing through thick and thin: Historiographical and philosophical implications

DNA sequencing has been characterised by scholars and life scientists as an example of 'big', 'fast' and 'automated' science in biology. This paper argues, however, that these characterisations are a product of a particular interpretation of what sequencing is, what I call 'thin sequencing'. The 'thin sequencing' perspective focuses on the determination of the order of bases in a particular stretch of DNA. Based upon my research on the pig genome mapping and sequencing projects, I provide an alternative 'thick sequencing' perspective, which also includes a number of practices that enable the sequence to travel across and be used in wider communities. If we take sequencing in the thin manner to be an event demarcated by the determination of sequences in automated sequencing machines and computers, this has consequences for the historical analysis of sequencing projects, as it focuses attention on those parts of the work of sequencing that are more centralised, fast (and accelerating) and automated. I argue instead that sequencing can be interpreted as a more open-ended process including activities such as the generation of a minimum tile path or annotation, and detail the historiographical and philosophical consequences of this move.

Contribution of the swine model in the study of human sexually transmitted infections

The pig has garnered more and more interest as a model animal to study various conditions in humans. The growing success of the pig as an experimental animal model is explained by its similarities with humans in terms of anatomy, genetics, immunology, and physiology, by their manageable behavior and size, and by the general public acceptance of using pigs for experimental purposes. In addition, the immunological toolbox of pigs has grown substantially in the last decade. This development led to a boost in the use of pigs as a preclinical model for various human infections including sexually transmitted diseases (STIs) like Chlamydia trachomatis. In the current review, we discuss the use of animal models for biomedical research on the major human STIs. We summarize results obtained in the most common animal models and focus on the contributions of the pig model towards the understanding of pathogenesis and the host immune response. In addition, we present the main features of the porcine model that are particularly relevant for the study of pathogens affecting human female and male genital tracts. We also inform on the technological advancements in the porcine toolbox to facilitate new discoveries in this biologically important animal model. There is a continued need for improvements in animal modeling for biomedical research inclusive STI research. With all its advantages and the highly improved toolbox, the porcine model can play a crucial role in STI research and open the door to new exciting discoveries.

The Oncopig Cancer Model: An Innovative Large Animal Translational Oncology Platform

Despite an improved understanding of cancer molecular biology, immune landscapes, and advancements in cytotoxic, biologic, and immunologic anti-cancer therapeutics, cancer remains a leading cause of death worldwide. More than 8.2 million deaths were attributed to cancer in 2012, and it is anticipated that cancer incidence will continue to rise, with 19.3 million cases expected by 2025. The development and investigation of new diagnostic modalities and innovative therapeutic tools is critical for reducing the global cancer burden. Toward this end, transitional animal models serve a crucial role in bridging the gap between fundamental diagnostic and therapeutic discoveries and human clinical trials. Such animal models offer insights into all aspects of the basic science-clinical translational cancer research continuum (screening, detection, oncogenesis, tumor biology, immunogenicity, therapeutics, and outcomes). To date, however, cancer research progress has been markedly hampered by lack of a genotypically, anatomically, and physiologically relevant large animal model. Without progressive cancer models, discoveries are hindered and cures are improbable. Herein, we describe a transgenic porcine model—the Oncopig Cancer Model (OCM)—as a next-generation large animal platform for the study of hematologic and solid tumor oncology. With mutations in key tumor suppressor and oncogenes, TP53^R167H and KRAS^G12D, the OCM recapitulates transcriptional hallmarks of human disease while also exhibiting clinically relevant histologic and genotypic tumor phenotypes. Moreover, as obesity rates increase across the global population, cancer patients commonly present clinically with multiple comorbid conditions. Due to the effects of these comorbidities on patient management, therapeutic strategies, and clinical outcomes, an ideal animal model should develop cancer on the background of representative comorbid conditions (tumor macro- and microenvironments). As observed in clinical practice, liver cirrhosis frequently precedes development of primary liver cancer or hepatocellular carcinoma. The OCM has the capacity to develop tumors in combination with such relevant comorbidities. Furthermore, studies on the tumor microenvironment demonstrate similarities between OCM and human cancer genomic landscapes. This review highlights the potential of this and other large animal platforms as transitional models to bridge the gap between basic research and clinical practice.

Engineered Swine Models of Cancer

Over the past decade, the technology to engineer genetically modified swine has seen many advancements, and because their physiology is remarkably similar to that of humans, swine models of cancer may be extremely valuable for preclinical safety studies as well as toxicity testing of pharmaceuticals prior to the start of human clinical trials. Hence, the benefits of using swine as a large animal model in cancer research and the potential applications and future opportunities of utilizing pigs in cancer modeling are immense. In this review, we discuss how pigs have been and can be used as a biomedical models for cancer research, with an emphasis on current technologies. We have focused on applications of precision genetics that can provide models that mimic human cancer predisposition syndromes. In particular, we describe the advantages of targeted gene-editing using custom endonucleases, specifically TALENs and CRISPRs, and transposon systems, to make novel pig models of cancer with broad preclinical applications.

A decade of pig genome sequencing: a window on pig domestication and evolution

Insight into how genomes change and adapt due to selection addresses key questions in evolutionary biology and in domestication of animals and plants by humans. In that regard, the pig and its close relatives found in Africa and Eurasia represent an excellent group of species that enables studies of the effect of both natural and human-mediated selection on the genome. The recent completion of the draft genome sequence of a domestic pig and the development of next-generation sequencing technology during the past decade have created unprecedented possibilities to address these questions in great detail. In this paper, I review recent whole-genome sequencing studies in the pig and closely-related species that provide insight into the demography, admixture and selection of these species and, in particular, how domestication and subsequent selection of Sus scrofa have shaped the genomes of these animals.

Stories and Challenges of Genome Wide Association Studies in Livestock - A Review

Undoubtedly livestock is one of the major contributors to the economy of any country. The economic value of livestock includes meat, dairy products, fiber, fertilizer etc. Understanding and identifying the associations of quantitative trait loci (QTL) with the economically important traits is believed to substantially benefit the livestock industry. The past two decades have seen a flurry of interest in mapping the QTL associated with traits of economic importance on the genome. With the availability of single nucleotide polymorphism chip of various densities it is possible to identify regions, QTL and genes on the genome that explain the association and its effect on the phenotype under consideration. Remarkable advancement has been seen in genome wide association studies (GWAS) since its inception till the present day. In this review we describe the progress and challenges of GWAS in various livestock species.

Bioinformatics Annotation of Human Y Chromosome-Encoded Protein Pathways and Interactions

We performed a comprehensive analysis of human Y chromosome-encoded proteins, their pathways, and their interactions using bioinformatics tools. From the NCBI annotation release 107 of human genome, we retrieved a total of 66 proteins encoded on Y chromosome. Most of the retrieved proteins were also matched with the proteins listed in the core databases of the Human Proteome Project including neXtProt, PeptideAtlas, and the Human Protein Atlas. When we examined the pathways of human Y-encoded proteins through KEGG database and Pathway Studio software, many of proteins fall into the categories related to cell signaling pathways. Using the STRING program, we found a total of 49 human Y-encoded proteins showing strong/medium interaction with each other. While using the Pathway studio software, we found that a total of 16 proteins interact with other chromosome-encoded proteins. In particular, the SRY protein interacted with 17 proteins encoded on other chromosomes. Additionally, we aligned the sequences of human Y-encoded proteins with the sequences of chimpanzee and mouse Y-encoded proteins using the NCBI BLAST program. This analysis resulted in a significant number of orthologous proteins between human, chimpanzee, and mouse. Collectively, our findings provide the scientific community with additional information on the human Y chromosome-encoded proteins.

Mining the pig genome to investigate the domestication process

Pig domestication began around 9000 YBP in the Fertile Crescent and Far East, involving marked morphological and genetic changes that occurred in a relatively short window of time. Identifying the alleles that drove the behavioural and physiological transformation of wild boars into pigs through artificial selection constitutes a formidable challenge that can only be faced from an interdisciplinary perspective. Indeed, although basic facts regarding the demography of pig domestication and dispersal have been uncovered, the biological substrate of these processes remains enigmatic. Considerable hope has been placed on new approaches, based on next-generation sequencing, which allow whole-genome variation to be analyzed at the population level. In this review, we provide an outline of the current knowledge on pig domestication by considering both archaeological and genetic data. Moreover, we discuss several potential scenarios of genome evolution under the complex mixture of demography and selection forces at play during domestication. Finally, we highlight several technical and methodological approaches that may represent significant advances in resolving the conundrum of livestock domestication.

Genome sequencing and analysis of Mangalica, a fatty local pig of Hungary

BACKGROUND

Mangalicas are fatty type local/rare pig breeds with an increasing presence in the niche pork market in Hungary and in other countries. To explore their genetic resources, we have analysed data from next-generation sequencing of an individual male from each of three Mangalica breeds along with a local male Duroc pig. Structural variations, such as SNPs, INDELs and CNVs, were identified and particular genes with SNP variations were analysed with special emphasis on functions related to fat metabolism in pigs.

RESULTS

More than 60 Gb of sequence data were generated for each of the sequenced individuals, resulting in 11× to 19× autosomal median coverage. After stringent filtering, around six million SNPs, of which approximately 10% are novel compared to the dbSNP138 database, were identified in each animal. Several hundred thousands of INDELs and about 1,000 CNV gains were also identified. The functional annotation of genes with exonic, non-synonymous SNPs, which are common in all three Mangalicas but are absent in either the reference genome or the sequenced Duroc of this study, highlighted 52 genes in lipid metabolism processes. Further analysis revealed that 41 of these genes are associated with lipid metabolic or regulatory pathways, 49 are in fat-metabolism and fatness-phenotype QTLs and, with the exception of ACACA, ANKRD23, GM2A, KIT, MOGAT2, MTTP, FASN, SGMS1, SLC27A6 and RETSAT, have not previously been associated with fat-related phenotypes.

CONCLUSIONS

Genome analysis of Mangalica breeds revealed that local/rare breeds could be a rich source of sequence variations not present in cosmopolitan/industrial breeds. The identified Mangalica variations may, therefore, be a very useful resource for future studies of agronomically important traits in pigs.

Generation and characterization of a transgenic pig carrying a DsRed-monomer reporter gene

BACKGROUND

Pigs are an optimal animal for conducting biomedical research because of their anatomical and physiological resemblance to humans. In contrast to the abundant resources available in the study of mice, few fluorescent protein-harboring porcine models are available for preclinical studies. In this paper, we report the successful generation and characterization of a transgenic DsRed-Monomer porcine model.

METHODS

The transgene comprised a CMV enhancer/chicken-beta actin promoter and DsRed monomeric cDNA. Transgenic pigs were produced by using pronuclear microinjection. PCR and Southern blot analyses were applied for identification of the transgene. Histology, blood examinations and computed tomography were performed to study the health conditions. The pig amniotic fluid progenitor/stem cells were also isolated to examine the existence of red fluorescence and differentiation ability.

RESULTS

Transgenic pigs were successfully generated and transmitted to offspring at a germ-line transmission rate of 43.59% (17/39). Ubiquitous expression of red fluorescence was detected in the brain, eye, tongue, heart, lung, liver, pancreas, spleen, stomach, small intestine, large intestine, kidney, testis, and muscle; this was confirmed by histology and western blot analyses. In addition, we confirmed the differentiation potential of amniotic fluid progenitor stem cells isolated from the transgenic pig.

CONCLUSIONS

This red fluorescent pig can serve as a host for other fluorescent-labeled cells in order to study cell-microenvironment interactions, and can provide optimal red-fluorescent-labeled cells and tissues for research in developmental biology, regenerative medicine, and xenotransplantation.

Evaluation of the gamma-H2AX assay for radiation biodosimetry in a swine model

There is a paucity of large animal models to study both the extent and the health risk of ionizing radiation exposure in humans. One promising candidate for such a model is the minipig. Here, we evaluate the minipig for its potential in γ-H2AX-based biodosimetry after exposure to ionizing radiation using both Cs137 and Co60 sources. γ-H2AX foci were enumerated in blood lymphocytes and normal fibroblasts of human and porcine origin after ex vivo γ-ray irradiation. DNA double-strand break repair kinetics in minipig blood lymphocytes and fibroblasts, based on the γ-H2AX assay, were similar to those observed in their human counterparts. To substantiate the similarity observed between the human and minipig we show that minipig fibroblast radiosensitivity was similar to that observed with human fibroblasts. Finally, a strong γ-H2AX induction was observed in blood lymphocytes following minipig total body irradiation. Significant responses were detected 3 days after 1.8 Gy and 1 week after 3.8 and 5 Gy with residual γ-H2AX foci proportional to the initial radiation doses. These findings show that the Gottingen minipig provides a useful in vivo model for validation of γ-H2AX biodosimetry for dose assessment in humans.

Large-scale sequencing based on full-length-enriched cDNA libraries in pigs: contribution to annotation of the pig genome draft sequence

BACKGROUND

Along with the draft sequencing of the pig genome, which has been completed by an international consortium, collection of the nucleotide sequences of genes expressed in various tissues and determination of entire cDNA sequences are necessary for investigations of gene function. The sequences of expressed genes are also useful for genome annotation, which is important for isolating the genes responsible for particular traits.

RESULTS

We performed a large-scale expressed sequence tag (EST) analysis in pigs by using 32 full-length-enriched cDNA libraries derived from 28 kinds of tissues and cells, including seven tissues (brain, cerebellum, colon, hypothalamus, inguinal lymph node, ovary, and spleen) derived from pigs that were cloned from a sow subjected to genome sequencing. We obtained more than 330,000 EST reads from the 5'-ends of the cDNA clones. Comparison with human and bovine gene catalogs revealed that the ESTs corresponded to at least 15,000 genes. cDNA clones representing contigs and singlets generated by assembly of the EST reads were subjected to full-length determination of inserts. We have finished sequencing 31,079 cDNA clones corresponding to more than 12,000 genes. Mapping of the sequences of these cDNA clones on the draft sequence of the pig genome has indicated that the clones are derived from about 15,000 independent loci on the pig genome.

CONCLUSIONS

ESTs and cDNA sequences derived from full-length-enriched libraries are valuable for annotation of the draft sequence of the pig genome. This information will also contribute to the exploration of promoter sequences on the genome and to molecular biology-based analyses in pigs.

Analyses of pig genomes provide insight into porcine demography and evolution

For 10,000 years pigs and humans have shared a close and complex relationship. From domestication to modern breeding practices, humans have shaped the genomes of domestic pigs. Here we present the assembly and analysis of the genome sequence of a female domestic Duroc pig (Sus scrofa) and a comparison with the genomes of wild and domestic pigs from Europe and Asia. Wild pigs emerged in South East Asia and subsequently spread across Eurasia. Our results reveal a deep phylogenetic split between European and Asian wild boars ∼1 million years ago, and a selective sweep analysis indicates selection on genes involved in RNA processing and regulation. Genes associated with immune response and olfaction exhibit fast evolution. Pigs have the largest repertoire of functional olfactory receptor genes, reflecting the importance of smell in this scavenging animal. The pig genome sequence provides an important resource for further improvements of this important livestock species, and our identification of many putative disease-causing variants extends the potential of the pig as a biomedical model.

[Manual annotation of the pig whole genomic sequence using Otterlace software]

In November 2009, scientists from the US, UK, and other countries announced the complete genome sequence draft of the domestic pig. With the release of improved versions of the pig genome assembly and the increase of correctly assembled sequenced fragments over the past two years, it is particularly urgent to have the pig genes annotated at whole-genome level. This article is aimed at introducing an excellent manual annotation tool, Otterlace software, developed by Sanger institute. We used CFL1 (Cofilin 1) gene as an example to expound the usage of the three main components of Otterlace, Zmap, Blixem, and Dotter tools, and developed a practical procedure for manual annotations. We have analyzed 243 immune-related genes, among which 180 genes have been completely or partially annotated, offering novel information to the porcine functional genomics.

Oligonucleotide microarray analysis of dietary-induced hyperlipidemia gene expression profiles in miniature pigs

Background

Hyperlipidemia animal models have been established, but complete gene expression profiles of the transition from normal lipid levels have not been obtained. Miniature pigs are useful model animals for gene expression studies on dietary-induced hyperlipidemia because they have a similar anatomy and digestive physiology to humans, and blood samples can be obtained from them repeatedly.

Methodology

Two typical dietary treatments were used for dietary-induced hyperlipidemia models, by using specific pathogen-free (SPF) Clawn miniature pigs. One was a high-fat and high-cholesterol diet (HFCD) and the other was a high-fat, high-cholesterol, and high-sucrose diet (HFCSD). Microarray analyses were conducted from whole blood samples during the dietary period and from white blood cells at the end of the dietary period to evaluate the transition of expression profiles of the two dietary models.

Principal Findings

Variations in whole blood gene expression intensity within the HFCD or the HFCSD group were in the same range as the controls provide with normal diet at all periods. This indicates uniformity of dietary-induced hyperlipidemia for our dietary protocols. Gene ontology- (GO) based functional analyses revealed that characteristics of the common changes between HFCD and HFCSD were involved in inflammatory responses and reproduction. The correlation coefficient between whole blood and white blood cell expression profiles at 27 weeks with the HFCSD diet was significantly lower than that of the control and HFCD diet groups. This may be due to the effects of RNA originating from the tissues and/or organs.

Conclusions

No statistically significant differences in fasting plasma lipids and glucose levels between the HFCD and HFCSD groups were observed. However, blood RNA analyses revealed different characteristics corresponding to the dietary protocols. In this study, whole blood RNA analyses proved to be a useful tool to evaluate transitions in dietary-induced hyperlipidemia gene expression profiles in miniature pigs.

Organization, complexity and allelic diversity of the porcine (Sus scrofa domestica) immunoglobulin lambda locus

We have characterized the organization, complexity, and expression of the porcine (Sus scrofa domestica) immunoglobulin lambda (IGL) light chain locus, which accounts for about half of antibody light chain usage in swine, yet is nearly totally unknown. Twenty-two IGL variable (IGLV) genes were identified that belong to seven subgroups. Nine genes appear to be functional. Eight possess stop codons, frameshifts, or both, and one is missing the V-EXON. Two additional genes are missing an essential cysteine residue and are classified as ORF (open reading frame). The IGLV genes are organized in two distinct clusters, a constant (C)-proximal cluster dominated by genes similar to the human IGLV3 subgroup, and a C-distal cluster dominated by genes most similar to the human IGLV8 and IGLV5 subgroups. Phylogenetic analysis reveals that the porcine IGLV8 subgroup genes have recently expanded, suggesting a particularly effective role in immunity to porcine-specific pathogens. Moreover, expression of IGLV genes is nearly exclusively restricted to the IGLV3 and IGLV8 genes. The constant locus comprises three tandem cassettes comprised of a joining (IGLJ) gene and a constant (IGLC) gene, whereas a fourth downstream IGLJ gene has no corresponding associated IGLC gene. Comparison of individual BACs generated from the same individual revealed polymorphisms in IGLC2 and several IGLV genes, indicating that allelic variation in IGLV further expands the porcine antibody light chain repertoire.

Evolution of the porcine (Sus scrofa domestica) immunoglobulin kappa locus through germline gene conversion

Immunoglobulin (IG) gene rearrangement and expression are central to disease resistance and health maintenance in animals. The IG kappa (IGK) locus in swine (Sus scrofa domestica) contributes to approximately half of all antibody molecules, in contrast to many other Cetartiodactyla, whose members provide the majority of human dietary protein and in which kappa locus utilization is limited. The porcine IGK variable locus is 27.9 kb upstream of five IG kappa J genes (IGKJ) which are separated from a single constant gene (IGKC) by 2.8 kb. Fourteen variable genes (IGKV) were identified, of which nine are functional and two are open reading frame (ORF). Of the three pseudogenes, IGKV3-1 contains a frameshift and multiple stop codons, IGKV7-2 contains multiple stop codons, and IGKV2-5 is missing exon 2. The nine functional IGKV genes are phylogenetically related to either the human IGKV1 or IGKV2 subgroups. IGKV2 subgroup genes were found to be dominantly expressed. Polymorphisms were identified on overlapping BACs derived from the same individual such that 11 genes contain amino acid differences. The most striking allelic differences are present in IGKV2 genes, which contain as many as 16 amino acid changes between alleles, the majority of which are in complementarity determining region (CDR) 1. In addition, many IGKV2 CDR1 are shared between genes but not between alleles, suggesting extensive diversification of this locus through gene conversion.

Characterization of porcine autism susceptibility candidate 2 as a candidate gene for the number of corpora lutea in pigs

In a previous study, we mapped two quantitative trait loci (QTL) approximately 50cM apart, both influencing the number of corpora lutea in pigs on chromosome 3. One locus included functional candidate genes for proteins related to specific aspects of fertility, such as the follicle-stimulating hormone receptor and the luteinizing hormone/choriogonadotropin receptor. However, specific genes related to the second locus have not yet been identified. This study aims to identify another candidate gene influencing the number of corpora lutea in pigs. Using 12 polymorphic markers, we fine-mapped a narrow region of pig chromosome 3 that had been shown to contain a QTL for corpora lutea. In the critical region, only 1 gene, autism susceptibility candidate 2 (AUTS2), was identified as a positional candidate. Our results demonstrate that the porcine AUTS2 gene consists of 19 exons with a complete open reading frame of 3768bp encoding an AUTS2 protein of 1256 amino acids. We screened the whole coding sequence and parts of the untranslated region for polymorphisms in an F(2) population of Duroc×Meishan crosses. We found 1 ins/del and 7 single nucleotide polymorphisms (SNP), including 2 nonsynonymous variants, c.943C>T in exon 7 and c.2828C>T in exon 19, resulting in P315S and A943V, respectively. The SNP c.943C>T within a proline-rich domain was genotyped in several breeds; the C allele occurred in all breeds, whereas the T allele occurred only in Meishan pigs. Using in situ hybridization, the mRNA expression of the AUTS2 gene was observed on granulosa cells in the porcine ovary and thus may be associated with hormone sensitivity.

Oligonucleotide microarray analysis of age-related gene expression profiles in miniature pigs

Miniature pigs are useful model animals for humans because they have similar anatomy and digestive physiology to humans and are easy to breed and handle. In this study, whole blood microarray analyses were conducted to evaluate variations of correlation among individuals and ages using specific pathogen-free (SPF) Clawn miniature pigs. Whole blood RNA is easy to handle compared to isolated white blood cell RNA and can be used for health and disease monitoring and animal control. In addition, whole blood is a heterogeneous mixture of subpopulation cells. Once a great change occurs in composition and expressing condition of subpopulations, their associated change will be reflected on whole blood RNA. From 12 to 30 weeks of age, fractions of lymphocytes, monocytes, neutrophils, eosinophils, and basophils in white blood cells showed insignificant differences with age as a result of ANOVA analysis. This study attempted to identify characteristics of age-related gene expression by taking into account the change in the number of expressed genes by age and similarities of gene expression intensity between individuals. As a result, the number of expressed genes was less in fetal stage and infancy period but increased with age, reaching a steady state of gene expression after 20 weeks of age. Variation in gene expression intensity within the same age was great in fetal stage and infancy period, but converged with age. The variation between 20 and 30 weeks of age was comparable to that among 30 weeks individuals. These results indicate that uniformity of laboratory animals is expected for miniature pigs after 20 weeks of age. Furthermore, a possibility was shown that whole blood RNA analysis is applicable to evaluation of physiological state.

Analysis of porcine MHC using microarrays

The major histocompatibility complex (MHC) in Mammals is one of the most gene dense regions of the genome and contains the polymorphic histocompatibility gene families known to be involved in pathogen response and control of auto-immunity. The MHC is a complex genetic system that provides an interesting model system to study genome expression regulation and genetic diversity at the megabase scale. The pig MHC or SLA (Swine Leucocyte Antigen) complex spans 2.4 megabases and 151 loci have been annotated. We will review key results from previous RNA expression studies using microarrays containing probes specific to annotated loci within SLA and in addition present novel data obtained using high-density tiling arrays encompassing the whole SLA complex. We have focused on transcriptome modifications of porcine peripheral blood mononuclear cells stimulated with a mixture of phorbol myristate acetate and ionomycin known to activate B and T cell proliferation. Our results show that numerous loci mapping to the SLA complex are affected by the treatment. A general decreased level of expression for class I and II genes and an up-regulation of genes involved in peptide processing and transport were observed. Tiling array-based experiments contributed to refined gene annotations as presented for one SLA class I gene referred to as SLA-11. In conclusion, high-density tiling arrays can serve as an excellent tool to draw comprehensive transcription maps, and improve genome annotations for the SLA complex. We are currently studying their relevance to characterize SLA genetic diversity in combination with high throughput next generation sequencing.

Contribution of large pig for renal ischemia-reperfusion and transplantation studies: the preclinical model

Animal experimentation is necessary to characterize human diseases and design adequate therapeutic interventions. In renal transplantation research, the limited number of in vitro models involves a crucial role for in vivo models and particularly for the porcine model. Pig and human kidneys are anatomically similar (characterized by multilobular structure in contrast to rodent and dog kidneys unilobular). The human proximity of porcine physiology and immune systems provides a basic knowledge of graft recovery and inflammatory physiopathology through in vivo studies. In addition, pig large body size allows surgical procedures similar to humans, repeated collections of peripheral blood or renal biopsies making pigs ideal for medical training and for the assessment of preclinical technologies. However, its size is also its main drawback implying expensive housing. Nevertheless, pig models are relevant alternatives to primate models, offering promising perspectives with developments of transgenic modulation and marginal donor models facilitating data extrapolation to human conditions.

An intronic polymorphism in the porcine IRF7 gene is associated with better health and immunity of the host during Sarcocystis infection, and affects interferon signalling

Interferon regulatory factor 7 (IRF7), as a key regulator of type I interferon response, plays an important role during innate response against viral infection. Although well conserved across species, the structure of IRF7 and its function during parasite infection are not well documented in farm animals, such as the pig. To bridge this gap, we have determined the porcine IRF7 gene structure and identified two intronic single nucleotide polymorphisms (SNPs), SNP g.748G>C and SNP g.761A>G, in commercial pig breeds. The distribution of SNP g.761A>G in multiple breeds suggested that it was in Hardy-Weinberg equilibrium and allowed us to map it at the top of SSC2. We found that during Sarcocystis miescheriana infection, the G allele was associated with high lymphocyte levels (P < 0.02), reduced drop in platelet levels (P < 0.002) and IgG1-Th2-dominated response (P < 0.05). This suggests that the G allele was associated with better health and immunity of the host during Sarcocystis infection. Furthermore, we have also provided suggestive evidence that the G allele of SNPc.761A>G enhances the transactivation activity of IRF7, possibly by improving IRF7 transcript splicing of intron-3. These findings would suggest that IRF7, as a transcriptional regulator, is involved in the defence mechanism against a larger spectrum of pathogens, and in more host species, than initially anticipated.

The pig genome project has plenty to squeal about

Significant progress on pig genetics and genomics research has been witnessed in recent years due to the integration of advanced molecular biology techniques, bioinformatics and computational biology, and the collaborative efforts of researchers in the swine genomics community. Progress on expanding the linkage map has slowed down, but the efforts have created a higher-resolution physical map integrating the clone map and BAC end sequence. The number of QTL mapped is still growing and most of the updated QTL mapping results are available through PigQTLdb. Additionally, expression studies using high-throughput microarrays and other gene expression techniques have made significant advancements. The number of identified non-coding RNAs is rapidly increasing and their exact regulatory functions are being explored. A publishable draft (build 10) of the swine genome sequence was available for the pig genomics community by the end of December 2010. Build 9 of the porcine genome is currently available with Ensembl annotation; manual annotation is ongoing. These drafts provide useful tools for such endeavors as comparative genomics and SNP scans for fine QTL mapping. A recent community-wide effort to create a 60K porcine SNP chip has greatly facilitated whole-genome association analyses, haplotype block construction and linkage disequilibrium mapping, which can contribute to whole-genome selection. The future 'systems biology' that integrates and optimizes the information from all research levels can enhance the pig community's understanding of the full complexity of the porcine genome. These recent technological advances and where they may lead are reviewed.

Identification of a second gene associated with variation in vertebral number in domestic pigs

Background

The number of vertebrae in pigs varies and is associated with body size. Wild boars have 19 vertebrae, but European commercial breeds for pork production have 20 to 23 vertebrae. We previously identified two quantitative trait loci (QTLs) for number of vertebrae on Sus scrofa chromosomes (SSC) 1 and 7, and reported that an orphan nuclear receptor, NR6A1, was located at the QTL on SSC1. At the NR6A1 locus, wild boars and Asian local breed pigs had the wild-type allele and European commercial-breed pigs had an allele associated with increased numbers of vertebrae (number-increase allele).

Results

Here, we performed a map-based study to define the other QTL, on SSC7, for which we detected genetic diversity in European commercial breeds. Haplotype analysis with microsatellite markers revealed a 41-kb conserved region within all the number-increase alleles in the present study. We also developed single nucleotide polymorphisms (SNPs) in the 450-kb region around the QTL and used them for a linkage disequilibrium analysis and an association study in 199 independent animals. Three haplotype blocks were detected, and SNPs in the 41-kb region presented the highest associations with the number of vertebrae. This region encodes an uncharacterized hypothetical protein that is not a member of any other known gene family. Orthologs appear to exist not only in mammals but also birds and fish. This gene, which we have named vertnin (VRTN) is a candidate for the gene associated with variation in vertebral number. In pigs, the number-increase allele was expressed more abundantly than the wild-type allele in embryos. Among candidate polymorphisms, there is an insertion of a SINE element (PRE1) into the intron of the Q allele as well as the SNPs in the promoter region.

Conclusions

Genetic diversity of VRTN is the suspected cause of the heterogeneity of the number of vertebrae in commercial-breed pigs, so the polymorphism information should be directly useful for assessing the genetic ability of individual animals. The number-increase allele of swine VRTN was suggested to add an additional thoracic segment to the animal. Functional analysis of VRTN may provide novel findings in the areas of developmental biology.

Pig genome sequence--analysis and publication strategy

Background

The pig genome is being sequenced and characterised under the auspices of the Swine Genome Sequencing Consortium. The sequencing strategy followed a hybrid approach combining hierarchical shotgun sequencing of BAC clones and whole genome shotgun sequencing.

Results

Assemblies of the BAC clone derived genome sequence have been annotated using the Pre-Ensembl and Ensembl automated pipelines and made accessible through the Pre-Ensembl/Ensembl browsers. The current annotated genome assembly (Sscrofa9) was released with Ensembl 56 in September 2009. A revised assembly (Sscrofa10) is under construction and will incorporate whole genome shotgun sequence (WGS) data providing > 30× genome coverage. The WGS sequence, most of which comprise short Illumina/Solexa reads, were generated from DNA from the same single Duroc sow as the source of the BAC library from which clones were preferentially selected for sequencing. In accordance with the Bermuda and Fort Lauderdale agreements and the more recent Toronto Statement the data have been released into public sequence repositories (Genbank/EMBL, NCBI/Ensembl trace repositories) in a timely manner and in advance of publication.

Conclusions

In this marker paper, the Swine Genome Sequencing Consortium (SGSC) sets outs its plans for analysis of the pig genome sequence, for the application and publication of the results.

MicroRNAome of porcine pre- and postnatal development

The domestic pig is of enormous agricultural significance and valuable models for many human diseases. Information concerning the pig microRNAome (miRNAome) has been long overdue and elucidation of this information will permit an atlas of microRNA (miRNA) regulation functions and networks to be constructed. Here we performed a comprehensive search for porcine miRNAs on ten small RNA sequencing libraries prepared from a mixture of tissues obtained during the entire pig lifetime, from the fetal period through adulthood. The sequencing results were analyzed using mammalian miRNAs, the precursor hairpins (pre-miRNAs) and the first release of the high-coverage porcine genome assembly (Sscrofa9, April 2009) and the available expressed sequence tag (EST) sequences. Our results extend the repertoire of pig miRNAome to 867 pre-miRNAs (623 with genomic coordinates) encoding for 1,004 miRNAs, of which 777 are unique. We preformed real-time quantitative PCR (q-PCR) experiments for selected 30 miRNAs in 47 tissue-specific samples and found agreement between the sequencing and q-PCR data. This broad survey provides detailed information about multiple variants of mature sequences, precursors, chromosomal organization, development-specific expression, and conservation patterns. Our data mining produced a broad view of the pig miRNAome, consisting of miRNAs and isomiRs and a wealth of information of pig miRNA characteristics. These results are prelude to the advancement in pig biology as well the use of pigs as model organism for human biological and biomedical studies.

Fine mapping of quantitative trait loci for meat color on Sus scrofa chromosome 6: analysis of the swine NUDT7 gene

In the livestock industry, meat color has become important because consumer acceptance is subject to the appearance of the product in the marketplace. Our previous analyses of a whole genome QTL scan for various meat qualities using 2 F(2) families from Japanese wild boar (known as a red meat) x Large White and from Duroc x Chinese Jinhua suggested that a meat color (heme content) QTL is located on SSC6. The objective of this study was to fine-map this SSC6 meat color QTL and subsequently investigate positional candidate genes for polymorphisms that may cause changes in meat color. Therefore, we conducted interval mapping on SSC6 using an additional 9 gene markers through combined analyses of the 2 F(2) families of Japanese wild boar x Large White (353 progeny) and Duroc x Chinese Jinhua (204 progeny). Comparative analysis with humans, mice, and cattle suggested that there were 10 functional genes in the region. Among these genes, we suggested that a novel pig gene encoding a nudix (nucleoside diphosphate linked moiety X)-type motif 7 (NUDT7, a member of the nudix hydrolases) is a strong candidate for the QTL because the mouse Nudt7 is reported to hydrolyze succinyl-CoA, a substrate of the reaction limiting the rate of heme biosynthesis. We therefore determined the pig NUDT7 gene sequence including the 5' promoter region and explored genetic polymorphisms between Japanese wild boar and Large White. We identified 116 polymorphisms within the NUDT7 CDS or in the 5' region. None of the AA substitutions were associated with the meat color QTL; however, 3 polymorphisms were found in putative transcription factor recognition sites. We then investigated the differential expression of NUDT7 in Japanese wild boar and Large White by allele-specific quantitative real-time PCR. The expression level of the Large White type allele was greater than that of the Japanese wild-boar-type allele. Consequently, we speculated that the difference in meat color between Japanese wild boar and Large White is caused partly by differential expression of this candidate gene. Upregulation of NUDT7 expression in muscle may reduce succinyl-CoA content and thus reduce the level of heme biosynthesis.

Application of massive parallel sequencing to whole genome SNP discovery in the porcine genome

BACKGROUND

Although the Illumina 1 G Genome Analyzer generates billions of base pairs of sequence data, challenges arise in sequence selection due to the varying sequence quality. Therefore, in the framework of the International Porcine SNP Chip Consortium, this pilot study aimed to evaluate the impact of the quality level of the sequenced bases on mapping quality and identification of true SNPs on a large scale.

RESULTS

DNA pooled from five animals from a commercial boar line was digested with DraI; 150-250-bp fragments were isolated and end-sequenced using the Illumina 1 G Genome Analyzer, yielding 70,348,064 sequences 36-bp long. Rules were developed to select sequences, which were then aligned to unique positions in a reference genome. Sequences were selected based on quality, and three thresholds of sequence quality (SQ) were compared. The highest threshold of SQ allowed identification of a larger number of SNPs (17,489), distributed widely across the pig genome. In total, 3,142 SNPs were validated with a success rate of 96%. The correlation between estimated minor allele frequency (MAF) and genotyped MAF was moderate, and SNPs were highly polymorphic in other pig breeds. Lowering the SQ threshold and maintaining the same criteria for SNP identification resulted in the discovery of fewer SNPs (16,768), of which 259 were not identified using higher SQ levels. Validation of SNPs found exclusively in the lower SQ threshold had a success rate of 94% and a low correlation between estimated MAF and genotyped MAF. Base change analysis suggested that the rate of transitions in the pig genome is likely to be similar to that observed in humans. Chromosome X showed reduced nucleotide diversity relative to autosomes, as observed for other species.

CONCLUSION

Large numbers of SNPs can be identified reliably by creating strict rules for sequence selection, which simultaneously decreases sequence ambiguity. Selection of sequences using a higher SQ threshold leads to more reliable identification of SNPs. Lower SQ thresholds can be used to guarantee sufficient sequence coverage, resulting in high success rate but less reliable MAF estimation. Nucleotide diversity varies between porcine chromosomes, with the X chromosome showing less variation as observed in other species.