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


     


  Add to CiteULike   Add to Complore   Add to Connotea   Add to Del.icio.us   Add to Digg   Add to Facebook   Add to Reddit   Add to Technorati   Add to Twitter
What's this?
This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Supplementary Data
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 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 Web of Science (2)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Zelaya, I. A.
Right arrow Articles by VanGessel, M. J.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Zelaya, I. A.
Right arrow Articles by VanGessel, M. J.
Agricola
Right arrow Articles by Zelaya, I. A.
Right arrow Articles by VanGessel, M. J.
Social Bookmarking
 Add to CiteULike   Add to Complore   Add to Connotea   Add to Del.icio.us   Add to Digg   Add to Facebook   Add to Reddit   Add to Technorati   Add to Twitter  
What's this?
(American Journal of Botany. 2007;94:660-673.)
© 2007 Botanical Society of America, Inc.


Reproductive Biology

Transfer of glyphosate resistance: evidence of hybridization in Conyza (Asteraceae)1

Ian A. Zelaya4, Micheal D. K. Owen and Mark J. VanGessel

2Iowa State University, Department of Agronomy, 2104 Agronomy Hall, Ames, Iowa 50011-1011 USA; 3University of Delaware, Department of Plant and Soil Sciences, Research and Education Center, 16684 County Seat Highway, Georgetown, Delaware 19947-9575 USA

Received for publication May 14, 2006. Accepted for publication February 11, 2007.

ABSTRACT

Transfer of herbicide resistance genes between crops and weeds is relatively well documented; however, far less information exists for weed-to-weed interactions. The hybridization between the weedy diploids Conyza canadensis (2n = 18) and C. ramosissima (2n = 18) was investigated by monitoring transmission of the allele conferring resistance to N-phosphonomethyl glycine (glyphosate). In a multivariate quantitative trait analysis, we described the phylogenic relationship of the plants, whereas we tested seed viability to assess potential postzygotic reproductive barriers (PZRB) thus affecting the potential establishment of hybrid populations in the wild. When inflorescences were allowed to interact freely, approximately 3% of C. ramosissima or C. canadensis ova were fertilized by pollen of the opposing species and produced viable seeds; >95% of the ova were fertilized under no-pollen competition conditions (emasculation). The interspecific Conyza hybrid (

Formula

) demonstrated an intermediate phenotype between the parents but superior resistance to glyphosate compared to the resistant C. canadensis parent. Inheritance of glyphosate resistance in the selfed

Formula

(

Formula

) followed the partially dominant nuclear, single-gene model;

Formula

backcrosses confirmed successful introgression of the resistance allele to either parent. Negligible PZRB were observed in the hybrid progenies, confirming fertility of the C. canadensis x C. ramosissima nothotaxa. The implications of introgressive hybridization for herbicide resistance management and taxonomy of Conyza are discussed.

Key Words: allogamy • Conyza canadensisConyza ramosissima • gene flow • herbicide resistance • interspecific hybridization • shikimic acid • transgressive segregation

Interspecific hybridization refers to the cross-fertilization between two species that produces a fertile or infertile progeny with phenotypic traits of both parents; this process of interspecific gene transfer promotes genetic diversity and genome evolution (Abbott, 1992 ; Barton, 2001 ). This natural process has been utilized in breeding efforts to improve crop traits (Anamthawat–Jónsson, 2001 ). This process, however, may also facilitate the rapid evolution and adaptation of introduced plant pathogens and contribute to the genetic diversity of crop pests, and thus it may increase production difficulties in current agroecosystems (Teal and Oostendorp, 1995 ; Schardl and Craven, 2003 ). Furthermore, interspecific hybridization may adversely affect crop production and weed management, as interspecific transfer of herbicide resistance and genetically engineered genes has been documented (Owen and Zelaya, 2005 ).

Much information exists regarding the transgene flow and transfer of herbicide resistance genes between crops and their wild relatives (Kwon and Kim, 2001 ; Ellstrand, 2003 ; Légère, 2005 ; Guadagnuolo et al., 2006 ; Reichman et al., 2006 ). However, far less attention has been focused on gene flow between weed species and the impact on dissemination of herbicide resistance alleles or the evolution of novel taxa with diverse "weedy" traits. Current estimates of weed-to-weed herbicide resistance transfer based on four genera vary from 0.15% to 85% as predicted by the frequency of resistant individuals in the interspecific hybrid progeny (Table 1). To address this disparity in knowledge, we investigated hybridization between the weedy diploids Conyza canadensis (L.) Cronq. (2n = 18) and C. ramosissima Cronq. (2n = 18) by monitoring transmission of the allele conferring resistance to the herbicide N-phosphonomethyl glycine (glyphosate).


View this table:
[in this window]
[in a new window]

 
Table 1. Cases of confirmed transfer of herbicide resistance gene(s) through introgressive hybridization in weeds.

 
Asteraceae is one of the largest and most diverse family within dicotyledonous plants, members of which are present in every global environment except for aquatic habitats (Cronquist, 1980 ). The genus Conyza Less., represented by approximately 60 species, is composed of annual herbaceous plants that prosper chiefly in the tropical and subtropical regions of the globe (Nesom, 1990 ). Species endemic to the United States include the winter or summer annual forbs C. ramosissima and C. canadensis (Cronquist, 1980 ).

Conzya ramosissima was originally described in Illinois as Erigeron divaricatus (Michaux, 1803 ); the species was later annexed to the genus Conyza based on the eligulate character of its multiseriate pistillate florets (Cronquist, 1980 ). Both C. ramosissima and C. canadensis are ubiquitous in disturbed habitats and play important roles in primary ecological successions; however, only C. canadensis is reported to reduce yields in row crops, serve as an alternate host for diverse pests, and limit grazing by reducing the palatability of pasture forages (Steyermark, 1963 ; Cronquist, 1980 ). Importantly, C. canadensis has evolved resistance to amide, bipyridilium, glycine, imidazolinone, sulfonylurea, and triazine herbicides in more than 10 countries worldwide and thus is considered one of the 10 most important herbicide-resistant weeds (Heap, 2006 ). Conzya ramosissima is found from North Dakota to Pennsylvania and south from New Mexico to Alabama in the United States, in southern Canada, and northern Mexico, while the more cosmopolitan C. canadensis is distributed throughout the Americas, West Indies, Europe, and Africa (Steyermark, 1963 ). Interestingly, the Conyza taxon represents the most successful case of intercontinental colonization of the Americas to the Old World, to the extent that C. canadensis and C. floribunda H.B.K. are probably the most widely distributed species throughout the world (Thébaud and Abbott, 1995 ; Pruski and Sancho, 2006 ).

Glyphosate is one of the world's most important herbicides; glyphosate-resistant crops provide farmers with a simple, economical, and effective tool to manage a diverse weed flora, hence favoring the rapid adoption of this technology in many United States crop production systems (Owen and Zelaya, 2005 ). The mechanism of glyphosate action is the competitive inhibition with respect to the phosphate moiety of phosphoenolpyruvate in the reaction mediated by 3-phosphoshikimate 1-carboxyvinyltransferase (EPSPS; EC 2.5.1.19) (Steinrücken and Amrhein, 1980 ; Holländer-Czytko and Amrhein, 1983 ). The unique mode of action and limited metabolism in plants are purported reasons for the low frequency of evolved glyphosate resistance compared to other herbicide chemistries (Jasieniuk, 1985 ; Bradshaw et al., 1997 ). Since the commercial introduction of glyphosate-resistant crops in 1996 and the accompanying ubiquitous glyphosate use in these production systems, 12 weed species resistant to glyphosate have been identified worldwide, including confirmation in 15 independent C. canadensis populations within the United States and populations in Brazil and China (Heap, 2006 ).

Transmission of glyphosate resistance via gene flow has been documented between several crop and weed systems (Watrud et al., 2004 ; Légère, 2005 ; Guadagnuolo et al., 2006 ; Reichman et al., 2006 ). Presently, however, limited information exists regarding the level of within-species glyphosate resistance transfer in self-incompatible, wind-pollinated weed species as Plantago lanceolata, Amaranthus palmeri, and A. tuberculatus or the between-species transfer in the interfertile Ambrosia artemisiifolia and A. trifida or Lolium rigidum and L. perenne (Vincent and Cappadocia, 1987 ; Tonsor, 1990 ; Balfourier et al., 1998 ; Wetzel et al., 1999 ). We previously reported that the within-species glyphosate resistance transfer, based on the proportion of resistant individuals in the progeny of C. canadensis plants, ranged from 0% to 14% and 92% to 100% under pollen competition and no-pollen competition studies, respectively (Zelaya et al., 2004 ). Introgressive hybridization (introgression) in Conyza is well documented in the European species; however, the existence of Conyza hybrid zones in the Americas is unknown (Knobloch, 1972 ; Stace, 1975 ; McClintock and Marshall, 1988 ; Thébaud and Abbott, 1995 ). Considering the importance of glyphosate as a global herbicide and the pervasive nature of Conyza worldwide, an investigation was undertaken to assess the potential for glyphosate resistance transfer from C. canadensis to C. ramosissima through hybridization. Postzygotic reproductive barriers and phenotypic characterization of the interspecific Conyza hybrid (

Formula

) were also determined because these factors may affect the fitness and potential of hybrids to establish as important agricultural weeds. A naturally occurring Conyza hybrid with introgressed resistance to glyphosate would likely have immediate and considerable economic implications to United States agriculture. The potential implications of hybridization on management of glyphosate resistance and taxonomy of Conyza are discussed.

MATERIALS AND METHODS

Plant materials
The glyphosate-susceptible C. ramosissima population (

Formula

) was collected from the Ontario Cemetery, in Ames, Iowa, in June 2003. According to city records, no glyphosate was used on the cemetery premises. The original C. canadensis glyphosate-resistant population (

Formula

) was that collected by Mark VanGessel in Delaware (VanGessel, 2001 ). We previously reported that glyphosate resistance in this population was governed by an incompletely dominant, single nuclear gene (R allele) (Zelaya et al., 2004 ). A stable, near-homozygous glyphosate-resistant C. canadensis population (RS2) was isolated in the greenhouse through two cycles of recurrent selection on

Formula

at the rate of 2.0 kg acid equivalent (AE) glyphosate/ha (Zelaya et al., 2004 ). Complete specimen sets of

Formula

, RS2, and the interspecific hybrid progenies (explained later) were placed in the Ada Hayden Herbarium (ISC; accession nos. 435636–435644) (Appendices S3–S9, see Supplemental Data accompanying online version of this article). Duplicate specimen sets were also deposited at the Botanical Research Institute of Texas (BRIT) and the New York Botanical Garden (NY).

Growth conditions
Twenty

Formula

plants were transplanted from the field into 12-cm diameter pots with a peat : perlite : loam (1 : 2 : 1) soil mix media. In contrast, 20 RS2 individuals originated from seeds that were germinated in flats containing the soil mix media and later transplanted to 12 cm diameter pots. All crosses were performed in growth cabinets set at a 16-h photoperiod, 35° : 25°C day : night, 70–90% relative humidity (RH), and 600 µmol·m–2·s–1 photosynthetic photon flux density (PPFD) conditions. The RS2,

Formula

, and hybrid progenies (explained next) were grown in a greenhouse set to 25°–35°C and 50–80% RH diurnal conditions and 20°–25°C and 50% RH nocturnal conditions; natural light was supplemented with 600–1000 µmol·m–2·s–1 PPFD artificial illumination and photoperiod set to 16 h. Pots were irrigated as needed and fertilized (Miracle Gro Excel, Scott-Sierra Horticultural Products, Marysville, Ohio, USA) 1 mo after establishment. Achene production per capitula and plant, in addition to the germination rate of seeds, were estimated after crosses and in subsequent generations (explained later).

Conyza interspecific hybridization studies
The gynomonoecious C. canadensis and C. ramosissima possess white pistillate ray florets in the capitulum periphery and yellow perfect disk florets in the capitulum core. Pollen release may occur prior to capitula opening; therefore, both pollen competition and no-pollen competition studies were conducted pre-anthesis.

Response of parents to glyphosate
Prior to conducting crosses, the phenotype of parents was assessed to verify that RS2 and

Formula

were resistant and susceptible to glyphosate, respectively. The RS2 parents were treated with 2.0 kg AE of glyphosate/ha (2.3 times the recommended rate), and survival was assessed 20 d after treatment (AT); resistant parents had no phytotoxicity symptoms. While the phenotype of susceptible parents was not confirmed prior to conducting crosses, Ames city records indicated that the

Formula

population had never been exposed to glyphosate. Furthermore, treatment of the progeny of selfed

Formula

plants with 0.85 kg AE of glyphosate/ha resulted in uniform plant death (data not shown).

Pollen competition studies
Reciprocal crosses were performed between C. canadensis and C. ramosissima plants (

Formula

x RS2, RS2 x

Formula

) by placing the intact inflorescences of both species inside a PQ218 DelNet bag (DelStar Technologies, Middletown, Delaware, USA) and fastening bags at the inflorescence base with a wire. The mesh size in DelNet bags allowed for airflow while preventing the entry of foreign pollen. Under these conditions, the receptive stigmas of each species could only be fertilized by pollen of either species, thus the term "pollen competition". Fertilized florets were permitted to mature on the mother plant, capitula were harvested, seeds germinated on flats, and the emerged seedlings transplanted to 12-cm pots. Pollen competition tests assessed the level of allogamy between the evaluated Conyza species under controlled conditions. Under field conditions, however, the allogamy levels estimated by this test may increase or decrease depending on factors such as plant vicinity, insect pollinators, or wind speed and direction.

No-pollen competition studies
Disk florets of

Formula

parent plants were manually removed with forceps (emasculation). Upon stigma protrusion of the remaining

Formula

pistillate florets, fertilization was accomplished by gently rubbing the intact capitula of the pollen donor RS2 parent plants. Crosses were performed in one direction with

Formula

and RS2 serving as the pollen receptor and pollen donor parent, respectively. Non-emasculated capitula in

Formula

plants were removed to prevent autogamy. Approximately 50 capitula per

Formula

parent plant were emasculated and fertilized daily for 1 wk. When the pappus became visible, capitula were harvested from plants, seeds germinated in flats, and the resultant seedlings were transplanted to 12-cm pots. This test assessed the level of compatibility between the studied Conyza species in the absence of pollen competition.

Assessment of emasculation efficiency
Ten capitula per

Formula

parent were emasculated and permitted to develop on inflorescences covered with DelNet bags, thus preventing fertilization from external pollen. Once the emasculated capitula reached maturity, capitula were harvested, seeds were planted on flats, and the germination was monitored. If emasculation was completely effective in preventing autogamy in Conyza, no viable seeds would be produced.

Crossing scheme
Ten RS2 and

Formula

plant pairs (families) were crossed in isolation for the reciprocal pollen competition studies, and an identical arrangement was used for the unidirectional no-pollen competition studies. Estimates of hybridization were based on characterization of the first filial generation (F1) and quantification of interspecific hybrid (

Formula

) frequencies. One

Formula

plant per family was then permitted to self-pollinate in isolation, and segregation of the R allele was monitored in the resultant generation (

Formula

); 297 individual

Formula

plants were evaluated in the segregation analysis to develop a genetic model. Finally, introgression of the R allele was confirmed by backcrossing one

Formula

plant per family to the original RS2 (

Formula

) and

Formula

(

Formula

) parents. The parents used in backcrosses originated from florets of the original

Formula

or RS2 parents that were allowed to self-pollinate and produce seeds. The phenotype of RS2 parents used in backcrosses was confirmed as indicated previously. The

Formula

parents were confirmed susceptible by treating plants with the sublethal glyphosate rate of 0.4 kg AE/ha; this caused ≤60% visual injuries but did not kill

Formula

plants, thus allowing for nondestructive verification of the susceptible phenotype.

Postzygotic reproductive barriers
Seed viability was tested according to the Association of Official Seed Analysts (AOSA) standard germination procedure recommended for Asteraceae (AOSA, 2003 ). One hundred seeds per each of the 10 parents (RS2 and

Formula

) or families (

Formula

,

Formula

,

Formula

, and

Formula

) were germinated (n = 1000) on deionized-water-moistened blue blotter circles (Anchor Paper, St. Paul, Minnesota, USA) inside plastic petri dishes in the aforementioned growth cabinet conditions. Dishes were monitored daily for 2 wk, and germinated seeds were counted and removed; the remnant seeds were then classified as dormant or nonviable based on tetrazolium test results (Moore, 1985 ). Seed viability experiments were repeated in time 1 wk after the initial assessment (n = 2000). For postzygotic reproductive barriers, hybrid inviability was tested by assessing viability of

Formula

seeds, hybrid sterility tested on the

Formula

seed, and hybrid breakdown on the viability of

Formula

and

Formula

seeds.

Characterization of interspecific hybrids
Ten randomly selected RS2 and

Formula

plants or one randomly selected

Formula

plant per each of the 10 families generated through no-pollen competition crosses were sampled to determine quantitative traits. Five measurements (samples) were determined per plant, and measurements were repeated in time (n = 100). Rosette measurements were taken 5–6 wk after emergence. These measurements included rosette diameter, rosette leaf number, leaf length, leaf width, leaf shape (leaf length ÷ leaf width), leaf dentation, adaxial leaf trichomes, leaf area, leaf mass, and the specific leaf area (SLA; leaf area ÷ leaf mass). Two weeks after stem differentiation but prior to anthesis, cauline measurements included leaf length, leaf width, leaf shape, leaf dentation, adaxial leaf trichomes, leaf area, leaf mass, and SLA, branch number per plant, and stem diameter at 2 cm from the soil or at the axis of the rosette. Quantitative traits were also recorded post-anthesis but prior to seed physiological maturity and included capitula length, capitula width, capitula shape (capitula length ÷ capitula width), the number of pistillate and perfect florets, and total florets per capitulum. Only the achene length parameter was measured at maturity. All dimensional measurements were recorded with a digital caliper; leaf area was estimated by outlining leaves on millimetric paper (2-mm divisions) and calculating the area therein. Trichome numbers were determined using a stereoscope by counting the pubescence of leaf disks (30 mm2) excised with a hole punch.

Response of interspecific Conyza hybrids to glyphosate
The recommended glyphosate rate for Conyza at the 10-cm diameter rosette stage is 0.85 kg AE/ha (Roundup UltraMAX, Monsanto, St. Louis, Missouri, USA) (Anonymous, 2004 ). The resistant (R) and intermediate-resistant (IR) phenotypes comprised those rosettes that developed ≤30% and 31–69% visual injury 20 d AT, respectively, when treated with 2.0 kg AE of glyphosate/ha (Zelaya et al., 2004 ). After treatment with glyphosate, both R and IR phenotypes reached reproductive stage; however, only the R phenotype had visual growth rates equivalent to those observed in untreated Conyza plants. The susceptible (S) phenotype had ≥70% visual injury and was killed at the 2.0 kg AE/ha glyphosate rate.

Whole-plant rate response
Plants were treated with deionized water (dH2O; control) or 0.5, 1, 2, 4, 8, or 16 times the recommended glyphosate rate of 0.85 kg AE/ha. Treatments were applied 30 cm above the plant canopy through an 80015-E nozzle (TeeJet Spraying Systems, Wheaton, Illinois, USA) in a CO2-powered spray chamber (SB5–66, DeVries Manufacturing, Hollandale, Minnesota, USA) delivering 187 L/ha at a pressure of 2.8 kg/cm2. Treatments had four replications per rate and were repeated once in time (n = 8). Glyphosate phytotoxicity symptoms included plant stunting, leaf chlorosis, and necrosis that developed from the meristems and leaf tips to the rest of the plant. At 20 d AT, glyphosate efficacy was visually estimated by comparing glyphosate-treated Conyza plants with the dH2O-treated control plants (0% = non-injured; 100% = completely necrotic). Plants were then cut at the soil surface, placed in paper bags, and dried at 80°C for 48 h. Biomass was estimated by weighing the individual Conyza sample per 12 cm diameter pot and used to determine the glyphosate rate that inhibited plant growth by 50% (GR50); meristems subsamples were also taken for shikimic acid determination based on a method previously reported for Conyza (Zelaya et al., 2004 ).

Statistical analysis
Statistical Analysis Software (SAS, 2000 ) was utilized to conduct data analyses. Seed viability tests, arranged in a complete randomized design (CRD), were subjected to analysis of variance (ANOVA; PROC GLM) as well as mean separation by Fisher's least significant difference (LSD{alpha}=0.05) when ANOVA identified significant taxon effects. Glyphosate rate response tests were analyzed as a randomized complete block (RCB). GR50, I50 (glyphosate rate resulting in 50% accumulation of the maximum estimable shikimic acid), LD50 (glyphosate rate inflicting 50% mortality within the population), |{lambda}50| (the absolute difference between two estimated GR50 values), and {chi}2 goodness-of-fit estimates were done as previously reported (Zelaya et al., 2004 ).

Quantitative trait data were tested for normality based on the univariate Shapiro-Wilk test and accepted if the P value for W100 was ≥ 0.05; otherwise, the variance ({sigma}2) of the data were normalized by natural-log transformation (Shapiro and Wilk, 1965 ). ANOVA was done on individual quantitative traits considering taxon and family nested within taxon as fixed and random effects, respectively. Fisher's LSD{alpha}=0.05 tested for differences between taxa means. When the normality assumption was not met, Kruskal-Wallis analyses (PROC NPAR1WAY) were conducted, and post hoc non-parametric mean separation was performed by Dunn's test (Bonferroni's method) (Conover, 1999 ). In addition, PROC VARCOMP was used to partition the total phenotypic {sigma}2 into the different ANOVA components for traits that converged to Shapiro-Wilk's normality assumption; these {sigma}2 components were then used to calculate the intraclass correlation coefficient (t):


Formula 1

(1)
where

Formula 1

and

Formula 1

corresponded to the {sigma}2 between and within taxon, respectively.

Phenotypic intermediacy of the

Formula 1

, as it related to the RS2 and

Formula 1

parents, was tested by a character count procedure using trait means based on an intermediate vs. non-intermediate one-tailed sign test (Wilson, 1992 ). Multivariate normality was tested by Mardia's kurtosis ({gamma}2) analysis (Mardia, 1970 ). PROC PRINCOMP was invoked to test whether the data covariance matrix was singular. If rejected, squared Mahalanobis distances (

Formula 1

) from the centroid were computed for chi-square ({chi}2) quantile–quantile (Q–Q) comparisons of multivariate normal data; 75% confidence intervals were constructed from standard deviation ({sigma}) estimates of g(z), as indicated by Chambers et al. (1983) . Quantitative traits that converged to Mardia's assumption were then subjected to canonical discriminant function (CDF) analysis utilizing PROC CANDISC (Thompson, 1984 ). Taxon clustering of group averages used Mahalanobis' distance matrix based on the unweighted pair-group method with arithmetic mean (UPGMA) (PROC CLUSTER); concomitantly, PROC TREE was used to display the phylogenic relationship between taxon.

RESULTS

The Conyza parental populations differed in their response to glyphosate
Nonsignificant lack-of-fit (LOF) tests and coefficients of determination estimates for biomass (F = 0.33; P = 0.98;

Formula 1

= 0.82) and shikimic acid (F = 1.36; P = 0.15;

Formula 1

= 0.92) confirmed suitability of the log-logistic model for describing plant response as a function of increasing glyphosate rates.

Formula 1

was killed at 0.85 kg AE of glyphosate/ha (Figs. 1 and 2). Compared to

Formula 1

, RS2 had a three fold and seven fold increase in estimated GR50 and LD50 values, respectively (Table 2). In addition, when comparing GR50 estimates for

Formula 1

and RS2, a statistically significant |{lambda}50| value of 1.60 (Fobs = 1.77; P < 0.01) was observed, reaffirming that the parental populations differed in their response to glyphosate. Shikimic acid determinations, which indirectly estimate the level of EPSPS inhibition by glyphosate (Harring et al., 1998 ), suggested that twice the rate of glyphosate was required to inhibit 50% of EPSPS in RS2 plants compared to

Formula 1

(Table 2). Untreated Conyza plants contained 7–14 µmol of shikimic acid/g dry mass, which increased sigmoidally to a maximum of 150 µmol at the two highest glyphosate rates (Fig. 1). Concurrently, endogenous shikimic acid levels correlated negatively with plant biomass and positively with visual injury, thus further supporting the differential response to glyphosate observed between the parental populations (Table 2).


Figure 1
Figure 1
Figure 1
Figure 1
Figure 1
Figure 1
View larger version (15K):
[in this window]
[in a new window]

 
Fig. 1. Main plot: N-phosphonomethyl glycine (glyphosate) whole-plant rate response of Conyza canadensis (RS2), C. ramosissima (), the interspecific hybrid (), and the hybrid progeny () 20 d after herbicide treatment. Insert: Endogenous shikimic acid accumulation in the tissue of treated plants. For either graph, data points represent the mean of four replications and two experiments (n = 8). Three randomly selected plants per each of the 10 families (n = 30) were used in and determinations. Extensions on symbols designate the standard error associated with individual means ({sigma}M).

 

Figure 2
Figure 2
Figure 2
Figure 2
Figure 2
Figure 2
Figure 2
Figure 2
Figure 2
View larger version (21K):
[in this window]
[in a new window]

 
Fig. 2. Main plot: Mortality observed 20 d after application of N-phosphonomethyl glycine (glyphosate) on Conyza canadensis (RS2), C. ramosissima (), the interspecific hybrid (), and the hybrid progeny (). Each symbol represents the mean of four replications and two experiments (n = 8); the and mortality was estimated from three randomly selected plants per each of the 10 generated families (n = 30). Solid and broken lines represent the mortality estimated by PROBIT and that expected considering partially dominant monofactorial inheritance in the (Tabashnik, 1991 ), respectively. Insert: Biomass of untreated (white bar), (gray bar), and RS2 (black bar) rosettes at the 10–12-cm diameter stage. Bars represent the mean of 10 randomly selected rosettes and determination repeated in time (n = 20). Letters above bars represent minimum statistical differences according to Fisher's LSD for comparisons between Conyza taxa means (LSD{alpha}=0.05 = 0.24 g); extensions on bars or symbols represent the standard error associated with individual means ({sigma}M).

 

View this table:
[in this window]
[in a new window]

 
Table 2. Response of Conyza canadensis (RS2), C. ramosissima (), the interspecific hybrid (), and the hybrid progeny () to glyphosate. Numbers in parentheses designate the 95% lower and upper confidence intervals (GR50; I50) or fiducial limits (LD50) for the preceding estimated parameter. The extent of association between shikimic acid levels and biomass or visual injury is reported according to Spearman's correlation (r2) analysis.

 
Levels of interspecific hybridization
Pistillate florets in emasculated

Formula 1

capitula contained <1% viable seeds in the absence of pollen (Table 3). Previously, we proposed that fertilization of pistillate florets under these conditions probably originated from pollen of perfect florets released prior to anthesis or from pollen of perfect florets incompletely excised during emasculation (Zelaya et al., 2004 ). Under greenhouse environments, approximately 60% of the total achenes on

Formula 1

plants produced viable seeds (Fig. 3). Therefore, we estimated that the emasculation method was 98% effective at preventing self-fertilization in

Formula 1

.


View this table:
[in this window]
[in a new window]

 
Table 3. Frequencies of hybridization in the Conyza canadensis x C. ramosissima cross. Emasculation efficiency was tested by evaluating the germination of seeds that developed in the absence of pollen in 10 emasculated capitula per each of 10 C. ramosissima () parents evaluated. Crosses between C. canadensis (RS2) and were conducted under pollen and no-pollen competition conditions. Hybridization was estimated from the frequency of Conyza plants with a susceptible (S), resistant (R), or hybrid (H) phenotype within the first filial generation (); data represent the sum of the 10 generated families (see Materials and Methods).

 

Figure 3
Figure 3
Figure 3
Figure 3
Figure 3
Figure 3
Figure 3
Figure 3
Figure 3
Figure 3
Figure 3
Figure 3
Figure 3
View larger version (20K):
[in this window]
[in a new window]

 
Fig. 3. Proportion of germinated (gray bars), dormant (dark gray bars), and nonviable (white bars) seeds among the Conyza canadensis (RS2), C. ramosissima (), interspecific hybrid (), hybrid progeny (), and the backcross to RS2 () or (). Bars represent the mean of two experiments comprised of 100 seeds per each of the 10 parents (RS2 and ) or 10 families (, , , and ) generated in the no-pollen competition studies (n = 2000); extensions above bars designate the standard error associated with individual means ({sigma}M). Letters above bars represent minimum statistical differences according to Fisher's LSD for comparisons within germinated (LSD{alpha}=0.05 = 12%), dormant (LSD{alpha}=0.05 = 7%), and nonviable (LSD{alpha}=0.05 = 14%) seeds.

 
Interspecific hybridization under pollen competition conditions ranged from 0% to 7% in the RS2 to

Formula 1

crosses and from 0% to 9% in the reciprocal

Formula 1

to RS2 crosses, although five in 10

Formula 1

and two in 10 RS2 families contained no identifiable Conyza hybrids (data not shown). In contrast, high hybridization levels (98.2%) were observed under no-pollen competition conditions when RS2 served as pollen donor to

Formula 1

(Table 3). The observed levels of glyphosate-susceptible individuals (<2%) in the no-pollen competition studies mirrored those of emasculation efficiency estimates (2%), suggesting that some level of self-fertilization occurred prior to the observed anthesis. Collectively, these data suggested that while the levels of hybridization were relatively low (2.6% to 4.1%) under pollen competition conditions, C. canadensis and C. ramosissima are genetically compatible (>95%), and thus the potential for interspecific hybridization in the field exists. Anecdotal herbaria records of "off-type" Conyza sp. suggest that hybridization does occur, albeit infrequently, in natural plant communities.

Phenotypic variance of parents and hybrid progenies
Differences among families within taxa were significant for approximately half of the quantitative traits evaluated (Table 4). Maternal family effects were significant (P < 0.05) in nine of 27 and 10 of 27 quantitative traits within the

Formula 1

and

Formula 1

populations, respectively, and in approximately half of the RS2 traits evaluated. Only the characters rosette leaf length, cauline leaf trichomes, and achene length had significant maternal family effects across RS2,

Formula 1

, and

Formula 1

(Table 4). Intraclass correlation coefficients (t) were generally (85%) different from zero and most (72%) were greater than 0.90 (Table 4). Combined with the ANOVA, data suggested that most of the observed phenotypic variance ({sigma}2) was accounted for differences between the studied Conyza taxa. Furthermore, we assumed that this {sigma}2 was primarily attributable to genetic differences between the studied taxa rather than to the interaction with the environment because the Conyza populations developed under controlled greenhouse conditions (Appendix S1, see Supplemental Data accompanying online version of this article).


View this table:
[in this window]
[in a new window]

 
Table 4. Partitioning of the total phenotypic variance ({sigma}2) into taxon and family effects for the evaluated quantitative traits. Numbers represent F values and probability (in parentheses) estimates for the analysis of variance (ANOVA). The total {sigma}2 between and within the Conyza taxon was utilized to calculate the intraclass correlation coefficients (t) for each trait; underlined coefficients correspond to traits with a significant family effect (P < 0.05) in the ANOVA for the individual taxa. Nomenclature: Conyza canadensis (RS2), C. ramosissima (), and the interspecific hybrid ().

 
Phenotype of

Formula 1


Statistical differences (P < 0.05) between the studied Conyza taxa were obtained for all 27 quantitative traits evaluated (Table 5). RS2 produced wider rosettes with fewer leaves than

Formula 1

. Furthermore, both RS2 leaf dimorphisms were longer and wider and possessed more dentations than

Formula 1

(Table 5). In contrast,

Formula 1

was profusely branched, had narrower stems, and developed more pubescent leaves with greater density (SLA) than RS2. Capitula dimensions and the number of individual and total florets were greater in RS2; however,

Formula 1

produced longer achenes (Table 5).


View this table:
[in this window]
[in a new window]

 
Table 5. Comparison of quantitative traits among Conyza canadensis (RS2), C. ramosissima (), and the interspecific hybrid (). Values represent the mean of five observations per each of 10 sampled plants and assessments repeated in time (n = 100); numbers in parentheses designate the standard error associated with means ({sigma}M). A one-tailed sign test of intermediate vs. non-intermediate was used to test the hypothesis of phenotypic intermediacy.

 
To investigate whether the phenotype of

Formula 1

plants was intermediate to both parents, a character count one-tailed sign test was conducted. Rosette leaf number, cauline branch number, cauline SLA, and capitula length

Formula 1

measurements tended to exceed those of either parent; however, only capitulum shape was significantly different (Table 5). Hence, five of the 27 evaluated quantitative traits did not meet the character count assumption that

Formula 1

was intermediate to both parents (D = 22). Because conditions for the DeMoivre–Laplace theorem were met (np and nq ≥ 5), normal approximation was used for probability estimation rather than a binomial distribution. The calculated one-tailed sign test value for phenotypic intermediacy was significant (z+ = 3.27; P = 0.001) and thus prompted rejection of the null hypothesis (H0) for equal sign proportions of the evaluated quantitative traits. The collective phenotypic data therefore suggested that

Formula 1

represented an intermediate between both parents.

Intermediacy of

Formula 1


Multivariate normality analysis determined that only 17 of the 27 evaluated morphometric traits converged to Mardia's statistics ({gamma}2 = 328.2; P = 0.08) and thus were combined for multivariate discriminant analysis; the pistillate florets trait was eliminated from the analysis as inclusion instigated a singular covariance matrix. The H0 that data arose from populations with a common distribution was tested by Q–Q plots; linearity of data points along the expected mean vectors and allocation within the estimated 75% confidence intervals suggested nondeparture from normality (Fig. 4). Most (81%) of intraclass correlation coefficients (t) in the multivariate normally distributed morphometric traits possessed values greater than 0.70, hence suggesting negligible interference of correlations among families that could potentially distort interpretation of the canonical discriminant function (CDF) analysis (Table 4).


Figure 4
Figure 4
Figure 4
View larger version (12K):
[in this window]
[in a new window]

 
Fig. 4. (A) Multivariate discriminant function analysis of canonical (CDF) variates 1 and 2 for 17 morphometric traits that converged to Mardia's multivariate normality assumption. Data points represent the score of individuals (n = 100) within the Conyza canadensis (RS2), C. ramosissima (), and interspecific hybrid () populations. Inserts: Chi-square ({chi}2) quantile–quantile (Q–Q) plots of taxon observations. The continuous and dashed lines represent expected mean vectors and 75% confidence intervals for the data, respectively. (B) Clustering analysis of Mahalanobis' distance matrix based on taxon averages by the unweighted pair-group method with arithmetic mean (UPGMA).

 
Univariate statistics for differences among taxa means were significant (P < 0.0001) for all traits, as were Wilks' exact multivariate statistics for comparison among taxa ({lambda} = 0.004; F = 229.9; P < 0.0001). The CDF analysis estimated that the first (Can1) and second (Can2) canonical variates accounted for 97% and 82% of the total quantitative trait {sigma}2, respectively. Taxa clustered in discrete sections within the canonical graph;

Formula 1

(6.8) and the

Formula 1

(1.3) were positive along the Can1 axis, while RS2 was diametrically negative (–8.2). Fewer taxa divergence was obtained along the Can2 axis (Fig. 4; Appendix S2, see Supplemental Data accompanying online version of this article).

Taxa separation was also discernable by group averages based on the UPGMA. Divergence in Mahalanobis distances was greatest among the RS2 and

Formula 1

parents (

Formula 1

= 226.5; F = 630.2; P < 0.0001), and approximately half and one-fourth that square distance was estimated for comparisons between the

Formula 1

and RS2 (

Formula 1

= 107.5; F = 299.1; P < 0.0001) and

Formula 1

and

Formula 1

(

Formula 1

= 53.7; F = 149.4; P < 0.0001), respectively. Hierarchical clustering formed a dendrogram with a common root node for the taxon (Fig. 4). The

Formula 1

and

Formula 1

leaves clustered within a single branch, whereas RS2 clustered to the common parent node. While the character count procedure lent credence to the thesis of an intermediate hybrid

Formula 1

phenotype, CDF analysis and UPGMA strongly emphasized that the

Formula 1

was more similar to

Formula 1

than to RS2.

Postzygotic reproductive barriers
Partial correlation estimates (r2 = 0.76; P = 0.001) associated with the sums of squares and crossproducts (SSCP) matrix suggested a strong relationship between the seed viability experiments repeated in time. Univariate (F = 0.14; P = 0.71) and multivariate (Wilks' {lambda} = 0.99; P = 0.71) tests for the between-time effects were not significant, therefore no difference in experiments repeated in time was inferred and data for the seed viability experiments were combined in further analyzes. Previous research reported near 100% C. canadensis germination under light and constant 28°C conditions (Shontz and Oosting, 1970 ). Viability of RS2 seeds under our conditions ranged from 17% to 71% within families with a mean of 46%, which was not statistically different (LSD{alpha}=0.05 = 12%) from the 57% mean viability estimated for

Formula 1

seeds (Fig. 3).

Evidence for hybrid inviability was apparent in the 23% increase in nonviable

Formula 1

seeds (Fig. 3). No-pollen competition tests estimated hybridization levels of 98% in the viable hybrid zygotes within families (Table 3). Viability tests nonetheless confirmed that approximately half of the available ova in

Formula 1

were not fertilized by RS2 pollen or ova were fertilized but aborted prematurely, an argument that some level of genetic incompatibility existed among the Conyza parents (Fig. 3). While physical damage of

Formula 1

capitula during emasculation could have contributed to the inviability of hybrid zygotes, epistatic, homeotic transformations, or interactions of complementary genes may have also resulted in embryo abortion of

Formula 1

zygotes.

The presence of a hybrid sterility reproductive barrier was disregarded because the

Formula 1

had embryos with approximately twice the level of viability of their progenitor (

Formula 1

) (Fig. 3). Estimates of

Formula 1

viability were significantly less than for either RS2 or

Formula 1

, an indication that some level of self-infertility remained within hybrids and evidence that the inviable

Formula 1

zygotes probably arose from genetic incompatibilities rather than the physical stress of emasculation (Fig. 3). Backcrosses of

Formula 1

to RS2 or

Formula 1

produced a viable progeny, 21% and 25%, respectively, confirming successful introgression of the R allele to either of the evaluated Conyza taxa through the intermediate hybrid (

Formula 1

). Establishment of hybrid interspecific populations in the environment is often limited by postzygotic reproductive barriers (Barton, 2001 ). While the germination of

Formula 1

seeds was less than that of the RS2 or

Formula 1

parents, the collective viability data established that the Conyza hybrid was fertile and probably capable of establishment in natural environments.

Transgressive segregation in

Formula 1


Most hybrid progenies demonstrate marked traits of transgression (Rieseberg et al., 1999 ). The majority of

Formula 1

plants appeared similar to either the canadensis or ramosissima epithets. However, specific characteristics attributable to transgression were observed, ranging from glabrous to pubescent leaves and stems, oblanceolate to subulate leaves with dentated to strait margins, green to purplish midribs, irregular arrangements of lateral leaf nervures, and profusely branched single-stem plants. Generally, flower morphology was a preserved trait among

Formula 1

individuals. Transgressive segregation was also discernable in whole-plant rate responses to glyphosate. Compared to the biomass accumulated by RS2 ({sigma}2 = 0.67),

Formula 1

({sigma}2 = 0.42), or the

Formula 1

({sigma}2 = 0.87),

Formula 1

rosettes demonstrated greater variance ({sigma}2 = 1.17) in response to glyphosate. Similarly, greater variations in the level of endogenous shikimic acid accumulation and visual injury were observed in the

Formula 1

(data not shown). The

Formula 1

also was characterized by the unexpected death of plants (<5%) at the early rosette stages. We theorize that these low lethal frequencies reflect possible deleterious homeotic transformations or epistatic interactions in

Formula 1

plants.

Phenotype of backcrosses
Visual observations confirmed that the progeny of

Formula 1

backcrosses to parents (

Formula 1

and

Formula 1

) were phenotypically similar to either C. canadensis or C. ramosissima. For example,

Formula 1

resembled RS2 at the rosette and cauline stages, except that the plant developed lateral branches below the center-main stem. In addition,

Formula 1

produced serrated and nonserrated leaves, while RS2 produced only the serrated dimorphism.

Formula 1

possessed a well-defined axis that was absent in

Formula 1

and generated slightly larger leaves and thicker stems than

Formula 1

.

Formula 1

was also precocious compared to

Formula 1

, completing the reproductive cycle in 2–3 mo.

Inheritance of glyphosate resistance
Application of 0.40 kg AE of glyphosate/ha to

Formula 1

resulted in visual injury levels of 30–60% and effectively delayed plant development compared to the dH2O-treated Conyza plants, although treated plants recovered from glyphosate injuries within 2–4 wk and reached reproductive stage. Rates of 0.85 kg AE of glyphosate/ha resulted in uniform kill of

Formula 1

or the selfed progeny of

Formula 1

, suggesting that the population was near-homozygous susceptible to glyphosate. We previously confirmed that RS2 represented a near-homozygous glyphosate-resistant population (Zelaya et al., 2004 ).

Overdominance of

Formula 1


The

Formula 1

produced larger rosettes compared to

Formula 1

and more and denser (SLA) leaves than RS2 (Table 5). The difference in these two leaf parameters resulted in an apparent heterotic response of

Formula 1

recorded at the 10–12-cm diameter rosettes stage (Fig. 2). Compared to

Formula 1

and RS2, greater glyphosate rates were required to reduce biomass accumulation or cause mortality in

Formula 1

(Table 2; Figs. 1 and 2). This divergence in response to glyphosate was confirmed by |{lambda}50| values for comparisons of

Formula 1

and RS2 (Fobs = 1.74; P < 0.01) or

Formula 1

and

Formula 1

(Fobs = 1.76; P < 0.01). Not only did

Formula 1

rosettes demonstrate a vigorous growth that probably required greater glyphosate rates to inhibit, but leaves produced more trichomes than RS2, which could have hindered glyphosate absorption into

Formula 1

plants (Table 5). The response of

Formula 1

to glyphosate therefore failed to obey the additive, dominant, or hybrid susceptibility models that describe resistance in hybrid populations; rather, the response was best explained by the hybrid resistance hypothesis, which predicts greater resistance in hybrid populations compared to their progenitors (Fritz et al., 1994 ).

Segregation of the R allele
Glyphosate resistance in RS2 is conferred by the incompletely dominant, nuclear R allele (Zelaya et al., 2004 ). This model for glyphosate resistance in the Conyza hybrid was tested by monitoring the segregation ratios of

Formula 1

,

Formula 1

, and

Formula 1

to glyphosate. Efficacy tests at 20 d AT with 2.0 kg AE of glyphosate/ha, a rate that differentiated R and S phenotypes in the parental populations, identified three distinct segregates in

Formula 1

families—R, IR, and S phenotypes. Exact goodness-of-fit (GOF) tests based on the H0 that the observed segregation ratios followed a 1 : 2 : 1 Mendelian distribution provided nonsignificant {chi}2 values for all

Formula 1

families, substantiating appropriateness of the partially dominant monogenic model (Table 6). Concomitantly, homogeneity analysis ({chi}2 = 5.72; P = 0.77) confirmed that the combined

Formula 1

data converged to the 1 : 2 : 1 genetic model.


View this table:
[in this window]
[in a new window]

 
Table 6. Segregation for glyphosate resistance of the R allele in the progeny () of one interspecific hybrid plant () per family allowed to self-pollinate in isolation. The 10 families evaluated arose from a single plant isolated from the Conzya canadensis (RS2) to C. ramosissima () no-pollen competition cross. For back-crosses, served as pollen donor to RS2 () or ().

 
Backcrossing of

Formula 1

to RS2 (

Formula 1

) and subsequent progeny treatment with 2.0 kg AE of glyphosate/ha identified only R and IR phenotypes (Table 6). The observed segregation ratios among

Formula 1

families were homogenous ({chi}2 = 4.40; P = 0.88), and the collective data were consistent with the expected 1 : 1 (R : IR) ratio. Similarly, analysis of

Formula 1

backcrosses to

Formula 1

(

Formula 1

) resulted in homogenous IR and S segregation ratios ({chi}2 = 5.81; P = 0.76) that obeyed the partially dominant model (Table 6). Because the backcross data followed the monofactorial model of inheritance, results corroborated our previous assertion that RS2 and

Formula 1

represented near-homozygous lineages with regard to their response to glyphosate.

Further substantiation of the incompletely dominant, monogenic model was obtained graphically from the pattern of observed

Formula 1

mortality by comparing to that mortality expected as suggested by Tabashnik (1991) : Y{chi} = WR (0.25) + WIR (0.50) + WS (0.25). Three distinct response phases were predicted based on partially dominant, monofactorial inheritance (Fig. 2). No mortality was recorded from 0.0 to 0.42 kg AE/ha, suggesting that both homozygous (RR and rr) and the heterozygous (Rr) genotypes were present at these glyphosate rates. A second segment was discernable at glyphosate rates of 0.85–3.38 kg AE/ha, which corresponded to approximately one-fourth (23–37%) mortality of the putative homozygous susceptible genotype (Fig. 2). The heterozygous (71%) and homozygous (100%) resistant genotypes were controlled at glyphosate rates of 6.77 and 13.54 kg AE/ha, respectively.

DISCUSSION

Hybridization of Conyza in nature
The phylogenic boundaries in Conyza are not clearly understood and considerable phenotypic variation has been reported, particularly in response to adverse environmental stimuli (Nesom, 1990 ). We initiated a project to assess potential hybridization of Conyza species and better understand the phylogenic relationship between C. canadensis and C. ramosissima, two weedy species in United States agroecosystems despite nothing described in the literature. Our work herein suggests that the studied taxa are genetically compatible, capable of transferring the R allele, and producing interspecific hybrid progenies that are vigorous and fertile. Given that C. canadensis has become a major economic weed problem in the Midwestern United States and the apparent vigor of the hybrid identified in this research, the ramifications of an interspecific Conyza hybrid that has resistance to glyphosate are potentially significant. Glyphosate-resistant crop systems are suggested to be simple and without great environmental consequences. However, we have demonstrated that there are major ecological and economic consequences from these presumed simple systems. New weeds typically evolve over a long period of time, and existing weeds tend to adapt to new agroecosystems slowly. We propose that if hybridization of new taxa with glyphosate resistance as a semi-dominant trait can occur with relative ease, the current agroecosystem is at considerable jeopardy. While the hybrid demonstrated in this research has not been specifically identified in the field, we have illustrated the potential for the occurrence.

We speculate that the relatively high genetic compatibility among the studied taxa is associated with the common diploidy structure (2n = 18) that allows for successful chromosome pairing during meiosis. Both C. canadensis and C. ramosissima represent sibling species as ascertained by nrDNA internal transcribed spacers (ITS) analysis, which clustered both taxa in a single branch within group VI of the Erigeron and allied Asteraceae cladogram (Noyes, 2000 ). This phylogenic analysis estimated a recent speciation event between the ramosissima and canadensis epithets and provides further support to our thesis of genetic compatibility between the species.

The likelihood of native hybridization in Conyza is probably low given the enclosed involucre arrangement in Conyza and the autogamous nature of the genus. Entomophily was reported in C. canadensis, although the relative importance to interspecific gene transfer remains unknown (Weaver, 2001 ). To our knowledge, there is no prior documentation of hybridization between C. canadensis and C. ramosissima. However, since early descriptions, plants with common traits between both species have been cited (Michaux, 1803 ).

Voss (1996) reported depauperate C. canadensis biotypes in the upper peninsula of Michigan with glabrous to pubescent patterns and ramification near or below the main stem, characteristics that we observed in the

Formula 1

. In Canada, Darbyshire (1990) described profusely branched C. canadensis biotypes that resembled our

Formula 1

, but attributed those abnormalities to the loss of apical dominance. Because environmental factors can impact phenotype, some argue that C. canadensis and C. ramosissima are conspecific and that C. ramosissima represents a depauperate extreme of the highly plastic canadensis epithet (S. Darbyshire, Agriculture and Agri-Food Canada, personal communication). However, because C. canadensis and C. ramosissima plants grown in the greenhouse preserved their distinctive phenotypes, we suggest that these extreme phenotypes possibly represent transgressive hybrids. Interspecific gene transfer and the interactions of genetic and environmental stimuli affecting the phenotypic expression of Conyza deserve further investigation (Voss, 1996 ).

Documentation of hybridization in Conyzinae is at present restricted to Europe. For example, in the British Isles hybridization between C. canadensis and Erigeron acer L. was reported to produce weak and apparently sterile plants with few capitula; the nothotaxa was classified as Conyzigeron x huelsenii (Vatke) Rauschert (Stace, 1975 ). Furthermore, hybrids of unknown fecundity, namely Conyza x flahaultiana (Thell.) Sennen and Conyza x daveauiana Sennen in Spain and France, originated from crosses between C. canadensis and C. bonariensis (L.) Cronq. and between C. bonariensis and C. sumatrensis (Retz.) E. Walker, respectively (McClintock and Marshall, 1988 ). In the Iberian Peninsula, Conyza x rouyana Sennen arose from the hybridization between C. albida Willd. ex Spreng. and C. canadensis (J. L. Carretero, Universidad Politécnica de Valencia, personal communication). Additionally, Thébaud and Abbott (1995) stated that in France, C. sumatrensis and C. blakei Cabr. cross under natural environments and produce hybrid progenies with moderate (30%) fertility. Importantly, significant phenotypic variability exists in C. blakei (Laínz, 2001 ). In Belgium, Verloove and Boullet (2001) reported that some xenophyte Conyza populations identified as C. floribunda were in fact C. canadensis x C. sumatrensis hybrids. Furthermore, Conyza x mixta Fouc. & Neyr. reportedly arose from crosses between C. canadensis and C. bonariensis in Belgium, France, Great Britain, and Portugal (F. Verloove, National Botanic Garden of Belgium, personal communication). More recently, Sída (2003) reported a putative hybrid between C. bonariensis and C. triloba Decne. in the Czech Republic.

Loss of vigor is apparently a common trait to European Conyza hybrids. We speculate that ploidy differences may be a significant barrier determining successful hybridization in Conyzinae. For example, more compatible and vigorous hybrids would be expected from crosses between the allopolyploids (2n = 54) C. sumatrensis, C. floribunda, and C. bonariensis compared to crosses with the diploid (2n = 18) C. canadensis. Just as geography delimits the major groups in Asteraceae, spatial isolation is probably the principal species barrier within Asteraceae (Nesom, 1990 ). Therefore, disturbances of these barriers will likely stimulate interactions between phylogenetically related taxa separated by geography (allopatry) and allow for the exchange of genetic material between unstructured random-mating (panmictic) Asteraceae populations.

Ploidy of the Conyza hybrid
Hybridization represents an important mechanism for plant speciation whereby fertile and stable taxa can arise by either chromosome number doubling (allopolyploidy) or recombinational speciation without polyploidization (homoploidy) (Rieseberg et al., 2000 ). In the absence of karyological studies, the ploidy of the C. canadensis x C. ramosissima hybrid remains unknown. Nonetheless, we propose that the segregation and compatibility data allude to the formation of an interspecific hybrid without an increase in ploidy.

Man-made allopolyploids often display homeotic transformations and aberrant chromosomal rearrangements that result in gene silencing, hybrid instability, and lethality (Comai, 2000 ). These phenomena were absent from the Conyza hybrid, as the interspecific hybrids demonstrated a vigorous growth and marginal (<5%)

Formula 1

mortality (Figs. 2 and 3). Additionally, homeologous recombination in neo-allopolyploids typically destabilizes genomes by chromosomal deletions or expansions, and homeologous pairing often hinders meiosis through the formation of univalent or trivalent chromosomes (Comai, 2000 ). The genome of

Formula 1

was apparently stable as hybrids demonstrated moderate to no postzygotic reproductive barriers and were genetically compatible, with capacity for introgression with either the C. canadensis or C. ramosissima parent (Table 6; Fig. 3). Finally, allopolyploids often preserve homeologous integrity by hindering intergenomic recombination, whereas in homoploids, numerous chromosome combinations form through chiasmata, thus allowing recombination of the parental genomes (Comai, 2000 ). Given that RS2 and

Formula 1

represented near-homozygous lineages and that

Formula 1

families segregated for glyphosate resistance, the C. canadensis and C. ramosissima genomes may have coalesced in the

Formula 1

. This characteristic of segregation is typically associated with homeologous recombination in homoploids.

Establishment of the Conyza hybrid in the environment
The low hybridization levels (0% to 12%) in the pollen competition studies and the estimated high genetic compatibility between C. canadensis and C. ramosissima (98%) suggest that the resultant hybrid will probably evolve in nature as a contiguous population in a common geographic range (parapatry) (Table 3). Under natural environments, hybridization may occur at different frequencies than those observed in the greenhouse, depending on insect- and wind-pollination levels and other environmental factors such as competition and resource availability, which can affect flowering time in Conyza (Thébaud et al., 1996 ). Regardless, the fact that

Formula 1

was fertile implied that fitness and habitat divergence may be two essential factors determining the adaptation of Conyza hybrids to the environment. Certainly, the overdominance for resistance to glyphosate provides the proposed Conyza hybrid considerable ecological fitness, given the widespread use of glyphosate in current agroecosystems.

Models for homoploid speciation suggest that superior hybrid competitiveness results in the rapid displacement of the parental taxa, while under a fitness disadvantage, either the hybrid taxa becomes extinct or it coexists with the parents, provided adequate niche differentiation in the environment (Rieseberg, 1997 ). Hybrids are rarely better fit than their well-adapted congeners. Low initial frequencies, reduced fertility and viability, and competitive disadvantage with respect to parents typically lead to hybrid extinction (Wolf et al., 2001 ).

Ascertaining fitness of the Conyza hybrid and the potential ecological displacement of other taxa in different environments (vicariation) would require further experimentation and is not within the scope of this study. The supposition that

Formula 1

hybrids are highly fit compared to the Conyza parents is provisional, pending proper field assessments of fitness. Nevertheless, we believe that several hybrid traits would serve as ecological advantages, at least under current agroecosystems. The

Formula 1

rosettes may be competitive by generating more and denser leaves than RS2 and by producing wider rosettes than

Formula 1

(Table 5), thus accumulating greater biomass than either parent (Fig. 2). In addition,

Formula 1

plants possessed leaves with higher trichome densities than RS2, a trait that may reduce herbicide uptake and hinder herbicide efficacy. This could be significant because the mechanism of glyphosate resistance in C. canadensis was attributed to an altered cellular distribution, resulting in reduced herbicide translocation (Feng et al., 2004 ; Koger and Reddy, 2005 ; Dinelli et al., 2006 ). Furthermore,

Formula 1

plants produced multiple branches, which may facilitate recovery from herbicide applications through meristematic regrowth (Table 5).

We therefore suggest that the Conyza hybrid may be well fitted to agroecosystems, because evaluation of the

Formula 1

evidenced heterosis and overdominance for glyphosate resistance compared to RS2 and

Formula 1

(Table 2; Fig. 2). The

Formula 1

developed an equal number, if not more, achenes than either parent (Table 5); this reproductive behavior could compensate for the estimated depression in

Formula 1

seed viability (Fig. 3). Achene adaptation for dispersion by wind (anemochory), the capacity to overwinter (therophyte), and a prolific achene production make C. canadensis extremely competitive under no-tillage agroecosystems. However, the species is poorly adapted to crop production systems where tillage is used as a management strategy (Brown and Whitwell, 1988 ). Hence, the larger

Formula 1

achenes may provide an ecological advantage over C. canadensis by allowing Conyza hybrid seeds to emerge from deeper soil profiles in tilled agroecosystems (Table 5).

According to Barton (2001) , reasons for higher hybrid fitness include the (1) occurrence of diverse phenotypes through transgression that may have ecological advantages under particular environments, (2) reconstruction of an "ancestral linkage" that was competitive in the past and is probably adapted to present environments, and (3) coalescence of previously disjunct gene sets that may have a positive impact on fitness. Hybridization may therefore explain some cases of niche differentiation and provide the raw materials for the adaptation of novel weeds to the environment (Rieseberg et al., 1999 ). It is important to acknowledge that the transgressive fitness with respect to glyphosate resistance resulted from the combination of phenotypes of the two Conyza populations herein studied; a different fitness may be observed in hybrids originating from crosses between other Conyza species.

Implication of hybridization on herbicide resistance management
The impact of hybridization on glyphosate resistance management would depend on the native introgression levels of the R allele and the stability, fertility, and fitness of the Conyza hybrid with respect to RS2 and

Formula 1

. Herein, we demonstrated that the Conyza hybrid is fertile and competitive (Figs. 2 and 3). In addition, successful R allele transfer occurred from C. canadensis to C. ramosissima through the

Formula 1

in backcrosses (Table 6). Furthermore, the

Formula 1

demonstrated many phenotypic traits that may facilitate the adaptation to different agroecological niches. These circumstances could potentially complicate the management of glyphosate resistant weed populations in glyphosate resistant crops and the containment of glyphosate resistance genes within these agroecosystems.

Absent from this hypothesis, however, are the natural hybridization frequencies in Conyza. Native hybridization frequencies in the European Conyza approximate 50 plants in thousands of individuals, although 60% of the resultant hybrids were infertile (Thébaud and Abbott, 1995 ). These hybridization frequencies are completely dependent on the ecological circumstances under which plants developed. However, we would expect higher hybridization frequencies than those reported in Europe since the parents evaluated in this study are compatible and the resultant Conyza hybrid demonstrated negligible postzygotic reproduction barriers. Even low initial hybrid frequencies within the population could increase in time because Conyza plants can produce more than 240 000 achenes per growing season, mostly viable, which are capable of dissemination to 30 m in 16 km/h wind (Muenscher, 1935 ; Dauer et al., 2006 ; Nandula et al., 2006 ).

Interspecific hybridization may therefore affect the extinction rates of weeds in the environment and serve as a mechanism for the dissemination of transgenic and herbicide resistance genes (Owen and Zelaya, 2005 ). Farmers should consider the potential for hybridization between weeds when developing programs aimed at managing herbicide-resistant weeds. Production systems that depend on a single herbicidal chemistry for weed control should be combined with alternative management tactics, thus mitigating the evolution of herbicide resistance and maintaining the sustainability of current agroecosystems. Further research is needed to monitor native gene flow levels between weeds and ascertain the potential for dissemination of herbicide resistance through introgressive hybridization.

FOOTNOTES

1 The authors thank J. Wendel for his valuable discussions on evolutionary biology; M. Graham, P. Tranel, J. Gressel, and R. Hartzler for their critical review of earlier versions of this manuscript; D. Lewis at the Ada Hayden Herbarium (ISC) for mounting the specimens submitted to ISC, the Botanical Research Institute of Texas (BRIT), and the New York Botanical Garden (NY); and P. Knosby, J. Ruhland, and R. van der Laat for assistance with greenhouse endeavors. Back

4 Author for correspondence (e-mail: iazelaya{at}iastate.edu ; current address: Syngenta Ltd., Weed Control Research, Jealott's Hill International Research Centre, Bracknell, Berkshire, RG42 6EY, UK Back

LITERATURE CITED

Abbott R. J.. 1992. Plant invasions, interspecific hybridization and the evolution of new plant taxa. Trends in Ecology and Evolution 7: 401-405..[CrossRef]

Anamthawat-Jónsson K.. 2001. Molecular cytogenetics of introgressive hybridization in plants. Methods in Cell Science 23: 139-148..[Medline]

Anonymous. 2004. Crop protection reference. C & P Press, New York, New York, USA..

AOSA. 2003. Rules for testing seeds. Association of Official Seed Analysts, Las Cruces, New Mexico, USA..

Balfourier F. Charmet G. Ravel C.. 1998. Genetic differentiation within and between natural populations of perennial and annual ryegrass (Lolium perenne and L. rigidum). Heredity 81: 100-110..[CrossRef][Web of Science]

Barton N. H.. 2001. The role of hybridization in evolution. Molecular Ecology 10: 551-568..[CrossRef][Medline]

Bradshaw L. D. Padgette S. R. Kimball S. L. Wells B. H.. 1997. Perspectives on glyphosate resistance. Weed Technology 11: 189-198..[Web of Science]

Brown S. M. Whitwell T.. 1988. Influence of tillage on horseweed, Conyza canadensis. Weed Technology 2: 269-270..[Web of Science]

Cavan G. Biss P. Moss S. R.. 1998. Herbicide resistance and gene flow in wild-oats (Avena fatua and Avena sterilis ssp. ludoviciana). Annals of Applied Biology 133: 207-217..[CrossRef][Web of Science]

Chambers J. M. Cleveland W. S. Kleiner B. Tukey P. A.. 1983. Graphical methods for data analysis. Duxbury, Boston, Massachusetts, USA..

Comai L.. 2000. Genetic and epigenetic interactions in allopolyploid plants. Plant Molecular Biology 43: 387-399..[CrossRef][Web of Science][Medline]

Conover W. J.. 1999. Practical nonparametric statistics. John Wiley, New York, New York, USA..

Cronquist A.. 1980. Vascular flora of the southeastern United States, vol. I. Asteraceae. University of North Carolina Press, Chapel Hill, North Carolina, USA..

Darbyshire S. J.. 1990. Conyza ramosissima, another new weed in the Ottawa District. Trail & Landscape 24: 94-98..

Dauer J. T. Mortensen D. A. Humston R.. 2006. Controlled experiments to predict horseweed (Conyza canadensis) dispersal distances. Weed Science 54: 484-489..[Web of Science]

Dinelli G. Marotti I. Bonetti A. Minelli M. Catizone P. Barnes J.. 2006. Physiological and molecular insight on the mechanisms of resistance to glyphosate in Conyza canadensis (L.) Cronq. biotypes. Pesticide Biochemistry and Physiology 86: 30-41..[CrossRef][Web of Science]

Ellstrand N. C.. 2003. Current knowledge of gene flow in plants: implications for transgene flow. Philosophical Transactions of the Royal Society of London, B, Biological Sciences 358: 1163-1170..[CrossRef]

Feng P. C. C. Tran M. Chiu T. Sammons R. D. Heck G. R. CaJacob C. A.. 2004. Investigations into glyphosate-resistant horseweed (Conyza canadensis): retention, uptake, translocation, and metabolism. Weed Science 52: 498-505..[CrossRef][Web of Science]

Franssen A. S. Skinner D. Z. Al-Khatib K. Horak M. J. Kulakow P. A.. 2001. Interspecific hybridization and gene flow of ALS resistance in Amaranthus species. Weed Science 49: 598-606..[CrossRef][Web of Science]

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

Guadagnuolo R. Clegg J. Ellstrand N. C.. 2006. Relative fitness of transgenic vs. non-transgenic maize x teosinte hybrids: a field evaluation. Ecological Applications 16: 1967-1974..[CrossRef][Web of Science][Medline]

Harring T. Streibig J. C. Husted S.. 1998. Accumulation of shikimic acid: a technique for screening glyphosate efficacy. Journal of Agricultural and Food Chemistry 46: 4406-4412..[CrossRef][Web of Science]

Heap I.. 2006. The international survey of herbicide resistant weeds [online]. Website http://www.weedscience.com/ [accessed 15 December 2006]..

Holländer-Czytko H. Amrhein N.. 1983. Subcellular compartmentation of shikimic acid and phenylalanine in buckwheat cell suspension cultures grown in the presence of shikimate pathway inhibitors. Plant Science Letters 29: 89-96..[CrossRef][Web of Science]

Jasieniuk M.. 1985. Constraints on the evolution of glyphosate resistance in weeds. Resistant Pest Management Newsletter 7: 31-32..

Knobloch I. W.. 1972. Intergeneric hybridization in flowering plants. Taxon 21: 97-103..[Medline]

Koger C. H. Reddy K. N.. 2005. Role of absorption and translocation in the mechanism of glyphosate resistance in horseweed (Conyza canadensis). Weed Science 53: 84-89..[CrossRef][Web of Science]

Kwon Y. W. Kim D. S.. 2001. Herbicide-resistant genetically-modified crop: its risks with an emphasis on gene flow. Weed Biology and Management 1: 42-52..[CrossRef]

Laínz M.. 2001. Conyza blakei (Cabrera) Cabrera (Compositae), ¿especie peninsular?. Anales Jardín Botánico de Madrid 59: 352-353..

Légère A.. 2005. Risks and consequences of gene flow from herbicide-resistant crops: canola (Brassica napus L.) as a case study. Pest Management Science 61: 292-300..[CrossRef][Web of Science][Medline]

Mardia K. V.. 1970. Measures of multivariate skewness and kurtosis with applications. Biometrika 57: 519-530..[Abstract/Free Full Text]

McClintock D. Marshall J. B.. 1988. On Conyza sumatrensis (Retz) E. Walker and certain hybrids in the genus. Watsonia 17: 172-173..

Michaux A.. 1803. Flora Boreali-Americana, vol. II. Parisiis et Argentorati, Apud fratres Levrault, Paris, France..

Moore R. P.. 1985. Handbook on tetrazolium testing. International Seed Testing Association, Zürich, Switzerland..

Muenscher W. C.. 1935. Weeds. Macmillan, New York, New York, USA..

Nandula V. K. Eubank T. W. Poston D. H. Koger C. H. Reddy K. N.. 2006. Factors affecting germination of horseweed (Conyza canadensis). Weed Science 54: 898-902..[CrossRef][Web of Science]

Natali L. Giordani T. Polizzi E. Pugliesi C. Fambrini M. Cavallini A.. 1998. Genomic alterations in the interspecific hybrid Helianthus annuus x Helianthus tuberosus. Theoretical and Applied Genetics 97: 1240-1247..[CrossRef][Web of Science]

Nesom G. L.. 1990. Further definition of Conyza (Asteraceae: Astereae). Phytologia 68: 229-233..

Noyes R. D.. 2000. Biogeographical and evolutionary insights on Erigeron and allies (Asteraceae) from ITS sequence data. Plant Systematics and Evolution 220: 93-114..[CrossRef][Web of Science]

Owen M. D. K. Zelaya I. A.. 2005. Herbicide-resistant crops and weed resistance to herbicides. Pest Management Science 61: 301-311..[CrossRef][Web of Science][Medline]

Pruski J. F. Sancho G.. 2006. Conyza sumatrensis var. leiotheca (Compositae: Astereae), a new combination for a common neotropical weed. Novon 16: 96-101..[CrossRef][Web of Science]

Reichman J. R. Watrud L. S. Lee E. H. Burdick C. A. Bollman M. A. Storm M. J. King G. A. Mallory-Smith C.. 2006. Establishment of transgenic herbicide-resistant creeping bentgrass (Agrostis stolonifera L.) in nonagronomic habitats. Molecular Ecology 15: 4243-4255..[CrossRef][Medline]

Rieseberg L. H.. 1997. Hybrid origins of plant species. Annual Review of Ecology and Systematics 28: 359-389..[CrossRef][Web of Science]

Rieseberg L. H. Archer M. A. Wayne R. K.. 1999. Transgressive segregation, adaptation and speciation. Heredity 83: 363-372..[CrossRef][Web of Science][Medline]

Rieseberg L. H. Baird S. J. E. Gardner K. A.. 2000. Hybridization, introgression, and linkage evolution. Plant Molecular Biology 42: 205-224..[CrossRef][Web of Science][Medline]

SAS. 2000. SAS procedures guide, version 8. SAS Institute, Cary, North Carolina, USA..

Schardl C. L. Craven K. D.. 2003. Interspecific hybridization in plant-associated fungi and oomycetes: a review. Molecular Ecology 12: 2861-2873..[CrossRef][Medline]

Shapiro S. S. Wilk M. B.. 1965. An analysis of variance test for normality. Biometrika 52: 591-611..[Free Full Text]

Shontz J. P. Oosting H. J.. 1970. Factors affecting interaction and distribution of Haplopappus divaricatus and Conyza canadensis in North Carolina old fields. Ecology 51: 780-793..[CrossRef][Web of Science]

Sída O.. 2003. Conyza triloba, new to Europe, and Conyza bonariensis, new to the Czech Republic. Preslia 75: 249-254..[Web of Science]

Stace C. A.. 1975. Hybridization and the flora of the British Isles. Academic Press, London, UK..

Steinrücken H. C. Amrhein N.. 1980. The herbicide glyphosate is a potent inhibitor of 5-enolpyruvyl-shikimic acid-3-phosphate synthase. Biochemical and Biophysical Research Communications 94: 1207-1212..[CrossRef][Web of Science][Medline]

Steyermark J. A.. 1963. Flora of Missouri. Iowa State University Press, Ames, Iowa, USA..

Tabashnik B. E.. 1991. Determining the mode of inheritance of pesticide resistance with backcross experiments. Journal of Economic Entomology 84: 703-712..[Web of Science][Medline]

Teal P. E. A. Oostendorp A.. 1995. Effect of interspecific hybridization between Heliothis virescens and H. subflexa (Lepidoptera: Noctuidae) on sex pheromone production by females. Journal of Insect Physiology 41: 519-525..[CrossRef][Web of Science]

Thébaud C. Abbott R. J.. 1995. Characterization of invasive Conyza species (Asteraceae) in Europe: quantitative trait and isoenzyme analysis. American Journal of Botany 82: 360-368..[CrossRef][Web of Science]

Thébaud C. Finzi A. C. Affre L. Debussche M. Escarre J.. 1996. Assessing why two introduced Conyza differ in their ability to invade Mediterranean old fields. Ecology 77: 791-804..[CrossRef][Web of Science]

Thompson B.. 1984. Canonical correlation analysis: uses and interpretation. Sage Publications, Thousand Oaks, California, USA..

Tonsor S. J.. 1990. Spatial patterns of differentiation for gene flow in Plantago lanceolata. Evolution 44: 1373-1378..[CrossRef][Web of Science]

Tranel P. J. Wassom J. J. Jeschke M. R. Rayburn A. L.. 2002. Transmission of herbicide resistance from a monoecious to a dioecious weedy Amaranthus species. Theoretical and Applied Genetics 105: 674-679..[CrossRef][Web of Science][Medline]

Trucco F. Jeschke M. R. Rayburn A. L. Tranel P. J.. 2005. Promiscuity in weedy amaranths: high frequency of female tall waterhemp (Amaranthus tuberculatus) x smooth pigweed (A. hybridus) hybridization under field conditions. Weed Science 53: 46-54..[CrossRef][Web of Science]

VanGessel M. J.. 2001. Glyphosate-resistant horseweed from Delaware. Weed Science 49: 703-705..[CrossRef][Web of Science]

Verloove F. Boullet V.. 2001. Conyza bonariensis en Conyza sumatrensis: recent ingeburgerd in België?. Dumortiera 77: 2-8..

Vincent G. Cappadocia M.. 1987. Interspecific hybridization between common ragweed (Ambrosia artemisiifolia) and giant ragweed (A. trifida). Weed Science 35: 633-636..[Web of Science]

Voss E. G.. 1996. Michigan flora, part III, Dicots concluded. Cranbrook Institute of Science, Bloomfield Hills, Michigan, USA..

Watrud L. S. Lee E. H. Fairbrother A. Burdick C. Reichman J. R. Bollman M. Storm M. King G. Van de Water P. K.. 2004. Evidence for landscape-level, pollen-mediated gene flow from genetically modified creeping bentgrass with CP4 EPSPS as a marker. Proceedings of the National Academy of Sciences, USA 101: 14533-14538..[Abstract/Free Full Text]

Weaver S. E.. 2001. The biology of Canadian weeds. 115. Conyza canadensis. Canadian Journal of Plant Science 81: 867-875..[Web of Science]

Wetzel D. K. Horak M. J. Skinner D. Z. Kulakow P. A.. 1999. Transferal of herbicide resistance traits from Amaranthus palmeri to Amaranthus rudis. Weed Science 47: 538-543..[Web of Science]

Wilson P.. 1992. On inferring hybridity from morphological intermediacy. Taxon 41: 11-23..[CrossRef][Web of Science]

Wolf D. E. Takebayashi N. Rieseberg L. H.. 2001. Predicting the risk of extinction through hybridization. Conservation Biology 15: 1039-1053..[CrossRef][Web of Science]

Zelaya I. A. Owen M. D. K.. 2006. Interspecific hybridization in Asteraceae: case studies in Ambrosia, Conyza and Helianthus. Proceedings of the Weed Science Society of America Meeting 46: 90 (abstract)..

Zelaya I. A. Owen M. D. K. VanGessel M. J.. 2004. Inheritance of evolved glyphosate resistance in Conyza canadensis (L.) Cronq. Theoretical and Applied Genetics 110: 58-70..[CrossRef][Web of Science][Medline]


Add to CiteULike CiteULike   Add to Complore Complore   Add to Connotea Connotea   Add to Del.icio.us Del.icio.us   Add to Digg Digg   Add to Facebook Facebook   Add to Reddit Reddit   Add to Technorati Technorati   Add to Twitter Twitter    What's this?



This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Supplementary Data
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 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 Web of Science (2)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Zelaya, I. A.
Right arrow Articles by VanGessel, M. J.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Zelaya, I. A.
Right arrow Articles by VanGessel, M. J.
Agricola
Right arrow Articles by Zelaya, I. A.
Right arrow Articles by VanGessel, M. J.
Social Bookmarking
 Add to CiteULike   Add to Complore   Add to Connotea   Add to Del.icio.us   Add to Digg   Add to Facebook   Add to Reddit   Add to Technorati   Add to Twitter  
What's this?


HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS