Mapping of fertility traits in Finnish Ayrshire by genome-wide association analysis

A whole-genome scan using single marker association was used to detect chromosome regions associated with seven female fertility traits in Finnish Ayrshire dairy cattle. The phenotypic data consisted of de-regressed estimated breeding values for 340 bulls which were estimated using a single trait model. Genotypes were obtained with the Illumina BovineSNP50 panel and a total of 35 630 informative, high-quality single nucleotide polymorphism (SNP) markers were used. The association analysis was performed using a mixed-model approach which fitted a fixed effect for each SNP and a random polygenic effect. We detected eleven genome-wide significant associations on eight different chromosomes. With at least chromosome-wise significance after Bonferroni correction, sixteen SNPs on nine chromosomes showed significant associations with one or more fertility traits. The results confirmed quantitative trait loci on three chromosomes (1, 2 and 20) for fertility traits previously reported for the same breed and one on chromosome four previously detected in Holstein cattle.     

Exploration of relationships between production and fertility traits in dairy cattle via association studies of SNPs within candidate genes derived by expression profiling

The objective of this work was to integrate findings from functional genomics studies with genome-wide association studies for fertility and production traits in dairy cattle. Association analyses of production and fertility traits with SNPs located within or close to 170 candidate genes derived from two gene expression studies and from the literature were performed. Data from 2294 Holstein bulls genotyped for 39557 SNPs were used. A total of 111 SNPs were located on chromosomal segments covered by a candidate gene. Allele substitution effects for each SNP were estimated using a mixed model with a fixed effect of marker and a random polygenic effect. Assumed covariance was derived either from marker or from pedigree information. Results from the analysis with the kinship matrix built from marker genotypes were more conservative than from the analysis with the pedigree-derived relationship matrix. From sixteen SNPs with significant effects on both classes of traits, ten provided evidence of an antagonistic relationship between productivity and fertility. However, we found four SNPs with favourable effects on fertility and on yield traits, one SNP with favourable effects on fertility and percentage traits, and one SNP with antagonistic effects on two fertility traits. While most quantitative genetic studies have proven genetic antagonisms between yield and functional traits, improvements in both production and functionality may be possible when focusing on a few relevant SNPs. Investigations combining input from quantitative genetics and functional genomics with association analysis may be applied for the identification of such SNPs.     

奶牛重要经济性状全基因组关联分析研究进展

产奶性状是奶牛最重要的生产性状,随着平衡育种理念的提出和发展,繁殖性状、体型性状、健康性状和长寿性等功能性状也逐渐被重视并纳入育种规划中。鉴定产奶性状和功能性状主效基因或遗传标记并将之应用于奶牛标记辅助选择可望加快遗传进展。随着高密度SNP标记的高通量检测技术的发展,全基因组关联分析已成为鉴定畜禽重要经济性状基因的重要途径。文章对奶牛产奶性状和功能性状全基因组关联分析研究进展进行综述。     

奶牛重要经济性状全基因组关联分析研究进展

产奶性状是奶牛最重要的生产性状, 随着平衡育种理念的提出和发展, 繁殖性状、体型性状、健康性状和长寿性等功能性状也逐渐被重视并纳入育种规划中。鉴定产奶性状和功能性状主效基因或遗传标记并将之应用于奶牛标记辅助选择可望加快遗传进展。随着高密度SNP标记的高通量检测技术的发展, 全基因组关联分析已成为鉴定畜禽重要经济性状基因的重要途径。文章对奶牛产奶性状和功能性状全基因组关联分析研究进展进行综述。     

Genome-wide associations for fertility traits in Holstein–Friesian dairy cows using data from experimental research herds in four European countries

Genome-wide association studies for difficult-to-measure traits are generally limited by the sample population size with accurate phenotypic data. The objective of this study was to utilise data on primiparous Holstein–Friesian cows from experimental farms in Ireland, the United Kingdom, the Netherlands and Sweden to identify genomic regions associated with traditional measures of fertility, as well as a fertility phenotype derived from milk progesterone profiles. Traditional fertility measures investigated were days to first heat, days to first service, pregnancy rate to first service, number of services and calving interval (CI); post-partum interval to the commencement of luteal activity (CLA) was derived using routine milk progesterone assays. Phenotypic and genotypic data on 37 590 single nucleotide polymorphisms (SNPs) were available for up to 1570 primiparous cows. Genetic parameters were estimated using linear animal models, and univariate and bivariate genome-wide association analyses were undertaken using Bayesian stochastic search variable selection performed using Gibbs sampling. Heritability estimates of the traditional fertility traits varied from 0.03 to 0.16; the heritability for CLA was 0.13. The posterior quantitative trait locus (QTL) probabilities, across the genome, for the traditional fertility measures were all <0.021. Posterior QTL probabilities of 0.060 and 0.045 were observed for CLA on SNPs each on chromosome 2 and chromosome 21, respectively, in the univariate analyses; these probabilities increased when CLA was included in the bivariate analyses with the traditional fertility traits. For example, in the bivariate analysis with CI, the posterior QTL probability of the two aforementioned SNPs were 0.662 and 0.123. Candidate genes in the vicinity of these SNPs are discussed. The results from this study suggest that the power of genome-wide association studies in cattle may be increased by sharing of data and also possibly by using physiological measures of the trait under investigation.     

Genome-wide association study for longevity in the Holstein cattle population

Longevity is one of the most important traits determining dairy cow profitability. In the last decades dairy cows suffered a lowering in the age at culling. With the aim to identify the genes involved in longevity, dates of birth, yields, dates of calving during lifespan and culling dates were collected for 946 culled cows which had been genotyped with the Bovine High Density panel. Using the GenABEL package in R, genome-wide association analysis was performed on three potential traits of longevity: (1) ‘days in production,’ (2) ‘days in herd,’ (3) number of calvings over lifespan.’ Five genome-wide significant single nucleotide polymorphisms (SNPs) associated with all three longevity traits were detected. Several consecutive SNPs identified on chromosomes 16 and 30 indicated the presence of two suggestive quantitative trait loci (QTL). The genes comprised in the QTL regions had biological functions related to fertility, reproductive disorders, heat stress and welfare of cows. These findings might contribute to improving breeding strategies to improve longevity.     

A genome-wide survey reveals a deletion polymorphism associated with resistance to gastrointestinal nematodes in Angus cattle

Gastrointestinal (GI) nematode infections are a worldwide threat to human health and animal production. In this study, we performed a genome-wide association study between copy number variations (CNVs) and resistance to GI nematodes in an Angus cattle population. Using a linear regression analysis, we identified one deletion CNV which reaches genome-wide significance after Bonferroni correction. With multiple mapped human olfactory receptor genes but no annotated bovine genes in the region, this significantly associated CNV displays high population frequencies (58.26 %) with a length of 104.8 kb on chr7. We further investigated the linkage disequilibrium (LD) relationships between this CNV and its nearby single nucleotide polymorphisms (SNPs) and genes. The underlining haplotype blocks contain immune-related genes such as ZNF496 and NLRP3. As this CNV co-segregates with linked SNPs and associated genes, we suspect that it could contribute to the detected variations in gene expression and thus differences in host parasite resistance.     

Non-synonymous mutations mapped to chromosome X associated with andrological and growth traits in beef cattle

Background

Previous genome-wide association analyses identified QTL regions in the X chromosome for percentage of normal sperm and scrotal circumference in Brahman and Tropical Composite cattle. These traits are important to be studied because they are indicators of male fertility and are correlated with female sexual precocity and reproductive longevity. The aim was to investigate candidate genes in these regions and to identify putative causative mutations that influence these traits. In addition, we tested the identified mutations for female fertility and growth traits.

Results

Using a combination of bioinformatics and molecular assay technology, twelve non-synonymous SNPs in eleven genes were genotyped in a cattle population. Three and nine SNPs explained more than 1% of the additive genetic variance for percentage of normal sperm and scrotal circumference, respectively. The SNPs that had a major influence in percentage of normal sperm were mapped to LOC100138021 and TAF7L genes; and in TEX11 and AR genes for scrotal circumference. One SNP in TEX11 was explained ~13% of the additive genetic variance for scrotal circumference at 12 months. The tested SNP were also associated with weight measurements, but not with female fertility traits.

Conclusions

The strong association of SNPs located in X chromosome genes with male fertility traits validates the QTL. The implicated genes became good candidates to be used for genetic evaluation, without detrimentally influencing female fertility traits.

Electronic supplementary material

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

The role of bovine causal genes underlying dairy traits in Spanish Churra sheep

In dairy cattle, quantitative trait nucleotides (QTNs) underlying quantitative trait loci (QTL) for milk production traits have been identified in bovine DGAT1, GHR and ABCG2 genes. The SPP1 gene has also been proposed to be a regulator of lactation. In sheep, QTL underlying milk production traits have been reported only recently, and no proven QTN has been identified. Taking into account the close phylogenetic relationship between sheep and cattle, this study examined the possible effects of the aforementioned genes on sheep milk production traits. We first studied the genetic variability of the DGAT1, GHR, ABCG2 and SPP1 genes in 15 rams of the Spanish Churra dairy sheep breed. Second, we performed an association analysis between SNPs identified in these genes and three milk production traits recorded in a commercial population of Churra sheep. This analysis revealed only three significant associations at the nominal level (P-value <0.05) involving allelic variants of the ABCG2 gene, whereas no significant association was found for the DGAT1, GHR and SPP1 genes. When the Bonferroni correction was applied to take into account the multiple tests performed, none of the associations identified at the nominal level remained significant. Nevertheless, taking into account the high level of false-negative findings that can arise when applying the stringent Bonferroni correction, we think that our results provide a valuable primary assessment of strong candidate genes for milk traits in sheep.     

A genome-wide association study of direct gestation length in US Holstein and Italian Brown populations

Direct gestation length influences economically important traits in dairy cattle that are related to birth and peri-natal survival of the calf. The objective of this study was to identify single nucleotide polymorphisms (SNPs) that are significantly associated with direct gestation length through a genome-wide association study. Data used in the analysis included 7,308,194 cow gestation lengths from daughters of 4743 United States Holstein sires in the Cooperative Dairy DNA Repository population and 580,157 gestation lengths from 749 sires in the Italian Brown population. Association analysis included 36,768 and 35,082 SNPs spanning all autosomes for Holstein and Brown Swiss, respectively. Multiple shrinkage Bayesian was employed. Estimates of heritability for both populations were moderate, with values of 0.32 (±0.03) and 0.29 (±0.02) for Holstein and Brown Swiss, respectively. A panel of SNPs was identified, which included SNPs that have significant effects on direct gestation length, of which the strongest candidate region is located on chromosome 18. Two regions not previously linked to direct calving ease and calf survival were identified on chromosome 7 and 28, corresponding to regions that contain genes related to embryonic development and foetal development. SNPs were also identified in regions that have been previously mapped for calving difficulty and longevity. This study identifies target regions for the investigation of direct foetal effects, which are a significant factor in determining the ease of calving.     

Genome‐wide association mapping for female fertility traits in Danish and Swedish Holstein cattle

A genome‐wide association study was conducted using a mixed model analysis for QTL for fertility traits in Danish and Swedish Holstein cattle. The analysis incorporated 2,531 progeny tested bulls, and a total of 36 387 SNP markers on 29 bovine autosomes were used. Eleven fertility traits were analyzed for SNP association. Furthermore, mixed model analysis was used for association analyses where a polygenic effect was fitted as a random effect, and genotypes at single SNPs were successively included as a fixed effect in the model. The Bonferroni correction for multiple testing was applied to adjust the significance threshold. Seventy‐four SNP‐trait combinations showed chromosome‐wide significance, and five of these were significant genome‐wide. Twenty‐four QTL regions on 14 chromosomes were detected. Strong evidence for the presence of QTL that affect fertility traits were observed on chromosomes 3, 5, 10, 13, 19, 20, and 24. The QTL intervals were generally smaller than those described in earlier linkage studies. The identification of fertility trait‐associated SNPs and mapping of the corresponding QTL in small chromosomal regions reported here will facilitate searches for candidate genes and candidate polymorphisms.     

Single nucleotide polymorphisms in the growth hormone and insulin-like growth factor-1 genes are associated with milk production, body condition score and fertility traits in dairy cows

The somatotrophic axis (GH-IGF) is a key regulator of animal growth and development, affecting performance traits that include milk production, growth rate, body composition, and fertility. The aim of this study was to quantify the association of previously identified SNPs in bovine growth hormone (GH1) and insulin-like growth factor 1 (IGF-1) genes with direct performance trait measurements of lactation and fertility in Holstein-Friesian lactating dairy cows. Sixteen SNPs in both IGF-1 and GH1 were genotyped across 610 cows and association analyses were carried out with traits of economic importance including calving interval, pregnancy rate to first service and 305-day milk production, using animal linear mixed models accounting for additive genetic effects. Two IGF-1 SNPs, IGF1i1 and IGF1i2, were significantly associated with body condition score at calving, while a single IGF-1 SNP, IGF1i3, was significantly associated with milk production, including milk yield (means ± SEM; 751.3 ± 262.0 kg), fat yield (21.3 ± 10.2 kg) and protein yield (16.5 ± 8.0 kg) per lactation. Only one GH1 SNP, GH33, was significantly associated with milk protein yield in the second lactation (allele substitution effect of 9.8 ± 5.0 kg). Several GH1 SNPs were significantly associated with fertility, including GH32, GH35 and GH38 with calving to third parity (22.4 ± 11.3 days) (GH32 and GH38 only), pregnancy rate to first service (0.1%) and overall pregnancy rate (0.05%). The results of this study demonstrate the effects of variants of the somatotrophic axis on milk production and fertility traits in commercial dairy cattle.     

Genome-wide association study identifies candidate genes for male fertility traits in humans

Despite the fact that hundreds of genes are known to affect fertility in animal models, relatively little is known about genes that influence natural fertility in humans. To broadly survey genes contributing to variation in male fertility, we conducted a genome-wide association study (GWAS) of two fertility traits (family size and birth rate) in 269 married men who are members of a founder population of European descent that proscribes contraception and has large family sizes. Associations between ～250,000 autosomal SNPs and the fertility traits were examined. A total of 41 SNPs with p ≤ 1 × 10(-4) for either trait were taken forward to a validation study of 123 ethnically diverse men from Chicago who had previously undergone semen analyses. Nine (22%) of the SNPs associated with reduced fertility in the GWAS were also associated with one or more of the ten measures of reduced sperm quantity and/or function, yielding 27 associations with p values < 0.05 and seven with p values < 0.01 in the validation study. On the basis of 5,000 permutations of our data, the probabilities of observing this many or more small p values were 0.0014 and 5.6 × 10(-4), respectively. Among the nine associated loci, outstanding candidates for male fertility genes include USP8, an essential deubiquitinating enzyme that has a role in acrosome assembly; UBD and EPSTI1, which have potential roles in innate immunity; and LRRC32, which encodes a latent transforming growth factor β (TGF-β) receptor on regulatory T cells. We suggest that mutations in these genes that are more severe may account for some of the unexplained infertility (or subfertility) in the general population.     

Genetic dissection of bull fertility in US Jersey dairy cattle

The service sire has been recognized as an important factor affecting herd fertility in dairy cattle. Recent studies suggest that genetic factors explain part of the difference in fertility among Holstein sires. The main objective of this study was to dissect the genetic architecture of sire fertility in US Jersey cattle. The dataset included 1.5 K Jersey bulls with sire conception rate (SCR) records and 96 K single nucleotide polymorphism (SNP) markers spanning the whole genome. The analysis included whole‐genome scans for both additive and non‐additive effects and subsequent functional enrichment analyses using KEGG Pathway, Gene Ontology (GO) and Medical Subject Headings (MeSH) databases. Ten genomic regions located on eight different chromosomes explained more than 0.5% of the additive genetic variance for SCR. These regions harbor genes, such as PKDREJ, EPB41L2, PDGFD, STX2, SLC25A20 and IP6K1, that are directly implicated in testis development and spermatogenesis, sperm motility and the acrosome reaction. In addition, the genomic scan for non‐additive effects identified two regions on BTA11 and BTA25 with marked recessive effects. These regions harbor three genes—FER1L5, CNNM4 and DNAH3—with known roles in sperm biology. Moreover, the gene‐set analysis revealed terms associated with calcium regulation and signaling, membrane fusion, sperm cell energy metabolism, GTPase activity and MAPK signaling. These gene sets are directly implicated in sperm physiology and male fertility. Overall, this integrative genomic study unravels genetic variants and pathways affecting Jersey bull fertility. These findings may contribute to the development of novel genomic strategies for improving sire fertility in Jersey cattle.     

Interactions of the bovine placental lactogen and prolactin receptor genes are associated with fertility traits in cattle

Decline in fertility of high-producing dairy cattle has become a global challenge to the dairy industry. Because of low heritability and complexity, it is difficult to find genetic markers for fertility traits in cattle. Here, we report the use of an in vitro fertilization (IVF) system and candidate gene approach to test genetic associations of a single-nucleotide polymorphism (SNP) in bovine placental lactogen (bPL), and its interactions with SNPs in the prolactin receptor (PRLR) and growth hormone receptor genes with fertility traits in an IVF system. The associations suggest a possible involvement of genetic interactions between bPL and PRLR in the fertilization and embryonic development processes, and the potential for application in a marker-assisted selection program.     

A Post-GWAS Replication Study Confirming the PTK2 Gene Associated with Milk Production Traits in Chinese Holstein

Our initial genome-wide association study (GWAS) demonstrated that two SNPs (ARS-BFGL-NGS-33248, UA-IFASA-9288) within the protein tyrosine kinase 2 (PTK2) gene were significantly associated with milk production traits in Chinese Holstein dairy cattle. To further validate if the statistical evidence provided in GWAS were true-positive findings, a replication study was performed herein through genotype-phenotype associations. The two tested SNPs were found to show significant associations with milk production traits, which confirmed the associations observed in the original study. Specifically, SNPs lying in the PTK2 gene were also detected by sequencing 14 unrelated sires in Chinese Holsteins and a total of thirty-three novel SNPs were identified. Thirteen out of these identified SNPs were genotyped and tested for association with milk production traits in an independent resource population. After Bonferroni correction for multiple testing, twelve SNPs were statistically significant for more than two milk production traits. Analyses of pairwise D’ measures of linkage disequilibrium (LD) between all SNPs were also explored. Two haplotype blocks were inferred and the association study at haplotype level revealed similar effects on milk production traits. In addition, the RNA expression analyses revealed that a non-synonymous coding SNP (g.4061098T>G) was involved in the regulation of gene expression. Thus the findings presented here provide strong evidence for associations of PTK2 variants with dairy production traits and may be applied in Chinese Holstein breeding program.     

A genome-wide association study of mare fertility in the Pura Raza Español horse

Despite the economic importance of fertility for the horse industry, few efforts have been made to achieve a better understanding of the genetic mechanisms underlying its control. This is probably due to the difficulty of obtaining reliable phenotypes and the complexity of modelling the environmental and management factors. This work is novel in that we propose to use reproductive efficiency (RE) as an indicator of mare fertility. To achieve this, we performed a genome-wide association study in the Pura Raza Español horse aimed at identifying genomic variants, regions, and candidate genes associated with fertility in mares. The dataset included 819 animals genotyped with the Affymetrix Axiom? Equine 670 K single-nucleotide polymorphisms (SNPs) Genotyping Array and the deregressed breeding values for RE trait, obtained using a ssBLUP model, employed as pseudo-phenotypic data. Our results showed 28 SNPs potentially associated with RE, which explained 87.19% of the genetic variance and 6.61% of the phenotypic variance. Those results were further validated in BayesB, showing a correlation between observed and predicted RE of 0.57. In addition, 15 candidate genes (HTRA3, SPIRE1, APOE, ERCC1, FOXA3, NECTIN-2, KLC3, RSPH6A, PDPK1, MEIOB, PAQR4, NM3, PKD1, PRSS21, IFT140) previously related to fertility in mammals were associated with the markers and genomic regions significantly associated with RE. To our knowledge, this is the first genome-wide association study performed on mare fertility.     

Genome-wide association mapping in Norwegian Red cattle identifies quantitative trait loci for fertility and milk production on BTA12

Reproductive performance is a critical trait in dairy cattle. Poor reproductive performance leads to prolonged calving intervals, higher culling rates and extra expenses related to multiple inseminations, veterinary treatments and replacements. Genetic gain for improved reproduction through traditional selection is often slow because of low heritability and negative correlations with production traits. Detection of DNA markers associated with improved reproductive performance through genome-wide association studies could lead to genetic gain that is more balanced between fertility and production. Norwegian Red cattle are well suited for such studies, as very large numbers of detailed reproduction records are available. We conducted a genome-wide association study for non-return rate, fertility treatments and retained placenta using almost 1 million records on these traits and 17 343 genome-wide single-nucleotide polymorphisms. Genotyping costs were minimized by genotyping the sires of the cows recorded and by using daughter averages as phenotypes. The genotyped sires were assigned to either a discovery or a validation population. Associations were only considered to be validated if they were significant in both groups. Strong associations were found and validated on chromosomes 1, 5, 8, 9, 11 and 12. Several of these were highly supported by findings in other studies. The most important result was an association for non-return rate in heifers in a region of BTA12 where several associations for milk production traits have previously been found. Subsequent fine-mapping verified the presence of a quantitative trait loci (QTL) having opposing effects on non-return rate and milk production at 18 Mb. The other reproduction QTL did not have pleiotropic effects on milk production, and these are therefore of considerable interest for use in marker-assisted selection.     

A Multiple-SNP Approach for Genome-Wide Association Study of Milk Production Traits in Chinese Holstein Cattle

The multiple-SNP analysis has been studied by many researchers, in which the effects of multiple SNPs are simultaneously estimated and tested in a multiple linear regression. The multiple-SNP association analysis usually has higher power and lower false-positive rate for detecting causative SNP(s) than single marker analysis (SMA). Several methods have been proposed to simultaneously estimate and test multiple SNP effects. In this research, a fast method called MEML (Mixed model based Expectation-Maximization Lasso algorithm) was developed for simultaneously estimate of multiple SNP effects. An improved Lasso prior was assigned to SNP effects which were estimated by searching the maximum joint posterior mode. The residual polygenic effect was included in the model to absorb many tiny SNP effects, which is treated as missing data in our EM algorithm. A series of simulation experiments were conducted to validate the proposed method, and the results showed that compared with SMMA, the new method can dramatically decrease the false-positive rate. The new method was also applied to the 50k SNP-panel dataset for genome-wide association study of milk production traits in Chinese Holstein cattle. Totally, 39 significant SNPs and their nearby 25 genes were found. The number of significant SNPs is remarkably fewer than that by SMMA which found 105 significant SNPs. Among 39 significant SNPs, 8 were also found by SMMA and several well-known QTLs or genes were confirmed again; furthermore, we also got some positional candidate gene with potential function of effecting milk production traits. These novel findings in our research should be valuable for further investigation.