PLCz Functional Haplotypes Modulating Promoter Transcriptional Activity Are Associated with Semen Quality Traits in Chinese Holstein Bulls

The sperm-specific phospholipase C zeta (PLCz) is a candidate sperm-borne oocyte-activating factor that triggers a characteristic series of physiological stimuli via cytoplasmic Ca2+ oscillations during fertilization. The molecular mechanisms involved in the regulation of PLCz gene expression remain largely unknown. To explore the genetic variations in the 5′-flanking region of the PLCz gene and their common haplotypes in Chinese Holstein bulls, as well as to determine whether these variations affect bovine semen quality traits and transcriptional activity, DNA samples were collected from Chinese Holstein bulls and sequenced for the identification of genetic variants in the 5′-flanking region of PLCz. Two genetic variants were identified, and their haplotypic profiles were constructed. The two novel genetic variations (g. −456 G>A and g. +65 T>C) were genotyped in 424 normal Chinese Holstein bulls. Bioinformatics analysis revealed that both loci are in transcription factor binding sites of the core promoter region. The association studies revealed that the two genetic variations and their haplotype combinations significantly affected semen quality traits. Using serially truncated constructs of the bovine PLCz promoters and the luciferase reporter, we found that a 726 bp (−641 nt to +112 nt) fragment constitutes the core promoter region. Furthermore, four haplotypes, H1H1 (GTGT), H2H2 (GCGC), H3H3 (ATAT), and H4H4 (ACAC), were significantly associated with semen quality traits and successfully transfected into MLTC-1 cell lines. The luciferase reporter assay showed that the different haplotypes exhibited distinct promoter activities. Maximal promoter activity was demonstrated by the H2H2 haplotypes, as compared with the other haplotypes. To the best of our knowledge, this study is the first report on genetic variants and their respective haplotypes in the 5′-flanking region of PLCz gene that can influence the semen quality of Chinese Holstein bulls as well as contribute to the transcriptional activity of the PLCz promoter.


Introduction
Fertility is one of the most economically important production traits in cattle improvement. However, fertility is greatly affected by the genetic make-up and the environment. Artificial insemination (AI) is known as one of the most valuable techniques for the genetic improvement of semen quality traits in dairy cattle herds [1]. However, the quantity and quality of semen from bull (including the ejaculate volume, sperm density, initial sperm motility, frozen semen motility, and the sperm deformity rate) are affected by heredity, breeding conditions, the physiological status of the parents, temperature, and other factors. Direct selection for semen quality traits is difficult because of their low heritability [2]. With the progress of genetic and molecular techniques, genetic variations at the DNA level are now key elements for the selection of animals via marker-assisted selection (MAS). Several reports have described the use of candidate genes for MAS to influence semen quality and fertility in boars [3][4][5], goats [6], and bulls [7][8][9][10][11][12].
Daghigh-Kia [28] reported a polymorphism specific to bulls at position 2749 (G.A) in intron 6 of the PLCz gene. Their association study between the genotypes and sperm quality traits showed that no significant genotype effects (P . 0.05) could be observed in black and red Holstein populations [28]. However, the genetic variants in the 59-flanking region of the PLCz gene have not yet been functionally characterized. Furthermore, little is known about the regulation of PLCz expression in bulls. The present study aimed to investigate potentially functional genetic variants in the 59-flanking region of the PLCz gene, to characterize the distribution in a population, and to analyze their effects on semen quality traits in Chinese Holstein bulls. The results indicated that novel molecular markers associated with semen quality traits can be used for MAS in bull breeding programs.

Animals and cell lines
Semen samples were collected from a total of 424 Chinese Holstein bulls from three bull stations, namely, 142 bulls from the Beijing Dairy Center, 218 bulls from the Shanghai Bright Dairy and Food Co., Ltd., and 64 bulls from the Shandong OX Bio-Technology Co., Ltd. Semen sample collection was permitted by the owners of the animals. Samples were collected by employees of the respective companies. This study was approved by the Bureau of Animal Husbandry and Veterinary and the Dairy Cattle Frozen Semen Quality Supervision Testing Center of the Chinese Ministry of Agriculture.
The 424 genotyped bulls were sons of 230 sires (each sire had 1 to 10 sons, average of 1.84 sons). Repeated measurements of the sperm quality traits from 2007 to 2011 were available for each bull. Data on 40121 ejaculates from the 424 bulls were used, with 4 to 379 ejaculates from each bull. Semen was collected from each bull at 3 d to 6 d intervals using an artificial vagina. Immediately after collection, the ejaculates were stored at 37uC in a water bath prior to evaluation of the fresh semen quality traits. These traits included the semen volume per ejaculate (in ml), sperm motility (in%), sperm concentration (expressed as 610 8 /ml), and the percentage of abnormal sperm (in%). The fresh semen was then diluted with a glycerol-egg yolk-citrate mixture in the Beijing and Shanghai bull stations or with Bioxcell (IMV Biotechnology, L'Aigle, France) in the Shandong bull station. The collected semen were packaged in 0.25 ml straws and cryopreserved. After storage in liquid nitrogen for 5 d to 7 d, two straws were randomly obtained from each ejaculate and thawed at 38uC for 20 s. These samples were immediately evaluated for the frozen semen quality traits, including post-thaw cryopreserved sperm motility and the percentage of abnormal sperm. The traits were measured using light microscopy, according to the guidelines of the World Health Organization.
Briefly, the ejaculate volume was obtained by dividing the weight of the ejaculate (the difference in weight of a semencollecting vial before and after sampling) by the semen density [7]. The motilities of the fresh and post-thaw cryopreserved sperm were viewed on a TV monitor connected to a camera mounted onto a phase-contrast microscope (Olympus-BX40; Optical Co., Ltd.) at 4006 magnification. A drop of semen was placed onto a pre-warmed (37uC) slide and overlaid with a cover slip. The sperm concentration was determined using a sperm densitometer (Accucell; IMV Biotechnology, L'Aigle, France) that was calibrated using the hemocytometer method. The percentage of viable sperm was calculated by examining more than 100 sperm cells stained by eosin Y-aniline blue from each sample at 4006 magnification [29]. The percentage of sperm deformities were similarly calculated at 4006 and 10006 magnification with Giemsa staining. To minimize sampling error, the quality trait assessment of all semen samples was performed by a single welltrained technician in each semen collection station.
The semen quality traits included the ejaculate volume, initial sperm motility, sperm density, frozen semen motility, and sperm deformity rate. The means and standard errors of the sperm quality traits investigated in the 424 Chinese Holstein bulls are shown in Table 1. The murine Leydig tumor cells (MLTC-1) were obtained from the Cell Culture Collection of the Chinese Academy of Sciences, Shanghai, China.
Genetic variants screening of 59-flanking region DNA was extracted from sperm using a high-concentration-salt protocol and subsequently stored at -20 uC prior to use. One primer pair (TA4 F/R, Table 2) was designed using primer PREMIER 5.0 to clone the 59-flanking region based on the GenBank reference sequence (Accession No. AC_000162.1). The primer pairs were synthesized by Shanghai Sangon Biological Engineering Co., Ltd. The DNA amplification fragments from 50 randomly selected samples were directly sequenced in both directions using an ABI PRISM 3730 DNA analyzer (Applied Biosystems, USA) following standard protocol. The sequencing results were analyzed using the DNASTAR 5.0 package (DNASTAR, Inc., USA) to detect genetic variations in the PLCz gene.
Nucleic acid sequences were analyzed using the accepted software formats.

Genotyping and association study
Sequencing revealed two novel genetic variations (g. 2456 G.A and g. +65 T.C) in the 59-flanking region of the bovine PLCz gene. The gene was subsequently genotyped using the polymerase chain reaction (PCR)-restriction fragment length polymorphism (RFLP). The corresponding restriction endonucleases (HphI and BsenI) were selected to digest PCR products according to the manufacturer's recommendations. The digested fragments were separated by polyacrylamide gel electrophoresis (PAGE) using 10% polyacrylamide gels  [30] of the sample population. The association of the identified genetic variations and haplotype combinations with the semen quality traits was analyzed by the least squares method as applied in the ''PROC GLM'' procedure of SAS software (SAS Institute Inc., Cary, NC, USA), according to the following linear model: where Y ijkl is the observed value of each semen quality trait, m is the overall mean, H i is the fixed effect of genotype or haplotype combinations, P j is the fixed effect of age (j = 2 to 10; classified as: (1) 2 y to 3 y; (2) 4 y to 5 y; (3) 6 y to 10 y); S k is the fixed effect of the origin of bull (k = 1 to 10; classified as: (1) 1 bull from one sire; (2) 2 bulls from one sire; (3) $3 bulls from one sire), M l is the effect of farm, and e ijkl is the random residual error. Values with P,0.05 and P,0.01 were regarded as significant. Multiple comparisons were performed using Duncan's test.

Cloning and construction of PLCz reporter plasmid
Bioinformatics analysis showed that the two genetic variations were located within the proximal region of the promoter. To evaluate the promoter activity of different parts of the 59-flanking region of the PLCz gene, we performed serial truncations of the PLCz promoter fragment, ranging from 22290 bp to +213 bp, and analyzed the activity of the reporter constructs. The four pairs of primers that were progressively located closer to the transcription starting site of the PLCz gene are listed in the Table 2. The forward and reverse primers contained restriction sites for KpnI and XhoI, respectively. The KpnI sites (GGGGTACC) were added to primers as recognition and protection sites. The reverse primer (in the 39 direction) contained an XhoI site (CCGCTCGAG). The amplified promoter fragments were purified, double-digested with the restriction enzymes, and cloned into the pGL3-Basic Luciferase Reporter Vector (Promega, Beijing, China). The recombinant plasmids were confirmed by sequencing, were confirmed by sequencing, namely, pGL3-TA1, -TA2, -TA3, and -TA4.
To examine the effect of different haplotypes on PLCz promoter activity, a series of reporter plasmids with different haplotypes encompassing the TA4 fragments were constructed. The resulting constructs were named pGL3-H1, pGL3-H2, pGL3-H3, and pGL3-H4. Each plasmid harbored a core promoter region with different haplotypes.

Transient transfection and luciferase reporter assays
The MLTC-1 cell line was cultured in RPMI-1640 medium (Sigma Co., St. Louis, MO, USA) with 10% fetal bovine serum (Invitrogen Life Technologies) containing 10 mg/L of penicillin and streptomycin (Invitrogen Life Technologies) at 37uC in 5% CO 2 . For the luciferase reporter assays, MLTC-1 cells were inoculated in 48-well plates and grown to 70% to 80% confluence. Transfection was performed using the Lipofectamine TM 2000 reagent (Invitrogen) according to the manufacturer's instructions. Cells were cotransfected with 50 ng of the pRL-TK vector DNA (Promega) and 400 ng of either the empty pGL3-Basic plasmid (a promoterless control from Promega) or one of the promoter constructs with different lengths and haplotypes of the PLCz promoter. The pRL-TK vector, which provided the constitutive expression of Renilla luciferase, was cotransfected as an internal control to correct for differences in transfection and harvesting efficiency. After 36 h of incubation, cells were harvested and analyzed for luciferase activity using the Dual-Luciferase Reporter Assay System (Promega). Promoter activity was reported in relative light units (RLUs) and normalized against the activity of the empty pGL3-Basic vector. All transfections were performed in triplicate and repeated at least thrice in independent experiments.

Identification of genetic variants within the 59-flanking region of PLCz gene
We sequenced a 1836 bp segment from the 59-flanking region of the PLCz gene in Chinese Holstein bulls. Compared with the PLCz gene sequence (GenBank accession no. NC_007303.4), two novel genetic variations (g. 2456 G.A and g. +65 T.C) were detected (Fig. 1). These two variations have been submitted to the NCBI database (submitted SNP numbers: g. 2456 G.A ss 478894128 and g. +65 T.C ss 478894130).

Prediction of the promoter region of PLCz gene
The bioinformatics analysis predictions indicated that the core promoter region of the PLCz gene is located in the 2461 bp to +100 bp region. Bovine PLCz genes do not have a typical TATA box in the promoter region. However, a putative TATA box with an ATCAATA sequence was present at 2289 bp to 2283 bp from the transcription start site (TSS; the nucleotide sequence numbered +1 is the first A of the TSS) [31]. The putative promoter contained six functional elements, including the AML-1a, CDP CR, P300, GATA-1, SRY, HSF, delta E, and Cdxa transcription factor binding sites (Fig.1).
The SNP g. 2456 G.A was found in the transcription factor binding site of the promoter core region. This SNP was in an activating protein-1 (CAP/AP1) transcription factor binding motif, which was eliminated by the presence of the G allele. In the locus g. +65T.C, the ''C'' allele caused the disappearance of the CdxA transcription factor binding site (Fig. 1). These results were used to investigate the potential effects of the two genetic variations on the regulation of bovine PLCz gene transcription. Genotyping the novel genetic variations by PCR-RFLP Two novel genetic variations g. 2456 G.A and g. +65 T.C were genotyped by PCR-RFLP in the 424 Chinese Holstein bulls. Digestion with HphI of the amplified PLCz gene g. 2456 G.A locus produced fragments of the following sizes: 318 and 40 bp for genotype AA; 358, 40, and 318 bp for genotype GA; and 358 bp for genotype GG (Fig. 2). Digestion with BsenI of the amplified g. +65 T.C locus produced the following fragment sizes: 507 and 145 bp for genotype TT; 652, 507, and 145 bp for genotype TC; and 652 bp for genotype CC (Fig. 2). The obtained genotypes were in agreement with the DNA sequencing results.

Genetic parameter analysis of two genetic variations
The genotypic and allelic frequencies of the g. 2456 G.A and g. +65 T.C loci as well as their genetic diversity are shown in Fig. 3. The A and C alleles are the dominant alleles of g. 2456 G.A and g. +65 T.C, respectively. Their high frequencies might be attributed to the long-term breeding of these populations for different purposes and selections. The PIC, H e , and Ne values as well as the x 2 -test results indicated that g. 2456 G.A was in Hardy-Weinberg equilibrium (x 2 2456 = 3.29, P.0.05) for the analyzed population. Therefore, the selection pressure on this SNP in the population was not powerful. By contrast, g. +65 T.C did not meet the criteria for Hardy-Weinberg equilibrium (x 2 +65 = 12.40, P,0.05).

Association between single variation and semen quality traits in Chinese Holstein bulls
The effects of the two novel genetic variations on semen quality traits (ejaculate volume, sperm density, initial sperm motility, frozen semen motility, and the deformity rate) of Chinese Holstein bulls are summarized in Table 3. g. 2456 G.A was associated with the ejaculate volume (P,0.01). This locus is situated within a putative binding site for the transcription factor CAP/AP1. Bulls with the genotype GG had higher ejaculate volume than those with genotype AA. Similarly, g. +65 T.C was significantly associated with the ejaculate volume in the analyzed population (P,0.05). The presence of the T allele at g. +65 T.C created a CdxA transcription factor binding motif that was eliminated upon its substitution by the C allele. Bulls with the genotype CC had higher ejaculate volume than those with genotypes CT and TT (P,0.05). The g. 2456 G.A and g. +65 T.C loci are within transcription factor binding sites, which implied their correlation with the transcriptional regulation of bovine PLCz via the binding of transcription factors. A number of significant associations were observed between the various haplotype combinations and bovine semen quality traits (Table 4), if the environment and peculiarities in the AI stations were disregarded. The ejaculate volume from bulls with the haplotype combinations H1H1 (P,0.05), H1H2 (P,0.01), H1H4 (P,0.01), H2H2 (P,0.05), H2H4 (P,0.01), and H4H4 (P,0.05) were significantly higher than those with the haplotype combina- tion H3H4. Bulls with the haplotype combination H4H4 demonstrated higher initial sperm motility by 6.15% and 7.23% than the animals with the haplotype combinations H1H3 (P,0.05) and H3H4 (P,0.05), respectively, but had lower motility than those with H2H2 (P,0.05). Bulls with the haplotype combinations H1H2 (P,0.05) and H2H2 (P,0.05) had significantly higher sperm density than the haplotype combination H3H4. Significant associations were not found among the different haplotype combinations for post-thaw cryopreserved sperm motility and the sperm deformity rate.

Promoter activity of the PLCz 59-flanking region
To determine whether the fragment 22290 nt to +111 nt constitutes an active promoter, we amplified a 2401 bp fragment containing the 59-flanking region of PLCz (Fig. 4A). Subsequently, we generated truncated constructs (TA1, TA2, TA3, and TA4) by progressive deletion of nucleotides from the 59-end and cloned  these fragments into the pGL3-Basic Luciferase vector (Fig. 4B). These constructs were transiently transfected into MLTC-1 cells and tested for luciferase activity to determine the shortest required sequence for the transcription of PLCz.
As shown in Fig. 4C, the promoter activity of all the constructs was significantly higher than the basal activity of the control pGL3 vector. Sequential additions caused the 59-end of the 726 bp TA1 fragment to generate different constructs, from TA1 to TA2 (1307 bp), TA3 (2191 bp), and TA4 (2401 bp), thereby decreasing the luciferase activity by approximately 74.7%, 90.8%, and 69.3%, respectively (relative to the shortest construct, TA1). This result indicated that the region between 2641 nt and +112 nt is responsible for most of the promoter activity in the PLCz gene. Therefore, the TA1 fragment of the proximal 59-flanking region of PLCz was recognized as the active core promoter. The TA1 fragment contained both g. 2456 G.A and g. +65 T.C. To investigate the effect of these potential functional genetic variations on the PLCz expression, the 5-flanking region of the PLCz gene containing these loci was subjected to further functional analysis of the promoter activity.

Differential transcriptional activity of different PLCz promoter haplotypes
To assess the effect of the potentially functional variants (Fig. 5A) on the transcriptional activity of the different PLCz promoter haplotypes, we generated four haplotype constructs (Fig. 5B), which were designated as H1 (GT), H2 (GC), H3 (AT), and H4 (AC). Different recombined haplotypes from the PLCz core promoter TA1 were obtained for transfection into MLTC-1 cells. The dual-luciferase activity of the pGL3-Basic vector alone was used as the control. As shown in Fig. 4C, all the constructs provided higher levels of luciferase expression than the basal activity (control) of the empty vector (pGL3) (P,0.05). However, the recombined haplotype constructs differed from each other in terms of their transcriptional activity. Maximal activity was observed for H2 (GC) haplotypes as compared with the others haplotypes. The H2 haplotypes containing the 2456 G and +65 C variants showed 80% and 53% higher transcriptional activity, as compared with the haplotypes H1 (GT; containing the 2456 G and +65 T variants) and H4 (AC; containing the 2456 A and +65 C variants), respectively (P,0.01; Fig. 5C). Given that the region   between 2641 nt and +112 nt represents the core promoter, g. 2456 G.A and g. +65T.C may be functional variants involved in the transcriptional regulation of PLCz. Taken together, our data indicated that the functional g. 2456 G.A and g. +65 T.C variants in the recombined haplotypes played significant roles in the transcriptional activity of the PLCz promoter.

Discussion
The heritability values of the ejaculation volume, sperm concentration, and sperm motility were 0.09, 0.16, and 0.22, respectively [32]. Correlations between breeding values of semen quality traits and routinely estimated breeding values of male fertility were low and ranged from 0.08 to 0.17 [33]. Therefore, single locus effects are expected to be minimal and would require the analysis of larger animal populations [34]. In this study, the g. 2456 G.A SNP in PLCz was in Hardy-Weinberg equilibrium in the population of 424 bulls, thereby implying that the tested population had constant allele and genotype frequencies.
PLCz, which is expressed in mammalian sperm, uniquely possesses all the essential properties of a sperm factor [16]. The bovine PLCz gene is composed of 14 exons and encodes a protein with 634 amino acids; PLCz is a key protein during sperm-egg fusion in mammals [35]. The 13 known mammalian PLCs are classified according to their structure into six types (Beta, Gamma, Delta, Epsilon, Zeta, and Eta), of which PLCz is the smallest [36]. PLCz contains the X and Y catalytic domains, which are associated with the catalytic activity common to all phosphoinositide-specific PLCs, two pairs of EF hand domains that are associated with Ca 2+ binding, the C2 domain that binds to phosphoinositde containing lipids [36], and a putative nuclear localization signal in the linker region between the X and Y domains. Polymorphisms in the 59 UTR of a gene may affect its gene product(s) by altering transcription factor binding or RNA stability. To date, little work has been done to investigate potentially functional genetic variants in the 59-flanking region of the PLCz gene and to characterize their distribution in Chinese Holstein bulls. To obtain more information on the regulatory polymorphisms in the 59-flanking region of the PLCz gene, we conducted a functional analysis of the 59-flanking region of PLCz. In the present study, 59-flanking sequences of the bovine PLCz gene were cloned and sequenced to determine the genetic variation in the population and to predict the core promoter and regulatory elements. Furthermore, the novel genetic variations g. 2456 G.A and g. +65 T.C were identified in the putative promoter region of the bovine PLCz gene.
In bulls, Daghigh-Kia [28] reported a polymorphism at +2749 (G .A) in intron 6 of the PLCz gene. DNA variations in introns and exons of PLCz were previously reported to have an important effect on male infertility [37]. In the present study, two SNPs were found by screening transcription factor binding sites in the 59flanking region of the PLCz locus. Genetic variations in transcription factor binding sites may cause significant potential phenotype diversity [38]. Numerous studies have reported that promoter region polymorphisms are associated with gene expression. For example, a promoter haplotype of the mannan-binding lectin (MBL) gene was associated with the dramatic decrease in the serum concentration of the gene product [39]. Similarly, a genetic variation (g. 2534 T.C) in the bovine peptidase S (PEPS) gene promoter was associated with fat percentage and SCS (somatic cell score) traits [40]. In this study, the g. 2456 G.A locus was in a transcription factor (CAP/AP1) binding site, which was present with the G allele but lost by the presence of the A allele. CAP/ AP1, as a crucial transcription factor, consists of homodimers and heterodimers of jun, fos, and other activating transcription factor proteins. Studies suggest that different AP-1 factors may regulate different target genes and consequently execute distinct biological functions [41]. Therefore, the g. 2456 G.A locus can be considered as a regulator of gene expression via transcription factors.
The g. 2456 G.A and g. +65 T.C loci with genotypes GG and CC, respectively, were significantly associated with bovine ejaculate volume traits. The haplotype combination H2H2 (H2 = GC) likewise presented relatively high ejaculate volume. These findings have possible interpretations. The ejaculate is composed of secretions from various sources. Approximately twothirds of the ejaculate volume is contributed by the seminal vesicles, whereas less than one-third is contributed by the prostate. Furthermore, up to 10% could have originated from the testicle and epididymis, with a small fraction from the bulbourethral glands [42]. Spermatogenesis is a complex process involving spermatogonial stem cells. As a surface protein, the equine PLCz expression is localized in the acrosome, equatorial segment, and head-midpiece junction, as well as the principal piece of the flagellum, in all epididymal, uncapacitated, and capacitated sperms by immunofluorescence [43]. Thus, PLCz might have important functions during spermatogenesis and sperm cell development in the seminiferous tubules and epididymis. Furthermore, the effect of the SNPs g. 2456 G.A and g. +65 T.C, specifically the haplotype combination H2H2, on the ejaculate volume traits requires further elucidation. With the application of molecular genetics, the PLCz gene associated with semen quality traits can be utilized in breeding programs through MAS. Therefore, the H2H2 haplotype may be used as a molecular marker for the selection of bulls with high ejaculate volume.
This study is the first report that analyzes the transcriptional activity of the 59-flanking regions of the bovine PLCz gene in MLTC-1. MLTC-1 was chosen as the cell line for our experiments because the predicted amino acid sequence of bovine PLCz had 69% homology with that of murine PLCz, as revealed by the multiple alignment using the ClustalW algorithm. Furthermore, our results support the notion that each fragment has very high transfection efficiency in MLTC-1 cells. The fragment TA1 had a significantly higher effect on the decreasing promoter activity, which is consistent with predicted results of the bioinformatics analysis. However, the TA2 and TA4 fragments had a significantly higher relative luciferase activity, as compared with the pGL3basis vector. The mechanism for the masking effect of these fragments remains unknown. However, we hypothesize that the complex upstream regulatory elements may interact with each other during the regulation of PLCz gene expression. The identification of other signals or extracellular stimuli that interact with the factors that induce changes in the PLCz promoter activity would be of interest.
Haplotypes, unlike single genetic variations, are more likely to have significant effects on traits [44,45]. In this study, the genetic variations g. 2456 G.A and g. +65 T.C are located at the core promoter region of the PLCz gene, thereby suggesting their important function in bovine PLCz expression. The significant associations of the eight haplotype combinations and the semen quality traits varied among the tested Holstein bulls. Furthermore, the haplotypes H1 (GT), H2 (GC), H3 (AT), and H4 (AC) were transfected into MLTC-1 cell lines to assess their effect on the transcriptional activity of the PLCz promoter. The reporter plasmid with pGL3-H2 had higher relative luciferase activity than those with pGL3-H4 and pGL3-H1 (P,0.001). Bulls with the H2H2 haplotype combination had notably higher ejaculate volume, initial sperm motility, sperm density, and frozen semen motility than the H4H4 and H1H1 haplotype combinations, as well as lower sperm deformity rates than H4H4 and H1H1. A possible explanation for these results is the localization of PLCz expression. PLCz expression was found to be localized in the postacrosomal region of murine, bovine, and human sperm [46]. The localization of PLCz in the post-acrosomal region at the head of the sperm is important because the head is rapidly exposed to the ooplasm after gamete fusion; PLCz is reported to have an oocytepenetrating effect in the sperm acrosome for successful fertilization [47]. Furthermore, PLCz mRNA was significantly expressed in the testis and epididymis of G-I boars as compared with G-II boars; G-I boars are noted for their relatively better sperm quality [48]. The association analysis of the four haplotypes confirmed the observed luciferase activity. The different haplotypes exhibited different promoter activities, which implied that genetic variants are likely to regulate PLCz expression and consequently influence the physiological function of PLCz on semen quality traits. For example, promoter region polymorphisms in the human MBL2 gene were found to control the baseline expression of MBL2 [39]. The haplotypes HY, LY, and LX were correlated with high, intermediate, and low MBL levels, respectively [39,49]. Therefore, future studies should identify natural genetic structures that could be changed or utilize artificial selection in Chinese Holstein bulls under various environmental or economic conditions. Thus, the frequency of favorable alleles in the population would be expanded, thereby gradually fixing the dominant alleles for highquality sperm characteristics in future generations of Chinese Holstein cattle breeds.
Our findings suggest that the functional genetic variants of PLCz and their respective haplotypes modulate the transcriptional activity of the PLCz promoter and may contribute to the semen quality traits of Chinese Holstein bulls. Moreover, promoter fragments with different haplotypes exhibited different promoter activities, which imply that these genetic variants are likely to regulate PLCz expression and consequently influence physiological functions of the gene.