Whole-genome sequencing identifies potential candidate genes for reproductive traits in pigs

Reproductive performance is a complex quantitative trait, that is determined by multiple genes, regulatory pathways and environmental factors. A list of major genes with large effect have been detected, although multiple QTLs are identified. To identify candidate genes for pig prolificacy, whole genome variants from five high- and five low-prolificacy Yorkshire sows were collected using whole-genome resequencing. A total of 13,955,609 SNPs and 2,666,366 indels were detected across the genome. Common differential SNPs and indels were identified between the two groups of sows. Genes encoding components of the TGF-beta signaling pathway were enriched with the variations, including BMP5, BMP6, BMP7, ACVR1, INHBA, ZFYVE9, TGFBR2, DCN, ID4, BAMBI, and ACVR2A. Several differential variants within these genes related to reproductive traits were identified to be associated with litter size. A comparison of selective regions and published QTL data suggests that NEDD9, SLC39A11, SNCA, and UNC5D are candidate genes for reproduction traits.     

QTL analysis of clinical-chemical traits in an F₂ intercross between Landrace and Korean native pigs

Clinical-chemical traits are essential when examining the health status of individuals. The aim of this study was to identify quantitative trait loci (QTL) and the associated positional candidate genes affecting clinical-chemical traits in a reciprocal F(2) intercross between Landrace and Korean native pigs. Following an overnight fast, 25 serum phenotypes related to clinical-chemical traits (e.g., hepatic function parameters, renal function parameters, electrolyte, lipids) were measured in >970 F(2) progeny. All experimental samples were subjected to genotyping analysis using 165 microsatellite markers located across the genome. We identified eleven genome-wide significant QTL in six chromosomal regions (SSC 2, 7, 8, 13, 14, and 15) and 59 suggestive QTL in 17 chromosomal regions (SSC 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 12, 13, 14, 15, 16, 17, and 18). We also observed significant effects of reciprocal crosses on some of the traits, which would seem to result from maternal effect, QTL on sex chromosomes, imprinted genes, or genetic difference in mitochondrial DNA. The role of genomic imprinting in clinical-chemical traits also was investigated. Genome-wide analysis revealed a significant evidence for an imprinted QTL in SSC4 affecting serum amylase levels. Additionally, a series of bivariate linkage analysis provided strong evidence that QTL in SSC 2, 13, 15, and 18 have a pleiotropic effect on clinical-chemical traits. In conclusion, our study detected both novel and previously reported QTL influencing clinical-chemical traits in pigs. The identified QTL together with the positional candidate genes identified here could play an important role in elucidating the genetic structure of clinical-chemical phenotype variation in humans and swine.     

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.     

猪的基因图谱及数量性状位点定位

在人类基因组计划的带动下,猪的遗传连锁图谱和细胞遗传学图谱有了较大的进步,利用目前猪基因组图谱的研究成果,通过基因组扫描法和候选基因法,可以对猪重要经济性状的主效基因位点进行区域定位,进而图位克隆,找到主效基因,为现代遗传育种奠定理论基础。     

Mapping quantitative trait loci affecting female reproductive traits on porcine chromosome 8

An understanding of the genetic control of porcine female reproductive performance would offer the opportunity to utilize natural variation and improve selective breeding programs through marker-assisted selection. The Chinese Meishan is one of the most prolific pig breeds known, farrowing three to five more viable piglets per litter than the European Large White breed. This difference in prolificacy is attributed to the Meishan's superior prenatal survival levels. The present study utilized a three-generation cross in which the founder grandparental animals were purebred Meishan and Large White pigs in a scan for quantitative trait loci (QTL) on porcine chromosome 8 (SSC8) associated with reproductive performance. Reproductive traits, including number of corpora lutea (ovulation rate), teat number, litter size, and prenatal survival, were recorded for as many as 220 F2 females. Putative QTL for the related traits of litter size and prenatal survival were identified at the distal end of the long arm of SSC8. A physiological candidate gene, SPP1, was found to lie within the 95% confidence interval of these QTL. A suggestive QTL for teat number was revealed on the short arm of SSC8. The present study demonstrates, to our knowledge, the first independent confirmation of QTL for fecundity on SSC8, and these QTL regions provide a crucial starting point in the search for the causal genetic variants.     

Polymorphisms on SSC15q21-q26 Containing QTL for reproduction in Swine and its association with litter size

Several quantitative trait loci (QTL) for important reproductive traits (ovulation rate) have been identified on the porcine chromosome 15 (SSC15). To assist in the selection of positional candidate swine genes for these QTL on SSC15, twenty-one genes had already been assigned to SSC15 in a previous study in our lab, by using the radiation hybrid panel IMpRH. Further polymorphism studies were carried out on these positional candidate genes with four breeds of pigs (Duroc, Erhualian, Dahuabai and Landrace) harboring significant differences in reproduction traits. A total of nineteen polymorphisms were found in 21 genes. Among these, seven in six genes were used for association studies, whereby NRP2 polymorphism was found to be significantly (p < 0.05) associated with litter-size traits. NRP2 might be a candidate gene for pig-litter size based on its chromosome location (Du et al., 2006), significant association with litter-size traits and relationships with Sema and the VEGF super families.     

QTL analysis of white blood cell, platelet and red blood cell-related traits in an F2 intercross between Landrace and Korean native pigs

Haematological traits play important roles in disease resistance and defence functions. The objective of this study was to locate quantitative trait loci (QTL) and the associated positional candidate genes influencing haematological traits in an F₂ intercross between Landrace and Korean native pigs. Eight blood‐related traits (six erythrocyte traits, one leucocyte trait and one platelet trait) were measured in 816 F₂ progeny. All experimental animals were genotyped with 173 informative microsatellite markers located throughout the pig genome. We report that nine chromosomes harboured QTL for the baseline blood parameters: genomic regions on SSC 1, 4, 5, 6, 8, 9, 11, 13 and 17. Eight of twenty identified QTL reached genome‐wide significance. In addition, we evaluated the KIT locus, an obvious candidate gene locus affecting variation in blood‐related traits. Using dense single nucleotide polymorphism marker data on SSC 8 and the marker‐assisted association test, the strong association of the KIT locus with blood phenotypes was confirmed. In conclusion, our study identified both previously reported and novel QTL affecting baseline haematological parameters in pigs. Additionally, the positional candidate genes identified here could play an important role in elucidating the genetic architecture of haematological phenotype variation in swine and in humans.     

Comparative mapping of bovine chromosome 27 with human chromosome 8 near a dairy form QTL in cattle

In the absence of a complete and annotated bovine genome sequence, detailed human-bovine comparative maps are one of the most effective tools for identification of positional candidate genes contributing to quantitative trait loci (QTL) in cattle. In the present study, eight genes from human chromosome 8 were selected for mapping in cattle to improve breakpoint resolution and confirm gene order on the comparative map near the 40 cM region of the BTA27 linkage map where a QTL affecting dairy form had previously been identified. The resulting map identified ADRB3 as a positional candidate gene for the QTL contributing to the dairy form trait based on its estimated position between 40 and 45 cM on the linkage map. It is also a functional candidate gene due to its role in fat metabolism, and polymorphisms in the ADRB3 gene associated with obesity and metabolic disease in humans, as well as, carcass fat in sheep. Further studies are underway to investigate the existence of polymorphisms in the bovine ADRB3 gene and their association with traits related to fat deposition in cattle.     

QTL modulating ear size and erectness in pigs

Ear size and erectness are important conformation measurements in pigs. An F(2) population established by crossing European Large White (small, erect ears) with Chinese Meishan (large, flop ears) was used to study the genetic influence of the two ear traits for the first time. A linkage map incorporating 152 markers on 18 autosomal chromosomes was utilised in a genome scan for QTL. Significant QTL were found on SSC1, 5, 7, 9 and 12 for the two traits. The QTL on SSC5 and SSC7 had major effects and were significant at the genome-wide level (P < 0.01). The QTL on SSC1 for ear erectness also had a major effect and was genome-wide significant (P < 0.01). The 95% confidence interval (CI) of the ear size QTL on SSC5 spanned only 4 cM. The QTL on SSC7 for the two ear traits each had a CI of <20 cM, and their positions overlapped with those of the major QTL affecting subcutaneous fat depths on the same chromosome. This study provides insights on the complex genetic influences underlying pig ear traits and will facilitate positional candidate gene analysis to identify causative DNA variants.     

A genome‐wide linkage analysis for reproductive traits in F2 Large White × Meishan cross gilts

Female reproductive performance traits in pigs have low heritabilities thus limiting improvement through traditional selective breeding programmes. However, there is substantial genetic variation found between pig breeds with the Chinese Meishan being one of the most prolific pig breeds known. In this study, three cohorts of Large White × Meishan F₂ cross‐bred pigs were analysed to identify quantitative trait loci (QTL) with effects on reproductive traits, including ovulation rate, teat number, litter size, total born alive and prenatal survival. A total of 307 individuals were genotyped for 174 genetic markers across the genome. The genome‐wide analysis of the trait‐recorded F₂ gilts in their first parity/litter revealed one QTL for teat number significant at the genome level and a total of 12 QTL, which are significant at the chromosome‐wide level, for: litter size (three QTL), total born alive (two QTL), ovulation rate (four QTL), prenatal survival (one QTL) and teat number (two QTL). Further support for eight of these QTL is provided by results from other studies. Four of these 12 QTL were mapped for the first time in this study: on SSC15 for ovulation rate and on SSC18 for teat number, ovulation rate and litter size.     

Genome‐wide association study to identify genomic regions affecting prolificacy in Lori‐Bakhtiari sheep

Several causative mutations in candidate genes affecting prolificacy have been detected in various sheep breeds. A genome‐wide association study was performed on estimated breeding values for litter size in Lori‐Bakhtiari sheep. Prolific ewes with twinning records and others with only singleton records were genotyped using the medium‐density Illumina Ovine SNP50 array. Four single nucleotide polymorphisms (SNPs) associated with litter size were identified on chromosomes 3, 6 and 22. The region on sheep chromosome 3 between 75 739 167 and 75 745 152 bp included two significant SNPs (s52383.1 and OAR3_80038014_X.1) in high linkage disequilibrium with each other. The region that surrounds these SNPs contains a novel putative candidate gene: luteinizing hormone/choriogonadotropin receptor (LHCGR), known to be involved in ovarian steroidogenesis and organism‐specific biosystem pathways in sheep. Known prolificacy genes BMPR1B, BMP15 and GDF9 were not associated with litter size in Lori‐Bakhtiari sheep, suggesting that other biological mechanisms could be responsible for the trait's variation in this breed.     

Indirect effect of IGF2 intron3 g.3072G>A mutation on prolificacy in sows

A QTL located in the paternally expressed insulin-like growth factor 2 (IGF2) gene is known to increase muscle growth and reduce fat deposition in pigs. This makes the QTL in IGF2 a good marker for use in pig breeding programmes. However, care has to be taken as it is postulated that increased leanness and lowered fat deposition may have a negative effect on the prolificacy and longevity of sows. Selection of sire and dam lines for different alleles of the mutation in the paternally imprinted IGF2 gene could actually provide a solution to this problem. Therefore, in this study, the effect of the IGF2 QTL on prolificacy-related traits in sows was investigated. It was found that the paternal IGF2 wild-type allele was associated with higher reproduction performance in the sow. Moreover, it was also examined whether the difference in prolificacy in sows could be a consequence of differential IGF2 expression in the ovarian follicles of the sow or whether it is mainly a secondary effect caused by differences in fatness traits. Therefore, IGF2 expression was measured in follicles of different sizes from sows with different genotypes for the paternal IGF2 allele. It was observed that, however, while the size of the follicles was associated with follicular IGF2 expression level, the IGF2 genotype was not. It could be concluded that the difference in prolificacy of sows with a different paternal IGF2 genotype could be a secondary effect, resulting from differences in fat deposition.     

Genotyping of Novel SNPs in BMPR1B,BMP15, and GDF9 Genes for Association with Prolificacy in Seven Indian Goat Breeds

Goats form the backbone of rural livelihood and financial security systems in India but their population is showing decreasing trend. Improvement of reproductive traits such as prolificacy offers a solution to stabilize the decreasing goat population and to meet the nutritional needs of growing human population. In the present study, six novel SNPs in three candidate genes for prolificacy (BMPR1B, BMP15, and GDF9) were genotyped in seven breeds of Indian goats to evaluate their association with litter size. Tetra primer ARMS-PCR and PCR-RFLP based protocols were developed for genotyping six novel SNPs, namely, T(-242)C in BMPR1B; G735A and C808G in BMP15; and C818T, A959C, and G1189A in GDF9 gene. The effect of breed was highly significant (p ≤ 0.01) on litter size but the effect of genotype was nonsignificant. The effect of parity on litter size was also significant in the prolific Black Bengal breed. The litter size differences observed between breeds are attributed to breed differences. Novel mutations observed at different loci in GDF9, BMP15, and BMPR1B genes do not contribute to the reproductive capability of the investigated breeds. Further studies with more number of breeds and animals exploring association of these novel SNPs with reproductive traits may be fruitful.     

SNP detection and genetic mapping of porcine genes encoding enzymes in hepatic metabolic pathways and evaluation of linkage with carcass traits

We have previously identified and mapped porcine expressed sequence tags (ESTs) derived from genes that are preferentially expressed in liver. The aim of the present study was to identify single nucleotide polymorphisms (SNPs) in porcine genes encoding enzymes in hepatic metabolic pathways and use the SNPs for mapping. Furthermore, these genes, which are involved in utilization and partitioning of nutrients, were examined for their effects on carcass and meat quality traits by linkage analyses. In total, 100 ESTs were screened for SNPs by single strand conformation polymorphism analyses across a diverse panel of animals with a 36% success rate. Twelve of 36 polymorphic loci segregated in a three-generation Duroc x Berlin Miniature Pig (F2) resource population, the DUMI resource population, and were genetically mapped. Interval mapping of the corresponding chromosomes was performed to verify mapping of the genes within quantitative trait loci (QTL) regions detected in this resource population. QTL with genome-wide significance were detected in the vicinity of GNMT, ESTL147 and HGD. These loci therefore are positional candidate genes.     

Quantitative trait loci analysis of phenotypic traits and principal components of maize tassel inflorescence architecture

Maize tassel inflorescence architecture is relevant to efficient production of F₁ seed and yield performance of F₁ hybrids. The objectives of this study were to identify genetic relationships among seven measured tassel inflorescence architecture traits and six calculated traits in a maize backcross population derived from two lines with differing tassel architectures, and identify Quantitative Trait Loci (QTL) involved in the inheritance of those tassel inflorescence architecture traits. A Principal Component (PC) analysis was performed to examine relationships among correlated traits. Traits with high loadings for PC1 were branch number and branch number density, for PC2 were spikelet density on central spike and primary branch, and for PC3 were lengths of tassel and central spike. We detected 45 QTL for individual architecture traits and eight QTL for the three PCs. For control of inflorescence architecture, important QTL were found in bins 7.02 and 9.02. The interval phi034—ramosa1 (ral) in bin 7.02 was associated with six individual architecture trait QTL and explained the largest amount of phenotypic variation (17.3%) for PC1. Interval bnlg344–phi027 in bin 9.02 explained the largest amount of phenotypic variation (14.6%) for PC2. Inflorescence architecture QTL were detected in regions with candidate genes fasciated ear2, thick tassel dwarf1, and ral. However, the vast majority of QTL mapped to regions without known candidate genes, indicating positional cloning efforts will be necessary to identify these genes.     

Impact of the ESR gene on litter size and production traits in Czech Large White pigs

To evaluate the effect of the PvuII polymorphism of the oestrogen receptor gene on litter size and production traits in Czech Large White swine, data from 1250 sows and 3600 litters were analysed with two four-trait animal models. The traits in the first model were number of piglets born alive in a sow's first litter, number of piglets born alive in second and subsequent litters, lifetime daily gain and lean meat percentage. The second model included number of piglets born, number of piglets born alive, number of piglets weaned and litter weight at weaning from first and subsequent litters. The oestrogen receptor (ESR) locus significantly affected prolicacy in the first parity and averaged over all parities (P < 0.05), with allele A superior to allele B. In the first parity, AA sows produced approximately 0.5 more live piglets per litter than BB sows. Averaged over all parities, this difference was c. 0.25 piglets. Results for total number of piglets born and number of piglets weaned were similar to results for numbers born alive. No significant dominance effect was found for prolificacy traits. For litter weight at weaning, no significant additive effect was observed at the ESR locus, but a significant negative dominance effect (-1.5 kg) was estimated averaged across parities (litters of AB sows were similar to litters of BB sows for this trait). No pleiotropic effect of the ESR polymorphism on average daily gain or lean meat percentage was found.