INCENP基因功能性单倍型可调控启动子活性并影响公牛精液品质

为研究内着丝粒蛋白(Inner centromere protein, INCENP)基因启动子区单核苷酸多态性(SNPs)与精液品质的相关性,本文利用聚合酶链反应-限制性片段长度多态性(PCR-RFLP)方法检测了250头中国荷斯坦公牛INCENP基因的基因型。在INCENP基因启动子区鉴定出两个SNPs (g.-556 G>T,rs 136823901和g.-692 C>T,rs 211010999),发现了3种单倍型(CG、TT、TG)。分析两个SNP位点的基因型频率和等位基因频率,各SNP及单倍型组合与中国荷斯坦公牛精液品质的相关性,结果表明SNP位点g.-556 G>T GT基因型个体的鲜精活力显著高于GG基因型个体(P<0.05),单倍型组合H1H1(CCGG)、H1H3(CTGT)、H2H3(TTGT)和H3H3(TTTT)个体的鲜精活力和冻精解冻后活力均显著高于H1H2个体(P<0.05)。为进一步研究g.-556 G>T和g.-692 C>T影响精液品质的可能机理,本文将3种单倍型质粒分别转染小鼠睾丸间质细胞(MLTC-1),结果显示含TG单倍型的载体荧光素酶活性最高。由此推测,g.-556 G>T和g.-692 C>T为启动子区功能性突变位点,可通过调节启动子活性来调控INCENP基因表达,进而影响精液品质。     

Genetic variations of HSBP1 gene and its effect on thermal performance traits in Chinese Holstein cattle

Molecular chaperones have been understood to be preferentially transcribed to prevent perturbations in response to various stresses. In this study, three single nucleotide polymorphisms (SNPs), g.324G>C, g.589C>T and g.651C>G in Heat shock factor binding protein 1 (HSBP1) gene were found and genotyped in 930 Chinese Holstein cattle. The results indicated that only g.589C>T polymorphism locus met Hardy–Weinberg equilibrium (P > 0.05). Pair linkage disequilibrium analysis and haplotype construction of HSBP1 gene were performed using SHEsis software. Seven haplotypes were constructed and fourteen haplotype combinations were found. Association analysis showed that H2H2 haplotype combination was advantageous for thermo tolerance breeding in Chinese Holstein. The cows with H2H2 haplotype combination have lower decrease rate of milk yield than those with H2H3 haplotype combination (P < 0.05) and lower potassium content in erythrocytes (PCE) than those with H2H5 (P < 0.05), H4H4 (P < 0.05) and H4H5 (P < 0.01) haplotype combination. The association between SNP and thermo tolerance traits showed that PCE of cows with GG genotype was lower than those with CG genotype at g.651C>G locus (P < 0.01). Pair linkage disequilibrium analysis revealed that the three loci were at a strong disequilibrium state. So we presumed that the effect of H2H2 haplotype combination on thermo tolerance traits major due to the SNP of g.651C>G.     

Six novel single-nucleotide polymorphisms in SPAG11 gene and their association with sperm quality traits in Chinese Holstein bulls

Sperm-associated antigen 11 (SPAG11) is predominant in the male reproductive tract. Similar to β-defensin, aside from its antibacterial activity, SPAG11 also has an important role in male reproductive function. In the present study, the association of bovine SPAG11 gene polymorphism with sperm quality traits was examined, including ejaculate volume, sperm concentration, fresh sperm motility, post-thaw cryopreserved sperm motility, and deformity rate of bull semen. Six novel single nucleotide polymorphisms (SNPs) of the SPGA11 gene were investigated in 426 normal mature Chinese Holstein bulls using polymerase chain reaction-single strand conformation polymorphism (PCR-SSCP), PCR-restriction fragment length polymorphism (PCR-RFLP), created restriction site-PCR (CRS-PCR), and DNA sequencing methods. Linkage disequilibrium analysis showed that g.1306G>A and g.1454G>A (SNP-1), and g.16904G>T, g.16974C>T, and g.17000A>G (SNP-2) are completely linked, respectively. Correlation analysis showed the SNP-2 marker had a marked effect on fresh sperm motility and sperm concentration (P<0.05). SNP-3 g.22696T>C had a marked effect on post-thaw cryopreserved sperm motility (P<0.05) and deformity rate (P<0.01). However, the presence of SNP-1 was not correlated with the sperm production traits (P>0.05). Furthermore, association analyses of the 8 haplotypes constructed from the 17 combined haplotypes and reproductive traits showed that the bulls with the combined haplotype H5H6 (GGT/TTC) have the highest ejaculate volumes and the bulls with combined haplotypes H1H1 (AAT/TTT) and H1H6 (AGT/TTC) had the highest fresh and post-thaw sperm motilities, respectively. These results indicate that new molecular markers associated with sperm quality traits can be used in marker-assisted selection in bull breeding programs.     

Genome‐wide association study for semen traits of the bulls in Chinese Holstein

A genome‐wide association study (GWAS) was performed to identify markers and candidate genes for five semen traits in the Holstein bull population in China. The analyzed dataset consisted of records from 692 bulls from eight bull stations; each bull was genotyped using the Illumina BovineSNP50 BeadChip. Association tests between each trait and the 41 188 informative high‐quality SNPs were achieved with gapit software. In total, 19 suggestive significant SNPs, partly located within the reported QTL regions or within or close to the reported candidate genes, associated with five semen traits were detected. By combining our GWAS results with the biological functions of these genes, eight novel promising candidate genes, including ETNK1, PDE3A, PDGFRB, CSF1R, WT1, DSCAML1, SOD1 and RUNX2, were identified that potentially relate to semen traits. Our findings may provide a basis for further research on the genetic mechanism of semen traits and marker‐assisted selection of such traits in Holstein bulls.     

Polymorphism in the exon 4 of β-lactoglobulin variant B precursor gene and its association with milk traits and protein structure in Chinese Holstein

β-lactoglobulin (β-LG) is the major whey protein in the milk. In order to investigate the polymorphism of β-LG variant B precursor (β-LG B*: GenBank accession no. DQ489319) gene and its effects on the milk traits, the single-strand conformation polymorphism method (PCR-SSCP) were adopted to analyze polymorphism between 5229th and 5476th bp in the β-LG B* gene in Chinese Holstein. Four genotypes were found (AA, AB, AC and ABC) and 3 single nucleotide polymorphisms (SNPs) were detected (g.5239C>A, g.5240A>C, g.5305C>T and mix type g.5305C/T) in the exon 4 of β-LG B* gene. It was also found that the protein contents of AB, AC and ABC dairy cows were higher than AA (P < 0.05), and AC cows were the highest among them. Three SNPs (g.5239C>A, g.5240A>C and g.5305C>T) might affect the milk trait and all of them were high polymorphism (0.5 < PIC < 1.0). In further researches, the three SNPs also caused amino acid change (Asp>Glu, Thr>Pro and Ala>Val) respectively, and the spatial secondary and tertiary structure forecasting result also showed that single amino acid change influence protein spatial structure change in Chinese Holstein. Taken together, it is suggested that these SNPs change β-LG B* gene structure and expression. The polymorphism possibly holds the secret of milk protein and fat contents in the milk of Chinese Holstein.     

Effect of FASN gene on milk yield and milk composition in the Chinese Holstein dairy population

Fatty acid synthase (FASN) is a multifunctional protein that catalyzes de novo synthesis of fatty acids in cells. It plays a key role in the lipid biosynthesis as well as in the general metabolism of all living animals. We herein investigated polymorphisms of FASN. As a result, six single nucleotide polymorphisms (SNPs) were found and then genotyped in 752 Chinese Holstein cows. It was found that g.17924A>G was non‐synonymous, g.13965 C>T, g.16907 T>C and g.18663T>C were synonymous mutations and two other two SNPs, g.8948 C>T (ss491228481) and g.14439T>C (rs133498277), were in intronic sequences of the gene. All such identified SNPs were found to be associated with milk yield and composition traits (P = 0.0441 to <0.0001). Significant additive and allele substitution effects were observed for three yield traits at all six loci as well (P < 0.05 to <0.01). Complete linkage disequilibrium among the five SNPs, with the exception of g.8948 C>T, was observed.     

A Function Accounting for Training Set Size and Marker Density to Model the Average Accuracy of Genomic Prediction

Prediction of genomic breeding values is of major practical relevance in dairy cattle breeding. Deterministic equations have been suggested to predict the accuracy of genomic breeding values in a given design which are based on training set size, reliability of phenotypes, and the number of independent chromosome segments (). The aim of our study was to find a general deterministic equation for the average accuracy of genomic breeding values that also accounts for marker density and can be fitted empirically. Two data sets of 5′698 Holstein Friesian bulls genotyped with 50 K SNPs and 1′332 Brown Swiss bulls genotyped with 50 K SNPs and imputed to ∼600 K SNPs were available. Different k-fold (k = 2–10, 15, 20) cross-validation scenarios (50 replicates, random assignment) were performed using a genomic BLUP approach. A maximum likelihood approach was used to estimate the parameters of different prediction equations. The highest likelihood was obtained when using a modified form of the deterministic equation of Daetwyler et al. (2010), augmented by a weighting factor (w) based on the assumption that the maximum achievable accuracy is . The proportion of genetic variance captured by the complete SNP sets () was 0.76 to 0.82 for Holstein Friesian and 0.72 to 0.75 for Brown Swiss. When modifying the number of SNPs, w was found to be proportional to the log of the marker density up to a limit which is population and trait specific and was found to be reached with ∼20′000 SNPs in the Brown Swiss population studied.     

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.     

Role of the DLGAP2 Gene Encoding the SAP90/PSD-95-Associated Protein 2 in Schizophrenia

Aberrant synaptic dysfunction is implicated in the pathogenesis of schizophrenia. The DLGAP2 gene encoding the SAP90/PSD-95-associated protein 2 (SAPAP2) located at the post-synaptic density of neuronal cells is involved in the neuronal synaptic function. This study aimed to investigate whether the DLGAP2 gene is associated with schizophrenia. We resequenced the putative promoter region and all the exons of the DLGAP2 gene in 523 patients with schizophrenia and 596 non-psychotic controls from Taiwan and conducted a case-control association analysis. We identified 19 known SNPs in this sample. Association analysis of 9 SNPs with minor allele frequency greater than 5% showed no association with schizophrenia. However, we found a haplotype (CCACCAACT) significantly associated with schizophrenia (odds ratio:2.5, p<0.001). We also detected 16 missense mutations and 1 amino acid-insertion mutation in this sample. Bioinformatic analysis showed some of these mutations were damaging or pathological to the protein function, but we did not find increased burden of these mutations in the patient group. Notably, we identified 5 private rare variants in 5 unrelated patients, respectively, including c.−69+9C>T, c.−69+13C>T, c.−69+47C>T, c.−69+55C>T at intron 1 and c.−32A>G at untranslated exon 2 of the DLGAP2 gene. These rare variants were not detected in 559 control subjects. Further reporter gene assay of these rare variants except c.−69+13C>T showed significantly elevated promoter activity than the wild type, suggesting increased DLGAP2 gene expression may contribute to the pathogenesis of schizophrenia. Our results indicate that DLGAP2 is a susceptible gene of schizophrenia.     

Single-nucleotide polymorphisms g.151435C>T and g.173057T>C in PRLR gene regulated by bta-miR-302a are associated with litter size in goats

Single-nucleotide polymorphisms (SNPs) located at microRNA-binding sites (miR-SNPs) can affect the expression of genes. This study aimed to identify the miR-SNPs associated with litter size. Guanzhong (n = 321) and Boer (n = 191) goat breeds were used to detect SNPs in the caprine prolactin receptor (PRLR) gene by DNA sequencing, primer-introduced restriction analysis-polymerase chain reaction, and polymerase chain reaction-restriction fragment length polymorphism. Three novel SNPs (g.151435C>T, g.151454A>G, and g.173057T>C) were identified in the caprine PRLR gene. Statistical results indicated that the g.151435C>T and g.173057T>C SNPs were significantly associated with litter size in Guanzhong and Boer goat breeds. Further analysis revealed that combinative genotype C6 (TTAACC) was better than the others for litter size in both goat breeds. Furthermore, the PRLR g.173057T>C polymorphism was predicted to regulate the binding activity of bta-miR-302a. Luciferase reporter gene assay confirmed that 173057C to T substitution disrupted the binding site for bta-miR-302a, resulting in the reduced levels of luciferase. Taken together, these findings suggested that bta-miR-302a can influence the expression of PRLR protein by binding with 3′untranslated region, resulting in that the g.173057T>C SNP had significant effects on litter size.     

Systematic Functional Study of Cytochrome P450 2D6 Promoter Polymorphisms in the Chinese Han Population

The promoter polymorphisms of drug-metabolizing genes can lead to interindividual differences in gene expression, which may result in adverse drug effects and therapeutic failure. Based on the database of CYP2D6 gene polymorphisms in the Chinese Han population established by our group, we functionally characterized the single nucleotide polymorphisms (SNPs) of the promoter region and corresponding haplotypes in this population. Using site-directed mutagenesis, all the five SNPs identified and ten haplotypes with a frequency equal to or greater than 0.01 in the population were constructed on a luciferase reporter system. Dual luciferase reporter systems were used to analyze regulatory activity. The activity produced by Haplo3(−2183G>A, −1775A>G, −1589G>C, −1431C>T, −1000G>A, −678A>G), Haplo8(−2065G>A, −2058T>G, −1775A>G, −1589G>C, −1235G>A, −678A>G) and MU3(−498C>A) was 0.7−, 0.7−, 1.2− times respectively compared with the wild type in human hepatoma cell lines(p<0.05). These findings might be useful for optimizing pharmacotherapy and the design of personalized medicine.     

Effect of CYP2B6 Gene Polymorphisms on Efavirenz Plasma Concentrations in Chinese Patients with HIV Infection

Objectives

The main aim of this study was to investigate the effect of CYP2B6 gene polymorphisms on efavirenz (EFV) plasma concentrations in Han Chinese patients with human immunodeficiency virus (HIV) infection.

Methods

In total, 322 patients were recruited for study. EFV plasma concentrations at steady-state were determined using high-performance liquid chromatography. Genotyping for seven single nucleotide polymorphisms (SNPs), including 171+967C>A, 171+3212C>T, 171+4335T>C, 516G>T, 785A>G, 1295-913G>A, and *1355A>G of CYP2B6, was performed using ligase detection reaction (LDR). SPSS 18.0 and Haploview 4.2 were applied for statistical analyses.

Results

The average EFV concentration of patients was 2.35±2.09 μg/mL. Overall, 22% patients displayed EFV concentrations out of the therapeutic range of 1–4 μg/mL (13.1% < 1 μg/mL, 9.3% > 4 μg/mL). We observed significant association of 171+967C>A, 171+4335T>C, 516G>T, 785A>G and *1355A>G with high plasma EFV levels (p<.01). The predictive accuracy values of 171+4335CC, 516TT and 785GG for EFV concentrations > 4 μg/mL were 56.7%, 56.7% and 60%, respectively. We observed strong linkage disequilibrium for 171+967C>A, 171+4335T>C, 516G>T and 785A>G, resulting in five haplotypes. The frequencies of the five haplotypes (high to low) were as follows: CCTG (0.328), ACTG (0.280), ACCT (0.189), ATTG (0.186) and ACCG (0.017). The frequency of CCTG (0.524) in patients with EFV plasma concentrations < 1 μg/mL was significantly higher than that in other patient groups, while that of ACCT (0.733) was significantly higher in patients with EFV concentrations > 4 μg/mL, relative to other patient groups. Average EFV concentrations of patients carrying ACTG (1.78 μg/mL), ACCT (7.50 μg/mL), and ATTG (1.92 μg/mL) haplotypes were markedly higher than those of patients carrying the CCTG haplotype. The predictive accuracy of ACCT for EFV > 4 μg/mL was 81%.

Conclusions

Chinese patients administered standard doses of EFV require therapeutic drug monitoring or personalized medication management. Based on the current findings, we propose that 171+4335T>C, 516G>T, 785A>G and haplotype ACCT may be effectively used as genomic markers for EFV, which should aid in improving the efficacy of EFV-containing treatments and reduce the incidence of adverse reactions.     

Polymorphisms in luteinizing hormone receptor and hypothalamic gonadotropin-releasing hormone genes and their effects on sperm quality traits in Chinese Holstein bulls

Genes of hypothalamic-pituitary-gonadal axis play a key role in male reproductive performance. This study evaluated the polymorphisms of luteinizing hormone receptor (LHR) and hypothalamic gonadotropin-releasing hormone (GnRH) genes and their effects on sperm quality traits including semen volume per ejaculate (VOL), sperm density (SD), fresh sperm motility (FSM), thawed sperm motility (TSM), acrosome integrity rate (AIR), and abnormal sperm rate (ASR) collected from 205 Chinese Hostein bulls. The study bulls consisted of 205 mature Chinese Holstein, 27 Simmental, 28 Charolais, and 14 German yellow cattle. One single nucleotide polymorphism (SNP) (A883G) in exon 2 of GnRH and two SNPs (A51703G and G51656T) in intron 9 of LHR were identified in 274 bulls. Analysis of variance in 205 Chinese Holstein bulls showed that age had significant effect on both SD and FSM (P < 0.01), and ASR (P < 0.05). With regards to genotype and its interaction with age, only the SNP of G51656T in LHR gene had significant effect on SD (P < 0.05, P < 0.01; respectively). The association result showed that bulls with AG genotype had higher FSM than bulls with AA and GG genotype in LHR at 51,703 locus (P < 0.10), and bulls with GG genotype had higher SD than bulls with TT genotype in LHR at G51656T locus (P < 0.10). Phenotypic correlation among the traits revealed that significant negative correlations were observed between ASR and AIR (r = -0.736, P < 0.01), ASR and AIR (r = -0.500, P < 0.01). There were moderate positive correlations between VOL and SD (r = 0.422, P < 0.01), as well as FSM (r = 0.411, P < 0.01). In conclusion, LHR may be a potential marker for sperm quality of SD and FSM.     

Sequence Variations in the Bovine Growth Hormone Gene Characterized by Single-Strand Conformation Polymorphism (Sscp) Analysis and Their Association with Milk Production Traits in Holsteins

Sequence variations in the bovine growth hormone (GH) gene were investigated by single strand conformation polymorphism (SSCP) analysis of seven amplified fragments covering almost the entire gene (2.7 kb). SSCPs were detected in four of these fragments and a total of six polymorphisms were found in a sample of 128 Holstein bulls. Two polymorphisms, a T->C transition in the third intron (designated GH4.1) and an A->C transversion in the fifth exon (designated GH6.2), were shown to be associated with milk production traits. GH4.1(c)/GH4.1(c) bulls had higher milk yield than GH4.1(c)/GH4.1(t) (P <= 0.005) and GH4.1(t)/GH4.1(t) (P <= 0.0022) bulls. GH4.1(c)/GH4.1(c) bulls had higher kg fat (P <= 0.0076) and protein (P <= 0.0018) than GH4.1(c)/GH4.1(t) bulls. Similar effects on milk production traits with the GH6.2 polymorphism were observed with the GH6.2(a) allele being the favorable allele. The average effects of the gene substitution for GH4.1 and GH6.2 are similar, with +/-300 kg for milk yield, +/-8 kg for fat content and +/-7 kg for protein content per lactation. The positive association of GH4.1(c) and GH6.2(a) with milk production traits may be useful for improving milk performance in dairy cattle.     

Association of IL10 Polymorphisms and Leprosy: A Meta-Analysis

Leprosy is a chronic infectious disease that depends on the interplay of several factors. Single nucleotide polymorphisms (SNPs) in host immune related genes have been consistently suggested as participants in susceptibility towards disease. Interleukin-10 (IL-10) is a crucial immunomodulatory cytokine in mycobacterial pathogenesis and especially the -819C>T SNP (rs1800871) has been tested in several case-control studies indicating association with leprosy risk, although a recent consensus estimate is still missing. In this study, we evaluated the association of the -819C>T SNP and leprosy in two new Brazilian family-based populations. Then, we performed meta-analysis for this polymorphism summarizing published studies including these Brazilian family-based groups. Finally, we also retrieved published studies for other distal and proximal IL10 polymorphisms: -3575 T>A (rs1800890), -2849 G>A (rs6703630), -2763 C>A (rs6693899), -1082 G>A (rs1800896) and -592 C>A (rs1800872). Results from meta-analysis supported a significant susceptibility association for the -819T allele, with pooled Odds Ratio of 1.22 (CI = 1.11–1.34) and P-value = 3x10^–5 confirming previous data. This result remained unaltered after inclusion of the Brazilian family-based groups (OR = 1.2, CI = 1.10–1.31, P-value = 2x10^–5). Also, meta-analysis confirmed association of -592 A allele and leprosy outcome (OR = 1.24, CI = 1.03–1.50, P-value = 0.02). In support of this, linkage disequilibrium analysis in 1000 genomes AFR, EUR, ASN and AMR populations pointed to r² = 1.0 between the -592C>A and -819C>T SNPs. We found no evidence of association for the other IL10 polymorphisms analyzed for leprosy outcome. Our results reinforce the role of the -819C>T as a tag SNP (rs1800871) and its association with leprosy susceptibility.     

Effects of <Emphasis Type="Italic">Mbo</Emphasis>II and <Emphasis Type="Italic">BspM</Emphasis>I polymorphisms in the gonadotropin releasing hormone receptor (<Emphasis Type="Italic">GnRHR</Emphasis>) gene on sperm quality in Holstein bulls

The hypothalamic gonadotropin-releasing hormone receptor (GnRHR) plays an essential physiological role in reproductive function, which triggers the synthesis and release of luteinizing hormone and follicle stimulating hormone in the pituitary. The objective of this study was to investigate the effects of polymorphisms of GnRHR gene on the quality of fresh and frozen semen in Holstein bulls. The PCR-RFLP method was applied to detect G286A and T340C transitions determining MboII and BspMI polymorphisms, respectively, in the exon I of bovine GnRHR gene and evaluated its associations with sperm quality traits in 131 Holstein bulls. In polymorphic locus 286, bulls with the GA genotype had significantly higher sperm motility in frozen semen (FMOT) than GG genotype (P < 0.01). In polymorphic locus 340, bulls with heterozygote CT genotype had significantly higher sperm motility (MOT), semen volume per ejaculate (VOL), and lower abnormal spermatozoa rate (ASR) than homozygote TT genotype (P < 0.05). Bulls contained one A allele or C allele had a favorable, positive effect on sperm quality traits. These results indicate that GnRHR gene can be a potential marker for improving sperm quality traits, and imply that bulls with GA or CT genotype should be selected in breeding program.     

A Nonsense Mutation in TMEM95 Encoding a Nondescript Transmembrane Protein Causes Idiopathic Male Subfertility in Cattle

Genetic variants underlying reduced male reproductive performance have been identified in humans and model organisms, most of them compromising semen quality. Occasionally, male fertility is severely compromised although semen analysis remains without any apparent pathological findings (i.e., idiopathic subfertility). Artificial insemination (AI) in most cattle populations requires close examination of all ejaculates before insemination. Although anomalous ejaculates are rejected, insemination success varies considerably among AI bulls. In an attempt to identify genetic causes of such variation, we undertook a genome-wide association study (GWAS). Imputed genotypes of 652,856 SNPs were available for 7962 AI bulls of the Fleckvieh (FV) population. Male reproductive ability (MRA) was assessed based on 15.3 million artificial inseminations. The GWAS uncovered a strong association signal on bovine chromosome 19 (P = 4.08×10⁻⁵⁹). Subsequent autozygosity mapping revealed a common 1386 kb segment of extended homozygosity in 40 bulls with exceptionally poor reproductive performance. Only 1.7% of 35,671 inseminations with semen samples of those bulls were successful. None of the bulls with normal reproductive performance was homozygous, indicating recessive inheritance. Exploiting whole-genome re-sequencing data of 43 animals revealed a candidate causal nonsense mutation (rs378652941, c.483C>A, p.Cys161X) in the transmembrane protein 95 encoding gene TMEM95 which was subsequently validated in 1990 AI bulls. Immunohistochemical investigations evidenced that TMEM95 is located at the surface of spermatozoa of fertile animals whereas it is absent in spermatozoa of subfertile animals. These findings imply that integrity of TMEM95 is required for an undisturbed fertilisation. Our results demonstrate that deficiency of TMEM95 severely compromises male reproductive performance in cattle and reveal for the first time a phenotypic effect associated with genomic variation in TMEM95.