A genome-wide association study in large white and landrace pig populations for number piglets born alive

Genome-Wide Association Study of Conformation Traits in Brazilian Holstein Cattle

The linear conformation of animals exerts an influence on health, reproduction, production, and welfare, in addition to longevity, which directly affects the profitability of milk-producing farms. The objectives of this study were (1) to perform genome-wide association studies (GWASs) of conformation traits, namely the Rump, Feet and Legs, Mammary System, Dairy Strength, and Final Classification traits, and (2) to identify genes and related pathways involved in physiological processes associated with conformation traits in Brazilian Holstein cattle. Phenotypic and genotypic data from 2339 Holstein animals distributed across the states of Rio Grande do Sul, Paraná, São Paulo, and Minas Gerais were used. The genotypic data were obtained with a 100 K SNP marker panel. The single-step genome-wide association study (ssGWAS) method was employed in the analyses. Genes close to a significant SNP were identified in an interval of 100 kb up- and downstream using the Ensembl database available in the BioMart tool. The DAVID database was used to identify the main metabolic pathways and the STRING program was employed to create the gene regulatory network. In total, 36 significant SNPs were found on 15 chromosomes; 27 of these SNPs were linked to genes that may influence the traits studied. Fourteen genes most closely related to the studied traits were identified, as well as four genes that showed interactions in important metabolic pathways such as myogenesis, adipogenesis, and angiogenesis. Among the total genes, four were associated with myogenesis (TMOD2, TMOD3, CCND2, and CTBP2), three with angiogenesis (FGF23, FGF1, and SCG3), and four with adipogenesis and body size and development (C5H12orf4, CCND2, EMILIN1, and FGF6). These results contribute to a better understanding of the biological mechanisms underlying phenotypic variability in conformation traits in Brazilian Holstein cattle.

Evaluation of genetic differentiation and search for candidate genes for reproductive traits in pigs

OBJECTIVE

The use of molecular genetic methods in pig breeding can significantly increase the efficiency of breeding and breeding work. We applied the Fst (fixsacion index) method, the main focus of the work was on the search for common options related to the number of born piglets and the weight of born piglets, since today the urgent task is to prevent a decrease in the weight of piglets at birth while maintaining high fertility of sows.

METHODS

One approach is to scan the genome, followed by an assessment of Fst and identification of selectively selected regions. We chose Large White sows (n = 237) with the same conditions of keeping and feeding. The data were collected from the sows across three farrowing. For genotyping, we used GeneSeek GGP Porcine HD Genomic Profiler v1, which included 68,516 single nucleotide polymorphisms evenly distributed with an average spacing of 25 kb (Illumina Inc, San Diego, CA, USA).

RESULTS

Based on the results of the Fst analysis, 724 variants representing selection signals for the signs BALWT, BALWT1, NBA, and TNB (weight of piglets born alive, average weight of the 1st piglets born alive, total number born alive, total number born). At the same time, 18 common variants have been identified that are potential markers for both the number of piglets at birth and the weight of piglets at birth, which is extremely important for breeding work to improve reproductive characteristics in sows.

CONCLUSION

Оur work resulted in identification of variants associated with the reproductive characteristics of pigs. Moreover, we identified, variants which are potential markers for both the number of piglets at birth and the weight of piglets at birth, which is extremely important for breeding work to improve reproductive performance in sows.

Genome-wide association studies for the number of piglets born alive and dead in Dongliao black pigs

Litter size (total number born) trait has a great impact on the economic success of pork production. The total number born consists of the number of piglets born alive and dead. To clarify the genetic background of litter size, genome-wide association studies were undertaken in the present study. Samples of DNA were collected and genotyped using the Porcine 50K BeadChip from 723 Dongliao Black sows. Using three different models (BLINK, FarmCPU, and MLM), a total of 155 significant SNPs were discovered, six of which had been reported in previous pig reproduction association studies. We suggest that rs81318434, located in the GLI3 gene, might be the promising candidate affecting litter size trait. Our findings may provide insights for uncovering the genetic mechanisms for the litter size of pigs.

Functional genomics of reproduction in pigs: Are we there yet?

Reproductive failure is the main reason for culling females in swine herds and is both a financial and sustainability issue. Because reproductive traits are complex and lowly to moderately heritable, genomic selection within populations can achieve substantial genetic gain in reproductive efficiency. A better understanding of the physiological components affecting the expression of these traits will facilitate greater understanding of the genes affecting reproductive traits and is necessary to improve and optimize management strategies to maximize reproductive success of gilts and sows. Large-scale genotyping with single-nucleotide polymorphism (SNP) arrays are used for genome-wide association studies (GWAS) and have facilitated identification of positional candidate genes. Transcriptomic data can be used to weight SNP for GWAS and could lead to previously unidentified candidate genes. Resequencing and fine mapping of candidate genes are necessary to identify putative functional variants and some of these have been incorporated into new genotyping arrays. Sequence imputation and genotype by sequence are newer strategies that could reveal novel functional mutations. In this study, these approaches are discussed. Advantages and limitations are highlighted where additional research is needed.

Allelic to genome wide perspectives of swine genetic variation to litter size and its component traits

Litter size is a complex and sex limited trait that depends on various biological, managemental and environmental factors. Owing to its low heritability it is inefficaciously selected by traditional methods. However, due to higher heritability of ovulation rate and embryo survival, selection based on component traits of litter size is advocated. QTL analysis and candidate gene approach are among the various supplementary/alternate strategies for selection of litter size. QTL analysis is aimed at identifying genomic regions affecting trait of interest significantly. Candidate gene approach necessitates identification of genes potentially affecting the trait. There are various genes that significantly affect litter size and its component traits viz. ESR, LEP, BF, IGFBP, RBP4, PRLR, CTNNAL1, WNT10B, TCF12, DAZ, and RNF4. These genes affect litter size in a complex interacting manner. Lately, genome wide association study (GWAS) have been utilized to unveil the genetic and biological background of litter traits, and elucidate the genes governing litter size. Favorable SNPs in these genes have been identified and offers a scope for inclusion in selection programs thereby increasing breeding efficiency and profit in pigs. The review provides a comprehensive coverage of investigations carried out globally to unravel the genetic variation in litter size and its component traits in pigs, both at allelic and genome wide level. It offers a current perspective on different strategies including the profiling of candidate genes, QTLs, and genome wide association studies as an aid to efficient selection for litter size and its component traits.

Identification of Genomic Regions and Candidate Genes for Litter Traits in French Large White Pigs Using Genome-Wide Association Studies

The reproductive traits of sows are one of the important economic traits in pig production, and their performance directly affects the economic benefits of the entire pig industry. In this study, a total of 895 French Large White pigs were genotyped by GeneSeek Porcine 50K SNP Beadchip and four phenotypic traits of 1407 pigs were recorded, including total number born (TNB), number born alive (NBA), number healthy piglets (NHP) and litter weight born alive (LWB). To identify genomic regions and genes for these traits, we used two approaches: a single-locus genome-wide association study (GWAS) and a single-step GWAS (ssGWAS). Overall, a total of five SNPs and 36 genomic regions were identified by single-locus GWAS and ssGWAS, respectively. Notably, fourof all five significant SNPs were located in 10.72–11.06 Mb on chromosome 7, were also identified by ssGWAS. These regions explained the highest or second highest genetic variance in the TNB, NBA and NHP traits and harbor the protein coding gene ENSSSCG00000042180. In addition, several candidate genes associated with litter traits were identified, including JARID2, PDIA6, FLRT2 and DICER1. Overall, these novel results reflect the polygenic genetic architecture of the litter traits and provide a theoretical reference for the following implementation of molecular breeding.

Effect of the polymorphism of genes related to lactogenesis in pigs on colostrum and milk composition and on rearing performance of piglets

The aim of the study was a search for new polymorphisms in the genes related to lactogenesis, including lactose synthesis in colostrum and milk, and determine their relationship with the content of colostrum and milk components and body weight of the piglets. Colostrum and milk were sampled during the second lactation on d 1, 7, 14 and 21 from 112 sows of Polish Large White (PLW) and Polish Landrace (PL). The piglets’ rearing performance was determined based on the number and body weight at 1, 7, 14 and 21 d of age. In the study six polymorphisms located in different genes were analysed. Among the six polymorphisms under analysis, those identified in the B4GALT1 genes had the most significant effect on colostrum and milk composition and the piglets’ body weight changes. Sows of the B4GALT1GG genotype, when compared to those of B4GALT1AA , produced colostrum and milk with a significantly lower lactose content and reared piglets with a substantially higher body weight (by 0.30 kg) at 21 d of age.

Genome-Wide Association Analysis Reveals Genetic Loci and Candidate Genes for Chest, Abdominal, and Waist Circumferences in Two Duroc Pig Populations

Chest circumference (CC), abdominal circumference (AC), and waist circumference (WC) are regarded as important indicators for improving economic traits because they can reflect the growth and physiological status in pigs. However, the genetic architecture of CC, AC, and WC is still elusive. Here, we performed single-trait and multi-trait genome-wide association studies (GWASs) for CC, AC, and WC in 2,206 American origin Duroc (AOD) and 2,082 Canadian origin Duroc (COD) pigs. As a result, one novel quantitative trait locus (QTL) on Sus scrofa chromosome (SSC) one was associated with CC and AC in COD pigs, which spans 6.92 Mb (from 170.06 to 176.98 Mb). Moreover, multi-trait GWAS identified 21 significant SNPs associated with the three conformation traits, indicating the multi-trait GWAS is a powerful statistical approach that uncovers pleiotropic locus. Finally, the three candidate genes (ITGA11, TLE3, and GALC) were selected that may play a role in the conformation traits. Further bioinformatics analysis indicated that the candidate genes for the three conformation traits mainly participated in sphingolipid metabolism and lysosome pathways. For all we know, this study was the first GWAS for WC in pigs. In general, our findings further reveal the genetic architecture of CC, AC, and WC, which may offer a useful reference for improving the conformation traits in pigs.

Genome-wide scan reveals genomic regions and candidate genes underlying direct and maternal effects of preweaning calf mortality in Nellore cattle

We conducted analysis to estimate genetic parameters and to identify genomic regions and candidate genes affecting direct and maternal effects of preweaning calf mortality (PWM) in Nellore cattle. Phenotypic records of 67,196 animals, and 8443 genotypes for 410,936 SNPs were used. Analysis were performed through the weighted single-step GBLUP approach and considering a threshold animal model via Bayesian Inference. Direct and maternal heritability estimates were of 0.2143 ± 0.0348 and 0.0137 ± 0.0066, respectively. The top 10 genomic regions accounted for 13.61 and 14.23% of the direct and maternal additive genetic variances and harbored a total of 63 and 91 positional candidate genes, respectively. Two overlapping regions on BTA2 were identified for both direct and maternal effects. Candidate genes are involved in biological mechanisms i.e. embryogenesis, immune response, feto-maternal communication, circadian rhythm, hormone alterations, myometrium adaptation, and milk secretion, which are critical for the successful calf growth and survival during preweaning period.

Potential functional variants in AHR signaling pathways are associated with age at puberty in swine

Puberty in female pigs is defined as age at first estrus and gilts that have an earlier age at puberty are more likely to have greater lifetime productivity. Because age at puberty is predictive for sow longevity and lifetime productivity, but not routinely measured in commercial herds, it would be beneficial to use genomic or marker-assisted selection to improve these traits. A GWAS at the US Meat Animal Research Center (USMARC) identified several loci associated with age at puberty in pigs. Candidate genes in these regions were scanned for potential functional variants using sequence information from the USMARC swine population founder animals and public databases. In total, 135 variants (SNP and insertion/deletions) in 39 genes were genotyped in 1284 phenotyped animals from a validation population sired by Landrace and Yorkshire industry semen using the Agena MassArray system. Twelve variants in eight genes were associated with age at puberty (P < 0.005) with estimated additive SNP effects ranging from 1.6 to 5.3 days. Nine of these variants were non-synonymous coding changes in AHR, CYP1A2, OR2M4, SDCCAG8, TBC1D1 and ZNF608, two variants were deletions of one and four codons in aryl hydrocarbon receptor, AHR, and the most significant SNP was near an acceptor splice site in the acetyl-CoA carboxylase alpha, ACACA. Several of the loci identified have a physiological and a genetic role in sexual maturation in humans and other animals and are involved in AHR-mediated pathways. Further functional validation of these variants could identify causative mutations that influence age at puberty in gilts and possibly sow lifetime productivity.

[Clinical, ethologic, endokrinologic, and metabolic aspects of the peripartal period in pigs]

The aim of this review article is to provide an overview of the literature relevant to the peripartal period in swine. As in all other mammals, the farrowing process in pigs is divided into 3 phases (I-III; opening, expulsion, and postnatal stage), during which various essential endocrine and metabolic mechanisms initiate or maintain parturition. These include the hormones progesterone, cortisol, prostaglandin F_2α, oxytocin, estradiol, relaxin as well as electrolytes, enzymes, and metabolites such as calcium, magnesium, inorganic phosphate, glucose, creatine kinase, lactate, non-esterified free fatty acids, and β-hydroxybutyrate. Exogenous or endogenous disruptive factors may result in a delay or even stagnation of labor. For example, the form of husbandry may represent a possible exogenous disruptive factor. Endogenous disruptive factors may arise from insufficient storage and/or distribution of the above-mentioned labor-associated parameters. Subsequent dystocia leads to temporary or permanent consequences for maternal reproductive fitness and impairs piglet vitality at the time of birth, possibly resulting in lower survival rates.

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.

Assessing the Effect of SNPs on Litter Traits in Pigs

The reproductive ability of sows is the principle of continuous and efficient production, based on such traits as the number of piglets, the total number of parities, and the period of economic use. Currently, SNPs associated with the TNB and NBA are presented in the PigQTLdb. The aim of this work was the assessment of the SNP effects on the litter traits in Large White (LW, n = 502) and Landrace (LN, n = 432) sow breeds in a farm in Russia. 9 SNPs (SNP_1: rs80956812; SNP_2: rs81471381; SNP_3: rs80891106; SNP_4: rs81399474; SNP_5: rs81421148; SNP_6: rs81242222; SNP_7: rs81319839; SNP_8: rs81312912; SNP_9: rs80962240) were selected for the study. Associative analysis was performed using the GLM procedure in R version 3.5.1. The analysis of reproductive traits was carried out according to the results of the first parity, the second and subsequent parities, and totals for lifetime of sows. The significant effect on litter traits in LW was determined for SNP rs80956812, SNP rs81471381, SNP rs81421148, and SNP rs81399474. The significant effect on litter traits in LN was determined for SNP rs81421148 and SNP rs81319839. AKT3 gene was identified as perspective candidate gene, whose biological functions, as well as the results obtained in our work and in other studies, indicate its potential role in the reproductive process regulation in pigs. In general, the data obtained help to explain the genetic mechanisms of reproductive traits.

miR-205 Expression Elevated With EDS Treatment and Induced Leydig Cell Apoptosis by Targeting RAP2B via the PI3K/AKT Signaling Pathway

The adult Leydig cells (ALCs), originated from stem Leydig cells (SLCs), can secrete testosterone which is essential for germ cell development and sexual behavior maintenance. As a synthetic compound, ethane dimethane sulfonate (EDS), a well-known alkylating agent, has been reported to specifically ablate ALCs. In this study, EDS was verified to ablate differentiated pig LCs by experiments. Subsequently, the primary isolated pig LCs (containing SLCs and differentiated LCs) and EDS-treated LCs (almost exclusively SLCs) were collected for RNA-seq 4,904 genes and 15 miRNAs were differently expressed between the two groups. Down-regulated genes in the EDS-treated group were mainly related to steroid hormone biosynthesis. The highest up-regulation miRNAs was miR-205 after EDS treatment. Additionally, miR-205 was expressed more highly in pig SLCs clones compared with differentiated LCs. Through qRT-PCR, western blot (WB), TUNEL, EDU and flow cytometry, miR-205 was found to induce cell apoptosis, but did not affect proliferation or differentiation in both TM3 and GC-1spg mouse cell lines. Through luciferase reporter assays and WB, RAP2B was identified as a target gene of miR-205. Besides, overexpression of miR-205 inhibited the expressions of PI3K, Akt and p-AKT. All these findings were helpful for elucidating the regulation mechanism in pig LCs.

Survey of SNPs Associated with Total Number Born and Total Number Born Alive in Pig

Reproductive productivity depend on a complex set of characteristics. The number of piglets at birth (Total number born, Litter size, TNB) and the number of alive piglets at birth (Total number born alive, NBA) are the main indicators of the reproductive productivity of sows in pig breeding. Great hopes are pinned on GWAS (Genome-Wide Association Studies) to solve the problems associated with studying the genetic architecture of reproductive traits of pigs. This paper provides an overview of international studies on SNP (Single nucleotide polymorphism) associated with TNB and NBA in pigs presented in PigQTLdb as "Genome map association". Currently on the base of Genome map association results 306 SNPs associated with TNB (218 SNPs) and NBA (88 SNPs) have been identified and presented in the Pig QTLdb database. The results are based on research of pigs such as Large White, Yorkshire, Landrace, Berkshire, Duroc and Erhualian. The presented review shows that most SNPs found in chromosome areas where candidate genes or QTLs (Quantitative trait locus) have been identified. Further research in the given direction will allow to obtain new data that will become an impulse for creating breakthrough breeding technologies and increase the production efficiency in pig farming.

Genetic and genomic characterization of vulva size traits in Yorkshire and Landrace gilts

BACKGROUND

Reproductive performance is critical for efficient swine production. Recent results indicated that vulva size (VS) may be predictive of reproductive performance in sows. Study objectives were to estimate genetic parameters, identify genomic regions associated, and estimate genomic prediction accuracies (GPA) for VS traits.

RESULTS

Heritability estimates of VS traits, vulva area (VA), height (VH), and width (VW) measurements, were moderately to highly heritable in Yorkshire, with 0.46 ± 0.10, 0.55 ± 0.10, 0.31 ± 0.09, respectively, whereas these estimates were low to moderate in Landrace, with 0.16 ± 0.09, 0.24 ± 0.11, and 0.08 ± 0.06, respectively. Genetic correlations within VS traits were very high for both breeds, with the lowest of 0.67 ± 0.29 for VH and VW for Landrace. Genome-wide association studies (GWAS) for Landrace, reveled genomic region associated with VS traits on Sus scrofa chromosome (SSC) 2 (154-157 Mb), 7 (107-110 Mb), 8 (4-6 Mb), and 10 (8-19 Mb). For Yorkshire, genomic regions on SSC 1 (87-91 and 282-287 Mb) and 5 (67 Mb) were identified. All regions explained at least 3.4% of the genetic variance. Accuracies of genomic prediction were moderate in Landrace, ranging from 0.30 (VH) to 0.61 (VA), and lower for Yorkshire, with 0.07 (VW) to 0.11 (VH). Between-breed and multi-breed genomic prediction accuracies were low.

CONCLUSIONS

Our findings suggest that VS traits are heritable in Landrace and Yorkshire gilts. Genomic analyses show that major QTL control these traits, and they differ between breed. Genomic information can be used to increase genetic gains for these traits in gilts. Additional research must be done to validate the GWAS and genomic prediction results reported in our study.

Association between single-nucleotide polymorphisms within candidate genes and fertility in Landrace and Duroc pigs

Finding effective predictors of traits related to boar fertility is essential for increasing the efficiency of artificial insemination systems in pig breeding. The objective of this study was to find associations between single-nucleotide polymorphisms (SNPs) within candidate genes and fertility in the breeds Landrace and Duroc. Animals with breeding values for total number of piglets born, were re-sequenced for exonic regions of 14 candidate genes related to male and female fertility using samples from 16 Landrace boars and 16 Duroc boars (four with high and four with low breeding value of total number of piglets born for each breed for male fertility, and the same for female fertility) to detect genetic variants. Genotyping for the detected SNPs was done in 619 Landrace boars and 513 Duroc boars. Two SNPs in BMPR1 and one SNP in COX-2 were found significantly associated with the total number of piglets born in Landrace. In Duroc, two SNPs in PLC_z, one SNP in VWF and one SNP in ZP3 were found significantly associated with total number of piglets born. These SNPs explained between 0.27% and 1.18% of the genetic variance. These effects are too low for being used directly for selection purposes but can be of interest in SNP-panels used for genomic selection.

A genomewide association study in divergently selected lines in rabbits reveals novel genomic regions associated with litter size traits

Uterine capacity (UC), defined as the total number of kits from unilaterally ovariectomized does at birth, has a high genetic correlation with litter size. The aim of our research was to identify genomic regions associated with litter size traits through a genomewide association study using rabbits from a divergent selection experiment for UC. A high-density SNP array (200K) was used to genotype 181 does from a control population, high and low UC lines. Traits included total number born (TNB), number born alive (NBA), number born dead, ovulation rate (OR), implanted embryos (IE) and embryo, foetal and prenatal survivals at second parity. We implemented the Bayes B method and the associations were tested by Bayes factors and the percentage of genomic variance (GV) explained by windows. Different genomic regions associated with TNB, NBA, IE and OR were found. These regions explained 7.36%, 1.27%, 15.87% and 3.95% of GV, respectively. Two consecutive windows on chromosome 17 were associated with TNB, NBA and IE. This genomic region accounted for 6.32% of GV of TNB. In this region, we found the BMP4, PTDGR, PTGER2, STYX and CDKN3 candidate genes which presented functional annotations linked to some reproductive processes. Our findings suggest that a genomic region on chromosome 17 has an important effect on litter size traits. However, further analyses are needed to validate this region in other maternal rabbit lines.

Genome-wide association studies for the number of animals born alive and dead in duroc pigs

Litter size is one of the most important economic traits for pig production as it is directly related to the production efficiency. As an important litter size trait in pigs, the number of piglets born alive at birth (NBA) receives widespread interests in the pig industry. However, traits of piglets born dead, including the number of stillborn piglets (NS) and the piglets mummified at birth (NM) should be noted to explain the loss of reproduction. Herein, in the present study, a total of 803 producing sows were sampled and 2807 farrowing records for NBA, NM, and NS traits were collected in a Duroc swine population. Subsequently, a genome-wide association study (GWAS) was performed for NBA, NS and NM in parity groups 1 to 5. In total, 10 putative regions were found associated with these traits. After stepwise conditional analyses around the putative regions, eight independent signals were ultimately identified for NBA, NS, and NM, and there were seven promising candidate genes related to these traits, including ARID1A, RXRG, NFATC4, ABTB2, GRAMD1B, NDRG1, and APC. Our findings contribute to the understanding of the significant genetic causes of piglets born alive and dead, and could have a positive effect on pig production efficiency and economic profits.

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.

A genome-wide single nucleotide polymorphism and copy number variation analysis for number of piglets born alive

Background

In this study we integrated the CNV (copy number variation) and WssGWAS (weighted single-step approach for genome-wide association) analyses to increase the knowledge about number of piglets born alive, an economically important reproductive trait with significant impact on production efficiency of pigs.

Results

A total of 3892 samples were genotyped with the Porcine SNP80 BeadChip. After quality control, a total of 57,962 high-quality SNPs from 3520 Duroc pigs were retained. The PennCNV algorithm identified 46,118 CNVs, which were aggregated by overlapping in 425 CNV regions (CNVRs) ranging from 2.5 Kb to 9718.4 Kb and covering 197 Mb (~ 7.01%) of the pig autosomal genome. The WssGWAS identified 16 genomic regions explaining more than 1% of the additive genetic variance for number of piglets born alive. The overlap between CNVR and WssGWAS analyses identified common regions on SSC2 (4.2–5.2 Mb), SSC3 (3.9–4.9 Mb), SSC12 (56.6–57.6 Mb), and SSC17 (17.3–18.3 Mb). Those regions are known for harboring important causative variants for pig reproductive traits based on their crucial functions in fertilization, development of gametes and embryos. Functional analysis by the Panther software identified 13 gene ontology biological processes significantly represented in this study such as reproduction, developmental process, cellular component organization or biogenesis, and immune system process, which plays relevant roles in swine reproductive traits.

Conclusion

Our research helps to improve the understanding of the genetic architecture of number of piglets born alive, given that the combination of GWAS and CNV analyses allows for a more efficient identification of the genomic regions and biological processes associated with this trait in Duroc pigs. Pig breeding programs could potentially benefit from a more accurate discovery of important genomic regions.

Electronic supplementary material

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

Genome-wide associations and detection of potential candidate genes for direct genetic and maternal genetic effects influencing dairy cattle body weight at different ages

Background

Body weight (BW) at different ages are of increasing importance in dairy cattle breeding schemes, because of their strong correlation with energy efficiency traits, and their impact on cow health, longevity and farm economy. In total, 15,921 dairy cattle from 56 large-scale test-herds with BW records were genotyped for 45,613 single nucleotide polymorphisms (SNPs). This dataset was used for genome-wide association studies (GWAS), in order to localize potential candidate genes for direct and maternal genetic effects on BW recorded at birth (BW0), at 2 to 3 months of age (BW23), and at 13 to 14 months of age (BW1314).

Results

The first 20 principal components (PC) of the genomic relationship matrix (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathbf{G}}$$\end{document}G) grouped the genotyped cattle into three clusters. In the statistical models used for GWAS, correction for population structure was done by including polygenic effects with various genetic similarity matrices, such as the pedigree-based relationship matrix (\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathbf{A}}$$\end{document}A), the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathbf{G}}$$\end{document}G-matrix, the reduced \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathbf{G}}$$\end{document}G-matrix LOCO (i.e. exclusion of the chromosome on which the candidate SNP is located), and LOCO plus chromosome-wide PC. Inflation factors for direct genetic effects using \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathbf{A}}$$\end{document}A and LOCO were larger than 1.17. For \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\mathbf{G}}$$\end{document}G and LOCO plus chromosome-wide PC, inflation factors were very close to 1.0. According to Bonferroni correction, ten, two and seven significant SNPs were detected for the direct genetic effect on BW0, BW23, and BW1314, respectively. Seventy-six candidate genes contributed to direct genetic effects on BW with four involved in growth and developmental processes: FGF6, FGF23, TNNT3, and OMD. For maternal genetic effects on BW0, only three significant SNPs (according to Bonferroni correction), and four potential candidate genes, were identified. The most significant SNP on chromosome 19 explained only 0.14% of the maternal de-regressed proof variance for BW0.

Conclusions

For correction of population structure in GWAS, we suggest a statistical model that considers LOCO plus chromosome-wide PC. Regarding direct genetic effects, several SNPs had a significant effect on BW at different ages, and only two SNPs on chromosome 5 had a significant effect on all three BW traits. Thus, different potential candidate genes regulate BW at different ages. Maternal genetic effects followed an infinitesimal model.

Electronic supplementary material

The online version of this article (10.1186/s12711-018-0444-4) contains supplementary material, which is available to authorized users.

Genome-Wide Association Mapping of Host-Plant Resistance to Soybean Aphid

Soybean aphid [ Matsumura (Hemiptera: Aphididae)] is the most damaging insect pest of soybean [ (L.) Merr.] in the Upper Midwest of the United States and is primarily controlled by insecticides. Soybean aphid resistance (i.e., genes) has been documented in some soybean accessions but more sources of resistance are needed. Incorporation of the resistance into marketed varieties has also been slow. Genome-wide association mapping can aid in identifying resistant accessions by correlating phenotypic data with single nucleotide polymorphisms (SNPs) across a genome. Aphid population measures from 2366 soybean accessions were collected from published studies screening cultivated soybean () and wild soybean ( Siebold & Zucc.) with aphids exhibiting Biotype 1, 2, or 3 characteristics. Genotypic data were obtained from the SoySNP50K high-density genotyping array previously used to genotype the USDA Soybean Germplasm Collection. Significant associations between SNPs and soybean aphid counts were found on 18 of the 20 soybean chromosomes. Significant SNPs were found on chromosomes 7, 8, 13, and 16 with known genes. SNPs were also significant on chromosomes 1, 2, 4 to 6, 9 to 12, 14, and 17 to 20 where genes have not yet been mapped, suggesting that many genes remain to be discovered. These SNPs can be used to determine accessions that are likely to have novel aphid resistance traits of value for breeding programs.

Genetic effect and combined genotype effect of ESR, FSHβ, CTNNAL1 and miR-27a loci on litter size in a Large White population

To select new Large White line with high number of piglets born, genotypes of estrogen receptor (ESR), the follicle stimulating hormone β subunit (FSHβ), catenin alpha like 1 (CTNNAL1) and miR-27a were tested in 472 Large White sows. The associations of different genotypes with litter size traits were also studied. The results showed ESR^BB and FSHβ^BB sows produced 0.41-1.49 more pigs per litter (p < .05) for total number born (TNB) and number born alive (NBA) than did other corresponding genotypes. TNB of CTNNAL1^CG sows is 0.50 more pigs per litter (p < .05) than that of CTNNAL1^GG sows with the dominance effect of 0.25 pigs per litter (p < .05). miR-27a^BB sows had a less estimated breeding value (EBV) to TNB and had a more number of mummified pigs (NM) than did miR-27a^AA or miR-27a^AB sows (p < .05). Therefore, ESR^B, FSHβ^B, CTNNAL1^G, miR-27a^A allele was favorable for litter size traits. Furthermore, combined genetic effect analysis showed ESR^AAFSHβ^BB, ESR^AACTNNAL1^CG, ESR^AAmiR-27a^AA, FSHβ^BBCTNNAL1^CC, FSHβ^BBmiR-27a^AA and CTNNAL1^CG miR-27a^AB was the favorable combined genotype for litter size traits. These results identified favorable alleles and genotypes for litter size traits and suggested a potential selection scheme for litter size in Large White pigs.

A novel synonymous SNP (A47A) of the <i>TMEM95</i> gene is significantly associated with the reproductive traits related to testis in male piglets

Abstract. Transmembrane protein 95 (TMEM95) is located on the acrosomal membrane of the sperm head involved in the acrosome reaction; thus, it is regarded as affecting spermatogenesis and reproduction traits. The aim of this study was to explore the novel single nucleotide polymorphisms (SNPs) within the pig TMEM95 gene as well as to evaluate their associations with the testicular sizes in male Landrace (LD) and Large White (LW) breeds. After pool sequencing and bioinformatics analysis, only one novel coding SNP was found in exon 1, namely NC_010454.3: g.341T > C, resulting in a synonymous mutation (A47A). This SNP could be genotyped using the StuI polymerase chain reaction–restriction fragment length polymorphism (PCR-RFLP) assay. The minor allelic frequencies (MAFs) were 0.259 and 0.480 in the LD and LW breeds. Their polymorphism information content (PIC) values were 0.310 and 0.375. The LW population was at the Hardy–Weinberg equilibrium (HWE) (p > 0.05), whereas the LD population was not (p < 0.05). Association analyses demonstrated that a significant relationship was found between this A47A polymorphism and testis weight at 40 days of age in the LW population (p  =  0.047), and the heterozygote individuals showed lower testis weight than those with other genotypes. Moreover, this SNP was significantly associated with three testis measurement traits at 15 days of age in the LW population (p < 0.05); the individuals with genotypes TT and TC showed consistently superior testis measurement traits than those with genotype CC. These findings demonstrate that the A47A polymorphism had a significant effect on testis measurement traits, suggesting that the TMEM95 gene could be a candidate gene associated with reproductive traits. These results could contribute to breeding and genetics programs in the pig industry via DNA marker-assisted selection (MAS).

Identification of Stem Leydig Cells Derived from Pig Testicular Interstitium

Stem Leydig cells (SLCs), located in the testicular interstitial compartment in the mammalian testes, are capable of differentiating to testosterone-synthesizing Leydig cells (LCs), thus providing a new strategy for treating testosterone deficiency. However, no previous reports have identified and cultured SLCs derived from the pig. The aim of the current study was to isolate, identify, and culture SLCs from pigs. Haematoxylin and eosin staining and immunochemical analysis showed that SLCs were present and that PDGFRα was mainly expressed in the pig testicular interstitium, indicating that PDGFRα was a marker for SLCs in the neonatal pig. In addition, reverse transcription-PCR results showed that SLC markers were expressed in primary isolated LCs, indicating that they were putative SLCs. The putative SLCs were subsequently cultured with a testicular fluid of piglets (pTF) medium. Clones formed after 7 days and the cells expressed PDGFRα. However, no clones grew in the absence of pTF, but the cells expressed CYP17A1, indicating that pTF could sustain the features of porcine SLCs. To summarize, we isolated porcine SLCs and identified their basic characteristics. Taken together, these results may help lay the foundation for research in the clinical application of porcine SLCs.

Can multi-subpopulation reference sets improve the genomic predictive ability for pigs?

In most countries and for most livestock species, genomic evaluations are obtained from within-breed analyses. To achieve reliable breeding values, however, a sufficient reference sample size is essential. To increase this size, the use of multibreed reference populations for small populations is considered a suitable option in other species. Over decades, the separate breeding work of different pig breeding organizations in Germany has led to stratified subpopulations in the breed German Large White. Due to this fact and the limited number of Large White animals available in each organization, there was a pressing need for ascertaining if multi-subpopulation genomic prediction is superior compared with within-subpopulation prediction in pigs. Direct genomic breeding values were estimated with genomic BLUP for the trait "number of piglets born alive" using genotype data (Illumina Porcine 60K SNP BeadChip) from 2,053 German Large White animals from five different commercial pig breeding companies. To assess the prediction accuracy of within- and multi-subpopulation reference sets, a random 5-fold cross-validation with 20 replications was performed. The five subpopulations considered were only slightly differentiated from each other. However, the prediction accuracy of the multi-subpopulations approach was not better than that of the within-subpopulation evaluation, for which the predictive ability was already high. Reference sets composed of closely related multi-subpopulation sets performed better than sets of distantly related subpopulations but not better than the within-subpopulation approach. Despite the low differentiation of the five subpopulations, the genetic connectedness between these different subpopulations seems to be too small to improve the prediction accuracy by applying multi-subpopulation reference sets. Consequently, resources should be used for enlarging the reference population within subpopulation, for example, by adding genotyped females.

Detection of one large insertion/deletion (indel) and two novel SNPs within the <i>SPEF2</i> gene and their associations with male piglet reproduction traits

Abstract. The sperm flagella 2 (SPEF2) gene is essential for normal sperm tail development and male fertility. To fully characterize the structure of the mutation and to further study the function of the pig SPEF2 gene, we explored the insertion/deletion (indel) and novel single-nucleotide polymorphisms (SNPs) within the pig SPEF2 gene, and tested their associations with the testicular sizes in male Large White (LW) and Landrace (LD) pigs from China. Herein, a large insertion located at the SPEF2 gene in chromosome 16 was found, and two alleles of "I" (insertion) and "D" (deletion) were designated. Allele "D" was dominant in all analyzed pigs. Two novel SNPs (namely (NC_010458) g.19642G > A, resulting in AfaI aCRS PCR–PFLP, and g.19886C > G, resulting in EcoRI aCRS PCR–PFLP) were found in LW and LD pigs. Association testing revealed that g.19886C > G was significantly associated with the testis long circumference (TLC) in LW pigs (P < 0.05), suggesting that this SNP would be the DNA marker for the marker-assisted selection (MAS) in reproduction traits. This preliminary result indicates that the pig SPEF2 gene had significant effects on male reproduction traits. These findings could not only extend the spectrum of genetic variations in the pig SPEF2 gene but also contribute to implementing MAS in genetics and breeding in pigs.

Genetic determinants of pig birth weight variability

Background

Piglet birth weight variability, a trait also known as the within-litter homogeneity of birth weight, reflects the sow’s prolificacy, because it is positively genetically correlated with preweaning mortality but negatively correlated with the mean growth of piglets during sucking. In addition, the maternal additive genetic variance and heritability has been found exist for this trait, thus, reduction in the variability of piglet birth weight to improve the sow prolificacy is possible by selective breeding.

Results

We performed a genome wide association study (GWAS) in 82 sows with extreme standard deviation of birth weights within the first parity to identify significant SNPs, and finally 266 genome-wide significant SNPs (p < 0.01) were identified. These SNPs were mainly enriched on chromosome 7, 1, 13, 14, 15 and 18. We further scanned genes of the top 50 SNPs with the lowest p values and found some genes involved in plasma glucose homeostasis (GLP1R) and lipid metabolism as well as maternal-fetal lipid transport (AACS, APOB, OSBPL10 and LRP1B) which may contribute to the birth weight variability trait.

Conclusions

Birth weight variability trait has a low heritability. It is not easy to get significant signal by GWAS using small sample size. Herein, we identified some candidate chromosome regions especially chromosome 7 and suggested five genes which may provide some information for the further study.

Electronic supplementary material

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

Genome-wide association analyses for boar taint components and testicular traits revealed regions having pleiotropic effects

Background

The aim of this study was to perform a genome-wide association analyses (GWAS) for androstenone, skatole and indole in different Pietrain sire lines and compare the results with previous findings in purebred populations. Furthermore, the genetic relationship of androstenone and skatole were investigated with respect to pleiotropy. In order to characterize the performance of intact boars, crossbred progenies of 136 Pietrain boars mated to crossbred sows from three different breeding companies were tested on four test stations. A total of 598 boars were performance tested according to the rules of stationary performance testing in Germany. Beside common fattening and carcass composition traits, the concentrations of the boar taint components and testicular size parameters were recorded. All boars were genotyped with the PorcineSNP60 Illumina BeadChip. The GWAS were performed using the whole data set as well as in sub groups according to the line of origin. Besides an univariate GWAS approach, principal component (PC) techniques were applied to identify common expression pattern affecting the biosynthesis and the metabolism of androstenone.

Results

In total, 33 SNPs were significantly associated with at least one of the boar taint components. Only one SNP was identified being significant in both subgroups. The analyses of the testes size parameters revealed 31 significant associations. The numbers of significant SNPs within the genetic groups evidenced the strong population specific effects. A multivariate approach using PC revealed 33 significant associations for five different PC.

Conclusions

Based on Pietrain sired cross bred boars, the mayor objective of our study was to identify QTL for boar taint components and to detect pleiotropy among boar taint and testes traits. The high number of identified QTL revealed that boar taint traits are influenced by a large number of loci. Analyzing pleiotropy allowed identifying a QTL affecting androstenone and the gonasomatic index. In this region, QTL for ovulation rate and age at puberty of sows have been described in literature. This supports the physiological findings that the androstenone level of boars and reproduction performance of sows might be linked by an antagonistic relationship.

Electronic supplementary material

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