Search

Abstract Background The application of DNA markers for the identification of biological samples from both human and non-human species is widespread and includes use in food authentication. In the food industry the financial incentive to substituting the true name of a food product with a higher value alternative is driving food fraud. This applies to British pork products where products derived from traditional pig breeds are of premium value. The objective of this study was to develop a genetic assay for regulatory authentication of traditional pig breed-labelled products in the porcine food industry in the United Kingdom. Results The dataset comprised of a comprehensive coverage of breed types present in Britain: 460 individuals from 7 traditional breeds, 5 commercial purebreds, 1 imported European breed and 1 imported Asian breed were genotyped using the PorcineSNP60 beadchip. Following breed-informative SNP selection, assignment power was calculated for increasing SNP panel size. A 96-plex assay created using the most informative SNPs revealed remarkably high genetic differentiation between the British pig breeds, with an average FST of 0.54 and Bayesian clustering analysis also indicated that they were distinct homogenous populations. The posterior probability of assignment of any individual of a presumed origin actually originating from that breed given an alternative breed origin was > 99.5% in 174 out of 182 contrasts, at a test value of log(LR) > 0. Validation of the 96-plex assay using independent test samples of known origin was successful; a subsequent survey of market samples revealed a high level of breed label conformity. Conclusion The newly created 96-plex assay using selected markers from the PorcineSNP60 beadchip enables powerful assignment of samples to traditional breed origin and can effectively identify mislabelling, providing a highly effective tool for DNA analysis in food forensics.

Abstract 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.

Abstract Background This work describes the first genome-wide analysis of the transcriptional landscape of the pig. A new porcine Affymetrix expression array was designed in order to provide comprehensive coverage of the known pig transcriptome. The new array was used to generate a genome-wide expression atlas of pig tissues derived from 62 tissue/cell types. These data were subjected to network correlation analysis and clustering. Results The analysis presented here provides a detailed functional clustering of the pig transcriptome where transcripts are grouped according to their expression pattern, so one can infer the function of an uncharacterized gene from the company it keeps and the locations in which it is expressed. We describe the overall transcriptional signatures present in the tissue atlas, where possible assigning those signatures to specific cell populations or pathways. In particular, we discuss the expression signatures associated with the gastrointestinal tract, an organ that was sampled at 15 sites along its length and whose biology in the pig is similar to human. We identify sets of genes that define specialized cellular compartments and region-specific digestive functions. Finally, we performed a network analysis of the transcription factors expressed in the gastrointestinal tract and demonstrate how they sub-divide into functional groups that may control cellular gastrointestinal development. Conclusions As an important livestock animal with a physiology that is more similar than mouse to man, we provide a major new resource for understanding gene expression with respect to the known physiology of mammalian tissues and cells. The data and analyses are available on the websites http://biogps.org and http://www.macrophages.com/pig-atlas.

Abstract Infection of sheep with the gastric nematode Teladorsagia circumcincta results in distinct Th2-type changes in the mucosa, including mucous neck cell and mast cell hyperplasia, eosinophilia, recruitment of IgA/IgE producing cells and neutrophils, altered T-cell subsets and mucosal hypertrophy. To address the protective mechanisms generated in animals on previous exposure to this parasite, gene expression profiling was carried out using samples of abomasal mucosa collected pre- and post- challenge from animals of differing immune status, using an experimental model of T. circumcincta infection. Recently developed ovine cDNA arrays were used to compare the abomasal responses of sheep immunised by trickle infection with worm-naïve sheep, following a single oral challenge of 50 000 T. circumcincta L3. Key changes were validated using qRT-PCR techniques. Immune animals demonstrated highly significant increases in levels of transcripts normally associated with cytotoxicity such as granulysin and granzymes A, B and H, as well as mucous-cell derived transcripts, predominantly calcium-activated chloride channel 1 (CLCA1). Challenge infection also induced up-regulation of transcripts potentially involved in initiating or modulating the immune response, such as heat shock proteins, complement factors and the chemokine CCL2. In contrast, there was marked infection-associated down-regulation of gene expression of members of the gastric lysozyme family. The changes in gene expression levels described here may reflect roles in direct anti-parasitic effects, immuno-modulation or tissue repair. (Funding; DEFRA/SHEFC (VT0102) and the BBSRC (BB/E01867X/1)).

Abstract Background Enterotoxigenic Escherichia coli (ETEC) that express the F4ab and F4ac fimbriae is a major contributor to diarrhoea outbreaks in the pig breeding industry, infecting both newborn and weaned piglets. Some pigs are resistant to this infection, and susceptibility is inherited as a simple dominant Mendelian trait. Indentifying the genetics behind this trait will greatly benefit pig welfare as well as the pig breeding industry by providing an opportunity to select against genetically susceptible animals, thereby reducing the number of diarrhoea outbreaks. The trait has recently been mapped by haplotype sharing to a 2.5 Mb region on pig chromosome 13, a region containing 18 annotated genes. Findings The coding regions of five candidate genes for susceptibility to ETEC F4ab/ac infection (TFRC, ACK1, MUC20, MUC4 and KIAA0226), all located in the 2.5 Mb region, were investigated for the presence of possible causative mutations. A total of 34 polymorphisms were identified in either coding regions or their flanking introns. The genotyping data for two of those were found to perfectly match the genotypes at the ETEC F4ab/ac locus, a G to C polymorphism in intron 11 of TFRC and a C to T silent polymorphism in exon 22 of KIAA0226. Transcriptional profiles of the five genes were investigated in a porcine tissue panel including various intestinal tissues. All five genes were expressed in intestinal tissues at different levels but none of the genes were found differentially expressed between ETEC F4ab/ac resistant and ETEC F4ab/ac susceptible animals in any of the tested tissues. Conclusions None of the identified polymorphisms are obvious causative mutations for ETEC F4ab/ac susceptibility, as they have no impact on the level of the overall mRNA expression nor predicted to influence the composition of the amino acids composition. However, we cannot exclude that the five tested genes are bona fide candidate genes for susceptibility to ETEC F4ab/ac infection since the identified polymorphism might affect the translational apparatus, alternative splice forms may exist and post translational mechanisms might contribute to disease susceptibility.

Abstract Background Genetic markers can be used to identify and verify the origin of individuals. Motivation for the inference of ancestry ranges from conservation genetics to forensic analysis. High density assays featuring Single Nucleotide Polymorphism (SNP) markers can be exploited to create a reduced panel containing the most informative markers for these purposes. The objectives of this study were to evaluate methods of marker selection and determine the minimum number of markers from the BovineSNP50 BeadChip required to verify the origin of individuals in European cattle breeds. Delta, Wright's FST, Weir & Cockerham's FST and PCA methods for population differentiation were compared. The level of informativeness of each SNP was estimated from the breed specific allele frequencies. Individual assignment analysis was performed using the ranked informative markers. Stringency levels were applied by log-likelihood ratio to assess the confidence of the assignment test. Results A 95% assignment success rate for the 384 individually genotyped animals was achieved with < 80, < 100, < 140 and < 200 SNP markers (with increasing stringency threshold levels) across all the examined methods for marker selection. No further gain in power of assignment was achieved by sampling in excess of 200 SNP markers. The marker selection method that required the lowest number of SNP markers to verify the animal's breed origin was Wright's FST (60 to 140 SNPs depending on the chosen degree of confidence). Certain breeds required fewer markers (< 100) to achieve 100% assignment success. In contrast, closely related breeds require more markers (~200) to achieve > 95% assignment success. The power of assignment success, and therefore the number of SNP markers required, is dependent on the levels of genetic heterogeneity and pool of samples considered. Conclusions While all SNP selection methods produced marker panels capable of breed identification, the power of assignment varied markedly among analysis methods. Thus, with effective exploration of available high density genetic markers, a diagnostic panel of highly informative markers can be produced.

Abstract 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.

Abstract Background Differences in gene expression may be caused by nearby DNA polymorphisms (cis regulation) or by interactions of gene control regions with polymorphic transcription factors (trans regulation). Trans acting loci are much harder to detect than cis acting loci and their effects are much more sensitive to genetic background. Results To quantify cis and trans regulation we correlated haplotype data with gene expression in two inbred mouse strains and two derived congenic lines. Upstream haplotype differences between the parental strains suggested that 30-43% of differentially expressed genes were differentially expressed because of cis haplotype differences. These cis regulated genes displayed consistent and relatively tissue-independent differential expression. We independently estimated from the congenic mice that 71-85% of genes were trans regulated. Cis regulated genes were associated with low p values (p < 0.005) for differential expression, whereas trans regulated genes were associated with values 0.005 < p < 0.05. The genes differentially expressed between congenics and controls were not a subset of those that were differentially expressed between the founder lines, showing that these were dependent on genetic background. For example, the cholesterol synthesis pathway was strongly differentially expressed in the congenic mice by indirect trans regulation but this was not observable in the parental mice. Conclusions The evidence that most gene regulation is trans and strongly influenced by genetic background, suggests that pathways that are modified by an allelic variant, may only exhibit differential expression in the specific genetic backgrounds in which they were identified. This has significant implications for the interpretation of any QTL mapping study.

Abstract Background Over the last decade, several studies have identified quantitative trait loci (QTL) affecting variation of immune related traits in mammals. Recent studies in humans and mice suggest that part of this variation may be caused by polymorphisms in genes involved in Toll-like receptor (TLR) signalling. In this project, we used a comparative approach to investigate the importance of TLR-related genes in comparison with other immunologically relevant genes for resistance traits in five species by associating their genomic location with previously published immune-related QTL regions. Results We report the genomic localisation of TLR1-10 and ten associated signalling molecules in sheep and pig using in-silico and/or radiation hybrid (RH) mapping techniques and compare their positions with their annotated homologues in the human, cattle and mouse whole genome sequences. We also report medium-density RH maps for porcine chromosomes 8 and 13. A comparative analysis of the positions of previously published relevant QTLs allowed the identification of homologous regions that are associated with similar health traits in several species and which contain TLR related and other immunologically relevant genes. Additional evidence was gathered by examining relevant gene expression and association studies. Conclusion This comparative genomic approach identified eight genes as potentially causative genes for variations of health related traits. These include susceptibility to clinical mastitis in dairy cattle, general disease resistance in sheep, cattle, humans and mice, and tolerance to protozoan infection in cattle and mice. Four TLR-related genes (TLR1, 6, MyD88, IRF3) appear to be the most likely candidate genes underlying QTL regions which control the resistance to the same or similar pathogens in several species. Further studies are required to investigate the potential role of polymorphisms within these genes.

No disponible

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.

Lawsonia intracellularis is the aetiological agent of proliferative enteropathy, an enteric disease endemic in swine. Survival in its intracellular niche of the ileum epithelial lining requires the capacity to subvert, repress or exploit the host immune response to create an environment conducive to bacterial propagation. To better understand how L. intracellularis survives in its intracellular niche, we have performed an investigation into the dynamic relationship between infection and the host autophagy response by immunohistochemistry in experimentally infected porcine ileum samples.Beclin1, a protein required early in the autophagy pathway was observed to be distributed with a basal to apical concentration gradient in the crypts of healthy piglets, whilst infected piglets were observed to have no gradient of distribution and an increase in the presence of Beclin1 in crypts with histological characteristics of L. intracellularis residence. Detecting microtubule-associated protein light chain 3 (LC3) is used as a method for monitoring autophagy progression as it associates with mature autophagosomes. For LC3 there was no notable change in signal intensity between crypts with characteristic L. intracellularis infection and healthy crypts of uninfected pigs. Finally, as p62 is degraded with the internal substrate of an autophagosome it was used to measure autophagic flux. There was no observed reduction or redistribution of p62.These preliminary results of the autophagy response in the ileum suggest that L. intracellularis affects autophagy. This disruption to host ileum homeostasis may provide a mechanism that assists in bacterial propagation and contributes to pathogenesis., (© 2021. The Author(s).)

There is very little information about how the genome is regulated in domestic pigs (Sus scrofa). This lack of knowledge hinders efforts to define and predict the effects of genetic variants in pig breeding programs. To address this knowledge gap, we need to identify regulatory sequences in the pig genome starting with regions of open chromatin. We used the "Improved Protocol for the Assay for Transposase-Accessible Chromatin (Omni-ATAC-Seq)" to identify putative regulatory regions in flash-frozen semitendinosus muscle from 24 male piglets. We collected samples from the smallest-, average-, and largest-sized male piglets from each litter through five developmental time points. Of the 4661 ATAC-Seq peaks identified that represent regions of open chromatin, >50% were within 1 kb of known transcription start sites. Differential read count analysis revealed 377 ATAC-Seq defined genomic regions where chromatin accessibility differed significantly across developmental time points. We found regions of open chromatin associated with downregulation of genes involved in muscle development that were present in small-sized fetal piglets but absent in large-sized fetal piglets at day 90 of gestation. The dataset that we have generated provides a resource for studies of genome regulation in pigs and contributes valuable functional annotation information to filter genetic variants for use in genomic selection in pig breeding programs., (© The Author(s) 2021. Published by Oxford University Press on behalf of Genetics Society of America.)

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., (© 2021. The Author(s).)

Background: The Pacific oyster (Crassostrea gigas) is a bivalve mollusc with vital roles in coastal ecosystems and aquaculture globally. While extensive genomic tools are available for C. gigas, highly contiguous reference genomes are required to support both fundamental and applied research. Herein we report the creation and annotation of a chromosome-level assembly for C. gigas., Findings: High-coverage long- and short-read sequence data generated on Pacific Biosciences and Illumina platforms were used to generate an initial assembly, which was then scaffolded into 10 pseudo-chromosomes using both Hi-C sequencing and a high-density linkage map. The assembly has a scaffold N50 of 58.4 Mb and a contig N50 of 1.8 Mb, representing a step advance on the previously published C. gigas assembly. Annotation based on Pacific Biosciences Iso-Seq and Illumina RNA-Seq resulted in identification of ∼30,000 putative protein-coding genes. Annotation of putative repeat elements highlighted an enrichment of Helitron rolling-circle transposable elements, suggesting their potential role in shaping the evolution of the C. gigas genome., Conclusions: This new chromosome-level assembly will be an enabling resource for genetics and genomics studies to support fundamental insight into bivalve biology, as well as for selective breeding of C. gigas in aquaculture., (© The Author(s) 2021. Published by Oxford University Press GigaScience.)

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.

The overall aim of the Ovine FAANG project is to provide a comprehensive annotation of the new highly contiguous sheep reference genome sequence ( Oar rambouillet v1.0 ). Mapping of transcription start sites (TSS) is a key first step in understanding transcript regulation and diversity. Using 56 tissue samples collected from the reference ewe Benz2616, we have performed a global analysis of TSS and TSS-Enhancer clusters using Cap Analysis Gene Expression (CAGE) sequencing. CAGE measures RNA expression by 5' cap-trapping and has been specifically designed to allow the characterization of TSS within promoters to single-nucleotide resolution. We have adapted an analysis pipeline that uses TagDust2 for clean-up and trimming, Bowtie2 for mapping, CAGEfightR for clustering, and the Integrative Genomics Viewer (IGV) for visualization. Mapping of CAGE tags indicated that the expression levels of CAGE tag clusters varied across tissues. Expression profiles across tissues were validated using corresponding polyA+ mRNA-Seq data from the same samples. After removal of CAGE tags with <10 read counts, 39.3% of TSS overlapped with 5' ends of 31,113 transcripts that had been previously annotated by NCBI (out of a total of 56,308 from the NCBI annotation). For 25,195 of the transcripts, previously annotated by NCBI, no TSS meeting stringent criteria were identified. A further 14.7% of TSS mapped to within 50 bp of annotated promoter regions. Intersecting these predicted TSS regions with annotated promoter regions (±50 bp) revealed 46% of the predicted TSS were "novel" and previously un-annotated. Using whole-genome bisulfite sequencing data from the same tissues, we were able to determine that a proportion of these "novel" TSS were hypo-methylated (32.2%) indicating that they are likely to be reproducible rather than "noise". This global analysis of TSS in sheep will significantly enhance the annotation of gene models in the new ovine reference assembly. Our analyses provide one of the highest resolution annotations of transcript regulation and diversity in a livestock species to date., (Copyright © 2020 Salavati, Caulton, Clark, Gazova, Smith, Worley, Cockett, Archibald, Clarke, Murdoch and Clark.)

More people globally depend on the water buffalo than any other domesticated species, and as the most closely related domesticated species to cattle they can provide important insights into the shared evolutionary basis of domestication. Here, we sequence the genomes of 79 water buffalo across seven breeds and compare patterns of between breed selective sweeps with those seen for 294 cattle genomes representing 13 global breeds. The genomic regions under selection between cattle breeds significantly overlap regions linked to stature in human genetic studies, with a disproportionate number of these loci also shown to be under selection between water buffalo breeds. Investigation of potential functional variants in the water buffalo genome identifies a rare example of convergent domestication down to the same mutation having independently occurred and been selected for across domesticated species. Cross-species comparisons of recent selective sweeps can consequently help identify and refine important loci linked to domestication.

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., (© The Author(s) 2020. Published by Oxford University Press.)

Lawsonia intracellularis is a Gram-negative obligate intracellular bacterium that is the aetiological agent of proliferative enteropathy (PE), a common intestinal disease of major economic importance in pigs and other animal species. To date, progress in understanding the biology of L. intracellularis for improved disease control has been hampered by the inability to culture the organism in vitro . In particular, our understanding of the genomic diversity and population structure of clinical L. intercellularis is very limited. Here, we utilized a metagenomic shotgun approach to directly sequence and assemble 21 L . intracellularis genomes from faecal and ileum samples of infected pigs and horses across three continents. Phylogenetic analysis revealed a genetically monomorphic clonal lineage responsible for infections in pigs, with distinct subtypes associated with infections in horses. The genome was highly conserved, with 94 % of genes shared by all isolates and a very small accessory genome made up of only 84 genes across all sequenced strains. In part, the accessory genome was represented by regions with a high density of SNPs, indicative of recombination events importing novel gene alleles. In summary, our analysis provides the first view of the population structure for L. intracellularis , revealing a single major lineage associated with disease of pigs. The limited diversity and broad geographical distribution suggest the recent emergence and clonal expansion of an important livestock pathogen.

The domestic pig ( Sus scrofa ) is both an economically important livestock species and a model for biomedical research. Two highly contiguous pig reference genomes have recently been released. To support functional annotation of the pig genomes and comparative analysis with large human transcriptomic data sets, we aimed to create a pig gene expression atlas. To achieve this objective, we extended a previous approach developed for the chicken. We downloaded RNAseq data sets from public repositories, down-sampled to a common depth, and quantified expression against a reference transcriptome using the mRNA quantitation tool, Kallisto. We then used the network analysis tool Graphia to identify clusters of transcripts that were coexpressed across the merged data set. Consistent with the principle of guilt-by-association, we identified coexpression clusters that were highly tissue or cell-type restricted and contained transcription factors that have previously been implicated in lineage determination. Other clusters were enriched for transcripts associated with biological processes, such as the cell cycle and oxidative phosphorylation. The same approach was used to identify coexpression clusters within RNAseq data from multiple individual liver and brain samples, highlighting cell type, process, and region-specific gene expression. Evidence of conserved expression can add confidence to assignment of orthology between pig and human genes. Many transcripts currently identified as novel genes with ENSSSCG or LOC IDs were found to be coexpressed with annotated neighbouring transcripts in the same orientation, indicating they may be products of the same transcriptional unit. The meta-analytic approach to utilising public RNAseq data is extendable to include new data sets and new species and provides a framework to support the Functional Annotation of Animals Genomes (FAANG) initiative., (Copyright © 2020 Summers, Bush, Wu, Su, Muriuki, Clark, Finlayson, Eory, Waddell, Talbot, Archibald and Hume.)

The domestic water buffalo ( Bubalus bubalis ) makes a major contribution to the global agricultural economy in the form of milk, meat, hides, and draught power. The global water buffalo population is predominantly found in Asia, and per head of population more people depend upon the buffalo than on any other livestock species. Despite its agricultural importance, there are comparatively fewer genomic and transcriptomic resources available for buffalo than for other livestock species. We have generated a large-scale gene expression atlas covering multiple tissue and cell types from all major organ systems collected from three breeds of riverine water buffalo (Mediterranean, Pandharpuri and Bhadawari) and used the network analysis tool Graphia Professional to identify clusters of genes with similar expression profiles. Alongside similar data, we and others have generated for ruminants as part of the Functional Annotation of Animal Genomes Consortium; this comprehensive transcriptome supports functional annotation and comparative analysis of the water buffalo genome.

One of the most significant physiological challenges to neonatal and juvenile ruminants is the development and establishment of the rumen. Using a subset of RNA-Seq data from our high-resolution atlas of gene expression in sheep ( Ovis aries ) we have provided the first comprehensive characterization of transcription of the entire gastrointestinal (GI) tract during the transition from pre-ruminant to ruminant. The dataset comprises 164 tissue samples from sheep at four different time points (birth, one week, 8 weeks and adult). Using network cluster analysis we illustrate how the complexity of the GI tract is reflected in tissue- and developmental stage-specific differences in gene expression. The most significant transcriptional differences between neonatal and adult sheep were observed in the rumen complex. Comparative analysis of gene expression in three GI tract tissues from age-matched sheep and goats revealed species-specific differences in genes involved in immunity and metabolism. This study improves our understanding of the transcriptomic mechanisms involved in the transition from pre-ruminant to ruminant by identifying key genes involved in immunity, microbe recognition and metabolism. The results form a basis for future studies linking gene expression with microbial colonization of the developing GI tract and provide a foundation to improve ruminant efficiency and productivity through identifying potential targets for novel therapeutics and gene editing., (Copyright © 2019 Bush et al.)

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., Competing Interests: No authors have competing interests. GAW and AEH are both employed by JSR Genetics Limited who historically marketed the genetic line used in this study as a sireline boar product to commercial farmers. The line has since has been culled and JSR Genetics Limited no longer receive any revenue from the product.

The Functional Annotation of ANimal Genomes (FAANG) project aims, through a coordinated international effort, to provide high quality functional annotation of animal genomes with an initial focus on farmed and companion animals. A key goal of the initiative is to ensure high quality and rich supporting metadata to describe the project's animals, specimens, cell cultures and experimental assays. By defining rich sample and experimental metadata standards and promoting best practices in data descriptions, deposition and openness, FAANG champions higher quality and reusability of published datasets. FAANG has established a Data Coordination Centre, which sits at the heart of the Metadata and Data Sharing Committee. It continues to evolve the metadata standards, support submissions and, crucially, create powerful and accessible tools to support deposition and validation of metadata. FAANG conforms to the findable, accessible, interoperable, and reusable (FAIR) data principles, with high quality, open access and functionally interlinked data. In addition to data generated by FAANG members and specific FAANG projects, existing datasets that meet the main-or more permissive legacy-standards are incorporated into a central, focused, functional data resource portal for the entire farmed and companion animal community. Through clear and effective metadata standards, validation and conversion software, combined with promotion of best practices in metadata implementation, FAANG aims to maximise effectiveness and inter-comparability of assay data. This supports the community to create a rich genome-to-phenotype resource and promotes continuing improvements in animal data standards as a whole., (© 2018 Stichting International Foundation for Animal Genetics.)

The human population is growing, and as a result we need to produce more food whilst reducing the impact of farming on the environment. Selective breeding and genomic selection have had a transformational impact on livestock productivity, and now transgenic and genome-editing technologies offer exciting opportunities for the production of fitter, healthier and more-productive livestock. Here, we review recent progress in the application of genome editing to farmed animal species and discuss the potential impact on our ability to produce food.

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., (Copyright © 2018 Burkard et al.)

Porcine proliferative ileitis is a major economic burden for the swine industry, affecting growing pigs and young adult pigs. In this study, the protective efficacy of an inactivated, injectable whole-cell bacteria vaccine against L. intracellularis - Porcilis® Ileitis was evaluated under field conditions. Eighty-five, three-week-old pigs on a commercial farrow-to-finish farm were vaccinated by the intramuscular route, either with a dose of injectable vaccine, or with saline. A subset of vaccinates and control pigs were necropsied at 21 days post-challenge. Incidence and severity of ileitis were evaluated by gross and microscopic observation of ileal tissues. Colonization of the gut after challenge was examined by L. intracellularis-specific immunohistochemistry, and qPCR of ileal scrapings. Integrity of the intestinal barrier was evaluated to quantify a range of intestinal markers including secreted mucin and intestinal alkaline phosphatase, and innate immune markers including Caspase-3 and Calprotectin. A second subset of pigs was monitored for fecal shedding of L. intracellularis, until resolution of shedding. Our investigation indicated that Porcilis Ileitis provided robust protection against ileitis, reduced bacterial shedding 15-fold (p < .05) and preserved normal gut barrier function in the face of an experimental challenge with virulent L. intracellularis., (Copyright © 2018 Merck Sharp & Dohme Corp, Authors. Published by Elsevier Ltd.. All rights reserved.)

Porcine reproductive and respiratory syndrome virus (PRRSv) infection alters the host's cellular and humoral immune response. Immunity against PRRSv is multigenic and vary between individuals. The aim of the present study was to compare several genes that encode for molecules involved in the immune response between two groups of vaccinated pigs that experienced short or long viremic periods after PRRSv challenge. These analyses include the sequencing of four SLA Class I, two Class II allele groups, and CD163, plus the analysis by quantitative realtime qRT-PCR of the constitutive expression of TLR2, TLR3, TLR4, TLR7, TLR8 and TLR9 mRNA and other molecules in peripheral blood mononuclear cells.

Balanced chromosomal aberrations have been shown to affect fertility in most species studied, often leading to hypoprolificacy (reduced litter size) in domestic animals such as pigs. With an increasing emphasis in modern food production on the use of a small population of high quality males for artificial insemination, the potential economic and environmental costs of hypoprolific boars, bulls, rams etc. are considerable. There is therefore a need for novel tools to facilitate rapid, cost-effective chromosome translocation screening. This has previously been achieved by standard karyotype analysis; however, this approach relies on a significant level of expertise and is limited in its ability to identify subtle, cryptic translocations. To address this problem, we developed a novel device and protocol for translocation screening using subtelomeric probes and fluorescence in situ hybridisation. Probes were designed using BACs (bacterial artificial chromosomes) from the subtelomeric region of the short (p-arm) and long (q-arm) of each porcine chromosome. They were directly labelled with FITC or Texas Red (p-arm and q-arm respectively) prior to application of a 'Multiprobe' device, thereby enabling simultaneous detection of each individual porcine chromosome on a single slide. Initial experiments designed to isolate BACs in subtelomeric regions led to the discovery of a series of incorrectly mapped regions in the porcine genome assembly (from a total of 82 BACs, only 45 BACs mapped correctly). Our work therefore highlights the importance of accurate physical mapping of newly sequenced genomes. The system herein described allows for robust and comprehensive analysis of the porcine karyotype, an adjunct to classical cytogenetics that provides a valuable tool to expedite efficient, cost effective food production., (© 2017 The Authors. Animal Genetics published by John Wiley & Sons Ltd on behalf of Stichting International Foundation for Animal Genetics.)

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.

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.

Lawsonia intracellularis is an obligate intracellular bacterial pathogen that causes proliferative enteropathy (PE) in pigs. L. intracellularis infection causes extensive intestinal crypt cell proliferation and inhibits secretory and absorptive cell differentiation. However, the affected host upstream cellular pathways leading to PE are still unknown. β-catenin/Wnt signalling is essential in maintaining intestinal stem cell (ISC) proliferation and self-renewal capacity, while Notch signalling governs differentiation of secretory and absorptive lineage specification. Therefore, in this report we used immunofluorescence (IF) and quantitative reverse transcriptase PCR (RTqPCR) to examine β-catenin/Wnt and Notch-1 signalling levels in uninfected and L. intracellularis infected pig ileums at 3, 7, 14, 21 and 28 days post challenge (dpc). We found that while the significant increase in Ki67+ nuclei in crypts at the peak of L. intracellularis infection suggested enhanced cell proliferation, the expression of c-MYC and ASCL2, promoters of cell growth and ISC proliferation respectively, was down-regulated. Peak infection also coincided with enhanced cytosolic and membrane-associated β-catenin staining and induction of AXIN2 and SOX9 transcripts, both encoding negative regulators of β-catenin/Wnt signalling and suggesting a potential alteration to β-catenin/Wnt signalling levels, with differential regulation of the expression of its target genes. We found that induction of HES1 and OLFM4 and the down-regulation of ATOH1 transcript levels was consistent with the increased Notch-1 signalling in crypts at the peak of infection. Interestingly, the significant down-regulation of ATOH1 transcript levels coincided with the depletion of MUC2 expression at 14 dpc, consistent with the role of ATOH1 in promoting goblet cell maturation. The lack of significant change to LGR5 transcript levels at the peak of infection suggested that the crypt hyperplasia was not due to the expansion of ISC population. Overall, simultaneous induction of Notch-1 signalling and the attenuation of β-catenin/Wnt pathway appear to be associated with the inhibition of goblet cell maturation and enhanced crypt cell proliferation at the peak of L. intracellularis infection. Moreover, the apparent differential regulation of apoptosis between crypt and lumen cells together with the strong induction of Notch-1 signalling and the enhanced SOX9 expression along crypts 14 dpc suggest an expansion of actively dividing transit amplifying and/or absorptive progenitor cells and provide a potential basis for understanding the development and maintenance of PE.

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.

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.

Background. Ruminants are successful herbivorous mammals, in part due to their specialized forestomachs, the rumen complex, which facilitates the conversion of feed to soluble nutrients by micro-organisms. Is the rumen complex a modified stomach expressing new epithelial (cornification) and metabolic programs, or a specialised stratified epithelium that has acquired new metabolic activities, potentially similar to those of the colon? How has the presence of the rumen affected other sections of the gastrointestinal tract (GIT) of ruminants compared to non-ruminants? Methods. Transcriptome data from 11 tissues covering the sheep GIT, two stratified epithelial and two control tissues, was analysed using principal components to cluster tissues based on gene expression profile similarity. Expression profiles of genes along the sheep GIT were used to generate a network to identify genes enriched for expression in different compartments of the GIT. The data from sheep was compared to similar data sets from two non-ruminants, pigs (closely related) and humans (more distantly related). Results. The rumen transcriptome clustered with the skin and tonsil, but not the GIT transcriptomes, driven by genes from the epidermal differentiation complex, and genes encoding stratified epithelium keratins and innate immunity proteins. By analysing all of the gene expression profiles across tissues together 16 major clusters were identified. The strongest of these, and consistent with the high turnover rate of the GIT, showed a marked enrichment of cell cycle process genes (P = 1.4 E-46), across the whole GIT, relative to liver and muscle, with highest expression in the caecum followed by colon and rumen. The expression patterns of several membrane transporters (chloride, zinc, nucleosides, amino acids, fatty acids, cholesterol and bile acids) along the GIT was very similar in sheep, pig and humans. In contrast, short chain fatty acid uptake and metabolism appeared to be different between the species and different between the rumen and colon in sheep. The importance of nitrogen and iodine recycling in sheep was highlighted by the highly preferential expression of SLC14A1-urea (rumen), RHBG-ammonia (intestines) and SLC5A5-iodine (abomasum). The gene encoding a poorly characterized member of the maltase-glucoamylase family (MGAM2), predicted to play a role in the degradation of starch or glycogen, was highly expressed in the small and large intestines. Discussion. The rumen appears to be a specialised stratified cornified epithelium, probably derived from the oesophagus, which has gained some liver-like and other specialized metabolic functions, but probably not by expression of pre-existing colon metabolic programs. Changes in gene transcription downstream of the rumen also appear have occurred as a consequence of the evolution of the rumen and its effect on nutrient composition flowing down the GIT.

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.

Genetic selection against boar taint, which is caused by high skatole and androstenone concentrations in fat, is a more acceptable alternative than is the current practice of castration. Genomic predictors offer an opportunity to overcome the limitations of such selection caused by the phenotype being expressed only in males at slaughter, and this study evaluated different approaches to obtain such predictors. Samples from 1000 pigs were included in a design which was dominated by 421 sib pairs, each pair having one animal with high and one with low skatole concentration (≥0.3 μg/g). All samples were measured for both skatole and androstenone and genotyped using the Illumina SNP60 porcine BeadChip for 62 153 single nucleotide polymorphisms. The accuracy of predicting phenotypes was assessed by cross-validation using six different genomic evaluation methods: genomic best linear unbiased prediction (GBLUP) and five Bayesian regression methods. In addition, this was compared to the accuracy of predictions using only QTL that showed genome-wide significance. The range of accuracies obtained by different prediction methods was narrow for androstenone, between 0.29 (Bayes Lasso) and 0.31 (Bayes B), and wider for skatole, between 0.21 (GBLUP) and 0.26 (Bayes SSVS). Relative accuracies, corrected for h(2) , were 0.54-0.56 and 0.75-0.94 for androstenone and skatole respectively. The whole-genome evaluation methods gave greater accuracy than using only the QTL detected in the data. The results demonstrate that GBLUP for androstenone is the simplest genomic technology to implement and was also close to the most accurate method. More specialised models may be preferable for skatole., (© 2015 The Authors. Animal Genetics published by John Wiley & Sons Ltd on behalf of Stichting International Foundation for Animal Genetics.)

Many applications of high throughput sequencing rely on the availability of an accurate reference genome. Variant calling often produces large data sets that cannot be realistically validated and which may contain large numbers of false-positives. Errors in the reference assembly increase the number of false-positives. While resources are available to aid in the filtering of variants from human data, for other species these do not yet exist and strict filtering techniques must be employed which are more likely to exclude true-positives. This work assesses the accuracy of the pig reference genome (Sscrofa10.2) using whole genome sequencing reads from the Duroc sow whose genome the assembly was based on. Indicators of structural variation including high regional coverage, unexpected insert sizes, improper pairing and homozygous variants were used to identify low quality (LQ) regions of the assembly. Low coverage (LC) regions were also identified and analyzed separately. The LQ regions covered 13.85% of the genome, the LC regions covered 26.6% of the genome and combined (LQLC) they covered 33.07% of the genome. Over half of dbSNP variants were located in the LQLC regions. Of copy number variable regions identified in a previous study, 86.3% were located in the LQLC regions. The regions were also enriched for gene predictions from RNA-seq data with 42.98% falling in the LQLC regions. Excluding variants in the LQ, LC, or LQLC from future analyses will help reduce the number of false-positive variant calls. Researchers using WGS data should be aware that the current pig reference genome does not give an accurate representation of the copy number of alleles in the original Duroc sow's 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.

The expression patterns of secreted (MUC2 and MUC5AC) and membrane-tethered (MUC1, MUC4, MUC12 and MUC13) mucins were monitored in healthy pigs and pigs challenged orally with Lawsonia intracellularis. These results showed that the regulation of mucin gene expression is distinctive along the GI tract of the healthy pig, and may reflect an association between the function of the mucin subtypes and different physiological demands at various sites. We identified a specific depletion of secreted MUC2 from goblet cells in infected pigs that correlated with the increased level of intracellular bacteria in crypt cells. We concluded that L. intracellularis may influence MUC2 production, thereby altering the mucus barrier and enabling cellular invasion., (Copyright © 2015 The Authors. Published by Elsevier B.V. All rights reserved.)

Whipworms (Trichuris spp.) infect a variety of hosts, including domestic animals and humans. Of considerable interest is the porcine whipworm, T. suis, which is particularly prevalent in outdoor production systems. High infection levels may cause growth retardation, anaemia and haemorrhagic diarrhoea. A significant proportion of the variation in Trichuris faecal egg count (FEC) has been attributed to the host's genetic make-up. The aim of the present study was to identify genetic loci associated with resistance to T. suis in pigs. We used single nucleotide polymorphism (SNP) markers to perform a whole-genome scan of an F1 resource population (n = 195) trickle-infected with T. suis. A measured genotype analysis revealed a putative quantitative trait locus (QTL) for T. suis FEC on chromosome 13 covering ∼ 4.5 Mbp, although none of the SNPs reached genome-wide significance. We tested the hypothesis that this region of SSC13 harboured genes with effects on T. suis burden by genotyping three SNPs within the putative QTL in unrelated pigs exposed to either experimental or natural T. suis infections and from which we had FEC (n = 113) or worm counts (n = 178). In these studies, two of the SNPs (rs55618716, ST) were associated with FEC (P < 0.01), thus confirming our initial findings. However, we did not find any of the SNPs to be associated with T. suis worm burden. In conclusion, our study demonstrates that genetic markers for resistance to T. suis as indicated by low FEC can be identified in pigs., (Copyright © 2015 Elsevier B.V. All rights reserved.)