A genome-wide association study of reproductive traits in a Yorkshire pig population

Genome-wide single nucleotide polymorphism (SNP) data reveal potential candidate genes for litter traits in a Yorkshire pig population

The litter trait is one of the most important economic traits, and increasing litter size is of great economic value in the pig industry. However, the molecular mechanisms underlying pig litter traits remain elusive. To identify molecular markers and candidate genes for pig litter traits, a genome-wide association study (GWAS) and selection signature analysis were conducted in a Yorkshire pig population. A total of 518 producing sows were genotyped with Illumina Porcine SNP 50 BeadChip, and 1969 farrowing records for the total number born (TNB), the number born alive (NBA), piglets born dead (PBD), and litter weight born alive (LWB) were collected. Then, a GWAS was performed for the four litter traits using a repeatability model. Based on the estimated breeding values (EBVs) of TNB, 15 high- and 15 low-prolificacy individuals were selected from the 518 sows to implement selection signature analysis. Subsequently, the selection signatures affecting the litter traits of sows were detected by using two methods including the fixation index (FST) and θ π . Combining the results of the GWAS and selection signature analysis, 20 promising candidate genes (NKAIN2, IGF1R, KISS1R, TYRO3, SPINT1, ADGRF5, APC2, PTBP1, CLCN3, CBR4, HPF1, FAM174A, SCP2, CLIC1, ZFYVE9, SPATA33, KIF5C, EPC2, GABRA2, and GABRA4) were identified. These findings provide novel insights into the genetic basis of pig litter traits and will be helpful for improving the reproductive performances of sows in pig breeding.

Elucidation of the genetic determination of clutch traits in Chinese local chickens of the Laiwu Black breed

BACKGROUND

Egg laying rate (LR) is associated with a clutch, which is defined as consecutive days of oviposition. The clutch trait can be used as a selection indicator to improve egg production in poultry breeding. However, little is known about the genetic basis of clutch traits. In this study, our aim was to estimate genetic parameters and identify quantitative trait single nucleotide polymorphisms for clutch traits in 399 purebred Laiwu Black chickens (a native Chinese breed) using a genome-wide association study (GWAS).

METHODS

In this work, after estimating the genetic parameters of age at first egg, body weight at first egg, LR, longest clutch until 52 week of age, first week when the longest clutch starts, last week when the longest clutch ends, number of clutches, and longest number of days without egg-laying until 52 week of age, we identified single nucleotide polymorphisms (SNPs) and potential candidate genes associated with clutch traits in Laiwu Black chickens. The restricted maximum likelihood method was used to estimate genetic parameters of clutch pattern in 399 Laiwu Black hens, using the GCTA software.

RESULTS

The results showed that SNP-based heritability estimates of clutch traits ranged from 0.06 to 0.59. Genotyping data were obtained from whole genome re-sequencing data. After quality control, a total of 10,810,544 SNPs remained to be analyzed. The GWAS revealed that 421 significant SNPs responsible for clutch traits were scattered on chicken chromosomes 1-14, 17-19, 21-25, 28 and Z. Among the annotated genes, NELL2, SMYD9, SPTLC2, SMYD3 and PLCL1 were the most promising candidates for clutch traits in Laiwu Black chickens.

CONCLUSION

The findings of this research provide critical insight into the genetic basis of clutch traits. These results contribute to the identification of candidate genes and variants. Genes and SNPs potentially provide new avenues for further research and would help to establish a framework for new methods of genomic prediction, and increase the accuracy of estimated genetic merit for egg production and clutch traits.

The impact of Box-Cox transformation on phenotypic and genomic characteristics of litter size variability in Landrace pigs

The genetic background of variability remains of interest especially in traits of high economic importance, e.g. litter size in pigs. It has been indicated that the data transformation can affect the variability phenotype. This study aims to evaluate the phenotypic and genomic background of variability of litter size obtained from data before and after the Box-Cox transformation. In total, 67 500 records on the total number born (TNB) in Landrace pig population were used. Since the data presented skewness, the decision was made to perform Box-Cox transformation on TNB and obtain bcTNB. Next, the phenotypic variability was estimated as log-transformed variance of residuals (LnVar) for both TNB (LnVar_TNB) and bcTNB (LnVar_bcTNB). The variability traits were further used in the genome-wide association study (GWAS) performed on 10 688 sows genotyped with Axiom porcine 660 K or imputed to 660 K SNP-chip. The substantial difference in skewness was observed after data transformation, represented as a change from -0.46 to -0.02. Heritability for TNB was 0.118 vs 0.125 for bcTNB. The heritability for LnVar_TNB was 0.0025 vs 0.0037 for LnVar_bcTNB. The change in the genetic variance was confirmed when genetic coefficients on SD level were compared: 2% for LnVar_TNB vs 4% for LnVar_bcTNB. In bivariate analysis, the genetic correlation between the additive genetic effects of the mean TNB and its variability changed from 0.38 to 0.63. The observed positive genetic correlations indicated that selection focused on increasing the litter size will simultaneously cause an increase in litter size variability. Based on GWAS, 14 SNPs were detected for LnVar_TNB and eight for LnVar_bcTNB, with two of them indicating the most promising candidate genes. First candidate gene located on Sus scrofa chromosome (SSC) 3 is STAG3, which plays an essential role in gametogenesis. Second gene located on SSC 10 is ESRRG, which affects placenta development. The additional post-GWAS analysis indicated even more candidate genes for LnVar_TNB and LnVar_bcTNB. The most promising candidate gene was located on SSC 13 - MFN1, which is involved in embryonic development. The results of this study indicated a substantial change in variance components for variability when the Box-Cox transformation was applied to data presenting skewness. Moreover, the data transformation changed the phenotype substantially enough that only part of SNP overlapped between two variability traits. Our investigation shows that it is essential to perform Box-Cox transformation for skewed data in order to properly describe phenotypic and genomic properties of litter size variability in Landrace pigs.

Further host-genomic characterization of total antibody response to PRRSV vaccination and its relationship with reproductive performance in commercial sows: genome-wide haplotype and zygosity analyses

Background

The possibility of using antibody response (S/P ratio) to PRRSV vaccination measured in crossbred commercial gilts as a genetic indicator for reproductive performance in vaccinated crossbred sows has motivated further studies of the genomic basis of this trait. In this study, we investigated the association of haplotypes and runs of homozygosity (ROH) and heterozygosity (ROHet) with S/P ratio and their impact on reproductive performance.

Results

There was no association (P-value ≥ 0.18) of S/P ratio with the percentage of ROH or ROHet, or with the percentage of heterozygosity across the whole genome or in the major histocompatibility complex (MHC) region. However, specific ROH and ROHet regions were significantly associated (P-value ≤ 0.01) with S/P ratio on chromosomes 1, 4, 5, 7, 10, 11, 13, and 17 but not (P-value ≥ 0.10) with reproductive performance. With the haplotype-based genome-wide association study (GWAS), additional genomic regions associated with S/P ratio were identified on chromosomes 4, 7, and 9. These regions harbor immune-related genes, such as SLA-DOB, TAP2, TAPBP, TMIGD3, and ADORA. Four haplotypes at the identified region on chromosome 7 were also associated with multiple reproductive traits. A haplotype significantly associated with S/P ratio that is located in the MHC region may be in stronger linkage disequilibrium (LD) with the quantitative trait loci (QTL) than the previously identified single nucleotide polymorphism (SNP) (H3GA0020505) given the larger estimate of genetic variance explained by the haplotype than by the SNP.

Conclusions

Specific ROH and ROHet regions were significantly associated with S/P ratio. The haplotype-based GWAS identified novel QTL for S/P ratio on chromosomes 4, 7, and 9 and confirmed the presence of at least one QTL in the MHC region. The chromosome 7 region was also associated with reproductive performance. These results narrow the search for causal genes in this region and suggest SLA-DOB and TAP2 as potential candidate genes associated with S/P ratio on chromosome 7. These results provide additional opportunities for marker-assisted selection and genomic selection for S/P ratio as genetic indicator for litter size in commercial pig populations.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12711-021-00676-5.

Genetic Analysis of Major Production and Reproduction Traits of Korean Duroc, Landrace and Yorkshire Pigs

The study aimed to investigate the genetic parameters of the production and reproduction traits of Korean Duroc, Landrace, and Yorkshire pigs. Three production traits, namely average daily gain (ADG), age at 105 kg body weight (DAYS105) and backfat thickness (BFT), and three reproduction traits, namely age at first farrowing (AFF), total number of piglets born (TNB) and number of piglets born alive (NBA), were analyzed. The reproduction dataset was based on first-parity gilts only. However, the production dataset comprised pigs of both sexes. Genetic parameters were estimated from individual datasets using a multiple-trait animal model in AIREMLF90 software. The heritability values of ADG, DAYS105 and BFT were 0.34-0.36, 0.41-0.44 and 0.38-0.48, respectively, across breeds. Heritability values for AFF, TNB and NBA were 0.07-0.14, 0.09-0.11 and 0.09-0.10, respectively. Strong genetic correlations existed between ADG and DAYS105 (-0.97) and between TNB and NBA (0.90 to 0.96). In line with breeding goals, all productive traits in Duroc pigs, and all reproduction traits except AFF in Landrace and Yorkshire pigs, also showed noticeable improvements in recent years. In conclusion, we believe that our findings on economic traits would greatly assist future pig breeding decisions in Korea.

Integrated analysis of lncRNA, miRNA and mRNA reveals novel insights into the fertility regulation of large white sows

Background

Improving sow fertility is extremely important as it can lead to increased reproductive efficiency and thus profitability for swine producers. There are considerable differences in fertility rates among individual animals, but the underlying molecular mechanisms remain unclear. In this study, by using different types of RNA libraries, we investigated the complete transcriptome of ovarian tissue during the luteal (L) and follicular (F) phases of the estrous cycle in Large White pigs with high (H) and low (L) fecundity, and performed a comprehensive analysis of long noncoding RNAs (lncRNAs), mRNAs and micro RNAs (miRNAs) from 16 samples by combining RNA sequencing (RNA-seq) with bioinformatics.

Results

In total, 24,447 lncRNAs, 27,370 mRNAs, and 216 known miRNAs were identified in ovarian tissues. The genomic features of lncRNAs, such as length distribution and number of exons, were further analyzed. We selected a threshold of P < 0.05 and |log₂ (fold change)| ≥ 1 to obtain the differentially expressed lncRNAs, miRNAs and mRNAs by pairwise comparison (LH vs. LL, FH vs. FL). Bioinformatics analysis of these differentially expressed RNAs revealed multiple significantly enriched pathways (P < 0.05) that were closely involved in the reproductive process, such as ovarian steroidogenesis, lysosome, steroid biosynthesis, and the estrogen and GnRH signaling pathways. Moreover, bioinformatics screening of differentially expressed miRNAs that share common miRNA response elements (MREs) with lncRNAs and their downstream mRNA targets were performed. Finally, we constructed lncRNA–miRNA–mRNA regulation networks. The key genes in these networks were verified by Reverse Transcription Real-time Quantitative PCR (RT-qRCR), which were consistent with the results from RNA-Seq data.

Conclusions

These results provide further insights into the fertility of pigs andcan contribute to further experimental investigation of the functions of these genes.

Transcriptomic comparison of liver tissue between Anqing six-end-white pigs and Yorkshire pigs based on RNA sequencing

Chinese indigenous pig and Western commercial pig breeds show different patterns of lipid metabolism, fat deposition, and fatty acid composition; for these reasons, they have become vitally important models of energy metabolism and obesity in humans. To compare the mechanisms underlying lipid metabolism between Yorkshire pigs (lean type) and Anqing six-end-white pigs (obese type), the liver transcriptomes of six castrated boars with a body weight of approximately 100 kg (three Yorkshire and three Anqing) were analyzed by RNA-seq. The total number of reads produced for each liver sample ranged from 47.05 to 62.6 million. Among 362 differentially expressed genes, 142 were up-regulated and 220 were down-regulated in Anqing six-end-white pigs. Based on these data, 79 GO terms were significantly enriched. The top 10 (the 10 with lowest corrected P-value) significantly enriched GO terms were identified, including lipid metabolic process and carboxylic acid metabolic process. Pathway analysis revealed three significantly enriched KEGG pathways including PPAR signaling pathway, steroid hormone biosynthesis, and retinol metabolism. Based on protein–protein interaction networks, multiple genes responsible for lipid metabolism were identified, such as PCK1, PPARA, and CYP7A1, and these were considered promising candidate genes that could affect porcine liver lipid metabolism and fat deposition. Our results provide abundant transcriptomic information that will be useful for animal breeding and biomedical research.

Effect of single-nucleotide polymorphisms on specific reproduction parameters in Hungarian Large White sows

The aim of this study was to reveal the effect of single-nucleotide polymorphisms (SNPs) on the total number of piglets born (TNB), the litter weight born alive (LWA), the number of piglets born dead (NBD), the average litter weight on the 21st day (M21D) and the interval between litters (IBL). Genotypes were determined on a high-density Illumina Porcine SNP 60K BeadChip. Data screening and data identification were performed by a multi-locus mixed-model. Statistical analyses were carried out to find associations between individual genotypes of 290 Hungarian Large White sows and the investigated reproduction parameters. According to the analysis outcome, three SNPs were identified to be associated with TNB. These loci are located on chromosomes 1, 6 and 13 (-log₁₀P = 6.0, 7.86 and 6.22, the frequencies of their minor alleles, MAF, were 0.298, 0.299 and 0.364, respectively). Two loci showed considerable association (-log₁₀P = 10.35 and 10.46) with LWA on chromosomes 5 and X, the MAF were 0.425 and 0.446, respectively. Seven loci were found to be associated with NBD. These loci are located on chromosomes 5, 6, 13, 14, 15, 16 and 18 (-log₁₀P = 10.95, 5.43, 8.29, 6.72, 6.81, 5.90, and 5.15, respectively). One locus showed association (-log₁₀P = 5.62) with M21D on chromosome 1 (the MAF was 0.461). Another locus was found to be associated with IBL on chromosome 8 (-log₁₀P = 7.56; the MAF was 0.438). The above-mentioned loci provide a straightforward possibility to assist selection by molecular tools and, consequently, to improve the competitiveness of the Hungarian Large White (HLW) breed.