Am. J. Bot. Visit Plant Cell Online
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Submit a response
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via ISI Web of Science (17)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Tovar-Sánchez, E.
Right arrow Articles by Oyama, K.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Tovar-Sánchez, E.
Right arrow Articles by Oyama, K.
Agricola
Right arrow Articles by Tovar-Sánchez, E.
Right arrow Articles by Oyama, K.
(American Journal of Botany. 2004;91:1352-1363.)
© 2004 Botanical Society of America, Inc.


Population Biology

Natural hybridization and hybrid zones between Quercus crassifolia and Quercus crassipes (Fagaceae) in Mexico: morphological and molecular evidence1

Efraín Tovar-Sánchez and Ken Oyama2

Centro de Investigaciones en Ecosistemas, UNAM (Universidad Nacional Autónama de México), Campus Morelia, Antigua Carretera a Pátzcuaro No. 8701, Col. San José de la Huerta, C.P. 58190, Morelia, Michoacán, México

Received for publication July 10, 2003. Accepted for publication May 20, 2004.


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Hybrid zones provide interesting systems to study genetic and ecological interaction between different species. The correct identification of hybrids is necessary to understand the evolutionary process involved in hybridization. An oak species complex occurring in Mexico formed by two parental species, Quercus crassifolia H. & B. and Q. crassipes H. & B., and their putative hybrid species, Q. dysophylla, was analyzed with molecular markers (random amplified polymorphic DNA [RAPDs]) and morphological tools in seven hybrid zones (10 trees per taxa in each hybrid zone) and two pure sites for each parental species (20 trees per site). We tested whether geographic proximity of hybrid plants to the allopatric site of a parental species increases its morphological and genetic similarity with its parent. Seventeen morphological traits were measured in 8700 leaves from 290 trees. Total DNA of 250 individuals was analyzed with six diagnostic RAPD primers. Quercus crassifolia differed significantly from Q. crassipes in all the examined characters. Molecular markers and morphological characters were highly coincident and support the hypothesis of hybridization in this complex, although both species remain distinct in mixed stands. Clusters and a hybrid index (for molecular and morphological data) showed that individuals from the same parental species were more similar among themselves than to putative hybrids, indicating occasional hybridization with segregation in hybrid types or backcrossing to parents. Evidence does not indicate a unidirectional pattern of gene flow.

Key Words: Fagaceae • hybrid zones • hybridization • leaf morphology • Mexico • Quercus • RAPDs


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Hybridization is a natural phenomenon that occurs frequently in plants and animals (Harrison, 1993 ). This process produces new genetic combinations by the introduction of semi-compatible genes into another genotype, upon which interaction of environment and genetic variation can isolate a novel taxon from parental types. Hybrids may be defined as species, subspecies, variants, or races, depending on the degree of divergence (Futuyma, 1998 ). Interspecific gene transfer is an important evolutionary force, because the genetic material introduced by introgression exceeds that which is produced directly by mutation (Anderson, 1949 ).

The genus Quercus (Fagaceae) is one of the most diversified groups of temperate trees with more than 500 species distributed worldwide (Nixon, 1993 ). Hybridization and hybrid zones are common among oaks (Trelease, 1924 ; Palmer, 1948 ; Muller, 1952 ; Tucker, 1961 ; Cottam et al., 1982 ; Jensen et al., 1993 ; Spellenberg, 1995 ; Howard et al., 1997 ; Ishida et al., 2003 ). However, despite the perception that hybrid zones are well documented among oaks, few comparative analyses of oak hybrid zones have used both morphological characters and genetic markers (Howard et al., 1997 ).

Hybridization in oaks was initially detected based on morphological characters (Stebbins et al., 1947 ; Barlett, 1951 ; Tucker, 1961 ; Benson et al., 1967 ; Hardin, 1975 ; Cottam et al., 1982 ; Rushton, 1993 ). Leaf morphology in particular has been useful to demonstrate hybridization (Bacon and Spellenberg, 1996 ). However, in some cases morphological characters alone do not confirm unequivocally the existence of hybridization (Bacilieri et al., 1995 ; Manos et al., 1999 ; Mayol and Rosselló, 2001 ) requiring other methods such as DNA markers (Crawford et al., 1993 ; Rieseberg and Ellstrand, 1993 ). Random amplified polymorphic DNA (RAPD) markers have been particularly successful in the detection of interspecific hybridization and introgression in plants (Arnold et al., 1991 ; Arnold, 1993 ; Crawford et al., 1993 ; Cruzan and Arnold, 1993 ; Marsolais et al., 1993 ; Fritz et al., 1994 ; Smith et al., 1996 ; Samuel, 1999 ).

Mexico is considered one of the centers of diversification of the genus Quercus (oaks) (Muller and McVaugh, 1972 ; Rzedowski, 1978 ; González, 1993 ; Nixon, 1993 ), with 135– 150 species that include 86 endemics (Nixon, 1993 ). However, hybridization has only recently been reported for some species of Mexican oaks (McVaugh, 1974 ; Boecklen and Spellenberg, 1990 ; Spellenberg, 1992 , 1995 ; Spellenberg and Bacon, 1996 ; González-Rodríguez et al., 2004 , in press ). The taxonomy and evolutionary relationships of Mexican oaks are currently being studied, and some species complexes formed by potential parental species and putative hybrids have been recently detected by oak specialists. We chose two red oak species (subg. Erythrobalanus) for this study, Q. crassifolia H. & B. and Q. crassipes H. & B. These species have noticeable differences in several morphological characters (Romero, 1993 ) when they form allopatric pure stands. However, intermediate trees with atypical leaf shapes are observed when both species occur in sympatry, suggesting that hybridization may explain the observed variations. It is important to indicate that other oak species that can be considered as reasonable putative parents as judged from the morphological features of Q. dysophylla do not occur in the area. The taxonomic status of Q. dysophylla is still under discussion; some authors consider this species as a hybrid formed by Q. crassifolia x Q. crassipes (K. Nixon and S. Valencia, Cornell University and UNAM [Universidad Nacional Autónama de México], respectively, personal communication), while others recognize it as a different species (Romero, 1993 ; Zavala Chávez, 1995).

In this paper, we describe and compare the patterns of morphological and genetic variation of the Q. crassifolia x Q. crassipes complex, document the structure of overlapping zones and hybridization in the Eje Neovolcánico and Central Mexico, and assess the taxonomic distinctness of the two species and the putative hybrid.


    MATERIALS AND METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Study sites
Quercus crassifolia is distributed throughout the Sierra Madre Occidental (SMOc, northwest Mexico), whereas Q. crassipes is distributed along the Sierra Madre Oriental (SMOr) on the opposite side of the country. The distribution of these species overlaps along the Eje Neovolcánico, a range of volcanic mountains that traverses central Mexico in an east–west direction, where individuals with intermediate leaf morphology are present (Romero, 1993 ). Seven sympatric zones were chosen in the Eje Neovolcánico and Sierra Madre Oriental, where Q. crassifolia and Q. crassipes overlap geographically: Cantera, Canalejas, Tlaxco, Acajete, Esperanza, Agua Blanca, and Palo Bendito. These localities were chosen because they contain the highest numbers of intermediate individuals of this complex. Two localities were also chosen for each parental species, where they are dominant and no hybrids were observed (Fig. 1). In general, the hybrids occur in a very low frequency sporadically interspersed and near to the putative parents. Intermediate trees occur in more disturbed habitats but we cannot assess the preference of hybrids for some type of habitats (e.g., type of soils, forest gaps). It is rare to find intermediate trees outside of the contact zones. No evidence exists on another oak species hybridizing with trees of the species studied. When other oak species were found, they were white oaks that cannot hybridize with red oak species.



View larger version (20K):
[in this window]
[in a new window]
 
Fig. 1. Map of sampled populations of Quercus crassifolia x Q. crassipes complex. The mixed stands are represented by numbers, 1 = Cantera, 2 = Canalejas, 3 = Tlaxco, 4 = Acajete, 5 = Esperanza (located at the Eje Neovolcánico), 6 = Agua Blanca, 7 = Palo Bendito (located at the Sierra Madre Oriental)

 
Morphological data
Seventeen morphological characters of leaves were measured (Table 1) in 290 trees. A total of 210 trees were sampled in the seven mixed stands, 10 trees for each taxa: Q. crassifolia, Q. crassipes, and the hybrid. Twenty trees were sampled in each pure site (two sites for each parental species). Thirty mature leaves, without any apparent damage, were randomly sampled in each tree. Seventeen characters were measured in a total of 8700 leaves.


View this table:
[in this window]
[in a new window]
 
Table 1. List of the leaf morphological characters examined for the Q. crassifolia x Q. crassipes complex in Mexico

 
Molecular data
We collected undamaged young leaves from 250 trees of Q. crassifolia (n = 90), Q. crassipes (n = 90), and hybrids (n = 70), distributed in seven mixed stands and 40 trees in two pure zones for each parental species. Total DNA was extracted using a DNAeasy Plant Mini Kit (Qiagen, Valencia, California, USA). A genetic analysis was performed using RAPD (Welsh and McClelland, 1990 ; Williams et al., 1990 ), and the polymerase chain reaction (PCR). Forty-eight 10-base pair (bp) primers of random sequence (Kits A, B, C; Operon Technologies, Alameda, California, USA) were tested to find specific, diagnostic markers for each parental species. Diagnostic markers are those that are present in all individuals of one species and absent in a second species, while species-specific markers are those unique to one species but not necessarily present in all individuals within the species (Howard et al., 1997 ). Lastly, 250 individuals were tested using six primers (four species-specific markers [two for Q. crassifolia, OPA-09 and OPA-13, and two for Q. crassipes, OPB-01 and OPC-06], and two diagnostic markers [Q. crassifolia, OPA-14, and Q. crassipes, OPB-18]). DNA fragments were separated by electrophoresis on 1.5% agarose gel, stained with ethidium bromide, and developed on a UV table. The molecular mass of the DNA fragments was estimated with a 1-kilo base (kb) DNA ladder.

Statistical analysis
Nested variance analyses were conducted (ANOVA) to determine the effects of oak species, locality, and individual (tree) on the morphological leaf variability of each of the 17 studied characters (Table 1). Hybrids were not included in this analysis. Trees were considered as a random factor nested within species, because they were representative of each population. Percentage data were corrected as X = arcsin (%)1/2;, and discontinuous data were transformed as X = (x)1/2; + 0.5 (Zar, 1999 ).

To quantify variation in leaf morphology among oak species and hybrids, we randomly selected 30 leaves from vouchered specimens, and all the 17 characters were measured. All morphological characters were used for discriminant analysis to assess the most useful character for taxonomic discrimination between Q. crassifolia and Q. crassipes and to determine how leaf morphological characters separate individuals into groups. Seventeen characters were measured for the character count procedure to determine hybrids with intermediate leaf morphology of the Q. crassifolia x Q. crassipes complex following the procedure of Wilson (1992) . Variation in size and shape due to shape alone was quantified using the ratios of the sums of eigen values from the discriminant analysis (Darroch and Mosimann, 1985 ).

In the seven mixed stands, the Anderson hybrid index (Anderson, 1949 , 1953 ) was used to identify intermediate individuals and possible backcrosses. A histogram for each mixed stand was obtained. The Anderson hybrid index was calculated using the 17 morphological characters because they demonstrated differences between parental species. The representative characters of Q. crassipes received a rank of 2; Q. crassifolia characters were assigned a rank of 0; while intermediate characters were assigned a rank of 1 (Wilson, 1992 ).

A general cluster diagram for all zones was obtained, including pure and mixed stands. STATISTICA 6.0 for Windows was used for all the statistical analyses (Statsoft, 1998 ).

The maximum likelihood (ML) hybrid index score from RAPD analysis was calculated using Hardig-Hybrid software (Hardig et al., 2000 ). This index is useful to identify intermediate individuals, showing backcrosses as well as the structures of hybrid swarms. The results were represented in a frequency histogram for each mixed stand and for the pure sites. A Mantel Z-test matrix and Tools for Population Genetic Analyses (TFPGA, version 1.3) were used to test isolation by distance, as well as to create a general cluster diagram for the mixed stands and the four pure zones (Miller, 1997 ).


    RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Morphological analysis for parental and hybrid plants
All of the examined characters in Quercus crassifolia differ significantly from those of Q. crassipes. Locality and tree also had highly significant effects on each of the measured leaf characters (Table 2). Three characters (NA, LL/MWL, and LLW/MWL) were clearly separated without overlap between the two species in the pure zones as well as in the mixed stands and can be considered as species-specific characters. Hybrids presented intermediate morphological leaf traits between Q. crassifolia and Q. crassipes. The binomial sign analysis showed that the deviation was highly significant (P < 0.001), accepting the hybridization hypothesis. Intermediate characters were found in the 17 examined characters of the samples from Palo Bendito ; 16 of 17 examined characters in the samples from Cantera and Esperanza; and 15 of 17 examined characters in the samples from Canalejas, Tlaxco, Agua Blanca, and Acajete (Table 3).


View this table:
[in this window]
[in a new window]
 
Table 2. Mean ± standard error and nested ANOVA results (F statistics) for all characters of Quercus crassifolia and Q. crassipes in Mexico. F significant at P < 0.001 (***); P < 0.01 (**); and P < 0.05 (*)

 

View this table:
[in this window]
[in a new window]
 
Table 3. Character count procedure from the seven hybrid zones of the Quercus crassifolia x Q. crassipes complex in Mexico. Abbreviations of characters are described in Table 1 ; + = the hybrid is intermediate between the parental species, – = the hybrid is not intermediate, SD = standard deviation

 

View this table:
[in this window]
[in a new window]
 
Table 3. Continued

 
Discriminant function analysis for size and shape variation showed that in all of the analyzed hybrid zones the hybrids had intermediate characters between their parents. The discriminant axis described a gradient of leaf shapes from elliptic narrow (Q. crassipes) to the hybrids and the obovate leaves of Q. crassifolia (Fig. 2).



View larger version (34K):
[in this window]
[in a new window]
 
Fig. 2. Discriminant function analysis for leaf morphology variation in Quercus crassifolia x Q. crassipes complex (17 measured characters) in seven hybrids zones in Mexico. See abbreviations in Table 1

 
For the Q. crassifolia x Q. crassipes comparison, one highly significant (P < 0.001) discriminant axis accounted from 81.23% (Esperanza) to 94.78% (Acajete) of the variation between taxa. Discriminant function analysis for leaf morphology variation showed that in general the MWL contributed with the highest value to the ordination model for the first axis in Cantera, Canalejas, and Agua Blanca while in Tlaxco and Palo Bendito it was LL, and in Esperanza HW% and in Acajete it was LLW. On the other hand, the character with the highest value for the second axis was LL in Cantera, Canalejas, Agua Blanca, Esperanza, and Acajete, LL/MWL in Palo Bendito, and TLL in Tlaxco (Table 4 ).


View this table:
[in this window]
[in a new window]
 
Table 4. Discriminant function analysis for leaf morphology variation from seven hybrid zones of the Quercus crassifolia x Q. crassipes complex in Mexico. Abbreviations of characters are described in Table 1

 

View this table:
[in this window]
[in a new window]
 
Table 4. Continued

 
The shape variable produced two significant discriminant axes, and the ratio of the sums of eigenvalues of the Cantera hybrid zone indicated that approximately 4.48/7.16 = 62.57% of shape alone was accountable for the total variation in size and shape. In the other hybrid zones, less than 50% of the total variation in size and shape was attributable to the shape alone: in Agua Blanca, 1.65/6.44 = 25.62 %; Canalejas, 2.48/ 8.14 = 30.47%; Palo Bendito, 2.88/8.52 = 33.80%; Tlaxco, 3.55/9.48 = 37.45%; Acajete, 2.28/5.17 = 44.10%; and Esperanza, 3.34/7.94 = 42.06%.

In a global cluster analysis based on morphological leaf traits for the seven mixed stands and four pure populations (two for Q. crassifolia and two for Q. crassipes), the populations having the same parental species were more similar to each other. For the five hybrid zones located on the Eje Neovolcánico (mixed stands from 1 to 5), nearness of the hybrid zones to the allopatric putative parent correlated with increasing similarity of the complex to the parental species. Lastly, mixed stands 6 and 7, located on the Sierra Madre Oriental, were more similar to Q. crassifolia (Fig. 3, left).



View larger version (8K):
[in this window]
[in a new window]
 
Fig. 3. Phenogram of similarity based on leaf morphology and RAPD data in seven hybrids zones of the Quercus crassifolia x Q. crassipes complex in Mexico. The mixed stands are represented by numbers, 1 = Cantera, 2 = Canalejas, 3 = Tlaxco, 4 = Acajete, 5 = Esperanza (located at the Eje Neovolcánico), 6 = Agua Blanca, 7 = Palo Bendito (located at the Sierra Madre Oriental)

 
We used 17 morphological characters to estimate the status of Q. crassifolia and Q. crassipes plants marked previously as "pure" or "hybrid" plants. The results for the Anderson hybrid index support the field identification of 250 plants (Fig. 4). In general, the frequency histogram (Anderson hybrid index) had a pattern similar to the global cluster analysis. In Eje Neovolcánico, the mixed stands (Cantera and Canalejas) nearer to SMOc presented unidirectional introgression towards Q. crassifolia. Subsequently, the site closest to central Mexico (Tlaxco) was the only site that showed bidirectional introgression, but the two mixed stands nearest to SMOr (Acajete and Esperanza) registered unidirectional introgression towards Q. crassipes. Finally, the two hybrid zones localized north of Tlaxco (Agua Blanca and Palo Bendito) on SMOr presented unidirectional introgression towards Q. crassifolia (Fig. 4).



View larger version (45K):
[in this window]
[in a new window]
 
Fig. 4. Frequency distribution of individuals vs. the Anderson hybrid index derived from 17 morphological characters in two pure and seven hybrid zones. The number of each plant evaluated is represented. The figure shows the cluster results of the pure or nearly pure parentals in four zones (index = 0 or 34), a cluster for probable F1 hybrids (index = 12.75–25.5), probable backcrosses towards Q. crassifolia (index = 8.5–10.63), and finally probable backcrosses toward Q. crassipes (index = 23.38–25.5). The plants misidentified in the field are indicated in boldface type

 
Individuals A10 (Cantera), A7, A8 (Tlaxco), and A9 (Agua Blanca) were classified as Q. crassifolia in the field, but the index showed that they were backcrosses towards Q. crassifolia. Individuals B4 and B8 (Tlaxco), B2 (Acajete), and B2 (Esperanza) were originally marked as Q. crassipes, but the index analysis confirmed that they were backcrosses toward Q. crassipes. Individuals AB1 and AB9 (Agua Blanca), as well as AB2 (Acajete), and AB2 and AB10 (Esperanza), were marked as hybrids, but the index showed that they were individuals of Q. crassifolia and Q. crassipes, respectively. Also, the individuals AB8, AB9, and AB10 (Cantera); AB1 and AB10 (Canalejas); AB8 and AB10 (Tlaxco); AB6 and AB10 (Agua Blanca); and AB3 and AB5 (Palo Bendito) were classified as hybrid plants, but the index analysis indicated that they were backcrosses toward Q. crassifolia. Finally, AB3 (Tlaxco), AB3, AB9, and AB10 (Acajete), and AB3, AB5, and AB9 (Esperanza) were marked as hybrids, but the data showed that they were backcrosses toward Q. crassipes (Fig. 4).

Genetic analysis (RAPDs) of parental and hybrid plants
We used six primers to estimate the genetic status of Q. crassifolia and Q. crassipes plants morphologically identified as "pure" or "hybrid" plants. These primers yielded 49 distinct markers (bands). The RAPD analysis proved to be a powerful tool for characterizing hybrid individuals between Q. crassifolia and Q. crassipes. The analysis of the ML hybrid index using six RAPD markers supported the field identification of 250 plants (Fig. 5). Individuals A10 (Cantera); A1, A7, and A8 (Tlaxco); and A5 and A6 (Palo Bendito) were classified as Q. crassifolia in the field, but the RAPD showed that they were backcrosses toward Q. crassifolia. Individuals B4 and B8 (Tlaxco) and B2 (Esperanza) were originally marked as Q. crassipes, but RAPDs analysis confirmed that they were backcrosses toward Q. crassipes. The individuals AB9 (Agua Blanca) and AB7 (Tlaxco) were marked as hybrids, but RAPDs and subsequent field examinations showed that they were individuals of Q. crassifolia and Q. crassipes, respectively. Individuals AB10 (Cantera), AB1 (Agua Blanca), and AB3, AB5, and AB7 (Palo Bendito) were classified as hybrid plants, but the molecular analysis indicated that they were backcrosses toward Q. crassifolia. Lastly, AB3, AB6, and AB9 (Tlaxco) ; AB4, AB5, AB8, AB9, and AB10 (Acajete); and AB2 and AB10 (Esperanza) were identified as hybrids, but RAPD data showed that they were backcrosses toward Q. crassipes (Fig. 5).



View larger version (42K):
[in this window]
[in a new window]
 
Fig. 5. Frequency distribution of individuals vs. the Hardig hybrid index derived from the RAPD band data using six primers. The number of each plant evaluated is represented. The figure shows cluster results of the pure or nearly pure parentals in four zones (index = 0 or 1), a cluster of probable F1 hybrids (index = 0.437–0.562), probable backcrosses toward Q. crassifolia (index = 0.687–0.750), and finally probable backcrosses towards Q. crassipes (index = 0.250–0.312). The plants misidentified in the field are indicated in boldface type

 
Few plants had perfect marker additivity as may be expected in F1 (12 plants), but 54 individuals were interpreted as F1, 26 deviated by only one character (0.437–0.562), and 16 plants deviated by only two characters (0.375–0.625). Twelve plants were interpreted as backcrosses toward Q. crassipes (0.250– 0.312) and 11 as backcrosses toward Q. crassifolia (0.687– 0.750) (Fig. 5).

The mean ML hybrid index for Q. crassifolia individuals was 0.85 (SD 0.07), for Q. crassipes was 0.08 (SD 0.06), and for the hybrids, 0.48 (SD 0.09).

The Hardig hybrid index showed the same general pattern as the Anderson hybrid index (see earlier). No introgression was registered only in Canalejas (Fig. 5).

The cluster analysis for genetic (RAPDs) and morphological data showed that these are very similar. The allopatric zones (two zones for Q. crassifolia and two for Q. crassipes) were located at the edges, while the seven hybrid zones were located between them (Fig. 3).

The Mantel Z-test matrix showed that no correlation exists between the geographic distances and the genetic distances for any of the species (Q. crassifolia, r = 0.37, P > 0.05 ; and Q. crassipes, r = 0.14, P > 0.05).


    DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Oaks frequently present complex patterns of variation leading to taxonomic problems in differentiating species (Burger, 1975 ). Interspecific hybridization and shared ancestral polymorphisms are two of the most common explanations for the observed pattern of variation (Jensen et al., 1993 ; Kleinschmit et al., 1995 ; Bruschi et al., 2000 ). Quercus crassifolia and Q. crassipes form hybrids, but they remain morphologically distinct in their allopatric and sympatric distributions. All the examined morphological leaf characters in these species differed significantly by localities, populations, and individuals. Relatively few diagnostic molecular markers differentiated between Q. crassifolia and Q. crassipes. However, these markers demonstrated geographic consistency in support of the morphological evidence, indicating that each species is distinct and that each has some degree of genetic cohesiveness.

Furthermore, the character count procedure has allowed us to confirm statistically that Q. dysophylla is the result of hybridization between Q. crassifolia and Q. crassipes in the seven hybrid zones. The ordination analysis also demonstrated that hybrids presented intermediate morphology between the parental species and that leaf shape explains a major percentage of variation.

Quercus crassifolia ranges from the Sierra Madre Occidental (SMOc) to the center of Mexico along the Eje Neovolcánico, whereas Q. crassipes ranges from the Sierra Madre Oriental (SMOr) to the Eje Neovolcánico, where both species overlap producing hybrid zones. The Eje Neovolcánico, an orographic system that traverses the central part of the country in an east–west direction, is considered geologically the youngest mountain range in Mexico and contains valleys higher than 2000 m in altitude and the tallest mountains in Mexico (Ferrusquia-Villafranca, 1993 ). Phylogeographic studies are in progress to understand the process of migration and the colonization routes of this oak complex, as has been done for other oak species (e.g., Dumolin-Lapègue et al., 1997 ).

The genetic results indicate that the introgression process is present in both species, but the direction changes depending on the localization of the hybrid zone. The hybrid zones closest to the SMOc (Cantera and Canalejas) registered unidirectional introgression towards Q. crassifolia, while the hybrids from the intermediate locality between the two mountain ranges (Tlaxco) showed bidirectional introgression, and the hybrids from the two closest localities to the SMOr registered unidirectional introgression towards Q. crassipes. These findings suggest that the closeness of hybrids to an allopatric site of either parental species is directly related to their similarity and vice versa. Thus, the Eje Neovolcánico acts as a corridor where the proximity to an allopatric site favors the introgression of the hybrid towards the parental species, increasing its variation from the species with which it is maintaining a genetic exchange, diluting the limits with the parental species in the allopatric site. Lastly, the two hybrid zones located north of Tlaxco (SMOr) showed unidirectional introgression towards Q. crassifolia. These results confirm that patterns of variation in oaks do not follow simple monotonic clines (e.g., Barton and Hewit, 1985 ) but form complex mosaic zones characterized by patches of pure populations and mixed populations scattered across a zone of overlap (Howard et al., 1997 ). A bidirectional hybrid zone was detected for Q. crassifolia and Q. crassipes.

In the seven hybrid zones studied along the Eje Neovolcánico, where intermediate plants are mixed with their parental species, hybrids are rare and they are in a narrow contact zone between well-differentiated taxa. The presence of hybrid individuals in the hybrid zones was very low (between 10 and 17 trees), requiring an extensive field search. Oak hybrids are produced in an isolated and sporadic manner and they may introgress with parental species (Bacon and Spellenberg, 1996 ). Hybrid zones with high levels of disturbance (i.e., Canalejas, Acajete, and Esperanza) were the ones with the highest number of hybrid individuals (mostly juveniles). Disturbances produced by human activities such as logging, deforestation, fires, and agriculture, may enhance the establishment of hybrids as they modify reproductive barriers (Arnold et al., 1990 ; Klier et al., 1991 ).

Our results suggest that the sympatric zones of Q. crassifolia and Q. crassipes are mosaic hybrid zones as proposed by Howard (1982 , 1986 ) and Harrison (1986 , 1990 ), because of the patchy distribution pattern of the parental species in sympatric and allopatric sites and the lack of a gradual transition from Q. crassipes to Q. crassifolia. It is important to mention that hybrid plants were less frequent than putative parents in the mixed stands and that Q. crassipes prefers drier habitats and lower sites than Q. crassifolia. Ecological divergence rather than genetic incompatibility may maintain hybrid zones (Jiggins and Mallet, 2000 ) by causing local adaptations to different environmental conditions (e.g., Howard et al., 1997 ).

In summary, we found that molecular markers (RAPD) and morphological leaf traits are highly coincident and support the phenomenon of hybridization between Q. crassifolia and Q. crassipes complex (Fig. 3). Inasmuch as hybridization was evident, both species remain distinct in mixed stands. We also observed that the Eje Neovolcánico acts as a corridor where proximity to an allopatric site favors the introgression of the hybrid with its parental species, increasing its divergence from the species with which it maintains a genetic exchange, and thus diluting the limits with parental species in the allopatric site. Hybrid plants constitute a heterogeneous group in which many individuals were F1 and others appeared as backcrosses of Q. crassifolia or Q. crassipes, depending on the locality.

Our data and field observations suggest that the sympatric zones of Q. crassifolia and Q. crassipes must be considered as mosaic hybrid zones (e.g., Howard, 1982 , 1986 ; Harrison, 1986 , 1990 ), because of the patchy distribution pattern of the parental species in sympatric and allopatric sites, and there is not a gradual transition from Q. crassipes to Q. crassifolia. Finally, we suggest that Q. dysophylla does not deserve the status of species but it must be recognized as an entity of potential evolutionary importance, named as Quercus x dysophylla Benth. pro sp.


    FOOTNOTES
 
1 The authors thank Richard Spellenberg and one anonymous reviewer whose comments improved previous versions of this manuscript. We also thank Antonio González, Susana Valencia, Patricia Mussali-Galante, Rocío Esteban-Jiménez, Marco Romero, and Nidia Pérez- Nasser for technical assistance. This research was supported by grants from CONACYT- Mexico (38550-V) and DGAPA, UNAM (IN212999) to K. O., a PAEP-UNAM grant and a CONACYT scholarship to E. T. S. Back

2 akoyama{at}oikos.unam.mx Back


    LITERATURE CITED
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Anderson E. 1949 Introgressive hybridization. John Wiley, New York, New York, USA

Anderson E. 1953 Introgressive hybridization. Biological Reviews 28: 280-307[CrossRef]

Arnold M. L. 1993 Iris nelsonii (Iridaceae): origin and genetic composition of a homoploid hybrid species. American Journal of Botany 80: 577-591[CrossRef][ISI]

Arnold M. L. C. M. Buckner J. J. Robinson 1991 Pollen-mediated introgression and hybrid speciation in Louisiana irises. Proceedings of the National Academy of Sciences, USA 80: 1398-1402

Arnold M. L. J. L. Hamrick B. D. Bennett 1990 Allozyme variation in Louisiana irises: a test for introgression and hybrid speciation. Heredity 65: 297-306[ISI]

Bacilieri R. A. Ducousso A. Kremer 1995 Genetical, morphological, ecological and phonological differentiation between Quercus petrea (Matt.) Liebl. and Quercus robur L. in a mixed stand of northwest France. Silvae Genetica 44: 1-10

Bacon J. R. R. Spellenberg 1996 Hybridization in two distantly related Mexican black oaks Quercus conzattii and Quercus eduardii (Fagaceae: Quercus: section Lobatae). SIDA 17: 17-41

Barlett H. H. 1951 Regression of Quercus deamii toward Quercus macrocarpa and Quercus muhlenbergii. Rhodora 53: 249-264

Barton N. H. G. M. Hewit 1985 Analysis of hybrid zones. Annual Review of Ecology and Systematics 16: 113-148[CrossRef][ISI]

Benson L. E. A. Phillips P. A. Wilder 1967 Evolutionary sorting of characters in a hybrid swarm. I. Direction of slope. American Journal of Botany 54: 1017-1026[CrossRef][ISI]

Boecklen W. R. Spellenberg 1990 Structure of herbivore communities in two oak (Quercus spp.) hybrid zones. Oecologia 85: 92-100[CrossRef][ISI]

Bruschi P. G. G. Vendramin F. Bussotti P. Grossoni 2000 Morphological and molecular differentiation between Quercus petrea (Matt.) Liebl. and Quercus pubescens Willd. (Fagaceae) in northern and central Italy. Annals of Botany 85: 325-333[Abstract/Free Full Text]

Burger W. C. 1975 The species concept in Quercus. Taxon 24: 45-50[CrossRef]

Cottam W. P. J. M. Tucker F. S. Santamour 1982 Oak hybridization at the University of Utah. State Arboretum of Utah, Salt Lake City, Utah, USA

Crawford D. J. S. Brauner M. B. Crosner T. F. Stuessy 1993 Use of RAPD markers to document the origin of the intergeneric hybrid x Margyracaena skottsbergii (Rosaceae) on the Juan Fernandez Island. American Journal of Botany 80: 89-92

Cruzan M. B. M. L. Arnold 1993 Ecological and genetic association in an Iris hybrid zone. Evolution 47: 1432-1445[CrossRef][ISI]

Darroch J. N. J. E. Mosimann 1985 Canonical and principal components of shape. Biometrika 72: 241-252[Abstract/Free Full Text]

Dumolin-Lapègue S. B. Demesure S. Fineschi V. Le Corre R. J. Petit 1997 Phylogeographic structure of white oaks throughout the European Continent. Genetics 146: 1475-1487[Abstract]

Ferrusquía-Villafranca I. 1993 Geology of Mexico: a synopsis. In T. P. Ramamoorthy, R. Bye, A. Lot, and J. Fa [eds.], Biological diversity of Mexico: origins and distribution, 3–107. Oxford University Press, New York, New York, USA

Fritz R. S. C. M. Nichols-Orians S. J. Brunsfeld 1994 Interspecific hybridization of plants and resistance to herbivores: hypotheses, genetics, and variable responses in a diverse herbivore community. Oecologia 97: 106-117[CrossRef][ISI]

Futuyma D. J. 1998 Evolutionary biology. Sinauer, Sunderland, Massachusetts, USA

González R. R. 1993 La diversidad de los encinos mexicanos. Revista de la Sociedad de Historia Natural 44: 125-142

González-Rodríguez A. D. M. Arias S. Valencia K. Oyama 2004 Morphological and RAPD analysis of hybridization between Quercus affinis and Quercus laurina (Fagaceae), two Mexican red oaks. American Journal of Botany 91: 401-409[Abstract/Free Full Text]

González-Rodríguez A. J. F. Bain J. L. Golden K. Oyama In press Chloroplast DNA variation in the Quercus affinis-Q. laurina complex in Mexico: geographic structure and associations with nuclear and morphological variation. Molecular Ecology in press

Hardig T. M. S. J. Brunsfeld R. S. Fritz M. Morgan C. M. Orians 2000 Morphological and molecular evidence for hybridization and introgression in a willow (Salix) hybrid zone. Molecular Ecology 9: 9-24[CrossRef][Medline]

Hardin J. W. 1975 Hybridization and introgression in Quercus alba. Journal of the Arnold Arboretum 56: 336-363[ISI]

Harrison R. G. 1986 Pattern and process in a narrow hybrid zone. Heredity 56: 337-149[ISI]

Harrison R. G. 1990 Hybrids zones: windows on evolutionary process. Oxford Surveys in Evolutionary Biology 7: 69-128

Harrison R. G. 1993 Hybrids and hybrid zones: historical perspective. In R. G. Harrison [ed.], Hybrid zones and the evolutionary process, 3–12. Oxford University Press, Oxford, UK

Howard D. J. 1982 Speciation and coexistence in a group of closely related ground crickets. Yale University Press, New Haven, Connecticut, USA

Howard D. J. 1986 A zone of overlap and hybridization between two ground cricket species. Evolution 40: 34-43[CrossRef][ISI]

Howard D. J. R. W. Preszler J. Williams S. Fenchel W. J. Boecklen 1997 How discrete are oaks species? Insights from a hybrid zone between Quercus grisea and Quercus gambelii. Evolution 51: 747-755[CrossRef][ISI]

Ishida T. A. K. Hattori H. Sato M. T. Kimura 2003 Differentiation and hybridization between Quercus crispula and Q. dentata (Fagaceae): insights from morphological traits, amplified fragment length polymorphism markers, and leaf miner composition. American Journal of Botany 90: 769-776[Abstract/Free Full Text]

Jensen R. J. S. C. Hokanson J. G. Isebrands J. F. Hancock 1993 Morphometric variation in oaks of the Apostle Islands in Wisconsin: evidence of hybridization between Quercus rubra and Q. ellipsoidalis (Fagaceae). American Journal of Botany 80: 1358-1366[CrossRef][ISI]

Jiggins C. D. J. Mallet 2000 Bimodal hybrid zones and speciation. Trends in Ecology and Evolution 15: 250-255

Kleinschmit J. R. G. R. Bacilieri A. Kremer A. Roloff 1995 Comparison of morphological and genetic traits of penduculate oak (Quercus robur (L.)) and sessile oak (Quercus petrea (Matt.)) Liebl). Silvae Genetica 44: 5-6

Klier K. M. J. Leoschke J. F. Wendel 1991 Hybridization and introgression in white and yellow ladyslipper orchids (Cypripedium candidum and C. pubescens). Journal of Heredity 82: 305-319[Abstract/Free Full Text]

Mayol M. J. A. Roselló 2001 Why nuclear ribosomal DNA spacers (ITS) tell different stories in Quercus. Molecular Phylogenetic and Evolution 19: 167-176[CrossRef]

Manos P. S. J. J. Doyle K. C. Nixon 1999 Phylogeny, biogeography and process of molecular differentiation in Quercus subgenus Quercus (Fagaceae). Molecular Phylogenetic and Evolution 12: 333-349[CrossRef]

Marsolais L. V. J. S. Pringle B. N. White 1993 Assessment of random amplified polymorphic DNA (RAPD) as genetic markers for determining the origin of interspecific lilac hybrids. Taxon 42: 531-537[CrossRef][ISI]

McVaugh R. 1974 Fagaceae. Flora Novo-Galiciana. Contributions from the University of Michigan Herbarium 12: 1-93

Miller M. P. 1997 Tools for population genetic analyses (TFPGA) 1.3: a windows program for the analysis of allozyme and molecular population genetic data. Computer software distributed by the author, Northern Arizona University, Flagstaff, Arizona, USA

Muller C. 1952 Ecological control of hybridization in Quercus: a factor in the mechanism of evolution. Evolution 6: 147-161[CrossRef][ISI]

Muller C. R. McVaugh 1972 The oaks (Quercus), with comments on related species. Contributions from the University of Michigan Herbarium 9: 507-522

Nixon K. C. 1993 The genus Quercus in Mexico. In K. C. Nixon [ed.], Biological diversity of Mexico: origins and distributions, 447–458. Oxford University Press, New York, New York, USA

Palmer E. J. 1948 Hybrid oaks of North America. Journal of the Arnold Arboretum 29: 1-48

Rieseberg L. N. C. Ellstrand 1993 What can molecular and morphological markers tell us about plant hybridization?. Critical Reviews in Plant Science 12: 213-241[CrossRef]

Romero R. S. 1993 El género Quercus (Fagaceae) en el Estado de México. M.Sc. thesis, Universidad Nacional Autónoma de México, Distrito Federal, México

Rushton B. S. 1993 Natural hybridization within the genus QuercusL. Annals of Science Forestry Supplement 50: 73-90

Rzedowski J. 1978 Vegetación de México. Limusa, Distrito Federal, México

Samuel R. 1999 Identification of hybrids between Quercus petrea and Q. robur (Fagaceae): results obtained with RAPD markers confirm allozyme studies based on the Got–2 locus. Plant Systematics and Evolution 217: 137-146[CrossRef][ISI]

Smith J. F. C. C. Burke W. L. Wagner 1996 Interspecific hybridization in natural populations of Cyrtandra (Gesneriaceae) on the Hawaiian Islands: evidence from RAPD markers. Plant Systematics and Evolution 200: 61-77[CrossRef][ISI]

Spellenberg R. 1992 A new species of black oak (Quercus, subg. Erythrobalanus, Fagaceae) from the Sierra Madre Occidental, México. American Journal of Botany 79: 1200-1206[CrossRef][ISI]

Spellenberg R. 1995 On the hybrid nature of Quercus basaseachicensis (Fagaceae: Sect. Quercus). Sida 16: 427-437

Spellenberg R. J. R. Bacon 1996 Taxonomy and distribution of a natural group of black oaks of Mexico (Quercus, section Lobatae, subsection Racemiflorae). Systematic Botany 21: 85-99

Statsoft. 1998 STATISTICA for Windows. Manual version 6.0. Statsoft, Tulsa, Oklahoma, USA

Stebbins G. L. E. B. Matzke C. Epling 1947 Hybridization in a population of Quercus marilandica and Q. ilifolia. Evolution 1: 79-88[CrossRef][ISI]

Trelease W. 1924 The American oaks. Memories of the National Academy of Science 20: 1-255

Tucker J. M. 1961 Studies in the Quercus undulate complex I. A preliminary statement. American Journal of Botany 48: 202-208[CrossRef][ISI]

Welsh J. M. McClelland 1990 Fingerprinting genomes using PCR with arbitrary primers. Nucleic Acids Research 18: 7213-7218[Abstract/Free Full Text]

Williams J. G. K. A. R. Kubelik K. L. Livak J. A. Rafalski S. V. Tingey 1990 DNA polymorphisms amplified by arbitrary primers are useful as genetics markers. Nucleic Acids Research 18: 6531-6535[Abstract/Free Full Text]

Wilson P. 1992 On inferring hybridity from morphological intermediacy. Taxon 41: 11-23

Zar J. H. 1999 Biostatistical analysis, 4th ed. Prentice Hall, Englewood Cliffs, New Jersey, USA

Zavala Chávez F. 1995 Encinos hidalguenses. Universidad Autónama de Chapingo, Chapingo, México




This article has been cited by other articles:


Home page
ANN BOT (LOND)Home page
J. Lihova, J. Kucera, M. Perny, and K. Marhold
Hybridization between Two Polyploid Cardamine (Brassicaceae) Species in North-western Spain: Discordance Between Morphological and Genetic Variation Patterns
Ann. Bot., June 1, 2007; 99(6): 1083 - 1096.
[Abstract] [Full Text] [PDF]


Home page
ANN BOT (LOND)Home page
F. Gugerli, J.-C. Walser, K. Dounavi, R. Holderegger, and R. Finkeldey
Coincidence of Small-scale Spatial Discontinuities in Leaf Morphology and Nuclear Microsatellite Variation of Quercus petraea and Q. robur in a Mixed Forest
Ann. Bot., April 1, 2007; 99(4): 713 - 722.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Submit a response
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via ISI Web of Science (17)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Tovar-Sánchez, E.
Right arrow Articles by Oyama, K.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Tovar-Sánchez, E.
Right arrow Articles by Oyama, K.
Agricola
Right arrow Articles by Tovar-Sánchez, E.
Right arrow Articles by Oyama, K.


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS