Research Article
Volume 1 Issue 5  2014
Determination of Cardinal Temperatures of Milk Thistle (Silybum marianum L.) Germination
Zohreh Heidari^{1}, Behnam Kamkar^{2}* and Jafar Masoud Sinaki^{1}
^{1}Agronomy Department, Damghan branch, Islamic Azad University, Damghan, Iran
^{2}Department of Agronomy, Gorgan University of Agricultural Sciences and Natural Resources (GUASNR), Iran
Received: September 16, 2014  Published: November 13, 2014
*Corresponding author: Behnam Kamkar, Department of Agronomy, Gorgan University of Agricultural Sciences and Natural Resources (GUASNR), Faculty of Plant Production, Basij Square, Pardis No. 2, 4918943464 Gorgan, Iran, Tel: +981714427060; Email: @
Citation: Heidari Z, Kamkar B, Sinaki J (2014) Determination of Cardinal Temperatures of Milk Thistle (
Silybum marianum L.) Germination. Adv Plants Agric Res 1(5): 00027. DOI:
10.15406/apar.2014.01.00027
Abstract
Plants response to temperature is an important issue. In order to determine cardinal temperatures and its related germination responses to temperature three varieties seeds of milk thistle (Shomal, Mollasani and Majarestan) were exposed to six constant temperatures (10, 15, 20, 25, 30 and 35 °C) in a four replicated factorial arranged, completely randomized design. Results indicated that the temperature and variety had a significant effect on the maximum germination percentage (MGR), germination uniformity (GU), germination rate (reciprocal time to 50% germination), and time to 5, 10, 50, 90 and 95% germination, but the interacted effect of temperature and variety was not significant. Assessment of three nonlinear regression models, including segmented, dentlike and beta models revealed that the beta model was an appropriate model for all varieties. Base and ceiling temperatures for all three varieties were estimated at near to 0 and 40 °C, respectively. Optimum temperature for Shomal, Mollasani and Majarestan varieties was estimated as (28.79, 29.97 and 29.58 °C), respectively. Biological time for 50% germination for Shomal, Mollasani and Majarestan varieties also was estimated at 43.61, 42.61 and 45.85 hours, correspondingly. It was concluded that these model can be used to quantify response of milk thistle germination to temperature and to obtain cardinal temperature of germination. These parameters are required to predict milk thistle germination and emergence.
Keywords: Beta model; Cardinal temperature; Biological time; Milk thistle; Phenology
Abbreviations
MGR: Maximum Germination Percentage; GU: Germination Uniformity; GR: Germination Rate; RMSE: Root Mean Square of Errors; CV: Coefficient of Variation; R^{2}: Coefficient of Determination; T_{b}: Base Temperature; T_{o}: Optimum Temperature; T_{o1}: Lower Limit of Optimum Temperature; T_{o2}: Upper Limit of Optimum Temperature; T_{c}: Ceiling Temperature; Fo: Biological Time; MGP: Maximum Germination Percentage
Introduction
Milk thistle (Silybum marianum L.) is one of the most important medicinal plants in the pharmaceutical industry worldwide, and is used in the production of flavonoids of the silymarin group (silybin, silidianin and silychristine). The origin of this plant has been reported to be the East Mediterranean region [1]. In this century, extensive research has been conducted on medicinal plants, and medicines that contain natural active ingredients have opened new horizons for the community of physicians and pharmacists engaged in research. Therefore, the pharmaceutical industry and research groups of many countries have turned their attention to the cultivation and production of medicinal plants [2].
Given the importance of these plants in the provision of medicines required by people, the collection, and the completion of, information concerning the ecophysiological behaviors of these plants becomes a necessity. This research was conducted to evaluate nonlinear regression models for describing germination responses of milk thistle to temperature, and to estimate the cardinal temperatures of its germination.
Germination is the first and the most important stage in plant development and is influenced by genetic and environmental factors [3,4]. Germination and establishment of seedlings are among the critical and important stages in the life cycle of plants [5]. Seed germination is the first stage in plant growth and consists of the three steps of water absorption, lag phase and root emergence. Enzymatic activity starts in the first and second steps. During the second step, respiration increases, degradation, and synthesis reactions are started and activated, enzymes break down stored tissues and transfer materials, and finally, in the third step, radicles appear [6].
Temperature and moisture can affect germination percentage and rate either singly or in combination [7]. If sufficient moisture and oxygen are provided, temperature determines the characteristics of mass germination (such as rate, uniformity and percentage of germination) of nondormant seeds [812], and the effects of temperature on plant development are the basis for models that predict germination timing [13]. Temperature is the most important driving force influencing crop development rate [14].
In 1860, it was found that plants have base or minimum, optimum and maximum or ceiling germination temperatures. The minimum or base temperature is the lowest temperature at which germination happens. The optimum temperature is the temperature with the highest germination rate and the maximum or ceiling temperature is the highest temperature at which seeds can germinate [12,15,16]. Cardinal temperatures may significantly differ in various species and varieties [17].
Each plant germinates in a specific temperature range known as cardinal temperatures [18], which includes minimum, optimum and maximum temperatures that have applications in presenting models for predicting germination [12]. Plant growth rate increases with temperature rises from the base to the optimum temperatures and declines from the optimum to the ceiling temperatures [19]. Temperature not only affects germination percentage but also germination rate [20]. Since temperature has considerable effects on germination characteristics, such as the start of percentage and germination percentage and rate, it is, therefore, the most critical factor determining the success or failure of plant establishment [21].
Seeds of different plants respond to temperature differently. Germination rate and percentage reach their maximums at the optimum temperature and decline to zero at the ceiling and the base temperatures [22]. Some studies suggest that, normally, germination rate increases linearly, at least in a suitable temperature range, but falls sharply at temperatures exceeding this range [23]. It is useful to study seed germination and emergence responses to temperature because it allows us to determine the cardinal temperatures for creating models that can predict seed germination and emergence and for selecting the suitable planting time. It also enables us to select species and genotypes that are able to tolerate low or high temperatures and to determine geographical regions where plant species and genotypes can germinate and establish successfully [11,12].
Models play an important role in structuring data and ideas, and can determine the weak areas, and the gaps, in our knowledge. A plant model is a mathematical description of our understanding concerning plant behavior and, since mathematics is used in the creation of the model, it must be completely clear and specified at all stages [24]. Therefore, some researchers have used these models to obtain the cardinal temperatures, that is, the base, the optimum, and the ceiling temperatures [23,25]. There are various functions for describing germination responses to temperature, among which three have been used more often: the segmented, the beta, and the dentlike functions [19,26,27].
Kamkar et al. [28], in their study of poppy by using nonlinear regression models, recommended the segmented model as the superior model for predicting germination rate. Kamkar et al. [14,15] also used segmented and logistic models to determine cardinal temperatures of germination in 3 millet varieties and emergence in wheat cv. “Tajan”. Jame and Cutforth [29] used the beta model in their study on quantification of germination and emergence of spring wheat. Khalili [30] also used nonlinear regression models to quantify germination rate responses of barley to temperature and to water potential, and recommended the beta model as the superior one. GhaderiFar et al. [31], in their evaluation of nonlinear regression models for the three medicinal plants of pumpkin, brago and black cumin, concluded that the dentlike model described germination rate responses to temperature better than other models.
Materials and Methods
This experiment was conducted at the laboratory of Seed Research of Gorgan University of Agriculture and Natural Resources in 2012 using a factorial arrangement of treatments and randomized complete block design in four replications. The experimental treatments included three cultivars of milk thistle (Shomal, Mollasani, and Majarestan) and six levels of temperature (10, 15, 20, 25, 30 and 35 degrees centigrade). Incubators with the accuracy of 0.5 degree centigrade were used for the temperature treatments. Each experimental unit consisted of a petri dish (9cm diameter) in which 50 milk thistle seeds were placed on top of two sheets of filter paper and were covered by a single sheet of filter paper. Subsequently, depending on the temperature used, the numbers of germinated seeds at different time intervals were counted. Obviously, this interval will be shorter for temperatures closer to the optimum one. At each interval, the germinated seeds were removed. The criterion for germination was a radical emergence of two millimeters or more. Counting was stopped when no seeds germinated within 48 hours, or when all the seeds had germinated. Each time the germinated seeds were counted, the exact time of counting was recorded. During the experiment, distilled water was added to the petri dish when required.
The GS2011 software (Kamkar, 2011) was used to extract the deciles of germinations. The software fitted seed germination data from these temperature regimes to the logistic function and performed the necessary calculations through interpolation. Moreover, germination rate (per hour) was determined by using formula 1 [19]:$GR=\frac{1}{{D}_{50}}$
...(1)
In the above formula, GR is the germination rate and D_{50} the time it takes for 50 percent of the seeds to germinate.
In this study, three temperature functions were compared to describe changes in germination rate against temperature. These functions were:
1. The Beta function, shown by the abbreviation B, is as follows [19]:$f(T)=\frac{[(\frac{T{T}_{b}}{{T}_{o}{T}_{b}})\ast (\frac{{T}_{c}T}{{T}_{c}{T}_{b}})(\frac{{T}_{c}{T}_{o}}{{T}_{o}{T}_{b}})]\alpha}{{f}_{o}}$...(2)
2. The Dentlike function, shown by the abbreviation D, is as follows [19]:$f\left(T\right)=\frac{\left(T{T}_{b}\right)}{\left({T}_{o1}{T}_{b}\right)}if{T}_{b}T\le {T}_{o1}$
$f\left(T\right)=\frac{\left(T{T}_{b}\right)}{\left({T}_{o1}{T}_{b}\right)}if{T}_{o2}T\le {T}_{c}$.....(3)
$f\left(T\right)=1if{T}_{o1}T\le {T}_{o2}$
$f\left(T\right)=0ifT\le {T}_{b}or\text{}T\ge {T}_{c}$
3. The Segmented function, shown by the abbreviation S, is as follows [19]:$f\left(T\right)=\frac{\left(T{T}_{b}\right)}{\left({T}_{o}{T}_{b}\right)}if{T}_{b}T{T}_{o}$
$f\left(T\right)=\left[\frac{1\left(T{T}_{o}\right)}{\left({T}_{c}{T}_{o}\right)}\right]if{T}_{o}T{T}_{c}$....(4)
$f\left(T\right)=0if\text{}T\le {T}_{b}or\text{}T\ge {T}_{c}$
In the above function, Tb is the base temperature, To1 the lower limit of optimum temperature, To2 the upper limit of optimum temperature, To the optimum temperature, Tc the ceiling temperature, Fo the biologic time, “a” a constant coefficient, and T the average daily temperature (the temperature at which the experiment was conducted).
Parameters of each model were estimated by a method of iterative optimization with the help of the PROC NLIN procedure of the SAS software [32]. In this method of iterative optimization, every time the initial values of the parameters are entered, the final value is estimated by the leastsquares method. Initial values are changed until the best estimate of the parameters is obtained. The root mean square of errors (RMSE) of germination time, the coefficient of determination (R2), the correlation coefficient (r), the simple linear regression coefficients a and b, and the relationship between the observed and the predicted germination rates were used to select the superior model from among the employed methods.
The higher the coefficient of determination (R2) is, the higher the percentage of changes in germination rate explained by the model. Coefficients “a” and “b” represent the degree of deviation of the regression line from the origin of coordinates and the bias of the regression line from the 1:1 line, respectively. The SAS software determined the numerical values of coefficients “a” and “b.” Significance of coefficient “a” was investigated by SAS and that of coefficient “b” by determining confidence limits. If coefficient “a” is significant, it means that the yintercept of the regression line and the yintercept of the 1:1 line do not coincide. If coefficient “b” is significant, it means that the slopes of the regression line and that of the 1:1 line are different and that the regression line is biased, in relation to the 1:1 line [33]. The less the points are scattered around the 1:1 line, the more efficient the model is. Given these criteria, the superior model was used to determine the cardinal germination temperatures for describing germination rates of various milk thistle varieties.
Results and Discussion
Results of analysis of variance of mean squares of maximum germination percentage (MGR), of germination uniformity (GU), of germination rate (GR), and of the time of reaching 5% (D5), 10% (D10), 50% (D50), 90% (D90) and 95% (D95) germination showed that the effects of temperature on all of the mentioned parameters were significant at the one percent probability level. Varieties did not affect the parameter of reaching 5% germination, but had significant effects on the time of reaching 50% and 90% germination at the one percent level of probability and exhibited significant effects on other parameters at the five percent level of probability. Moreover, the mutual effects of variety and temperature on the parameter of germination uniformity at the one percent level of probability, and on the parameter of reaching 90% germination at the five percent level of probability, were significant (but these mutual effects did not significantly influence the other parameters) (Table 1). Nadjafi et al. [34] studied 11 medicinal plant species and reported that the maximum percentage and rate of milk thistle germination were obtained in the temperature range of 1030 degrees centigrade. The reduction in germination percentage and rate at nonoptimal temperatures can be attributed to the reduction or inhibition of, enzymatic activity and, hence, to a reduction in the required rates of biosynthetic reactions that necessary for germination at these temperatures [35]. It is expected that high temperatures cause seed deterioration in addition to reducing germination rate [25].
Source of variation 
DF 
MGR 
GU 
GR 
D_{5} 
D_{10} 
D_{50} 
D_{90} 
D_{95} 
Variety 
2 
158.16* 
765.47* 
0.000009* 
10.50^{n.s} 
22.58* 
269.38** 
1046.77** 
1206.06* 
Temperature 
5 
4461.96** 
4810.2** 
0.0005** 
25751.97** 
26224.03** 
30934.68** 
36655.20** 
36046.54** 
Variety*Temperature 
10 
46.03^{n.s} 
453.62** 
0.00000^{n.s} 
7.07^{n.s} 
6.81^{n.s} 
16.44^{n.s} 
407.81* 
454.85^{n.s} 
Errors 
54 
44.46 
163.44 
0.000003 
5.54 
7.50 
48.49 
196.05 
330.38 
Table 1: Results of analysis of variance of mean squares of maximum germination percentage (MGR), of germination uniformity (GU), of germination rate (GR), and time to 5% (D5), 10% (D10), 50% (D50), 90% (D90) and 95% (D95) germination.
The symbols ** and * indicate significance at the 1 and 5 percent level of probability, respectively. Values of mean squares and those of F are used for error and for other sources, respectively
Shomal variety 
RMSE 
CV 
R^{2} 
a 
b 
Beta 
0.0006 
3.90 
0.99 
0.0003^{n.s} 
0.9755^{n.s} 
Segmented 
0.0013 
8.79 
0.97 
0.0008* 
0.9511^{n.s} 
Dentlike 
0.0016 
7.71 
0.96 
0.0012^{n.s} 
0.9254^{n.s} 
Mollasani variety 





Beta 
0.0003 
3.72 
0.99 
0.00007^{n.s} 
0.9960^{n.s} 
Segmented 
0.0009 
8.22 
0.98 
0.0004* 
0.9778^{n.s} 
Dentlike 
0.0010 
5.08 
0.98 
0.0005^{n.s} 
0.9727^{n.s} 
Majarestan variety 





Beta 
0.0004 
3.27 
0.99 
0.000006^{n.s} 
0.9994^{n.s} 
Segmented 
0.0006 
7.003 
0.99 
0.0002^{n.s} 
0.9877^{n.s} 
Dentlike 
0.0013 
6.53 
0.97 
0.0009^{n.s} 
0.9440^{n.s} 
Table 2: Root mean square of errors (RMSE), coefficient of variation (CV), coefficient of determination (R2) and coefficient of regression (a and b) for different regression models for Shomal, Mollasani, and Majarestan varieties in the milk thistle.
The symbols * and ** accompanying “a” show significant differences with zero, and in “b” significant differences with one
The beta and the segmented models estimate the base, the optimum, and the ceiling temperatures and the biological time, while the dentlike model estimates the range of optimum temperatures (from the lower limit of optimum temperature to the upper limit of optimum temperature). The criteria used in comparing different models for selecting the superior one to be used in describing germination responses to temperature in various varieties of milk thistle include root mean square of error (RMSE), the coefficient of determination (R2), the regression coefficients “a” and “b”, and the relationship between the observed germination rate and the predicted one (Table 2).
Results obtained from comparison of these criteria are as follow:
In the Shomal variety, the “b” coefficient (the bias of the regression line from the 1:1 line) was not significant in any of the three models, and its value in the beta model (0.97) was closer to one compared to the other two models. Moreover, the “a” coefficient (the degree of deviation of the regression line from the origin of the coordinates) in the beta and dentlike models was not significant, but it was significant in the segmented model at the five percent probability level. The coefficient “a” in the beta model (0.0003) was closer to zero compared to the other two models. The coefficient of determination (R2) in the beta model (99%) was higher compared to the other two models. The root mean square error (RMSE) and the coefficient of variation (CV) in the beta model (0.0006 and 3.9 respectively) were less compared to the other two models (Table 2). The points around the 1:1 line were less scattered in the beta model compared to the other two models. Taking the above points into consideration, and since the “a” coefficient was significant in the segmented method, this model was eliminated in the selection of the superior model for the Shomal variety. The root of mean square error, the coefficient of variation, and the scattering of the points around the 1:1 line were lower in the beta model compared to the dentlike model. Therefore, the beta model was selected as the superior model for the Shomal variety and was used in later stages to determine the cardinal temperatures (Figure 1).
In the Mollasani variety, the “b” coefficient (the bias of the regression line from the 1:1 line) was not significant in any of the three models, and its value in the beta model (0.99) was closer to one compared to the other two models. Moreover, the “a” coefficient (the degree of deviation of the regression line from the origin of the coordinates) in the beta and dentlike models was not significant, but it was significant in the segmented model at the five percent level of probability. The “a” coefficient in the beta model was 0.00007 and closer to zero than in the other two models. The coefficient of determination (R2) in the beta model (99%) was higher compared to the other two models (Table 2). The root mean square error (RMSE) and the coefficient of variation (CV) in the beta model (0.0003 and 3.72, respectively) were less compared to the other two models. The points around the 1:1 line in the beta model was less scattered compared to the other two models. Considering the above points, and because the “a” coefficient was significant in the segmented model, this model was eliminated in the selection of the superior model. The root of mean square error, the coefficient of variation, and the scattering of points around the 1:1 line were less, and the coefficient of determination was higher, in the beta model compared to the dentlike model. Therefore, the beta model was selected as the superior model for the Mollasani variety and was used in later stages to determine the cardinal temperatures (Figure 2).
In the Majarestan variety, the “b” coefficient (the bias of the regression line from the 1:1 line) in none of the three models was significant, and its value in the beta model (0.99) was closer to one compared to the other models. Moreover, the “a” coefficient (the degree of deviation of the regression line from the origin of the coordinates) was not significant in any of three models, and its value in the beta model ( 0.000006) was closer to zero compared to the other two models. The coefficient of determination (R2) in the beta and segmented model, 99%, was higher compared to the dentlike model. The root mean square error (RMSE) and the coefficient of variation (CV) in the beta model (0.0004 and 3.27 respectively) were less compared to the other two models (Table 2). The points around the 1:1 in the beta model were less scattered compared to the other two models. The root mean square error, the coefficient of variation, and the degree of the scattering of the points around the 1:1 lines were less, and the coefficient of determination was higher in the beta model compared to the other two models. Therefore, considering what was earlier said, the beta model was selected as the superior model for the Majarestan variety too, and was used in later stages to determine the cardinal temperatures (Figure 3).
Since the beta model was selected as the superior model for all three studied varieties, its parameters including the base, the optimum, and the ceiling temperatures were used for all three varieties. Based on obtained results:
The estimated base temperatures were 0.00000001, 0.0602, and 0.00000001 degree centigrade for the Shomal, the Mollasani, and the Majarestan varieties, respectively. The estimated optimum temperature for the Shomal variety was 28.79 degrees centigrade, and 29.97 and 29.58 degrees centigrade for the Mollasani and the Majarestan varieties, respectively. The estimated ceiling temperature for all three varieties was 40 degrees centigrade. Moreover, the biological time for 50% germination was equivalent to 43.61, 42.61, and 45.85 hours for the Shomal, Mollasani, and Majarestan varieties, respectively (Table 3).
Researchers have reported the cardinal temperatures of many plants. Nadjafi et al. [34], in a study on cardinal temperatures of germination of 11 medicinal plants, reported that the base, optimum, and ceiling temperatures for milk thistle (Silybum marianum) were 2.7, 18 and 40.3 degrees centigrade, respectively. In cornflower (Centaurea benedicta), they were 0, 13, and 35.5, in hyssop (Hyssopus officinalis) 0, 17.8 and 44.6 in marrubium (Marrubium vulgare) 4, 17.8 and 42.6, and in raceme catnip (Nepeta racemosa), 1.3, 24.1 and 42.6 degrees centigrade. The cardinal temperatures were 3.4, 22.5 and 42.5 degrees centigrade in common evening primrose (Oenothera biennis), 4, 23.2 and 41.3 in basil (Ocimum basilicum), 1, 18.4 and 40.1 in marjoram (Origanum majorana). The cardinal temperatures were 0, 21 and 40.3 degrees centigrade in sage (Salvia sclarea), 0, 17 and 41 in sage (Salvia nemorosa), and 3.1, 22.7 and 42.4 degrees centigrade in thyme (Thymus daenensis subsp. Daenensis). Kamkar et al. [28], in a study they conducted on poppy, found its base, optimum and ceiling temperatures for germination were 3.02, 27.36 and 36.31 degrees centigrade, respectively. Jami AlAhmadi and Kafi [21], in their research on finding the cardinal temperatures of germination in the species kochia (Kochia scoparia), reported that it could germinate in a wide range of temperatures from 3.5 (the base temperature) to 50 (the ceiling temperature), with its optimum germination temperature being 24 degrees centigrade. Zeinali et al. [36] studied germination in 12 wheat varieties and reported their optimum temperatures varied from 31.8 to 32.4, and their ceiling temperatures from 38.1 to 42.1 degrees centigrade. Jame and Cutforth [29] showed in their study that the base, optimum and ceiling temperatures for wheat were 0, 30 and 42 degrees centigrade, respectively. Khalili [30] studied germination in barley and found the base, optimum, and ceiling temperatures were 0.09, 27.5 and 35.7 degrees centigrade, respectively.
Figure 1: The relationship between germination speed and temperature (in degrees centigrade), The predicted values for 50% germination obtained by using the beta (B), the dentlike (D), and the segmented (S) functions are presented against the observed value in the Shomal variety.
Figure 2: The relationship between germination speed and temperature (in degrees centigrade), The predicted values for 50% germination obtained by using the beta (B), the dentlike (D), and the segmented (S) functions are presented against the observed value in the Mollasani variety.
Figure 3: The relationship between germination speed and temperature (in degrees centigrade), The predicted values for 50% germination obtained by using the beta (B), the dentlike (D), and the segmented (S) functions are presented against the observed value in the Majarestan variety.
Beta 
A 
T_{b} 
T_{O} 
T_{C} 
Fo 
Shomal variety 
2.06 
0.00000001 
28.79 
40.00 
43.6165 
Mollasani variety 
1.86 
0.0602 
29.97 
40.00 
42.6132 
Majarestan variety 
1.95 
0.00000001 
29.58 
40.00 
45.8576 
Segmented 





Shomal variety 
 
2.44 
32.44 
40.00 
36.8229 
Mollasani variety 
 
2.48 
33.13 
40.00 
36.0030 
Majarestan variety 
 
2.69 
32.81 
40.00 
38.7295 
Dentlike 
T_{b} 
T_{O1} 
T_{O2} 
T_{C} 
Fo 
Shomal variety 
4.29 
24.40 
35.00 
40.00 
47.1840 
Mollasani variety 
4.34 
24.75 
35.00 
40.00 
45.8527 
Majarestan variety 
3.84 
24.24 
35.00 
40.00 
49.4035 
Table 3: Estimated of constant coefficient (a) (for beta model), base temperature (Tb), optimum temperature (To), lower limit of optimum temperature (To1), upper limit of optimum temperature (To2), ceiling temperature (Tc) and biological time (Fo) by using the nonlinear regression models (beta, segmented and dentlike) for Shomal, Mollasani, and Majarestan varieties in the milk thistle.
Conclusion
Our research showed that of the three common nonlinear regression models for quantifying the effects of temperature on germination rate in the three varieties of Shomal, Mollasani, and Majarestan (the beta, the dentlike, and the segmented models), the beta model fitted the observed data best. Therefore, it was selected as the superior model and was used in the determination the cardinal temperatures of germination (the base, the optimum, and the base temperatures). Consequently, the beta model and the parameters estimated by it may be utilized in creating and evaluating predictive models of germination. These parameters can be saved in data banks of various varieties to be used in determining their cardinal temperatures of germination. Results of our research indicated that, in the studied varieties, there were no great differences between the curves showing responses to temperature, and this makes us more confident in using the obtained values, coefficients and biological time.
References
 Keville K (1991) The Illustrated Herb Encyclopedia; A Complete Culinary, Cosmetic, Medicinal, and Ornamental Guide to Hers. Simon a Schuster Australia. East Roseville, New South Wales.
 Omidbaigi R (2009) Production and Processing of Medicinal Plants. Vol 2, Astan Quds Razavi Publications. p.439 (In Farsi).
 Foley ME, Fennimore SA (1998) Genetic basis for seed dormancy. Seed Sci Res 8(2): 173182.
 Meyer SE, Pendleton RL (2000) Genetic regulation of seed dormancy in purshia tridentata (Rosaceae). Ann Bot 85(4): 521529.
 Windauer L, Altuna A, Arnold RB (2007) Hydrotime analysis of Lesquerella fendleri seed germination responses to priming treatments. Indt Crop Prod 25(1): 7074.
 McDonald MB (2000) Seed priming: 287–525. In: Black M and Bewley JD (Eds.), Seed Technology and Its Biological Basis. Sheffield Academic press, Shettield, Uk, pp.428.
 Bloomberg M, Sedcole JR, Mason EG, Buchan G (2009) Hydrothermal time germination models for radiata pine (Pinus radiata D. Don). Seed Sci. Res. 19(3): 171182.
 Dubetz S, Russell V, Anderson DT (1962) Effect of soil temperature on seedling emergence. Can J Plant Sci 42(3): 481487.
 Heydecker W (1977) Stress and seed germination. An agronomic view In: (Eds.), Khan AA. The physiology and chemistry of seed dormancy and germination. 237276. North Holland. Oxford. Biomedical Press.
 Huidobro JG, Monteith JL, Saquire GR (1982) Time, temperature and germination of pearl millet (Pennisetum typhoides S. and H.) I. Constant temperature. J Exp Bot 33(2): 288296.
 Mwale SS, AzamAli SN, Clark JA, Bradley RG, Chatha MR (1994) Effect of temperature on the germination of sunflower. Seed Sci Technol 22: 565571.
 Ramin AA (1997) The influence of temperature on germination taree Irani. Seed Sci Technol 25(3): 419426.
 Bewley JD, Black M (1994) Seeds: Physiology of Development and Germination, 2^{nd} (edn), Plenum Press, New York, London, USA.
 Kamkar B, Ahmadi M, Soltani A, Zeinali E (2008) Evaluating nonlinear regression models to describe response of wheat emergence rate to temperature. Seed Sci Technol 2: 5357.
 Kamkar B, Koocheki AR, Nassiri Mahallati M, Rezvani Moghaddam P (2006) Cardinal temperature for germination in three millet species (Panicum miliaceum, Pennisetum glaucumand Setaria italica). Asian J Plant Sci 5: 316319.
 Alvarado V, Bradford KJ (2002) A hydrothermal time model explains the cardinal temperatures for seed germination. Plant Cell Environ 25(8): 10611069.
 Phartyal SS, Thapliyal RC, Nayal JS, Rawat MS, Joshi G (2003) The influences of tempeloprasum on seed germination rate in Himalayan elm (Ulmus wallichiana). Seed Sci Technol 31: 8393.
 Bewley JD (1997) Seed germination and dormancy. Plant Cell. 9(7): 10551066.
 Soltani A, Robertson MJ, Torabi B, YousefiDaz M, Sarparast R (2006) Modeling seedling emergence in chickpea as influenced by temperature and sowing depth. Agr Forest Meteorol 138(14): 156167.
 Riemens MM, Scheepens PC, Van der Weide RY (2004) Dormancy, germination and emergence of weed seeds , with emphasis on influence of light. Plant Research International BV 302: 12.
 AlAhmadi JM, Kafi M (2007) Cardinal temperatures for germination of Kochia scoparia(L). Journal of Arid Environments 68(2): 308314.
 Mahmoodi A, Soltani E, Barani H (2008) Germination response to temperature of snail medic (Medicago scutellata L.). Elect J Crop prod 1(1): 5463 (In Farsi).
 Hardegree SP, Winstral AH (2006) Predicting Germination Response to Temperature. II. Threedimensional Regression, Statistical Gridding and Iterativeprobit Optimization Using Measured and Interpolatedsubpopulation Data. Ann Bot 98(2): 403410.
 Banayan M (2002) Development and applications of simulation models in agriculture. Ferdowsi University Press, Iran, pp.208 (in Farsi).
 Hardegree SP (2006) Predicting germination response to temperature. I. Cardinal temperature models and subpopulationspecific regression. Ann. Bot. 97(6): 11151125.
 AkramGhaderi F (2008) The study of seed quality development, germination, longevity and deterioration in some medicinal plants: medicinal pumpkin (Cucurbita pepo L. subsp. Pepo. Convar: Pepo var. styriaca Gerb), cumin blank (Nigella sativa L.) and borago (Borago officinalis L.). Ph.D. Thesis, Gorgan. Univer. Agric Sci Natur Resour p.180.
 Soltani E, Ghaderi FA, Soltani A (2008) Applications of germination modeling on the response to temperature and water potential in seed science research. 1st National Conference of Seed Sci and Technology in Iran. Gorgan, Iran. Plenum Press, p.445.
 Kamkar B, Jami AlAhmadi M, Mahdavi Damghani A, Villalobos FJ (2012) Quantification of the cardinal temperatures and thermal time requirement of opium poppy (Papaver somniferum L.) seeds to germinate using nonlinear regression models. Ind Crop Prod 35(1): 192198.
 Jame YW, Cutforth HW (2004) Simulating the effects of temperature and seeding depth on germination and emergence of spring wheat. Agr Forest Meteorol 124(34): 207218.
 Khalili N (2013) Predicting barley (Hordeum vulgare L.) emergence as affected by temperature, moisture and sowing depth. Gorgan Univer Agric Sci Natur Resour p.60.
 GhaderiFar F, Soltani A, Sadeghipour HR (2009) Evaluation of nonlinear regression models in quantifying germination rate of medicinal pumpkin (Cucurbita pepo L. subsp. Pepo. Convar: Pepo var. styriaca Gerb), borago (Borago officinalis L.) and black cumin (Nigella sativa L.) to temperature. J Agric Sci Nat Resour 16(4): 120.
 Soltani A (2007) Appilication of SAS in statistical analysis. Jihaddaneshgahi; Press, Mashhad, Iran, pp.180 (In Farsi).
 Ahmadi M, Kamkar B, Soltani A, Zeinali E (2010) Evaluation of nonLinear regression models to predict stem elongation rate of wheat (Tajan cultivar) in response to temperature and Photoperiod. Elect J Crop Prod 2(4): 3954.
 Nadjafi F, Tabrizi L, Shabahang J, Damghani AMM (2009) Cardinal Germination Temperatures of Some Medicinal Plant Species. Seed Sci Technol 31(2): 156165.
 Kamaha C, Magure YD (1992) Effect of temperature on germination of six winter wheat cultivars. Seed Sci Technol 20: 181185.
 Zeinali E, Soltani A, Galeshi S, Sadati SJ (2010) Cardinal temperatures, response to temperature and range of thermal tolerance for seed germination in wheat (Triticum aestivum L.) cultivars. Elect J Crop Prod 3(3): 2342.