|
|
||||||||
Population Biology |
2Departamento de Ecología, Genética y Evolución, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, 1428 Buenos Aires, Argentina; 3Instituto de Botánica Darwinion, Labardén 200, 1642 San Isidro, Buenos Aires, Argentina
Received for publication July 16, 2002. Accepted for publication August 5, 2003.
| ABSTRACT |
|---|
|
|
|---|
Key Words: Acacia Argentina Fabaceae genetic variability isozymes mating system parameters population structure
| INTRODUCTION |
|---|
|
|
|---|
In Argentina the genus is represented by 21 species. Two of them, A. aroma Hook. et Arn. and A. macracantha Willd., are widely distributed in northern and central Argentina (Cialdella, 1984
) because they can tolerate dryness and eroded soils. Both are economically important and promising multipurpose species for agroforesty programs (Dimitri and Biloni, 1973
; Allen and Allen, 1981
; Cialdella, 1984
; Martínez et al., 1996
). The only study on population genetics of Argentinean species of Acacia indicated that A. aroma and A. macracantha are closely related (Casiva et al., 2002
). While data on their genetic structure is just beginning to emerge, little information is avaliable on their mating systems. Mating system analysis in other species of Acacia have shown high levels of outcrossing, sometimes together with self-incompatibility systems, and some species are able to form a hybrid complex (Ali and Qaiser, 1980
; Bernhardt et al., 1984
; Kenrick and Knox, 1985
; Kenrick et al., 1986
; Moran et al., 1989b
; Sedgley et al., 1992
). Studies on the reproductive biology of A. caven from Chile and Argentina showed that this species is protogynous and has a self-incompatibility system (Peralta et al., 1992
; Baranelli et al., 1995
). Analyses of floral morphology and breeding systems in A. aroma and A. macracantha suggested that they are also outcrossers and self-incompatible, when bees, Trigona testaceicornis, Metadontia sp., and Allograpta sp. are considered as the main pollen vector (Zapata and Arroyo, 1978
; Bernhardt et al., 1984
) but mating system parameters have not been studied in these species.
An efficient method to evaluate mating system parameters in natural populations is based on the multilocus mixed mating model and the estimation procedure of Ritland and Jain (1981)
. The method requires grouping the population samples in family arrays and using neutral codominant loci that segregate independently. In the present case, we chose isozyme markers because they fulfill these requirements and give information about the genetic structure and distribution of genetic variability at population and species levels.
Here we tested the hypothesis that A. aroma and A. macracantha are outcrossers, estimated the distribution of genetic diversity, and evaluated the effect of sympatry on the genetic differentiation between populations of different species.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Three populations of A. aroma and two populations of A. macracantha were collected in Argentina (Fig. 1, Table 1). Thirty to 50 pods were collected from at least 10 trees separated from each other by more than 50 m following the method of Vilardi et al. (1988)
and Saidman and Vilardi (1993)
for legumes of genus Prosopis. All seeds (300500) from each shrub were placed in a single bag from which they were chosen at random for isozyme analysis of progeny arrays.
|
|
We used progeny arrays of open-pollinated families to infer genetic control of each enzyme system. Alleles within loci were indicated with superscript numbers with respect to their mobility relative to bromophenol blue, with "1" denoting the fastest allozyme. Bands that did not conform to Mendelian-segregation patterns were omitted from the analysis. For mating system analysis, we chose five polymorphic loci (SOD-1, SOD-3, SOD-4, 6PGD-2, and SKD-1) because they could be evaluated simultaneously in each individual (seed).
Data analysis
We used the program Biosys-1 (version 1.7) (Swofford and Selander, 1981
) to calculate allelic frequencies, mean number of alleles per locus (A), percentage of polymorphic loci (Pp, 0.95 criterion), and unbiased expected heterozygosity (He; Nei, 1978
) in the populations sampled at Tipas. The estimates of He, A, and Pp of the populations sampled at Campo Quijano and Huilla Catina had been previously reported in Casiva et al. (2002)
. Wright's (1951)
fixation index (FIS) was estimated for all five populations using the same program. We used this parameter to estimate the bias from Hardy-Weinberg expectations. Nei's (1978)
genetic distances among all populations were represented in a UPGMA phenogram (Sneath and Sokal, 1973
). Using the method of Sokal and Rohlf (1962)
, we estimated the cophenetic correlation coefficient. Bootstrap support values of each phenogram branch were obtained using PHYLIP (Felsenstein, 1993
). The phenetic relationships among populations were also represented by multidimensional scaling using the program STATISTICA (StatSoft, 1995
).
The mean expected heterozygosity was calculated at the population level (HP), species level (HC), and for the whole sample (HT). We performed a hierarchical analysis of population structure (Wright, 1978
), considering variation within populations (1 FST), among populations, within species (FST FCT) and between species (FCT). These statistics were estimated according to Nei (1987)
as FST = 1 HS/HT; FCT = 1 HC/HT; FSC = 1 HS/HC.
To test isolation by distance between populations and species, we compared the matrix of Nei's genetic distances with that of geographic distances among populations. The latter was estimated using the measurement tool of Encarta World Atlas (Microsoft 19951996
). The significance of that correlation was estimated with Mantel's (1976)
test using the ISOLDE program in the GENEPOP package (Raymond and Rousset, 1995
).
Using MLTR 2.2 (Ritland, 2002
), we calculated the following mating system parameters: multilocus outcrossing rate (tm), single locus outcrossing rate (ts), correlation of outcrossed paternity (rp), correlation of tm between progeny arrays (rt), and fixation index of maternal parents (FISM). Because the mixed mating model assumes independent segregation of alleles at different marker loci, we tested possible genotypic association among loci using GENEPOP (Raymond and Rousset, 1995
). This method performs a probability test (or Fisher exact test) for each contingency table using a Markov chain.
| RESULTS |
|---|
|
|
|---|
|
|
|
|
|
|
|
The single locus inbreeding coefficient of maternal parents (FISM) was in all cases lower than that estimated for their offspring (Tables 3, 7). In A. aroma (Tipas) FISM was negative, indicating some heterozygosity excess in the maternal plants. For the other four populations, FISM was near zero, suggesting no Hardy-Weinberg deviations in the populations.
In most populations, significant differences were observed for several loci between pollen and ovule allelic frequencies (confidence intervals of estimated frequencies do not overlap) (Table 8).
|
| DISCUSSION |
|---|
|
|
|---|
In the hierarchical analysis, we found no differences among species, and the most variability occurred within populations. These results agree with the predictions of Hamrick and Godt (1990)
that plant species, on average, maintain high levels of allozyme variation within populations. This same pattern was observed by Playford et al. (1993)
in populations of Acacia melanoxylon, for which most of the genetic diversity was found within populations. However, they found a high level of genetic patterning in the species and a strong differentiation among populations and geographic zones. Although in other studies on Acacia geographical patterns of differentiation of populations was found (Brain, 1986
; Joly et al., 1992
; Playford et al., 1993
), we found no geographic isolation and no genetic differentiation between species, suggesting a wide area of distribution of the allelic variants among all populations of these species. In all the studied populations of A. aroma and A. macracantha, polymorphic loci had the same allelic variants. Seed dispersal by livestock and/or human activity may explain the lack of correlation between genetic distances and geographic distances.
The FCT index indicated no significant differences between species. The high similarity was also revealed by Nei's genetic distances. The UPGMA phenogram and the MDS plot showed a similar trend, i.e., the populations are not grouped as expected according to morphological determination. The small genetic distance between these two taxa was also noticed previously by Casiva et al. (2002)
. In fact, the distances recorded between populations belonging to these species are as low as those corresponding to conspecific populations of other Acacia species (e.g., A. caven; Casiva et al., 2002
) or populations of some species of the mimosoid genus Prosopis that belong to a syngameon (see Saidman and Vilardi, 1987
).
The levels and distribution of genetic variability may also be related to mating system (Hamrick and Godt, 1990
). Mating system analysis showed high levels of outcrossing in all populations studied in our work. On average, only 2% of progeny would result from selfing. The small differences between tm and ts and the differences between pollen and ovule allelic frequencies are consistent with the lack of biparental inbreeding. The rp values suggest that individuals within progeny arrays are full rather than half sibs. This pattern could be a consequence of hierarchical mating, previously noted for Acacia species (Muona et al., 1991
). This is the result of the composite pollen grains, which are distributed in polyads. The number of pollen grains per polyad varies among Acacia species, with four, eight, 16, and 32 grains per polyad (Vassal, 1972
; Guinet and Vassal, 1978
; Cialdella, 1984
). A correlation between polyad grain number (pgn) and maximum pod seed number (mpsn) was suggested by Kenrick and Knox (1982)
, who observed that pgn
mpsn and proposed that polyads in Acacia could be a mechanism helping to ensure seed set following a single pollination event (a single polyad fertilizing all available ovules at a single ovary). However, A. macracantha might be an exception to the predictions of Kenrick and Knox (1982)
. In fact, although both A. aroma and A. macracantha have 16 pollen grains per polyad, A. aroma presents 616 seeds per fruit and A. macracantha pods may have as many as 20 seeds (Cialdella, 1984
). Therefore, in A. macracantha two pollination events might be involved in some single pods.
Nonrandom pollination between pods may also occur because flowers are grouped into clusters, which often form parts of a complex raceme. Because all flowers within a cluster often open simultaneously, the probability that the same insect pollinates several flowers with polyads from the same paternal tree may be high (Muona et al., 1991
). Because we sampled randomly seeds from different pods, the high proportion of full sibs in our sample suggests that pollinating insects tend to visit many flowers in the same shrub before moving to a different plant.
Based on the results from our study of mating system parameters, we can reject selfing as a mechanism of isolation between these sympatric entities. The high variability observed by Casiva et al. (2002)
and the present estimates of He and tm indicate that these species are outcrossers.
The actual relationships between A. aroma and A. macracantha are under debate. The morphological criteria differ between authors. According to Cialdella (1984
, 1997
) and Casiva et al. (2002)
, the most important characters for morphological differentiation are the shape and length of the thorny stipules. On the other hand, Ebinger et al. (2000)
hold that A. aroma and A. macracantha could be differentiated by fruit and petiolar gland morphology. However, in Argentinean material we observed variation in fruit morphology within individuals, the range of peduncle length in both species overlaps, and both species have sessile petiolar glands. The biochemical evidence is also intriguing. Lamarque et al. (2000)
found that seed chemical components (moisture, fat, protein, ash, and fatty acid profiles) permit separation of A. aroma from A. macracantha. On these grounds, they support recognition of two distinct species. However, the associations among species obtained from morphology by Cialdella (1984
, 1997)
and from isozymes and RAPDs by Casiva et al. (2002)
differ from those obtained by Lamarque et al. (2000)
.
The high genetic similarity between A. aroma and A. macracantha is not expected for different biological species. The hierarchical analysis indicated absence of differentiation between them. The low estimates of genetic distances are compatible with very short divergence between these entities (see Nei, 1987
). Indirect estimations of gene flow suggested that these species are able to exchange genes (Casiva et al., 2001a
, 2002
). The outcrossing system is not compatible with sympatric speciation if the species undergo significant gene flow.
Casiva et al. (2002)
advanced three possible hypothesis: (1) A. aroma and A. macracantha may be varieties of a single species, (2) they may be members of a single species where thorny stipule shape and length are determined by a few linked loci, or (3) they are different species that diverged very recently and have yet to accumulate significant genetic differences. The geographical range together with the present results on the mating system constitute evidence against the possibility that A. aroma and A. macracantha are varieties or subspecies unless they could be isolated by different pollen vectors. In Venezuelan populations of A. macracantha, Trigona testaceicornis, Metadontia sp., and Allograpta sp. are described as pollinators (Zapata and Arroyo, 1978
), but no information is available for A. aroma. Entomological studies in Argentinean populations showed no differences in the insect fauna hosted by both species (Di Iorio, in press
). According to Cialdella (1984)
and to our field observations, they are sympatric and their flowering and fruiting periods overlap. Because interspecific Nm estimates are high and the hierarchical analysis does not provide evidence of isolation between these entities, the present results rule out hypothesis 1. Because the morphological differences between these species refer to a single trait (thorny stipule), the condition resembles that of Prosopis glandulosa for which a polymorphism was described for the presence of the nectarium (Golubov et al., 1999
). Although hypothesis 3 cannot be ruled out, the lack of genetic divergence between A. aroma and A. macracantha favors the hypothesis that they may be a single polymorphic species.
| FOOTNOTES |
|---|
| LITERATURED CITED |
|---|
|
|
|---|
Allen O. N. E. K. Allen 1981 The Leguminosae. A source book of characteristics, uses and nodulation. University of Wisconsin, Madison, Wisconsin, USA
Baranelli J. L. A. A. Cocucci A. M. Anton 1995 Reproductive biology in Acacia caven (Mol.) Mol. (Leguminosae) in the central region of Argentina. Botanical Journal of the Linnean Society 119: 65-76[CrossRef]
Bernhardt P. J. Kenrick R. B. Knox 1984 Pollination biology and the breeding system of Acacia retinodes (Leguminosae: Mimosoideae). Annals of the Missouri Botanical Garden 71: 17-29[CrossRef][ISI]
Brain P. 1986 Leaf peroxidase types in Acacia karroo. Geographical distribution and influence of the environment. South African Journal of Botany 26: 365-379
Butcher P. A. G. F. Moran H. D. Perkins 1998 RFLP diversity in the nuclear genome of Acacia mangium. Heredity 81: 205-213[CrossRef][ISI]
Casiva P. 2001 Utilización de marcadores moleculares y morfológicos para el estudio de la diferenciación y variabilidad en poblaciones naturales argentinas del género Acacia (Leguminosae). Tesis de Licenciatura, University of Buenos Aires, Buenos Aires, Argentina
Casiva P. B. O. Saidman A. M. Cialdella J. C. Vilardi 2001a Escasa diferenciación interespecífica entre Acacia aroma y A. macracantha (Fabaceae). Journal of Basic and Applied Genetics 2: 142-143
Casiva P. B. Saidman J. Vilardi A. Cialdella 2001b Caracterización morfométrica en cuatro especies argentinas de Acacia (Morphometric characterization in four argentinean species of Acacia). Boletin de la Sociedad Argentina de Botánica 36: (Supplement 2001) 11.
Casiva P. V. B. O. Saidman J. C. Vilardi A. M. Cialdella 2002 First comparative phenetic studies of Argentinean species of Acacia (Fabaceae), using morphometric, isozymal, and RAPD approaches. American Journal of Botany 89: 843-853
Cialdella A. M. 1984 El género Acacia (Leguminosae) en la Argentina. Darwiniana 25: 59-111
Cialdella A. M. 1997 Acacia. In A. T. Hunziker [ed.], Flora fanerogámica Argentina, fascicle 25, 321. Proflora CONICET, Córdoba, Argentina
Di Iorio O. In press Torneutini (Coleoptera: Cerambycidae) of Argentina. Part 1b. New records, larval and adult biologies, and comparisons with Prioninae and Trachyderini. Giornale Italiani di Entomologia
Dimitri M. J. J. S. Biloni 1973 El libro del árbol, I. Celulosa Argentina S.A. Buenos Aires, Argentina
Ebinger J. E. D. Seigler H. D. Clarke 2000 Taxonomic revision of South American species of the genus Acacia (Fabaceae: Mimosoideae). Systematic Botany 25: 588-617[CrossRef][ISI]
Fagg C. W. R. D. Barnes C. Marunda 1997 African Acacia Trials Network: a seed collection of six species for provenance-progeny test: held at the Oxford Forestry Institute. FAOFOREST Genetic Resources 25, Oxford, UK
Felsenstein J. 1993 PHYLIP (Phylogenetic inference package), version 3.5c. Department of Genetics, University of Washington, Seattle, Washington, USA
Golubov J. L. E. Eguiarte M. C. Mandujano J. López-Portillo C. Montana 1999 Why be a honeyless honey mesquite? Reproduction and mating system of nectarful and nectarless individuals. American Journal of Botany 86: 955-963
Guinet P. J. Vassal 1978 Hypotheses on the differentiation of the major groups in the genus Acacia (Leguminosae). Kew Bulletin 32: 509-527[CrossRef]
Hamrick J. C. J. W. Godt 1990 Allozyme diversity in plant species. In H. D. Brown, M. T. Clegg, A. L. Kahler, and B. S. Weir [eds.], Plant population genetics, breeding, and genetic resources, 4365, Sinauer, Sunderland, Massachusetts, USA
IUPAC-IUB (International Union of Pure and Applied Chemistry and the International Union of Biochemistry) 1979 Enzyme nomenclature: recommendations of the nomenclature committee of the International Union of Biochemistry. Academic Press, New York, New York, USA
Joly H. I. M. Zeh-Nlo P. Danthu C. Aygalent 1992 Population genetics of an African acacia, Acacia albida I. Genetic diversity of populations from West Africa. Australian Journal of Botany 40: 59-73[CrossRef][ISI]
Kenrick J. V. Kaul E. G. Williams 1986 Self-incompatibility in Acacia retinodes: site of pollen-tube arrest in the nucellus. Planta 169: 245-250[CrossRef][ISI]
Kenrick J. R. B. Knox 1982 Function of the polyad in reproduction of Acacia. Annals of Botany 50: 721-727
Kenrick J. R. B. Knox 1985 Self-incompatibility in the nitrogen-fixing tree, Acacia retinodes: quantitative cytology of pollen tube growth. Theoretical and Applied Genetics 69: 481-488[CrossRef][ISI]
Lamarque A. L. R. H. Fortunato D. M. Maestri C. Guzmán 2000 Seed components and taxonomy of some Acacia species. Biochemical Systematic and Ecology 28: 53-60
Mantel N. 1976 The detection of disease clustering and a generalized regression approach. Cancer Research 27: 209-220
Martínez M. G. León de Pinto C. Rivas E. Ocando 1996 Chemical and spectroscopic studies of the gum polysaccharide from Acacia macracantha. Carbohydrate Polymers 29: 247-252[CrossRef][ISI]
McGranahan M. J. C. Bell G. F. Moran M. Slee 1997 High genetic divergence between geographic regions in the highly out crossing species Acacia aulacocarpa (Cunn. ex Benth). Forest Genetics 4: 1-13
Microsoft. 19951996 Microsoft Encarta World Atlas (Spanish Version). Microsoft Corporation. Redmond, Washington, USA
Moran G. F. 1992 Patterns of genetic diversity in Australian tree species. New Forest 6: 49-66
Moran G. F. O. Muona J. C. Bell 1989a Acacia mangium: a tropical forest tree of the coastal lowlands with low genetic diversity. Evolution 43: 231-235[CrossRef][ISI]
Moran G. F. O. Muona J. C. Bell 1989b Breeding systems and genetic diversity in Acacia auriculiformis and Acacia crassicarpa. Biotropica 21: 250-256[CrossRef][ISI]
Muona O. G. F. Moran J. C. Bell 1991 Hierarchical patterns of correlated mating in Acacia melanoxylon. Genetics 127: 619-626[Abstract]
Nei M. 1978 Estimation of average heterozygosity and genetic distance from a small number of individuals. Genetics 83: 583-590
Nei M. 1987 Molecular evolutionary genetics. Columbia University Press, New York, New York, USA
Peralta I. J. G. Rodríguez M. T. Kalin Arroyo 1992 Breeding system and aspects of pollination in Acacia caven (Mol.) Mol. (Leguminosae: Mimosoideae) in the mediterranean-type climate zone of central Chile. Botanische Jahrbücher für Systematik 114: 297-314
Playford J. R. Appels B. R. Baum 1992 The 5S DNA units of Acacia species (Mimosaceae). Plant Systematics and Evolution 183: 235-247[CrossRef][ISI]
Playford J. J. C. Bell G. F. Moran 1993 A major disjunction in genetic diversity over the geographic range of Acacia melanoxylon R. Br. Australian Journal of Botany 41: 355-368[CrossRef]
Raymond M. F. Rousset 1995 GENEPOP version 1.2: population genetic software for exact tests and ecumenicism. Journal of Heredity 86: 248-249
Ritland K. 2002 Extensions of models for the estimation of mating systems using n independent loci. Heredity 88: 221-228[CrossRef][ISI][Medline]
Ritland K. S. Jain 1981 A model for the estimation of outcrossing rate and gene frequencies using n independent loci. Heredity 47: 35-52[ISI]
Ross J. H. 1981 An analysis of the African Acacia species: their distribution, possible origins and relationships. Bothalia 13: 389-413
Saidman B. O. 1985 Estudio de la variación alozímica en el género Prosopis. Ph.D. thesis, University of Buenos Aires, Buenos Aires, Argentina
Saidman B. O. J. C. Vilardi 1987 Analysis of the genetic similarities among seven species of Prosopis (Leguminosae: Mimosoideae). Theoretical and Applied Genetics 75: 109-116[ISI]
Saidman B. O. J. C. Vilardi 1993 Genetic variability and germplasm conservation in the genus Prosopis. In S. Puri [ed.], Nursery technology of forest tree species of arid and semiarid regions, 187198. Winrock-Oxford & IBH Publishing, New Delhi, India
Sedgley M. J. Harbard R. M. M. Smith R. Wickneswari A. R. Griffin 1992 Reproductive biology and interspecific hybridization of Acacia mangium and A. auriculiformis A. Cunn. ex Benth. (Leguminosae: Mimosoideae). Australian Journal of Botany 40: 37-48[CrossRef]
Sneath P. H. A. R. R. Sokal 1973 Numerical taxonomy: the principles and practice of numerical classification. Freeman, San Francisco, California, USA
Sokal R. R. F. J. Rohlf 1962 The comparisons of dendograms by objective methods. Taxon 11: 33-43
StatSoft. 1995 STATISTICA for Windows. StatSoft, Tulsa, Oklahoma, USA
Swofford D. L. R. B. Selander 1981 BIOSYS-1: a FORTRAN program for the comprehensive analysis of electrophoretic data in population genetics and systematics. Journal of Heredity 72: 281-283
Vassal J. 1972 Apport des recherches ontogéniques et séminologiques a l=étude morphologique, taxonomique et phylogénique des genre Acacia. Bulletin de la Societé d=Histoire Naturelle de Toulouse 108: 1-115
Vilardi J. C. B. O. Saidman R. A. Palacios 1988 Muestreo según variabilidad. In Prosopis en Argentina, 119124. Documento preliminar elaborado para el I Taller Internacional sobre Recurso Genético y Conservación de Germoplasma en Prosopis. Facultad de Ciencias Agrarias (U. N. Córdoba), FAOPIRB, Córdoba, Argentina
Wright S. 1951 The genetical structure of populations. Annals of Eugenics 15: 323-354[ISI]
Wright S. 1978 Evolution and the genetics of populations, vol. 4, Variability within and among natural populations. University of Chicago Press, Chicago, Illinois, USA
Zapata T. M. Arroyo 1978 Plant reproductive ecology of a secondary deciduous tropical forest in Venezuela. Biotropica 10: 221-230
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |