Investigative power of Genomic Informational Field Theory (GIFT) relative to GWAS for genotype-phenotype mapping.

Identifying associations between phenotype and genotype is the fundamental basis of genetic analyses. Inspired by frequentist probability and the work of R.A. Fisher, genome-wide association studies (GWAS) extract information using averages and variances from genotype-phenotype datasets. Averages and variances are legitimated upon creating distribution density functions obtained through the grouping of data into categories. However, as data from within a given category cannot be differentiated, the investigative power of such methodologies is limited. Genomic Informational Field Theory (GIFT) is a method specifically designed to circumvent this issue. The way GIFT proceeds is opposite to that of GWAS. Whilst GWAS determines the extent to which genes are involved in phenotype formation (bottom-up approach), GIFT determines the degree to which the phenotype can select microstates (genes) for its subsistence (top-down approach). Doing so requires dealing with new genetic concepts, a.k.a. genetic paths, upon which significance levels for genotype-phenotype associations can be determined. By using different datasets obtained in ovis aries related to bone growth (Dataset-1) and to a series of linked metabolic and epigenetic pathways (Dataset-2), we demonstrate that removing the informational barrier linked to categories enhances the investigative and discriminative powers of GIFT, namely that GIFT extracts more information than GWAS. We conclude by suggesting that GIFT is an adequate tool to study how phenotypic plasticity and genetic assimilation are linked.

nf-core/isoseq: simple gene and isoform annotation with PacBio Iso-Seq long-read sequencing.

MOTIVATION: Iso-Seq RNA long-read sequencing enables the identification of full-length transcripts and isoforms, removing the need for complex analysis such as transcriptome assembly. However, the raw sequencing data need to be processed in a series of steps before annotation is complete. Here, we present nf-core/isoseq, a pipeline for automatic read processing and genome annotation. Following nf-core guidelines, the pipeline has few dependencies and can be run on any of platforms. AVAILABILITY AND IMPLEMENTATION: The pipeline is freely available online on the nf-core website (https://nf-co.re/isoseq) and on GitHub (https://github.com/nf-core/isoseq) under MIT License (DOI: 10.5281/zenodo.7116979).

Stem cell-derived porcine macrophages as a new platform for studying host-pathogen interactions.

BACKGROUND: Infectious diseases of farmed and wild animals pose a recurrent threat to food security and human health. The macrophage, a key component of the innate immune system, is the first line of defence against many infectious agents and plays a major role in shaping the adaptive immune response. However, this phagocyte is a target and host for many pathogens. Understanding the molecular basis of interactions between macrophages and pathogens is therefore crucial for the development of effective strategies to combat important infectious diseases. RESULTS: We explored how porcine pluripotent stem cells (PSCs) can provide a limitless in vitro supply of genetically and experimentally tractable macrophages. Porcine PSC-derived macrophages (PSCdMs) exhibited molecular and functional characteristics of ex vivo primary macrophages and were productively infected by pig pathogens, including porcine reproductive and respiratory syndrome virus (PRRSV) and African swine fever virus (ASFV), two of the most economically important and devastating viruses in pig farming. Moreover, porcine PSCdMs were readily amenable to genetic modification by CRISPR/Cas9 gene editing applied either in parental stem cells or directly in the macrophages by lentiviral vector transduction. CONCLUSIONS: We show that porcine PSCdMs exhibit key macrophage characteristics, including infection by a range of commercially relevant pig pathogens. In addition, genetic engineering of PSCs and PSCdMs affords new opportunities for functional analysis of macrophage biology in an important livestock species. PSCs and differentiated derivatives should therefore represent a useful and ethical experimental platform to investigate the genetic and molecular basis of host-pathogen interactions in pigs, and also have wider applications in livestock.

Balancing selection at a premature stop mutation in the myostatin gene underlies a recessive leg weakness syndrome in pigs.

Balancing selection provides a plausible explanation for the maintenance of deleterious alleles at moderate frequency in livestock, including lethal recessives exhibiting heterozygous advantage in carriers. In the current study, a leg weakness syndrome causing mortality of piglets in a commercial line showed monogenic recessive inheritance, and a region on chromosome 15 associated with the syndrome was identified by homozygosity mapping. Whole genome resequencing of cases and controls identified a mutation causing a premature stop codon within exon 3 of the porcine Myostatin (MSTN) gene, similar to those causing a double-muscling phenotype observed in several mammalian species. The MSTN mutation was in Hardy-Weinberg equilibrium in the population at birth, but significantly distorted amongst animals still in the herd at 110 kg, due to an absence of homozygous mutant genotypes. In heterozygous form, the MSTN mutation was associated with a major increase in muscle depth and decrease in fat depth, suggesting that the deleterious allele was maintained at moderate frequency due to heterozygous advantage (allele frequency, q = 0.22). Knockout of the porcine MSTN by gene editing has previously been linked to problems of low piglet survival and lameness. This MSTN mutation is an example of putative balancing selection in livestock, providing a plausible explanation for the lack of disrupting MSTN mutations in pigs despite many generations of selection for lean growth.

Illuminating the dark side of the human transcriptome with long read transcript sequencing.

BACKGROUND: The human transcriptome annotation is regarded as one of the most complete of any eukaryotic species. However, limitations in sequencing technologies have biased the annotation toward multi-exonic protein coding genes. Accurate high-throughput long read transcript sequencing can now provide additional evidence for rare transcripts and genes such as mono-exonic and non-coding genes that were previously either undetectable or impossible to differentiate from sequencing noise. RESULTS: We developed the Transcriptome Annotation by Modular Algorithms (TAMA) software to leverage the power of long read transcript sequencing and address the issues with current data processing pipelines. TAMA achieved high sensitivity and precision for gene and transcript model predictions in both reference guided and unguided approaches in our benchmark tests using simulated Pacific Biosciences (PacBio) and Nanopore sequencing data and real PacBio datasets. By analyzing PacBio Sequel II Iso-Seq sequencing data of the Universal Human Reference RNA (UHRR) using TAMA and other commonly used tools, we found that the convention of using alignment identity to measure error correction performance does not reflect actual gain in accuracy of predicted transcript models. In addition, inter-read error correction can cause major changes to read mapping, resulting in potentially over 6 K erroneous gene model predictions in the Iso-Seq based human genome annotation. Using TAMA's genome assembly based error correction and gene feature evidence, we predicted 2566 putative novel non-coding genes and 1557 putative novel protein coding gene models. CONCLUSIONS: Long read transcript sequencing data has the power to identify novel genes within the highly annotated human genome. The use of parameter tuning and extensive output information of the TAMA software package allows for in depth exploration of eukaryotic transcriptomes. We have found long read data based evidence for thousands of unannotated genes within the human genome. More development in sequencing library preparation and data processing are required for differentiating sequencing noise from real genes in long read RNA sequencing data.

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.

A high resolution atlas of gene expression in the domestic sheep (Ovis aries).

Sheep are a key source of meat, milk and fibre for the global livestock sector, and an important biomedical model. Global analysis of gene expression across multiple tissues has aided genome annotation and supported functional annotation of mammalian genes. We present a large-scale RNA-Seq dataset representing all the major organ systems from adult sheep and from several juvenile, neonatal and prenatal developmental time points. The Ovis aries reference genome (Oar v3.1) includes 27,504 genes (20,921 protein coding), of which 25,350 (19,921 protein coding) had detectable expression in at least one tissue in the sheep gene expression atlas dataset. Network-based cluster analysis of this dataset grouped genes according to their expression pattern. The principle of 'guilt by association' was used to infer the function of uncharacterised genes from their co-expression with genes of known function. We describe the overall transcriptional signatures present in the sheep gene expression atlas and assign those signatures, where possible, to specific cell populations or pathways. The findings are related to innate immunity by focusing on clusters with an immune signature, and to the advantages of cross-breeding by examining the patterns of genes exhibiting the greatest expression differences between purebred and crossbred animals. This high-resolution gene expression atlas for sheep is, to our knowledge, the largest transcriptomic dataset from any livestock species to date. It provides a resource to improve the annotation of the current reference genome for sheep, presenting a model transcriptome for ruminants and insight into gene, cell and tissue function at multiple developmental stages.

Pigs Lacking the Scavenger Receptor Cysteine-Rich Domain 5 of CD163 Are Resistant to Porcine Reproductive and Respiratory Syndrome Virus 1 Infection.

Porcine reproductive and respiratory syndrome virus (PRRSV) has a narrow host cell tropism, limited to cells of the monocyte/macrophage lineage. CD163 protein is expressed at high levels on the surface of specific macrophage types, and a soluble form is circulating in blood. CD163 has been described as a fusion receptor for PRRSV, with the scavenger receptor cysteine-rich domain 5 (SRCR5) region having been shown to be the interaction site for the virus. As reported previously, we have generated pigs in which exon 7 of the CD163 gene has been deleted using CRISPR/Cas9 editing in pig zygotes. These pigs express CD163 protein lacking SRCR5 (ΔSRCR5 CD163) and show no adverse effects when maintained under standard husbandry conditions. Not only was ΔSRCR5 CD163 detected on the surface of macrophage subsets, but the secreted, soluble protein can also be detected in the serum of the edited pigs, as shown here by a porcine soluble CD163-specific enzyme-linked immunosorbent assay (ELISA). Previous results showed that primary macrophage cells from ΔSRCR5 CD163 animals are resistant to PRRSV-1 subtype 1, 2, and 3 as well as PRRSV-2 infection in vitro Here, ΔSRCR5 pigs were challenged with a highly virulent PRRSV-1 subtype 2 strain. In contrast to the wild-type control group, ΔSRCR5 pigs showed no signs of infection and no viremia or antibody response indicative of a productive infection. Histopathological analysis of lung and lymph node tissue showed no presence of virus-replicating cells in either tissue. This shows that ΔSRCR5 pigs are fully resistant to infection by the virus.IMPORTANCE Porcine reproductive and respiratory syndrome (PRRS) virus (PRRSV) is the etiological agent of PRRS, causing late-term abortions, stillbirths, and respiratory disease in pigs, incurring major economic losses to the worldwide pig industry. The virus is highly mutagenic and can be divided into two species, PRRSV-1 and PRRSV-2, each containing several subtypes. Current control strategies mainly involve biosecurity measures, depopulation, and vaccination. Vaccines are at best only partially protective against infection with heterologous subtypes and sublineages, and modified live vaccines have frequently been reported to revert to virulence. Here, we demonstrate that a genetic-control approach results in complete resistance to PRRSV infection in vivo CD163 is edited so as to remove the viral interaction domain while maintaining protein expression and biological function, averting any potential adverse effect associated with protein knockout. This research demonstrates a genetic-control approach with potential benefits in animal welfare as well as to the pork industry.

Precision engineering for PRRSV resistance in pigs: Macrophages from genome edited pigs lacking CD163 SRCR5 domain are fully resistant to both PRRSV genotypes while maintaining biological function.

Porcine Reproductive and Respiratory Syndrome (PRRS) is a panzootic infectious disease of pigs, causing major economic losses to the world-wide pig industry. PRRS manifests differently in pigs of all ages but primarily causes late-term abortions and stillbirths in sows and respiratory disease in piglets. The causative agent of the disease is the positive-strand RNA PRRS virus (PRRSV). PRRSV has a narrow host cell tropism, limited to cells of the monocyte/macrophage lineage. CD163 has been described as a fusion receptor for PRRSV, whereby the scavenger receptor cysteine-rich domain 5 (SRCR5) region was shown to be an interaction site for the virus in vitro. CD163 is expressed at high levels on the surface of macrophages, particularly in the respiratory system. Here we describe the application of CRISPR/Cas9 to pig zygotes, resulting in the generation of pigs with a deletion of Exon 7 of the CD163 gene, encoding SRCR5. Deletion of SRCR5 showed no adverse effects in pigs maintained under standard husbandry conditions with normal growth rates and complete blood counts observed. Pulmonary alveolar macrophages (PAMs) and peripheral blood monocytes (PBMCs) were isolated from the animals and assessed in vitro. Both PAMs and macrophages obtained from PBMCs by CSF1 stimulation (PMMs) show the characteristic differentiation and cell surface marker expression of macrophages of the respective origin. Expression and correct folding of the SRCR5 deletion CD163 on the surface of macrophages and biological activity of the protein as hemoglobin-haptoglobin scavenger was confirmed. Challenge of both PAMs and PMMs with PRRSV genotype 1, subtypes 1, 2, and 3 and PMMs with PRRSV genotype 2 showed complete resistance to viral infections assessed by replication. Confocal microscopy revealed the absence of replication structures in the SRCR5 CD163 deletion macrophages, indicating an inhibition of infection prior to gene expression, i.e. at entry/fusion or unpacking stages.

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.

Normalized long read RNA sequencing in chicken reveals transcriptome complexity similar to human.

BACKGROUND: Despite the significance of chicken as a model organism, our understanding of the chicken transcriptome is limited compared to human. This issue is common to all non-human vertebrate annotations due to the difficulty in transcript identification from short read RNAseq data. While previous studies have used single molecule long read sequencing for transcript discovery, they did not perform RNA normalization and 5'-cap selection which may have resulted in lower transcriptome coverage and truncated transcript sequences. RESULTS: We sequenced normalised chicken brain and embryo RNA libraries with Pacific Bioscience Iso-Seq. 5' cap selection was performed on the embryo library to provide methodological comparison. From these Iso-Seq sequencing projects, we have identified 60 k transcripts and 29 k genes within the chicken transcriptome. Of these, more than 20 k are novel lncRNA transcripts with ~3 k classified as sense exonic overlapping lncRNA, which is a class that is underrepresented in many vertebrate annotations. The relative proportion of alternative transcription events revealed striking similarities between the chicken and human transcriptomes while also providing explanations for previously observed genomic differences. CONCLUSIONS: Our results indicate that the chicken transcriptome is similar in complexity compared to human, and provide insights into other vertebrate biology. Our methodology demonstrates the potential of Iso-Seq sequencing to rapidly expand our knowledge of transcriptomics.

Quasispecies evolution of the prototypical genotype 1 porcine reproductive and respiratory syndrome virus early during in vivo infection is rapid and tissue specific.

Porcine reproductive and respiratory syndrome virus (PRRSV) is a major infectious threat to the pig industry worldwide. Increasing evidence suggests that microevolution within a quasispecies population can give rise to high sequence heterogeneity in PRRSV; potentially impacting the pathogenicity of the virus. Here, we report on micro-evolutionary events taking place within the viral quasispecies population in lung and lymph node 3 days post infection (dpi) following experimental in vivo infection with the prototypical Lelystad PRRSV (LV). Sequence analysis revealed 16 high frequency single nucleotide variants (SNV) or differences from the reference LV genome which are assumed to be representative of the consensus inoculum genome. Additionally, 49 other low frequency SNVs were also found in the inoculum population. At 3 dpi, a total of 9 and 10 SNVs of varying frequencies could already be detected in the LV population infecting the lung and lymph nodes, respectively. Interestingly, of these, three and four novel SNVs emerged independently in the two respective tissues when compared to the inoculum. The remaining variants, though already present at lower frequencies in the inoculum, were positively selected and their frequency increased within the quasispecies population. Hence, we were able to determine directly from tissues infected with PRRSV the repertoire of genetic variants within the viral quasispecies population. Our data also suggest that microevolution of these variants is rapid and some may be tissue-specific.

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.

Characterization of the Interactome of the Porcine Reproductive and Respiratory Syndrome Virus Nonstructural Protein 2 Reveals the Hyper Variable Region as a Binding Platform for Association with 14-3-3 Proteins.

Porcine reproductive and respiratory syndrome virus (PRRSV) is a major threat to the swine industry worldwide and hence global food security, exacerbated by a newly emerged highly pathogenic (HP-PRRSV) strain from China. PRRSV nonstructural protein 2 (nsp2) is a multifunctional polypeptide with strain-dependent influences on pathogenicity. A number of discrete functional regions have been identified on the protein. Quantitative label free proteomics was used to identify cellular binding partners of nsp2 expressed by HP-PRRSV. This allowed the identification of potential cellular interacting partners and the discrimination of nonspecific interactions. The interactome data were further investigated and validated using biological replicates and also compared with nsp2 from a low pathogenic (LP) strain of PRRSV. Validation included both forward and reverse pulldowns and confocal microscopy. The data indicated that nsp2 interacted with a number of cellular proteins including 14-3-3, CD2AP, and other components of cellular aggresomes. The hyper-variable region of nsp2 protein was identified as a binding platform for association with 14-3-3 proteins.

Distinct functional enrichment of transcriptional signatures in pigs with high and low IFN-gamma responses after vaccination with a porcine reproductive and respiratory syndrome virus (PRRSV).

Little is known about the host factor in the response to PRRSV vaccination. For this purpose, piglets were immunized with a commercial PRRSV-live vaccine and classified as high responders (HR) or low responders (LR) as regards to the frequencies of virus-specific IFN-γ-secreting cells. Six weeks post vaccination, PBMCs isolated from three individuals with the most extreme responses in each HR and LR groups and 3 unvaccinated controls, were either stimulated with phytohaemagglutinin, challenged with the vaccine or mock treated for 24 h, prior conducting transcriptional studies, gene ontology and pathway analyses. The LR group had very low neutralizing antibody levels and showed a higher number of down-regulated transcripts compared with the HR group (FDR < 0.2, P < 0.001). Down-regulated genes encoded chemoattractants, proinflammatory cytokines and the interferon-inducible GBP family, and showed enrichment in wounding (FDR < 3.6E-13), inflammation (FDR < 8E-12), defence (FDR < 8.7E-09) and immunity (FDR < 7.6E-08), suggesting immune response impairment. In the HR group, down-regulated genes were involved in protein transport (FDR < 4.77E-03), locomotory behavior (FDR < 5.47E-3), regulation of protein localization (FDR < 1.02E-02), and regulation of TNF superfamily member 15 and miR181. In contrast, the HR group presented up-regulated transcripts associated with wounding (FDR < 4.95). Moreover, IFN-γ was predicted to be an inhibited upstream regulator since IFN-γ pathways were associated with higher number of down-regulated genes in the LR (n = 40) than the HR (n = 10). Divergent responses to PRRSV-vaccination may be the result of the genetic background of the host.

Genome-wide association reveals QTL for growth, bone and in vivo carcass traits as assessed by computed tomography in Scottish Blackface lambs.

BACKGROUND: Improving meat quality including taste and tenderness is critical to the protection and development of markets for sheep meat. Phenotypic selection for such measures of meat quality is constrained by the fact that these parameters can only be measured post-slaughter. Carcass composition has an impact on meat quality and can be measured on live animals using advanced imaging technologies such as X-ray computed tomography (CT). Since carcass composition traits are heritable, they are potentially amenable to improvement through marker-assisted and genomic selection. We conducted a genome-wide association study (GWAS) on about 600 Scottish Blackface lambs for which detailed carcass composition phenotypes, including bone, fat and muscle components, had been captured using CT and which were genotyped for ~40,000 single nucleotide polymorphisms (SNPs) using the Illumina OvineSNP50 chip. RESULTS: We confirmed that the carcass composition traits were heritable with moderate to high (0.19-0.78) heritabilities. The GWAS analyses revealed multiple SNPs and quantitative trait loci (QTL) that were associated with effects on carcass composition traits and were significant at the genome-wide level. In particular, we identified a region on ovine chromosome 6 (OAR6) associated with bone weight and bone area that harboured SNPs with p values of 5.55 × 10(-8) and 2.63 × 10(-9), respectively. The same region had effects on fat area, fat density, fat weight and muscle density. We identified plausible positional candidate genes for these OAR6 QTL. We also detected a SNP that reached the genome-wide significance threshold with a p value of 7.28 × 10(-7) and was associated with muscle density on OAR1. Using a regional heritability mapping approach, we also detected regions on OAR3 and 24 that reached genome-wide significance for bone density. CONCLUSIONS: We identified QTL on OAR1, 3, 24 and particularly on OAR6 that are associated with effects on muscle, fat and bone traits. Based on available evidence that indicates that these traits are genetically correlated with meat quality traits, these associated SNPs have potential applications in selective breeding for improved meat quality. Further research is required to determine whether the effects associated with the OAR6 QTL are caused by a single gene or several closely-linked genes.

GO-FAANG meeting: a Gathering On Functional Annotation of Animal Genomes.

The Functional Annotation of Animal Genomes (FAANG) Consortium recently held a Gathering On FAANG (GO-FAANG) Workshop in Washington, DC on October 7-8, 2015. This consortium is a grass-roots organization formed to advance the annotation of newly assembled genomes of domesticated and non-model organisms (www.faang.org). The workshop gathered together from around the world a group of 100+ genome scientists, administrators, representatives of funding agencies and commodity groups to discuss the latest advancements of the consortium, new perspectives, next steps and implementation plans. The workshop was streamed live and recorded, and all talks, along with speaker slide presentations, are available at www.faang.org. In this report, we describe the major activities and outcomes of this meeting. We also provide updates on ongoing efforts to implement discussions and decisions taken at GO-FAANG to guide future FAANG activities. In summary, reference datasets are being established under pilot projects; plans for tissue sets, morphological classification and methods of sample collection for different tissues were organized; and core assays and data and meta-data analysis standards were established.

Signatures of diversifying selection in European pig breeds.

Following domestication, livestock breeds have experienced intense selection pressures for the development of desirable traits. This has resulted in a large diversity of breeds that display variation in many phenotypic traits, such as coat colour, muscle composition, early maturity, growth rate, body size, reproduction, and behaviour. To better understand the relationship between genomic composition and phenotypic diversity arising from breed development, the genomes of 13 traditional and commercial European pig breeds were scanned for signatures of diversifying selection using the Porcine60K SNP chip, applying a between-population (differentiation) approach. Signatures of diversifying selection between breeds were found in genomic regions associated with traits related to breed standard criteria, such as coat colour and ear morphology. Amino acid differences in the EDNRB gene appear to be associated with one of these signatures, and variation in the KITLG gene may be associated with another. Other selection signals were found in genomic regions including QTLs and genes associated with production traits such as reproduction, growth, and fat deposition. Some selection signatures were associated with regions showing evidence of introgression from Asian breeds. When the European breeds were compared with wild boar, genomic regions with high levels of differentiation harboured genes related to bone formation, growth, and fat deposition.

Identification and annotation of conserved promoters and macrophage-expressed genes in the pig genome.

BACKGROUND: The FANTOM5 consortium used Cap Analysis of Gene Expression (CAGE) tag sequencing to produce a comprehensive atlas of promoters and enhancers within the human and mouse genomes. We reasoned that the mapping of these regulatory elements to the pig genome could provide useful annotation and evidence to support assignment of orthology. RESULTS: For human transcription start sites (TSS) associated with annotated human-mouse orthologs, 17% mapped to the pig genome but not to the mouse, 10% mapped only to the mouse, and 55% mapped to both pig and mouse. Around 17% did not map to either species. The mapping percentages were lower where there was not clear orthology relationship, but in every case, mapping to pig was greater than to mouse, and the degree of homology was also greater. Combined mapping of mouse and human CAGE-defined promoters identified at least one putative conserved TSS for >16,000 protein-coding genes. About 54% of the predicted locations of regulatory elements in the pig genome were supported by CAGE and/or RNA-Seq analysis from pig macrophages. CONCLUSIONS: Comparative mapping of promoters and enhancers from humans and mice can provide useful preliminary annotation of other animal genomes. The data also confirm extensive gain and loss of regulatory elements between species, and the likelihood that pigs provide a better model than mice for human gene regulation and function.