HOME JOURNALS CONTACT

Pakistan Journal of Biological Sciences

Year: 2003 | Volume: 6 | Issue: 13 | Page No.: 1163-1175
DOI: 10.3923/pjbs.2003.1163.1175
F1 and F2 Backcrosses in the Hybrids Between Two Unnamed Genetically Distinct Species of Tropical Sea Urchins, Echinometra sp. A and Echinometra sp. C
M. Aminur Rahman and Tsuyoshi Uehara

Abstract: Experiments on backcrosses using the gametes of the reciprocal F1 hybids and the conspecific controls of the two sympatric species of tropical sea urchins, Echinometra sp. A (Ea) and Echinometra sp. C (Ec) were conducted in the laboratory. The ova from the female hybrid of Ea x Ec and Ec x Ea yielded a higher percent of fertilization with Ea x Ea sperm than with Ec x Ec sperms. Conversely backcrosses by hybrid males of either crosses yielded a higher percentage of fertilization with Ec x Ec ova than with Ea x Ea ova, indicating that Ea ova appeared to be more descriminating than Ec ova due to the sequence differences in their gamete recognition alleles. On the other hand, the higher fertilization rates between the same types of hybrids versus the different types of hybrids indicated the presence of a complex sorting of gamete compatibility genes. However, the higher fertilization rates as well as higher survival rates of larvae, juveniles and adults of the F2 hybrids eliminates the possibility that hybrid inviability/sterility is a postzygotic mechanism of reproductive isolation. In adults, EaEa x EaEc and EaEa x EcEa juveniles were consistantly larger than those of equivalent ages of the other crosses in terms of relative test dimesions and growth performances; largest to smallest were EaEa x EcEa, EaEa x EaEc, EcEa x EaEa, EaEc x EaEa, EcEa x EcEc, EcEa x EcEc, EaEc x EcEc, EcEc x EcEa and EcEc x EaEc, respectively. These results of growth performances however, indicated that the F2 hybrids were viable in the lab-reared conditions. This experiment represents the first successful production of F2 hybrids (Progeny of backcrosses) in the laboratory. Although, coloration patterns of F2 hybrids tended to be maternal, other characters such as test sizes and spine lengths, growth performances and phenotypic characteristics (such as tubefoot and gonad spicules, pedicellaria valve length and gamete sizes) tended to be intermediate but closer to their maternal F1. Despite these findings, hybrids with these morphological characters were not found in the field, indicating possible prezygotic isolating mechanism(s) that either singlely or incombination with others creating barriers to hybridize in the field. The results of this study further indicate that considerable morphological and genetic integrity is maitained between Echinometra sp. A (Ea) and Echinometra sp. C (Ec) which strongly supports their recognition as distinct evolutionary species despite the fact that they produce F1 and F2 hybrids under laboratory conditions.

Fulltext PDF Fulltext HTML

How to cite this article
M. Aminur Rahman and Tsuyoshi Uehara, 2003. F1 and F2 Backcrosses in the Hybrids Between Two Unnamed Genetically Distinct Species of Tropical Sea Urchins, Echinometra sp. A and Echinometra sp. C. Pakistan Journal of Biological Sciences, 6: 1163-1175.

Keywords: reproductive isolation, Echinometra sp. Gamete compatibility, speciation, sea urchins and F1 and F2 backcrosses

Introduction

A detailed studies on physiological adaptation and reproduction, external features of eggs and sperms, differences in cross-fertilization rates, larval and adult morphologies, karyotypes, spawning seasons and distribution patterns, it has been concluded that four closely related sympatric sea urchins belonging to the genus Echinometra in Okinawa should be recognized as four different species, distinguished as Echinometra spp. A, B, C and D (Ea, Eb, Ec and Ed) (Shingaki and Uehara, 1984; Uehara and Shingaki, 1985; Uehara and Taira, 1987; Arakaki and Uehara, 1991; Uehara, 1990 and Uehara et al., 1986, 1990, 1991). Recent biochemical studies on enzyme electrophoresis (Matsuoka and Hatanaka, 1991); microhabitats (Nishihira et al., 1991); DNA analysis (Metz et al., 1991; Palumbi and Metz, 1991; Palumbi, 1996 and Palumbi et al., 1997) and gamete recognition protein binding (Metz and Plumbi, 1996 and Palumbi, 1998) also suggested these types are separate but closely related species which have recently been speciated. Consequently, mitrochondrial DNA sequence data from the four species of Echinometra show that genetic distance between these species are very small. This result as well as calibration of the rate of mtDNA across sea urchins suggested that the Echinometra in the central and west Pacific diverged over the past 1-3 million years (Palumbi, 1996). Although the four Echinometra are recognized as four distinct species, valid names for these species have been debated (Palumbi, 1996; Palumbi et al., 1997). E. sp. B is now recognized as Echinometra mathaei; sensu stricto (Arakaki et al., 1998); while E. sp. D belongs in the Echinometra oblonga species complex, which may include a cryptic species composed of at least three species (Arakaki and Uehara, 1999). The other two, "sp. A" and "sp. C" have yet to be described and named.

Among broadcast free-spawning marine invertebrates including sea urchins gametic incompatibility is considered as particularly an important mechanism of reproductive isolation, creating barriers to gene flow (Dobzhansky et al., 1977; O'Rand, 1988; Minor et al., 1989; Lessios and Cunningham, 1990; Vacquier et al., 1990, 1995; Palumbi and Metz, 1991; Byrne and Anderson, 1994; Metz et al., 1994 and Palumbi, 1998). Molecules with species-specific composition that attract only conspecific sperm (Ward et al., 1985) and others that do not permit primary binding of heterospecific sperm (Summers and Hylander, 1975, 1976; Bellet et al., 1977; Glabe and Vacquier, 1977, 1978; Glabe and Lennarz, 1979; Vacquier, 1980; Metz et al., 1994; Vacquier et al., 1995 and Metz and Palumbi, 1996) have been found in the egg membrane of echinoids. Though the effectiveness of the barriers in relation to phylogenetic affinity between echinoid species has not yet been determined, Strathmann (1981) has found unidirectional gametic incompatibility between closely related sympatric echinoid species Strongylocentrotus droebachiensis and S. pallidus. The single species spawning have been reported for echinoderms and very recently multiple species spawning of echinoderms has been reported in which participants are sympatric congeners (Pearse et al., 1988 and Byrne and Barker, 1991). In these multiple species events heterospecific gamete fusion may occur, producing viable hybrids (Aslan and Uehara, 1997; Rahman, 1997, 2001 and Rahman et al., 2001) thereby allowing gene flow through hybridization or wasting gametes through formation of inviable hybrids (Pearse et al., 1988).

Hybridization have long been of interest of evolutionary biologists because they offer excellent opportunities to study speciation processes which contribute to reproductive isolation by either pre- or post-zygotic isolating mechanisms (e.g., Strathmann, 1981; Byrne and Anderson, 1994; Palumbi, 1994; Aslan and Uehara, 1997 and Rahman et al., 2001). This interbreeding brings together different genomes, introduces new sources of genetic variation (Willis et al., 1997) and experiences variability in morphological characters (e.g. Wallace and Willis, 1994; Willis et al., 1997; Aslan and Uehara, 1997 and Rahman et al., 2001). Most hybrids are usually morphologically intermediate including in sea urchins in respect of majority of the characters examined (Mayr, 1979; Campton, 1987; 1990; Uehara et al., 1991; Byrne and Anderson, 1994; Wallace and Willis, 1994; Aslan and Uehara, 1997 and Rahman et al., 2001). Laboratory hybridization experiments have been utilized extensively to confirm the probable hybrid nature to certain individuals by demonstrating that two taxa will interbreed when provided with the opportunity to do so, or that gametes from two taxa can be artificially cross-fertilized. Such information provides a valuable contribution to the natural marine hybridization. Laboratory crosses have also demonstrated that there exists considerable potential for hybridization in echinoids which live in sympatry (Uehara et al., 1991; Byrne and Anderson, 1994; Aslan and Uehara, 1997 and Rahman et al., 2001). But none of them had succeeded to find out hybrids form the field collected suspected specimens by using the hybrid characteristics as identifying markers.

Whether all Okinawan Echinometra diversed in allopatry or sympatry is not known, also what mechanisms, pre-or postzygotic maintain their genetic identity today is unsure. Numerous papers have been published on artificial crosses of sea urchin showing that prezygotic isolating mechanisms for instance gamete incompatibility, usually take place (Strathmann, 1981; Uehara et al. 1990 and Metz et al., 1994) revealed by the substantial genetic divergence of sea urchin species. Up to now, only Strathmann (1981) has reported postzygotic isolating mechanisms between hybrid produced experimentally from the two species of sea urchin Strongylocentrotus droebachiensis and S. pallidus because of difficulties of rearing sea urchins (Metz et al., 1994). Conversely, reciprocal adult hybrids produced by using E. sp. A and E. sp. D in the laboratory revealed that the two species are isolated by prezygotic barriers in the field (Aslan and Uehara, 1997).

In our previous experiment we reported that reciprocal F1 hybrids, produced experimentally through cross-fertilization between E. sp. A and E. sp. C were viable, fertile and morphologically intermediate between two of these species (Rahman et al., 2001). To determine the gametic compatibility and developmental compatibility of the F2 hybrids and to assess the degree to which such hybridization may be responsible for genetic introgresion across species borders, a series of cross-fertilization experiments were conducted using the gametes of F1 and F2 hybrids and their conspecific controls.

Materials and Methods

Samples collection and maintenance: Healthy mature adults of F1 conspecifics, Ea x Ea (brownish dark test and white-tipped spines), Ec x Ec (greenish test, spines without white tips, but with white basal ring) and F1 hybrids, Ea x Ec (dark brownish test and spines with a translucent white basal ring), Ec x Ea (uniformly deep greenish test and the spines had barely detectable white tips and basal white rings), produced through cross-fertilization between Echinometra sp. A and Echinometra sp. C (Rahman, 1997) were collected from the culture tanks of Sesoko Marine Research Center and were transported to the laboratory at the Dept. Marine and Environmental Sciences, University of the Ryukyus and maintained in closed aquarium before use for in vitro crosses. The experiment was conducted during two consecutive breeding seasons from May to October, 1998 and from May to September, 1999 in laboratory conditions.

Backcross-fertilization of gametes: Cross-fertilization among the gametes of the F1 conspecifics and the reciprocal F1 hybrids of two Echinometra species, Ea and Ec were conducted using all possible combinations of ova and sperm at room temperature (26-28°C). For simplicity, when referring to the backcrosses, the maternal species is named first. For example, F2 hybrid produced through backcross between the F1 female of Ea x Ea and F1 male of Ea x Ec is denoted as EaEa x EaEc. Fertilization was done by mixing two drops of a diluted sperm into petridishes containing 15 ml egg suspensions. Sperm concentration was generally 10-4 to 10-3 of "dry sperm" (Uehara et al., 1990 and Rahman, 1997). Sperms were allowed to remain with the eggs for ten minutes, excess sperms were then removed with three consecutive washes with FSW (filtered sea water). In each fertilization, a conspecific fertilization by use of ova from the same female was also conducted as a control. The egg were then layered on the bottom of petridishes and incubated at ambient room temperature (26-28°C) for one hour. The first 100 eggs encountered were classifed as "fertilized" if they had reached the 2-4 cell stage. The fertilized eggs were then transferred in glass beakers and incubated in FSW at ambient room temperature until they attained free-swimming blastula stage.

Sperm dilution experiments: The protocol and techniques used in these experiments were different from those employed for production of F2 hybrids. To determine the fertilization rate for both F2 conspecific and backcrosses, a 0.1 ml aliquot of diluted egg suspensions (350-400 eggs) was placed in a small vial with 0.8 ml of FSW. Fresh "dry" sperm was first adjusted to 10°C and then quickly diluted in a series of 7, 10-fold dilutions. A 0.1 ml aliquot from one of these sperm solutions was then placed into the vial, to bring the final volumes to 1 ml. For example, mixing a 0.1 ml aliquot of 10o undiluted sperm with 0.9 ml egg suspensions in a vial was called 10-1 diluted concentration of sperm. This procedure was followed through a series until 10-8diluted concentration of sperm was made. After 5-10 min of gamete mixing, excess sperm were removed by 4-5 consecutive washes with FSW and the eggs were then resuspended in 5 ml of FSW into the vial for incubation. Eighteen replicate crosses were performed with gametes from new individuals each time. Fertilization rate was estimated 1.25-1.5 h after gamete mixing by counting the number of eggs reaching 2-4 cell stages among the first 100 eggs observed.

Larval rearing and culture of juveniles and adults: Only the embryos and larvae of F2 conspecific controls (EaEa x EaEa and EcEc x EcEc) and backcrossed F2 hybrids (EaEc x EaEa, EaEa x EaEc, EcEa x EaEa, EaEa x EcEa, EaEc x EcEc, EcEc x EaEc, EcEa x EcEc and EcEc x EcEa) produced from the above fertilization experiments were continued to rear through metamorphosis, as described by Rahman et al. (2000). After 22-24 days of rearing, the competent larvae were placed in small (25 x 20 x 10 cm) aquaria with aerated filtered sea water and pieces of coralline red algal skeletons were added to induce settlement and metamorphosis (Rahman and Uehara, 2001). Sea water was partially changed once a week with fresh filtered sea water. This was continued for up to three months, by which time the juveniles were 6.0-7.0 mm in test diameter. The juveniles were then transferred to plastic aquaria (36 x 45 x 18 cm) supplied with aerated flow through sea water at Sesoko Marine Research Centre. Coral skeletons covered with encrusting coralline algae, served as food. The cultures were continued for one year by which time the urchins attained sexual maturity. The performances of larvae, juvenile and adults were examined and compared among the F2 hybrid groups and their parental species controls.

Morphological characteristics: Detailed morphological characteristics for describing the differences were recorded or measured from the above one year old F2 conspecifics, EaEa x EaEa, EcEc x EcEc and their F2 hybrid groups were: test size and spine lengths, spicules in the tubefeet and gonads, pedicellaria valve length, color patterns of oral and aboral spines and test and gamete sizes.

The tests of echinoids in the family Echinometridae are oblong so after removing the spines, length, width and height were measured with calipers. For the measurement of spine length, 30 spines were randomly selected from the equator of the test.

Phenotypic coloration patterns of the body and spine of both aboral and oral view of adult urchins were observed by the color book of Kornerup and Wanscher (1978).

The spicules in gonad and tubefoot of the urchins were thoroughly examined. A small piece of the gonadal tissue was clipped off using a forceps and immersed in distilled water. The samples were first treated with 10% KOH and then squashed on a slide glass with a cover slip and finally observed under an objective microscope (10 x 10). Several tubefeet were clipped off with a forceps and the morphology of the spicules was investigated by the same methods as above.

In order to collect the pedicellarias easily, all spines around the body were removed. The pedicellaria were plucked form the test near peristome and ambital area. Soft tissues from the pedicellaria were removed by treating with 10% sodium hypochloride solution (household bleaching). After bleaching, the pedicellaria were rinsed with three changes of distilled water. The valve length (VL) of pedicellaria were then measured under a compound microscope (10 x 10) with presetting micrometer.

Gamete sizes (egg diameter and sperm-head length) were also measured from the sexually matured hybrids and their conspecifics using a differential microscope (eggs at 400x in a slide well, sperm at 1000x on a flat slide), following the methods of Amy (1983) and Rahman et al. (2001, 2002).

F2 backcrosses: After 1 year of culturing, the F2 conspecifics and F2 hybrids reached sexual maturity and contained mature gametes. Gametes were obtained, following the similar method of gamete shedding experiment. Their gametes were then reciprocally backcrossed, following the methods of F1 backcrosses.

Hybrids in nature: To examine the occurrence of natural hybridization between Ea and Ec, field surveys were conducted along Sunabe coast of Okinawa and coast of Sesoko Island, where both of the species occur nearly sympatrically in adjacent microhabitats. On the basis of the phenotypic color patterns of the spines and tests, about 400 suspected hybrids were collected from Sunabe and Sesoko coasts and compared to the laboratory cultured hybrids with respect to coloration patterns of the body and the spines, spicule characteristics and other possible morphological features.

Data analysis: Percentage data where statistical analysis conducted were arcsine transformed. Those replicates in which none or all eggs fertilized, especially in fertilization experiments were given a value of 1/4n and 1-1/4n (n = number of observations) to improve the arcsine transformation (Zar, 1996). This transformation helped to normalized the data and reduce heterogeneity in variances. A "Bartlett's test" was used to analyze the homogeneity of variances (Bartlett, 1937). When variances were not significantly heterogeneous and no major departures from normality, a one way analysis of variance (ANOVA) was done followed by Tukey's multiple comparison test. Data that did not meet the normality assumption of parametric analysis were analyzed using non-parametric statistics. This was done by transforming values to ranks and then applying one way ANOVA followed by Tukey's multiple comparison test. The level for statistical significance was set at 0.05. Untransformed data are presented in tables and figures.

Results

Fertilization rates in F1 backcrosses: Fertilization rates of backcrosses using the gametes of F1 hybrids and their conspecific parental species under various sperm concentrations are depicted in Fig. 1. It is evident that gametes from both the hybrids were completely fertile and reciprocally compatible with the gametes of their F1 conspecifics and vise versa, though the percentages of fertilization showed asymmetrical values in some crosses under lower sperm concentrations (10-8 to 10-3 dilutions) and the identical values for all crosses under higher sperm concentrations (10-2 to 10-1 dilutions) (Fig. 1). Under limited sperm concentration (10-5 dilution), eggs from both Ea x Ec and Ec x Ea yielded higher percentages of fertilization with Ea x Ea sperms (82.27 and 84.11%) than with Ec x Ec sperms (74.56 and 80.39%) (Table 1), similar to the findings with parental crosses where Ea sperm more readily fertilized Ec eggs than Ec sperms fertilized Ea eggs (Rahman et al., 2001). The results obtained from backcrosses among the F1 conspecifics and F1 hybrids however, indicated that gene flow through female hybrids should run predominantly back to E. sp. A rather than to E. sp. C but backcrosses by hybrid males altered the direction of gene flow between these species (i.e., back to E. sp. C rather than to E. sp. A). However, the higher fertilization rates in all backcrosses eliminates the possibility that hybrid infertility/sterility is a postzygotic mechanism of reproductive isolation of these two species. Moreover, backcrosses by sperms from males of Ea x Ec and Ec x Ea yielded higher percentages of fertilization with Ec x Ec ova (86.94 and 93.44%) than with Ea x Ea ova (74 and 74.83%) (Table 1); Ea ova appears to be more discriminating than Ec ova. The fertilization rates between the same types of hybrids versus the different types of hybrids (85.89% versus 76.39% and 83.94% versus 80%) indicate a complex sorting of gamete recognition genes (Table 1). The higher fertilization rates in F1 backcrosses, however, indicated the closer genetic affinity between Ea and Ec.

Table 1:
Percentage of eggs fertilized in backcrosses among laboratory-reared F1 generation of conspecifics and hybrids of Echinometra sp. A (Ea) and Echinometra sp. C (Ec) at limited sperm concentration (10-5 dilution). Counts made 1.5 hours after gamete mixing of number of eggs out of 100 that had reached 2-4 cell stage (fertilized). Each value represents 18 replicate crosses with gametes from new individuals in each replicate; mean±SD, ranges in parentheses

Table 2:
Comparison of larval, Juvenile and adult performances of F2 hybrids, produced experimentally through backcrossing among F1 conspecifics and the reciprocal hybrids of Echinometra sp. A (Ea) and Echinometra sp. C (Ec). Five replicate experiments were conducted in each cross for each type of performance. All values represent mean±SD with ranges in parentheses. Mean values within each experiment with common superscripts in the same column are not significantly different ((Tukey's test, P > 0.05)
*Matured larvae that were deemed competent for metamorphosis after a 22-24 days culture period in laboratory condition.
**Three months old juvenile urchins that were transferred to flow-through sea water system for advanced culture.

Table 3:
Comparison of test sizes and spine lengths of F2 hybrids, produced experimentally through backcrossing among F1 conspecifics and the reciprocal F1 hybrids of Echinometra sp. A (Ea) and Echinometra sp. C (Ec) one year after metamorphosis. A total of 20 adult individuals were measured for each cross. All values represent (mean±SD) in mm with ranges in parentheses
*Mean values within each experiment with common superscripts in the same column are not significantly different (Tukey's test, P > 0.05)

Table 4:
Comparison of types and percentages of tubefoot spicules of F2 hybrids, produced experimentally through backcrossing among F1 conspecifics and the reciprocal F1 hybrids of Echinometra sp. A (Ea) and Echinometra sp. C (Ec) one year after metamorphosis. A total of 20 adult individuals were measured for each treatment. All values represent (mean±SD) in percentage with ranges in parentheses
*Mean values within each experiment with common superscripts in the same column are not significantly different (Tukey's test, P > 0.05)

Table 5:
Comparison of types and percentages of gonad spicules of F2 hybrids, produced experimentally through backcrossing among F1 conspecifics and the reciprocal F1 hybrids of Echinometra sp. A (Ea) and Echinometra sp. C (Ec) one year after metamorphosis. A total of 20 adult individuals were measured for each treatment. All values represent (mean±SD) in percentage with ranges in parentheses
*Mean values within each experiment with common superscripts in the same column are not significantly different ((Tukey's test, P > 0.05).

Larval, juvenile and adult performances: Despite, there were slight differences in the survival of matured larvae, the values did not differed significantly (Tukey's test, P > 0.05) among the F2 conspecifcs and F2 hybrids. Metamorphosis, recovery and adult survival were also followed the same trends as larval survival (Table 2). The mean live weight attained by the one year old adult EaEa x EcEa, EcEa x EaEa, EaEa x EaEc and EaEc x EaEa was non-significantly very closer to their F2 parental EaEa x EaEa, but all of them were differed significantly from slow-growing EcEc x EcEc. Whereas, the same for EcEc x EcEa, EcEa x EcEc, EcEc x EaEc and EaEc x EcEc was closer to their F2 conspecific EcEc x EcEc, though the values of these crosses were differed significantly (Tukey's test, P < 0.05) (Table 2). Similar trends were observed in test sizes (length, width and height) and spine length of adult F2 hybrids and their conspecific controls (Table 3). Therefore, the F2 hybrids in repect of all growth performances were tended to be intermediate but very closer to the maternal character of their F2 conspecifics. Similar results were also obtained in the F1 hybrids between Echinometra sp. A and Echinometra sp. C (Rahman et al., 2000, 2001).

Table 6:
Valve length (VL) of four types of pedicellariae in F2 hybrids, produced experimentally through backcrossing among F1 conspecifics and the reciprocal F1 hybrids of Echinometra sp. A (Ea) and Echinometra sp. C (Ec) one year after metamorphosis. A total of twenty adult individuals were examined from each cross with 10 pedicellariae of each type from each individual. All values represent mean±SD in Fm followed by the ranges in parentheses
*Mean values within each experiment with common superscripts in the same column are not significantly different ((Tukey's test, P > 0.05).

Table 7:
Gamate sizes of F2 hybrids, produced experimentally through backcrossing among F1 conspecifics and the reciprocal F1 hybrids of Echinometra sp. A (Ea) and Echinometra sp. C (Ec) one year after metamorphosis. Twenty adult individuals were examined from each cross with 25 eggs and 25 sperm from each individual. All values represent mean±SD in Fm followed by the ranges in parentheses
*Mean values within each experiment with common superscripts in the same column are not significantly different ((Tukey's test, P > 0.05)

Table 8:
Percentage of eggs fertilized in backcrosses among laboratory-reared F2 generation of conspecifics and hybrids of Echinometra sp. A (Ea) and Echinometra sp. C (Ec) at limited sperm concentration (10-5 dilution). Counts made 1.5 hours after gamete mixing of number of eggs out of 100 that had reached 2-4 cell stage (fertilized). Each value represents 9 replicate crosses with gametes from new individuals in each replicate. All values represent mean in percentage with ranges in parentheses

Fig. 1:
Mean percentages of fertilization in F1 backcrosses using the gametes from laboratory cultured Ea x Ea, Ec x Ec and their reciprocal F1 hybrids under a series of sperm concentrations ("dry" sperm dilution). A total of 18 replicate crosses were done for each combination using the gametes from new individuals in each time. Standard deviations were not included, as these values were usually smaller than the symbols

Comparison of phenotypic characteristics: Major phenotypic color patterns of the F2 hybrids and their conspecific parents were examined at the end of the experiment. From aboral view, the spine and test of EaEa x EaEa were brownish dark and clear white ring at the base of each white tipped spine, whereas EcEc x EcEc were uniformly green without white tip but each spine had fadded white basal ring and the test color was dominated by deep green. EaEa x EcEa, EaEa x EaEc, EaEc x EaEa and EaEc x EcEc hybrids were more similar to EaEa x EaEa in having brownish dark test and white tipped spines with a clear white basal ring. On the other hand, EcEc x EcEa, EcEa x EcEc, EcEa x EaEa and EcEc x EaEc hybrids were more closer to EcEc x EcEc, having uniformly deep green test and spines and each spine had a barely detectable white tips and basal white rings. At oral coloration, EaEa x EaEa had white tipped spines around the mouth and brownish dark test, whereas EcEc x EcEc had yellowish green spines around the mouth and a greenish test color. EaEa x EcEa, EaEa x EaEc, EaEc x EaEa and EaEc x EcEc hybrids were more similar to EaEa x EaEa, whereas EcEc x EcEa, EcEa x EcEc, EcEa x EaEa and EcEc x EaEc hybrids were more similar to EcEc x EcEc. Therefore, majority of the F2 hybrid urchins showed the coloration patterns more closer to their maternal species. This may be due to the fact that maternal genomes are prefentially expressed in their F2 progenies. A few of them also experienced intermediate colors between their parental progenies, may be because of the combine effects of their parental genomes. Similar color patterns were also observed in the F1 hybrids between Echinometra sp. A and Echinometra sp. C (Rahman et al., 2001) and between Echinometra sp. A and Echinometra oblonga (Aslan, 1995).

Tubefoot spicules in EaEa x EaEa were always bihamate (100%), whereas EcEc x EcEc were triradiate (100%) (Table 4). Although tubefoot spicules of EaEa x EcEa, EcEa x EaEa, EaEa x EaEc and EaEc x EaEa hybrids were bihamate (65.73-82.62%), bihamate-like (6.25-29.08%), triradiate-bihamate (7.42-13.69%) and triradiate (1.01-3.12%) but they were dominated by bihamate type (Table 4), whereas those of EcEc x EcEa, EcEa x EcEc, EcEc x EaEc and EaEc x EcEc were bihamate (13.09-25.44%), bihamate-like (15.82-27.57%), triradiate-bihamate (8.07-19.59%) and triradiate (38.92-51.50%) but dominated by triradiate type (Table 4). The F2 hybrids and their conspecific F2 parents differed significantly (Tukey's test, P < 0.05) in respect of corresponding similar types of spicules and no intermediate types were found in conspecifics (Table 4). Therefore, the tubefoot spicule morphologies of hybrids were intermediate and tended towards maternal affinities.

The spicules of gonad in EaEa x EaEa were almost all spindle-shaped (98.4%); other spicules seen: bihamete (0.6%), Irregular (1.0%), no triradiates were found (Table 5). In contrast, those in Ec x Ec gonads were nearly all triradiate (94.0%); other spicules seen: spindle (4.9%) and bihamete (1.1%) (Table 5). Gonads in EaEa x EcEa, EcEa x EaEa, EaEa x EaEc and EaEc x EaEa hybrids had high proportions of spindle-shaped spicules (58.4-82.1%) with spindle-like (6.7-27.8%), spindle-triradiate (6.3-10.1%), triradiate (1.9-3.9%) and bihamate (0.8-1.0%) and irregular (0.8-1.1%) in smaller proportions whereas, those in EcEc x EcEa, EcEa x EcEc, EcEc x EaEc and EaEc x EcEc hybrids had high proportion of triradiate-shaped spicules (38.7-49.7%) followed by spindle-like (17.5-29.9%), spindle (12.6-24.9%) and spindle-triradiate (5.0-18.5%) with a very small proportion of bihamate (0.8-1.1%) and irregular (0.7-1.0%) types. The very fewer proportions of bihamate and irregular types were not differed significantly (Tukey's test, P > 0.05) among the conspecific and hybrid groups (Table 5). Although, both the hybrids showed various proportions of major types of spicules, the values of each types were significantly different (Tukey's test, P > 0.05). Though, some significant differences were recognized among the spicule types and proportions, hybrids of both directions seemed to be more closer to those of their maternal types (Table 5).

The pedicellaria found in both F2 conspecifics and their hybrids were tridentate, globiferous, ophiocephalous and triphyllus. Only the valve length (VL) of all four types were measured and compared among conspecifics and their hybrid groups (Table 6). As shown in Table 6, all four types of pedicellaria VL of EaEa x EaEa were significantly (Tukey's test, P < 0.05) larger than those of their corresponding types of EcEc x EcEc. All the hybrid groups experienced intermediate sizes, but non-significantly closer to their maternal controls (Table 6).

Egg diameters of EcEc x EcEc control were significantly (Tukey's test, P < 0.05) larger than the EaEa x EaEa eggs. Hybrids contained intermediate sized eggs which were non-significantly closer to their maternal conspecifics (Table 7). The size of sperm heads were also smallest in EaEa x EaEa and were significantly different (Tukey's test, P < 0.05) from all the crosses. Though, the values differed among the hybrid groups, all were intermediate between values of Ea x Ea (smallest) and Ec x Ec (largest) (Table 7).

Existence of natural hybrids: Five hundred individuals with more or less similar coloration to the laboratory reared F2 hybrids were collected from intertidal zone of the sea where their parental Ea and Ec inhabit at closer proximity. However, detailed comparisons of the above morphological characters revealed that none of these individuals actually had common character combinations to the experimentally obtained F2 hybrids. Similarly, Aslan and Uehara (1997) and Rahman et al. (2001) did not find any natural hybrids between Ea and Ed and between Ea and Ec in the field even though they were able to produce sexually matured hybrids of these species in the laboratory.

Fertilization rates in F2 backcrosses: Fertilization rates in F2 backcrosses using the gametes of F2 hybrids and their conspecific F2 parental species were also conducted under various sperm concentrations. Only the fertilization rates at limited sperm concentration (10-5 dilution of "dry" sperm) are shown in Table 8. At this concentration of "dry" sperm, ova from all F2 hybrid groups (aa x ca, ca x aa, aa x ac, ac x aa, cc x ca, ca x cc, ac x cc and cc x ac) exhibited similar high rates of fertilization (86.9-88.8%) with the sperm from F2 conspecific aa x aa, whereas ova of aa x ca, ca x aa, aa x ac and ac x aa hybrids produced very low percentage of fertilization (43.1-46.2%) than the ova of cc x ca, ca x cc, ac x cc and cc x ac (98.3-100%) with the same sperm from F2 conspecific cc x cc (Table 8). Consequently, ova from F2 conspecific cc x cc produced higher percentages of fertilization (83.8-100%) with the sperm from each of all hybrid groups than the percent fertilization exhibited by the ova from conspecific aa x aa (21.9-89.1%) with the sperm from each groups of F2 hybrids (Table 8). However, the sperm of F2 hybrids appeared to be more discriminating with ova of F2 conspecific aa x aa than with F2 conspecific cc x cc. Thus, the differences in fertilization rates of the gametes among the F2 conspecifics and their F2 hybrids clearly indicates the sequence differences in their gamete recognition alleles. Similarly, the discrimination of fertilization between and among the hybrid groups also indicates the presence of polymorphic gamete recognition genes in their binding alleles. However, the F2 hybrids were completely fertile, indicating that Ea and Ec are genetically very close.

Discussion

Hybridization was asymmetrical however and fertilization rate was generally lower in some backcrosses. This reduction, at least among some heterogametic crosses, indicates the presence of a gamete recognition protein binding system, as reported by Metz et al. (1994), Metz and Palumbi (1996) and Palumbi (1998) which might eventually lead to gametic incompatibility and reproductive isolation. However, the high fertilization rates in many of the heterogametic crosses suggested that gametic incompatibility among these species and their hybrids is unlikely to provide a mechanism for maintaining species integrity in these species (Metz and Palumbi, 1996; Aslan and Uehara, 1997; Rahman, 1997; Rahman; 2000 and Rahman et al., 2001). The compatibility of the gametes of the F1 and F2 hybrids of E. sp. A and E. sp. C demonstrates that if gamete recognition molecules are involved in fertilization in these species, they are not strongly species specific (Byrne and Anderson, 1994). Moreover, the similarly higher larval survival, metamorphosis and recovery rates of F2 hybrids in laboratory conditions eliminate the possibilities that hybrid inviability, sterility or breakdown is the postygotic mechanism of reproductive isolation.

The growth patterns of all the ten combinations showed similar trends in triplicate aquaria. It was observed that one year old adult of EaEa x EcEa, EcEa x EaEa, EaEa x EaEc and EaEc x EaEa while being very closer to their F2 parental EaEa x EaEa, registered similar faster growth than the other groups. Whereas, the same for EcEc x EcEa, EcEa x EcEc, EcEc x EaEc and EaEc x EcEc was slower but close to their F2 conspecific EcEc x EcEc, the slowest. The growth of F2 hybrids, produced through backcrossing among the ova of F1 hybrids and sperm of their F1 conspecifics were almost the similar during the culture period of one year. Subsequently, the growth of EaEa x EcEa, EcEa x EaEa, EaEa x EaEc and EaEc x EaEa was surpassed by EcEc x EcEa, EcEa x EcEc, EcEc x EaEc and EaEc x EcEc, mostly appears to be intermediate between that of their parental species. The first-growing F1 parents partially transmitting this trait to their F2 hybrids. The growth and survival data however, confirmed that hybridization among F1 generation of E. sp. A and E. sp. C has been successful in laboratory-rearing conditions and the resulting F2 hybrids are as viable as conspecifics and showed evidence of parental heterosis. Rahman et al. (2000, 2001) also observed the similar phenomena in the F1 hybrids between Echinometra sp. A and Echinometra sp. C.

Although the expression of an intermediate phenotypes by the laboratory-reared F1 (previous Study) and F2 hybrids (present study) has assisted to find out hybrid genotyes in the field (Aslan and Uehara, 1997; Rahman, 1997 and Rahman et al., 2001). The coloration patterns of both the F1 and F2 hybrids tended to be maternal. Other remarkable characters such as test sizes, growth performances, spine lengths, gonad and tubefoot spicules, gamete sizes tended to be intermediate. In our previous experiment (Rahman, 1997), we had done an intensive investigation to discover the F1 hybrids in the field by using these morphological characters, but we were failed to find any evidence of natural hybridization. This may be due to the fact that there might be some sorts of isolating mechanisms (such as gametic incompatibility, habitat segregassion and gamete competition) prevent these two species from hybridizing in the field and so that no introgression takes place despite their sympatric existence (Rahman, 2000; Rahamn et al., 2001). Though we were not able to find any hybrids in the field, the higher fertilization and survival rates of the F1 and F2 hybrids could open the door to do further investigations by using molecular and genetic markers. If genetic analyses, using allozymes and DNA markers also fails to find evidence of hybrids (F1 or F2) in the field, there must be some effective isolating mechanisms that separates these two congeners into distinct species.

A good result regarding the fertilization rates through backcrossing between the hybrid and parental progeny suggests that these two species may share the same gene pool with their hybrid progeny. In other words, these species are genetically very close to each other and perhaps no effective introgression takes place. The recent hybridization experiment conducted by Lessios and Pearse (1996) also revealed that gene introgression among the three tropical Indo-Pacific species of the genus Diadema (D. paucispinum, D. savignyi and D. setosum) in Okinawa is limited.

Whether the Echinometra sp. A and Echinometra sp. C are distinct species or not, the results from our hybridization experiments (F1 and F2) indicate that the two species are very close in respect to their genetic similarity and phylogenetic relationships and it seems almost certain that the two species are the recent derivatives from one ancestral species and also sufficient for their speciation as suggested by Matsuoka and Hatanaka (1991) and Palumbi and Metz (1991). Whether incipient gametic incompatibility plays an important role in speciation of these species as suggested by Metz et al. (1994), Metz and Palumbi (1996) and Aslan and Uehara, (1997) or if some other mechanism (s) is involved deserve further study. We believe that the genetic integrity maintained between Echinometra sp. A and Echinometra sp. C warrants their recognition as distinct evolutionary species despite the viable F1 and F2 hybrids were produced experimentally through cross-fertilization in the laboratory.

Acknowledgments

We are much indebted to the Director and staff of the Sesoko Station of the Tropical Biosphere Research Center, University of the Ryukyus, Okinawa, who provided space and physical facilities for culturing sea urchins. The first author also wish to extend his appreciation and grateful thanks to Rotary Yoneyama Memorial Foundation, Naha Club, Okinawa for scholarship grants during the period of this experiment.

REFERENCES

  • Amy, R.L., 1983. Gamete sizes and developmental time tables of five tropical sea urchins. Bull. Mar. Sci., 33: 173-176.


  • Arakaki, Y. and T. Uehara, 1999. Morphological comparison of black Echinometra individuals among those in the Indo-west pacific. Zool. Sci., 16: 551-558.


  • Arakaki, Y., T. Uehara and I. Fagoone, 1998. Comparative studies of the genus Echinometra from okinawa and mauritius. Zool. Sci., 15: 159-168.
    CrossRef    


  • Aslan, L. M. and T. Uehara, 1997. Hybridization and F1 backcrosses between two closely related tropical species of sea urchins (genus Echinometra) in Okinawa. Invert. Reprod. Dev., 31: 319-324.
    CrossRef    


  • Bartlett, M.S., 1937. Some examples of statistical methods of research in agriculture and applied biology. Suppl. J. R. Stat. Soc., 4: 137-183.
    CrossRef    Direct Link    


  • Bellet, N.F., J.P. Vacquier and V.D. Vacquier, 1977. Characterization and comparison of bindin isolated from sperm of two species of sea urchins. Biotechnol. Bioph. Res. Commun., 79: 159-165.


  • Byrne, M. and M.F. Barker, 1991. Embryogenesis and larval development of the asteroid Patiriella regularis viewed by light and scanning electron microscopy. Biol. Bull., 180: 332-345.


  • Byrne, M. and M.J. Anderson, 1994. Hybridiation of sympatric Patiriela species (Echinodermata: Asteroidea) in New South Wales. Evolution, 48: 564-576.


  • Cameron, R.A. and S. Schroeter, 1980. Sea urchin recruitment: Effect of substrate selection on juvenile distribution. Mar. Ecol. Prog. Ser., 2: 243-247.


  • Glabe, C.G. and V.D. Vacquier, 1977. Species-specific aglutination of eggs by bindin isolated from sea urchin sperm. Nature, 267: 836-838.


  • Glabe, C.G. and V.D. Vacquier, 1978. Egg surface glycoprotein receptor for sea urchin sperm bindin. Proc. Nat. Acad. Sci. USA., 75: 881-885.


  • Glabe, C.G. and W.J. Lennarz, 1979. Species-specific sperm adhesion in sea urchins: A quantitaive investigation of bindin-mediated egg agglutination. J. Cell Biol., 83: 595-604.


  • Kornerup, A. and J.H. Wanscher, 1978. Methuen Handbook of Colour. 3rd Edn., Eyre Methuen Ltd., London, pp: 259


  • Lessios, H.A. and C.W. Cunningham, 1990. Gametic incompatibility between species of the sea urchin genus, Echinometra on the two sides of the Isthmus of Panama. Evolution, 44: 933-941.


  • Lessios, H.A. and J.S. Pearse, 1996. Hybridization and introgression between Indo-Pacific species of Diadema. Mar. Biol., 126: 715-723.


  • Matsuoka, N. and T. Hatanaka, 1991. Molecular evidence for the existence of four sibling species within the sea urchin, Echinometra mathaei in Japanese waters and their evulutionary relationships. Zool. Sci., 8: 121-133.
    Direct Link    


  • Mayr, E., 1979. Animal species and Evolution. 6th Edn., Belknap Press, Cambridge, Massachusetts, USA


  • Metz, E.C. and S.R. Palumbi, 1996. Positive selection and sequence rearrangements generate extensive polymorphism in the gamete recognition protein bindin. Mol. Biol. Evol., 13: 397-406.


  • Metz, E.C., R.E. Kane, H. Yanagimachi and S.R. Palumbi, 1994. Fertilization between closely related sea urchins is blocked by incompatibiliies during sperm-egg attachment and early stags of fusion. Biol. Bull., 187: 23-24.
    Direct Link    


  • O'Rand, M.G., 1988. Sperm-egg recognition and barriers to interspecies fertilization. Gamete Res., 19: 315-328.


  • Palumbi, S.R., 1994. Genetic divergence, reproductive isolation and marine speciation. Annu. Rev. Ecol. Syst., 25: 547-572.
    Direct Link    


  • Palumbi, S.R., 1996. What can molecular genetics contribute to marine biogeography an urchini's tale. J. Exp. Mar. Biol. Ecol., 203: 75-92.


  • Palumbi, S.R. and E.C. Metz, 1991. Strong reproductive isolation between closely related tropical sea urchins (genus Echinometra). Mol. Biol. Evol., 8: 227-239.


  • Palumbi, S.R., G. Grabowsky, T. Duda, L. Geyer and N. Tachino, 1997. Speciation and population genetic structure in tropical Pacific sea urchins. Evolution, 5: 1506-1517.


  • Pearse, J.S., D.J. McClay, M.A. Sewell, W.C. Austin, A. Prez-Ruzafa and M. Byrne, 1988. Simultaneous spawning of six species of echinoderms in Barkley Sound, British Columbia. Invert. Reprod., 14: 279-288.


  • Rahman, M.A., 1997. Hybridization between two closely related species of sea urchins (genus Echinometra) with special reference to their phenotypes and F1 backcrosses. M.Sc. Thesis, Department of Biology, University of the Ryukyus, Japan, pp: 101.


  • Rahman, M.A., 2000. Comparative Studies on two closely related species of sea urchins (genus Echinometra) in Okinawa: Fertilization, hybridization and larval and juvenile growth with respect to species integrity and aquaculture potential. M.Sc. Thesis, University of the Ryukyus, Okinawa, Japan, pp: 151.


  • Rahman, M.A., T. Uehara and L.M. Aslan, 2000. Comparative viability and growth of hybrids between two sympatric species of sea urchins (Genus Echinometra) in Okinawa. Aquaculture, 183: 45-56.
    CrossRef    


  • Rahman, M.A. and T. Uehara, 2001. Induction of metamorphosis and substratum preference in four sympatric and closely related species of sea urchins (genus Echinometra) in Okinawa. Zool. Stud., 40: 29-43.


  • Rahman, M.A., T. Uehara and J.S. Pearse, 2001. Hybrids of two closely related tropical sea urchins (genus Echinometra): Evidence against postzygotic isolating mechanisms. Biol. Bull., 200: 97-106.
    PubMed    


  • Rahman, M.A., T. Uehara and S.M. Rahman, 2002. Efeects of egg size on fertilization, fecundity and offspring performance: A comparative study between two sibling species of tropical sea urchins. Pak. J. Biol. Sci., 5: 114-121.


  • Shingaki, M. and T. Uehara, 1984. Chromosome studies in the several species of sea urchin from Okinawa. Zool. Sci., 1: 1008-1008.


  • Strathmann, R.R., 1981. On barriers to hybridization between Strongylocentrotus droebachiensis (O.F. Muller) and S. pallidus (G.O. Sars). J. Exp. Mar. Biol. Ecol., 55: 39-97.


  • Summers, R.G. and B.L. Hylander, 1975. Species-specifity of acrosome reaction and primary gamete binding in echinoids. Exp. Cell Res., 96: 63-68.


  • Summers, R.G. and B.L. Hylander, 1976. Primary gamete binding: Quantitative determination of its specificity in echinoid fertilization. Exp. Cell Res., 100: 190-194.


  • Uehara, T., 1990. Speciation in Echinometra mathaei. Iden, 44: 47-53.


  • Uehara, T. and M. Shingaki, 1985. Taxonomic studies in the sea urchin, Echinometra mathaei from Okinawa, Japan. Zool. Sci., 3: 1114-1114.


  • Uehara, T., M. Shingaki and K. Taira, 1986. Taxonomic studies in the sea urchin, genus Echinometra from Okinawa and Hawaii. Zool. Sci., 3: 1114-1114.


  • Uehara, T. and K. Taira, 1987. Heteromorphyic chromosomes in sea urchin, Type B of Echinometra mathaei, from Okinawa. Zool. Sci., 4: 1001-1001.


  • Vacquier, V.D., K.R. Carner and C.D. Stout, 1990. Species-specific sequences of abalone lysin, the sperm protein that creates a hole in the egg envelope. Proc. Nat. Acad. Sci., 87: 5792-5796.


  • Vacquier, V.D., W.J. Swanson and M.E. Hellberg, 1997. What have we learned about sea urchin sperm bindin?. Dev. Growth Diff., 37: 1-10.


  • Wallace, C.C. and B.L. Willis, 1994. Systematics of the coral genus Acropora: Implications of new biological findings for species concepts. Annu. Rev. Ecol. Syst., 25: 237-262.


  • Ward, G.E., C.J. Brokaw, D.L. Garbers and V.D. Vacquier, 1985. Chemotaxis of Arbachia punctulata spermatozoa to resact, a peptide from the egg jelly layer. J. Cell Biol., 101: 2324-2329.


  • Willis, B.L., R.C. Babcock, P.L. Harrison and C.C. Wallace, 1997. Experimental hybridization and breeding incompatibilities within the mating systems of mass spawning reef corals. Coral Reefs, 16: 553-565.


  • Zar, J.H., 1996. Biostatistical Analysis. 3rd Edn., Prentice Hall International, Inc., Upper Saddle River, New Jersy


  • Aslan, L.M., 1995. Hybridization between two closely related tropical species of sea urchins (genus Echinometra) in Okinawa. M.Sc. Thesis, Department of Biology, University of the Ryukyus, Japan, pp: 59.


  • Dobzhansky, T., F.J. Ayala, G.L. Stebbins and J.W. Valentine, 1977. Evolution. W.H. Freeman, San Francisco


  • Arakaki, Y. and T. Uehara, 1991. Physiological Adaptation and Reproduction of the Four Types of Echinometra mathaei (Blainville). In: Biology of Echinodermata, Yanagisawa, T., I. Yasumasu, C. Oguro, N. Suzuki and T. Motokawa (Eds.). A.A. Balkema, Rotterdam, pp: 105-111


  • Uehara, T., H. Asakura and Y. Arakaki, 1990. Fertilization Blockage and Hybridization among Species of Sea Urchins. In: Advances in Invertebrate Reproduction, Hoshi, M. and O. Yamashita (Eds.). Elsevier, Amsterdam, pp: 305-310


  • Uehara, T., H. Asakura and Y. Arakaki, 1990. Fertilization Blockage and Hybridization among Species of Sea Urchins. In: Advances in Invertebrate Reproduction, Hoshi, M. and O. Yamashita (Eds.). Elsevier, Amsterdam, pp: 305-310


  • Uehara, T., M. Shingaki, K. Taira, Y. Arakaki and H. Nakatomi, 1991. Chromosome Studies in Eleven Okinawan Species of Sea Urchins, with Special Reference to four Species of the Indo-Pacific Echinometra. In: Biology of Echinodermata, Yanagisawa, T., I. Yasumasu, C. Oguro, N. Suzuki and T. Motokawa (Eds.). A.A. Balkema, Rotterdam, pp: 119-129


  • Vacquier, V.D., 1980. The Isolation of Gamete Surface Components Involved in the Adhesion of Sperm to Eggs During Sea Urchin Fertilization. In: Advances in Invertebrates Reproduction Clark, W.H. and T.S. Adams (Eds.). Elsevier, North Holland, New York, pp: 215-224


  • Campton, D.E., 1987. Natural Hybridization and Introgression in Fishes. Methods of Detection and Genetic Interpretations. In: Population Genetics and Fishery Management Ryman, N. and F. Utter (Eds.). University of Washington Press, Seattle, Washington, USA., pp: 161-192


  • Campton, D.E., 1990. Application of Biochemical and Molecular Genetic Markers to Analysis of Hybridization. In: Electrophoretic and Isoelectric Focussing Techniques in Fisheries Management. Whitmore, D.H. (Ed.). CRC Press, Boca Raton, Florida, USA., pp: 241-263


  • Metz, E.C., H. Yanagimachi and S.R. Palumbi, 1991. Gamete Compatibility and Reproductive Isolation of Closely Related Indo-Pacific Sea Urchins, Genus Echinometra. In: Biol. Echinodermata, Yanagisawa, T., I. Yasumasu, C. Oguro, N. Suzuki and T. Motokawa (Eds.). A.A. Balkema, Rotterdam, pp: 131-137


  • Minor, J.E., B.B. Gao and E. Davidson, 1989. The Molecular Biology of Bindin. In: The Molecular Biology of Fertiliztion. Schatten, H. and G. Schatten (Eds.). Academic Press, San Diego, pp: 73-88


  • Nishihira, M., Y. Sato, Y. Arakaki and M. Tsuchiya, 1991. Eological Distribution and Habitat Preference of Four Types of Echinometra Mathaei on Okinawan Coral Reef. In: Biology of Echinodermata, Yanagisawa, T., I. Yasumasu, C. Oguro, N. Suzuki and T. Motokawa (Eds.). A.A. Balkema, Rotterdam, pp: 91-104


  • Palumbi, S.R., 1998. Species Formation and the Evolution of Gamete Recognition Loci. In: Endless Forms: Species and Speciation. Howard, D.J. and S.H. Berlocher (Eds.). Oxford University Press, New York, pp: 271-278

  • © Science Alert. All Rights Reserved