|
|
||||||||
Population Biology |
2Department of Biology, Faculty of Science, Chiba University, Yayoi, Inage, Chiba 263-8522, Japan; 3Department of Forest Genetics, Forestry and Forest Products Research Institute, Kukizaki, Ibaraki 305-8687, Japan
Received for publication April 1, 2003. Accepted for publication August 8, 2003.
| ABSTRACT |
|---|
|
|
|---|
Key Words: EST primers hybrid index hybrid zone introgression Japan Pinaceae Pinus SSCP
| INTRODUCTION |
|---|
|
|
|---|
Hybrid zones between Pinus pumila (Pallas) Regel and P. parviflora Siebold et Zucc. var. pentaphylla (Mayr) Henry have several interesting characteristics suitable for the study of interspecific genetic exchange and its evolutionary significance. First, Pinus species exhibit a paternal chloroplast inheritance and a maternal mitochondrial inheritance (Neale and Sederoff, 1989
; Mogensen, 1996
). Therefore, two cytoplasmic genomes can transgress independently across species boundaries. In fact, previous papers on the hybrid zones between P. pumila and P. parviflora var. pentaphylla indicated that the paternal chloroplast DNA (cpDNA) flowed from P. parviflora var. pentaphylla to P. pumila, and, in contrast, the maternal mitochondrial DNA (mtDNA) flowed from P. pumila to P. parviflora var. pentaphylla (Watano et al., 1995
, 1996
; Senjo et al., 1999
). The independent introgression of the cpDNA and mtDNA provides a unique opportunity to study the relationship among modes of inheritance and the levels and patterns of interspecific gene flow. Two distinct means of genetic movement, pollen and seed, are employed in seed plants. The paternal cpDNA is transmitted by both pollen and seed, and in contrast, the maternal mtDNA is carried only by seed. Therefore, the patterns of introgression of the cpDNA and mtDNA can be also viewed as a window to assess pollen- and seed-mediated gene flow across species boundary. Second, several hybrid zones in Japan seem to have formed independently. Pinus pumila is a creeping shrub, which dominates the vegetation zone above the forest limit of the high mountains in Japan, while P. parviflora var. pentaphylla is a tall tree in subalpine to montane zones. Hybrid zones are observed in the ecotonal zone of the two species in several mountains. Because the populations of P. pumila have an island-like geographic distribution due to the alpine habitat of the species, hybrid zones in different mountain regions seems to have been formed independently. A comparison of independent hybrid zones will allow us to elucidate the necessary genetic consequence of hybridization. Finally, the genomics project for Pinus taeda L. is ongoing (Brown et al., 2001
; Temesgen et al., 2001
), and the database of expressed sequence tags (ESTs) is available to the public (http://dendrome.ucdavis.edu/). The ESTs developed in P. taeda are also available in other pine species (Brown et al., 2001
).
In the previous studies on the introgressive patterns of cpDNA and mtDNA (Watano et al., 1995
, 1996
; Senjo et al., 1999
), we examined cpDNA and mtDNA haplotypes and morphology in each individual sampled and then assessed the direction and level of cpDNA and mtDNA introgressions based on morphological criteria. Although P. pumila and P. parviflora var. pentaphylla can be discriminated by many morphological characters such as cones, seeds, and needle anatomy (Isii, 1941
), the previous studies used only needle anatomical characters. This is because plants that have cones with seeds are relatively rare in nature. In the present study, we developed codominant nuclear markers, then elucidated the genetic structure of the hybrid zones by estimating the molecular hybrid index (MHI) for each individual. The MHI measures are much more unbiased criteria for the genetic makeup of hybrids than the morphological one with limited characters (Carney et al., 2000
). The combination of the MHI and cytoplasmic genome composition allowed us to compare the introgressive patterns of three independently inherited genomes: the biparental nuclear, paternal chloroplast, and maternal mitochondrial DNAs.
To develop nuclear DNA markers in our materials, we used polymerase chain reaction (PCR) primers developed from P. taeda ESTs (Temesgen et al., 2001
). First, we tested the cross-amplification of P. taeda PCR primers in P. pumila and P. parviflora var. pentaphylla, which belong to a different subgenus from P. taeda (Price et al., 1998
). Next, the allelic variation of the amplified DNA fragments was analyzed using the single-strand conformation polymorphism (SSCP) method. The PCR-SSCP of nuclear DNAs generates codominant markers (Bodenes et al., 1996
; Ishikawa et al., 2002
). Diagnostic or mostly diagnostic markers of P. pumila and P. parviflora var. pentaphylla were selected, and the multilocus genotype of each sample were determined. Finally, the genetic structure of the hybrid zones was described by comparing the MHI of each sample to previous data from the cytoplasmic genomes.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Development of diagnostic nuclear markers by PCR-SSCP
DNA from each plant was extracted from fresh leaves by the cetyltrimethylammonium bromide (CTAB) method (Doyle and Doyle, 1987
). To test cross-amplification of the P. taeda primers developed by Temesgen et al. (2001)
in our materials, we carried out PCR amplifications using 93 primer sets with 12 DNA templates (six P. pumila and six P. parviflora var. pentaphylla). The PCR reaction was conducted in a total volume of 25 µL containing 25 ng of template DNA, 0.2 µmol/L of each primer, 0.2 mmol/L of each dNTP, 1x PCR buffer, 1.5 mmol/L MgCl2, and 0.625 unit of TaKaRa ExTaq DNA polymerase (TaKaRa Bio., Tokyo, Japan). The thermal profile for PCR was as follows: initial 3 min denaturation at 95°C, three cycles at 95°C for 1 min, 58°C for 1 min, 72°C for 2 min, followed by three cycles at 95°C for 1 min, 55°C for 1 min, 72°C for 2 min, followed by 34 cycles at 95°C for 45 s, 52°C for 45 s, 72°C for 1 min 30 s, and a final extension at 72°C for 10 min. The samples obtained from the PCR reaction were resolved by 2% agarose gel, and the primer sets that generated the expected size of DNA fragments with no artificial bands were screened. The sequence variation of PCR products amplified in P. pumila and P. parviflora var. pentaphylla were then analyzed by SSCP. A portion of PCR sample (1 µL) was mixed with 19 µL of formamide-dye solution (90% formamide, 0.005% bromophenol blue, 8% glycerol), and then denatured for 3 min at 95°C. The denatured samples were cooled on ice, and 4 µL of the sample was loaded on a nondenatured 50% MDE gel (TaKaRa Bio., Shiga, Japan) with dimensions of 135 mm wide, 130 mm long, and 0.75 mm thick. Electrophoresis was carried out at two electrophoretic conditions (gels with 2% and 5% glycerol) at 20°C in 50% TBE (50 mmol/L Tris, 41.5 mmol/L boric acid, and 1 mmol/L EDTA-2Na) using an electrophoretic apparatus with a thermostat-controlled cooled water circulator (AE-6290, ATTO, Tokyo, Japan). The DNA bands were visualized using a DNA Silver Staining Kit (Amersham Biosciences, Piscataway, New Jersey, USA). Finally, we selected primer sets for which the PCR products had diagnostic SSCP band patterns for P. pumila and P. parviflora var. pentaphylla.
Molecular hybrid index
Some of the marker loci used in this study were not perfectly diagnostic. To calculate the hybrid index of each plant, therefore, we employed the maximum-likelihood method of Rieseberg et al. (1998)
. Here, we hypothesize that a hybrid individual has a hybrid index h, which represents the overall proportion of alleles in the hybrid individual inferred to be derived from P. pumila. Under this hypothesis, the probability that the hybrid individual has an allele x is
![]() | (1) |
![]()
| RESULTS |
|---|
|
|
|---|
|
|
|
Genetic makeup of hybrid zones
Nuclear genotypes at eight marker loci were determined for 46 individuals from Mt. Asahidake and 69 from Mt. Higashiazuma. The multilocus genotypes are shown in Appendix 1 (see Supplemental Data accomanying the online version of this article) together with data on the altitude, horizontal distance from the first sample collected, cytoplasmic haplotypes, and needle morphology. No apparent difference of introgressive patterns was observed among eight marker loci (Appendix 1). The changes in the MLEs of the hybrid index and in the haplotypes of the cpDNA and mtDNA along the sampling routes are illustrated in Fig. 2. Because the MHI represents the overall proportion of alleles in the hybrid individual inferred to be derived from P. pumila, a high MHI value corresponds to a P. pumila-like plant and a low MHI to a P. parviflora-like plant.
|
It is clear that the values of the hybrid index at Mt. Higashiazuma (Fig. 2b) rapidly changed between sample numbers AZ43 (hybrid index = 0.016) and AZ44 (0.875). In this mountain, a gap in the distribution of five-needle pines existed between AZ44 and AZ45: 125 m in altitude and 350 m in distance. This gap was covered with dense forest of Abies mariesii Masters and Tsuga diversifolia Masters. The values of the hybrid index below the Abies and Tsuga forest zone varied from 0 to 0.203 with a mean of 0.075 except for AZ44, while the hybrid index above the zone varied from 0.173 to 1 with a mean of 0.945. Only one individual (AZ58) at the top of mountain had a low value of the hybrid index (0.173). That individual was specifically collected because it was the only individual that had an erect stem, which is a key character of P. parviflora var. pentaphylla.
Comparison of the levels of nuclear introgression with cpDNA and mtDNA introgression
To compare nuclear introgression with cpDNA and mtDNA, we divided samples into several zones mainly based on altitudinal discontinuity and calculated the mean value of the hybrid index and the frequencies of the cpDNA and mtDNA haplotypes in each zone (Table 3). As mentioned in the introduction, the introgression of cytoplasmic genomes are unidirectional in hybrid zones of P. pumila and P. parviflora var. pentaphylla: cpDNA flowed from P. parviflora var. pentaphylla to P. pumila (from the bottom to the top), and, in contrast, mtDNA flowed from P. pumila to P. parviflora var. pentaphylla (from the top to the bottom) (Watano et al., 1996
; Senjo et al., 1999
).
|
On Mt. Higashiazuma, in contrast to Mt. Asahidake, the introgression of P. pumila mtDNA to the lower zones was prominent; the frequencies of the P. pumila mtDNA haplotype was much higher than the values of the hybrid index in Bottom 1 and 2. The level of introgression of P. parviflora var. pentaphylla cpDNA to the highest zone (Top) was low (0.120) and was nearly identical to the level of the nuclear portion of P. parviflora var. pentaphylla (0.055 = 1 0.945).
Molecular hybrid index and needle morphology
Watano et al. (1995)
and Senjo et al. (1999)
classified each sample into three categories, the P. pumila type, the P. parviflora type, and the intermediate type, based on the anatomical characters of needles. The characters used were (1) the number and position of resin canals, (2) the presence or absence of sclerenchymatous cells above and below a fibrovascular bundle, and (3) the presence or absence of idioblasts. Table 4 summarizes the relationship among morphological categories and the hybrid index.
|
As for Mt. Higashiazuma, morphological data were available for all samples (Senjo et al., 1999
). Sixteen intermediate types were classified into two different types. The type with an intermediate resin canal arrangement had low values of hybrid index with the mean of 0.084, which did not significantly deviate from the mean value (0.051) of the P. parviflora type (P = 0.17). All samples of this type, except AZ58, were below the gap of distribution between AZ44 and AZ45 (Fig. 1b). The other intermediate type had high values of the hybrid index, with the mean of 0.950, and all samples except AZ44 grew above the gap.
| DISCUSSION |
|---|
|
|
|---|
Although Mt. Asahidake lacks a subalpine coniferous forest, a gap in the pine distribution exists at the lower region, from 1300 to 1600 m in elevation. Nuclear and mtDNA introgression from P. pumila were stopped by the gap, while nuclear and cpDNA introgression from P. parviflora var. pentaphylla transgressed the gap. It is interesting to note that the effects of the distribution gap on gene flow are different in the two hybrid zones. One possible explanation is the difference in the altitudinal locations of the distribution gaps on these two mountains. The gap at the high elevation (Mt. Higashiazuma) may be more effective in stopping the upward gene flow from P. parviflora var. pentaphylla, and the gap at the low elevation (Mt. Asahidake) may be a more effective barrier for the downward gene flow from P. pumila.
The seeds of P. pumila have no wings and are dispersed mainly by the Eurasian Nutcracker (Nucifraga caryocatactes) (Kajimoto, 2002
). The dispersal of pine seeds by nutcrackers and the range of seed dispersal may be critical to understanding detailed gene flow patterns in hybrid zones. In this respect, it should be noted that only AZ58 had a low hybrid index (0.173) in the top zone of Mt. Higashiazuma. All plants surrounding AZ58 had high hybrid indices over 0.9. Because the hybrid index of the seeds produced by the plants in the top zone should be higher than 0.45, AZ45 cannot be the direct offspring of the plants in the top zone. AZ58 probably originated from seed dispersed from the zones below the gap. This indicates that the transport of pine seeds by nutcrackers may occasionally become a medium for gene flow between the zones above and below the gap. As a whole, however, the contrasting hybrid index values above and below the gap in Mt. Higashiazuma suggest that such occasional seed flow has been ineffective in mixing the two zones.
Nuclear vs. cytoplasmic introgression
The paternal cpDNA had a higher level of introgression than the nuclear DNAs at Mt. Asahidake, while the maternal mtDNA introgression exceeded nuclear introgression at Mt. Higashiazuma (Fig. 2, Table 3). These results indicated that the levels of introgression of the two differently inherited cytoplasmic genomes were largely affected by unique situations on each mountain. As discussed earlier, the nature and the location of the distribution gap in the hybrid zone could be one of the factors determining the introgressive pattern of each genome. In addition to the present situation on each mountain, historical factors should also be considered. An excess of cpDNA introgression over nuclear introgression has been reported in many angiosperm species (Rieseberg and Wendel, 1993
). Because cpDNA is maternally inherited in most angiosperms, the examples of extensive cpDNA introgression in angiosperms may correspond to extensive mtDNA introgression, which is the case with Mt. Higashiazuma in our study. Senjo et al. (1999)
examined mtDNA haplotypes in many populations throughout the distribution range of P. parviflora var. pentaphylla and found that the populations that had high frequencies of the mtDNA haplotype of P. pumila were clustered within the southern and middle parts of the Ohu Mountains in Tohoku. Mt. Higashiazuma is located in the southernmost position of an extensive mtDNA introgression area. Several authors (Dong and Wagner, 1993
; Martinsen et al., 2001
) discussed that such extensive introgression of cytoplasmic genomes may reflect ancient introgression, namely, a footprint of the ancient location of the hybrid zone. It is not certain that the ancient introgression hypothesis holds true in our study. However, it may be a good working hypothesis for further studies.
Molecular hybrid index and morphology
Based on the comparison of the hybrid index for categories of needle morphology, needle morphology was a poor indicator for assessing the nuclear genetic constitution of individuals in hybrid zones. For example, on Mt. Asahidake, there were no significant differences in mean values of hybrid index among the three morphological categories for samples above 1600 m in elevation (Table 4). Carney et al. (2000)
examined both the morphological and molecular hybrid index in a hybrid population between Helianthus annuus and H. bolanderi. The hybrid population was divided into western and eastern halves by a central patch of grass. In the eastern half, there was no significant correlation between the morphological and molecular hybrid index measures, and morphologically pure H. annuus-like plants contained a considerable portion of H. bolanderi genes (approximately 20%). The authors suggested that this might be the result of the ecological selection on hybrids for H. annuus morphological traits. Although our morphological data are limited, the situation observed in the hybrid population of Helianthus may be very similar to that in the pine hybrid zone on Mt. Asahidake. Pinus pumila is adapted to severe environments in the alpine zone. Although the pines of the highest region on Mt. Asahidake had a highly introgressed nuclear genome (25%), they still survive in a "pure" P. pumila habitat. It is possible that a strong conflict between gene flow and natural selection occurred in the past and may still be ongoing. A detailed description of the morphological and physiological characters of this mountain population is needed.
Summary and conclusions
In this study, we developed codominant nuclear DNA markers diagnostic or mostly diagnostic for P. pumila and P. parviflora var. pentaphylla. The molecular hybrid index measures calculated by these nuclear markers clarified the genetic structure of the hybrid zones on two mountains, Mt. Higashiazuma and Mt. Asahidake. On Mt. Higashiazuma, a P. parviflora-like (low hybrid index) zone and a P. pumila-like (high hybrid index) zone were clearly recognized below and above the subalpine coniferous forest zone located in the altitudinal range from 1800 to 1900 m. Mt. Asahidake lacks the subalpine coniferous forest, and the hybrid zone at this mountain showed a gradual increase in the hybrid index from the bottom to the top. The contrasting genetic structure of the two mountains suggested that even a small distribution gap (approximately 100 m in elevation) could play a role in creating an effective barrier to nuclear gene flow in the hybrid zone. As for the relationship among inheritance modes of genomes and the levels of interspecific gene flow, the two hybrid zones had different trends. The paternal cpDNA introgression was more extensive than the nuclear introgression on Mt. Asahidake, but not on Mt. Higashiazuma. In contrast, the maternal mtDNA introgression exceeded the nuclear introgression on Mt. Higashiazuma, but not on Mt. Asahidake. It is possible that the level and pattern of introgression of each genome could be affected by unique circumstances on each mountain, such as the size and position of the distribution gap.
| FOOTNOTES |
|---|
4 watano{at}faculty.chiba-u.jp ![]()
| LITERATURE CITED |
|---|
|
|
|---|
Banks M. A. W. Eichert 2000 WHICHRUN (version 3.2): a computer program for population assignment of individuals based on multilocus genotype data. Journal of Heredity 91: 87-89
Barton N. H. K. S. Gale 1993 Genetic analysis of hybrid zones. In G. G. Harrison [ed.], Hybrid zones and the evolutionary process. Oxford University Press, New York, New York, USA
Bodenes C. F. Laigret A. Kremer 1996 Inheritance and molecular variations of PCR-SSCP fragments in pedunculate oak (Quercus robur L). Theoretical and Applied Genetics 93: 348-354[CrossRef][ISI]
Brown G. R. E. E. Kadel D. L. Bassoni K. L. Kiehne B. Temesgen J. P. van Buijtenen M. M. Sewell K. A. Marshall D. B. Neale 2001 Anchored reference loci in loblolly pine (Pinus taeda L.) for integrating pine genomics. Genetics 159: 799-809
Carney S. E. K. A. Gardner L. H. Rieseberg 2000 Evolutionary changes over the fifty-year history of a hybrid population of sunflowers (Helianthus). Evolution 54: 462-474[CrossRef][ISI][Medline]
Dong J. D. B. Wagner 1993 Taxonomic and population differentiation of mitochondrial diversity in Pinus banksiana and Pinus contorta. Theoretical and Applied Genetics 86: 573-578[CrossRef][ISI]
Doyle J. J. J. L. Doyle 1987 A rapid DNA isolation procedure for small quantities of fresh leaf material. Phytochemical Bulletin 19: 11-15
Grant V. 1981 Plant speciation. Columbia University Press, New York, New York, USA
Ishikawa H. Y. Watano K. Kano M. Ito S. Kurita 2002 Development of primer sets for PCR amplification of PgiC gene in ferns. Journal of Plant Research 115: 65-70[CrossRef][ISI][Medline]
Isii S. 1941 On the various forms of Pinus pumila and other northern Japanese soft pines with special reference to their distribution (IV). Journal of Japanese Forestry Society 23: 47-55
Kaji M. 1982 Studies on the ecological geography of subalpine conifers; distribution pattern of Abies mariesii in relation to the effect of climate in the postglacial warm period. Bulletin of the Tokyo University Forests 72: 31-120
Kajimoto T. 2002 Factors affecting seedling recruitment and survivorship of the Japanese subalpine stone pine, Pinus pumila, after seed dispersal by nutcrackers. Ecological Research 17: 481-491[CrossRef][ISI]
Martinsen G. D. T. G. Whitham R. J. Turek P. Keim 2001 Hybrid populations selectively filter gene introgression between species. Evolution 55: 1325-1335[CrossRef][ISI][Medline]
Mogensen H. L. 1996 The hows and whys of cytoplasmic inheritance in seed plants. American Journal of Botany 83: 383-404[CrossRef][ISI]
Neale D. B. R. R. Sederoff 1989 Paternal inheritance of chloroplast DNA and maternal inheritance of mitochondrial DNA in loblolly pine. Theoretical and Applied Genetics 77: 212-216[ISI]
Price R. A. A. Liston S. Strauss 1998 Phylogeny and systematics of Pinus. In D. M. Richardson [ed.], Ecology and biogeography of Pinus, 4968. Cambridge University Press, Cambridge, UK
Rattenbury J. A. 1962 Cyclic hybridization as a survival mechanism in the New Zealand forest flora. Evolution 16: 348-363[CrossRef][ISI]
Rieseberg L. H. 1995 The role of hybridization in evolution: old wine in new skins. American Journal of Botany 82: 944-953[CrossRef][ISI]
Rieseberg L. H. S. J. E. Baird A. M. Desrochers 1998 Patterns of matings in wild sunflower hybrid zones. Evolution 52: 713-726[CrossRef][ISI]
Rieseberg L. H. J. F. Wendel 1993 Introgression and its consequences in plants. In R. G. Harrison [ed.], Hybrid zones and the evolutionary process, 70109. Oxford University Press, New York, New York, USA
Rundel P. W. B. J. Yoder 1998 Ecophysiology of Pinus. In D. M. Richardson [ed.], Ecology and biogeography of Pinus, 296323. Cambridge University Press, Cambridge, UK
Senjo M. K. Kimura Y. Watano K. Ueda T. Shimizu 1999 Extensive mitochondrial introgression from Pinus pumila to P. parviflora var. pentaphylla (Pinaceae). Journal of Plant Research 112: 97-105[CrossRef][ISI]
Temesgen B. G. R. Brown D. E. Harry C. S. Kinlaw M. M. Sewell D. B. Neale 2001 Genetic mapping of expressed sequence tag polymorphism (ESTP) markers in loblolly pine (Pinus taeda). Theoretical and Applied Genetics 102: 664-675[CrossRef][ISI]
Watano Y. M. Imazu T. Shimizu 1995 Chloroplast DNA typing by PCR-SSCP in the Pinus pumila-P. parviflora var. pentaphylla complex (Pinaceae). Journal of Plant Research 108: 493-499[CrossRef][ISI]
Watano Y. M. Imazu T. Shimizu 1996 Spatial distribution of cpDNA and mtDNA haplotypes in a hybrid zone between Pinus pumila and P. parviflora var. pentaphylla (Pinaceae). Journal of Plant Research 109: 403-408[CrossRef][ISI]
This article has been cited by other articles:
![]() |
G. Aldridge and D. R. Campbell Asymmetrical pollen success in Ipomopsis (Polemoniaceae) contact sites Am. J. Botany, June 1, 2006; 93(6): 903 - 909. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Ebihara, H. Ishikawa, S. Matsumoto, S.-J. Lin, K. Iwatsuki, M. Takamiya, Y. Watano, and M. Ito Nuclear DNA, chloroplast DNA, and ploidy analysis clarified biological complexity of the Vandenboschia radicans complex (Hymenophyllaceae) in Japan and adjacent areas Am. J. Botany, September 1, 2005; 92(9): 1535 - 1547. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |