Molecular advances in QTL discovery and application in pig breeding

The Identification of Goat KCNJ15 Gene Copy Number Variation and Its Association with Growth Traits

(1) Background: Copy number variation (CNV) is a critical component of genome structural variation and has garnered significant attention. High-throughput screening of the KCNJ15 gene has revealed a correlation between the CNV region and the growth traits of goats. We aimed to identify the CNV of the KCNJ15 gene in five goat breeds and analyze its association with growth characteristics. (2) Methods: We utilized 706 goats from five breeds: Guizhou black goat (GZB), Guizhou white goat (GZW), Bohuai goat (BH), Huai goat (HH), and Taihang goat (TH). To evaluate the number of copies of the KCNJ15 gene using qPCR, we analyzed the correlation between the CNV and growth characteristics and then used a universal linear model. The findings revealed variations in the distribution of different copy number types among the different goat breeds. (3) Results: Association analysis revealed a positive influence of the CNV in the KCNJ15 gene on goat growth. In GZB, individuals with duplication types exhibited superior performance in terms of cannon bone circumference (p < 0.05). In HH, individuals with duplication types exhibited superior performance in terms of body slanting length (p < 0.05). Conversely, normal TH demonstrated better body height and body weight (p < 0.05), while in GZW, when CN = 3, it performed better than other types in terms of body weight and chest circumference (p < 0.05). However, in BH, it had no significant effect on growth traits. (4) Conclusions: We confirmed that the CNV in the KCNJ15 gene significantly influences the growth characteristics of four distinct goat breeds. The correlation between KCNJ15 gene CNVs and goat growth traits offers valuable insights to breeders, enabling them to employ precise and efficient breeding methods that enhance livestock welfare, productivity, and overall economic benefits in the industry.

Leveraging molecular quantitative trait loci to comprehend complex diseases/traits from the omics perspective

Comprehending the molecular basis of quantitative genetic variation is a principal goal for complex diseases or traits. Molecular quantitative trait loci (molQTLs) have made it possible to investigate the effects of genetic variants hiding behind large-scale omics data. A deeper understanding of molQTL is urgently required in light of the multi-dimensionalization of omics data to more fully elucidate the pertinent biological mechanisms. Herein, we reviewed molQTLs with the corresponding resource from the omics perspective and further discussed the integrative strategy of GWAS-molQTL to infer their causal effects. Subsequently, we described the opportunities and challenges encountered by molQTL. The case studies showed that molQTL is essential for complex diseases and traits, whether single- or multi-omics QTLs. Overall, we highlighted the functional significance of genetic variants to employ the discovery of molQTL in complex diseases and traits.

Identification of gut microbes associated with feed efficiency by daily-phase feeding strategy in growing-finishing pigs

Feed efficiency is one of the most important issues for sustainable pig production. Daily-phase feeding (DPF) is a form of precision feeding that could improve feed efficiency in pigs. Gut microbiota can regulate host nutrient digestion, absorption, and metabolism. However, which key microbes may play a vital role in improving the feed efficiency during DPF remains unclear. In the present study, we used a DPF program compared to a three-phase feeding (TPF) program in growing-finishing pigs to investigate the effects of gut microbiota on feed efficiency. A total of 204 Landrace × Yorkshire pigs (75 d) were randomly assigned into 2 treatments. Each treatment was replicated 8 times with 13 to 15 pigs per replicate pen. Pigs in the TPF group were fed with a commercial feeding program that supplied fixed feed for phases I, II, and III, starting at 81, 101, and 132 d of age, respectively, and pigs in the DPF group were fed a blend of adjacent phase feed from 81 to 155 d at a gradual daily ratio and phase III feed from 155 to 180 d of age. Daily feed intake and body weight were recorded by a computerized device in the feeders. Feces and blood samples were collected from 1 pig per replicate at 155 and 180 d of age. The results showed that the DPF program remarkably improved the feed efficiency at 155 d (P < 0.001) and 180 d of age (P < 0.001), with a significant reduction of the intake of crude protein (P < 0.01), net energy (P < 0.001), crude fiber (P < 0.001), ether extract (P < 0.01), and ash (P < 0.001). The daily-phase feeding program increased the abundance of Prevotella copri (P < 0.05) and Paraprevotella clara (P < 0.05), while it decreased the abundance of Ocilibacter (P < 0.05) at 155 d of age. The results of correlation analysis indicated that the differentially abundant microbiota communities were closely associated with 20 metabolites which enriched amino acid and phenylalanine metabolism. Our results suggest that 2 key microbes may contribute to feed efficiency during daily-phase feeding strategies in pigs.

A composite strategy of genome-wide association study and copy number variation analysis for carcass traits in a Duroc pig population

BACKGROUND

Carcass traits are important in pig breeding programs for improving pork production. Understanding the genetic variants underlies complex phenotypes can help explain trait variation in pigs. In this study, we integrated a weighted single-step genome-wide association study (wssGWAS) and copy number variation (CNV) analyses to map genetic variations and genes associated with loin muscle area (LMA), loin muscle depth (LMD) and lean meat percentage (LMP) in Duroc pigs.

RESULTS

Firstly, we performed a genome-wide analysis for CNV detection using GeneSeek Porcine SNP50 Bead chip data of 3770 pigs. A total of 11,100 CNVs were detected, which were aggregated by overlapping 695 CNV regions (CNVRs). Next, we investigated CNVs of pigs from the same population by whole-genome resequencing. A genome-wide analysis of 21 pigs revealed 23,856 CNVRs that were further divided into three categories (851 gain, 22,279 loss, and 726 mixed), which covered 190.8 Mb (~ 8.42%) of the pig autosomal genome. Further, the identified CNVRs were used to determine an overall validation rate of 68.5% for the CNV detection accuracy of chip data. CNVR association analyses identified one CNVR associated with LMA, one with LMD and eight with LMP after applying stringent Bonferroni correction. The wssGWAS identified eight, six and five regions explaining more than 1% of the additive genetic variance for LMA, LMD and LMP, respectively. The CNVR analyses and wssGWAS identified five common regions, of which three regions were associated with LMA and two with LMP. Four genes (DOK7, ARAP1, ELMO2 and SLC13A3) were highlighted as promising candidates according to their function.

CONCLUSIONS

We determined an overall validation rate for the CNV detection accuracy of low-density chip data and constructed a genomic CNV map for Duroc pigs using resequencing, thereby proving a value genetic variation resource for pig genome research. Furthermore, our study utilized a composite genetic strategy for complex traits in pigs, which will contribute to the study for elucidating the genetic architecture that may be influenced and regulated by multiple forms of variations.

Genes Related to Fat Metabolism in Pigs and Intramuscular Fat Content of Pork: A Focus on Nutrigenetics and Nutrigenomics

Simple Summary

The intramuscular fat (IMF) or marbling is an essential pork sensory quality that influences the preference of the consumers and premiums for pork. IMF is the streak of visible fat intermixed with the lean within a muscle fibre and determines sensorial qualities of pork such as flavour, tenderness and juiciness. Fat metabolism and IMF development are controlled by dietary nutrients, genes, and their metabolic pathways in the pig. Nutrigenetics explains how the genetic make-up of an individual pig influences the pig’s response to dietary nutrient intake. Differently, nutrigenomics is the analysis of how the entire genome of an individual pig is affected by dietary nutrient intake. The knowledge of nutrigenetics and nutrigenomics, when harmonized, is a powerful tool in estimating nutrient requirements for swine and programming dietary nutrient supply according to an individual pig’s genetic make-up. The current paper aimed to highlight the roles of nutrigenetics and nutrigenomics in elucidating the underlying mechanisms of fat metabolism and IMF deposition in pigs. This knowledge is essential in redefining nutritional intervention for swine production and the improvement of some economically important traits such as growth performance, backfat thickness, IMF accretion, disease resistance etc., in animals.

Abstract

Fat metabolism and intramuscular fat (IMF) are qualitative traits in pigs whose development are influenced by several genes and metabolic pathways. Nutrigenetics and nutrigenomics offer prospects in estimating nutrients required by a pig. Application of these emerging fields in nutritional science provides an opportunity for matching nutrients based on the genetic make-up of the pig for trait improvements. Today, integration of high throughput “omics” technologies into nutritional genomic research has revealed many quantitative trait loci (QTLs) and single nucleotide polymorphisms (SNPs) for the mutation(s) of key genes directly or indirectly involved in fat metabolism and IMF deposition in pigs. Nutrient–gene interaction and the underlying molecular mechanisms involved in fatty acid synthesis and marbling in pigs is difficult to unravel. While existing knowledge on QTLs and SNPs of genes related to fat metabolism and IMF development is yet to be harmonized, the scientific explanations behind the nature of the existing correlation between the nutrients, the genes and the environment remain unclear, being inconclusive or lacking precision. This paper aimed to: (1) discuss nutrigenetics, nutrigenomics and epigenetic mechanisms controlling fat metabolism and IMF accretion in pigs; (2) highlight the potentials of these concepts in pig nutritional programming and research.

miR-152 targets pyruvate kinase to regulate the glycolytic activity of pig skeletal muscles and affects pork quality

As a type of non-coding RNA, microRNAs are widely involved in the biological processes of animals. In the present study, the expression of miR-152 in glycolytic muscle fibers (Longissimus thoracis, LT) was lower than that of oxidative muscle fibers (Psoas major, PM). Using dual luciferase assay, miR-152 was shown to target muscle pyruvate kinase (PKM) to perform biological functions. Moreover, overexpression of miR-152 in primary porcine cells inhibited PKM gene expression and reduced lactic acid production in cells, whereas inhibition of miR-152 expression promoted PKM gene expression and increased lactic acid production. Correlation analysis showed that the expression of miR-152 was significantly positively correlated with the ultimate pH of LT after slaughter, while the expression of the PKM gene was significantly negatively correlated with the final pH of LT. In vivo and in vitro experiments discussed herein suggest that miR-152 may affect muscle pH by targeting the expression of the PKM gene. Our findings enrich the understanding of the genetic regulatory network that influences pork quality.

Evidence for and localization of proposed causative variants in cattle and pig genomes

BACKGROUND

This paper reviews the localization of published potential causative variants in contemporary pig and cattle reference genomes, and the evidence for their causality. In spite of the difficulties inherent to the identification of causative variants from genetic mapping and genome-wide association studies, researchers in animal genetics have proposed putative causative variants for several traits relevant to livestock breeding.

RESULTS

For this review, we read the literature that supports potential causative variants in 13 genes (ABCG2, DGAT1, GHR, IGF2, MC4R, MSTN, NR6A1, PHGK1, PRKAG3, PLRL, RYR1, SYNGR2 and VRTN) in cattle and pigs, and localized them in contemporary reference genomes. We review the evidence for their causality, by aiming to separate the evidence for the locus, the proposed causative gene and the proposed causative variant, and report the bioinformatic searches and tactics needed to localize the sequence variants in the cattle or pig genome.

CONCLUSIONS

Taken together, there is usually good evidence for the association at the locus level, some evidence for a specific causative gene at eight of the loci, and some experimental evidence for a specific causative variant at six of the loci. We recommend that researchers who report new potential causative variants use referenced coordinate systems, show local sequence context, and submit variants to repositories.

A further survey of the quantitative trait loci affecting swine body size and carcass traits in five related pig populations

Breeding for good meat quality performance while maintaining large body size and desirable carcass traits has been the major challenge for modern swine selective breeding. To address this goal, in the present work we studied five related populations produced by two commercial breeds (Berkshire and Duroc) and two Chinese breeds (Licha black pig and Lulai black pig). A single-trait GWAS performed on 20 body size and carcass traits using a self-developed China Chip-1 porcine SNP50K BeadChip identified 11 genome-wide significant QTL on nine chromosomes and 22 suggestive QTL on 15 chromosomes. For the 11 genome-wide significant QTL, eight were detected in at least two populations, and the rest were population-specific and only mapped in Shanxia black pig. Most of the genome-wide significant QTL were pleiotropic; for example, the QTL around 75.65 Mb on SSC4 was associated with four traits at genome-wide significance level. After screening the genes within 50 kb of the top SNP for each genome-wide significant QTL, NR6A1 and VRTN were chosen as candidate genes for vertebrae number; PLAG1 and BMP2 were identified as candidate genes for body size; and MC4R was the strong candidate gene for body weight. The four genes have been reported as candidates for thoracic vertebrae number, lumbar vertebrae number, carcass length and body weight respectively in previous studies. The effects of VRTN on thoracic vertebrae number, carcass length and body length have been verified in Shanxia black pig. Therefore, the VRTN genotype could be used in gene-assisted selection, and this could accelerate genetic improvement of body size and carcass traits in Shanxia black pig.

[Whole-genome association studies of distribution of developmental abnormalities and other breeding-valuable qualitative traits in offspring of the Russian large-white boars]

Выявление областей генома, прямо или опосредованно связанных с признаками пороков развития у домашних свиней, может способствовать идентификации генетических мишеней, используемых в качестве биомаркеров индивидуальных особенностей формирования экстерьера, их метаболического статуса, а также подверженности генетическим заболеваниям. Такие исследования напрямую связаны с повышением экономической эффективности, поскольку позволяют выявлять и исключать из селекционного процесса животных-носителей нежелательных генов, фенотип которых может не проявляться. В данной работе проведен поиск подобных целевых генов и геномных регионов с помощью полногеномных ассоциативных исследований (GWAS) с использованием ДНК-чипов PorcineSNP60K BeadChips (Illumina, San Diego, USA). Проанализировано 48 хряков свиней крупной белой породы селекционно-гибридного центра «Знаменский» Орловской области по 21 недостатку экстерьера и дефектам развития у 39 153 их потомков. Расчеты производили по линейной модели смешанного типа в программном пакете GEMMA. Из изначального сета в 61 000 SNP были отобраны 36 704 полиморфных SNP, в которых найдены 24 полиморфизма, входящих в 11 генов (P < 0.1), статистически значимо коррелирующих с признаками аномалий развития в геноме свиней, такими как атрезия анального отверстия (ARMC7, FANCC, RND3, ENSSSCG00000017216), проблемы с конечностями (PAWR, NTM, OPCML, ENSSSCG00000040250, ENSSSCG00000017018) и тремор поросят (RIC3, ENSSSCG00000032665). Также была выявлена коэкспрессия генов NTM, OPCML и RND3, участвующих в регуляции клеточной адгезии. Проведенная работа подтвердила актуальность применения подобного подхода в полногеномно-ассоциативных исследованиях для детектирования единичных SNP, связанных с отдельными признаками, даже для небольших выборок.

The Association of an SNP in the EXOC4 Gene and Reproductive Traits Suggests Its Use as a Breeding Marker in Pigs

In mammals, the exocyst complex component 4 (EXOC4) gene has often been reported to be involved in vesicle transport. The SNP rs81471943 (C/T) is located in the intron of porcine EXOC4, while six quantitative trait loci (QTL) within 5-10 Mb around EXOC4 are associated with ovary weight, teat number, total offspring born alive, and corpus luteum number. However, the molecular mechanisms between EXOC4 and the reproductive performance of pigs remains to be elucidated. In this study, rs81471943 was genotyped from a total of 994 Duroc sows, and the genotype and allele frequency of SNP rs81471943 (C/T) were statistically analyzed. Then, the associations between SNP rs81471943 and four reproductive traits, including number of piglets born alive (NBA), litter weight at birth (LWB), number of piglets weaned (NW), and litter weight at weaning (LWW), were determined. Sanger sequencing and PCR restriction fragment length polymorphism (PCR-RFLP) were utilized to identify the rs81471943 genotype. We found that the genotype frequency of CC was significantly higher than that of CT and TT, and CC was the most frequent genotype for NBA, LWB, NW, and LWW. Moreover, 5'-deletion and luciferase assays identified a positive transcription regulatory element in the EXOC4 promoter. After exploring the EXOC4 promoter, SNP -1781G/A linked with SNP rs81471943 (C/T) were identified by analysis of the transcription activity of the haplotypes, and SNP -1781 G/A may influence the potential binding of P53, E26 transformation specific sequence -like 1 transcription factor (ELK1), and myeloid zinc finger 1 (MZF1). These findings provide useful information for identifying a molecular marker of EXOC4-assisted selection in pig breeding.

The Effect of Dietary Lycopene Supplementation on Drip Loss during Storage of Lamb Meat by iTRAQ Analysis

This study was designed to investigate the impact of dietary lycopene (antioxidant extracted from tomato) supplementation on postmortem antioxidant capacity, drip loss and protein expression profiles of lamb meat during storage. Thirty male Hu lambs were randomly divided into three treatment groups and housed in individual pens and received 0, 200 or 400 mg·kg⁻¹ lycopene in their diet, respectively. All lambs were slaughtered after 3 months of fattening, and the longissimus thoracis (LT) muscle was collected for analyses. The results indicated that drip loss of LT muscle increased with storage days (P < 0.05). After storage for 7 days, significantly lower drip loss of meat was found in fed the lycopene-supplemented diet (P < 0.05). Dietary lycopene supplementation increased the activity of antioxidant enzymes (total antioxidant capacity (T-AOC), superoxide dismutase (SOD), glutathione peroxidase (GSH-Px), catalase (CAT)) (P < 0.05) and decreased the thiobarbituric acid reactive substance (TBARS) and carbonyl contents (P < 0.05). During the storage period (days 0, 5 and 7), a number of differentially abundant proteins (DAPs), including oxidases, metabolic enzymes, calcium channels and structural proteins, were identified based on iTRAQ data, with roles predominantly in carbon metabolism, oxidative phosphorylation, cardiac muscle contraction and proteasome pathways, and which contribute to decreased drip loss of lamb meat during storage. It can be concluded that dietary lycopene supplementation increased antioxidant capacity after slaughter, and the decreased drip loss during postmortem storage might occur by changing the expression of proteins related to enzyme activity and cellular structure in lamb muscle.

GALLO: An R package for genomic annotation and integration of multiple data sources in livestock for positional candidate loci

BACKGROUND

The development of high-throughput sequencing and genotyping methodologies has enabled the identification of thousands of genomic regions associated with several complex traits. The integration of multiple sources of biological information is a crucial step required to better understand patterns regulating the development of these traits.

FINDINGS

Genomic Annotation in Livestock for positional candidate LOci (GALLO) is an R package developed for the accurate annotation of genes and quantitative trait loci (QTLs) located in regions identified in common genomic analyses performed in livestock, such as genome-wide association studies and transcriptomics using RNA sequencing. Moreover, GALLO allows the graphical visualization of gene and QTL annotation results, data comparison among different grouping factors (e.g., methods, breeds, tissues, statistical models, studies), and QTL enrichment in different livestock species such as cattle, pigs, sheep, and chickens.

CONCLUSIONS

Consequently, GALLO is a useful package for annotation, identification of hidden patterns across datasets, and data mining previously reported associations, as well as the efficient examination of the genetic architecture of complex traits in livestock.

Genetic parameters and associated genomic regions for global immunocompetence and other health-related traits in pigs

The inclusion of health-related traits, or functionally associated genetic markers, in pig breeding programs could contribute to produce more robust and disease resistant animals. The aim of the present work was to study the genetic determinism and genomic regions associated to global immunocompetence and health in a Duroc pig population. For this purpose, a set of 30 health-related traits covering immune (mainly innate), haematological, and stress parameters were measured in 432 healthy Duroc piglets aged 8 weeks. Moderate to high heritabilities were obtained for most traits and significant genetic correlations among them were observed. A genome wide association study pointed out 31 significantly associated SNPs at whole-genome level, located in six chromosomal regions on pig chromosomes SSC4, SSC6, SSC17 and SSCX, for IgG, γδ T-cells, C-reactive protein, lymphocytes phagocytic capacity, total number of lymphocytes, mean corpuscular volume and mean corpuscular haemoglobin. A total of 16 promising functionally-related candidate genes, including CRP, NFATC2, PRDX1, SLA, ST3GAL1, and VPS4A, have been proposed to explain the variation of immune and haematological traits. Our results enhance the knowledge of the genetic control of traits related with immunity and support the possibility of applying effective selection programs to improve immunocompetence in pigs.

Phenotyping of the Visceral Adipose Tissue Using Dual Energy X-Ray Absorptiometry (DXA) and Magnetic Resonance Imaging (MRI) in Pigs

The objective of this study was to phenotype visceral adipose tissue (VAT) in pigs. In this context, the ability to detect VAT by using the DXA CoreScan mode within the enCORE software, version 17 (GE Healthcare) was evaluated in comparison with MRI measurements (Siemens Magnetom C!) of the same body region. A number of 120 crossbred pigs of the F1 and F2 generation, with the parental breeds Large White, Landrace, Piétrain, and Duroc, were examined at an age of 150 days. A whole-body scan in two different modes ("thick", "standard") was carried out by a GE Lunar iDXA scanner. Very strong relationships (R² = 0.95, RMSE = 175cm³) were found for VAT between the two DXA modes. The comparison of VAT measured by MRI and DXA shows high linear relationships ("thick": R² = 0.76, RMSE = 399.25cm³/"standard": R² = 0.71, RMSE = 443.42cm³), but is biased, according to the Bland-Altman analysis. A variance analysis of VAT shows significant differences for both DXA modes and for MRI between male and female pigs, as well as between F1 and F2. In conclusion, DXA "CoreScan" has the ability to estimate VAT in pigs with a close relationship to MRI but needs bias correction.

Genome-wide association analyses identify known and novel loci for teat number in Duroc pigs using single-locus and multi-locus models

BACKGROUND

More teats are necessary for sows to nurse larger litters to provide immunity and nutrient for piglets prior to weaning. Previous studies have reported the strong effect of an insertion mutation in the Vertebrae Development Associated (VRTN) gene on Sus scrofa chromosome 7 (SSC7) that increased the number of thoracic vertebrae and teat number in pigs. We used genome-wide association studies (GWAS) to map genetic markers and genes associated with teat number in two Duroc pig populations with different genetic backgrounds. A single marker method and several multi-locus methods were utilized. A meta-analysis that combined the effects and P-values of 34,681 single nucleotide polymorphisms (SNPs) that were common in the results of single marker GWAS of American and Canadian Duroc pigs was conducted. We also performed association tests between the VRTN insertion and teat number in the same populations.

RESULTS

A total of 97 SNPs were found to be associated with teat number. Among these, six, eight and seven SNPs were consistently detected with two, three and four multi-locus methods, respectively. Seven SNPs were concordantly identified between single marker and multi-locus methods. Moreover, 26 SNPs were newly found by multi-locus methods to be associated with teat number. Notably, we detected one consistent quantitative trait locus (QTL) on SSC7 for teat number using single-locus and meta-analysis of GWAS and the top SNP (rs692640845) explained 8.68% phenotypic variance of teat number in the Canadian Duroc pigs. The associations between the VRTN insertion and teat number in two Duroc pig populations were substantially weaker. Further analysis revealed that the effect of VRTN on teat number may be mediated by its LD with the true causal mutation.

CONCLUSIONS

Our study suggested that VRTN insertion may not be a strong or the only candidate causal mutation for the QTL on SSC7 for teat number in the analyzed Duroc pig populations. The combination of single-locus and multi-locus GWAS detected additional SNPs that were absent using only one model. The identified SNPs will be useful for the genetic improvement of teat number in pigs by assigning higher weights to associated SNPs in genomic selection.

DNA N6-Methyladenosine modification role in transmitted variations from genomic DNA to RNA in Herrania umbratica

Background

DNA methylation is an important epigenetic modification. Recently the developed single-molecule real-time (SMRT) sequencing technology provided an efficient way to detect DNA N6-methyladenine (6mA) modification that played an important role in epigenetic and positively regulated gene expression. In addition, the gene expression was also regulated by genetic variation. However, the relationship between DNA 6mA modification and variation is still unknown.

Results

We collected the SMRT long-reads DNA, Illumina short reads DNA and RNA datasets from the young leaves of Herrania umbratica, and used them to detect 35,654 6mA modification sites, 829,894 DNA variations and 60,672 RNA variations respectively, among which, there are 303 DNA variations and 19 RNA variations with 6mA modification, and 57,468 transmitted genetic variations from DNA to RNA. The results illustrated that the genes with 6mA modification were significant disadvantage to mutate than those genes without modification (p-value< 4.9e-08). And result from the linear regression model showed the 6mA densities of genes were associated with the transmitted variations type 0/1 to 1/1 (p-value < 0.001).

Conclusions

The variations of DNA and RNA in genes with 6mA modification were significant less than those in unmodified genes. Furthermore, the variations in 6mA modified genes were easily transmitted from DNA to RNA, especially the transmitted variation from DNA heterozygote to RNA homozygote.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5776-0) contains supplementary material, which is available to authorized users.

Genetic control of longissimus dorsi muscle gene expression variation and joint analysis with phenotypic quantitative trait loci in pigs

Background

Economically important growth and meat quality traits in pigs are controlled by cascading molecular events occurring during development and continuing throughout the conversion of muscle to meat. However, little is known about the genes and molecular mechanisms involved in this process. Evaluating transcriptomic profiles of skeletal muscle during the initial steps leading to the conversion of muscle to meat can identify key regulators of polygenic phenotypes. In addition, mapping transcript abundance through genome-wide association analysis using high-density marker genotypes allows identification of genomic regions that control gene expression, referred to as expression quantitative trait loci (eQTL). In this study, we perform eQTL analyses to identify potential candidate genes and molecular markers regulating growth and meat quality traits in pigs.

Results

Messenger RNA transcripts obtained with RNA-seq of longissimus dorsi muscle from 168 F2 animals from a Duroc x Pietrain pig resource population were used to estimate gene expression variation subject to genetic control by mapping eQTL. A total of 339 eQTL were mapped (FDR ≤ 0.01) with 191 exhibiting local-acting regulation. Joint analysis of eQTL with phenotypic QTL (pQTL) segregating in our population revealed 16 genes significantly associated with 21 pQTL for meat quality, carcass composition and growth traits. Ten of these pQTL were for meat quality phenotypes that co-localized with one eQTL on SSC2 (8.8-Mb region) and 11 eQTL on SSC15 (121-Mb region). Biological processes identified for co-localized eQTL genes include calcium signaling (FERM, MRLN, PKP2 and CHRNA9), energy metabolism (SUCLG2 and PFKFB3) and redox hemostasis (NQO1 and CEP128), and results support an important role for activation of the PI3K-Akt-mTOR signaling pathway during the initial conversion of muscle to meat.

Conclusion

Co-localization of eQTL with pQTL identified molecular markers significantly associated with both economically important phenotypes and gene transcript abundance. This study reveals candidate genes contributing to variation in pig production traits, and provides new knowledge regarding the genetic architecture of meat quality phenotypes.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5386-2) contains supplementary material, which is available to authorized users.

Genome-wide association studies and meta-analysis uncovers new candidate genes for growth and carcass traits in pigs

Genome-wide association studies (GWAS) have been widely used in the genetic dissection of complex traits. As more genomic data is being generated within different commercial or resource pig populations, the challenge which arises is how to collectively investigate the data with the purpose to increase sample size and implicitly the statistical power. This study performs an individual population GWAS, a joint population GWAS and a meta-analysis in three pig F2 populations. D1 is derived from European type breeds (Piétrain, Large White and Landrace), D2 is obtained from an Asian breed (Meishan) and Piétrain, and D3 stems from a European Wild Boar and Piétrain, which is the common founder breed. The traits investigated are average daily gain, backfat thickness, meat to fat ratio and carcass length. The joint and the meta-analysis did not identify additional genomic clusters besides the ones discovered via the individual population GWAS. However, the benefit was an increased mapping resolution which pinpointed to narrower clusters harboring causative variants. The joint analysis identified a higher number of clusters as compared to the meta-analysis; nevertheless, the significance levels and the number of significant variants in the meta-analysis were generally higher. Both types of analysis had similar outputs suggesting that the two strategies can complement each other and that the meta-analysis approach can be a valuable tool whenever access to raw datasets is limited. Overall, a total of 20 genomic clusters that predominantly overlapped at various extents, were identified on chromosomes 2, 7 and 17, many confirming previously identified quantitative trait loci. Several new candidate genes are being proposed and, among them, a strong candidate gene to be taken into account for subsequent analysis is BMP2 (bone morphogenetic protein 2).

The effect of QTL-rich region polymorphisms identified by targeted DNA-seq on pig production traits

The aim of the present study was to analyse the effect of PLCD4, PECR, FN1 and PNKD mutations on pig productive traits and tested the usefulness of targeted enrichment DNA sequencing method as tool for preselection of genetic markers. The potential genetic markers for pig productive traits were identified by using targeted enrichment DNA sequencing of chromosome 15 region that is QTL-rich. The selected mutations were genotyped by using HRM, Sanger sequencing and PCR-ACRS methods. The association study was performed by using GLM model in SAS program and included over 500 pigs of 5 populations maintained in Poland. The variation (C/T) of PLCD4 gene affected feed conversion, intramuscular fat and water exudation. The PNKD mutations were associated with texture parameters measured after cooking. In turn, the variation rs792423408 (C/T) in the FN1 gene influenced toughness measured in semimembranosus muscle and growth traits that was observed particularly in Duroc breed. Summarizing, the investigated gene variants delivered valuable information that could be used during developing SNP microarray for genomic estimated breeding value procedure in pigs. Moreover, the study showed that the TEDNA-seq method could be used to preselect the molecular markers associated with pig traits. However, the further association study that included large number animal populations is necessary.

Systems genomics study reveals expression quantitative trait loci, regulator genes and pathways associated with boar taint in pigs

Boar taint is an offensive odour and/or taste from a proportion of non-castrated male pigs caused by skatole and androstenone accumulation during sexual maturity. Castration is widely used to avoid boar taint but is currently under debate because of animal welfare concerns. This study aimed to identify expression quantitative trait loci (eQTLs) with potential effects on boar taint compounds to improve breeding possibilities for reduced boar taint. Danish Landrace male boars with low, medium and high genetic merit for skatole and human nose score (HNS) were slaughtered at ~100 kg. Gene expression profiles were obtained by RNA-Seq, and genotype data were obtained by an Illumina 60K Porcine SNP chip. Following quality control and filtering, 10,545 and 12,731 genes from liver and testis were included in the eQTL analysis, together with 20,827 SNP variants. A total of 205 and 109 single-tissue eQTLs associated with 102 and 58 unique genes were identified in liver and testis, respectively. By employing a multivariate Bayesian hierarchical model, 26 eQTLs were identified as significant multi-tissue eQTLs. The highest densities of eQTLs were found on pig chromosomes SSC12, SSC1, SSC13, SSC9 and SSC14. Functional characterisation of eQTLs revealed functions within regulation of androgen and the intracellular steroid hormone receptor signalling pathway and of xenobiotic metabolism by cytochrome P450 system and cellular response to oestradiol. A QTL enrichment test revealed 89 QTL traits curated by the Animal Genome PigQTL database to be significantly overlapped by the genomic coordinates of cis-acting eQTLs. Finally, a subset of 35 cis-acting eQTLs overlapped with known boar taint QTL traits. These eQTLs could be useful in the development of a DNA test for boar taint but careful monitoring of other overlapping QTL traits should be performed to avoid any negative consequences of selection.

Squalene epoxidase plays a critical role in determining pig meat quality by regulating adipogenesis, myogenesis, and ROS scavengers

In mammals, Squalene epoxidase (SQLE) is an enzyme that converts squalene to 2,3-oxidosqualene, in the early stage of cholesterol generation. Here, we identified single nucleotide polymorphisms (SNPs) in the SQLE gene (c.2565 G > T) by RNA Sequencing from the liver tissue of Berkshire pigs. Furthermore, we found that homozygous GG pigs expressed more SQLE mRNA than GT heterozygous and TT homozygous pigs in longissimus dorsi tissue. Next, we showed that the SNP in the SQLE gene was associated with several meat quality traits including backfat thickness, carcass weight, meat colour (yellowness), fat composition, and water-holding capacity. Rates of myogenesis and adipogenesis induced in C2C12 cells and 3T3-L1 cells, respectively, were decreased by Sqle knockdown. Additionally, the expression of myogenic marker genes (Myog, Myod, and Myh4) and adipogenic marker genes (Pparg, Cebpa, and Adipoq) was substantially downregulated in cells transfected with Sqle siRNA. Moreover, mRNA expression levels of ROS scavengers, which affect meat quality by altering protein oxidation processes, were significantly downregulated by Sqle knockdown. Taken together, our results suggest the molecular mechanism by which SNPs in the SQLE gene can affect meat quality.

Pilot study of large-scale production of mutant pigs by ENU mutagenesis

N-ethyl-N-nitrosourea (ENU) mutagenesis is a powerful tool to generate mutants on a large scale efficiently, and to discover genes with novel functions at the whole-genome level in Caenorhabditis elegans, flies, zebrafish and mice, but it has never been tried in large model animals. We describe a successful systematic three-generation ENU mutagenesis screening in pigs with the establishment of the Chinese Swine Mutagenesis Consortium. A total of 6,770 G1 and 6,800 G3 pigs were screened, 36 dominant and 91 recessive novel pig families with various phenotypes were established. The causative mutations in 10 mutant families were further mapped. As examples, the mutation of SOX10 (R109W) in pig causes inner ear malfunctions and mimics human Mondini dysplasia, and upregulated expression of FBXO32 is associated with congenital splay legs. This study demonstrates the feasibility of artificial random mutagenesis in pigs and opens an avenue for generating a reservoir of mutants for agricultural production and biomedical research.

DOI:http://dx.doi.org/10.7554/eLife.26248.001

Joint QTL mapping and gene expression analysis identify positional candidate genes influencing pork quality traits

Meat quality traits have an increasing importance in the pig industry because of their strong impact on consumer acceptance. Herewith, we have combined phenotypic and microarray expression data to map loci with potential effects on five meat quality traits recorded in the longissimus dorsi (LD) and gluteus medius (GM) muscles of 350 Duroc pigs, i.e. pH at 24 hours post-mortem (pH₂₄), electric conductivity (CE) and muscle redness (a*), lightness (L*) and yellowness (b*). We have found significant genome-wide associations for CE of LD on SSC4 (~104 Mb), SSC5 (~15 Mb) and SSC13 (~137 Mb), while several additional regions were significantly associated with meat quality traits at the chromosome-wide level. There was a low positional concordance between the associations found for LD and GM traits, a feature that reflects the existence of differences in the genetic determinism of meat quality phenotypes in these two muscles. The performance of an eQTL search for SNPs mapping to the regions associated with meat quality traits demonstrated that the GM a* SSC3 and pH₂₄ SSC17 QTL display positional concordance with cis-eQTL regulating the expression of several genes with a potential role on muscle metabolism.

Genomic structure, expression and association study of the porcine FSD2

The fibronectin type III and SPRY domain containing 2 (FSD2) on porcine chromosome 7 is considered a candidate gene for pork quality, since its two domains, which were present in fibronectin and ryanodine receptor. The fibronectin type III and SPRY domains were first identified in fibronectin and ryanodine receptor, respectively, which are candidate genes for meat quality. The aim of this study was to elucidate the genomic structure of FSD2 and functions of single nucleotide polymorphisms (SNPs) within FSD2 that are related to meat quality in pigs. Using a bacterial artificial chromosome clone sequence, we revealed that porcine FSD2 consisted of 13 exons encoding 750 amino acids. In addition, FSD2 was expressed in heart, longissimus dorsi muscle, psoas muscle, and tendon among 23 kinds of porcine tissues tested. A total of ten SNPs, including four missense mutations, were identified in the exonic region of FSD2, and two major haplotypes were obtained based on the SNP genotypes of 633 Berkshire pigs. Both haplotypes were associated significantly with intramuscular fat content (IMF, P < 0.020) and moisture percentage (MP, P < 0.002). Moreover, haplotype 2 was associated with meat color, affecting yellowness (P = 0.002). These haplotype effects were further supported by the alteration of putative protein structures with amino acid substitutions. Taken together, our results suggest that FSD2 haplotypes are involved in regulating meat quality including IMF, MP, and meat color in pigs, and may be used as meaningful molecular makers to identify pigs with preferable pork quality.

Integration of flow studies for robust selection of mechanoresponsive genes

Blood flow is an essential contributor to plaque growth, composition and initiation. It is sensed by endothelial cells, which react to blood flow by expressing > 1000 genes. The sheer number of genes implies that one needs genomic techniques to unravel their response in disease. Individual genomic studies have been performed but lack sufficient power to identify subtle changes in gene expression. In this study, we investigated whether a systematic meta-analysis of available microarray studies can improve their consistency. We identified 17 studies using microarrays, of which six were performed in vivo and 11 in vitro. The in vivo studies were disregarded due to the lack of the shear profile. Of the in vitro studies, a cross-platform integration of human studies (HUVECs in flow cells) showed high concordance (> 90 %). The human data set identified > 1600 genes to be shear responsive, more than any other study and in this gene set all known mechanosensitive genes and pathways were present. A detailed network analysis indicated a power distribution (e. g. the presence of hubs), without a hierarchical organisation. The average cluster coefficient was high and further analysis indicated an aggregation of 3 and 4 element motifs, indicating a high prevalence of feedback and feed forward loops, similar to prokaryotic cells. In conclusion, this initial study presented a novel method to integrate human-based mechanosensitive studies to increase its power. The robust network was large, contained all known mechanosensitive pathways and its structure revealed hubs, and a large aggregate of feedback and feed forward loops.

Differential Gene Expression across Breed and Sex in Commercial Pigs Administered Fenbendazole and Flunixin Meglumine

Characterizing the variability in transcript levels across breeds and sex in swine for genes that play a role in drug metabolism may shed light on breed and sex differences in drug metabolism. The objective of the study is to determine if there is heterogeneity between swine breeds and sex in transcript levels for genes previously shown to play a role in drug metabolism for animals administered flunixin meglumine or fenbendazole. Crossbred nursery female and castrated male pigs (n = 169) spread across 5 groups were utilized. Sires (n = 15) of the pigs were purebred Duroc, Landrace, Yorkshire or Hampshire boars mated to a common sow population. Animals were randomly placed into the following treatments: no drug (control), flunixin meglumine, or fenbendazole. One hour after the second dosing, animals were sacrificed and liver samples collected. Quantitative Real-Time PCR was used to measure liver gene expression of the following genes: SULT1A1, ABCB1, CYP1A2, CYP2E1, CYP3A22 and CYP3A29. The control animals were used to investigate baseline transcript level differences across breed and sex. Post drug administration transcript differences across breed and sex were investigated by comparing animals administered the drug to the controls. Contrasts to determine fold change were constructed from a model that included fixed and random effects within each drug. Significant (P-value <0.007) basal transcript differences were found across breeds for SULT1A1, CYP3A29 and CYP3A22. Across drugs, significant (P-value <0.0038) transcript differences existed between animals given a drug and controls across breeds and sex for ABCB1, PS and CYP1A2. Significant (P <0.0038) transcript differences across breeds were found for CYP2E1 and SULT1A1 for flunixin meglumine and fenbendazole, respectively. The current analysis found transcript level differences across swine breeds and sex for multiple genes, which provides greater insight into the relationship between flunixin meglumine and fenbendazole and known drug metabolizing genes.

Copy number variation-based genome wide association study reveals additional variants contributing to meat quality in Swine

Pork quality is important both to the meat processing industry and consumers' purchasing attitude. Copy number variation (CNV) is a burgeoning kind of variants that may influence meat quality. In this study, a genome-wide association study (GWAS) was performed between CNVs and meat quality traits in swine. After false discovery rate (FDR) correction, a total of 8 CNVs on 6 chromosomes were identified to be significantly associated with at least one meat quality trait. All of the 8 CNVs were verified by next generation sequencing and six of them were verified by qPCR. Only the haplotype block containing CNV12 is adjacent to significant SNPs associated with meat quality, suggesting the effects of those CNVs were not likely captured by tag SNPs. The DNA dosage and EST expression of CNV12, which overlap with an obesity related gene Netrin-1 (Ntn1), were consistent with Ntn1 RNA expression, suggesting the CNV12 might be involved in the expression regulation of Ntn1 and finally influence meat quality. We concluded that CNVs may contribute to the genetic variations of meat quality beyond SNPs, and several candidate CNVs were worth further exploration.

Population structure and history of the Welsh sheep breeds determined by whole genome genotyping

Background

One of the most economically important areas within the Welsh agricultural sector is sheep farming, contributing around £230 million to the UK economy annually. Phenotypic selection over several centuries has generated a number of native sheep breeds, which are presumably adapted to the diverse and challenging landscape of Wales. Little is known about the history, genetic diversity and relationships of these breeds with other European breeds. We genotyped 353 individuals from 18 native Welsh sheep breeds using the Illumina OvineSNP50 array and characterised the genetic structure of these breeds. Our genotyping data were then combined with, and compared to, those from a set of 74 worldwide breeds, previously collected during the International Sheep Genome Consortium HapMap project.

Results

Model based clustering of the Welsh and European breeds indicated shared ancestry. This finding was supported by multidimensional scaling analysis (MDS), which revealed separation of the European, African and Asian breeds. As expected, the commercial Texel and Merino breeds appeared to have extensive co-ancestry with most European breeds. Consistently high levels of haplotype sharing were observed between native Welsh and other European breeds. The Welsh breeds did not, however, form a genetically homogeneous group, with pairwise F_ST between breeds averaging 0.107 and ranging between 0.020 and 0.201. Four subpopulations were identified within the 18 native breeds, with high homogeneity observed amongst the majority of mountain breeds. Recent effective population sizes estimated from linkage disequilibrium ranged from 88 to 825.

Conclusions

Welsh breeds are highly diverse with low to moderate effective population sizes and form at least four distinct genetic groups. Our data suggest common ancestry between the native Welsh and European breeds. These findings provide the basis for future genome-wide association studies and a first step towards developing genomics assisted breeding strategies in the UK.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-015-0216-x) contains supplementary material, which is available to authorized users.

Polymorphisms and expression analysis of SOX-6 in relation to porcine growth, carcass, and meat quality traits

The aim of the study was to investigate single nucleotide polymorphisms (SNPs) and expression of SOX-6 to support its candidacy for growth, carcass, and meat quality traits in pigs. The first SNP, rs81358375, was associated with pH 45 min post mortem in loin (pH1L), the thickness of backfat and side fat, and carcass length in Pietrain (Pi) population, and related with backfat thickness and daily gain in Duroc × Pietrain F2 (DuPi) population. The other SNP, rs321666676, was associated with meat colour in Pi population. In DuPi population, the protein, not mRNA, level of SOX-6 in high pH1L pigs was significantly less abundant compared with low pH1L pigs, where microRNAs targeting SOX-6 were also differently regulated. This paper shows that SOX-6 could be a potential candidate gene for porcine growth, carcass, and meat quality traits based on genetic association and gene expression.

Single-Nucleotide Polymorphisms in Pig EPHX1 Gene are Associated with Pork Quality Traits

Epoxide hydrolase 1 (EPHX1) plays an important role in both the activation and detoxification of exogenous chemicals. Reverse transcription-polymerase chain reaction (RT-PCR) analysis showed that the highest level of EPHX1 expression occurred in Berkshire liver, which is an organ that plays a key role in detoxification. We examined EPHX1 SNPs to analyze effect on increased expression of EPHX1 gene in Berkshire liver by total of 192 pigs of a pure Berkshire line (males = 97; females = 95). As a result, two nonsynonymous SNPs (nsSNPs) of EPHX1 were found from c.685T>G and c.776C > T, and located in 5th and 6th exons, respectively, which constitute the A/b hydrolase 1 domain of epoxide hydrolase. The nsSNP c.685T > G was significant differences in meat color, protein content, collagen content, and pH24 hr. Especially, T and G alleles of the nsSNP c.685T > G were significantly associated with CIE a*/CIE b* and protein content/pH24 hr, respectively. The nsSNP c.776C > T was significant differences in drip loss and protein content. Among meat quality traits to associate with SNPs, the protein content was only significantly associated with sex. Therefore, it is suggested that nsSNP c.685T > G in EPHX1 gene is a potential to apply as appropriate DNA markers for improvement of porcine economic traits.

Copy number variations in high and low fertility breeding boars

Background

In this study we applied the extreme groups/selective genotyping approach for identifying copy number variations in high and low fertility breeding boars. The fertility indicator was the calculated Direct Boar Effect on litter size (DBE) that was obtained as a by-product of the national genetic evaluation for litter size (BLUP). The two groups of animals had DBE values at the upper (high fertility) and lower (low fertility) end of the distribution from a population of more than 38,000 boars. Animals from these two diverse phenotypes were genotyped with the Porcine SNP60K chip and compared by several approaches in order to prove the feasibility of our CNV analysis and to identify putative markers of fertility.

Results

We have identified 35 CNVRs covering 36.5 Mb or ~1.3% of the porcine genome. Among these 35 CNVRs, 14 were specific to the high fertility group, while 19 CNVRs were specific to the low fertility group which overlap with 137 QTLs of various reproductive traits. The identified 35 CNVRs encompassed 50 genes, among them 40 were specific to the low fertility group, seven to the high fertility group, while three were found in regions that were present in both groups but with opposite gain/loss status. A functional analysis of several databases revealed that the genes found in CNVRs from the low fertility group have been significantly enriched in members of the innate immune system, Toll-like receptor and RIG-I-like receptor signaling and fatty acid oxidation pathways.

Conclusions

We have demonstrated that our analysis pipeline could identify putative CNV markers of fertility, especially in case of low fertility boars.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1473-9) contains supplementary material, which is available to authorized users.

From association to prediction: statistical methods for the dissection and selection of complex traits in plants

Quantification of genotype-to-phenotype associations is central to many scientific investigations, yet the ability to obtain consistent results may be thwarted without appropriate statistical analyses. Models for association can consider confounding effects in the materials and complex genetic interactions. Selecting optimal models enables accurate evaluation of associations between marker loci and numerous phenotypes including gene expression. Significant improvements in QTL discovery via association mapping and acceleration of breeding cycles through genomic selection are two successful applications of models using genome-wide markers. Given recent advances in genotyping and phenotyping technologies, further refinement of these approaches is needed to model genetic architecture more accurately and run analyses in a computationally efficient manner, all while accounting for false positives and maximizing statistical power.

A retrospective analysis of allele frequency changes of major genes during 20 years of selection in the Italian Large White pig breed

In this study, we investigated whether a selection programme based on boar genetic evaluation obtained with a classical BLUP animal model can change allele frequencies in a pig population. All Italian Large White boars born from 1992 to 2012 with estimated breeding value reliability >0.85 (n = 200) were selected among all boars of this breed. Boars were genotyped with markers in major genes (IGF2 intron3-g.3072G>A, MC4R p.D298N, VRTN PRE1 insertion, PRKAG3 p.I199V and FTO g.276T>G). Genotyping data were analysed grouping boars in eight classes according to their year of birth. To evaluate the influence of time on allele frequencies of the genotyped markers, multinomial logistic regression models were computed. Four of five polymorphic sites (IGF2, MC4R, VRTN and FTO) showed significant (p < 0.01) changes in allele frequencies over time due to a progressive and continuous increase of one allele (associated with higher lean meat content, higher average daily gain and favourable feed: gain ratio) and, consequently, decrease of the other one, following the directional selection of the selection programme of this pig breed. The retrospective analysis that was carried out in Italian Large White boars suggests that selection based on methodologies assuming the infinitesimal model is able to modify in a quite short period of time allele frequencies in major genes, increasing the frequency of alleles explaining a relevant (non-infinitesimal) fraction of the overall genetic variability for production traits.

A hidden Markov approach for ascertaining cSNP genotypes from RNA sequence data in the presence of allelic imbalance by exploiting linkage disequilibrium

BACKGROUND

Allelic specific expression (ASE) increases our understanding of the genetic control of gene expression and its links to phenotypic variation. ASE testing is implemented through binomial or beta-binomial tests of sequence read counts of alternative alleles at a cSNP of interest in heterozygous individuals. This requires prior ascertainment of the cSNP genotypes for all individuals. To meet the needs, we propose hidden Markov methods to call SNPs from next generation RNA sequence data when ASE possibly exists.

RESULTS

We propose two hidden Markov models (HMMs), HMM-ASE and HMM-NASE that consider or do not consider ASE, respectively, in order to improve genotyping accuracy. Both HMMs have the advantages of calling the genotypes of several SNPs simultaneously and allow mapping error which, respectively, utilize the dependence among SNPs and correct the bias due to mapping error. In addition, HMM-ASE exploits ASE information to further improve genotype accuracy when the ASE is likely to be present. Simulation results indicate that the HMMs proposed demonstrate a very good prediction accuracy in terms of controlling both the false discovery rate (FDR) and the false negative rate (FNR). When ASE is present, the HMM-ASE had a lower FNR than HMM-NASE, while both can control the false discovery rate (FDR) at a similar level. By exploiting linkage disequilibrium (LD), a real data application demonstrate that the proposed methods have better sensitivity and similar FDR in calling heterozygous SNPs than the VarScan method. Sensitivity and FDR are similar to that of the BCFtools and Beagle methods. The resulting genotypes show good properties for the estimation of the genetic parameters and ASE ratios.

CONCLUSIONS

We introduce HMMs, which are able to exploit LD and account for the ASE and mapping errors, to simultaneously call SNPs from the next generation RNA sequence data. The method introduced can reliably call for cSNP genotypes even in the presence of ASE and under low sequencing coverage. As a byproduct, the proposed method is able to provide predictions of ASE ratios for the heterozygous genotypes, which can then be used for ASE testing.

Feed intake, average daily gain, feed efficiency, and real-time ultrasound traits in Duroc pigs: II. Genomewide association

Efficient use of feed resources has become a clear challenge for the U.S. pork industry as feed costs continue to be the largest variable expense. The availability of the Illumina Porcine60K BeadChip has greatly facilitated whole-genome association studies to identify chromosomal regions harboring genes influencing those traits. The current study aimed at identifying genomic regions associated with variation in feed efficiency and several production traits in a Duroc terminal sire population, including ADFI, ADG, feed conversion ratio, residual feed intake (RFI), real-time ultrasound back fat thickness (BF), ultrasound muscle depth, intramuscular fat content (IMF), birth weight (BW at birth), and weaning weight (BW at weaning). Single trait association analyses were performed using Bayes B models with 35,140 SNP on 18 autosomes after quality control. Significance of nonoverlapping 1-Mb length windows (n = 2,380) were tested across 3 QTL inference methods: posterior distribution of windows variances from Monte Carlo Markov Chain, naive Bayes factor, and nonparametric bootstrapping. Genes within the informative QTL regions for the traits were annotated. A region ranging from166 to 140 Mb (4-Mb length) on SSC 1, approximately 8 Mb upstream of the MC4R gene, was significantly associated with ADFI, ADG, and BF, where SOCS6 and DOK6 are proposed as the most likely candidate genes. Another region affecting BW at weaning was identified on SSC 4 (84-85 Mb), harboring genes previously found to influence both human and cattle height: PLAG1, CHCHD7, RDHE2 (or SDR16C5), MOS, RPS20, LYN, and PENK. No QTL were identified for RFI, IMF, and BW at birth. In conclusion, we have identified several genomic regions associated with traits affecting nutrient utilization that could be considered for future genomic prediction to improve feed utilization.

Genome-wide QTL mapping of nine body composition and bone mineral density traits in pigs

Background

Since the pig is one of the most important livestock animals worldwide, mapping loci that are associated with economically important traits and/or traits that influence animal welfare is extremely relevant for efficient future pig breeding. Therefore, the purpose of this study was a genome-wide mapping of quantitative trait loci (QTL) associated with nine body composition and bone mineral traits: absolute (Fat, Lean) and percentage (FatPC, LeanPC) fat and lean mass, live weight (Weight), soft tissue X-ray attenuation coefficient (R), absolute (BMC) and percentage (BMCPC) bone mineral content and bone mineral density (BMD).

Methods

Data on the nine traits investigated were obtained by Dual-energy X-ray absorptiometry for 551 pigs that were between 160 and 200 days old. In addition, all pigs were genotyped using Illumina’s PorcineSNP60 Genotyping BeadChip. Based on these data, a genome-wide combined linkage and linkage disequilibrium analysis was conducted. Thus, we used 44 611 sliding windows that each consisted of 20 adjacent single nucleotide polymorphisms (SNPs). For the middle of each sliding window a variance component analysis was carried out using ASReml. The underlying mixed linear model included random QTL and polygenic effects, with fixed effects of sex, housing, season and age.

Results

Using a Bonferroni-corrected genome-wide significance threshold of P < 0.001, significant peaks were identified for all traits except BMCPC. Overall, we identified 72 QTL on 16 chromosomes, of which 24 were significantly associated with one trait only and the remaining with more than one trait. For example, a QTL on chromosome 2 included the highest peak across the genome for four traits (Fat, FatPC, LeanPC and R). The nearby gene, ZNF608, is known to be associated with body mass index in humans and involved in starvation in Drosophila, which makes it an extremely good candidate gene for this QTL.

Conclusions

Our QTL mapping approach identified 72 QTL, some of which confirmed results of previous studies in pigs. However, we also detected significant associations that have not been published before and were able to identify a number of new and promising candidate genes, such as ZNF608.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-014-0068-2) contains supplementary material, which is available to authorized users.

Identification of species-specific novel transcripts in pig reproductive tissues using RNA-seq

Although structural properties of the porcine reproductive system are shared by many placental mammals, some combination of these properties is unique to pigs. To explore whether genomic elements specific to pigs could potentially underlie this uniqueness, we made the first step to identify novel transcripts in two representative pig reproductive tissues by the technique of massively parallel sequencing. To automate the whole process, we built a computational pipeline, which can also be easily extended for similar studies in other species. In total, 5516 and 9061 novel transcripts were found, and 159 and 252 novel transcripts appear to be specific to pigs for the placenta and testis respectively. Furthermore, these novel transcripts were found to be enriched in quantitative trait loci (QTL) regions for reproduction traits in pigs. We validated eight of these novel transcripts by quantitative real-time PCR. With respect to their genomic organization and their functional relationship to reproduction, these transcripts need to be further validated and explored in various pig breeds to better comprehend the relevant aspects of pig physiology that contribute to reproductive performance.