- Research article
- Open Access
Two-stage genome-wide association study identifies integrin beta 5 as having potential role in bull fertility
BMC Genomics volume 10, Article number: 176 (2009)
Fertility is one of the most critical factors controlling biological and financial performance of animal production systems and genetic improvement of lines. The objective of this study was to identify molecular defects in the sperm that are responsible for uncompensable fertility in Holstein bulls. We performed a comprehensive genome wide analysis of single nucleotide polymorphisms (SNP) for bull fertility followed by a second-stage replication in additional bulls for a restricted set of markers.
In the Phase I association study, we genotyped the genomic sperm DNA of 10 low-fertility and 10 high-fertility bulls using Bovine SNP Gene Chips containing approximately 10,000 random SNP markers. In these animals, 8,207 markers were found to be polymorphic, 97 of which were significantly associated with fertility (p < 0.01). In the Phase II study, we tested the four most significant SNP from the Phase I study in 101 low-fertility and 100 high-fertility bulls, with two SNPs (rs29024867 and rs41257187) significantly replicated. Rs29024867 corresponds to a nucleotide change of C → G 2,190 bp 3' of the collagen type I alpha 2 gene on chromosome 4, while the rs41257187 (C → T) is in the coding region of integrin beta 5 gene on chromosome 1. The SNP rs41257187 induces a synonymous (Proline → Proline), suggesting disequilibrium with the true causative locus (i), but we found that the incubation of bull spermatozoa with integrin beta 5 antibodies significantly decreased the ability to fertilize oocytes. Our findings suggest that the bovine sperm integrin beta 5 protein plays a role during fertilization and could serve as a positional or functional marker of bull fertility.
We have identified molecular markers associated with bull fertility and established that at least one of the genes harboring such variation has a role in fertility. The findings are important in understanding mechanisms of uncompensatory infertility in bulls, and in other male mammals. The findings set the stage for more hypothesis-driven research aimed at discovering the role of variation in the genome that affect fertility and that can be used to identify molecular mechanisms of development.
Fertilization is a critical event at the onset of mammalian development. The widespread use of artificial insemination has revealed great variation in fertility among sires . Some males display reduced fertility that can be overcome with higher semen volume for insemination, called compensable infertility, while others show an uncompensable infertility [2, 3].
Uncompensable infertility defects may result from molecular defects caused by abnormalities in spermatozoal DNA, RNA, or proteins, which impair the ability of spermatozoa to interact with oocytes and induce embryonic development [4–6]. The quality of nuclear vacuoles, DNA integrity, and chromatin structure have been proposed as potential causes of uncompensable fertility defects [7–10]. However, most causes of bull subfertility are still unknown and are likely multigenic.
Recent advances in animal genome sequencing and associated technologies are providing new insights into the genomics study of gametes and embryos [11–14]. For instance, high-throughput technologies, including massively parallel expression and protein quantification, have revealed numerous differences between the spermatozoa of subinfertile and fertile males [15–17]. These phenotypes reflect, among other things, the genetic differences among the various sires. Single nucleotide polymorphisms (SNPs) which represent the most abundant genomic variation, have proved useful in studies of genes associated with human diseases (e.g., cancer, stroke, and diabetes) [18–21] and economically important traits in livestock (e.g., horse, pig, and cattle) [12, 22–29]. The previous use of SNPs for fertility studies has been limited to a few markers, and their implication in male infertility has not yet been fully demonstrated [19, 30–33].
The objective of the present study was to use a high-throughput and a high-density SNP array to conduct a near-genome-wide association study of bull fertility. Spermatozoa DNA were isolated from well-characterized low fertility (n = 10) and high fertility (n = 10) bulls (Phase I study) and examined for approximately 10,000 SNPs, followed by the screening of the four most significant SNPs in a larger population (101 low- and 100 high-fertility bulls; Phase II study).
Pure Holstein bulls were selected based on their fertility as previously described by Peddinti et al. . Briefly, the progeny test program from Alta Genetics Inc. (Alta Advantage® program) involving approximately 180 farms milking an average of 850 cows each was used to select the bulls (Alta Genetics Inc; Calgary, Alberta, Canada). This program provides certain benefits, including DNA verification of the paternity of offspring and pregnancy diagnoses by veterinary palpation, instead of relying solely on non-return rates 60–90 days after breeding. This depth of data allows an accurate determination of both male and female fertility traits. In addition, this program provides unique research materials to identify molecular markers associated with fertility.
Definition of low- and high-fertility bulls
The fertility of bulls (Holstein), drawn from a total population of 874 bulls, is based on an average of 788 breeding, ranging from 101 to 11,997. We used the model described by Zwald et al. , which takes in account the breeding event, environmental factors and herd management factors that influence fertility performance of sires (i.e., effects of herd-year-month, parity, cow, days in milk, sire proven status). All these factors were adjusted using threshold models similar to the previously published models . The fertility prediction of each bull was obtained using the Probit.F90 software  and expressed as the percent deviation (Table 1) of its conception from the average conception of all bulls. For the QTL analysis, the fertility was converted to a Z-score.
Isolation of pure sperm cells and DNA extraction
Alta Genetics Inc. (Watertown, WI) provided frozen semen straws of selected bulls. Thawed spermatozoa were then purified through a Percoll gradient, washed, counted, and pelleted for DNA isolation . DNA was extracted from a pool of three different ejaculates of spermatozoa using the Puregene DNA isolation kit (Qiagen, Valencia, CA), with minor modifications. Spermatozoa were homogenized in the lysis buffer (containing 60 mM DTT and 60 μg proteinase K), incubated for 60 minutes at 55°C, and treated with RNase A (12 μg). Proteins were sedimented and DNA was subsequently precipitated using isopropanol. DNA was washed in ethanol, dissolved in TE buffer (pH 8.0), and quantified using the NanoDrop ND-1000 spectrophotometer (NanoDrop Technologies). DNA integrity was verified on an electrophoresis gel stained-agarose. DNA samples with high purity (A260/A280 ≥ 1.8) and no degradation were used for the Phase I and Phase II studies.
SNP genotyping (Phase I study)
The 10 K SNP Bovine Gene Chip (Affymetrix/ParAllele GeneChip; Affymetrix Inc., Santa Clara, CA) was used to genotype DNA samples (250 ng/μl) of 10 low-fertility and 10 high-fertility bulls. The experiment was carried out at Baylor College of Medicine (Houston, TX), and assays utilized molecular inversion probe (MIP) technology, allowing the multiplex detection of single base variants using a 4-color array hybridization assay [37, 38]. The hybridization, washing, staining, and chip scanning procedures were performed using the standard protocol recommended by the manufacturer (Affymetrix Inc., Santa Clara, CA) of 9,919 SNPs analyzed in all DNA samples, the pass and completeness rates were 94.04% and 98.2%, respectively.
Allelic discrimination analysis (Phase II study)
The allelic discrimination technique, based on the TaqMan technology (ABI Prism 7000 Sequence Detection System, Applied Biosystems, Foster City, CA) was performed by scientists at SeqWright, Inc. seqwright.com on DNA samples (4 ng/μl) of 101 low-fertility and 100 high-fertility bulls to validate the Phase I association study. Primer and probe sets were designed on a sequence of 501 nucleotides containing the SNP (250 nucleotides downstream and upstream of the SNP allele) using the Primer Express software (Applied Biosystems, Foster City, CA). The probes were designed and labeled with FAM (6-carboxy-fluorescein) or VIC fluorescent dyes (Applied Biosystems, Foster City, CA) to match perfectly either one of the alleles (Table 1).
PCR reactions were carried out in a total volume of 25 μl as recommended by the manufacturer. Each reaction consisted of 5 μl of sperm DNA in 20 μl Master Mix solution (Applied Biosystems, TaqMan Universal PCR Master Mix) containing primers (900 nM) and probes (200 nM). DNA samples were amplified by 40 times (2 min-50°C, 10 min-95°C, 15 sec-92°C, and 1 min 60°C). For each PCR run, negative (no-template) and positive (oligos) controls were added, and each run was preceded or followed by 1-minute incubation at 60°C to determine the background, or final levels of fluorescence, in each reaction. The base calls were made by examining all samples on the allelic discrimination graph, and the Sequence Detection Software (Applied Biosystems, Foster City, CA) was used to determine the homozygosis or heterozygosis of alleles. All samples were run in triplicates for each SNP allele.
The major repository of SNP data in the National Center for Biotechnology Information (NCBI) database (dbSNP), combined with Entrez Genome (NCBI) as well as the Ensembl automatic analysis pipeline, were used to characterize and annotate the markers (SNPs and genes) based on the cattle genome assembly version 3.1. For each SNP, a sequence of 501 nucleotides consisting of the 250 nucleotides immediately upstream and downstream of the SNP allele, were extracted and the best hit in the Bovine 3.1 genome identified using BLASTN. We considered only the hits (i) found on the same chromosome as the query, (ii) with E-values equal or close to 0.00 and (iii) presenting at least 94% alignment with the full length of the query sequence. Furthermore, the same query sequences were used to search for putative transcription factor binding sites through the Transfac-blastX (TRANSFAC Database 7.0 for searching eukaryotic transcription factors).
Functional analysis of integrin beta 5 (ITGB5)
Sperm preparation and treatment
Thawed spermatozoa were purified through a Percoll gradient as previously reported . Motile spermatozoa were resuspended in the fertilization medium containing Heparin and PHE. The sperm concentration was adjusted to 50 × 106/ml and incubated with or without the integrin beta 5 antibody (ITGB5; sc-5401: 5 and 20 μg/ml). As a control for structural effects of the antibody, spermatozoa (50 × 106/ml) were also incubated with a nonmammalian protein (BIT) antibody (sc-33757: 20 μg/ml). After a 2-hour incubation at 38°C under 5% CO2 in air, spermatozoa were washed twice in fertilization medium and used to fertilize the oocytes. The motility of spermatozoa was comparable before and after the incubation period.
Oocyte maturation and in vitro fertilization
Bovine oocytes were purchased from Bomed Inc. (Madison, WI), and maturation took place during the transportation. These oocytes were washed and fertilized with treated and untreated spermatozoa at a final concentration of 106 spermatozoa/ml. After 18 hours of co-incubation, oocytes were collected, denuded, washed, fixed, and placed on slides for staining with 0.1% aceto-orcein on slides. The nuclear status of oocytes was observed under a microscope, and the fertilized oocytes (two extruded polar bodies and/or two pronuclei) were counted to assess the fertilization rates.
In the Phase I study, the 20 bulls of varied fertility were typed for 9,919 SNP using Affymetrix 10 K Xba 142 2.0 array. A total of 1,712 markers fixed for a single genotype in this sample were excluded from further analysis. Markers were tested for Hardy-Weinberg Equilibrium using a 1 degree of freedom chi-square test. Markers with a p value < 0.05 were flagged as potentially out of HWE, but were used for association analysis.
The fertility of the bulls was converted to Z-scores and used as the quantitative trait for the association analysis. Analysis was conducted only at the markers' positions; no attempt at interval mapping was made. At each marker, the mean fertility, scored as a Z-score, was compared between genotypes using single marker regression. If the n for a genotype was 1, the samples were collapsed into the heterozygote. Benjamini- Hochberg FDR was calculated at each marker, adjusting 8,207 tests. A p-value less than 0.01 was set as the threshold of a significant association between the SNP marker and bull fertility; this corresponds to an FDR of 84%. Similar analyses were conducted in Phase II data.
Results and discussion
The lack of methods to accurately predict sire fertility obliges the artificial insemination (AI) industry to keep and test hundreds of bulls. The selection of such bulls, whose fertility has been evaluated with progeny tests prior to their use in large-scale breeding programs, is costly and can take several years. Recent advances in cattle genome projects and molecular genetic technologies have increased the likelihood of identifying uncompensatory defects impairing the functions of spermatozoa. High-throughput technologies may help pinpoint relationships between a single DNA marker (i.e., SNP) and economically important traits in dairy cattle. Information about specific SNPs can enhance the efficiency of genetic selection, especially for traits that are (i) too difficult or expensive to measure in all animals and/or (ii) expressed after a long field trial (i.e., fertility).
Determination of bulls with uncompensatory infertility
The most fertile and the least fertile pure Holstein bulls were selected from a pool of 874 available bulls with at least 300 breeding records. These were defined as the representative outliers for their corresponding groups and the scaling of fertility as the deviation from the population average fertility. The average breeding and fertility rates of the 10 low- and 10 high-fertility bulls used in the Phase I association study were 1,974 and -9.2% versus 3,540 and 6.2%, respectively. In the Phase II study, these values were 904 and -4.2% for low-fertility (100), and 994 and 3.4%, for high (101) fertility bulls (Table 2). The average differences in fertility rates between low and high fertile groups were 15.4% and 7.6% in Phase I and Phase II studies, respectively. These differences were considered significant by Amann and Hammersted, 2002 .
Identification of SNPs associated with fertility using DNA microarrays (Phase I study)
The recent progress in genomics and automation has rendered the SNP genotyping a promising technology for genetic studies of the cow genome, which contains approximately one SNP every 252 base pairs . Using bovine SNP genechip microarrays (Affymetrix Inc.), we successfully genotyped 9,919 SNP markers using bull spermatozoa genomic DNA isolated from 10 low-fertility and 10 high-fertility Holstein bulls (Figure 1). Approximately 50% (4,963) of SNPs were physically located to chromosomes using the bovine genome assembly Btau-3.1 (National Center for Biotechnology Information, 2007).
No significant difference was found between the average of SNP call rates in both groups (95.8% ± 4.8% versus 98.6% ± 0.7%, respectively; p-value = 0.15). The HWE test revealed an FDR of 61% corresponding to 790 SNP markers that failed the test (versus the 496 expected at random), but were used for QTL analysis as mentioned above (Statistical analysis section). A total of 1,712 SNP markers were excluded from analysis because they were not polymorphic. In the Phase I samples, the fertility (Z-score) of bulls was used as the quantitative trait for the association analysis of the remaining 8,207 markers. The p-value distribution of the 8,207 markers' test showed a slight deviation from the uniform to left end (near 0) of the distribution, indicating the presence of more putative fertility-associated SNP markers in the data (Figure 2) than expected at random. Additional file 1 shows the Z-score data obtained from the four SNP markers with the highest association (p ≤ 10-4), and Additional file 2 shows the 97 significant markers p < 10-2. These most significant markers were selected to test in the Phase II samples.
High-throughput SNP arrays are powerful tools commonly used in humans to track population history and genes associated with diseases [19–21, 42, 43] or fertility [32, 44]. The application of a such high-throughput technology has been limited in livestock [24, 28, 29]. The 10 K SNP Bovine Gene Chip commercialized by Affymetrix Inc. ensures a limited coverage of the cow genome (≈ 2.8 Gb, based on the Btau 3.1 genome assembly), with one SNP every 300 to 350 Kb. However, the SNP location on the array seems evenly distributed across the genome (1.18 to 2.24 SNP/cM) .
The extent of linkage disequilibrium (LD) around the SNPs in the Holstein population is not clear yet, but we believe the LD may extend a good distance (≤ 0.5 Mb) around the markers due to the extensive inbreeding in the population [46–48]. For example, one specific bull not in this study was found to be an ancestor 163 times, often several times in a single bull's ancestors, within 5 generations of the bulls in our study. As a result, there may be some founding effects within this study and large stretches of LD, but it is doubtful we were able to fully cover the genome as close to the maximal Fisher's Information given the marker density and coupled with the low power of this study. This suggests there are additional fertility loci to be found in the genome. The need of higher density SNP arrays (50,000 to 300,000 SNPs) has been suggested for power association and fine mapping studies in cattle [41, 47, 49].
Genotyping of large numbers of bulls for the most significant SNP markers (Phase II study)
Large population sizes, often too large to be easily collected, are generally required for defining association between a given SNP and a trait with reasonable power. Thus, the use of several stages of analysis with progressively fewer markers typed in each replication population has become typical in human genome-wide association studies. The key reliability of a significant association result is replication in a second (or third) population. For example, a genome-wide-association (GWA) study of macular degeneration was successfully replicated when the initial study contained only 96 cases and 50 controls . Here, we selected the four most significant SNPs from the genotyping (p ≤ 10-4; Phase I) study (Table 3) for replication in a larger cohort of 100 low- and 101 high-fertility bulls (Figure 3). Interestingly, the marker rs29015574 out of HWE in the Phase I study became in HWE in Phase II, as well as three other markers (HWE test > 0.05; Table 4). The fertility (quantitative trait) analysis confirmed the significant association of two of the four SNPs (rs29024867 and rs41257187; P < 0.05), while the rs29016875 tended to be significant (P = 0.09). Our calculated FDR rate for these markers (26%–63%) in the Phase I data is very close to our observed replication rate of 50%, which is known as very high for a genome wide association study [51, 52].
Bioinformatics on the four highly associated SNPs
Bioinformatics were restricted on the four highest significantly associated SNPs (rs29016875, rs29015574, rs29024867, and rs41257187). Only three SNPs gave valuable information (Table 5). The SNP rs29024867 appeared as a positional marker for the collagen I alpha 2 gene, located at 2,190 base pairs of the 5' side of the SNP. The SNP rs29015574 had a single hit corresponding to a sequence similar to Rab3A-interacting molecule, which might be a potential gene candidate for fertility. The SNP rs41257187 is a synonymous polymorphism located in the exon 16 of integrin beta 5 gene (ITGB5). Each of these SNPs was located in a region (≤ 1 Mb) with low density (1–3 SNPs). Using similar SNP genechip array, Daetwyler et al., have reported a SNP frequency (SNP/cM) of 1.55, 1.75, 1.84 and 2.05 on Bos taurus autosome 9, 1, 4 and 10, respectively (Confer Table 3 for corresponding SNPs) . This observation makes the four SNP as strong positional candidate gene markers.
Putative binding sites for transcription factors with roles in embryogenesis were found within the nucleotide sequences flanking the SNPs rs29015574 (MEF-2D and Irx for, myocyte enhancer factor-2D and Iroquois homeobox, respectively) and rs29024867 (RAR β3 and Ceh-24 for, retinoic acid receptor and Caenorhabditis elegans-24, respectively). Hence, our findings suggest potential roles of the SNP mutations on regulation of gene expressions as already suggested in other studies [53–56].
Taken together, our results provide novel loci candidates whose associations with the bull fertility have not previously been reported. Because of the great interest of the non-synonymous SNP rs41257187 associated with the ITGB5, we performed a Six Frame Translation of ITGB5 mRNA reference sequence reported on NCBI (NM_174679.2) using the Baylor College of Medicine's HGSC Search Launcher to find the reading frame leading to the amino acid change. We found that the complete ITGB5 protein could be obtained from two different reading frames (+3 or +1) when using the whole length of RNA sequence or the length from the start codon (ATG), respectively. Contrary to the NCBI report, we found that the SNP mutation site (C/T) induced a synonymous amino acid change of a Proline (CCC) to another Proline (CCT) at the position 778 (P778P). The same SNP (rs41257187) inducing a synonymous amino acid change (Proline to Proline) is also reported on the Ensembl database.
The interest of integrin beta 5 in our study resided in the role(s) played by the integrin family members during fertilization and embryogenesis. Integrins are known to be expressed in a variety of tissues, including reproductive tissues of mice, humans, pigs and cattle [57–62].
Functional analysis of integrin beta 5
The biological function of the integrin beta 5 subunit is not clearly defined. We tested the potential involvement of ITGB5 isoform on sperm-egg interaction, and the results showed that the percentages of matured oocytes (metaphase II) undergoing a successful fertilization (1–2 PN) with spermatozoa pre-exposed to anti-integrin beta 5 antibody were decreased: 74, 60, and 47% normally fertilized oocytes in the presence of 0, 5, and 20 μg ITGB5/ml, respectively. This dose-dependent inhibition was not attributed to the antibody itself since the pre-exposure of spermatozoa to the non-mammalian protein (BIT) antibody had no effect on the ability of spermatozoa to fertilize the oocytes (Table 6).
We ruled out the likelihood of non-specific interactions of the anti-ITGB5 antibody and other integrin beta subunits. Indeed, the antibody is raised against a specific N-terminal epitope, which amino acid sequence is not found in other subunits of integrin beta protein. Therefore, our results contribute to the growing body of reports supporting the presence of αvβ integrins on sperm membranes [62, 63]. Their putative differential expression in subfertile and fertile bulls could be used as markers for fertility, as already suggested for the αv6β3 integrin in humans . Furthermore, our results contrast with the main body of literature reporting the presence of integrins mainly on the oocyte membrane while their ligands, ADAM family proteins, are on the sperm membrane [57, 58, 60, 61, 65]. In addition to this structural receptor function of integrins at fertilization (participation in cell-cell and cell-matrix interactions), integrin beta 5 might serve as a signaling receptor that induces serial events (such as inositol lipid turnover and protein phosphorylation) in the sperm, affecting fertilization and early embryo development.
We have employed a high-density SNP genome association study to identify loci that may play a role in fertility in dairy cattle and significant results in two of four loci were replicated in a second population. Functional studies, including one for genes harboring a replicated locus, the ITGB5 gene, suggest that it may play a role in sperm-egg interaction. These results provide a foundation for more hypothesis-driven research in genome biology, evolution, QTL mapping, and discovering genes and genomic regions for selectable traits.
Killian G: High-fertility proteins enhance reproduction rates in Dairy cattle. NRI Research Highlights, Competitive grants program, Factsheet, USDA. 1999
Den Daas JHG, De Jong G, Lansbergen LMTE, Van Wagtendonk-De Leeuw AM: The Relationship Between the Number of Spermatozoa Inseminated and the Reproductive Efficiency of Individual Dairy Bulls. Journal of Dairy Science. 1998, 81 (6): 1714-1723.
Watson PF: The causes of reduced fertility with cryopreserved semen. Animal Reproduction Science. 2000, 60–61: 481-492. 10.1016/S0378-4320(00)00099-3.
Lucy MC: Reproductive loss in high-producing dairy cattle: where will it end?. Journal of dairy science. 2001, 84 (6): 1277-1293.
Aitken RJ: Sperm function tests and fertility. International Journal of Andrology. 2006, 29 (1): 69-75. 10.1111/j.1365-2605.2005.00630.x.
Braundmeier AG, Miller DJ: Invited Review: The Search is on: Finding Accurate Molecular Markers of Male Fertility. Journal of Dairy Science. 2001, 84 (9): 1915-1925.
Ballachey BE, Evenson DP, Saacke RG: The sperm chromatin structure assay. Relationship with alternate tests of semen quality and heterospermic performance of bulls. Journal of Andrology. 1988, 9 (2): 109-115.
Dejarnette JM: The effect of semen quality on reproductive efficiency. Vet Clin North Am Food Anim Pract. 2005, 21 (2): 409-418. 10.1016/j.cvfa.2005.02.011.
Saacke RG, Dejarnette JM, Bame JH, Karabinus DS, Whitman SS: Can spermatozoa with abnormal heads gain access to the ovum in artificially inseminated super- and single-ovulating cattle?. Theriogenology. 1998, 50 (1): 117-128. 10.1016/S0093-691X(98)00119-8.
Saacke RG, Nadir S, Nebel RL: Relationship of semen quality to sperm transport, fertilization, and embryo quality in ruminants. Theriogenology. 1994, 41 (1): 45-50. 10.1016/S0093-691X(05)80047-0.
Feugang JM, Rodriguez-Osorio N, Memili E: Advances and applications of genomics and proteomics in gamete and embryo science. Acta Scientiae Veterinariae. 2007, 35 (Supl. 3): s847-s856.
Georges M: Mapping, Fine Mapping, and Molecular Dissection of Quantitative Trait Loci in Domestic Animals. Annual Review of Genomics and Human Genetics. 2007, 8 (1): 131-162. 10.1146/annurev.genom.8.080706.092408.
Sellner EM, Kim JW, McClure MC, Taylor KH, Schnabel RD, Taylor JF: BOARD-INVITED REVIEW: Applications of genomic information in livestock. Journal of Animal Science. 2007, 85 (12): 3148-3158. 10.2527/jas.2007-0291.
Misirlioglu M, Page G, Sagirkaya H, Kaya A, Parrish J, First N, Memili E: Dynamics of global transcriptome in bovine matured oocytes and preimplantation embryos. Proc Natl Acad Sci USA. 2006, 103 (50): 18905-18910. 10.1073/pnas.0608247103.
Li L-W, Fan L-Q, Zhu W-B, Nie H-C, Sun B-L, Luo K-L, Liao T-T, Tang L, Lu G-X: Establishment of a high-resolution 2-D reference map of human spermatozoal proteins from 12 fertile sperm-bank donors. Asian Journal of Andrology. 2007, 9 (3): 321-329. 10.1111/j.1745-7262.2007.00261.x.
Martínez-Heredia J, Maria Estanyol J, Luis Ballescà J, Oliva R: Proteomic identification of human sperm proteins. PROTEOMICS. 2006, 6 (15): 4356-4369. 10.1002/pmic.200600094.
Miller D, Ostermeier GC: Spermatozoal RNA: why is it there and what does it do?. Gynecologie Obstetrique & Fertilite. 2006, 34 (9): 840-846. 10.1016/j.gyobfe.2006.07.013.
Craig DW, Stephan DA: Applications of whole-genome high-density SNP genotyping. Expert Review of Molecular Diagnostics. 2005, 5 (2): 159-170. 10.1586/14737220.127.116.11.
Newton-Cheh C, Hirschhorn JN: Genetic association studies of complex traits: design and analysis issues. Mutat Res. 2005, 573 (1–2): 54-69.
Suh Y, Vijg J: SNP discovery in associating genetic variation with human disease phenotypes. Mutation Research/Fundamental and Molecular Mechanisms of Mutagenesis. 2005, 573 (1–2): 41-53. 10.1016/j.mrfmmm.2005.01.005.
Tishkoff SA, Reed FA, Ranciaro A, Voight BF, Babbitt CC, Silverman JS, Powell K, Mortensen HM, Hirbo JB, Osman M, et al: Convergent adaptation of human lactase persistence in Africa and Europe. Nat Genet. 2007, 39 (1): 31-40. 10.1038/ng1946.
Amills M, Vidal O, Varona L, Tomas A, Gil M, Sanchez A, Noguera JL: Polymorphism of the pig 2,4-dienoyl CoA reductase 1 gene (DECR1) and its association with carcass and meat quality traits. J Anim Sci. 2005, 83 (3): 493-498.
Brym P, Kaminski S, Wojcik E: Nucleotide sequence polymorphism within exon 4 of the bovine prolactin gene and its associations with milk performance traits. J Appl Genet. 2005, 46 (2): 179-185.
Haegeman A, Williams JL, Law A, Van Zeveren A, Peelman LJ: Mapping and SNP analysis of bovine candidate genes for meat and carcass quality. Animal Genetics. 2003, 34 (5): 349-353. 10.1046/j.1365-2052.2003.01008.x.
Heaton MP, Harhay GP, Bennett GL, Stone RT, Grosse WM, Casas E, Keele JW, Smith TP, Chitko-McKown CG, Laegreid WW: Selection and use of SNP markers for animal identification and paternity analysis in U.S. beef cattle. Mamm Genome. 2002, 13 (5): 272-281. 10.1007/s00335-001-2146-3.
Vallet JL, Freking BA, Leymaster KA, Christenson RK: Allelic variation in the erythropoietin receptor gene is associated with uterine capacity and litter size in swine. Anim Genet. 2005, 36 (2): 97-103. 10.1111/j.1365-2052.2005.01233.x.
Vychodilova-Krenkova L, Matiasovic J, Horin P: Single nucleotide polymorphisms in four functionally related immune response genes in the horse: CD14, TLR4, Cepsilon, andFcepsilon R1 alpha. Int J Immunogenet. 2005, 32 (5): 277-283. 10.1111/j.1744-313X.2005.00522.x.
Werner FA, Durstewitz G, Habermann FA, Thaller G, Kramer W, Kollers S, Buitkamp J, Georges M, Brem G, Mosner J, et al: Detection and characterization of SNPs useful for identity control and parentage testing in major European dairy breeds. Anim Genet. 2004, 35 (1): 44-49. 10.1046/j.1365-2052.2003.01071.x.
Michal JJ, Zhang ZW, Gaskins CT, Jiang Z: The bovine fatty acid binding protein 4 gene is significantly associated with marbling and subcutaneous fat depth in Wagyu × Limousin F2 crosses. Anim Genet. 2006, 37 (4): 400-402. 10.1111/j.1365-2052.2006.01464.x.
Aoki VW, Christensen GL, Atkins JF, Carrell DT: Identification of novel polymorphisms in the nuclear protein genes and their relationship with human sperm protamine deficiency and severe male infertility. Fertility and Sterility. 2006, 86 (5): 1416-1422. 10.1016/j.fertnstert.2006.04.033.
Tanaka H, Hirose M, Tokuhiro K, Matsuoka Y, Miyagawa Y, Tsujimura A, Okuyama A, Nishimune Y: Single nucleotide polymorphisms: discovery of the genetic causes of male infertility. Soc Reprod Fertil Suppl. 2007, 65: 531-534.
Tanaka H, Miyagawa Y, Tsujimura A, Matsumiya K, Okuyama A, Nishimune Y: Single nucleotide polymorphisms in the protamine-1 and -2 genes of fertile and infertile human male populations. Molecular Human Reproduction. 2003, 9 (2): 69-73. 10.1093/molehr/gag010.
Tschanter P, Kostova E, Luetjens CM, Cooper TG, Nieschlag E, Gromoll J: No association of the A260G and A386G DAZL single nucleotide polymorphisms with male infertility in a Caucasian population. Human Reproduction. 2004, 19 (12): 2771-2776. 10.1093/humrep/deh522.
Peddinti D, Nanduri B, Kaya A, Feugang JM, Burgess SC, Memili E: Comprehensive proteomic analysis of bovine spermatozoa of varying fertility rates and identification of biomarkers associated with fertility. BMC Systems Biology. 2008, 2: 19-10.1186/1752-0509-2-19.
Zwald NR, Weigel KA, Chang YM, Welper RD, Clay JS: Genetic Selection for Health Traits Using Producer-Recorded Data. I. Incidence Rates, Heritability Estimates, and Sire Breeding Values. Journal of Dairy Science. 2004, 87 (12): 4287-4294.
Chang YM, Gianola D, Heringstad B, Klemetsdal G: Effects of trait definition on genetic parameter estimates and sire evaluation for clinical mastitis with threshold models. Animal science. 2004, 79: 3555-3364.
Hardenbol P, Baner J, Jain M, Nilsson M, Namsaraev EA, Karlin-Neumann GA, Fakhrai-Rad H, Ronaghi M, Willis TD, Landegren U, et al: Multiplexed genotyping with sequence-tagged molecular inversion probes. Nat Biotech. 2003, 21 (6): 673-678. 10.1038/nbt821.
Hardenbol P, Yu F, Belmont J, MacKenzie J, Bruckner C, Brundage T, Boudreau A, Chow S, Eberle J, Erbilgin A, et al: Highly multiplexed molecular inversion probe genotyping: Over 10,000 targeted SNPs genotyped in a single tube assay. Genome Research. 2005, 15 (2): 269-275. 10.1101/gr.3185605.
Sagirkaya H, Misirlioglu M, Kaya A, First NL, Parrish JJ, Memili E: Developmental potential of bovine oocytes cultured in different maturation and culture conditions. Animal Reproduction Science. 2007, 101 (3–4): 225-240. 10.1016/j.anireprosci.2006.09.016.
Amann RP, Hammerstedt RH: Detection of Differences in Fertility. Journal of Andrology. 2002, 23 (3): 317-325.
Khatkar MS, Nicholas FW, Collins AR, Zenger KR, Cavanagh JA, Barris W, Schnabel RD, Taylor JF, Raadsma HW: Extent of genome-wide linkage disequilibrium in Australian Holstein-Friesian cattle based on a high-density SNP panel. BMC genomics. 2008, 9: 187-10.1186/1471-2164-9-187.
Healy DG, Abou-Sleiman PM, Goldstein DB, Wood NW: Population genetics for target identification. Drug Discovery Today: Technologies. 2004, 1 (1): 69-74. 10.1016/j.ddtec.2004.08.012.
Taramelli R, Acquati F: The human genome project and the discovery of genetic determinants of cancer susceptibility. European Journal of Cancer. 2004, 40 (17): 2537-2543. 10.1016/j.ejca.2004.07.030.
Iguchi N, Yang S, Lamb DJ, Hecht NB: An SNP in protamine 1: a possible genetic cause of male infertility?. Journal of Medical Genetics. 2006, 43 (4): 382-384. 10.1136/jmg.2005.037168.
Daetwyler HD, Schenkel FS, Sargolzaei M, Robinson JA: A genome scan to detect quantitative trait loci for economically important traits in Holstein cattle using two methods and a dense single nucleotide polymorphism map. Journal of dairy science. 2008, 91 (8): 3225-3236. 10.3168/jds.2007-0333.
Farnir F, Coppieters W, Arranz J-J, Berzi P, Cambisano N, Grisart B, Karim L, Marcq F, Moreau L, Mni M, et al: Extensive Genome-wide Linkage Disequilibrium in Cattle. Genome Research. 2000, 10 (2): 220-227. 10.1101/gr.10.2.220.
McKay SD, Schnabel RD, Murdoch BM, Matukumalli LK, Aerts J, Coppieters W, Crews D, Dias Neto E, Gill CA, Gao C, et al: Whole genome linkage disequilibrium maps in cattle. BMC genetics. 2007, 8: 74-10.1186/1471-2156-8-74.
Odani M, A N, Watanabe T, Yokouchi K, Sugimoto Y, Fujita T, Oguni T, Matsumoto M, Sasaki Y: Genome-wide linkage disequilibrium in two Japanese beef cattle breeds. Animal Genetics. 2006, 37 (2): 139-144. 10.1111/j.1365-2052.2005.01400.x.
Khatkar MS, Zenger KR, Hobbs M, Hawken RJ, Cavanagh JA, Barris W, McClintock AE, McClintock S, Thomson PC, Tier B, et al: A primary assembly of a bovine haplotype block map based on a 15,036-single-nucleotide polymorphism panel genotyped in holstein-friesian cattle. Genetics. 2007, 176 (2): 763-772. 10.1534/genetics.106.069369.
Klein RJ, Zeiss C, Chew EY, Tsai J-Y, Sackler RS, Haynes C, Henning AK, SanGiovanni JP, Mane SM, Mayne ST, et al: Complement Factor H Polymorphism in Age-Related Macular Degeneration. Science. 2005, 308 (5720): 385-389. 10.1126/science.1109557.
Scott EM, Carter AM, Grant PJ: Association between polymorphisms in the Clock gene, obesity and the metabolic syndrome in man. Int J Obes (Lond). 2007
McPherson R, Pertsemlidis A, Kavaslar N, Stewart A, Roberts R, Cox DR, Hinds DA, Pennacchio LA, Tybjaerg-Hansen A, Folsom AR, et al: A Common Allele on Chromosome 9 Associated with Coronary Heart Disease. Science. 2007, 316 (5830): 1488-1491. 10.1126/science.1142447.
Miyagawa Y, Nishimura H, Tsujimura A, Matsuoka Y, Matsumiya K, Okuyama A, Nishimune Y, Tanaka H: Single-Nucleotide Polymorphisms and Mutation Analyses of the TNP1 and TNP2 Genes of Fertile and Infertile Human Male Populations. Journal of Andrology. 2005, 26 (6): 779-786. 10.2164/jandrol.05069.
Anjos SM, Shao W, Marchand L, Polychronakos C: Allelic effects on gene regulation at the autoimmunity-predisposing CTLA4 locus: a re-evaluation of the 3[prime] +6230G>A polymorphism. Genes Immun. 2005, 6 (4): 305-311. 10.1038/sj.gene.6364211.
Leyva-Baca I, Schenkel F, Martin J, Karrow NA: Polymorphisms in the 5' Upstream Region of the CXCR1 Chemokine Receptor Gene, and Their Association with Somatic Cell Score in Holstein Cattle in Canada. Journal of Dairy Science. 2008, 91 (1): 407-417. 10.3168/jds.2007-0142.
Mahr S, Burmester GR, Hilke D, Göbel U, Grützkau A, Häupl T, Hauschild M, Koczan D, Krenn V, Neidel J, et al: Cis- and trans-acting gene regulation is associated with osteoarthritis. Am J Hum Genet. 2006, 78 (5): 793-803. 10.1086/503849.
Goncalves RF, Wolinetz CD, Killian GJ: Influence of arginine-glycine-aspartic acid (RGD), integrins ([alpha]V and [alpha]5) and osteopontin on bovine sperm-egg binding, and fertilization in vitro. Theriogenology. 2007, 67 (3): 468-474. 10.1016/j.theriogenology.2006.08.013.
Campbell S, Swann HR, Seif MW, Kimber SJ, Aplin JD: Integrins and adhesion mlecules: Cell adhesion molecules on the oocyte and preimplantation human embryo. Human Reproduction. 1995, 10 (6): 1571-1578.
Bloor DJ, Metcalfe AD, Rutherford A, Brison DR, Kimber SJ: Expression of cell adhesion molecules during human preimplantation embryo development. Molecular Human Reproduction. 2002, 8 (3): 237-245. 10.1093/molehr/8.3.237.
Zuccotti M, Rossi PG, Fiorillo E, Garagna S, Forabosco A, Redi CA: Timing of Gene Expression and Oolemma Localization of Mouse [alpha]6 and [beta]1 Integrin Subunits during Oogenesis. Developmental Biology. 1998, 200 (1): 27-34. 10.1006/dbio.1998.8923.
Miller DJ, Shi X, Burkin H: Molecular Basis of Mammalian Gamete Binding. Recent Progress in Hormone Research. 2002, 57 (1): 37-73. 10.1210/rp.57.1.37.
Glander H, Schaller J: Beta 1-integrins of spermatozoa: a flow cytophotometric analysis. Int J Androl. 1993, 16 (2): 105-111. 10.1111/j.1365-2605.1993.tb01162.x.
Barraud-Lange V, Naud-Barriant N, Saffar L, Gattegno L, Ducot B, Drillet A, Bomsel M, Wolf J, Ziyyat A: Alpha6beta1 integrin expressed by sperm is determinant in mouse fertilization. BMC Dev Biol. 2007, 7: 102-10.1186/1471-213X-7-102.
Reddy K, Meherji P, Shahani S: Integrin cell adhesion molecules on human spermatozoa. Indian J Exp Biol. 1998, 36 (5): 456-463.
Pate Barry, White Kenneth, Winger Quinton, Rickords Lee, Aston Kenneth, Sessons Benjamin, Li Guang-Peng, Campbell Kenneth, Weimer Bart, Bunch TD: Specific integrin subunits in bovine oocytes, including novel sequences for alpha 6 and beta 3 subunits. Molecular Reproduction and Development. 2007, 74 (5): 600-607. 10.1002/mrd.20649.
This project was funded by the Mississippi Agricultural and Forestry Experiment Station of Mississippi State University (J-11262), and by Alta Genetics Inc.
JMF carried out the sperm DNA isolation, performed the bioinformatics and functional analyses, and participated in the manuscript writing. AK and ET carried out the collection and analysis of fertility data and contributed to the study design and manuscript writing. GPP, LC, and TM participated in the design of the study and performed the statistical analysis, the SNP genotyping, and allelic discrimination data, and participated in manuscript writing. KH, LN, and RG performed the SNP genotyping and contributed to the manuscript writing. EM participated in the design of the study, coordinated the study and participated in manuscript writing.
Electronic supplementary material
Additional File 1: Table S1. Mean fertility (in Z-score), standard deviation, minimum and maximum values for the four markers from the Phase II study. (DOC 46 KB)
About this article
Cite this article
Feugang, J.M., Kaya, A., Page, G.P. et al. Two-stage genome-wide association study identifies integrin beta 5 as having potential role in bull fertility. BMC Genomics 10, 176 (2009). https://doi.org/10.1186/1471-2164-10-176
- Single Nucleotide Polymorphism
- Integrin Beta
- Single Nucleotide Polymorphism Marker
- Single Nucleotide Polymorphism Rs41257187
- Single Nucleotide Polymorphism Location