Genome wide analysis reveals single nucleotide polymorphisms associated with fatness and putative novel copy number variants in three pig breeds

Identifying Significant SNPs of the Total Number of Piglets Born and Their Relationship with Leg Bumps in Pigs

The aim of this study was to identify genetic variants and pathways associated with the total number of piglets born and to investigate the potential negative consequences of the intensive selection for reproductive traits, particularly the formation of bumps on the legs of pigs. We used genome-wide association analysis and methods for identifying selection signatures. As a result, 47 SNPs were identified, localized in genes that play a significant role during sow pregnancy. These genes are involved in follicle growth and development (SGC), early embryonic development (CCDC3, LRRC8C, LRFN3, TNFRSF19), endometrial receptivity and implantation (NEBL), placentation, and embryonic development (ESRRG, GHRHR, TUSC3, NBAS). Several genes are associated with disorders of the nervous system and brain development (BCL11B, CDNF, ULK4, CC2D2A, KCNK2). Additionally, six SNPs are associated with the formation of bumps on the legs of pigs. These variants include intronic variants in the CCDC3, ULK4, and MINDY4 genes, as well as intergenic variants, regulatory region variants, and variants in the exons of non-coding transcripts. The results suggest important biological pathways and genetic variants associated with sow fertility and highlight the potential negative impacts on the health and physical condition of pigs.

Identification of a novel QTL for lean meat percentage using imputed genotypes

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.

Additive and Dominance Genome-Wide Association Studies Reveal the Genetic Basis of Heterosis Related to Growth Traits of Duhua Hybrid Pigs

Heterosis has been extensively used for pig genetic breeding and production, but the genetic basis of heterosis remains largely elusive. Crossbreeding between commercial and native breeds provides a good model to parse the genetic basis of heterosis. This study uses Duhua hybrid pigs, a crossbreed of Duroc and Liangguang small spotted pigs, as materials to explore the genetic basis underlying heterosis related to growth traits at the genomic level. The mid-parent heterosis (MPH) analysis showed heterosis of this Duhua offspring on growth traits. In this study, we examined the impact of additive and dominance effects on 100 AGE (age adjusted to 100 kg) and 100 BF (backfat thickness adjusted to 100 kg) of Duhua hybrid pigs. Meanwhile, we successfully identified SNPs associated with growth traits through both additive and dominance GWASs (genome-wide association studies). These findings will facilitate the subsequent in-depth studies of heterosis in the growth traits of Duhua pigs.

Risk of Sperm Disorders and Impaired Fertility in Frozen-Thawed Bull Semen: A Genome-Wide Association Study

Cryopreservation is a widely used method of semen conservation in animal breeding programs. This process, however, can have a detrimental effect on sperm quality, especially in terms of its morphology. The resultant sperm disorders raise the risk of reduced sperm fertilizing ability, which poses a serious threat to the long-term efficacy of livestock reproduction and breeding. Understanding the genetic factors underlying these effects is critical for maintaining sperm quality during cryopreservation, and for animal fertility in general. In this regard, we performed a genome-wide association study to identify genomic regions associated with various cryopreservation sperm abnormalities in Holstein cattle, using single nucleotide polymorphism (SNP) markers via a high-density genotyping assay. Our analysis revealed a significant association of specific SNPs and candidate genes with absence of acrosomes, damaged cell necks and tails, as well as wrinkled acrosomes and decreased motility of cryopreserved sperm. As a result, we identified candidate genes such as POU6F2, LPCAT4, DPYD, SLC39A12 and CACNB2, as well as microRNAs (bta-mir-137 and bta-mir-2420) that may play a critical role in sperm morphology and disorders. These findings provide crucial information on the molecular mechanisms underlying acrosome integrity, motility, head abnormalities and damaged cell necks and tails of sperm after cryopreservation. Further studies with larger sample sizes, genome-wide coverage and functional validation are needed to explore causal variants in more detail, thereby elucidating the mechanisms mediating these effects. Overall, our results contribute to the understanding of genetic architecture in cryopreserved semen quality and disorders in bulls, laying the foundation for improved animal reproduction and breeding.

Finding Predictors of Leg Defects in Pigs Using CNV-GWAS

One of the most important areas of modern genome research is the search for meaningful relationships between genetic variants and phenotypes. In the livestock field, there has been research demonstrating the influence of copy number variants (CNVs) on phenotypic variation. Despite the wide range in the number and size of detected CNVs, a significant proportion differ between breeds and their functional effects are underestimated in the pig industry. In this work, we focused on the problem of leg defects in pigs (lumps/growths in the area of the hock joint on the hind legs) and focused on searching for molecular genetic predictors associated with this trait for the selection of breeding stock. The study was conducted on Large White pigs using three CNV calling tools (PennCNV, QuantiSNP and R-GADA) and the CNVRanger association analysis tool (CNV-GWAS). As a result, the analysis identified three candidate CNVRs associated with the formation of limb defects. Subsequent functional analysis suggested that all identified CNVs may act as potential predictors of the hock joint phenotype of pigs. It should be noted that the results obtained indicate that all significant regions are localized in genes (CTH, SRSF11, MAN1A1 and LPIN1) responsible for the metabolism of amino acids, fatty acids, glycerolipids and glycerophospholipids, thereby related to the immune response, liver functions, content intramuscular fat and animal fatness. These results are consistent with previously published studies, according to which a predisposition to the formation of leg defects can be realized through genetic variants associated with the functions of the liver, kidneys and hematological characteristics.

A Comprehensive Genomic Analysis of Chinese Indigenous Ningxiang Pigs: Genomic Breed Compositions, Runs of Homozygosity, and Beyond

Ningxiang pigs are a renowned indigenous pig breed in China, known for their meat quality, disease resistance, and environmental adaptability. In recent decades, consumer demand for meats from indigenous breeds has grown significantly, fueling the selection and crossbreeding of Ningxiang pigs (NXP). The latter has raised concerns about the conservation and sustainable use of Ningxiang pigs as an important genetic resource. To address these concerns, we conducted a comprehensive genomic study using 2242 geographically identified Ningxiang pigs. The estimated genomic breed composition (GBC) suggested 2077 pigs as purebred Ningxiang pigs based on a ≥94% NXP-GBC cut-off. The remaining 165 pigs were claimed to be crosses, including those between Duroc and Ningxiang pigs and between Ningxiang and Shaziling pigs, and non-Ningxiang pigs. Runs of homozygosity (ROH) were identified in the 2077 purebred Ningxiang pigs. The number and length of ROH varied between individuals, with an average of 32.14 ROH per animal and an average total length of 202.4 Mb per animal. Short ROH (1-5 Mb) was the most abundant, representing 66.5% of all ROH and 32.6% of total ROH coverage. The genomic inbreeding estimate was low (0.089) in purebred Ningxiang pigs compared to imported western pig breeds. Nine ROH islands were identified, pinpointing candidate genes and QTLs associated with economic traits of interest, such as reproduction, carcass and growth traits, lipid metabolism, and fat deposition. Further investigation of these ROH islands and candidate genes is anticipated to better understand the genomics of Ningxiang pigs.

Detection of Selection Signatures in Anqing Six-End-White Pigs Based on Resequencing Data

As a distinguished Chinese indigenous pig breed that exhibits disease resistance and high meat quality, the Anqing six-end-white (AQ) pig represents a valuable germplasm resource for improving the quality of the pig breeding industry. In this study, 24 AQ pigs that were distantly blood-related and 6 Asian Wild Boar (AWB) were selected for 10× deep-genome resequencing. The signatures of the selection were analyzed to explore the genetic basis of their germplasm characteristics and to identify excellent germplasm-related functional genes based on NGS data. A total of 49,289,052 SNPs and 6,186,123 indels were detected across the genome in 30 pigs. Most of the genetic variations were synonym mutations and existed in the intergenic region. We identified 275 selected regions (top 1%) harboring 85 genes by applying a crossover approach based on genetic differentiation (F_ST) and polymorphism levels (π ratio). Some genes were found to be positively selected in AQ pigs' breeding. The SMPD4 and DDX18 genes were involved in the immune response to pseudorabies virus (PRV) and porcine reproductive and respiratory syndrome virus (PRRSV). The BCL6 and P2RX6 genes were involved in biological regulation of immune T cells and phagocytes. The SLC7A4 and SPACA4 genes were related to reproductive performance. The MSTN and HIF1A genes were related to fat deposition and muscle development. Moreover, 138 overlapping regions were detected in selected regions and ROH islands of AQ pigs. Additionally, we found that the QTLs with the most overlapping regions were related to back fat thickness, meat color, pH value, fatty acid content, immune cells, parasitic immunity, and bacterial immunity. Based on functional enrichment analysis and QTLs mapping, we conducted further research on the molecular genetic basis of germplasm traits (disease resistance and excellent meat quality). These results are a reliable resource for conserving germplasm resources and exploiting molecular markers of AQ pigs.

Lost pigs of Angola: Whole genome sequencing reveals unique regions of selection with emphasis on metabolism and feed efficiency

Angola, in the western coast of Africa, has been through dramatic social events that have led to the near-disappearance of native swine populations, and the recent introduction of European exotic breeds has also contributed to the erosion of this native swine repertoire. In an effort to investigate the genetic basis of native pigs in Angola (ANG) we have generated whole genomes from animals of a remote local pig population in Huambo province, which we have compared with 78 genomes of European and Asian pig breeds as well as European and Asian wild boars that are currently in public domain. Analyses of population structure showed that ANG pigs grouped within the European cluster and were clearly separated from Asian pig breeds. Pairwise F_ST ranged from 0.14 to 0.26, ANG pigs display lower levels of genetic differentiation towards European breeds. Finally, we have identified candidate regions for selection using a complementary approach based on various methods. All results suggest that selection towards feed efficiency and metabolism has occurred. Moreover, all analysis identified CDKAL1 gene, which is related with insulin and cholesterol metabolism, as a candidate gene overlapping signatures of selection unique to ANG pigs. This study presents the first assessment of the genetic relationship between ANG pigs and other world breeds and uncovers selection signatures that may indicate adaptation features unique to this important genetic resource.

Genome-Wide Association Study Reveals Additive and Non-Additive Effects on Growth Traits in Duroc Pigs

Growth rate plays a critical role in the pig industry and is related to quantitative traits controlled by many genes. Here, we aimed to identify causative mutations and candidate genes responsible for pig growth traits. In this study, 2360 Duroc pigs were used to detect significant additive, dominance, and epistatic effects associated with growth traits. As a result, a total number of 32 significant SNPs for additive or dominance effects were found to be associated with various factors, including adjusted age at a specified weight (AGE), average daily gain (ADG), backfat thickness (BF), and loin muscle depth (LMD). In addition, the detected additive significant SNPs explained 2.49%, 3.02%, 3.18%, and 1.96% of the deregressed estimated breeding value (DEBV) variance for AGE, ADG, BF, and LMD, respectively, while significant dominance SNPs could explain 2.24%, 13.26%, and 4.08% of AGE, BF, and LMD, respectively. Meanwhile, a total of 805 significant epistatic effects SNPs were associated with one of ADG, AGE, and LMD, from which 11 sub-networks were constructed. In total, 46 potential genes involved in muscle development, fat deposition, and regulation of cell growth were considered as candidates for growth traits, including CD55 and NRIP1 for AGE and ADG, TRIP11 and MIS2 for BF, and VRTN and ZEB2 for LMD, respectively. Generally, in this study, we detected both new and reported variants and potential candidate genes for growth traits of Duroc pigs, which might to be taken into account in future molecular breeding programs to improve the growth performance of pigs.

Copy Number Variations Contribute to Intramuscular Fat Content Differences by Affecting the Expression of PELP1 Alternative Splices in Pigs

Simple Summary

Copy number variation (CNV) is a type of variant that may influence meat quality of, for example intramuscular fat (IMF). In this study, a genome-wide association study (GWAS) was then performed between CNVs and IMF in a pig F2 resource population. A total of 19 CNVRs were found to be significantly associated with IMF. RNA-seq and qPCR validation results indicated that CNV150, which is located on the 3′UTR end of the proline, as well as glutamate and the leucine rich protein 1 (PELP1) gene may affect the expression of PELP1 alternative splices. We infer that the CNVR may influence IMF content by regulating the alternative splicing of the PELP1 gene and ultimately affects the structure of the PELP1 protein. These findings suggest a novel mechanistic approach for meat quality improvement in animals and the potential treatment of insulin resistance in human beings.

Abstract

Intramuscular fat (IMF) is a key meat quality trait. Research on the genetic mechanisms of IMF decomposition is valuable for both pork quality improvement and the treatment of obesity and type 2 diabetes. Copy number variations (CNVs) are a type of variant that may influence meat quality. In this study, a total of 1185 CNV regions (CNVRs) including 393 duplicated CNVRs, 432 deleted CNVRs, and 361 CNVRs with both duplicated and deleted status were identified in a pig F2 resource population using next-generation sequencing data. A genome-wide association study (GWAS) was then performed between CNVs and IMF, and a total of 19 CNVRs were found to be significantly associated with IMF. QTL colocation analysis indicated that 3 of the 19 CNVRs overlapped with known QTLs. RNA-seq and qPCR validation results indicated that CNV150, which is located on the 3′UTR end of the proline, as well as glutamate and the leucine rich protein 1 (PELP1) gene may affect the expression of PELP1 alternative splices. Sequence alignment and Alphafold2 structure prediction results indicated that the two alternative splices of PELP1 have a 23 AA sequence variation and a helix-fold structure variation. This region is located in the region of interaction between PELP1 and other proteins which have been reported to be significantly associated with fat deposition or insulin resistance. We infer that the CNVR may influence IMF content by regulating the alternative splicing of the PELP1 gene and ultimately affects the structure of the PELP1 protein. In conclusion, we found some CNVRs, especially CNV150, located in PELP1 that affect IMF. These findings suggest a novel mechanistic approach for meat quality improvement in animals and the potential treatment of insulin resistance in human beings.

Fine Mapping of a Major Backfat QTL Reveals a Causal Regulatory Variant Affecting the CCND2 Gene

Backfat is an important trait in pork production, and it has been included in the breeding objectives of genetic companies for decades. Although adipose tissue is a good energy storage, excessive fat results in reduced efficiency and economical losses. A large QTL for backfat thickness on chromosome 5 is still segregating in different commercial pig breeds. We fine mapped this QTL region using a genome-wide association analysis (GWAS) with 133,358 genotyped animals from five commercial populations (Landrace, Pietrain, Large White, Synthetic, and Duroc) imputed to the porcine 660K SNP chip. The lead SNP was located at 5:66103958 (G/A) within the third intron of the CCND2 gene, with the G allele associated with more backfat, while the A allele is associated with less backfat. We further phased the QTL region to discover a core haplotype of five SNPs associated with low backfat across three breeds. Linkage disequilibrium analysis using whole-genome sequence data revealed three candidate causal variants within intronic regions and downstream of the CCND2 gene, including the lead SNP. We evaluated the association of the lead SNP with the expression of the genes in the QTL region (including CCND2) in a large cohort of 100 crossbred samples, sequenced in four different tissues (lung, spleen, liver, muscle). Results show that the A allele increases the expression of CCND2 in an additive way in three out of four tissues. Our findings indicate that the causal variant for this QTL region is a regulatory variant within the third intron of the CCND2 gene affecting the expression of CCND2.

Increased accuracy of genomic predictions for growth under chronic thermal stress in rainbow trout by prioritizing variants from GWAS using imputed sequence data

Through imputation of genotypes, genome-wide association study (GWAS) and genomic prediction (GP) using whole-genome sequencing (WGS) data are cost-efficient and feasible in aquaculture breeding schemes. The objective was to dissect the genetic architecture of growth traits under chronic heat stress in rainbow trout (Oncorhynchus mykiss) and to assess the accuracy of GP based on imputed WGS and different preselected single nucleotide polymorphism (SNP) arrays. A total of 192 and 764 fish challenged to a heat stress experiment for 62 days were genotyped using a customized 1 K and 26 K SNP panels, respectively, and then, genotype imputation was performed from a low-density chip to WGS using 102 parents (36 males and 66 females) as the reference population. Imputed WGS data were used to perform GWAS and test GP accuracy under different preselected SNP scenarios. Heritability was estimated for body weight (BW), body length (BL) and average daily gain (ADG). Estimates using imputed WGS data ranged from 0.33 ± 0.05 to 0.55 ± 0.05 for growth traits under chronic heat stress. GWAS revealed that the top five cumulatively SNPs explained a maximum of 0.94%, 0.86% and 0.51% of genetic variance for BW, BL and ADG, respectively. Some important functional candidate genes associated with growth-related traits were found among the most important SNPs, including signal transducer and activator of transcription 5B and 3 (STAT5B and STAT3, respectively) and cytokine-inducible SH2-containing protein (CISH). WGS data resulted in a slight increase in prediction accuracy compared with pedigree-based method, whereas preselected SNPs based on the top GWAS hits improved prediction accuracies, with values ranging from 1.2 to 13.3%. Our results support the evidence of the polygenic nature of growth traits when measured under heat stress. The accuracies of GP can be improved using preselected variants from GWAS, and the use of WGS marginally increases prediction accuracy.

An association study on imputed whole-genome resequencing from high-throughput sequencing data for body traits in crossbred pigs

Body traits are important economic factors in the pig industry. Genome-wide association studies (GWASs) have been widely applied using high-density genotype data to detect QTL in pigs. The aim of the present study was to detect the genetic variants significantly associated with body traits in crossbred pigs using the Illumina Porcine SNP50 BeadChip and imputed whole-genome sequence data. A set of seven body traits - body length, body height, chest circumference, cannon bone circumference, leg buttock circumference, back fat thickness and loin muscle depth - were measured. Moderate to high heritabilities were obtained for most traits (from 0.14 to 0.46), and significant genetic and phenotypic correlations among them were observed. GWAS identified 714 significantly associated SNPs located at 39 regions on all autosomes for body traits, and a total of seven functionally related candidate genes: PIK3CD, HOXA, PCGF2, CHST11, COL2A1, BMI1 and OSR2. Functional enrichment analysis revealed that candidate genes were enriched in the estrogen signaling pathway, embryonic skeletal system morphogenesis and embryonic skeletal system development. These results aim to uncover the genetic mechanisms underlying body development and marker-assisted selection programs focusing on body traits in pigs.

Mitochondrially localized MPZL3 emerges as a signaling hub of mammalian physiology

MPZL3 is a nuclear-encoded, mitochondrially localized, immunoglobulin-like V-type protein that functions as a key regulator of epithelial cell differentiation, lipid metabolism, ROS production, glycemic control, and energy expenditure. Recently, MPZL3 has surfaced as an important modulator of sebaceous gland function and of hair follicle cycling, an organ transformation process that is also governed by peripheral clock gene activity and PPARγ. Given the phenotype similarities and differences between Mpzl3 and Pparγ knockout mice, we propose that MPZL3 serves as a signaling hub that is regulated by core clock gene products and/or PPARγ to translate signals from these nuclear transcription factors to the mitochondria to modulate circadian and metabolic regulation. Conservation between murine and human MPZL3 suggests that human MPZL3 may have similarly complex functions in health and disease. We summarize current knowledge and discuss future directions to elucidate the full spectrum of MPZL3 functions in mammalian physiology.

Genome-wide association study reveals a quantitative trait locus and two candidate genes on Sus scrofa chromosome 5 affecting intramuscular fat content in Suhuai pigs

Intramuscular fat content (IFC) is an essential quantitative trait of meat, affecting multiple meat quality indicators. A certain amount of IFC could not only improve the sensory score of pork but also increase the flavour, tenderness, juiciness and shelf-life. To dissect the genetic determinants of IFC, two methods, including genome-wide efficient mixed-model analysis (GEMMA) and linkage disequilibrium adjusted kinships (LDAKs), were used to carry out genome-wide association studies for IFC in Suhuai pig population. A total of 14 and 18 significant single nucleotide polymorphisms (SNPs) were identified by GEMMA and LDAK, respectively. The results of these two methods were highly consistent and all 14 significant SNPs in GEMMA were detected by LDAK. Seven of the 18 SNPs reached the genome-wide significance level (P < 9.85E-07) while 11 cases reached the suggestive significance level (P < 1.77E-05). These significant SNPs were mainly distributed on Sus scrofa chromosome (SSC) 5, 3, and 7. Moreover, one locus resides in a 2.27 Mb (71.37-73.64 Mb) region on SSC5 harbouring 13 significant SNPs associated with IFC, and the lead SNP (rs81302978) also locates in this region. Linkage disequilibrium (LD) analysis showed that there were four pairs of complete LD (r² = 1) among these 13 SNPs, and the remaining 9 SNPs with incomplete LD (r² ≠ 1) were selected for subsequent analyses of IFC. Association analyses showed that 7 out of 9 SNPs were significantly associated with IFC (P < 0.05) in 330 Suhuai pigs, and the other 2 SNPs tended to reach a significant association level with IFC (P < 0.1). The phenotypic variance explained (PVE) range of these 9 SNPs was 0.92-3.55%. Meanwhile, the lead SNP was also significantly associated (rs81302978) with IFC (P < 0.05) in 378 commercial hybrid pigs (Pietrain × Duroc) × (Landrace × Yorkshire) (PDLY), and the PVE was 1.38%. Besides, two lipid metabolism-relevant candidate genes, the leucine rich repeat kinase 2 (LRRK2) and PDZ domain containing ring finger 4 (PDZRN4) were identified in the 2.27 Mb region on SSC5. In conclusion, our results may provide a set of markers useful for genetic improvement of IFC in pigs and will advance the genome selection process of IFC on pig breeding programmes.

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.

Genome-Wide Detection of Copy Number Variations and Their Association With Distinct Phenotypes in the World's Sheep

Copy number variations (CNVs) are a major source of structural variation in mammalian genomes. Here, we characterized the genome-wide CNV in 2059 sheep from 67 populations all over the world using the Ovine Infinium HD (600K) SNP BeadChip. We tested their associations with distinct phenotypic traits by conducting multiple independent genome-wide tests. In total, we detected 7547 unique CNVs and 18,152 CNV events in 1217 non-redundant CNV regions (CNVRs), covering 245 Mb (∼10%) of the whole sheep genome. We identified seven CNVRs with frequencies correlating to geographical origins and 107 CNVRs overlapping 53 known quantitative trait loci (QTLs). Gene ontology and pathway enrichment analyses of CNV-overlapping genes revealed their common involvement in energy metabolism, endocrine regulation, nervous system development, cell proliferation, immune, and reproduction. For the phenotypic traits, we detected significantly associated (adjusted P < 0.05) CNVRs harboring functional candidate genes, such as SBNO2 for polycerate; PPP1R11 and GABBR1 for tail weight; AKT1 for supernumerary nipple; CSRP1, WNT7B, HMX1, and FGFR3 for ear size; and NOS3 and FILIP1 in Wadi sheep; SNRPD3, KHDRBS2, and SDCCAG3 in Hu sheep; NOS3, BMP1, and SLC19A1 in Icelandic; CDK2 in Finnsheep; MICA in Romanov; and REEP4 in Texel sheep for litter size. These CNVs and associated genes are important markers for molecular breeding of sheep and other livestock species.

Transcriptomics analysis of differentially expressed genes in subcutaneous and perirenal adipose tissue of sheep as affected by their pre- and early postnatal malnutrition histories

BACKGROUND

Early life malnutrition is known to target adipose tissue with varying impact depending on timing of the insult. This study aimed to identify differentially expressed genes in subcutaneous (SUB) and perirenal (PER) adipose tissue of 2.5-years old sheep to elucidate the biology underlying differential impacts of late gestation versus early postnatal malnutrition on functional development of adipose tissues. Adipose tissues were obtained from 37 adult sheep born as twins to dams fed either NORM (fulfilling energy and protein requirements), LOW (50% of NORM) or HIGH (110% of protein and 150% of energy requirements) diets in the last 6-weeks of gestation. From day 3 to 6 months of age, lambs were fed high-carbohydrate-high-fat (HCHF) or moderate low-fat (CONV) diets, and thereafter the same moderate low-fat diet.

RESULTS

The gene expression profile of SUB in the adult sheep was not affected by the pre- or early postnatal nutrition history. In PER, 993 and 186 differentially expressed genes (DEGs) were identified in LOW versus HIGH and NORM, respectively, but no DEG was found between HIGH and NORM. DEGs identified in the mismatched pre- and postnatal nutrition groups LOW-HCHF (101) and HIGH-HCHF (192) were largely downregulated compared to NORM-CONV. Out of 831 DEGs, 595 and 236 were up- and downregulated in HCHF versus CONV, respectively. The functional enrichment analyses revealed that transmembrane (ion) transport activities, motor activities related to cytoskeletal and spermatozoa function (microtubules and the cytoskeletal motor protein, dynein), and responsiveness to the (micro) environmental extracellular conditions, including endocrine and nervous stimuli were enriched in the DEGs of LOW versus HIGH and NORM. We confirmed that mismatched pre- and postnatal feeding was associated with long-term programming of adipose tissue remodeling and immunity-related pathways. In agreement with phenotypic measurements, early postnatal HCHF feeding targeted pathways involved in kidney cell differentiation, and mismatched LOW-HCHF sheep had specific impairments in cholesterol metabolism pathways.

CONCLUSIONS

Both pre- and postnatal malnutrition differentially programmed (patho-) physiological pathways with implications for adipose functional development associated with metabolic dysfunctions, and PER was a major target.

Whole genome variants across 57 pig breeds enable comprehensive identification of genetic signatures that underlie breed features

Background

A large number of pig breeds are distributed around the world, their features and characteristics vary among breeds, and they are valuable resources. Understanding the underlying genetic mechanisms that explain across-breed variation can help breeders develop improved pig breeds.

Results

In this study, we performed GWAS using a standard mixed linear model with three types of genome variants (SNP, InDel, and CNV) that were identified from public, whole-genome, sequencing data sets. We used 469 pigs of 57 breeds, and we identified and analyzed approximately 19 million SNPs, 1.8 million InDels, and 18,016 CNVs. We defined six biological phenotypes by the characteristics of breed features to identify the associated genome variants and candidate genes, which included coat color, ear shape, gradient zone, body weight, body length, and body height. A total of 37 candidate genes was identified, which included 27 that were reported previously (e.g., PLAG1 for body weight), but the other 10 were newly detected candidate genes (e.g., ADAMTS9 for coat color).

Conclusion

Our study indicated that using GWAS across a modest number of breeds with high density genome variants provided efficient mapping of complex traits.

Supplementary Information

Supplementary information accompanies this paper at 10.1186/s40104-020-00520-8.

Genome-wide association study and pathway analysis for fat deposition traits in nellore cattle raised in pasture-based systems

Genome-wide association study (GWAS) is a powerful tool to identify candidate genes and genomic regions underlying key biological mechanisms associated with economically important traits. In this context, the aim of this study was to identify genomic regions and metabolic pathways associated with backfat thickness (BFT) and rump fat thickness (RFT) in Nellore cattle, raised in pasture-based systems. Ultrasound-based measurements of BFT and RFT (adjusted to 18 months of age) were collected in 11,750 animals, with 39,903 animals in the pedigree file. Additionally, 1,440 animals were genotyped using the GGP-indicus 35K SNP chip, containing 33,623 SNPs after the quality control. The single-step GWAS analyses were performed using the BLUPF90 family programs. Candidate genes were identified through the Ensembl database incorporated in the BioMart tool, while PANTHER and REVIGO were used to identify the key metabolic pathways and gene networks. A total of 18 genomic regions located on 10 different chromosomes and harbouring 23 candidate genes were identified for BFT. For RFT, 22 genomic regions were found on 14 chromosomes, with a total of 29 candidate genes identified. The results of the pathway analyses showed important genes for BFT, including TBL1XR1, AHCYL2, SLC4A7, AADAT, VPS53, IDH2 and ETS1, which are involved in lipid metabolism, synthesis of cellular amino acids, transport of solutes, transport between Golgi Complex membranes, cell differentiation and cellular development. The main genes identified for RFT were GSK3β, LRP1B, EXT1, GRB2, SORCS1 and SLMAP, which are involved in metabolic pathways such as glycogen synthesis, lipid transport and homeostasis, polysaccharide and carbohydrate metabolism. Polymorphisms located in these candidate genes can be incorporated in commercial genotyping platforms to improve the accuracy of imputation and genomic evaluations for carcass fatness. In addition to uncovering biological mechanisms associated with carcass quality, the key gene pathways identified can also be incorporated in biology-driven genomic prediction methods.

Genomics of Heat Tolerance in Reproductive Performance Investigated in Four Independent Maternal Lines of Pigs

Improving swine climatic resilience through genomic selection has the potential to minimize welfare issues and increase the industry profitability. The main objective of this study was to investigate the genetic and genomic determinism of tolerance to heat stress in four independent purebred populations of swine. Three female reproductive traits were investigated: total number of piglets born (TNB), number of piglets born alive (NBA) and average birth weight (ABW). More than 80,000 phenotypic and 12,000 genotyped individuals were included in this study. Genomic random-regression models were fitted regressing the phenotypes of interest on a set of 95 environmental covariates extracted from public weather station records. The models yielded estimates of (genomic) reactions norms for individual pigs, as indicator of heat tolerance. Heat tolerance is a heritable trait, although the heritabilities are larger under comfortable than heat-stress conditions (larger than 0.05 vs. 0.02 for TNB; 0.10 vs. 0.05 for NBA; larger than 0.20 vs. 0.10 for ABW). TNB showed the lowest genetic correlation (-38%) between divergent climatic conditions, being the trait with the strongest impact of genotype by environment interaction, while NBA and ABW showed values slightly negative or equal to zero reporting a milder impact of the genotype by environment interaction. After estimating genetic parameters, a genome-wide association study was performed based on the single-step GBLUP method. Heat tolerance was observed to be a highly polygenic trait. Multiple and non-overlapping genomic regions were identified for each trait based on the genomic breeding values for reproductive performance under comfortable or heat stress conditions. Relevant regions were found on chromosomes (SSC) 1, 3, 5, 6, 9, 11, and 12, although there were important regions across all autosomal chromosomes. The genomic region located on SSC9 appears to be of particular interest since it was identified for two traits (TNB and NBA) and in two independent populations. Heat tolerance based on reproductive performance indicators is a heritable trait and genetic progress for heat tolerance can be achieved through genetic or genomic selection. Various genomic regions and candidate genes with important biological functions were identified, which will be of great value for future functional genomic studies.

Genome-wide detection of copy number variants in European autochthonous and commercial pig breeds by whole-genome sequencing of DNA pools identified breed-characterising copy number states

In this study, we identified copy number variants (CNVs) in 19 European autochthonous pig breeds and in two commercial breeds (Italian Large White and Italian Duroc) that represent important genetic resources for this species. The genome of 725 pigs was sequenced using a breed-specific DNA pooling approach (30-35 animals per pool) obtaining an average depth per pool of 42×. This approach maximised CNV discovery as well as the related copy number states characterising, on average, the analysed breeds. By mining more than 17.5 billion reads, we identified a total of 9592 CNVs (~683 CNVs per breed) and 3710 CNV regions (CNVRs; 1.15% of the reference pig genome), with an average of 77 CNVRs per breed that were considered as private. A few CNVRs were analysed in more detail, together with other information derived from sequencing data. For example, the CNVR encompassing the KIT gene was associated with coat colour phenotypes in the analysed breeds, confirming the role of the multiple copies in determining breed-specific coat colours. The CNVR covering the MSRB3 gene was associated with ear size in most breeds. The CNVRs affecting the ELOVL6 and ZNF622 genes were private features observed in the Lithuanian Indigenous Wattle and in the Turopolje pig breeds respectively. Overall, the genome variability unravelled here can explain part of the genetic diversity among breeds and might contribute to explain their origin, history and adaptation to a variety of production systems.

CNV analysis of Meishan pig by next-generation sequencing and effects of AHR gene CNV on pig reproductive traits

Background

Reproductive performance of livestock is an economically important aspect of global food production. The Chinese Meishan pig is a prolific breed, with an average of three to five more piglets per litter than European breeds; however, the genetic basis for this difference is not well understood.

Results

In this study, we investigated copy number variations (CNVs) of 32 Meishan pigs and 29 Duroc pigs by next-generation sequencing. A genome-wide analysis of 61 pigs revealed 12,668 copy number variable regions (CNVRs) that were further divided into three categories based on copy number (CN) of the whole population, i.e., gain (n = 7,638), and loss (n = 5,030) CNVRs. We then compared Meishan and Duroc pigs and identified 17.17 Mb of 6,387 CNVRs that only existing in Meishan pigs CNVRs that overlapped the reproduction-related gene encoding the aryl hydrocarbon receptor (AHR) gene. We found that normal AHR CN was more frequent than CN loss in four different pig breeds. An association analysis showed that AHR CN had a positive effect on litter size (P < 0.05) and that a higher CN was associated with higher total number born (P < 0.05), number born alive (P < 0.05), number of weaned piglets, and birth weight.

Conclusions

The present study provides comprehensive CNVRs for Meishan and Duroc pigs through large-scale population resequencing. Our results provide a supplement for the high-resolution map of copy number variation in the porcine genome and valuable information for the investigation of genomic structural variation underlying traits of interest in pig. In addition, the association results provide evidence for AHR as a candidate gene associated with reproductive traits that can be used as a genetic marker in pig breeding programs.

Molecular mapping of major QTL conferring resistance to orange wheat blossom midge (Sitodiplosis mosellana) in Chinese wheat varieties with selective populations

Two novel midge resistance QTL were mapped to a 4.9-Mb interval on chromosome arm 4AL based on the genetic maps constructed with SNP markers. Orange wheat blossom midge (OWBM) is a devastating insect pest affecting wheat production. In order to detect OWBM resistance genes and quantitative trait loci (QTL) for wheat breeding, two recombinant inbred line (RIL) populations were established and used for molecular mapping. A total of seven QTL were detected on chromosomes 2D, 4A, 4D and 7D, respectively, of which positive alleles were all from the resistant parents except for the QTL on 7D. Two stable QTL (QSm.hbau-4A.2-1 and QSm.hbau-4A.2-2) were detected in both populations with the LOD scores ranging from 5.58 to 29.22 under all three environments, and they explained a combined phenotypic variation of 24.4-44.8%. These two novel QTL were mapped to a 4.9-Mb physical interval. The single-nucleotide polymorphism (SNP) markers AX-109543456, AX-108942696 and AX-110928325 were closely linked to the QTL and could be used for marker-assisted selection (MAS) for OWBM resistance in wheat breeding programs.

Impact of merging commercial breeding lines on the genetic diversity of Landrace pigs

BACKGROUND

The pig breeding industry has undergone a large number of mergers in the past decades. Various commercial lines were merged or discontinued, which is expected to reduce the genetic diversity of the pig species. The objective of the current study was to investigate the genetic diversity of different former Dutch Landrace breeding lines and quantify their relationship with the current Dutch Landrace breed that originated from these lines.

RESULTS

Principal component analysis clearly divided the former Landrace lines into two main clusters, which are represented by Norwegian/Finnish Landrace lines and Dutch Landrace lines. Structure analysis revealed that each of the lines that are present in the Dutch Gene bank has a unique genetic identity. The current Dutch Landrace breed shows a high level of admixture and is closely related to the six former lines. The Dumeco N-line, which is conserved in the Dutch Gene bank, is poorly represented in the current Dutch Landrace. All seven lines (the six former and the current line) contribute almost equally to the genetic diversity of the Dutch Landrace breed. As expected, the current Dutch Landrace breed comprises only a small proportion of unique genetic diversity that was not present in the other lines. The genetic diversity level, as measured by Eding's core set method, was equal to 0.89 for the current Dutch Landrace breed, whereas total genetic diversity across the seven lines, measured by the same method, was equal to 0.99.

CONCLUSIONS

The current Dutch Landrace breed shows a high level of admixture and is closely related to the six former Dutch Landrace lines. Merging of commercial Landrace lines has reduced the genetic diversity of the Landrace population in the Netherlands, although a large proportion of the original variation is maintained. Thus, our recommendation is to conserve breeding lines in a gene bank before they are merged.

Associations of ORMDL1 gene copy number variations with growth traits in four Chinese sheep breeds

Copy number variations (CNVs) are gains and losses of genomic sequence of more than 50 bp between two individuals of a species. Also, CNV is considered to be one of the main elements affecting the phenotypic diversity and evolutionary adaptation of animals. ORMDL sphingolipid biosynthesis regulator 1 (ORMDL1) is a protein-coding gene associated with diseases and development. In our study, the polymorphism of ORMDL1 gene copy numbers in four Chinese sheep breeds (abbreviated CK, HU, STH, and LTH) was detected. In addition, we analyzed the transcriptional expression level of ORMDL1 gene in different tissues of sheep and examined the association of ORMDL1 CNV with growth traits. The statistical analysis revealed that ORMDL1 CNV was remarkably correlated with body height, heart girth, and circumference of cannon bone in HU sheep (P<0.05), and there are significant effects on body weight, body height, body length, chest depth, and height of hip cross in STH sheep (P<0.05). In conclusion, our results provide a basis for the relationship between CNV of ORMDL1 gene and sheep growth traits, suggesting that ORMDL1 CNV may be considered a promising marker for the molecular breeding of Chinese sheep.

Copy number variation of MYLK4 gene and its growth traits of Capra hircus (goat)

Copy number variation (CNV) is a form of genetic variation caused by genome rearrangement, with abnormal fragments ranging from 50 bp to Mb. And, CNV is closely related to disease, growth and reproductive shape of livestock. As a member of myosin light chain kinase (MYLK) family with serine/threonine specificity, MYLK4 belongs to an enzyme encoded by MYLK4 gene. Although MYLK4 is a recognized kinase, its function has yet to be revealed in subsequent studies. This study aims to analyze CNV and genetic effects of MYLK4 gene in goats. We used qPCR to detect CNV of MYLK4 gene in African Nubian goat (n = 32), Guizhou black goat (n = 196) and Guizhou white goat (n = 95), respectively, and correlated CNV data of MYLK4 gene with goat growth traits in Chinese goats. The results showed that the effect of MYLK4 gene CNV on body weight, body length and body height of goats had significantly different (p < 0.05, Q < 0.05), in which CNV showed better growth traits in type of deletion. Therefore, CNV of MYLK4 gene can be used as a molecular marker for assisted selection of goat growth traits, which provides a theoretical basis for the genetic improvement of goat breeds in China.

About the existence of common determinants of gene expression in the porcine liver and skeletal muscle

BACKGROUND

The comparison of expression QTL (eQTL) maps obtained in different tissues is an essential step to understand how gene expression is genetically regulated in a context-dependent manner. In the current work, we have compared the transcriptomic and eQTL profiles of two porcine tissues (skeletal muscle and liver) which typically show highly divergent expression profiles, in 103 Duroc pigs genotyped with the Porcine SNP60 BeadChip (Illumina) and with available microarray-based measurements of hepatic and muscle mRNA levels. Since structural variation could have effects on gene expression, we have also investigated the co-localization of cis-eQTLs with copy number variant regions (CNVR) segregating in this Duroc population.

RESULTS

The analysis of differential expresssion revealed the existence of 1204 and 1490 probes that were overexpressed and underexpressed in the gluteus medius muscle when compared to liver, respectively (|fold-change| > 1.5, q-value < 0.05). By performing genome scans in 103 Duroc pigs with available expression and genotypic data, we identified 76 and 28 genome-wide significant cis-eQTLs regulating gene expression in the gluteus medius muscle and liver, respectively. Twelve of these cis-eQTLs were shared by both tissues (i.e. 42.8% of the cis-eQTLs identified in the liver were replicated in the gluteus medius muscle). These results are consistent with previous studies performed in humans, where 50% of eQTLs were shared across tissues. Moreover, we have identified 41 CNVRs in a set of 350 pigs from the same Duroc population, which had been genotyped with the Porcine SNP60 BeadChip by using the PennCNV and GADA softwares, but only a small proportion of these CNVRs co-localized with the cis-eQTL signals.

CONCLUSION

Despite the fact that there are considerable differences in the gene expression patterns of the porcine liver and skeletal muscle, we have identified a substantial proportion of common cis-eQTLs regulating gene expression in both tissues. Several of these cis-eQTLs influence the mRNA levels of genes with important roles in meat (CTSF) and carcass quality (TAPT1), lipid metabolism (TMEM97) and obesity (MARC2), thus evidencing the practical importance of dissecting the genetic mechanisms involved in their expression.

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

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

Distribution and association study in copy number variation of KCNJ12 gene across four Chinese cattle populations

Copy number variation is a large genome variation which usually happens in the noncoding-region, and it may occur at the locus associated with the functional gene to further influence the phenotype. Potassium inwardly-rectifying channel, subfamily J 12 (KCNJ12) gene expressed widely in cardiomyocytes and neurons, plays an important role in tumor therapy and muscle movement regulation. In this study, we detected the distribution of CNVs for KCNJ12 gene in 404 individuals belonging to four Chinese cattle breeds (NY, JX, JA and GF). We also investigated the KCNJ12 gene expression in different tissues of JX cattle. Additionally, we examined the association of two CNV regions (CNV1: 1,600 bp, intron 1; CNV2: 4,800 bp, intergenic) with growth traits. The statistical analyses indicated that the CNV1 is associated with the body length, rump length and weight in JX cattle population (P < 0.05); and there has a significant association with the body length, chest circumference, and body weight in GF cattle (P < 0.05).The CNV2 had a significant effect on the body length and body weight in JX cattle (P < 0.05); the body length, chest circumference, rump length and body weight in GF cattle (P < 0.01 or P < 0.05). The copy numbers of KCNJ12 gene presented the negative correlations with the transcript level of gene in skeletal muscles (P < 0.05). Our results provide evidence that CNV1 and CNV 2 in KCNJ12 are associated with growth traits in two cattle populations and may be used as candidates for marker-assisted selection and breeding management in cattle.

RAPID COMMUNICATION: A haplotype information theory method reveals genes of evolutionary interest in European vs. Asian pigs

Asian and European wild boars were independently domesticated ca. 10,000 yr ago. Since the 17th century, Chinese breeds have been imported to Europe to improve the genetics of European animals by introgression of favorable alleles, resulting in a complex mosaic of haplotypes. To interrogate the structure of these haplotypes further, we have run a new haplotype segregation analysis based on information theory, namely compression efficiency (CE). We applied the approach to sequence data from individuals from each phylogeographic region (n = 23 from Asia and Europe) including a number of major pig breeds. Our genome-wide CE is able to discriminate the breeds in a manner reflecting phylogeography. Furthermore, 24,956 nonoverlapping sliding windows (each comprising 1,000 consecutive SNP) were quantified for extent of haplotype sharing within and between Asia and Europe. The genome-wide distribution of extent of haplotype sharing was quite different between groups. Unlike European pigs, Asian pigs haplotype sharing approximates a normal distribution. In line with this, we found the European breeds possessed a number of genomic windows of dramatically higher haplotype sharing than the Asian breeds. Our CE analysis of sliding windows captures some of the genomic regions reported to contain signatures of selection in domestic pigs. Prominent among these regions, we highlight the role of a gene encoding the mitochondrial enzyme LACTB which has been associated with obesity, and the gene encoding MYOG a fundamental transcriptional regulator of myogenesis. The origin of these regions likely reflects either a population bottleneck in European animals, or selective targets on commercial phenotypes reducing allelic diversity in particular genes and/or regulatory regions.

A genome-wide association study of growth and fatness traits in two pig populations with different genetic backgrounds

Improvement in growth and fatness traits are the main objectives in pig all breeding programs. Tenth rib backfat thickness (10RIBBFT) and days to 100 kg (D100), which are good predictors of carcass lean content and growth rate, respectively, are economically important traits and also main breeding target traits in pigs. To investigate the genetic mechanisms of 10RIBBFT and D100 of pigs, we sampled 1,137 and 888 pigs from 2 Yorkshire populations of American and British origin, respectively, and conducted genome-wide association study (GWAS) through combined analysis and meta-analysis, to identify SNPs associated with 10RIBBFT and D100. A total of 11 and 7 significant SNPs were identified by combined analysis for 10RIBBFT and D100, respectively. And in meta-analysis, 8 and 7 significant SNPs were identified for 10RIBBFT and D100, respectively. Among them, 6 and 5 common significant SNPs in two analysis results were, respectively, identified associated with 10RIBBFT and D100, and correspondingly explained 2.09% and 0.52% of the additive genetic variance of 10RIBBFT and D100. Further bioinformatics analysis revealed 10 genes harboring or close to these common significant SNPs, 5 for 10RIBBFT and 5 for D100. In particular, Gene Ontology analysis highlighted 6 genes, PCK1, ANGPTL3, EEF1A2, TNFAIP8L3, PITX2, and PLA2G12, as promising candidate genes relevant with backfat thickness and growth. PCK1, ANGPTL3, EEF1A2, and TNFAIP8L3 could influence backfat thickness through phospholipid transport, regulation of lipid metabolic process through the glycerophospholipid biosynthesis and metabolism pathway, the metabolism of lipids and lipoproteins pathway. PITX2 has a crucial role in skeletal muscle tissue development and animal organ morphogenesis, and PLA2G12A plays a role in the lipid catabolic and phospholipid catabolic processes, which both are involved in the body weight pathway. All these candidate genes could directly or indirectly influence fat production and growth in Yorkshire pigs. Our findings provide novel insights into the genetic basis of growth and fatness traits in pigs. The candidate genes for D100 and 10RIBBFT are worthy of further investigation.

Genome-wide association studies for seven production traits highlight genomic regions useful to dissect dry-cured ham quality and production traits in Duroc heavy pigs

Protected designation of origin dry-cured hams are obtained from heavy pigs (slaughtered at about 160 kg of live weight). A specific breeding program designed to improve meat quality for this production has included as key traits the level of intermuscular fat between the leg muscles and ham weight loss during the seasoning period together with a balance between fat and lean cuts. In this study we carried out genome-wide association studies for seven traits used in the genetic merit of Italian Duroc heavy pigs, five related to meat and carcass quality traits (visible intermuscular fat, ham weight loss at first salting, backfat thickness, ham weight and lean cuts), and two related to performance and efficiency traits (average daily gain and feed : gain ratio). A total of 573 performance-tested pigs were genotyped with the Illumina PorcineSNP60 BeadChip and genome-wide association analyses were carried out using the Bayes B approach with the 1 Mb window option of GenSel and random residuals for each of the seven traits. Detected windows were supported by independent single nucleotide polymorphism analyses with a linear mixed model (LMM) approach on the same animals for the same traits. A total of 30 windows identifying different quantitative trait loci (QTL) were detected and among those, 27 were confirmed by LMM in one of these traits. Among the confirmed windows, three QTL were reported for visible intermuscular fat, seven for ham weight loss at first salting and five and four for backfat thickness and lean cut, respectively. A total of eight QTL were detected for the other production traits. No overlapping QTL were reported except for one window on porcine chromosome 10 between lean cuts and ham weight that contained the CACNB2 gene that has been already associated with loin marbling score in other Duroc pigs. Several regions contained genes that have been already associated with production traits in other pig breeds, including Duroc lines, related to fat deposition or muscle structure. This work reports, for the first time, genome-wide association study results for several traits in Italian Duroc heavy pigs. These results will be useful to dissect the genetic basis for dry-cured ham production traits that determine the total genetic merit index of Italian Duroc pigs.

Detection of breed-specific copy number variations in domestic chicken genome

Copy number variations (CNVs) are important large-scale variants. They are widespread in the genome and may contribute to phenotypic variation. Detection and characterization of CNVs can provide new insights into the genetic basis of important traits. Here, we perform whole-genome short read sequence analysis to identify CNVs in two indigenous and commercial chicken breeds to evaluate the impact of the identified CNVs on breed-specific traits. After filtration, a total of 12 955 CNVs spanning (on average) about 9.42% of the chicken genome were found that made up 5467 CNV regions (CNVRs). Chicken quantitative trait loci (QTL) datasets and Ensembl gene annotations were used as resources for the estimation of potential phenotypic effects of our CNVRs on breed-specific traits. In total, 34% of our detected CNVRs were also detected in earlier CNV studies. These CNVRs partly overlap several previously reported QTL and gene ontology terms associated with some important traits, including shank length QTL in Creeper-specific CNVRs and body weight and egg production characteristics, as well as muscle and body organ growth, in the Arian commercial breed. Our findings provide new insights into the genomic structure of the chicken genome for an improved understanding of the potential roles of CNVRs in differentiating between breeds or lines.

Identification of new single nucleotide polymorphisms affecting total number born and candidate genes related to ovulation rate in Chinese Erhualian pigs

The Chinese Erhualian pig has the highest record for litter size in the world. However, the genetic mechanism of its high prolificacy remains poorly understood. In our study, large phenotypic variations in litter size were found among Erhualian sows. Significant differences in total number born (TNB) and corpora lutea numbers were observed between sows with high and low estimated breeding values (EBVs) for TNB. To identify single nucleotide polymorphisms (SNPs) associated with TNB, a selective genomic scan was conducted on 18 sows representing the top 10% and 18 sows representing the bottom 10% of EBVs of 177 sows using Illumina Porcine SNP60 genotype data. Genome-wide fixation coefficient (F_ST ) values were calculated for each SNP between the high- and low-EBV groups. A total of 154 SNPs were significantly differentiated loci between the two groups. Of the top 10 highest F_ST SNPs, rs81399474, rs81400131 and rs81405013 on SSC8 and rs81434499 and rs81434489 on SSC 12 corresponded to previously reported QTL for litter size. The other five SNPs, rs81367039 on SSC2, rs80891106 on SSC7, rs81477883 on SSC12 and rs80938898 and rs80971725 on SSC14, appeared to be novel QTL for TNB. Significant associations between rs81399474 on SSC8 and TNB were confirmed in 313 Erhualian sows. Forty genes were identified around the top 10 highest F_ST SNPs, of which UCHL1, adjacent to rs81399474, and RPS6KB1 and CLTC, adjacent to rs81434499, have been reported to affect the ovulation rate in pig. The findings can advance understanding of the genetic variations in litter size of pigs.

Transcriptome analysis revealed chimeric RNAs, single nucleotide polymorphisms and allele-specific expression in porcine prenatal skeletal muscle

Prenatal skeletal muscle development genetically determines postnatal muscle characteristics such as growth and meat quality in pigs. However, the molecular mechanisms underlying prenatal skeletal muscle development remain unclear. Here, we performed the first genome-wide analysis of chimeric RNAs, single nuclear polymorphisms (SNPs) and allele-specific expression (ASE) in prenatal skeletal muscle in pigs. We identified 14,810 protein coding genes and 163 high-confidence chimeric RNAs expressed in prenatal skeletal muscle. More than 94.5% of the chimeric RNAs obeyed the canonical GT/AG splice rule and were trans-splicing events. Ten and two RNAs were aligned to human and mouse chimeric transcripts, respectively. We detected 106,457 high-quality SNPs (6,955 novel), which were mostly (89.09%) located within QTLs for production traits. The high proportion of non-exonic SNPs revealed the incomplete annotation status of the current swine reference genome. ASE analysis revealed that 11,300 heterozygous SNPs showed allelic imbalance, whereas 131 ASE variants were located in the chimeric RNAs. Moreover, 4 ASE variants were associated with various economically relevant traits of pigs. Taken together, our data provide a source for studies of chimeric RNAs and biomarkers for pig breeding, while illuminating the complex transcriptional events underlying prenatal skeletal muscle development in mammals.

SNP- and haplotype-based genome-wide association studies for growth, carcass, and meat quality traits in a Duroc multigenerational population

Background

The aim of the present study was to compare the power of single nucleotide polymorphism (SNP)-based genome-wide association study (GWAS) and haplotype-based GWAS for quantitative trait loci (QTL) detection, and to detect novel candidate genes affecting economically important traits in a purebred Duroc population comprising seven-generation pedigree. First, we performed a simulation analysis using real genotype data of this population to compare the power (based on the null hypothesis) of the two methods. We then performed GWAS using both methods and real phenotype data comprising 52 traits, which included growth, carcass, and meat quality traits.

Results

In total, 836 animals were genotyped using the Illumina PorcineSNP60 BeadChip and 14 customized SNPs from regions of known candidate genes related to the traits of interest. The power of SNP-based GWAS was greater than that of haplotype-based GWAS in a simulation analysis. In real data analysis, a larger number of significant regions was obtained by SNP-based GWAS than by haplotype-based GWAS. For SNP-based GWAS, 23 genome-wide significant SNP regions were detected for 17 traits, and 120 genome-wide suggestive SNP regions were detected for 27 traits. For haplotype-based GWAS, 6 genome-wide significant SNP regions were detected for four traits, and 11 genome-wide suggestive SNP regions were detected for eight traits. All genome-wide significant SNP regions detected by haplotype-based GWAS were located in regions also detected by SNP-based GWAS. Four regions detected by SNP-based GWAS were significantly associated with multiple traits: on Sus scrofa chromosome (SSC) 1 at 304 Mb; and on SSC7 at 35–39 Mb, 41–42 Mb, and 103 Mb. The vertnin gene (VRTN) in particular, was located on SSC7 at 103 Mb and was significantly associated with vertebrae number and carcass lengths. Mapped QTL regions contain some candidate genes involved in skeletal formation (FUBP3; far upstream element binding protein 3) and fat deposition (METTL3; methyltransferase like 3).

Conclusion

Our results show that a multigenerational pig population is useful for detecting QTL, which are typically segregated in a purebred population. In addition, a novel significant region could be detected by SNP-based GWAS as opposed to haplotype-based GWAS.

Electronic supplementary material

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

Genetics of Adiposity in Large Animal Models for Human Obesity-Studies on Pigs and Dogs

The role of domestic mammals in the development of human biomedical sciences has been widely documented. Among these model species the pig and dog are of special importance. Both are useful for studies on the etiology of human obesity. Genome sequences of both species are known and advanced genetic tools [eg, microarray SNP for genome wide association studies (GWAS), next generation sequencing (NGS), etc.] are commonly used in such studies. In the domestic pig the accumulation of adipose tissue is an important trait, which influences meat quality and fattening efficiency. Numerous quantitative trait loci (QTLs) for pig fatness traits were identified, while gene polymorphisms associated with these traits were also described. The situation is different in dog population. Generally, excessive accumulation of adipose tissue is considered, similar to humans, as a complex disease. However, research on the genetic background of canine obesity is still in its infancy. Between-breed differences in terms of adipose tissue accumulation are well known in both animal species. In this review we show recent advances of studies on adipose tissue accumulation in pigs and dogs, and their potential importance for studies on human obesity.

Identification of Copy Number Variations in Xiang and Kele Pigs

Xiang and Kele pigs are two well-known local Chinese pig breeds that possess rich genetic resources and have enormous economic and scientific value. We performed a comprehensive genomic analysis of the copy number variations (CNVs) in these breeds. CNVs are one of the most important forms of genomic variation and have profound effects on phenotypic variation. In this study, PorcineSNP60 genotyping data from 98 Xiang pigs and 22 Kele pigs were used to identify CNVs. In total, 172 candidate CNV regions (CNVRs) were identified, ranging from 3.19 kb to 8175.26 kb and covering 80.41 Mb of the pig genome. Approximately 56.40% (97/172) of the CNVRs overlapped with those identified in seven previous studies, and 43.60% (75/172) of the identified CNVRs were novel. Of the identified CNVRs, 82 (47 gain, 33 loss, and two gain-loss events that covered 4.58 Mb of the pig genome) were found only in a Xiang population with a large litter size. In contrast, 13 CNVRs (8 gain and 5 loss events) were unique to a Xiang population with small litter sizes, and 30 CNVRs (14 loss and 16 gain events) were unique to Kele pigs. The CNVRs span approximately 660 annotated Sus scrofa genes that are significantly enriched for specific biological functions, such as sensory perception, cognition, reproduction, ATP biosynthetic processes, and neurological processes. Many CNVR-associated genes, particularly the genes involved in reproductive traits, differed between the Xiang populations with large and small litter sizes, and these genes warrant further investigation due to their importance in determining the reproductive performance of Xiang pigs. Our results provide meaningful information about genomic variation, which may be useful in future assessments of the associations between CNVs and important phenotypes in Xiang and Kele pigs to ultimately help protect these rare breeds.

Genome-wide Association Study to Identify Quantitative Trait Loci for Meat and Carcass Quality Traits in Berkshire

Meat and carcass quality attributes are of crucial importance influencing consumer preference and profitability in the pork industry. A set of 400 Berkshire pigs were collected from Dasan breeding farm, Namwon, Chonbuk province, Korea that were born between 2012 and 2013. To perform genome wide association studies (GWAS), eleven meat and carcass quality traits were considered, including carcass weight, backfat thickness, pH value after 24 hours (pH24), Commission Internationale de l’Eclairage lightness in meat color (CIE L), redness in meat color (CIE a), yellowness in meat color (CIE b), filtering, drip loss, heat loss, shear force and marbling score. All of the 400 animals were genotyped with the Porcine 62K SNP BeadChips (Illumina Inc., USA). A SAS general linear model procedure (SAS version 9.2) was used to pre-adjust the animal phenotypes before GWAS with sire and sex effects as fixed effects and slaughter age as a covariate. After fitting the fixed and covariate factors in the model, the residuals of the phenotype regressed on additive effects of each single nucleotide polymorphism (SNP) under a linear regression model (PLINK version 1.07). The significant SNPs after permutation testing at a chromosome-wise level were subjected to stepwise regression analysis to determine the best set of SNP markers. A total of 55 significant (p<0.05) SNPs or quantitative trait loci (QTL) were detected on various chromosomes. The QTLs explained from 5.06% to 8.28% of the total phenotypic variation of the traits. Some QTLs with pleiotropic effect were also identified. A pair of significant QTL for pH24 was also found to affect both CIE L and drip loss percentage. The significant QTL after characterization of the functional candidate genes on the QTL or around the QTL region may be effectively and efficiently used in marker assisted selection to achieve enhanced genetic improvement of the trait considered.

Genome-Wide Copy Number Variations Using SNP Genotyping in a Mixed Breed Swine Population

Copy number variations (CNVs) are increasingly understood to affect phenotypic variation. This study uses SNP genotyping of trios of mixed breed swine to add to the catalog of known genotypic variation in an important agricultural animal. PorcineSNP60 BeadChip genotypes were collected from 1802 pigs that combined to form 1621 trios. These trios were from the crosses of 50 boars with 525 sows producing 1621 piglets. The pigs were part of a population that was a mix of ¼ Duroc, ½ Landrace and ¼ Yorkshire breeds. Merging the overlapping CNVs that were observed in two or more individuals to form CNV regions (CNVRs) yielded 502 CNVRs across the autosomes. The CNVRs intersected genes, as defined by RefSeq, 84% of the time – 420 out of 502. The results of this study are compared and contrasted to other swine studies using similar and different methods of detecting CNVR. While progress is being made in this field, more work needs to be done to improve consistency and confidence in CNVR results.

Genome-wide association study on legendre random regression coefficients for the growth and feed intake trajectory on Duroc Boars

Background

Feed intake and growth are economically important traits in swine production. Previous genome wide association studies (GWAS) have utilized average daily gain or daily feed intake to identify regions that impact growth and feed intake across time. The use of longitudinal models in GWAS studies, such as random regression, allows for SNPs having a heterogeneous effect across the trajectory to be characterized. The objective of this study is therefore to conduct a single step GWAS (ssGWAS) on the animal polynomial coefficients for feed intake and growth.

Results

Corrected daily feed intake (DFI_Adj) and average daily weight measurements (DBW_Avg) on 8981 (n = 525,240 observations) and 5643 (n = 283,607 observations) animals were utilized in a random regression model using Legendre polynomials (order = 2) and a relationship matrix that included genotyped and un-genotyped animals. A ssGWAS was conducted on the animal polynomials coefficients (intercept, linear and quadratic) for animals with genotypes (DFI_Adj: n = 855; DBW_Avg: n = 590). Regions were characterized based on the variance of 10-SNP sliding windows GEBV (WGEBV). A bootstrap analysis (n =1000) was conducted to declare significance. Heritability estimates for the traits trajectory ranged from 0.34-0.52 to 0.07-0.23 for DBW_Avg and DFI_Adj, respectively. Genetic correlations across age classes were large and positive for both DBW_Avg and DFI_Adj, albeit age classes at the beginning had a small to moderate genetic correlation with age classes towards the end of the trajectory for both traits. The WGEBV variance explained by significant regions (P < 0.001) for each polynomial coefficient ranged from 0.2-0.9 to 0.3-1.01 % for DBW_Avg and DFI_Adj, respectively. The WGEBV variance explained by significant regions for the trajectory was 1.54 and 1.95 % for DBW_Avg and DFI_Adj. Both traits identified candidate genes with functions related to metabolite and energy homeostasis, glucose and insulin signaling and behavior.

Conclusions

We have identified regions of the genome that have an impact on the intercept, linear and quadratic terms for DBW_Avg and DFI_Adj. These results provide preliminary evidence that individual growth and feed intake trajectories are impacted by different regions of the genome at different times.

Electronic supplementary material

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

A genomewide association mapping study using ultrasound-scanned information identifies potential genomic regions and candidate genes affecting carcass traits in Nellore cattle

The aim of this study was to identify candidate genes and genomic regions associated with ultrasound-derived measurements of the rib-eye area (REA), backfat thickness (BFT) and rumpfat thickness (RFT) in Nellore cattle. Data from 640 Nellore steers and young bulls with genotypes for 290 863 single nucleotide polymorphisms (SNPs) were used for genomewide association mapping. Significant SNP associations were explored to find possible candidate genes related to physiological processes. Several of the significant markers detected were mapped onto functional candidate genes including ARFGAP3, CLSTN2 and DPYD for REA; OSBPL3 and SUDS3 for BFT; and RARRES1 and VEPH1 for RFT. The physiological pathway related to lipid metabolism (CLSTN2, OSBPL3, RARRES1 and VEPH1) was identified. The significant markers within previously reported QTLs reinforce the importance of the genomic regions, and the other loci offer candidate genes that have not been related to carcass traits in previous investigations.

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.

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.

Copy number variants in Italian Large White pigs detected using high-density single nucleotide polymorphisms and their association with back fat thickness

The aim of this study was to identify copy number variants (CNVs) in Italian Large White pigs and test them for association with back fat thickness (BFT). Within a population of 12 000 performance-tested pigs, two groups of animals with extreme and divergent BFT estimated breeding values (EBVs; 147 with negative and 150 with positive EBVs) were genotyped with the Illumina Porcine SNP60 BeadChip. CNVs were detected with PENNCNV software. We identified a total of 4146 CNV events in 170 copy number variation regions (CNVRs) located on 15 porcine autosomes. Validation of detected CNVRs was carried out (i) by comparing CNVRs already detected by other studies and (ii) by semiquantitative fluorescent multiplex (SQFM) PCR of a few CNVRs. Most of CNVRs detected in Italian Large White pigs (71.2%) were already reported in other pig breeds/populations, and 82.1% of the CNV events detected by PENNCNV were confirmed by SQFM PCR. For each CNVR, we compared the occurrence of CNV events between the pigs of the high and low BFT EBV tails. Sixteen regions showed significance at P < 0.10, and seven were significant at P < 0.05 but were not significant after Bonferroni correction (Fisher's exact test). These results indicated that CNVs could explain a limited fraction of the genetic variability of fat deposition in Italian Large White pigs. However, it was interesting to note that one of these CNVRs encompassed the ZPLD1 gene. In humans, a rare CNV event including this gene is associated with obesity. Studies identifying CNVs in pigs could assist in elucidating the genetic mechanisms underlying human obesity.