Sperm morphometric subpopulations and hetero morphism in Bos taurus indicus bull semen and their influence on in vitro fertilization were investigated. Area, center of mass and width obtained on 200 digitized cells/animal (n=9) were used to determine the subpopulations. The analysis of Ward’s grouping with Euclidean distance was used to determine the subpopulations, and number of subpopulations was defined according to the behavior of distance level and similarity. Fertility was evaluated by Pro-Nuclei Formation Rate (PFR) after In Vitro Fertilization (IVF). Sperm hetero morphism in one hundred cells considered to be morphologically normal was assessed by Harmonic Fourier Descriptors and Principal Component Analysis (PCA). Four subpopulations (SP1 to SP4) were determined, exhibiting asymmetrical distribution. SP4 was composed by 43.0% of the cells that showed the greatest area (39.3 ± 1.8 µm 2 ), while SP1 (35.3±2.4 µm 2 ), SP2 (35.9 ± 2.5 µm 2 ), and SP3 (36.9 ± 3.0 µm 2 ) presented similar distribution (21.0%; 17.5%; 18.5%; respectively). A negative correlation (r= -0.86; P=0.01) was observed between the SP4 and PFR. Hetero morphism rate did not influence PFR. In conclusion, semen of Bos taurus indicus bulls exhibited morphological heterogeneous structure and cells exhibiting large head area presented lower capacity of IVF.

Introduction

Methodologies of sperm cell evaluation have experienced advances in the last years, but it is still not possible to precisely predict the fertility of a semen straw [1]. One reason for this constraint should be the attempt to consider the ejaculate as a homogeneous solution of spermatozoa, but cells present distinct patterns of movement and morphometry [1,2] and different clusters defined by these parameters were described in various species [3-9]. Besides, the use of mean values of seminal characteristics to correlate with fertility may be no representative [2], and currently, a subpopulation structure based on principal components and cluster analysis has been introduced [10].

The semen of the bull exhibits a heterospermic template. Sperm morphometry analysis has shown that variations within and among individuals influence fertility of the ejaculates [11]. Studies have investigated the phenotype of the sperm head and its possible relation to fertility [12,13], but those evidencing influence of morphological sperm subpopulations and fertility are speculative [10].

Fertility may be estimated in vivo or by the capacity of spermatozoa to penetrate oocytes in vitro [14-16]. However, only one study has addressed the correlation between the influence of subpopulations and fertility [3]. The present study aimed to identify sperm subpopulations and morphological parameters in frozen semen of Bos taurus indicus bulls related to in vitro fertility.

Material and Methods

Experiment 1: Determining Sperm Subpopulations and Semen Hetero morphism

Semen samples

One semen straw of nine sexually mature bulls of the Nelore breed (11.0 ± 3.1 years old) were used. Semen was collected via artificial vagina and frozen according standard protocols in an artificial insemination center. The experiment was conducted in the Veterinary Department - UFLA, and EMBRAPA, Brazil, being approved by the ethics committee on the use of animals (CEUA-UFLA 073/15).

Semen processing

Straws were thawed (37 °C/30s) and motility was assessed using computer-assisted sperm analysis (IVOS-Ultimate 12’s, Hamilton Thorne Biosciences, Beverly, MA, USA), configured for bovine spermatozoa: number of frames: 30; frames/s: 60Hz; minimum contrast: 50; minimum cell size: 6 pixels; contrast to static cells: 30; straightness: 60%; average path velocity cutoff: 30 µm/s; minimum average path velocity: 40 μm/s; straight-line velocity cutoff: 20 μm/s; cell intensity: 80; static head size: 0.23-1.91; static head intensity: 0.56-1.20; static elongation: 8-92; magnification: X 1.89 and temperature: 37 °C. For the evaluation, 9 μL of semen diluted in TALP at a concentration of 25´106 cells/mL were placed in a chamber (Makler, Santa Ana, CA, USA). Three fields were selected for analysis. Parameters evaluated were Total Motility (TM), Progressive Motility (PM) and Linearity (LIN) in %; Average Path Velocity (VAP); Straight Line Velocity (VSL) and Curvilinear Velocity (VCL) in micrometers/second; Lateral Head Amplitude (ALH) in micrometers and Beat Cross Frequency (BCF) in Hz.

A semen aliquot of 50 µL was added to a tube containing 100 µL of formaldehyde saline solution. Ten microliters were examined and 200 cells considered to be morphologically normal per straw were selected. Morphology was assessed by light microscopy, Harmonic Fourier Descriptors [17] and Principal Component Analysis [18], using 100 cells per bull in the latter analyses. A phase contrast microscope (Olympus BX41) with a green filter (IF550) coupled to a camera (Charge Coupled Detector - SDC-415, Samsung) was used to acquire the images, stored as .bmp files. Spermatozoa were centered on the screen at 400X magnification and it was assured that the cells were in the same focus plane with none artifacts.

Image processing

Sperm images analyses were performed using the Image J platform (National Institutes of Health, USA). Images were transformed to 8-bit and digitized using limits around 65-76 in gray scale (Threshold tool). After adjustments, the spermatozoon head was cut and transferred to an image with white background. The head contour line was identified (Wand tool) creating a Region of Interest (ROI) to be analyzed. Parameters evaluated were area (µm2), centroid (X and Y), center of mass (XM and YM), perimeter, rectangle, ellipse axis, width, and length (µm). To avoid redundancy [19], after statistical analysis, only three measurements were used to determine the subpopulations.

Determination of subpopulations

Sperm subpopulations were defined by multivariate analysis of hierarchical and nonhierarchical grouping. Samples were evaluated independently and together. The multivariate cluster analysis K-means was performed to classify spermatozoa in optimum reduced number of subpopulations, according to the head measurements, after applying the Ward’s grouping analysis, using Euclidean distance. Cells with similar morphometry were clustered in the same group and different clusters were formed according to the observed mean values of each group. Subpopulation numbers was based on the previous analysis of hierarchical dendrogram and calculated on the estimate of the level of distance and similarity among groups.

Experiment 2: Fertility Evaluation

Semen samples

Semen of three out of the nine bulls were excluded from in vitro fertility evaluation, because significant differences (P<0.05) for area after categorization in clusters were not observed, total motility was below 30% and because one was considered an outlier in the residual analysis.

Bovine oocytes collection and maturation

All reagents used in IVF were obtained from Sigma-Aldrich Química S.A (São Paulo, SP, Brazil).

Ovaries from slaughterhouses were transported in saline solution with 100 μg/mL of streptomycin and 100 UI/mL of penicillin G at 37 °C. Cumulus-Oocyte Complexes (COCs) aspirated from follicles from 3 to 8 mm were deposited in tubes of 50 mL for sedimentation. The supernatant was removed, precipitate was poured in a 100´200 mm cell culture plate containing TALP-HEPES for COCs evaluation. COCs with three or more layers of cumulus cells with homogeneous cytoplasm were selected for In Vitro Maturation (IVM). Groups of 25-30 COCs were placed in 400 μL of medium for IVM compounded by TCM 199 with Earle salts, L-glutamine, 2.2 g/L of sodium bicarbonate, 10% of estrus cow serum, 20 mg/mL of follicle stimulating hormone, 49.4 mg/mL of sodium pyruvate (0.45 mM), 100 UI/mL of penicillin and 0.1mg/mL of streptomycin in four-well plates and cultivated in a chamber during 22-24h at 38.5 °C, with 5% of CO 2.

Preparing semen and in vitro fertilization

After maturation, COCs were washed in Phosphate-Buffered Saline (PBS) and transferred into drops of 70 μL of the fertilization medium (FERT-TALP plus heparin, 10 UI/mL). One semen straw of the same ejaculate of each bull used to determine the subpopulations was used for IVF. Straws were thawed at 37 °C/30s and semen motility was evaluated under light microscopy.

A semen aliquot of 250 μL was centrifuged during 7min at 300´g at room temperature. Supernatant was removed and the sediment was re-suspended in 400 μL of medium FERT-TALP and centrifuged again (5min; 300´g). Supernatant was removed and the sediment was re-suspended in 70 μL of medium FERT-TALP. Motility and vigor were evaluated once more. Spermatozoa were counted in a Neubauer chamber and diluted in fertilization medium for a final concentration of 2´106 sperm cells/mL. A semen aliquot of 30 μL was added to each drop of fertilization, achieving a volume of 100 μL. Plates were incubated at 38.5 ºC/20h in chambers with 5% of CO 2.

Evaluation of pronuclei formation rate

Possible zygotes were transferred to 1 mL of TALP 20h after insemination and stirred during 5min to remove cumulus cells, being washed in PBS with polyvinyl alcohol (PVA; 1 mg/mL), fixed during 30min in paraformaldehyde 4% and washed again in PBS-PVA. Specimens were stained with Hoechst 33342 (5 µg/mL), placed under a coverslip and PFR was evaluated [3] under 400X magnification in epifluorescence microscopy (Motic BA400, excitation of 330-385 nm and emission of 420-490 nm).

Statistical Analyses

Principal Component Analysis (PCA) was done using scripts of the Matlab software version R2012 (8.0.783). Statistical analyses of subpopulations were performed in software Minitab 17.0.1. Regression analyses and residues were made in R software version 3.4.2, package MVar.pt version 1.9.9. Data previously standardized because metrics were at different scales were tested regarding normality of residues by the Shapiro-Wilk test. According to Martinez-Pastor, et al. [19], data were submitted to variance group analysis and one parameter of each formed group was selected, and selected parameters were submitted to Pearson’s correlation test. Only those presenting r values <0.3 and P<0.05 were used.

Hierarchical and nonhierarchical methods were used for grouping 1800 spermatozoa from one straw per bull evaluated together. For each bull 200 cells were analyzed. A multivariate analysis of Ward’s hierarchical grouping based on Euclidean distance was performed to create a hierarchical alignment cluster and then to obtain a dendrogram. The optimum number of clusters was determined by the behavior of similarity level and the distance in the steps of grouping. During the analysis to establish the groups, similarity among conglomerates decreases and distance increases. Sharp variations of these measurements indicate the ideal number of clusters. After defining the number of groups, the K-means method was applied by a nonhierarchical multivariate analysis to determine the proportion of cells belonging to each subpopulation. Validation of groups to verify the mean values difference among each subpopulation was made by the Tukey test.

Each straw was considered one independent observation. Motility at the time of IVF, PFR and proportion of subpopulations were tested regarding normality, homoscedasticity, and independence by means of the Shapiro-Wilk, Bartlett and Box-Pierce tests, respectively. The relationship between the proportion of each subpopulation in the semen samples and IVF was evaluated by regression analysis according to a previous work [3], as well as the relationship between sperm motility before and at the time of IVF and PFR. Differences were indicated by a probability of P < 0.05.

Results

To determine subpopulations, using the variance grouping analysis and Pearson’s correlation, it was selected area, center of mass and width. These variables presented the lower correlation values (r<0.3; P<0.05).

Four different sperm subpopulations were defined on the basis of the difference in the index obtained by grouping similarity and distance levels (Figure 1). When the cells of all animals were analyzed together, spermatozoa head area presented the greatest difference among the four subpopulations (P<0.05). Subpopulation 1 was composed by sperm cells with smaller area (35.3 2.4 µm2), as well as the SP2 (35.9 2.5 µm2), representing 21.3% and 17.5% of the total cells, respectively. Spermatozoa of these two subpopulations were defined as short and tapered. SP3 was composed by average spermatozoa head area of 36.9 3.0 µm2 and has 18.5% of the total cells, being defined as short and wider. SP4 was composed by the largest areas (39.3 1.8 µm2), representing 42.7% of the total of cells, being defined as large and elongated (Table 1).

SP4 presented the smallest variations of the center of mass, measured by Coefficients of Variation (CVs), while SP2 showed the greatest. Distribution of subpopulations among all animals was asymmetrical. Four out of nine animals presented normal distribution among the subpopulations. SP1 exhibited normal distribution in six animals, SP2 in seven and SP3 and SP4 in eight of the nine bulls. Six bulls had SP4 as the most representative cluster, SP2 in two and SP1 in one individual (Table 2).

Morphology was evaluated in the nine bulls by PCA and the first two components revealed a variance ranging from 45.5 to 95.6% (Figure 2). No correlation was observed between hetero morphism rate and PFR.

After processing for IVF, all straws had motility above 70% (Table 3). The fertility rate of six straws used to IVF ranged from 54.4% to 64.7% and a negative correlation (r= -0.86, P<0.05) was observed between the number of cells in SP4 and PFR. The regression equation (PFR= 70.8-0.2893SP4) showed R2=0.73 and the model demonstrated that when the relative proportion of SP4 cells increased, PFR decreased. There was no significant correlation between sperm motility before (total and progressive) and at the time of IVF and PFR (R2=0.14, Figure 3). Negative correlations between PFR and VCL (r=-0.27), VAP (r=-0.35) and VSL (r=-0.50) were observed (P<0.05).

Discussion

As it was observed previously in Bos taurus taurus [1,2,10], frozen semen of Bos taurus indicus also exhibits morphological heterogeneity and four sperm subpopulations defined by different head area were verified in our study. In Bos taurus taurus three to four subpopulations based on morphometric head parameters were observed [8,10,20] and four subpopulations were identified by kinetic attributes [3,5,6]. According to our results, the number of subpopulations determined by sperm head morphometric parameters ranges from three to four between the two subspecies.

Most studies on subpopulations showed a different cluster distribution among the animals involved [4,10,21], like data of the present study, and differences among ejaculates of the same individual [10]. These differences are considered to be a result of genetic and endocrine changes, related to adaptations to the fertilization environment [2,22], that should influence the fertility of the animal.

The average head area of the cells observed in this experiment, varying from 35.3 to 39.3 µm2, is comparable to other data [12,13,23,24]. Our results demonstrated that the higher proportion of SP4 decreases PFR. Estimation of semen fertility involves several factors [25], and IVF has been used as an indicator for it [3,26]. One study showed that motility could be one factor responsible for PFR [3], and positive correlations among rapid and progressive cells and the number of spermatozoa linked to the pellucid zone (r=0.79), penetration rate (r=0.66) and PFR (r=0.63) were verified. In cervids [27], high fertility males have ejaculates with high percentages of spermatozoa exhibiting fast and linear movements. Morphologically, these patterns of motility were correlated with a large proportion of spermatozoa having small and elongated heads.

In the present experiment, semen samples with total motility <30% were removed from IVF and statistical analysis. Variations in the relative proportions of cell clusters did not interfere with semen kinetics. No correlation was verified between motility ranging from 39.0% to 69.0% (≥70% after Percoll treatment) and PFR. Moreover, and negative correlations between PFR and VCL, VAP and VSL were observed. These results suggest that large and wide cells were responsible for the differences observed in fertility among the bulls. However, in the sheep, a positive correlation between fertility and proportion of wide and long cells was observed [28], while a negative correlation was observed between fertility and percentage of spermatozoa with short and wide heads [29].

Some efforts were made to investigate sperm cell morphometric parameters and fertility [12,13,30]. In the filtered and extended semen of the stallion, average measurements of length, area and perimeter of the spermatozoa head were higher in sub fertile animals [30]. In the frozen semen of the bull, no differences were observed among the mean values of length, width, area and perimeter, and the non-return to estrus rate [13]. Nonetheless, these studies used mean values of sperm head variables and heterogeneity of spermatozoa was not considered. In the fresh semen of highly fertile bulls it was observed that small variations in size of the sperm head can exist not interfering with the fertilization capacity of the sperm cells [12].

Sperm morphology [31] was assessed comparing frozen semen of bulls divided in high and low fertility groups. Bulls of high-fertility have sperm nuclear shape more elongated and tapered than those of lower fertility. Our approach consisted in cluster categorization and semen in vitro fertility evaluated by PFR, and a negative correlation between cells allocated in SP4 (long and wide) and IVF was observed. High fertility rates among rams [32] were associated to a subpopulation with short and elongated heads, while no relationship between fertility and average values of sperm head dimensions were found. Thus, cluster analysis is relevant in semen evaluation.

Morphology of sperm cells is sensible to physical [23] and physiological stress [24]. PCA in four bulls showed that the first two components explained more than 81.3% of the variance among sperm cells, thus cell morphology in these ejaculates was very similar. In the other five individuals, this variation ranged from 45.5 to 74.5%, which implies a greater sperm cell hetero morphism. In addition, the center of mass CVs indicated that some cell clusters are morphologically more stable than the others. No relationship was observed between hetero morphism rate and PFR. Number of cells showing normal morphology and morphometry increase after Percoll treatment [33], but how these parameters influence the fertilization process and if only the sperm dimensions are involved remain uncertainty. Nonetheless, sperm head morphometry appears to be an independent semen parameter and normal morphology is not necessarily an indicative of morphometric attributes.

The association found in the present study between poor fertility and large sperm head size may be related to chromatin condensation, as it was observed previously [34,35]. Bovine spermatozoa with chromatin condensation abnormalities presented alterations that include larger head area [31,36,37] and bulls that possess spermatozoa with high levels of DNA fragmentation and poor protamination also have lower fertility [35]. Lower rates of PFR verified in this experiment found in semen with more large-head cells could be a consequence of inadequate protamination, resulting in sperm heads exhibiting larger areas.

Conclusion

In conclusion, heterogeneity among sperm cell clusters must be considered in an accurate semen analysis. The frozen semen of Bos taurus indicus bulls is composed by four sperm cell subpopulations, defined by head morphometric parameters. Morphological subpopulation structure and distribution is different among individuals, and a high relative proportion of cells in a cluster exhibiting large sperm head area is related to poor in vitro fertility. The heterospermic pattern in cells considered to be normal in light microscopy observed in the ejaculate does not influence fertility evaluated by pro-nuclei formation rate.