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.

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.

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.

Background: Meat quality traits are important in pig breeding programs, but they are difficult to include in a traditional selection program. Marker assisted selection (MAS) of meat quality traits is therefore of interest in breeding programs and a Quantitative Trait Locus (QTL) analysis is the key to identifying markers that can be used in MAS. In this study, Landrace and Hampshire intercross and backcross families were used to investigate meat quality traits. Hampshire pigs are commonly used as the sire line in commercial pig breeding. This is the first time a pedigree including Hampshire pigs has been used for a QTL analysis of meat quality traits., Results: In total, we analyzed 39 meat quality traits and identified eight genome-wide significant QTL peaks in four regions: one on chromosome 3, two on chromosome 6 and one on chromosome 16. At least two of the QTLs do not appear to have been detected in previous studies. On chromosome 6 we identified QTLs for water content in M. longissimus dorsi (LD), drip loss in LD and post mortem pH decline in LD. On chromosomes 3 and 16 we identified previously undetected QTLs for protein content in LD and for freezing and cooking loss respectively., Conclusion: We identified at least two new meat quality trait QTLs at the genome-wide significance level. We detected two QTLs on chromosome 6 that possibly coincide with QTLs detected in other studies. We were also able to exclude the C1843T mutation in the ryanodine receptor (RYR1) as a causative mutation for one of the chromosome 6 QTLs in this cross.

We report on the construction of a high-resolution comparative map of porcine chromosome 17 (SSC17) focusing on evolutionary breakpoints with human chromosomes. The comparative map shows high homology with human chromosome 20 but suggests more limited homologies with other human chromosomes. SSC17 is of particular interest in studies of chromosomal organization due to the presence of QTLs that affect meat quality and carcass composition. A total of 158 pig ESTs available in databases or developed by the Sino-Danish Pig Genome Sequencing Consortium were mapped using the INRA-University of Minnesota porcine radiation hybrid panel. The high-resolution map was further anchored by fluorescence in situ hybridization. This study confirmed the extensive conservation between SSC17 and HSA20 and enabled the gene order to be determined. The homology of the SSC17 pericentromeric region was extended to other human chromosomes (HSA4, HSA8) and the chromosomal breakpoint boundaries were accurately defined. In total 15 breakpoints were identified.