Swine global genomic resources: insights into wild and domesticated populations.

Suids, both domesticated and wild, are found on all continents except for Antarctica and provide valuable food resources for humans in addition to serving as important models for biomedical research. Continuing advances in genome sequencing have allowed researchers to compare the genomes from diverse populations of suids helping to clarify their evolution and dispersal. Further analysis of these samples may provide clues to improve disease resistance/resilience and productivity in domestic suids as well as better ways of classifying and conserving genetic diversity within wild and captive suids. Collecting samples from diverse populations of suids is resource intensive and may negatively impact endangered populations. Here we catalog extensive tissue and DNA samples from suids in collections in both Europe and North America. We include samples that have previously been used for whole genome sequencing, targeted DNA sequencing, RNA sequencing, and reduced representation bisulfite sequencing (RRBS). This work provides an important centralized resource for researchers who wish to access published databases.

Swine models for translational oncological research: an evolving landscape and regulatory considerations.

Swine biomedical models have been gaining in popularity over the last decade, particularly for applications in oncology research. Swine models for cancer research include pigs that have severe combined immunodeficiency for xenotransplantation studies, genetically modified swine models which are capable of developing tumors in vivo, as well as normal immunocompetent pigs. In recent years, there has been a low success rate for the approval of new oncological therapeutics in clinical trials. The two leading reasons for these failures are either due to toxicity and safety issues or lack of efficacy. As all therapeutics must be tested within animal models prior to clinical testing, there are opportunities to expand the ability to assess efficacy and toxicity profiles within the preclinical testing phases of new therapeutics. Most preclinical in vivo testing is performed in mice, canines, and non-human primates. However, swine models are an alternative large animal model for cancer research with similarity to human size, genetics, and physiology. Additionally, tumorigenesis pathways are similar between human and pigs in that similar driver mutations are required for transformation. Due to their larger size, the development of orthotopic tumors is easier than in smaller rodent models; additionally, porcine models can be harnessed for testing of new interventional devices and radiological/surgical approaches as well. Taken together, swine are a feasible option for preclinical therapeutic and device testing. The goals of this resource are to provide a broad overview on regulatory processes required for new therapeutics and devices for use in the clinic, cross-species differences in oncological therapeutic responses, as well as to provide an overview of swine oncology models that have been developed that could be used for preclinical testing to fulfill regulatory requirements.

Future of biomedical, agricultural, and biological systems research using domesticated animals.

Increased knowledge of reproduction and health of domesticated animals is integral to sustain and improve global competitiveness of U.S. animal agriculture, understand and resolve complex animal and human diseases, and advance fundamental research in sciences that are critical to understanding mechanisms of action and identifying future targets for interventions. Historically, federal and state budgets have dwindled and funding for the United States Department of Agriculture (USDA) National Institute of Food and Agriculture (NIFA) competitive grants programs remained relatively stagnant from 1985 through 2010. This shortage in critical financial support for basic and applied research, coupled with the underappreciated knowledge of the utility of non-rodent species for biomedical research, hindered funding opportunities for research involving livestock and limited improvements in both animal agriculture and animal and human health. In 2010, the National Institutes of Health and USDA NIFA established an interagency partnership to promote the use of agriculturally important animal species in basic and translational research relevant to both biomedicine and agriculture. This interagency program supported 61 grants totaling over $107 million with 23 awards to new or early-stage investigators. This article will review the success of the 9-year Dual Purpose effort and highlight opportunities for utilizing domesticated agricultural animals in research.

Transarterial Embolization of Liver Cancer in a Transgenic Pig Model.

PURPOSE: To develop and characterize a porcine model of liver cancer that could be used to test new locoregional therapies. MATERIALS AND METHODS: Liver tumors were induced in 18 Oncopigs (transgenic pigs with Cre-inducible TP53R167H and KRASG12D mutations) by using an adenoviral vector encoding the Cre-recombinase gene. The resulting 60 tumors were characterized on multiphase contrast-enhanced CT, angiography, perfusion, micro-CT, and necropsy. Transarterial embolization was performed using 40-120 µm (4 pigs) or 100-300 µm (4 pigs) Embosphere microspheres. Response to embolization was evaluated on imaging. Complications were determined based on daily clinical evaluation, laboratory results, imaging, and necropsy. RESULTS: Liver tumors developed at 60/70 (86%) inoculated sites. Mean tumor size was 2.1 cm (range, 0.3-4 cm) at 1 week. Microscopically, all animals developed poorly differentiated to undifferentiated carcinomas accompanied by a major inflammatory component, which resembled undifferentiated carcinomas of the human pancreatobiliary tract. Cytokeratin and vimentin expression confirmed epithelioid and mesenchymal differentiation, respectively. Lymph node, lung, and peritoneal metastases were seen in some cases. On multiphase CT, all tumors had a hypovascular center, and 17/60 (28%) had a hypervascular rim. After transarterial embolization, noncontrast CT showed retained contrast medium in the tumors. Follow-up contrast-enhanced scan showed reduced size of tumors after embolization using either 40-120 µm or 100-300 µm Embosphere microspheres, while untreated tumors showed continued growth. CONCLUSIONS: Liver tumors can be induced in a transgenic pig and can be successfully treated using bland embolization.

Correction to: Genome-wide SNP data unveils the globalization of domesticated pigs.

An amendment to this paper has been published and can be accessed via the original article.

Characterization of an Inducible Alcoholic Liver Fibrosis Model for Hepatocellular Carcinoma Investigation in a Transgenic Porcine Tumorigenic Platform.

PURPOSE: This study used the Oncopig Cancer Model (OCM) to develop alcohol-induced fibrosis in a porcine model capable of developing hepatocellular carcinoma. MATERIALS AND METHODS: Liver injury was induced in 8-week-old Oncopigs (n = 10) via hepatic transarterial infusion of 0.75 mL/kg ethanol-ethiodized oil (1:3 v/v). Feasibility was assessed in an initial Oncopig cohort (n = 5) by histologic analysis at 8 weeks after induction, and METAVIR results were compared to age- and sex-matched healthy controls (n = 5). Liver injury was then induced in a second OCM cohort (n = 5) for a time-course study, with post-induction disease surveillance via biweekly physical exam, lab analysis, and liver biopsies until 20 weeks after induction. RESULTS: In Cohort 1, 8-week post-induction liver histologic analysis revealed median METAVIR F3 (range, F3-F4) fibrosis, A2 (range, A2-A3) inflammation, and 15.3% (range, 5.0%-22.9%) fibrosis. METAVIR and inflammation scores were generally elevated compared to healthy controls (F0-F1, P = 0.0013; A0-A1, P = .0013; median percent fibrosis 8.7%, range, 5.8%-12.1%, P = .064). In Cohort 2, histologic analysis revealed peak fibrosis severity of median METAVIR F3 (range, F2-F3). However, lack of persistent alcohol exposure resulted in liver recovery, with median METAVIR F2 (range, F1-F2) fibrosis at 20 weeks after induction. No behavioral or biochemical abnormalities were observed to indicate liver decompensation. CONCLUSIONS: This study successfully validated a protocol to develop METAVIR F3-F4 fibrosis within 8 weeks in the OCM, supporting its potential to serve as a model for hepatocellular carcinoma in a fibrotic liver background. Further investigation is required to determine if repeated alcohol liver injury is required to develop an irreversible METAVIR grade F4 porcine cirrhosis model.

Genome-wide SNP data unveils the globalization of domesticated pigs.

BACKGROUND: Pigs were domesticated independently in Eastern and Western Eurasia early during the agricultural revolution, and have since been transported and traded across the globe. Here, we present a worldwide survey on 60K genome-wide single nucleotide polymorphism (SNP) data for 2093 pigs, including 1839 domestic pigs representing 122 local and commercial breeds, 215 wild boars, and 39 out-group suids, from Asia, Europe, America, Oceania and Africa. The aim of this study was to infer global patterns in pig domestication and diversity related to demography, migration, and selection. RESULTS: A deep phylogeographic division reflects the dichotomy between early domestication centers. In the core Eastern and Western domestication regions, Chinese pigs show differentiation between breeds due to geographic isolation, whereas this is less pronounced in European pigs. The inferred European origin of pigs in the Americas, Africa, and Australia reflects European expansion during the sixteenth to nineteenth centuries. Human-mediated introgression, which is due, in particular, to importing Chinese pigs into the UK during the eighteenth and nineteenth centuries, played an important role in the formation of modern pig breeds. Inbreeding levels vary markedly between populations, from almost no runs of homozygosity (ROH) in a number of Asian wild boar populations, to up to 20% of the genome covered by ROH in a number of Southern European breeds. Commercial populations show moderate ROH statistics. For domesticated pigs and wild boars in Asia and Europe, we identified highly differentiated loci that include candidate genes related to muscle and body development, central nervous system, reproduction, and energy balance, which are putatively under artificial selection. CONCLUSIONS: Key events related to domestication, dispersal, and mixing of pigs from different regions are reflected in the 60K SNP data, including the globalization that has recently become full circle since Chinese pig breeders in the past decades started selecting Western breeds to improve local Chinese pigs. Furthermore, signatures of ongoing and past selection, acting at different times and on different genetic backgrounds, enhance our insight in the mechanism of domestication and selection. The global diversity statistics presented here highlight concerns for maintaining agrodiversity, but also provide a necessary framework for directing genetic conservation.

Evolutionary patterns of Toll-like receptor signaling pathway genes in the Suidae.

BACKGROUND: The Toll-like receptor (TLR) signaling pathway constitutes an essential component of the innate immune system. Highly conserved proteins, indicative of their critical roles in host survival, characterize this pathway. Selective constraints could vary depending on the gene's position within the pathway as TLR signaling is a sequential process and that genes downstream of the TLRs may be more selectively constrained to ensure efficient immune responses given the important role of downstream genes in the signaling process. Thus, we investigated whether gene position influenced protein evolution in the TLR signaling pathway of the Suidae. The members of the Suidae examined included the European Sus scrofa (wild boar), Asian Sus scrofa (wild boar), Sus verrucosus, Sus celebensis, Sus scebifrons, Sus barbatus, Babyrousa babyrussa, Potamochoerus larvatus, Potamochoerus porcus and Phacochoerus africanus. RESULTS: A total of 33 TLR signaling pathway genes in the Suidae were retrieved from resequencing data. The evolutionary parameter ω (dn/ds) had an overall mean of 0.1668 across genes, indicating high functional conservation within the TLR signaling pathway. A significant relationship was inferred for the network parameters gene position, number of protein-protein interactions, protein length and the evolutionary parameter dn (nonsynonymous substitutions) such that downstream genes had lower nonsynonymous substitution rates, more interactors and shorter protein length than upstream genes. Gene position was significantly correlated with the number of protein-protein interactions and protein length. Thus, the polarity in the selective constraint along the TLR signaling pathway was due to the number of molecules a protein interacted with and the protein's length. CONCLUSION: Results indicate that the level of selective constraints on genes within the TLR signaling pathway of the Suidae is dependent on the gene's position and network parameters. In particular, downstream genes evolve more slowly as a result of being highly connected and having shorter protein lengths. These findings highlight the critical role of gene network parameters in gene evolution.

Impact of neonatal iron deficiency on hippocampal DNA methylation and gene transcription in a porcine biomedical model of cognitive development.

BACKGROUND: Iron deficiency is a common childhood micronutrient deficiency that results in altered hippocampal function and cognitive disorders. However, little is known about the mechanisms through which neonatal iron deficiency results in long lasting alterations in hippocampal gene expression and function. DNA methylation is an epigenetic mark involved in gene regulation and altered by environmental factors. In this study, hippocampal DNA methylation and gene expression were assessed via reduced representation bisulfite sequencing and RNA-seq on samples from a previous study reporting reduced hippocampal-based learning and memory in a porcine biomedical model of neonatal iron deficiency. RESULTS: In total 192 differentially expressed genes (DEGs) were identified between the iron deficient and control groups. GO term and pathway enrichment analysis identified DEGs associated with hypoxia, angiogenesis, increased blood brain barrier (BBB) permeability, and altered neurodevelopment and function. Of particular interest are genes previously implicated in cognitive deficits and behavioral disorders in humans and mice, including HTR2A, HTR2C, PAK3, PRSS12, and NETO1. Altered genome-wide DNA methylation was observed across 0.5 million CpG and 2.4 million non-CpG sites. In total 853 differentially methylated (DM) CpG and 99 DM non-CpG sites were identified between groups. Samples clustered by group when comparing DM non-CpG sites, suggesting high conservation of non-CpG methylation in response to neonatal environment. In total 12 DM sites were associated with 9 DEGs, including genes involved in angiogenesis, neurodevelopment, and neuronal function. CONCLUSIONS: Neonatal iron deficiency leads to altered hippocampal DNA methylation and gene regulation involved in hypoxia, angiogenesis, increased BBB permeability, and altered neurodevelopment and function. Together, these results provide new insights into the mechanisms through which neonatal iron deficiency results in long lasting reductions in cognitive development in humans.

Evidence for adaptation of porcine Toll-like receptors.

Naturally endemic infectious diseases provide selective pressures for pig populations. Toll-like receptors (TLRs) represent the first line of immune defense against pathogens and are likely to play a crucial adaptive role for pig populations. This study was done to determine whether wild and domestic pig populations representing diverse global environments demonstrate local TLR adaptation. The genomic sequence encoding the ectodomain, responsible for interacting with pathogen ligands of bacterial (TLR1, TLR2 and TLR6) and viral (TLR3, TLR7 and TLR8) receptors, was obtained. Mitochondrial D-loop region sequences were obtained and a phylogenetic analysis using these sequences revealed a clear separation of animals into Asian (n = 27) and European (n = 40) clades. The TLR sequences were then analyzed for population-specific positive selection signatures within wild boars and domesticated pig populations derived from Asian and European clades. Using within-population and between-population tests for positive selection, a TLR2-derived variant 376A (126Thr), estimated to have arisen in 163,000 years ago with a frequency of 83.33% within European wild boars, 98.00% within domestic pig breeds of European origin, 40.00% within Asian wild boars, and 11.36% within Asian domestic pigs, was identified to be under positive selection in pigs of European origin. The variant is located within the N terminal domain of the TLR2 protein 3D crystal structure and could affect ligand binding. This study suggests the TLR2 gene contributing to responses to bacterial pathogens has been crucial in adaptation of pigs to pathogens.

Peripheral viral infection induced microglial sensome genes and enhanced microglial cell activity in the hippocampus of neonatal piglets.

Although poorly understood, early-life infection is predicted to affect brain microglial cells, making them hypersensitive to subsequent stimuli. To investigate this, we assessed gene expression in hippocampal tissue obtained from a previously published study reporting increased microglial cell activity and reduced hippocampal-dependent learning in neonatal piglets infected with porcine reproductive and respiratory syndrome virus (PRRSV), a virus that induces interstitial pneumonia. Infection altered expression of 455 genes, of which 334 were up-regulated and 121 were down-regulated. Functional annotation revealed that immune function genes were enriched among the up-regulated differentially expressed genes (DEGs), whereas calcium binding and synaptic vesicle genes were enriched among the down-regulated DEGs. Twenty-six genes encoding part of the microglia sensory apparatus (i.e., the sensome) were up-regulated (e.g., IL1R1, TLR2, and TLR4), whereas 15 genes associated with the synaptosome and synaptic receptors (e.g., NPTX2, GABRA2, and SLC5A7) were down-regulated. As the sensome may foretell microglia reactivity, we next inoculated piglets with culture medium or PRRSV at PD 7 and assessed hippocampal microglia morphology and function at PD 28 when signs of infection were waning. Consistent with amplification of the sensome, microglia from PRRSV piglets had enhanced responsiveness to chemoattractants, increased phagocytic activity, and secreted more TNFα in response to lipopolysaccharide and Poly I:C. Immunohistochemical staining indicated PRRSV infection increased microglia soma length and length-to-width ratio. Bipolar rod-like microglia not evident in hippocampus of control piglets, were present in infected piglets. Collectively, this study suggests early-life infection alters the microglia sensome as well as microglial cell morphology and function.

Adult porcine genome-wide DNA methylation patterns support pigs as a biomedical model.

BACKGROUND: Pigs (Sus scrofa) provide relevant biomedical models to dissect complex diseases due to their anatomical, genetic, and physiological similarities with humans. Aberrant DNA methylation has been linked to many of these diseases and is associated with gene expression; however, the functional similarities and differences between porcine and human DNA methylation patterns are largely unknown. METHODS: DNA and RNA was isolated from eight tissue samples (fat, heart, kidney, liver, lung, lymph node, muscle, and spleen) from the adult female Duroc utilized for the pig genome sequencing project. Reduced representation bisulfite sequencing (RRBS) and RNA-seq were performed on an Illumina HiSeq2000. RRBS reads were aligned using BSseeker2, and only sites with a minimum depth of 10 reads were used for methylation analysis. RNA-seq reads were aligned using Tophat, and expression analysis was performed using Cufflinks. In addition, SNP calling was performed using GATK for targeted control and whole genome sequencing reads for CpG site validation and allelic expression analysis, respectively. RESULTS: Analysis on the influence of DNA variation in methylation calling revealed a reduced effectiveness of WGS datasets in covering CpG rich regions, as well as the usefulness of a targeted control library for SNP detection. Analysis of over 500,000 CpG sites demonstrated genome wide methylation patterns similar to those observed in humans, including reduced methylation within CpG islands and at transcription start sites (TSS), X chromosome inactivation, and anticorrelation of TSS CpG methylation with gene expression. In addition, a positive correlation between TSS CpG density and expression, and a negative correlation between TSS TpG density and expression were demonstrated. Low but non-random non-CpG methylation (<1%) was also detected in all non-neuronal somatic tissues, with differences in tissue clustering observed based on CpG and non-CpG methylation patterns. Finally, allele specific expression analysis revealed enrichment of genes involved in metabolic and regulatory processes. DISCUSSION: These results provide transcriptional and DNA methylation datasets for the biomedical community that are directly relatable to current genomic resources. In addition, the correlation between TSS CpG density and expression suggests increased mutation rates at CpG sites play a significant role in adaptive evolution by reducing CpG density at TSS over time, resulting in higher methylation levels in these regions and more permanent changes to lower gene expression. This is proposed to occur predominantly through deamination of 5-methylcytosine to thymidine, resulting in the replacement of CpG with TpG sites in these regions, as indicated by the increased TSS TpG density observed in non-expressed genes, resulting in a negative correlation between expression and TSS TpG density. CONCLUSIONS: This study provides baseline methylation and gene transcription profiles for a healthy adult pig, reports similar patterns to those observed in humans, and supports future porcine studies related to human disease and development. Additionally, the observed reduced CpG and increased TpG density at TSS of lowly expressed genes suggests DNA methylation plays a significant role in adaptive evolution through more permanent changes to lower gene expression.

Regions of homozygosity in the porcine genome: consequence of demography and the recombination landscape.

Inbreeding has long been recognized as a primary cause of fitness reduction in both wild and domesticated populations. Consanguineous matings cause inheritance of haplotypes that are identical by descent (IBD) and result in homozygous stretches along the genome of the offspring. Size and position of regions of homozygosity (ROHs) are expected to correlate with genomic features such as GC content and recombination rate, but also direction of selection. Thus, ROHs should be non-randomly distributed across the genome. Therefore, demographic history may not fully predict the effects of inbreeding. The porcine genome has a relatively heterogeneous distribution of recombination rate, making Sus scrofa an excellent model to study the influence of both recombination landscape and demography on genomic variation. This study utilizes next-generation sequencing data for the analysis of genomic ROH patterns, using a comparative sliding window approach. We present an in-depth study of genomic variation based on three different parameters: nucleotide diversity outside ROHs, the number of ROHs in the genome, and the average ROH size. We identified an abundance of ROHs in all genomes of multiple pigs from commercial breeds and wild populations from Eurasia. Size and number of ROHs are in agreement with known demography of the populations, with population bottlenecks highly increasing ROH occurrence. Nucleotide diversity outside ROHs is high in populations derived from a large ancient population, regardless of current population size. In addition, we show an unequal genomic ROH distribution, with strong correlations of ROH size and abundance with recombination rate and GC content. Global gene content does not correlate with ROH frequency, but some ROH hotspots do contain positive selected genes in commercial lines and wild populations. This study highlights the importance of the influence of demography and recombination on homozygosity in the genome to understand the effects of inbreeding.

Quantification of the relative roles of niche and neutral processes in structuring gastrointestinal microbiomes.

The theoretical description of the forces that shape ecological communities focuses around two classes of models. In niche theory, deterministic interactions between species, individuals, and the environment are considered the dominant factor, whereas in neutral theory, stochastic forces, such as demographic noise, speciation, and immigration, are dominant. Species abundance distributions predicted by the two classes of theory are difficult to distinguish empirically, making it problematic to deduce ecological dynamics from typical measures of diversity and community structure. Here, we show that the fusion of species abundance data with genome-derived measures of evolutionary distance can provide a clear indication of ecological dynamics, capable of quantifying the relative roles played by niche and neutral forces. We apply this technique to six gastrointestinal microbiomes drawn from three different domesticated vertebrates, using high-resolution surveys of microbial species abundance obtained from carefully curated deep 16S rRNA hypervariable tag sequencing data. Although the species abundance patterns are seemingly well fit by the neutral theory of metacommunity assembly, we show that this theory cannot account for the evolutionary patterns in the genomic data; moreover, our analyses strongly suggest that these microbiomes have, in fact, been assembled through processes that involve a significant nonneutral (niche) contribution. Our results demonstrate that high-resolution genomics can remove the ambiguities of process inference inherent in classic ecological measures and permits quantification of the forces shaping complex microbial communities.

Oncopig bladder cancer cells recapitulate human bladder cancer treatment responses in vitro.

Introduction: Bladder cancer is a common neoplasia of the urinary tract that holds the highest cost of lifelong treatment per patient, highlighting the need for a continuous search for new therapies for the disease. Current bladder cancer models are either imperfect in their ability to translate results to clinical practice (mouse models), or rare and not inducible (canine models). Swine models are an attractive alternative to model the disease due to their similarities with humans on several levels. The Oncopig Cancer Model has been shown to develop tumors that closely resemble human tumors. However, urothelial carcinoma has not yet been studied in this platform. Methods: We aimed to develop novel Oncopig bladder cancer cell line (BCCL) and investigate whether these urothelial swine cells mimic human bladder cancer cell line (5637 and T24) treatment-responses to cisplatin, doxorubicin, and gemcitabine in vitro. Results: Results demonstrated consistent treatment responses between Oncopig and human cells in most concentrations tested (p>0.05). Overall, Oncopig cells were more predictive of T24 than 5637 cell therapeutic responses. Microarray analysis also demonstrated similar alterations in expression of apoptotic (GADD45B and TP53INP1) and cytoskeleton-related genes (ZMYM6 and RND1) following gemcitabine exposure between 5637 (human) and Oncopig BCCL cells, indicating apoptosis may be triggered through similar signaling pathways. Molecular docking results indicated that swine and humans had similar Dg values between the chemotherapeutics and their target proteins. Discussion: Taken together, these results suggest the Oncopig could be an attractive animal to model urothelial carcinoma due to similarities in in vitro therapeutic responses compared to human cells.

Compression-based distance (CBD): a simple, rapid, and accurate method for microbiota composition comparison.

BACKGROUND: Perturbations in intestinal microbiota composition have been associated with a variety of gastrointestinal tract-related diseases. The alleviation of symptoms has been achieved using treatments that alter the gastrointestinal tract microbiota toward that of healthy individuals. Identifying differences in microbiota composition through the use of 16S rRNA gene hypervariable tag sequencing has profound health implications. Current computational methods for comparing microbial communities are usually based on multiple alignments and phylogenetic inference, making them time consuming and requiring exceptional expertise and computational resources. As sequencing data rapidly grows in size, simpler analysis methods are needed to meet the growing computational burdens of microbiota comparisons. Thus, we have developed a simple, rapid, and accurate method, independent of multiple alignments and phylogenetic inference, to support microbiota comparisons. RESULTS: We create a metric, called compression-based distance (CBD) for quantifying the degree of similarity between microbial communities. CBD uses the repetitive nature of hypervariable tag datasets and well-established compression algorithms to approximate the total information shared between two datasets. Three published microbiota datasets were used as test cases for CBD as an applicable tool. Our study revealed that CBD recaptured 100% of the statistically significant conclusions reported in the previous studies, while achieving a decrease in computational time required when compared to similar tools without expert user intervention. CONCLUSION: CBD provides a simple, rapid, and accurate method for assessing distances between gastrointestinal tract microbiota 16S hypervariable tag datasets.

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.

A high density recombination map of the pig reveals a correlation between sex-specific recombination and GC content.

BACKGROUND: The availability of a high-density SNP genotyping chip and a reference genome sequence of the pig (Sus scrofa) enabled the construction of a high-density linkage map. A high-density linkage map is an essential tool for further fine-mapping of quantitative trait loci (QTL) for a variety of traits in the pig and for a better understanding of mechanisms underlying genome evolution. RESULTS: Four different pig pedigrees were genotyped using the Illumina PorcineSNP60 BeadChip. Recombination maps for the autosomes were computed for each individual pedigree using a common set of markers. The resulting genetic maps comprised 38,599 SNPs, including 928 SNPs not positioned on a chromosome in the current assembly of the pig genome (build 10.2). The total genetic length varied according to the pedigree, from 1797 to 2149 cM. Female maps were longer than male maps, with a notable exception for SSC1 where male maps are characterized by a higher recombination rate than females in the region between 91-250 Mb. The recombination rates varied among chromosomes and along individual chromosomes, regions with high recombination rates tending to cluster close to the chromosome ends, irrespective of the position of the centromere. Correlations between main sequence features and recombination rates were investigated and significant correlations were obtained for all the studied motifs. Regions characterized by high recombination rates were enriched for specific GC-rich sequence motifs as compared to low recombinant regions. These correlations were higher in females than in males, and females were found to be more recombinant than males at regions where the GC content was greater than 0.4. CONCLUSIONS: The analysis of the recombination rate along the pig genome highlighted that the regions exhibiting higher levels of recombination tend to cluster around the ends of the chromosomes irrespective of the location of the centromere. Major sex-differences in recombination were observed: females had a higher recombination rate within GC-rich regions and exhibited a stronger correlation between recombination rates and specific sequence features.

Association of the porcine transforming growth factor beta type I receptor (TGFBR1) gene with growth and carcass traits.

BACKGROUND: Growth and carcass traits are of great economic importance in livestock production. A large number of quantitative trait loci (QTL) have been identified for growth and carcass traits on porcine chromosome one (SSC1). A key positional candidate for this chromosomal region is TGFBR1 (transforming growth factor beta type I receptor). This gene plays a key role in inherited disorders at cardiovascular, craniofacial, neurocognitive, and skeletal development in mammals. RESULTS: In this study, 27 polymorphic SNPs in the porcine TGFBR1 gene were identified on the University of Illinois Yorkshire × Meishan resource population. Three SNPs (SNP3, SNP43, SNP64) representing major polymorphic patterns of the 27 SNPs in F1 and F0 individuals of the Illinois population were selected for analyses of QTL association and genetic diversity. An association analysis for growth and carcass traits was completed using these three representative SNPs in the Illinois population with 298 F2 individuals and a large commercial population of 1008 animals. The results indicate that the TGFBR1 gene polymorphism (SNP64) is significantly associated (p < 0.05) with growth rates including average daily gains between birth and 56 kg (p = 0.049), between 5.5 and 56 kg (p = 0.024), between 35 and 56 kg (p = 0.021). Significant associations (p < 0.05) were also identified between TGFBR1 gene polymorphisms (SNP3/SNP43) and carcass traits including loin-eye-area (p = 0.022) in the Illinois population, and back-fat thickness (p = 0.0009), lean percentage (p = 0.0023) and muscle color (p = 0.021) in the commercial population. These three SNPs were also used to genotype a diverse panel of 130 animals representing 11 pig breeds. Alleles SNP3_T and SNP43_G were fixed in Pietrain and Sinclair pig breeds. SNP64_G allele was uniquely identified in Chinese Meishan pigs. Strong evidence of association (p < 0.01) between both SNP3 and SNP64 alleles and reproductive traits including gestation length and number of corpora lutea were also observed in the Illinois population. CONCLUSION: This study gives the first evidence of association between the porcine TGFBR1 gene and traits of economic importance and provides support for using TGFBR1 markers for pig breeding and selection programs. The genetic diversities in different pig breeds would be helpful to understand the genetic background and migration of the porcine TGFBR1 gene.