|
|
||||||||
1Department of Biological Sciences, Wayne State University, Detroit, Michigan 48202; 2USDA Forest Service Shrub Sciences Laboratory, 735N 500E, Provo, Utah 84606; and 4Institute of Neuroscience, 1254 University of Oregon, Eugene, Oregon 97403
Received for publication January 6, 1998. Accepted for publication January 14, 1999.
| ABSTRACT |
|---|
|
|
|---|
Key Words: Artemisia tridentata Asteraceae biological absorption coefficients elemental concentrations hybrid zone niche differentiation reciprocal transplant experiment sagebrush
| INTRODUCTION |
|---|
|
|
|---|
For the past several decades, considerable attention has been given to explaining the persistence of stable hybrid zones, because their existence violates adaptive speciation theory and calls the biological species concept into question (see Harrison [1993]
for a discussion of the importance of stable hybrid zones). While three models have been proposed to account for the existence of stable hybrid zones, they fall into two basic groups. Both the dynamic equilibrium and mosaic hybrid zone models assume that hybrids are relatively unfit when compared to the parental taxa (Barton and Hewitt, 1985
; Harrison, 1986
, 1990
; Howard, 1986
; Hewitt, 1988
; Harrison and Rand, 1989
; Arnold, 1997
). This unfitness is due to the disruption of coadapted gene complexes and is thus independent of the environment. These models differ from one another in that the mosaic hybrid zone model assumes that the parental taxa have adapted to different environments, whereas the dynamic equilibrium model is ecologically neutral (Arnold, 1997
). In contrast, the bounded hybrid superiority model assumes that hybrids are more fit than both parental taxa within the hybrid zone, but less fit than the parental taxa when raised in either parental zone. Accordingly, genotype-by-environment interactions are believed to stabilize the hybrid zone, with the constraint that hybrids have higher fitness only in the hybrid zone. (Note that genotype-by-environment interactions could cause decreased hybrid fitness as presumed by the other two models.)
The factors determining the width and position of the hybrid zone also differ between the two sets of models. According to the dynamic equilibrium and mosaic hybrid zone models, the balance between endogenous selection against hybrids and gene flow determines the width of the zone. Because the dynamic equilibrium model lacks any ecological term there is nothing to link the zone to a particular habitat. Indeed, the zone is believed to move until it encounters an area of low population density where it becomes trapped (Barton and Hewitt, 1985
). In order for the zone to move, the parental taxa must be adapted to the new environment. Thus, this model predicts complete niche overlap by the parental taxa at least within the vicinity of the hybrid zone. This provides another prediction that can be examined independent of hybrid fitness considerations and is the subject of this study. The mosaic hybrid zone assumes only that the parental taxa are adapted to their respective environments, thus the hybrid zone must fall between them with the width being determined by the balance between endogenous selection and gene flow (see Arnold, 1997
).
In contrast, the bounded hybrid superiority model assumes that ecological selection via genotype-by-environment interactions determines both the position and width of the zone. Consequently, it is the nature of the environment that determines the width of the hybrid zone (for a discussion of these models, see Moore, 1977
; Moore and Buchanan, 1985
; Moore and Price, 1993
; Arnold, 1997
). Accordingly, the parental taxa and the hybrids are adapted to their respective environments and must show niche differentiation.
The critical issues that distinguish these models concern the relative fitness of the hybrids (particularly in the hybrid zone), the nature of selection constraining gene flow, and the role of the environment in determining the position and width of the hybrid zone. We have previously resolved fitness issues for the big sagebrush hybrid zone in northern Utah. Using reciprocal transplant experiments, we found significant genotype-by-environment interactions for major life history features (germination, survivorship, growth, and fecundity) in our study of the big sagebrush hybrid zone (Graham, Freeman, and McArthur, 1995
; Wang et al., 1997
). Hybrids were significantly more fit than parental taxa within the hybrid zone and less fit than the native parental taxon when transplanted outside the hybrid zone (Wang et al., 1997
). These findings are clearly consistent with the bounded hybrid superiority model but conflict with the predictions of the other models.
Using the elemental concentration of leaves as indicators of niche, we test another prediction of these models, that of niche overlap/differentiation. The bounded hybrid superiority model posits that each parental taxon is adapted to its native environment, and the hybrids are adapted to the hybrid zone, i.e., they have different niches, or at least exhibit some degree of niche separation (Moore, 1977
; Moore and Buchanan, 1985
; Moore and Price, 1993
). In contrast, the dynamic equilibrium model assumes that within the vicinity of the hybrid zone, the taxa do not differentially respond to environmental heterogeneity, i.e., the parental taxa and hybrids should have complete niche overlap. The mosaic hybrid zone presumes that the parental taxa have different niches, but makes no direct prediction about their niche overlap with the hybrids, and thus is not a subject of this study.
Plants make ideal experimental subjects to test these models because all plants require the same 17 essential nutrients (Salisbury and Ross, 1992
). Most of these nutrients are absorbed directly from the soil (Foth, 1992
). Elemental concentrations found in various plant parts can serve as indicators of plant responses to edaphic differences (Shkolnik, 1984
; Jeffrey, 1987
; Rendig and Taylor, 1989
), i.e., they can be used to determine whether there is niche separation or overlap. If elemental concentrations in plants raised in a common garden do not differ significantly among the taxa, then niche overlap is indicated, as suggested by the dynamic equilibrium model. However, if elemental concentrations in plants do differ significantly, particularly when raised in the same garden, then this would indicate at least a partial niche separation, as assumed by the bounded hybrid superiority model.
We examined elemental concentrations of the leaves of basin big sagebrush (Artemisia tridentata ssp. tridentata), mountain big sagebrush (A. tridentata. ssp. vaseyana), and their hybrids in native stands and reciprocally transplanted seedlings raised in common gardens. In particular, we address two questions: (1) Does a taxon have the same elemental concentrations when grown in different gardens, i.e., does it respond to the edaphic differences among gardens? and (2) Do the different taxa exhibit the same elemental concentrations when raised in the same garden? The null hypotheses correspond to the predictions of the dynamic equilibrium model. Significant genotype-by-environment interactions allow us to reject both null hypotheses and, thus, accept the alternative bounded hybrid superiority model. Note that we can accept the bounded hybrid superiority model in this case only because we have previously shown that there are significant genotype-by-environment interactions for fitness components and that hybrids have the highest relative fitness in the hybrid zone (Wang et al., 1997
).
| MATERIALS AND METHODS |
|---|
|
|
|---|
Sample collection
Leaves were collected from a total of 130 plant samples: 40 from the five native stands (eight samples per stand) and 90 from the seedling transplant experiments (six leaf samples per seedling source per common garden). For the native plants (hereafter referred to as site indigenous),
100 g of leaves were collected from each of eight adult plants per stand. To minimize the impact on the transplanted seedlings, only 5 g of leaves were collected from one seedling within each half-sib family in the common gardens.
Sample treatment and digestion
Samples were air dried for at least 30 d, then pulverized using a grinding mill, and temporarily stored in paper bags. A modified dry-ashing procedure (Munter and Grande, 1983
) was used: all pulverized samples were redried in an oven at 85°C for 10 h and stored in a desiccator until weighed. Three-tenths of a gram of each sample were transferred to a 10-mL crucible and ashed in a furnace at 485°C for 14 h. The ashed sample was first soaked with 1 mL of concentrated nitric acid (75%) for 1 h, then transferred to a 50-mL centrifuge tube and diluted to 15-mL with diluted nitric acid (2%). This solution was used to determine elemental concentrations of B, Ca, Cu, Fe, K, Mg, Mn, Mo, Na, and Zn, using a Leeman Lab inductively coupled plasma emission spectrophotometer.
Previously, we had determined the concentrations of eight of these elements (Ca, Cu, Fe, K, Mg, Mo, Na, and Zn) in soils from the same area (Wang, 1996
; Wang et al., 1998
). The majority of these elements (B, Ca, Cu, Fe, K, Mg, Mn, Mo, and Zn) are essential to angiosperms. Sodium is often required by plants in arid regions (Brownell, 1979
; Salisbury and Ross, 1992
).
Analytic precision
Duplicate digestions were conducted for each plant sample, and the concentrations of each element were measured twice per digestion. Further, one empty control was done for each batch of digestions. The corresponding control value (if any) was subtracted from each measurement. All concentrations are expressed in milligrams per kilogram of dry matter.
Data analysis
Because we collected multiple data points from the same plant, we used a MANOVA (Zar, 1984
; Norusis, 1988
). We used univariate ANOVAs and post hoc Student-Newman-Keuls multiple comparison tests to analyze for differences among means. We examined the relationships between elemental concentrations of plants and soils in two ways. First, correlation analyses were conducted between the concentrations of the same element in plants and soils, and secondly, we calculated the Biological Absorption Coefficient (BAC) for each element (Brooks, 1983
). BAC = Cp /Cs, where Cp is the concentration of the element in plants and Cs is the concentration of the same element in soils. BACs have been used to examine the differential response of plants (Shkolnik, 1984
; Jeffrey, 1987
) and thus are appropriate for our purposes.
Correspondence analysis (CA) (CANOCO; Ter Braak, 1992
; Randerson, 1993
) was used to analyze the pooled elemental concentrations and BACs, following a square-root transformation (to improve normality), for leaves from both site-indigenous plants and seedlings in the reciprocal transplant gardens. We used CA rather than principal components analysis because CA better reflects changes across gradients (Randerson, 1993
). Sample scores of the first two axes from CA were used for ordination. To conserve space we only report the result for the BACs. The patterns for the elemental concentrations accounted for a lower percentage of the variance, but were remarkably similar to those for the BACs (Wang, 1996
).
| RESULTS |
|---|
|
|
|---|
|
|
|
|
Biological Absorption Coefficients (BACs) of elements by site-indigenous plants and transplanted seedlings
The biological absorption coefficients of the elements in the site-indigenous taxa and transplanted seedlings have the following order: BACMo > BACNa
BACK > BACZn
BACMg
BACFe2 > BACCu > BACCa. BACs of the eight elements are all greater than one, suggesting that the big sagebrush taxa concentrate these elements. The MANOVA showed a significant site effect (Wilks' lambda = 0.0006, P << 0.001). Univariate analyses showed that BACK, BACMg, BACFe, BACCu, BACMo, and BACNa differed significantly among plants from the five native sites, while BACCa and BACZn did not (Table 3).
For the transplanted seedlings, the MANOVA showed that BACs differed significantly among seedling origins (Wilks' lambda = 0.296, P < 0.001) and among the gardens (Wilks' lambda = 0.005, P < 0.001). There was also a significant garden by seedling origin interaction Wilks' lambda = 0.193, P < 0.001). Univariate analyses showed that all BACs differed significantly among the gardens, while BACCa, BACFe, BACZn, and BACNa varied significantly among seedling sources and BACFe and BACNa had significant garden site-by-seedling source interactions (Fig. 2, Table 4).
|
|
|
|
|
|
Correlations between elemental concentrations and plant size
We previously measured a number of parameters on the same plants that we used in this study and thus can correlate the elemental concentrations reported here with these prior measurements. However, in doing so we are limited by the degrees of freedom. Consequently, we report only correlations with plant size (the product of height and crown diameter) as we have the most degrees of freedom for that analysis. Overall, Fe (R = -0.25, P < 0.05) and Mn (R = -0.35, P < 0.05) were negatively correlated with plant size, while B (R = 0.17, P < 0.053) and Cu (R = 0.15, P < 0.073) showed weak positive correlations.
| DISCUSSION |
|---|
|
|
|---|
Given the essential nature of the elements studied, it is not surprising that BACs are all greater than one, suggesting that big sagebrush taxa concentrate them, but they do so differently. Thus, it is probably not surprising that linear correlations between soil and leaf concentrations were not strong (only four of the eight correlations were significant, and these did not account for a high fraction of the variance). Nevertheless, the concentrations of elements in the soils probably contribute to the structuring of the hybrid zone, because the relative contributions of the eight elements to the first axes in both the soil and leaf analyses are highly correlated.
Based upon our examination of life history features (Wang et al., 1997
), each parental taxon was most fit within its native habitat, and hybrids were overwhelmingly more fit than either parental taxon in the hybrid zone. Thus, our life history data also show significant genotype-by-environment interactions. Here, we found that the same taxon responded differentially when transplanted into different gardens and that the different taxa concentrate minerals to different extents. Thus, it is not surprising that the different taxa responded differentially when transplanted into the same garden. These findings indicate that genotype-by-environment interactions influence the concentrations and BACs of essential elements from the soil. However, because of our limited degrees of freedom we have not identified which elements are responsible for the fitness differentials. Nevertheless, we find it most interesting that Fe showed significant genotype-by-environment interactions, was a major component of the CA axes, and was negatively correlated with plant size. However, it may be naive to suppose that there should be a simple correlation between fitness and elemental concentration or BACs. Often it is the ratios or proportions among elements that influence the outcome of competition (Tilman, 1985
). While we have not identified which elements, or combination of elements, limit the fitness of each taxon in each environment and believe such studies would be very useful and informative, we do not think they are necessary to establish niche differentiation. What is necessary is to show that the taxa utilize essential aspects of the environment differently (see, for example, Lack, 1944
) which we have done.
We have shown that the parental and hybrid taxa have different niches as they differentially use available resources. That such microevolution (the width of the zone is
1 km) can occur despite high gene flow (big sagebrush is wind pollinated with wind-dispersed seeds) is testament to the strong selection required to stabilize this zone. We estimate selection coefficients on the order of 10-1 (Wang et al., 1997
). We raised seedlings from three sites within the hybrid zone in our common gardens. In no case was the common garden in the hybrid site >250 m from the native location of any hybrid seedling. Yet, the hybrids exhibited some degree of niche separation among plants from the three hybrid sites. This indicates that extremely powerful local selection is operating to structure this hybrid zone.
Our evidence indicates that each taxon is adapted to the local environment in which it occurs naturally. We have tried to identify the nature of that adaptation. Accordingly, we have also examined water stress and respiration rates of plants in the common gardens (McArthur et al., 1998a
). In that study we found significant differences among the gardens and taxa, but not genotype-by-environment interactions. Thus, the physiological evidence, at this time, indicates that adaptations to the different environments encountered across the hybrid zone probably involves the active uptake of the elements described in this study. Because fitness does not readily map onto the elemental concentrations or BACs, further work is required to elucidate the precise adaptive features.
Our study of soil chemistry (Wang et al., 1998
) indicates that the sites have markedly different soil chemistry and, within the hybrid zone, soils are spatially more variable than corresponding sites in either parental zone. It appears that hybridization provides the enhanced genetic variance required to successfully occupy such sites (see Lewontin and Birch, 1966
; Grant and Grant, 1994
).
While the taxa do differ for a number of traits (e.g., morphology, secondary compounds, and PCR markers; (McArthur, Welch, and Sanderson, 1988
; Freeman et al., 1991
; McArthur et al., 1998b
), they apparently have not diverged to the point that hybridization disrupts coadaptation (see Graham, 1992
). Freeman et al. (1995)
found no evidence of enhanced developmental instability among native hybrids compared to the parental taxa. Thus, it appears that strong exogenous selection stabilizes the sagebrush hybrid zone and sifts among the phenotypic variants, adapting hybrid populations to unique microhabitats within the hybrid zone.
| FOOTNOTES |
|---|
3 Author for correspondence (cfreeman@sun.science.wayne.edu). ![]()
| LITERATURE CITED |
|---|
|
|
|---|
Barton, N. H., and G. M. Hewitt. 1985 Analysis of hybrid zones. Annual Review of Ecology and Systematics 16: 113148.[CrossRef][ISI]
Brooks, R. R. 1983 Biological methods of prospecting for minerals. John Wiley & Sons, New York, NY.
Brownell, P. E. 1979 Sodium as an essential micro-nutrient element for plants and its possible role in metabolism. Advances in Botanical Research 7: 117224.
Burke, J. M., S. E. Carney, and M. L. Arnold. 1998 Hybrid fitness in the Louisiana irises: analysis of parental and F1 performance. Evolution 52: 3743.
Burton, R. S. 1990a Hybrid breakdown in physiological responses: a mechanistic approach. Evolution 44: 18061813.[CrossRef][ISI]
. 1990b Hybrid breakdown in developmental time in the copepod Tigriopus californicus. Evolution 44: 18141822.
Emms, S. K., and M. L. Arnold. 1997 The effect of habitat on parental and hybrid fitness: reciprocal transplant experiments with Louisiana iris. Evolution 51: 11121119.[CrossRef][ISI]
Fisher, R. A. 1930 Genetical theory of natural selection. Clarendon Press, Oxford.
Foth, H. D. 1990 Fundamentals of soil sciences, 8th ed. John Wiley & Sons, New York, NY.
Freeman, D. C., W. A. Turner, E. D. McArthur, and J. H. Graham. 1991 Characterization of a narrow hybrid zone between two subspecies of big sagebrush (Artemisia tridentata: Asteraceae). American Journal of Botany 78: 805815.[CrossRef][ISI]
, J. H. Graham, D. W. Byrd, E. D. McArthur, and W. A. Turner. 1995 Narrow hybrid zone between two subspecies of big sagebrush (Artemisia tridentata: Asteraceae). III. Developmental instability. American Journal of Botany 82: 11441152.[CrossRef][ISI]
Graham, J. H. 1992 Genomic coadaptation and developmental stability in hybrid zones. Acta Zoologica Fennica 191: 121132.
, D. C. Freeman, and E. D. McArthur. 1995 Narrow hybrid zone between two subspecies of big sagebrush (Artemisia tridentata: Asteraceae). II. Selection gradients and hybrid fitess. American Journal of Botany 82: 709716.[CrossRef][ISI]
Grant, P. R., and B. R. Grant. 1994 Phenotypic and genetic effects of hybridization in Darwin's finches. Evolution 48: 297316.[CrossRef][ISI]
Harrison, R. G. 1986 Patterns and process in a narrow hybrid zone. Heredity 56: 337349.[ISI]
. 1990 Hybrid zones: windows on evolutionary process. Oxford Surveys in Evolutionary Biology 7: 69128.
. 1993 Hybrids and hybrid zones: historical perspectives. In R. G. Harrison [ed.], Hybrid zones and the evolutionary process, 312. Oxford University Press, Oxford.
, and J. Arnold. 1982 A narrow hybrid zone between closely related cricket species. Evolution 36: 535552.[CrossRef][ISI]
, and D. M. Rand. 1989 Mosaic hybrid zones and the nature of species boundaries. In D. Otte and J. A. Endler [eds.], Speciation and its consequences, 11311133. Sinauer, Sunderland, MA.
Hewitt, G. M. 1988 Hybrid zonesnatural laboratories for evolutionary studies. Trends in Ecology and Evolution 3: 158167.[CrossRef]
Howard, D. J. 1986 A zone of overlap and hybridization between two ground cricket species. Evolution 40: 3443.[CrossRef][ISI]
Jeffrey, D. W. 1987 Soil-plant relationships: an ecological approach. Timber Press, Portland, OR.
Lack, D. 1944 Symposium on "The ecology of closely allied species." Journal of Animal Ecology 13: 176177.
Lewontin, R. C., and L. C. Birch. 1966 Hybridization as a source of variation for adaptation to new environments. Evolution 20: 315336.[CrossRef][ISI]
Mallet, J., and N. H. Barton. 1989a Strong natural selection in a warning-color hybrid zone. Evolution 43: 421431.[CrossRef][ISI]
McArthur, E. D., D. C. Freeman, J. H. Graham, H. Wang, S. C. Sanderson, T. A. Monaco, and B. N. Smith. 1998a Narrow hybrid zone between two subspecies of big sagebrush (Artemisia tridentata: Asteraceae) VI. Respiration and water potential. Canadian Journal of Botany 76: 567574.[CrossRef]
, J. E. Mudge, R. Van Buren, W. R. Anderson, S. C. Sanderson, and D. G. Babbel. 1998b Random amplified polymorphic DNA analysis (RAPD) of Artemisia subgenus Tridentatae species and hybrids. Great Basin Naturalist 58: 1227.
, B. L. Welch, and S. C. Sanderson. 1988 Natural and artificial hybridization between big sagebrush (Artemisia tridentata) subspecies. Journal of Heredity 79: 268276.
McKenzie, J. A., and K. O. Farrell. 1993 Modification of developmental instability and fitness: malathion-resistance in the Australian sheep blowfly, Lucilia cuprina. In T. A. Markow [ed.], Developmental instability, its origins and evolutionary implications, 6078. Kluwer Academic Publishers, Dordrecht, The Netherlands.
Moore, W. S. 1977 An evaluation of narrow hybrid zones in vertebrates. Quarterly Review of Biology 52: 263277.[CrossRef][ISI]
, and D. B. Buchanan. 1985 Stability of the Northern Flicker hybrid zone in historical times: implications for adaptive speciation theory. Evolution 39: 135151.
, and J. T. Price. 1993 Nature of selection in the northern flicker hybrid zone and its implications for speciation theory. In R. G. Harrison [ed.], Hybrid zones and the evolutionary process, 196225. Oxford University Press, Oxford.
Munter, R. C., and R. A. Grande. 1983 Plant tissue and soil extract analysis by ICP-atomic emission spectrometry. In R. M. Barnes [ed.], Developments in atomic plasma spectrochemical analysis: Proceedings of international winter conference, San Juan, Puerto Rico, 711 January, 1980, 653672. John Wiley & Sons, New York, NY.
Norusis, M. J. 1988 SPSS/PC+ V2.0. Base manual, advanced statistics. SPSS Inc., Chicago, IL.
Randerson, P. F. 1993 Ordination. In J. C. Fry [ed.], Biological data analysis: a practical approach, 173217. Oxford University Press, Oxford.
Rendig, V. V., and H. M. Taylor. 1989 Principles of soil-plant interrelationships. McGraw-Hill, New York, NY.
Rice, W. R. 1989 Analysis tables of statistical tests. Evolution 43: 223225.[CrossRef][ISI]
Salisbury, F. B., and C. W. Ross. 1992 Plant physiology (4th ed.). Wadsworth, Belmount, CA.
Shkolnik, M. Ya. 1984 Trace elements in plants. Elsevier, Amsterdam.
Ter Braak, C. J. F. 1992 CANOCOa FORTRAN program for canonical community ordination. Microcomputer Power, Ithaca, NY.
Tilman, D. 1985 The resource ratio hypothesis of succession. American Naturalist 125: 827852.[CrossRef][ISI]
Wang, H. 1996 Ecological analysis of a stable hybrid zone between two big sagebrush subspecies (Artemisia tridentata ssp. tridentata and A. t ssp. vaseyana. Ph.D. dissertation, Wayne State University, Detroit, MI.
, D. W. Byrd, J. L. Howard, E. D. McArthur, J. H. Graham, and D. C. Freeman. 1998 Narrow hybrid zone between two subspecies of big sagebrush (Artemisia tridentata: Asteraceae) VI. Soil properties. International Journal of Plant Science 159: 139147.[CrossRef]
, E. D. McArthur, S. C. Sanderson, J. H. Graham, and D. C. Freeman. 1997 Narrow hybrid zone between two subspecies of big sagebrush (Artemisia tridentata: Asteraceae). V. Reciprocal transplant experiments. Evolution 51: 95102.[CrossRef][ISI]
Young, N. D. 1996 An analysis of the causes of genetic isolation in two Pacific Coast iris hybrid zones. Canadian Journal of Botany 74: 20062013.[CrossRef]
Zar, J. H. 1984 Biostatistical analysis, 2nd ed. Prentice-Hall, Englewood Cliffs, NJ.
This article has been cited by other articles:
![]() |
K. J. Miglia, E. D. McArthur, R. S. Redman, R. J. Rodriguez, J. C. Zak, and D. C. Freeman Genotype, soil type, and locale effects on reciprocal transplant vigor, endophyte growth, and microbial functional diversity of a narrow sagebrush hybrid zone in Salt Creek Canyon, Utah Am. J. Botany, March 1, 2007; 94(3): 425 - 436. [Abstract] [Full Text] [PDF] |
||||
![]() |
F. Fritsche and O. Kaltz Is the Prunella (Lamiaceae) hybrid zone structured by an environmental gradient? Evidence from a reciprocaltransplant experiment Am. J. Botany, July 1, 2000; 87(7): 995 - 1003. [Abstract] [Full Text] |
||||
![]() |
E. D. McArthur and S. C. Sanderson Cytogeography and chromosome evolution of subgenus Tridentatae of Artemisia (Asteraceae) Am. J. Botany, December 1, 1999; 86(12): 1754 - 1775. [Abstract] [Full Text] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |