Search

The Farm Animal Genotype-Tissue Expression (FarmGTEx) project has been established to develop a public resource of genetic regulatory variants in livestock, which is essential for linking genetic polymorphisms to variation in phenotypes, helping fundamental biological discovery and exploitation in animal breeding and human biomedicine. Here we show results from the pilot phase of PigGTEx by processing 5,457 RNA-sequencing and 1,602 whole-genome sequencing samples passing quality control from pigs. We build a pig genotype imputation panel and associate millions of genetic variants with five types of transcriptomic phenotypes in 34 tissues. We evaluate tissue specificity of regulatory effects and elucidate molecular mechanisms of their action using multi-omics data. Leveraging this resource, we decipher regulatory mechanisms underlying 207 pig complex phenotypes and demonstrate the similarity of pigs to humans in gene expression and the genetic regulation behind complex phenotypes, supporting the importance of pigs as a human biomedical model.

The functional annotation of livestock genomes is crucial for understanding the molecular mechanisms that underpin complex traits of economic importance, adaptive evolution and comparative genomics. Here, we provide the most comprehensive catalogue to date of regulatory elements in the pig (Sus scrofa) by integrating 223 epigenomic and transcriptomic data sets, representing 14 biologically important tissues. We systematically describe the dynamic epigenetic landscape across tissues by functionally annotating 15 different chromatin states and defining their tissue-specific regulatory activities. We demonstrate that genomic variants associated with complex traits and adaptive evolution in pig are significantly enriched in active promoters and enhancers. Furthermore, we reveal distinct tissue-specific regulatory selection between Asian and European pig domestication processes. Compared with human and mouse epigenomes, we show that porcine regulatory elements are more conserved in DNA sequence, under both rapid and slow evolution, than those under neutral evolution across pig, mouse, and human. Finally, we provide biological insights on tissue-specific regulatory conservation, and by integrating 47 human genome-wide association studies, we demonstrate that, depending on the traits, mouse or pig might be more appropriate biomedical models for different complex traits and diseases.

Feed efficiency (FE) in pigs is an important factor in the profitability of pig farming operations. It refers to the ability of a pig to convert the feed it consumes into body weight. We used two metrics to measure FE: feed conversion ratio and average daily residual feed intake. A previous genome-wide association study and transcriptome study in crossbred pigs identified two QTL regions on SSC9 associated with residual feed intake and pointed out two candidate genes of interest: (a) the gene encoding the Aryl Hydrocarbon Receptor gene (AHR) transcription factor; and (b) the Dynein, Axonemal, Heavy Polypeptide 11 gene (DNAH11). The previous study identified missense mutations in both genes leading to a conservative substitution of glycine to cysteine in AHR (AHR_rs339939442) and two non-conservative substitutions in DNAH11, where arginine is replaced by threonine (DNAH11_rs325475644) and alanine is replaced by threonine (DNAH11_rs346074031). We have now genotyped the missense mutations in independent cohorts of 107 Duroc, 155 Landrace and 160 Yorkshire pigs to substantiate further if these variants directly impact FE-related phenotypes. We verified that allele T of AHR_rs339939442 in AHR improves FE in Yorkshire pigs. Genotype GG of AHR_rs339939442 was fixed in Duroc pigs. We also confirmed that the variants rs325475644 and rs346074031 in DNAH11 did not affect FE. The findings contribute valuable insights into the genetic mechanisms governing FE in pigs, potentially offering contributions for future enhancements of FE., (© 2024 Stichting International Foundation for Animal Genetics.)

Lean meat percentage is a critical production trait in pig breeding systems with direct implications for the sustainability of the industry. In this study, we conducted a genome-wide association study for lean meat percentage using a cohort of 850 Duroc × (Landrace × Yorkshire) crossbred pigs and we identified QTL on SSC3 and SSC18. Based on the predicted effect of imputed variants and using the PigGTEx database of molecular QTL, we prioritized candidate genes and SNPs located within the QTL regions, which may be involved in the regulation of porcine leanness. Our results indicate that a nonsense mutation in ZC3HAV1L on SSC18 has a direct effect on lean meat percentage., (© 2024 The Authors. Animal Genetics published by John Wiley & Sons Ltd on behalf of Stichting International Foundation for Animal Genetics.)

Post-weaning diarrhea in pigs is a considerable challenge in the pig farming industry due to its effect on animal welfare and production costs, as well as the large volume of antibiotics, which are used to treat diarrhea in pigs after weaning. Previous studies have revealed loci on SSC6 and SSC13 associated with susceptibility to specific diarrhea causing pathogens. This study aimed to identify new genetic loci for resistance to diarrhea based on phenotypic data. In depth clinical characterization of diarrhea was performed in 257 pigs belonging to two herds during the first 14 days post weaning. The daily diarrhea assessments were used for the classification of pigs into case and control groups. Pigs were assigned to case and control groups based only on the incidence of diarrhea in the second week of the study in order to differentiate between differences in etiology. Genome-wide association studies and metabolomics association analysis were performed in order to identify new biological determinants for diarrhea susceptibility. With the present work, we revealed a new locus for diarrhea resistance on SSC16. Furthermore, studies of metabolomics in the same pigs revealed one metabolite associated with diarrhea., (© 2024 The Authors. Animal Genetics published by John Wiley & Sons Ltd on behalf of Stichting International Foundation for Animal Genetics.)

Cystinuria is a genetic disease that can lead to cystine urolith formation. The English bulldog is the dog breed most frequently affected. In this breed, three missense mutations have been suggested to be associated with cystinuria: c.568A>G and c.2086A>G in SLC3A1 and c.649G>A in SLC7A9. In this study, the occurrence of these three mutations in the Danish population of English bulldogs was investigated. Seventy-one English bulldogs were genotyped using TaqMan assays. The dogs' owners were given questionnaires concerning the medical histories of their dogs. Allele frequencies of 0.40, 0.40, and 0.52 were found for the mutant alleles in the three loci: c.568A>G, c.2086A>G, and c.649G>A, respectively. For both mutations in SLC3A1, a statistically significant association was found between cystinuria and homozygosity for the G allele among male, English bulldogs. For the mutation in SLC7A9, there was no statistically significant association between homozygosity for the mutant allele and cystinuria. Due to high allele frequencies, limited genetic diversity, continued uncertainty about the genetic background of cystinuria, and more severe health problems in the breed, selection based on genetic testing for the mutations in SLC3A1 cannot be recommended in the Danish population of English bulldogs. However, results of the genetic test may be used as a guide to recommend prophylactic treatment., (© 2023 The Authors. Animal Genetics published by John Wiley & Sons Ltd on behalf of Stichting International Foundation for Animal Genetics.)

Whole-genome sequencing of 217 animals from three Danish commercial pig breeds (Duroc, Landrace [LL], and Yorkshire [YY]) was performed. Twenty-six million single-nucleotide polymorphisms (SNPs) and 8 million insertions or deletions (indels) were uncovered. Among the SNPs, 493,099 variants were located in coding sequences, and 29,430 were predicted to have a high functional impact such as gain or loss of stop codon. Using the whole-genome sequence dataset as the reference, the imputation accuracy for pigs genotyped with high-density SNP chips was examined. The overall average imputation accuracy for all biallelic variants (SNP and indel) was 0.69, while it was 0.83 for variants with minor allele frequency > 0.1. This study provides whole-genome reference data to impute SNP chip-genotyped animals for further studies to fine map quantitative trait loci as well as improving the prediction accuracy in genomic selection. Signatures of selection were identified both through analyses of fixation and differentiation to reveal selective sweeps that may have had prominent roles during breed development or subsequent divergent selection. However, the fixation indices did not indicate a strong divergence among these three breeds. In LL and YY, the integrated haplotype score identified genomic regions under recent selection. These regions contained genes for olfactory receptors and oxidoreductases. Olfactory receptor genes that might have played a major role in the domestication were previously reported to have been under selection in several species including cattle and swine., (© The Author(s) 2020. Published by Oxford University Press on behalf of the American Society of Animal Science. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

Reprogramming of adipocyte function in obesity is implicated in metabolic disorders like type 2 diabetes. Here, we used the pig, an animal model sharing many physiological and pathophysiological similarities with humans, to perform in-depth epigenomic and transcriptomic characterization of pure adipocyte fractions. Using a combined DNA methylation capture sequencing and Reduced Representation bisulfite sequencing (RRBS) strategy in 11 lean and 12 obese pigs, we identified in 3529 differentially methylated regions (DMRs) located at close proximity to-, or within genes in the adipocytes. By sequencing of the transcriptome from the same fraction of isolated adipocytes, we identified 276 differentially expressed transcripts with at least one or more DMR. These transcripts were over-represented in gene pathways related to MAPK, metabolic and insulin signaling. Using a candidate gene approach, we further characterized 13 genes potentially regulated by DNA methylation and identified putative transcription factor binding sites that could be affected by the differential methylation in obesity. Our data constitute a valuable resource for further investigations aiming to delineate the epigenetic etiology of metabolic disorders., (Copyright © 2019 Jacobsen, Havgaard, Anthon, Mentzel, Cirera, Krogh, Pundhir, Karlskov-Mortensen, Bruun, Lesnik, Guerin, Gorodkin, Jørgensen, Fredholm and Barrès.)

We have established a pig resource population specifically designed to elucidate the genetics involved in development of obesity and obesity related co-morbidities by crossing the obesity prone Göttingen Minipig breed with two lean production pig breeds. In this study we have performed genome wide association (GWA) to identify loci with effect on blood lipid levels. The most significantly associated single nucleotide polymorphisms (SNPs) were used for linkage disequilibrium (LD) and haplotype analyses. Three separate haploblocks which influence the ratio between high density lipoprotein cholesterol and total cholesterol (HDL-C/CT), triglycerides (TG) and low density lipoprotein cholesterol (LDL-C) levels respectively were identified on Sus Scrofa chromosome 3 (SSC3). Large additive genetic effects were found for the HDL-C/CT and LDL-C haplotypes. Haplotypes segregating from Göttingen Minipigs were shown to impose a positive effect on blood lipid levels. Thus, the genetic profile of the Göttingen Minipig breed seems to support a phenotype comparable to the metabolic healthy obese (MHO) phenotype in humans.

The pig is a well-known animal model used to investigate genetic and mechanistic aspects of human disease biology. They are particularly useful in the context of obesity and metabolic diseases because other widely used models (e.g. mice) do not completely recapitulate key pathophysiological features associated with these diseases in humans. Therefore, we established a F2 pig resource population (n = 564) designed to elucidate the genetics underlying obesity and metabolic phenotypes. Segregation of obesity traits was ensured by using breeds highly divergent with respect to obesity traits in the parental generation. Several obesity and metabolic phenotypes were recorded (n = 35) from birth to slaughter (242 ± 48 days), including body composition determined at about two months of age (63 ± 10 days) via dual-energy x-ray absorptiometry (DXA) scanning. All pigs were genotyped using Illumina Porcine 60k SNP Beadchip and a combined linkage disequilibrium-linkage analysis was used to identify genome-wide significant associations for collected phenotypes. We identified 229 QTLs which associated with adiposity- and metabolic phenotypes at genome-wide significant levels. Subsequently comparative analyses were performed to identify the extent of overlap between previously identified QTLs in both humans and pigs. The combined analysis of a large number of obesity phenotypes has provided insight in the genetic architecture of the molecular mechanisms underlying these traits indicating that QTLs underlying similar phenotypes are clustered in the genome. Our analyses have further confirmed that genetic heterogeneity is an inherent characteristic of obesity traits most likely caused by segregation or fixation of different variants of the individual components belonging to cellular pathways in different populations. Several important genes previously associated to obesity in human studies, along with novel genes were identified. Altogether, this study provides novel insight that may further the current understanding of the molecular mechanisms underlying human obesity.

Obesity is a complex condition that increases the risk of life threatening diseases such as cardiovascular disease and diabetes. Studying the gene regulation of obesity is important for understanding the molecular mechanisms behind the obesity derived diseases and may lead to better intervention and treatment plans. MicroRNAs (miRNAs) are short non-coding RNAs regulating target mRNA by binding to their 3'UTR. They are involved in numerous biological processes and diseases, including obesity. In this study we use a mixed breed pig model designed for obesity studies to investigate differentially expressed miRNAs in subcutaneous adipose tissue by RNA sequencing (RNAseq). Both male and female pigs are included to explore gender differences. The RNAseq study shows that the most highly expressed miRNAs are in accordance with comparable studies in pigs and humans. A total of six miRNAs are differentially expressed in subcutaneous adipose tissue between the lean and obese group of pigs, and in addition gender specific significant differential expression is observed for a number of miRNAs. The differentially expressed miRNAs have been verified using qPCR. The results of these studies in general confirm the trends found by RNAseq. Mir-9 and mir-124a are significantly differentially expressed with large fold changes in subcutaneous adipose tissue between lean and obese pigs. Mir-9 is more highly expressed in the obese pigs with a fold change of 10 and a p-value < 0.001. Mir-124a is more highly expressed in the obese pigs with a fold change of 114 and a p-value < 0.001. In addition, mir-124a is significantly higher expressed in abdominal adipose tissue in male pigs with a fold change of 119 and a p-value < 0.05. Both miRNAs are also significantly higher expressed in the liver of obese male pigs where mir-124a has a fold change of 12 and mir-9 has a fold change of 1.6, both with p-values < 0.05.

Background: Peripheral nerve sheath tumors (PNSTs) are frequently found in Danish cattle at slaughter. Bovine PNSTs share several gross and histopathological characteristics with the PNSTs in humans with heritable neurofibromatosis syndromes. The aim of the present study was to investigate a possible hereditary disposition to PNSTs in dairy cattle by statistical analysis performed on data from 567 cattle with PNSTs. Furthermore, a preliminary genome-wide association study (GWAS) was performed on DNA isolated from 28 affected and 28 non-affected Holstein cows to identify loci in the bovine genome involved in the development of PNSTs., Results: PNSTs were significantly more common in the Danish Holstein breed than in other breeds with 0.49% of Danish Holsteins slaughtered during an eight-year-period having PNSTs. PNSTs also occurred significantly more frequently in the offspring of some specific Holstein sires. Examination of three generation pedigrees showed that these sires were genetically related through a widely used US Holstein sire. The PNSTs included in GWAS were histologically classified as neurofibroma-schwannoma (43%), schwannoma (36%) and neurofibroma (21%) and derived from Holstein cows with multiple PNSTs. A single SNP on chromosome 27 reached genome-wide significance., Conclusions: Gross and histological characteristics of bovine PNSTs are comparable to PNSTs in humans (schwannomatosis). Danish Holsteins are genetically disposed to develop PNSTs but the examined materials are insufficient to allow determination of the mode of inheritance.

Obesity is a rising worldwide public health problem. Difficulties to precisely measure various obesity traits and the genetic heterogeneity in human have been major impediments to completely disentangle genetic factors causing obesity. The pig is a relevant model for studying human obesity and obesity-related (OOR) traits. Using founder breeds divergent with respect to obesity traits we have created an F2 pig resource population (454 pigs), which has been intensively phenotyped for 36 OOR traits. The main rationale for our study is to characterize the genetic architecture of OOR traits in the F2 pig design, by estimating heritabilities, genetic, and phenotypic correlations using mixed- and multi-trait BLUP animal models. Our analyses revealed high coefficients of variation (15-42%) and moderate to high heritabilities (0.22-0.81) in fatness traits, showing large phenotypic and genetic variation in the F2 population, respectively. This fulfills the purpose of creating a resource population divergent for OOR traits. Strong genetic correlations were found between weight and lean mass at dual-energy x-ray absorptiometry scanning (0.56-0.97). Weight and conformation also showed strong genetic correlations with slaughter traits (e.g., r g between abdominal circumference and leaf fat at slaughtering: 0.66). Genetic correlations between fat-related traits and the glucose level vary between 0.35 and 0.74 and show a strong correlation between adipose tissue and impaired glucose metabolism. Our power calculations showed a minimum of 80% power for QTL detection for all phenotypes. We revealed genetic correlations at population level, for the first time, for several difficult to measure and novel OOR traits and diseases. The results underpin the potential of the established F2 pig resource population for further genomic, systems genetics, and functional investigations to unravel the genetic background of OOR traits.

Our previously published second generation genetic map for the American mink (Neovison vison) has been used and redesigned in its best for genome-wide studies with maximum of efficiency. A number of 114 selected markers, including 33 newly developed microsatellite markers from the CHORI-231 mink Bacterial Artificial Chromosome (BAC) library, have been genotyped in a two generation population composed of 1200 individuals. The outcome reassigns the position of some markers on the chromosomes and it produces a more reliable map with a convenient distance between markers. A total of 104 markers mapped to 14 linkage groups corresponding to the mink autosomes. Six markers are unlinked and four markers are allocated to the X chromosome by homology but no linkage was detected. The sex-average linkage map spans 1192 centiMorgans (cM) with an average intermarker distance of 11.4cM and 1648cM when the ends of the linkage groups and the autosomal unlinked markers are added. Sex-specific genetic linkage maps were also generated. The male sex-specific map had a total length of 1014.6cM between the linked markers and an average inter-marker interval of 9.7cM. The female map has a corresponding length of 1378.6cM and an average inter-marker interval of 13.3cM. The study is complemented with additional anchorage for most of the chromosomes of the map by BAC in situ hybridization with clones containing microsatellites strategically selected from the various parts of the genome. This map provides an improved tool for genetic mapping and comparative genomics in mink, also useful for the future assembly of the mink genome sequence when this will be taken forward., (Copyright © 2012 Elsevier B.V. All rights reserved.)

Herniation of the intervertebral disk is a common cause of neurological dysfunction in the dog, particularly in the Dachshund. Using the Illumina CanineHD BeadChip, we have previously identified a major locus on canine chromosome 12 nucleotide positions 36,750,205-38,524,449 that strongly associates with intervertebral disk calcification in Danish wire-haired Dachshunds. In this study, targeted resequencing identified two synonymous variants in MB21D1 and one in the 5'-untranslated region of KCNQ5 that associates with intervertebral disk calcification in an independent sample of wire-haired Dachshunds. Haploview identified seven linkage disequilibrium blocks across the disease-associated region. The effect of haplotype windows on disk calcification shows that all haplotype windows are significantly associated with disk calcification. However, our predictions imply that the causal variant(s) are most likely to be found between nucleotide 36,750,205-37,494,845 as this region explains the highest proportion of variance in the dataset. Finally, we develop a risk prediction model for wire-haired Dachshunds. We validated the association of the chromosome 12 locus with disk calcification in an independent sample of wire-haired Dachshunds and identify potential risk variants. Additionally, we estimated haplotype effects and set up a model for prediction of disk calcifications in wire-haired Dachshunds based on genotype data. This genetic prediction model may prove useful in selection of breeding animals in future breeding programs.

Inflammatory diseases in mouse models are under strong impact from the gut microbiota. Therefore increased interindividual gut microbiota similarity may be seen as a way to reduce group sizes in mouse experiments. The composition of the gut microbiota is to a high extent defined by genetics, and it is known that selecting siblings as mothers even in inbred colonies may increase the gut microbiota similarity among the mice with 3-4%. We therefore hypothesized that selective breeding of mice aiming at a high similarity in the gut microbiota would increase the interindividual similarity of the gut microbiota. BALB/cCrl mice were, however, found to have a mean heterozygosity of only 0.8% in their genome, and selection of breeders with a high similarity in the gut microbiota for three generations did not change the overall gut microbiota similarity, which was 66% in the P generation and 66%, 64% and 63% in the F1, F2 and F3 generations, respectively. Increased gut microbiota similarity in closely related mice in inbred mouse colonies is, therefore, more likely to be caused by other factors, such as imprinting or different intrauterine conditions, rather than by residual heterozygosity.

Helminths almost invariably have an over-dispersed distribution in the host population. Human and animal studies have provided evidence suggesting that a large part of this variation is due to host genetic factors. Recently, the heritability for roundworm (Ascaris suum) infection levels in pigs was estimated to be 0.45. We used single nucleotide polymorphism markers to perform a whole-genome scan on 195 pigs experimentally infected with A. suum. A putative quantitative trait locus for worm burden on chromosome 4 covering 2.5 Mbp was identified by measured genotype analysis, although none of the SNPs reached genome-wide significance. To validate the putative quantitative trait locus, we genotyped two of the SNPs within the region in unrelated, informative animals exposed to experimental or natural infections and from which we had worm counts and/or faecal egg counts; the validation studies showed that one of the SNPs (TXNIP) was associated with total worm burden (P < 0.001) and adult worm burden(P < 0.0001), whereas the other SNP (ARNT) was associated with adult worm burden (P < 0.025) in these populations. We were thus able to confirm the existence of the quantitative trait locus on chromosome 4.This is to our knowledge the first report of a quantitative trait locus associated with helminth burden in pigs.

Epilepsy is the most common neurological disorder in dogs, with an incidence ranging from 0.5% to up to 20% in particular breeds. Canine epilepsy can be etiologically defined as idiopathic or symptomatic. Epileptic seizures may be classified as focal with or without secondary generalization, or as primary generalized. Nine genes have been identified for symptomatic (storage diseases) and one for idiopathic epilepsy in different breeds. However, the genetic background of common canine epilepsies remains unknown. We have studied the clinical and genetic background of epilepsy in Belgian Shepherds. We collected 159 cases and 148 controls and confirmed the presence of epilepsy through epilepsy questionnaires and clinical examinations. The MRI was normal while interictal EEG revealed abnormalities and variable foci in the clinically examined affected dogs. A genome-wide association study using Affymetrix 50K SNP arrays in 40 cases and 44 controls mapped the epilepsy locus on CFA37, which was replicated in an independent cohort (81 cases and 88 controls; combined p = 9.70×10⁻¹⁰, OR = 3.3). Fine mapping study defined a ∼1 Mb region including 12 genes of which none are known epilepsy genes or encode ion channels. Exonic sequencing was performed for two candidate genes, KLF7 and ADAM23. No variation was found in KLF7 but a highly-associated non-synonymous variant, G1203A (R387H) was present in the ADAM23 gene (p = 3.7×10⁻⁸, OR = 3.9 for homozygosity). Homozygosity for a two-SNP haplotype within the ADAM23 gene conferred the highest risk for epilepsy (p = 6.28×10⁻¹¹, OR = 7.4). ADAM23 interacts with known epilepsy proteins LGI1 and LGI2. However, our data suggests that the ADAM23 variant is a polymorphism and we have initiated a targeted re-sequencing study across the locus to identify the causative mutation. It would establish the affected breed as a novel therapeutic model, help to develop a DNA test for breeding purposes and introduce a novel candidate gene for human idiopathic epilepsies.

Fertility is one of the most important traits in dairy cattle, and has been steadily declining over the last decades. We herein use state-of-the-art genomic tools, including high-throughput SNP genotyping and next-generation sequencing, to identify a 3.3 Kb deletion in the FANCI gene causing the brachyspina syndrome (BS), a rare recessive genetic defect in Holstein dairy cattle. We determine that despite the very low incidence of BS (<1/100,000), carrier frequency is as high as 7.4% in the Holstein breed. We demonstrate that this apparent discrepancy is likely due to the fact that a large proportion of homozygous mutant calves die during pregnancy. We postulate that several other embryonic lethals may segregate in livestock and significantly compromise fertility, and propose a genotype-driven screening strategy to detect the corresponding deleterious mutations.

The aim of the study was to investigate the different phases of the immune response after DNA immunization with the hemagglutinin and nucleoprotein genes from canine distemper virus (CDV). Although attenuated live CDV vaccines have effectively reduced the incidence of disease, canine distemper is still a problem worldwide. The broad host range of CDV creates a constant viral reservoir among wildlife animals. Our results demonstrated early humoral and cell-mediated immune responses (IFN-gamma) in DNA vaccinated mink compared to mock-vaccinated mink after challenge with a Danish wild-type CDV. The DNA vaccine-induced immunity protected the natural host against disease development.