Genome-wide selective signatures mining the candidate genes for egg laying in goose

BACKGROUND

Improving the egg production of goose is a crucial goal of breeding, because genetics is the key factor affecting egg production. Thus, we sequenced the genomes of 55 Chinese indigenous geese from six breeds, which were divided into the high egg-laying group (ZE, HY, and SC) and low egg-laying group (ZD, LH, and ST). Based on the results of the inter-population selection signal analysis, we mined the selected genome regions in the high egg-laying germplasm population to identify the key candidate genes affecting the egg-laying traits.

RESULTS

According to the whole-genome sequencing data, the average sequencing depth reached 11.75X. The genetic relationships among those six goose breeds coincided with the breed's geographical location. The six selective signal detection results revealed that the most selected regions were located on Chr2 and Chr12. In total, 12,051 single-nucleotide polymorphism (SNP) sites were selected in all six methods. Using the enrichment results of candidate genes, we detected some pathways involved in cell differentiation, proliferation, and female gonadal development that may cause differences in egg production. Examples of these pathways were the PI3K-Akt signaling pathway (IGF2, COMP, and FGFR4), animal organ morphogenesis (IGF2 and CDX4), and female gonad development (TGFB2).

CONCLUSION

On analyzing the genetic background of six local goose breeds by using re-sequencing data, we found that the kinship was consistent with their geographic location. 107 egg-laying trait-associated candidate genes were mined through six selection signal analysis. Our study provides a critical reference for analyzing the molecular mechanism underlying differences in reproductive traits and molecular breeding of geese.

CYP24A1 is associated with fetal mummification in pigs

Mummified piglets are among the leading causes of fertility loss and severely hamper reproductive performance in pigs. However, the contributions of genomic variation to the emergence of mummified piglets (MUM) have rarely been studied. This study aims to (1) elucidate the genetic architecture of MUM in sows of parity 1 - 3 using a single-step genome-wide association study (ssGWAS). The ssGWAS involved genotyping-by-sequencing of Large White and Landrace pig breeds. (2) Explore the biological role of the candidate genes at the cellular level. A total of 185 and 48 genome-wide significant SNPs are associated with MUM in Large White and Landrace pigs, explaining 0.01-36.52% genetic variance for different significant loci, respectively. All the significant SNPs are parity-specific, and the numerous, consecutive significant loci likely generated the nine significant peaks in different parities. Multiple candidate genes (including CYP24A1, FBXO30, and ARHGEF28) are associated with fetal congenital and maternal diseases. Collectively, CYP24A1 regulation contributes to steady-state levels of embryo development genes. CYP24A1 is involved in reproduction and, immune and gestational disorders. Thus, it is associated with known newborn death traits and MUM in Large White sows. Altogether, these results improve the current understanding of the genetic architecture of MUM and expand the knowledge on genetic variations for selecting against mummified piglets in pig breeding.

Genome-Wide Association Study of Potential Meat Quality Trait Loci in Ducks

With continuously increasing living standards and health requirements of consumers, meat quality is becoming an important consideration while buying meat products. To date, no genome-wide association study (GWAS) for copy number variants (CNVs) and single nucleotide polymorphisms (SNPs) has been conducted to reveal the genetic effects on meat quality in ducks. This study analyzed the phenotypic correlation and heritability of fat, water, collagen, and protein content of duck breast muscle. To identify the candidate variants for meat quality, we performed a GWAS using 273 ducks from an F2 population. The results of the SNP GWAS showed that the BARHL2, COPS7B, and CCDC50 genes were associated with fat content; BLM, WDR76, and EOMES with water content; CAMTA1, FGD5, GRM7, and RAPGEF5 with collagen production; and RIMS2, HNRNPU, and SPTBN1 with protein content. Additionally, 3, 7, 1, and 3 CNVs were associated with fat, water, collagen, and protein content, respectively, in duck breast muscle. The genes identified in this study can serve as markers for meat quality. Furthermore, our findings may help devise effective breeding plans and selection strategies to improve meat quality.

Insulin-like growth factor 2 (IGF2) gene: molecular characterization, expression analysis and association of polymorphisms with egg production traits in goose (Anser cygnoides)

Insulin-like growth factor 2 (IGF2) belongs to the member of the insulin-like growth factors family, which plays key roles in animal growth, differentiation and proliferation, as well as reproduction and the regulation of ovarian follicle development. However, little is known about the goose IGF2 gene. In this study, a 1879 bp fragment that covered the complete coding region (CDS) of goose IGF2 cDNA was identified for the first time. The cDNA consists of an open reading frame of 574 nucleotides with the capacity to encode a prepro-IGF-II protein of 187 amino acids, which comprises a signal peptide (24 residues), IGF-II peptide (67 residues), and C-terminal peptide (96 residues), and is closely related to that of chicken. qPCR indicates that the goose IGF2 mRNA is differentially expressed in all examined tissues of fertilized eggs (28 days) and laying Zhedong White geese (270 days). Two novel single nucleotide polymorphisms (SNPs) were detected in exon 1 (G63A, Chr2: G26541617A) and intron 1 (G38A, Chr2: G26541479A) regions, and the synonymous mutation G63A showed a significant association with egg numbers (E180d) of Sanhua goose population (p < 0.05). All the information derived from this study could be valuable and facilitate further studies on the functions of goose IGF2 gene.

Associations of IGF2 and DRD2 polymorphisms with laying traits in Muscovy duck

Insulin-like growth factor 2 (IGF2) and dopamine receptor 2 (DRD2) play important roles in ovarian follicular development. In this study, we analyzed tissue-specific expression of the Muscovy duck IGF2 and DRD2 genes and cloned those genes transcripts. Polymorphisms in these genes were tightly linked with egg production traits and both genes were highly expressed in the ovary. Moreover, we identified five single nucleotide polymorphisms (SNPs) for IGF1 and 28 for DRD2. Mutations A-1864G and C-1704G of IGF2 were positively correlated with increased egg laying at 59 weeks (E59W) (P < 0.05). The C+7T and C+364G mutations of DRD2 were highly and significantly associated with first-egg age (FEA) and egg numbers at 300 days (E300D) (P < 0.01). Moreover, C+3301G and C+3545G of DRD2 were highly significantly associated with FEA, E59W and E300D (P < 0.01). Other mutations were positively associated with FEA or E300D or E59W (P < 0.05). These data suggest specific roles for IGF1 and DRD2 polymorphisms in egg production in Muscovy ducks.

Exploring transcriptomic diversity in muscle revealed that cellular signaling pathways mainly differentiate five Western porcine breeds

Background

Among transcriptomic studies, those comparing species or populations can increase our understanding of the impact of the evolutionary forces on the differentiation of populations. A particular situation is the one of short evolution time with breeds of a domesticated species that underwent strong selective pressures. In this study, the gene expression diversity across five pig breeds has been explored in muscle. Samples came from: 24 Duroc, 33 Landrace, 41 Large White dam line, 10 Large White sire line and 39 Piétrain. From these animals, 147 muscle samples obtained at slaughter were analyzed using the porcine Agilent 44 K v1 microarray.

Results

A total of 12,358 genes were identified as expressed in muscle after normalization and 1,703 genes were declared differential for at least one breed (FDR < 0.001). The functional analysis highlighted that gene expression diversity is mainly linked to cellular signaling pathways such as the PI3K (phosphoinositide 3-kinase) pathway. The PI3K pathway is known to be involved in the control of development of the skeletal muscle mass by affecting extracellular matrix - receptor interactions, regulation of actin cytoskeleton pathways and some metabolic functions. This study also highlighted 228 spots (171 unique genes) that differentiate the breeds from each other. A common subgroup of 15 genes selected by three statistical methods was able to differentiate Duroc, Large White and Piétrain breeds.

Conclusions

This study on transcriptomic differentiation across Western pig breeds highlighted a global picture: mainly signaling pathways were affected. This result is consistent with the selection objective of increasing muscle mass. These transcriptional changes may indicate selection pressure or simply breed differences which may be driven by human selection. Further work aiming at comparing genetic and transcriptomic diversities would further increase our understanding of the consequences of human impact on livestock species.

Electronic supplementary material

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

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

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

Epigenetic consequences of artificial reproductive technologies to the bovine imprinted genes SNRPN, H19/IGF2, and IGF2R

Animal breeders have made widespread use of assisted reproductive technologies to accelerate genetic improvement programs aimed at obtaining more, better and cheaper food products. Selection approaches have traditionally focused on Mendel’s laws of inheritance using parental phenotypic characteristics and quantitative genetics approaches to choose the best parents for the next generation, regardless of their gender. However, apart from contributing DNA sequence variants, male and female gametes carry parental-specific epigenetic marks that play key roles during pre- and post-natal development and growth of the offspring. We herein review the epigenetic anomalies that are associated with artificial reproductive technologies in current use in animal breeding programs. For instance, we demonstrate that bovine embryos and fetuses derived by in vitro culture and somatic cell nuclear transfer show epigenetic anomalies in the differentially methylated regions controlling the expression of some imprinted genes. Although these genomic imprinting errors are undetected in the somatic tissues after birth, further research is warranted to examine potential germ cell transmission of epimutations and the potential risks of reproducing cattle using artificial reproductive technologies.

Genetic analysis of reproductive traits and antibody response in a PRRS outbreak herd

Porcine reproductive and respiratory syndrome (PRRS) is the most economically significant disease impacting pig production in North America, Europe, and Asia, causing reproductive losses such as increased rates of stillbirth and mummified piglets. The objective of this study was to explore the genetic basis of host response to the PRRS virus (PRRSV) in a commercial multiplier sow herd before and after a PRRS outbreak, using antibody response and reproductive traits. Reproductive data comprising number born alive (NBA), number alive at 24 h (NA24), number stillborn (NSB), number born mummified (NBM), proportion born dead (PBD), number born dead (NBD), number weaned (NW), and number of mortalities through weaning (MW) of 5,227 litters from 1,967 purebred Landrace sows were used along with a pedigree comprising 2,995 pigs. The PRRS outbreak date was estimated from rolling averages of farrowing traits and was used to split the data into a pre-PRRS phase and a PRRS phase. All 641 sows in the herd during the outbreak were blood sampled 46 d after the estimated outbreak date and were tested for anti-PRRSV IgG using ELISA (sample-to-positive [S/P] ratio). Genetic parameters of traits were estimated separately for the pre-PRRS and PRRS phase data sets. Sows were genotyped using the PorcineSNP60 BeadChip, and genome-wide association studies (GWAS) were performed using method Bayes B. Heritability estimates for reproductive traits ranged from 0.01 (NBM) to 0.12 (NSB) and from 0.01 (MW) to 0.12 (NBD) for the pre-PRRS and PRRS phases, respectively. S/P ratio had heritability (0.45) and strong genetic correlations with most traits, ranging from -0.72 (NBM) to 0.73 (NBA). In the pre-PRRS phase, regions associated with NSB and PBD explained 1.6% and 3% of the genetic variance, respectively. In the PRRS phase, regions associated with NBD, NSB, and S/P ratio explained 0.8%, 11%, and 50.6% of the genetic variance, respectively. For S/P ratio, 2 regions on SSC 7 (SSC7) separated by 100 Mb explained 40% of the genetic variation, including a region encompassing the major histocompatibility complex, which explained 25% of the genetic variance. These results indicate a significant genomic component associated with PRRSV antibody response and NSB in this data set. Also, the high heritability and genetic correlation estimates for S/P ratio during the PRRS phase suggest that S/P ratio could be used as an indicator of the impact of PRRS on reproductive traits.

Endometrial gene expression profiling in pregnant Meishan and Yorkshire pigs on day 12 of gestation

Background

Litter size in pigs is a major factor affecting the profitability in the pig industry. The peri-implantation window in pigs is characterized by the coordinated interactions between the maternal uterine endometrium and the rapidly elongating conceptuses and represents a period of time during which a large percentage of the developing conceptuses are lost. However, the gene expression and regulatory networks in the endometrium contributing to the establishment of the maternal: placental interface remain poorly understood.

Results

We characterized the endometrial gene expression profile during the peri-implantation stage of development by comparing two breeds that demonstrate very different reproductive efficiencies. We employed the porcine Affymetrix GeneChip® to assay the transcriptomic profiles of genes expressed in the uterine endometrium obtained from Meishan and Yorkshire gilts (n = 4 for each breed) on day 12 of gestation (M12 and Y12, respectively). Total of 17,076 probesets were identified as "present" in at least two arrays. A mixed model-based statistical analysis predicted a total of 2,656 (q < 0.1) transcripts as differentially expressed between Meishan and Yorkshire pigs. Eighteen differentially expressed transcripts of interest were validated by quantitative real-time PCR. Gene ontology (GO) annotation revealed that the known functions of the differentially expressed genes were involved in a series of important biological processes relevant to early pregnancy establishment in the pig.

Conclusions

The results identified endometrial gene expression profiles of two breeds differing in litter size and identified candidate genes that are related to known physiological pathways related to reproductive prolificacy. These findings provide a deeper understanding of molecular pathways differing between two breeds at the critical peri-implantation stage of pregnancy, which can be utilized to better understand the events contributing to pregnancy establishment in the pig.

Lack of the associations of the polymorphisms in IGF2, MC4R and GNAS genes with reproduction traits in pigs and imprinting analysis of IGF2 gene in ovary and cornus uteri

In recent years, intensive attention has been put on improving reproductive performance of pigs. Several experiments aimed to identify markers associated with prolificacy, but this issue still remains open. In our study, we investigated associations between polymorphisms in IGF2, GNAS and MC4R genes with reproductive traits of Polish Landrace and Large White pigs. We did not find any significant associations for g. GNAS314T > 324C, IGF2 intron3-g.3072G > A or g. MC4R 1426G > A in Polish Landrace and Large White pigs. In the case of IGF2 intron3-g.3072G > A, this information is of great importance, because this marker is widely implemented in pigs breeding and previous experiments suggested its role in prolificacy of pigs. We also investigated expression of IGF2 gene and showed that this gene is monoallelically expressed in reproductive organs (ovary and cornus uteri).

The imprinted gene DIO3 is a candidate gene for litter size in pigs

Genomic imprinting is an important epigenetic phenomenon, which on the phenotypic level can be detected by the difference between the two heterozygote classes of a gene. Imprinted genes are important in both the development of the placenta and the embryo, and we hypothesized that imprinted genes might be involved in female fertility traits. We therefore performed an association study for imprinted genes related to female fertility traits in two commercial pig populations. For this purpose, 309 SNPs in fifteen evolutionary conserved imprinted regions were genotyped on 689 and 1050 pigs from the two pig populations. A single SNP association study was used to detect additive, dominant and imprinting effects related to four reproduction traits; total number of piglets born, the number of piglets born alive, the total weight of the piglets born and the total weight of the piglets born alive. Several SNPs showed significant (q-value < 0.10) additive and dominant effects and one SNP showed a significant imprinting effect. The SNP with a significant imprinting effect is closely linked to DIO3, a gene involved in thyroid metabolism. The imprinting effect of this SNP explained approximately 1.6% of the phenotypic variance, which corresponded to approximately 15.5% of the additive genetic variance. In the other population, the imprinting effect of this QTL was not significant (q-value > 0.10), but had a similar effect as in the first population. The results of this study indicate a possible association between the imprinted gene DIO3 and female fertility traits in pigs.

Single nucleotide polymorphisms in the imprinted bovine insulin-like growth factor 2 receptor gene (IGF2R) are associated with body size traits in Irish Holstein-Friesian cattle

The regulation of the bioavailability of insulin-like growth factors (IGFs) is critical for normal mammalian growth and development. The imprinted insulin-like growth factor 2 receptor gene (IGF2R) encodes a transmembrane protein receptor that acts to sequester and degrade excess circulating insulin-like growth factor 2 (IGF-II) - a potent foetal mitogen - and is considered an important inhibitor of growth. Consequently, IGF2R may serve as a candidate gene underlying important growth- and body-related quantitative traits in domestic mammalian livestock. In this study, we have quantified genotype-phenotype associations between three previously validated intronic bovine IGF2R single nucleotide polymorphisms (SNPs) (IGF2R:g.64614T>C, IGF2R:g.65037T>C and IGF2R:g.86262C>T) and a range of performance traits in 848 progeny-tested Irish Holstein-Friesian artificial insemination sires. Notably, all three polymorphisms analysed were associated (P ≤ 0.05) with at least one of a number of performance traits related to animal body size: angularity, body depth, chest width, rump width, and animal stature. In addition, the C-to-T transition at the IGF2R:g.65037T>C polymorphism was positively associated with cow carcass weight and angularity. Correction for multiple testing resulted in the retention of two genotype-phenotype associations (animal stature and rump width). None of the SNPs analysed were associated with any of the milk traits examined. Analysis of pairwise r(2) measures of linkage disequilibrium between all three assayed SNPs ranged between 0.41 and 0.79, suggesting that some of the observed SNP associations with performance may be independent. To our knowledge, this is one of the first studies demonstrating associations between IGF2R polymorphisms and growth- and body-related traits in cattle. These results also support the increasing body of evidence that imprinted genes harbour polymorphisms that contribute to heritable variation in phenotypic traits in domestic livestock species.

Association analysis between variants in bovine progesterone receptor gene and superovulation traits in Chinese Holstein cows

The objective of this study was to identify a predictor to forecast superovulation response on the basis of associations between superovulation performance and gene polymorphism. The PCR-RFLP method was applied to detect two reported single nucleotide polymorphisms (SNPs) of G59752C and T81637C (rs41614030) located in introns 3 and 4 of the bovine progesterone receptor (PGR) gene in 171 Chinese Holstein cows treated for superovulation and evaluate its associations with superovulation traits. In polymorphic locus 81637, all cows without superovulation response were g.81637TC and g.81637TT genotypes. Association analysis showed that these two SNPs had significant effects on the total number of ova (TNO) (p<0.05), and the T81637C polymorphism was significantly associated with the number of transferable embryos (p<0.05). In addition, significant additive effects (p<0.05) on TNO were detected in the polymorphisms of G59752C and T81637C. These results showed for the first time that the G59752C and T81637C polymorphisms in PGR gene were associated with superovulation traits and indicated that PGR gene can be used as a predictor for superovulation in Chinese Holstein cows.

DNA sequence polymorphisms within the bovine guanine nucleotide-binding protein Gs subunit alpha (Gsα)-encoding (GNAS) genomic imprinting domain are associated with performance traits

BACKGROUND

Genes which are epigenetically regulated via genomic imprinting can be potential targets for artificial selection during animal breeding. Indeed, imprinted loci have been shown to underlie some important quantitative traits in domestic mammals, most notably muscle mass and fat deposition. In this candidate gene study, we have identified novel associations between six validated single nucleotide polymorphisms (SNPs) spanning a 97.6 kb region within the bovine guanine nucleotide-binding protein Gs subunit alpha gene (GNAS) domain on bovine chromosome 13 and genetic merit for a range of performance traits in 848 progeny-tested Holstein-Friesian sires. The mammalian GNAS domain consists of a number of reciprocally-imprinted, alternatively-spliced genes which can play a major role in growth, development and disease in mice and humans. Based on the current annotation of the bovine GNAS domain, four of the SNPs analysed (rs43101491, rs43101493, rs43101485 and rs43101486) were located upstream of the GNAS gene, while one SNP (rs41694646) was located in the second intron of the GNAS gene. The final SNP (rs41694656) was located in the first exon of transcripts encoding the putative bovine neuroendocrine-specific protein NESP55, resulting in an aspartic acid-to-asparagine amino acid substitution at amino acid position 192.

RESULTS

SNP genotype-phenotype association analyses indicate that the single intronic GNAS SNP (rs41694646) is associated (P ≤ 0.05) with a range of performance traits including milk yield, milk protein yield, the content of fat and protein in milk, culled cow carcass weight and progeny carcass conformation, measures of animal body size, direct calving difficulty (i.e. difficulty in calving due to the size of the calf) and gestation length. Association (P ≤ 0.01) with direct calving difficulty (i.e. due to calf size) and maternal calving difficulty (i.e. due to the maternal pelvic width size) was also observed at the rs43101491 SNP. Following adjustment for multiple-testing, significant association (q ≤ 0.05) remained between the rs41694646 SNP and four traits (animal stature, body depth, direct calving difficulty and milk yield) only. Notably, the single SNP in the bovine NESP55 gene (rs41694656) was associated (P ≤ 0.01) with somatic cell count--an often-cited indicator of resistance to mastitis and overall health status of the mammary system--and previous studies have demonstrated that the chromosomal region to where the GNAS domain maps underlies an important quantitative trait locus for this trait. This association, however, was not significant after adjustment for multiple testing. The three remaining SNPs assayed were not associated with any of the performance traits analysed in this study. Analysis of all pairwise linkage disequilibrium (r2) values suggests that most allele substitution effects for the assayed SNPs observed are independent. Finally, the polymorphic coding SNP in the putative bovine NESP55 gene was used to test the imprinting status of this gene across a range of foetal bovine tissues.

CONCLUSIONS

Previous studies in other mammalian species have shown that DNA sequence variation within the imprinted GNAS gene cluster contributes to several physiological and metabolic disorders, including obesity in humans and mice. Similarly, the results presented here indicate an important role for the imprinted GNAS cluster in underlying complex performance traits in cattle such as animal growth, calving, fertility and health. These findings suggest that GNAS domain-associated polymorphisms may serve as important genetic markers for future livestock breeding programs and support previous studies that candidate imprinted loci may act as molecular targets for the genetic improvement of agricultural populations. In addition, we present new evidence that the bovine NESP55 gene is epigenetically regulated as a maternally expressed imprinted gene in placental and intestinal tissues from 8-10 week old bovine foetuses.

Single nucleotide polymorphisms at the imprinted bovine insulin-like growth factor 2 (IGF2) locus are associated with dairy performance in Irish Holstein-Friesian cattle

The imprinted insulin-like growth factor 2 gene (IGF2) encodes a fetal mitogenic hormone protein (IGF-II) and has previously been shown to be associated with performance in dairy cattle. In this study we assessed genotype-phenotype associations between four single nucleotide polymorphisms (SNPs) located within the bovine IGF2 locus on chromosome 29 and a range of performance traits related to milk production, animal growth and body size, fertility and progeny survival in 848 progeny-tested Irish Holstein-Friesian sires. Two of the four SNPs (rs42196909 and IGF2.g-3815A>G), which were in strong linkage disequilibrium (r2 = 0·995), were associated with milk yield (P ≤ 0·01) and milk protein yield (P ≤ 0·05); the rs42196901 SNP was also associated (P ≤ 0·05) with milk fat yield. Associations (P ≤ 0·05) with milk fat percentage and milk protein percentage were observed at the rs42196901 and IGF2.g-3815A>G SNPs, respectively. The rs42196909 and IGF2.g-3815A>G SNPs were also associated with progeny carcass conformation (P ≤ 0·05), while an association (P ≤ 0·01) with progeny carcass weight was observed at the rs42194733 SNP locus. None of the four SNPs were associated with body size, fertility and progeny survival. These findings support previous work which suggests that the IGF2 locus is an important biological regulator of milk production in dairy cattle and add to an accumulating body of research showing that imprinted genes influence many complex performance traits in cattle.