CC BY
32
UDC 57.08:636:31 doi: 10.15407/biotech9.02.070
POSSIBILITY OF COMPUTER EXPERIMENT IN STUDY OF ANIMAL SPERMATOZOA HETEROGENEITY
L. V. Gorbunov National Technical University
Y. M. Mazharova «Kharkiv Polytechnic Institute», Ukraine
E-mail: 4glv@i.ua
Received 01.03.2016
A simulation model for evaluating the survival and fertilizing capacity of animal spermatozoa was developed, taking into account the initial condition of the sperm and the effectiveness of cryopreservation stages. The model is based on an analytical expression that reflects the main reasons for the survival of reproductive cells in onto-, techno- and phylogenesis. The decrease in spermatozoa resistance depends on a number of biological factors — the animal species, physiological conditions of sperm donor and recipient, the ejaculate quality, and technological factors — the effectiveness of the methods of cell cryopreservation and egg insemination. The discrepancy between the results of cell motility obtained by calculation and experimental methods amounted to less than 2% as a result of our own experiments and to less than 5% for the data taken from literature. A feature of the model is the complete independence of the effectiveness of studied techniques from the heterogeneity of animal sperm.
The conducted computer experiment showed that the difference between the values of initial motility and fertilizing capacity of sperm varies from 50 to 100% depending on the difference of biological parameters, while the index of the effectiveness of selected technique creates an error of about 1%. Comparative analysis of alternative technologies of spermatozoa cryopreservation showed the maximum efficiency of the stages of cryoprotectant use, freeze mode, survival and fertilizing capacity of the object. The use of computer modeling allows to greatly reduce the spread in spermatozoa preservation values that were obtained in different experiments, and thus to reduce the time and costs it takes to obtain reliable results.
Key words: system biology, cryopreservation, heterogeneity, spermatozoa.
The use of mathematical models when conducting experiments makes it possible to not only detect, but also to explain the obtained regularities, which allows scientists to create entirely new ways for cryopreservation of biological objects. Modern scientific methodology identifies three stages of knowledge mathemization: statistical processing of empirical data, experiment modeling on the basis of regression equations and the existence of relatively complete mathematical theories represented by analytical expressions [1]. Achievements in the field of cryobiology allow to widely use analytical methods in mathematical modeling and to explore the potential of the use of statistical methods for the analysis and optimization of multi-factorial experiments.
To date, there have been a number of studies aimed at improving the motility of animal sperm by optimizing the parameters
that define the freezing mode of biological objects [2-10]. It is found that the success of cryopreservation of biological objects, including sperm, mainly depends on three tied components: proper selection of cryoprotectant and its concentration, composition of a diluent and the freeze-thaw mode characteristics. Analysis of the cryopreservation results [2, 3, 6-8, 10-15] shows that bull and carp sperm motility rate lowers during sperm freezing and thawing.
The result of animal sperm cryopreservation depends on the variation of the initial state of the ejaculate. Consequently, the initial state of the biological object and features of technological operations that were conducted prior to cryopreservation may have an impact on the mobility of deconserved material. Therefore, to study patterns of repeatability of cryopreservation results, it is necessary to develop a mathematical model
for evaluating the effectiveness of the chosen method of cryopreservation that takes into account the initial state of biological objects and the features of different stages of the chosen cryopreservation method.
The aim of the study is to create a simulation model to describe the patterns of influence of animal heterogeneity and the efficiency of steps taken during semen cryopreservation on the fertilizing capacity of sperm.
Materials and Methods
Sperm of breeding bulls of black-and-white and Simmental breeds was used when conducting experiments. Breeding bulls from the experimental-methodological breeding enterprise of the 1st category of the Department Biotechnology at the Institute of Animal Science of the National Academy of Agrarian Sciences of Ukraine (NAAS) were used. When conducting computer simulation, we used the data obtained from the semen of one bull to ensure the homogeneity of the biomaterial. All technological methods for bull semen receiving, evaluation and dilution were performed in accordance with [9, 16]. Fresh sperm with a motility index not lower than 70% was used. The sperm was stored at a temperature of 30 °C no longer than 30 min after it has been obtained.
Upon receipt and evaluation semen was diluted at 1:1 ration with sterile lactose-vitelline-glycerol (LVG) medium (11% solution of lactose — 62 ml, yolk — 30 ml, glycerol — 7 ml). Then it was diluted again with lactose-citrate-glycerol (LCG) (lactose — 6 g, sodium citrate — 1.4 g, glycerol — 5 ml, distilled water — 100 ml) to a concentration of spermatozoa of 15 million/ml. Then sperm was placed into cone-shaped containers of 0.25 ml in volume. Sperm equilibration was performed at 2-4 °C for 3-4 hours with constant stirring.
Sperm of Nivky scaly carp was obtained from sexually mature males (Cyprinus carpio). The fish were caught in a pond during carp spawning at an ambient temperature of 20 °C. Fish that matured naturally without hormonal treatment was used for the experiment [12, 17]. Sperm was obtained separately from each male at 5th stage of maturity, which is important for achieving high cryoresistance. Carp sperm was obtained, assessed and prepared for freezing according to [11, 12].
Cryoprotectant was prepared by dissolving 15 g of Tris-HCl-buffer in 1 l of distilled water, adding hydrochloric acid to pH 8.
To make 1 l of cryoprotectant, 7,3 g NaCl, 0,04 g KCl, 2,3 g NaHCO3, 1,15 g sucrose, 0,52 g MgSO4 and 0,15 g CaCl2 were alternately dissolved in 600-700 ml of Tris-HCl-buffer.
150 g of egg yolk and 100-120 g of ethylene glycol was added to the resulting solution 30 minutes before using it. Then Tris-HCl-buffer was added to the solution to get the total volume of 1 l [11, 12].
Before the experiment, all the dishes, sperm, cryoprotectant and activators (pond water) were cooled to 15 °C. Cryoprotectant was added to the sperm at 1:1 volume ratio while constantly stirring. Then the diluted carp sperm was cooled at the speed of 1-2 °C/min to a temperature of 2-3 °C and equilibrated for 1 hour. Fresh sperm with a mobility index no lower than 70% was used in experiments. Selected ejaculates were placed in cone-shaped containers with tapered initiation of crystallization, sealed and placed in a refrigerator at 2-3 °C for further testing in various freeze modes. Each container was 0.25 ml in volume.
Ejaculate volume, motility and concentration of cells were determined before the experiment using a photoelectro-colorimeter. Motility index was used to assess deconserved carp sperm, motility, survivability and fertilizing capacity indexes were used to assess the condition of deconserved bull sperm. Fertilizing capacity was determined by first insemination of Simmental heifers.
Sperm motility of carp was determined as follows. A drop of water was placed on a glass slide with a pipette, then a small amount of sperm on the tip of a needle was added. Total duration of these manipulations did not exceed 10 seconds. When determining the motility of bovine spermatozoa, biomaterial was applied on a glass warmed to 38 °C, under a microscope. The motility of sperm in each sample was determined at least three times and the average result was calculated.
Microscopes of MBI or MBR brands with magnification of 300-600 were used during the experiment. Bull and carp sperm motility was determined by the ratio of spermatozoa with rectilinear translational motion to the total number of sperm in the field of view of the microscope, expressed as a percentage (Fig. 1).
Sperm cryopreservation was carried out in devices that were developed by us earlier. These devises use passive cooling of the fuser in the neck of Dewar vessels, as a basic process [18]. To make ultrahigh freeze-thaw speeds possible, envelopes made of polyethylene
and aluminum foil were used. Before the actual freezing, the device is cooled in liquid nitrogen. During thawing, it was heated in a water bath to 40 °C. Freezing (thawing) rate in this case was 15 000 (12 000) °C/min and 24 000 (22 200) °C/min for envelopes made of polyethylene and aluminum foil, respectively [18].
The condition of a biological object is characterized by several criteria: the motility of the sperm determines whether the object is suitable for further use in the study; it is evaluated by morphological parameters of the object; fertilizing capacity characterizes the probability of successful development of a biological object further in the experiment; the effectiveness of cryopreservation stage gives an indication of how the condition of a biological object changes as a result of the cryopreservation stage alone; reproducibility characterizes the dispersion of the results obtained by using a particular method of cryopreservation; comparability enables us to compare results obtained using different methods of cryopreservation.
To improve the reproducibility of obtained results, the study was performed using ejaculates of the same quality in each group, while determining the motility and fertilizing capacity of sperm. Animal sperm motility — the S f(1) is defined as the ratio of the number of moving sperm — n, divided by the total number of sperm in the sample — n0:
S = n/no.
(1)
Survival rate was also evaluated for deconserved bull sperm 38 °C. Survival rate was measured as the amount of time during which all spermatozoa lost their rectilinear translational movement.
The initial sperm motility S0 was determined experimentally, after semen was obtained and diluted — S f(1).
Fertilizing capacity of sperm — P f(2) was calculated as a ratio of the number of pregnancies — n to the total number of fertilized cows (for fish — embryos) — n0:
P = n/no,
(2)
where n0 — the initial number of samples used in the experiment.
Insemination of eggs and incubation was carried out at a temperature of 18-20 °C. To eliminate the influence of differences in fertilization and cultivation methods for fresh and deconserved sperm on evaluation
of effectiveness of cryopreservation of carp sperm, the following experiments were conducted. 0.02 ml of fresh diluted sperm was added to each Petri dish in control group, and 0.02 ml of deconserved sperm was added to each Petri dish in experiment group. Semen (0.02 ml) was mixed with eggs (50-60 eggs) and activator (0.25 ml 0,3% NaHCO3) in Petri dishes (10 cm in diameter), stirred thoroughly with a feather for 2 min and quenched with water [19, 20]. When assessing the fertilizing capacity of fresh and deconserved sperm, ejaculate of the same animal was used for 1 h after freezing. Experiments were conducted in the laboratory; inseminated eggs were incubated in Petri dishes. Embryos that were stuck to the bottom of the cup embryos were left to incubate at 18-20 °C, water was changed every 3-4 hours. To remove the sticky layer from the surface of eggs, pond water was added to eggs 5 minutes after they were mixed with activated sperm. After that, eggs and water were stirred for 30-40 minutes to remove the adhesive layer from egg surface. Fertilization percentage was determined by 4-cell blastomere stage, and calculated using the formula:
Pc = .100%, no
Pd = npL .100%, no
(2 а) (2 b)
where Pc and Pd — probabilities of fertilization using native and deconserved sperm; npc, npd and n0 — number of cells that were fertilized with fresh, deconserved sperm, and total number of eggs in the sample.
Statistical analysis of the results obtained in the experiment was carried out using standard formulas for qualitative and quantitative analysis [21]. Standard and developed by us software was used for the calculations. Reliability of differences in control and experiment groups was assessed using Student's t-test by comparing the average values of motility, viability and efficiency that were obtained for the experiment and control group. Differences in conjugate parameters for control and experiment groups were averaged to take into account individual characteristics of biological objects [21]. Reproducibility of cryopreservation results was determined using the value of the standard deviation and the coefficient of variation C,.
Results and Discussion
The need to create a simulation model lies in the need to identify values that cannot be identified directly when conducting a physical experiment. We have proposed a method for calculating parameters of sperm resistance for sperm that was obtained from a particular donor and fertilized recipient (Fig. 1, rd and rp), and for determining the efficiency of technologies that were used (Wd and Wp). These figures summarize the impact of biological (F, A, U, H) and technological factors (a, e, z, T). Calculated parameters (rd, rp, Wd, Wp) determine the conditions of the particular experiment. The model allows to calculate the fertilizing capacity of sperm cells when conducting a computer experiment using the given parameters (S0, rd, rp, Wd, Wp). Simulation modeling ensures comparability of empirical data (S, P) that was obtained under different conditions of the study (S, rd, rp, Wd, Wp). Calculated efficiency indexes of different technological stages (a, t, z, T) are independent from the condition of each biological object (rd and rp), which greatly improves the reproducibility of the results. The model allows one to determine the heterogeneity of biological characteristics (S0, S, A, rd, rp, F, P, U, H), given the fact that each of them depends on many interrelated parameters.
A logistic function was used as a basis for the proposed mathematical model to describe spermatzoa resistance when obtaining sperm, cryopreserving it and inseminating the recipient:
— = r • S(1 - S/K),
dt '
(3)
where S — spermatozoa motility, calculated in terms of f(1);
t — time, expressed in hours; r — coefficient that reflects the rate at which spermatozoa condition changes depending on the physiological condition of the donor;
K — coefficient that corresponds to the maximum value of sperm motility.
The solution of the differential equation f(3) gives a numerical value to describe the survival rate of sperm cells at physiological temperature:
S (t ) =
S0 • K • exp(r • t) K - S0 + S0 • exp(r • t)
(4)
where S0 — the value that corresponds to the value of the initial sperm motility.
The coefficient, which reflects the rate at which the condition of spermatozoa changes, S(t) in this analytical expression, is determined using computer simulation by analyzing the values of sperm motility f(1) that were obtained in the experiment.
According to Bayes' theorem, motility of deconserved spermatozoa Sd(t) can be expressed as the sum of conditionally independent probabilities (expressed in units) that correspond to fresh sperm motility S0(t) minus the probability of the death of some
Fig. 1. Chart of computer study of spermatozoa heterogeneity in cryopreservation and artificial insemination:
S0 , Sd and Р0 , Pd — motility and fertilizing capacity of fresh, deconserved sperm; rd and r„ — resistance of the donor and of fertilized recipients to sperm; Wd and Wp — efficiency of sperm cryopreservation and artificial insemination techniques
cells. This value depends on the effectiveness of the cryopreservation technique that was
used Wd:
Sd(t) = S0(t)-S0 • (1 - Wd).
(5)
In order to improve the accuracy of the analytical dependence S(t) f(4), we created a regression equation that connects the change in spermatozoa motility S to time. The method of least squares was used for this (Fig. 1-5). The values of cryopreservation efficiency Wd were calculated using the proposed model f(3-5) on the basis of spermatozoa motility data obtained in the experiment S f(1).
The error in determining empirically derived sperm motility index S f(1) is determined by assessing the state of a biological object of specified quality. The error is 5% (0.5 points) when using the visual method, and 1% (0.1 points) when using computer video fixation [18]. When using the conventional method for evaluating cryopreservation that takes into account the average value of sperm motility, depending on the initial spermatozoa condition, in different experiments the difference may be more than 10%. The need for such a transformation is due to the fact that average values do not carry information about the condition of each individual ejaculate, while the proposed model S(t) f(5), in conjunction with the empirically derived function S f(1) reflect the individual characteristics of a particular ejaculate.
Fertilizing capacity of fresh sperm of given quality P(S) depends on their motility — S f(1), on the physiological condition of the recipient — r f(7), and on the value that reflects the effectiveness of the technique used for artificial insemination of cows Wp f(8):
P( S ) =
Pm • K • exp(r • S) K - Pm + Pm • exp(r • S y
(7)
where K — coefficient that reflects the maximum value of egg fertilization, and the minimum — Pm, respectively;
r — coefficient that corresponds to the physiological condition of the recipient. It determines the probability of fertilization.
Analytical expressions S(t) and P(S) are derived from f(4) and f(7), using the regression dependence of sperm motility values that were obtained experimentally during sperm cultivation f(1) and insemination f(2) (Fig. 2-6).
The probability of fertilizing capacity of fresh and deconserved sperm P(S) is expressed through the product of dependent parameters of the fertilizing capacity of fresh sperm P0(S) f(7), the effectiveness of cryopreservation techniques Wd and artificial insemination techniques Wp, expressed in units:
P(S) = Po(S)•Wd •Wp
(8)
To improve the accuracy of determining the cryopreservation efficiency by reducing the total variance [18] this value is represented the product of the efficiency of the following technological stages: collection and dilution of the ejaculate Wa, equilibration in cryoprotectant and temperature adaptation We, selection of the freeze-thaw mode W2:
wd = wa w W.
(9)
The following system of equations is used to calculate the parameters of effectiveness:
Sa (t ) = So(t ) •Wa ; Sae (t ) = So(t ) •Wa .ffe
Sd (t ) = So(t ) • Wd.
(10)
Software was developed in Microsoft Excel in order to simplify the procedure for calculating the main biotechnological parameters. A graphical representation of calculated results is shown in Fig. 2-6.
To test the proposed simulation model f(3-10), we analyzed the data for fresh and deconserved bull sperm motility at 38 °C that was obtained at the Livestock Institute of National Academy of Agrarian Sciences of Ukraine [10]. Spermatozoa were frozen in plastic containers (V = 0.75 ml) in a device that was developed by us [15].
It was established experimentally [15], that when frozen with crystallization initiation, the motility index of deconserved bull sperm increased — from 38.6 ± 1.83% to 52.3 ± 1.64%, respectively. Survival time of deconserved bull sperm that was frozen with the initiation of crystal growth was 9 hours, and without the initiation of crystal growth — 7 hours (Fig. 2).
Regression dependence of experimental f(1) and model f(4, 5) values of spermatozoa motility in time allowed to determine approximation coefficients for fresh diluted
Fig. 2. Dependence of bull sperm motility from cultivation time at a temperature of 38 °С for fresh diluted cells — S0(t), deconserved cells in the first group — Sd1(t) and in the second — Sd2(t), respectively:
The effectiveness of cryopreservation techniques amounted to Wd1 = 62% and Wd2 = 50%. Coefficients that reflect sperm motility reduction, depending on the viability of the donor, r1 = r2 = -0,4. Coefficients that reflect the values of initial spermatozoa motility amounted to S01 = S02 = 80, and maximum motility —
Ks1 = Ks2 = 86.
cells S0(t) R2 = 0,985, and for deconserved cells in experimental group (with crystal growth initiation) — Sd1(t) R2 = 0,987 and in control (without initiation) — Sd2(t) R2 = 0,970, respectively.
Obtained results indicate that seeding has a positive effect on deconserved spermatozoa motility, which increases by 10-15%, and on the effectiveness of cryopreservation technique, which increases by 12% due to reduction in supercooling and, consequently, significantly reduced temperature change during freezing. The difference between control and experiment is reliable at the level of reliability of Р > 0,95.
To evaluate the performance of the proposed model we used a computer experiment that assesses motility of deconserved and fresh spermatozoa obtained from different ejaculates (Fig. 3). Difference in initial sperm motility was 20% and the error of efficiency estimation did not exceed 1%.
When the proposed simulation model was tested using data that was published by us [15], we received results that were highly similar to experimental results that were obtained when the experiment was conducted — R2 = 0,999 for values of fertilizing capacity of deconserved bull sperm (Fig. 4) that was preserved using the proposed technique. This biological quality indicator is a comprehensive characteristic of reproductive cells that were
frozen and thawed given a high preservation safety rate of heifer gametes.
Studies have shown that functions that tie fertilization probability to spermatozoa motility Pd(S) f(7) for fresh and deconserved cells coincide with high reliability probability R2 = 0,99. Effectiveness of cryopreservation amounted to Wd1 = Wd2 = 100%. Consequently, the effect of temperature lowering regimes on the fertilizing capacity of Simmental bull sperm depends solely on sperm motility and fertilization method. When determining the fertilizing capacity, sperm doses that were conserved according to Kharkiv technology [9] were used as a control.
Indicators of sperm motility, survival and absolute survival of deconserved sperm in all experiments conformed to the figures, outlined in GOST 26030-83 "Frozen bovine semen. Technical specifications" for heifer insemination. All frozen semen doses were studied in SEC "Ukraine" of Zmiev district to determine the effectiveness of insemination.
In publications [13, 17, 19, 20] it has been shown that the likelihood of fish embryo development depends on the genotype and physiological characteristics of the female and male that were used to obtain eggs and sperm. The authors found that fertilizing capacity depends not only on the initial state of fresh sperm and eggs, but also on many other factors: cryoprotectant type, freezing-thawing
v = <иХЮк* 10.003 - 0.03 Д*4 ♦ 0.12 S* * • 0,604*J. 2, S JO»« SO .84
&
S S
я О N
О «
Ё U О.
ей
3 4 5 6 Time, hours
10
Fig. 3. Dependence of bull spermatozoa motility on cultivation time at 38 °С from fresh diluted sperm in the first group — S01(t) and in the second — S02(t), deconserved spermatozoa — Sd1(t) и Sd2(t), respectively:
effectiveness of cryopreservation techniques amounted to: Wd1 = 60% and Wd2 = 60%. Coefficients that reflect the decrease in sperm motility, depending on donor viability amounted to r1 = r2 = -0,4. Coefficients that reflect the value of initial sperm motility — S01 = 70 and S02 = 90,
and maximum — Ks1 = 86 and Ks2 = 96
Fig. 4. Dependence of fertilizing capacity of cows on bull sperm motility for fresh — P0(S) and deconserved spermatozoa in the first group — Pd1(S), and in the second — Pd2(S), respectively:
effectiveness of cryopreservation techniques amounted to: W„1 = Wp2 = 95%. Coefficients that reflect the probability of fertilization, depending on the recipient's physiological condition — r1 = r2 = 1,0. Coefficients that reflect the maximum value of animal fertilization Kp1 = Kp2 = 100, and minimum — Pm1 = Pm2 = 3, respectively
mode, method of fertilization and embryo cultivation [13, 17, 19, 20]. A particularly important factor is the individual fertilizing ability of both fresh and deconserved spermatozoa obtained from a particular male [13, 17, 19]. In this regard, a common practice in industry is to fertilize eggs using a mixture of sperm that was obtained from at least three manufacturers [11, 12, 17, 19].
The effect of different cryopreservation techniques on fresh and deconserved carp sperm motility is shown in Fig. 5. The error in determining the effectiveness of cryopreservation technique Wd1 = 65% and Wd2 = 50% amounted to less than 1%.
Fertilizing capacity of spermatozoa depends on the number of viable cells in the sample which, consequently, depends on their concentration and motility [13, 17, 19, 20]. Results of fertility evaluation of ejaculates taken from two males with initial mobility 80,0 ± 3,6% and 70,0 ± 3,5% are shown in Fig. 6.
Fertilizing capacity of carp sperm in the two groups was not significantly different and was 23,3-31,2% for deconserved object (the accuracy of the difference Р > 0,90) and 51-61% for the freshly obtained object. Fertilizing capacity in the first group was higher than in the second, in addition, sperm motility value was higher in the first group for both fresh and deconserved sperm. When
assessing the effectiveness of insemination techniques, the obtained values were: for the first group — 51,0% and for the second — 45,8%. The difference in values for two groups decreased by 1,6-2,0 times due to a decrease in the impact of the initial condition of carp sperm and eggs which depends on their individual fertilizing ability and specific features of insemination and cultivation techniques that are used.
When comparing values for fertilizing capacity of fresh and deconserved carp sperm that were obtained by us with those from the literature [3, 13, 17], it should be noted that the obtained absolute values of the fertilizing capacity of fresh sperm can vary over a wide range from 15 to 80%, and for deconserved spermatozoa — 5-65%. However, values for fertilizing capacity after cryopreservation in publication [13] correspond to those obtained by us, while those in publications [3, 17] are higher, which seems to be associated with the use of mixed sperm from several males during fertilization.
Obtained values of deconserved sperm motility (30-65%) correspond to the data correspond to the data presented in the literature for carp spermatozoa [13], which makes it possible to have fertilizing capacity of deconserved carp spermatozoa ranging from 5 to 65%. It is shown that in the fertilizing
0 20 40 60 80 100 120
Time, sec
Fig. 5. Dependence of carp spermatozoa motility on cultivation time at room temperature for fresh cells — S0(t), deconserved spermatozoa in the first group — Sd1(t) and in the second — Sd2(t), respectively:
effectiveness of cryopreservation techniques amounted to Wd1 = 65% and Wd2 = 50%. Coefficients that reflect the decrease in spermatozoa motility, depending on donor viability — r1 = r2 = -0,04. Coefficients that reflect the minimum value of spermatozoa motility S01 = S02 = 90 and maximum — Ks1 = Ks2 = 96
Fig. 6. Dependence fertilizing ability of carp eggs on spermatozoa motility for fresh cells in the first group — P01(S) and in the second — P02(S), deconserved cells — Pd1(S), rd2(S), respectively:
Cryopreservation effectiveness amounted to Wd1 = 54 and Wd2 = 63%, and insemination effectiveness — Wp1 = Wp2 = 53%. Coefficients that reflect fertilization probability depending on the recipient's physiological state — r1 = r2 = 1,0. Coefficients that reflect maximum values for egg fertilization Кр1 = Кр2 = 100
and minimum — Pm1 = Pm2 = 3, respectively
capacity is significantly influenced by the condition of fresh sperm and eggs, and this in turn depends on the physiological state of the donor, the methods of hormonal treatment and superovulation induction [3, 13, 17].
The use of the proposed model reduces the influence of the condition of fresh carp sperm and eggs, conditions of fertilization and cultivation by a factor of two or more, and allows us to determine the effectiveness of the chosen method of cryopreservation, assessed in terms of fertilizing capacity (54-63%).
Along with the assessment of differences of averages of deconserved carp sperm motility, we also analyzed the variations in the current values of mobility. When transitioning to model values, the variation decreased 2 to 13 times (Fig. 5). This indicates an increase in the accuracy of the proposed method for evaluating cryobiological operations by half or more by reducing the influence of the initial state of carp sperm and related factors that determine the individual characteristics of the selected method of cryopreservation. Such a significant discrepancy between the obtained values can be tied to the effect of "noise" — a difference in the initial condition of sperm and eggs on the one hand, and features of the cryopreservation and fertilization methods — on the other, on the useful "signal" — the condition of the cryopreservation system being examined.
To evaluate the effectiveness of different methods of carp sperm cryopreservation
and their constituent stages, the following computer experiment was performed. Maximum deconserved sperm mobility — 51.7 ± 1.7% (cryopreservation effectiveness 68.9 ± 2.2%) is obtained using a cryopreservation method that involves slow two-step freezing of spermatozoa in a polyethylene envelope (table). As a result of rapid freezing, mobility rates 35.1 ± 5.0% (43.8 ± 1.0%) are close to those rates 36.9 ± 1.2% (50.0 ± 1.3%) that were obtained when spermatozoa were thawed slowly Во = 725 °C/min. The reasons for this — similar values of cryoprotectant effectiveness (92%) and freeze-thaw conditions (60-68%).
The range of motility rates 35-52% obtained by using different methods of cryopreservation of carp sperm (table) gives a similar discrepancy in motility values as when using ejaculates of different quality (Fig. 5) with the error in measurement of less than 1%.
The values of cryopreservation effectiveness 44-67% and the mistake of their measurement also have close values (1%). The decrease in differences of values when transitioning to model values is 1.2-2.7 times. This is indicative of the fact that variation in the initial condition of carp sperm and features of the chosen cryopreservation method have the same effect on deconserved sperm motility. Such a significant discrepancy between the obtained values is due to the influence of different initial condition of sperm samples, on
The effectiveness of different methods of carp sperm cryopreservation Wd and their stages, evaluated in terms of the use of cryoprotectant Wt and freeze-thaw mode Wz
Container Freezing rates В1, В2 (В3 — thawing), °C/min Motility Sd,% Step effectiveness rates,%
cryopreservation, Wd equilibration, We freezing, Wz
Test tube 4, 25 (318) 36.9 ± 1.5a 50.0 ± 1.3a 92.0 + 1.5 68.0 + 1.2a
Straw 4, 25 (725) 44.6 + 1.5b 56.9 + 1.6b 75.2 + 1.5b
Envelope 4, 25 (4780) 51.7 + 1.7с 68.9 + 1.2с 88.6 + 1.6с
Envelope 4500 (12000) 35.1 ± 1.5a 43.8 ± 1.0d 59.8 + 1.6d
Note: B1 and B2 — freezing rate during the first stage from 5 to -20 °C; the second stage from -20 to -60 °C; B3 — thawing rate.
the one hand, and features of cryopreservation methods that were used — on the other hand, on the condition of the deconserved biological object.
It was found that the initial state of the biological object in some cases affects his mobility to a greater extent than the features of the cryopreservation method that was used. In the case when reproducibility of carp sperm motility rates is connected to different semen quality, it reaches 96136%. If the reproducibility depends on the cryopreservation method, it is 94-137%. When using effectiveness rates, it is possible to increase the reproducibility of research results 1,2-4,7 times. To achieve this, it is necessary to take into account the impact of quality of the object on the study and to evaluate the effectiveness of the method cryopreservation method that was used and of its stages. Effectiveness rates of different technological stages that were obtained by us allow to compare results of studies that were conducted using different methods. The most effective method of freezing sperm carp, ensuring sperm safety rate of 45-52%, is based on the use of two-stage freezing. Cryoprotector effectiveness rate is 92%, freeze mode — 75-89%.
Thus, the developed simulation model describes the dependence of deconserved spermatozoa motility on their initial condition and on the effectiveness of cryopreservation technology stages. The use of the model greatly increases the reproducibility of the results of the study and provides conditions for their compatibility when using different technologies and heterogeneous biomaterial.
Comparative analysis of the effectiveness of different stages of bull and carp sperm cryopreservation process indicates that in each case the obtained values depend more on the experimental conditions than on the initial sperm motility. When transitioning to the calculated values of efficiency, it is possible to determine the impact of each studied factor on sperm motility, evaluated at different stages of cryopreservation procedures, while reducing the number of experiments by a factor of three or more.
We believe that the proposed model makes it possible to determine the conditions needed to solve a problem at the stage of exploratory experiments. The proposed method of biotechnological computer experiment can be used in studies that aim to identify hidden patterns of influence that the physiological state of the donor and recipient, different technologies of cryopreservation and artificial insemination have on the fertilizing capacity of deconserved sperm.
The model is a method that ensures the comparability of heterogeneous data obtained empirically under different experimental conditions, which makes it possible to quantitatively explain the strength of the influence that a number of biological and technological factors have on the probability of obtaining offspring. The same model is necessary for multi-factor research aimed at the intensification of existing technologies and the creation of new ones.
Authors would like to express gratitude to the employees at the Institute of Animal Science NAAS and at the Institute of Fisheries NAAS for helping to conduct these studies.
REFERENCES
1. Gorbunov L. Mathematical modeling of the condition of a depreserved biological object. LAPLAMBERT Academic Publishing. Germany. 2014, 116 p. (In Russian).
2. Forero-Gonzalez R. A., Celeghini E. C. C., Raphael C. F., Andrade A. F. C., Bressan F. F., Arruda R. P. Effects of bovine sperm cryopreservation using different freezing techniques and cryoprotective agents on plasma, acrosomal and mitochondrial membranes. Andrologia. 2012, V. 44, P. 154-159. doi: 10.1111/j. 1439-0272.2010.01154.x
3. Irawan H, Vuthiphandchai V., Nimrat S. The Effect of Extenders, Cryoprotectants and Cryopreservation Methods on Common Carp (Cyprinus Carpio) Sperm. Animal Reproduction Science. 2010, 122 (3-4), 236243. doi:10.1016/j.anireprosci.2010.08.017.
4. Celeghini E. C, Arruda R. P., Andrade A. F., Nascimento J., Raphael C. F., Rodrigues P. H. Effects That Bovine Sperm Cryopreservation Using Two Different Extenders Has on Sperm Membranes and Chromatin. Animal Reproduction Scienc. 2008, 104 (2-4), 119131. doi:10.1016/j.anireprosci. 2007.02.001.
5. Critser, J. K, Mobraoten L. E. Cryopreservation of Murine Spermatozoa. ILAR J. 2000, 41(4), 197-206. doi:10.1093/ilar.41.4.197.
6. Pacey, A. A., Tomlinson M. J. The History of Sperm Cryopreservation. Sperm Banking: Theory and Practice. Cambridge, UK: Cambridge UP, 2009.
7. Behma A., Behma L. Energy of spermatozoa, as a criterion in forecasting the fertilization ability of sperm. Biolohiia tvaryn. 2000, 2 (10), 119-124. (In Ukrainian).
8. Behma L., Behma A. Methods of sperm quality evaluation in cryopreservation. Naukovyi zbirnyk. Metodyky naukovykh doslidzhen iz selektsii, henetyky ta biotekhnolohii u tvarynnytstvi. Kyiv: Ahrarna nauka. 2005, P. 161-172. (In Ukrainian).
9. Kharkiv technology of aceptic collection and cryopreservation of sperm from sires: guidelines. Redaktor Ostashko F. Kharkiv. 1990, P. 47. (In Russian).
10. Bugrov A. Development of theory and practice of sire sperm cryopreservation. M. S. thesis in biological sciences. Kharkiv, Ukraina. 1983, P. 19-24. (In Russian).
11. Kopeyka E. Guidelines on low-temperature carp sperm preservation. Moskva. 1986, 9 p. (In Russian).
12. Tsvetkova L. Methodical manual on sperm cryopreservation of carp, salmon and sturgeon fish species. Moskva:VNIIPRKH. 1997, 11 p. (In Russian).
13. Bibenko O. Influence of cryopreservation factors on the genetic apparatus of fish spermatozoa. Ph. D. dissertation in biological sciences. 1993, 16 p. (In Russian).
14. Pavlenko B. Results of studying the cryoconvector method of freezing sperm. Mizhvidomchyi tematychnyi naukovyi zbirnyk. Rozvedennia i henetyka tvaryn. Instytut rozvedennia i henetyky tvaryn. 2002, V. 36, P. 132-133. (In Ukrainian).
15. Mishchenko A., Gorbunov L. Accelerated modes of freezing of biological material. LAPLAMBERTAcademic Publishing. Germany. 2014, 136 p. (In Russian).
16. DSTU 3535-97. Native bull sperm. Technical specifications. (In Ukrainian).
17. Gerasimov Yu. Fundamentals of fish economy: textbook. Samara: Izdatelstvo "Samarskiy universitet". 2003, 108 p. (In Russian).
18. Gorbunov L., Buchackiy L. Cryopreservation of reproductive cells and embryos. Kyiv: Izdatelsko-poligraficheskiy tsentr "Kievskiy universitet". 2005, 325 p. (In Russian).
19. Dzyuba B. Activation, reactivation and fertilization ability of white carp spermatozoa after dilution with cryoprotector environment. Problemy kriobiologii i kriomedetsiny. 1997, N. 3, P. 62-63. (In Russian).
20. Ginzburg A. Detlaf T. Development of sturgeon fish. The maturation of eggs and embryogenesis. Moskva: Nauka. 1969, 134 p. (In Russian).
21. Lakin B. Biometrics. Moskva: Vysshaya shkola. 1990, 254 p. (In Russian).
МОЖЛИВОСТ1 КОМП'ЮТЕРНОГО ЕКСПЕРИМЕНТУ В ДОСЛ1ДЖЕНН1 ГЕТЕРОГЕННОСТ1 СПЕРМАТОЗО1Д1В ТВАРИН
Л. В. Горбунов G. М. Мажарова
Нащональний техшчний унiверситет «Харшвський полггехшчний шститут», Украша
E-mail: 4glv@i.ua
Розроблено iмiтацiйну модель оцiнки рух-ливост й заплщнювально! здатностi сперма-тозоЩв тварин з урахуванням 1х початкового стану та ефективност етапiв кр^консерву-вання. В основу моделi покладено анал^ич-ний вираз, що вщображае основнi причини виживаностi репродуктивних кл^ин в онто-, техно- i ф^огенезь Зниження резистентно-стi сперматозоЗдав залежить в^ низки бш-логiчних (виду тварини, фiзiологiчного стану донора i реципieнта, якостi еякуляту) та технолопчних (ефективностi способiв кр^-консервування i заплiднення яйцеклiтин тварин) чиннишв. Розбiжнiсть мiж значеннями рухливосм клiтин, отриманими розрахунко-вим i експериментальним способами, стано-вила не б^ьше 2% у результатi проведення власних дослав i менше 5% — для даних, наведених у л^ературь Особлив^тю моделi е повна незалежшсть показникiв ефективност дослiджуваних технологш вiд гетерогенност еякуляту тварин.
Проведений комп'ютерний експеримент показав, що розбiжнiсть початково1 рухли-востi й заплщнювально! здатностi спермато-зоЩв залежно вiд варiацil б^лопчних пара-метрiв змiнюеться вiд 50 до 100%, при цьому похибка визначення показника ефективнос-т обрано! технологи становить близько 1%. Порiвняльний аналiз альтернативних тех-нологiй кршконсервування сперматозоЩв показав максимальну ефективнiсть еташв застосування крiопротектора, режиму замо-рожування, виживаност та заплщнювально! здатностi. Застосування моделювання дае змогу багаторазово знизити розкид мiж значеннями рухливост сперматозоЩв, отриманими у рiзних дослiдах, i тим самим б^ьш нiж на порядок скоротити час та грошовi ви-трати, необхщш для одержання достовiрного результату.
Ключовi слова: системна бмлопя, крмкон-сервацiя, гетерогеншсть, сперматозо!ди.
ВОЗМОЖНОСТИ КОМПЬЮТЕРНОГО ЭКСПЕРИМЕНТА В ИССЛЕДОВАНИИ ГЕТЕРОГЕННОСТИ СПЕРМАТОЗОИДОВ ЖИВОТНЫХ
Л. В. Горбунов Е. Н. Мажарова
Национальный технический университет «Харьковский политехнический институт», Украина
E-mail: 4glv@i.ua
Разработана имитационная модель оценки подвижности и оплодотворяющей способности сперматозоидов животных с учетом их начального состояния и эффективности этапов криоконсервирования. В основу модели положено аналитическое выражение, отражающее основные причины выживаемости репродуктивных клеток в онто-, техно- и филогенезе. Снижение резистентности сперматозоидов зависит от ряда биологических (вида животного, физиологического состояния донора и реципиента, качества эякулята) и технологических (эффективности способов криоконсервирова-ния и оплодотворения яйцеклеток животных) факторов. Расхождение между значениями подвижности клеток, полученными расчетным и экспериментальным способами, составило не более 2% в результате проведения собственных опытов и менее 5% — для данных, представленных в литературе. Особенностью модели является полная независимость показателей эффективности исследуемых технологий от гетерогенности эякулята животных.
Проведенный компьютерный эксперимент показал, что расхождение начальной подвижности и оплодотворяющей способности сперматозоидов в зависимости от вариации биологических параметров изменяется от 50 до 100%, при этом погрешность определения показателя эффективности выбранной технологии составляет около 1%. Сравнительный анализ альтернативных технологий криоконсервирования сперматозоидов показал максимальную эффективность этапов применения криопротек-тора, режима замораживания, выживаемости и оплодотворяющей способности. Применение моделирования дает возможность многократно снизить разброс между значениями подвижности сперматозоидов, полученными в разных опытах, и тем самым более чем на порядок сократить время и денежные затраты, необходимые для получения достоверного результата.
Ключевые слова: системная биология, сперматозоиды, гетерогенность, криоконсервиро-вание.
