|
|
||||||||
a Department of Botany and Zoology, The University of Hong Kong, Pokfulam Road, Hong Kong; and b School of Applied Science, Griffith University Gold Coast, Gold Coast MC PMB 50, Queensland 9726, Australia
| ABSTRACT |
|---|
|
|
|---|
Key Words: Kandelia mangrove outcrossing rate pollination biology population genetic structure Rhizophoraceae
| INTRODUCTION |
|---|
|
|
|---|
One of the important attributes of many mangroves is the production of viviparous offspring. Vivipary in flowering plants is defined as the precocious and continuous growth of the offspring when still attached to the maternal parent. In many viviparous mangrove species, the embryo penetrates through the fruit pericarp and grows to a considerable size before dispersal, such as in Kandelia candel and many other species in the Rhizophoraceae. Vivipary could be an adaptive feature for the plant to overcome the harsh tidal environment for seedling establishment in the parental sites, but it is not considered adaptive for dispersal either in time or space (Elmqvist and Cox, 1996). The rare successful long-distance dispersal event and rare opportunity for seedling establishment would cause founder effects in the resultant populations. The subsequent migration rate through dispersal of the viviparous seedlings at a microgeographical scale would determine the extent of genetic differentiation among local populations. However, the effect of this unique dispersal mechanism on the level of genetic variation within species and its distribution among populations have not previously been reported for any viviparous mangrove species.
In addition to the rate of migration or gene flow through seed or seedling dispersal, the pollination mechanism and mating system of a species also play an important role in determining population genetic structure. Although the mating systems of mangrove species may be inferred from their reproductive features, genetic marker analysis can provide more accurate information on the actual level of outcrossing in different populations. A few previous studies have used the offspring segregation ratios of chlorophyll-deficient heterozygotes to estimate outcrossing rates in natural populations of Rhizophora mangle (Lowenfeld and Klekowski, 1992; Klekowski, Lowenfeld, and Hepler, 1994) and Kandelia candel (Chen, Lin, and Lin, 1996). However, allozyme marker-based multilocus estimates of mating-system parameters have apparently not been obtained for any mangrove species.
The objective of this study was to investigate population genetic attributes of a typical species of mangrove, Kandelia candel (L.) Druce (Rhizophoraceae), in relation to its reproductive biology and viviparous seedling dispersal. Kandelia candel is one of the dominant species of mangroves along the Chinese coast (Li and Lee, 1997). Although considered to be of limited abundance throughout its range (Tomlinson, 1986), K. candel demonstrates reproductive habits typical of the Rhizophoraceae. Information on the mating system parameters will help explain the pattern of genetic diversity in the species. Knowledge of population genetic structure of this important mangrove species will also contribute to its conservation and management. In Hong Kong, mangrove stands have been greatly reduced to ~15% of their original distribution over the past few decades, mainly due to reclamation and urban development (Yipp, Hau, and Walthew, 1995). Up to the present, only limited ecological surveys on these mangroves have been conducted (e.g., Tam and Wong, 1997). No genetic information is available for incorporation in developing conservation strategies. This study provides needed genetic data for guiding conservation and management practice that would maximize genetic diversity in the species.
| MATERIALS AND METHODS |
|---|
|
|
|---|
Pollination ecology and bagging experiment
For studying population genetic structure, 13 populations of K. candel throughout Hong Kong were sampled in May and June 1997. Flower buds or the youngest leaf buds were collected from ~30 individuals per population site (see Fig. 1 for population names and locations). For estimating outcrossing rates, two of the population sites, Site A and Maipo, were sampled. These sites were chosen for having abundant propagules still attached to the trees at the time of collection. An average of eight propagules per tree and ~30 trees per site were sampled to provide materials for progeny arrays used in mating system analysis. The propagules collected from the same tree were kept separate on a family basis and half-immersed in water with a salinity similar to that at the sites of collection, until the completion of isozyme electrophoresis.
|
Isozyme electrophoresis
Flower buds or the youngest leaf buds were collected from the field and kept cool and moist until electrophoresis, normally conducted on the day following the collection. The fresh materials were ground with an extraction buffer as described in Sun and Ganders (1990). Two buffer systems were used to survey 14 enzyme systems in 12.5% starch gels. The "H" buffer system (histidine-tris-citrate system, pH 7.0; Sun and Corke, 1992) was used to resolve isocitrate dehydrogenase (IDH; E.C. [Enzyme Commission number] 1.1.1.42), malate dehydrogenase (MDH; E.C. 1.1.1.40), malic enzyme (ME; 1.1.1.40), phosphoglucomutase (PGM; E.C. 2.7.5.1), 6-phosphogluconate dehydrogenase (6PGD; E.C. 1.1.1.44), and shikimate dehydrogenase (SKDH; E.C. 1.1.1.25). The "L" buffer system (pH 7.8; Shields, Orton, and Stuber, 1983) was used to resolve aspartate aminotransferase (AAT; E.C. 2.6.1.1), aconitase (ACO; E.C. 4.2.1.3), esterase (EST; E.C. 3.1.1.1), leucine aminopeptidase (LAP; E.C. 3.4.11.1), phosphoglucoisomerase (PGI; E.C. 5.3.1.9), and triose-phosphate isomerase (TPI; E. C. 5.3.1.1). Fluorescent esterase (F-EST; E.C. 3.1.1.1) and ß-glucosidase (GLU; E.C.3.2.1.21) were resolved by either of the buffer systems. The protocols of Wendel and Weeden (1989) were used for enzyme activity staining. Genetic interpretation of band patterns followed standard principles (Weeden and Wendel, 1989; Wendel and Weeden, 1989).
Data analysis
Measures of within-population genetic variation include polymorphism, allelic diversity, and heterozygosity. The percentage of polymorphic loci (P ) was calculated based on the 99% criterion. The P values were not significantly different when the 95% criterion was used. Allelic diversity was measured as the number of alleles per locus (A). Gene diversity (H), the mean expected heterozygosity assuming Hardy-Weinberg equilibrium, was calculated as H = 1 - 1/n
j
ipij, where pij is the frequency of the ith allele at the jth locus, and n is the total number of isozyme loci surveyed. Heterozygosities observed (Ho) and expected (He) were calculated only for two commonly polymorphic loci, Skdh-1 and Tpi-2, in these populations.
The proportion of genetic variation distributed within (HS) vs. between populations (DST) was measured using Nei's (1973) gene diversity statistics. The coefficient of genetic differentiation among populations, GST, was used to estimate the level of gene flow, Nm (the number of migrants exchanged between local populations per generation), based on the relationship GST = 1/(4Nm+1) where GST is Nei's (1973) estimator of FST (Wright, 1951). The UPGMA (unweighted pair-group arithmetic averaging) dendrogram of Nei's (1972) genetic distance was constructed using the GDD computer program provided by K. Ritland.
For the estimation of outcrossing rates, two polymorphic allozyme loci, Skdh-1 and Tpi-2, were used as gene markers in Site A, and two more loci, Mdh-2 and Pgm-1, were sufficiently polymorphic as markers in Maipo. Single-locus Mendelian inheritance at these loci was verified by the conformity of progeny array patterns to Mendelian predictions, a useful method when controlled crosses are not performed (Kephart, 1990). Mating-system parameters were estimated using Ritland's MLTR program following Sun and Ritland (1998). The mating-system parameters estimated included: (a) the multilocus population outcrossing rate tm (Ritland and Jain, 1981); (b) the average of single-locus population outcrossing rate ts; (c) the inbreeding coefficient of maternal parents F; (d) the correlated selfing rate rs; and (e) the correlation of paternal genotypes between outcrossed sibs rp. After it was determined that pollen and ovule population gene frequencies did not significantly differ, these frequencies were pooled to increase statistical power to estimate mating system parameters. The Newton-Raphson method was used to find the maximum likelihood estimates. Standard errors of population estimates, including tm - ts, were calculated using 100 bootstraps, where the unit of resampling was the progeny array.
| RESULTS |
|---|
|
|
|---|
Of the 20 inflorescence branches bagged in late June, we found only two fruits in early September, one each from two bags. The sample estimate of the number of flowers per bagged branch is 18.6 ± 5.0, and thus the fruit set in the bagged inflorescences was only ~0.54%. A sample of 28 open-pollinated inflorescences on the same set of trees produced an average of 4.54 ± 2.0 fruits each, giving an estimate of 24.41% natural fruit set at the site. A common feature among the viviparous mangroves is that only one seedling per fruit emerges while still attached to the mother tree, and there are six ovules per flower, so the actual fertility of the bagged flowers is as low as 0.09%. These data suggest that automatic self-fertilization events are rare in the species, and natural fruit setting is almost completely pollinator dependent.
The mating system analysis showed that the inbreeding coefficient of maternal parents (F) was not significant in the populations (Table 1). Multilocus outcrossing rates ranged from 70 to 80%, and averaged 75 ± 7% for the two populations studied. The single-locus outcrossing estimates varied within populations, and the standard error (SE) estimates at each marker locus appear to be affected by the allele frequency distribution (Table 2). The marker loci having skewed allele frequencies, such as Mdh-2 and Pgm-1, tend to have much higher SE values than the markers having more even allele frequencies, such as Skdh-1 and Tpi-2. It is apparent that more accurate outcrossing estimates could be obtained from these more informative marker loci. Genetic substructuring, in terms of biparental inbreeding (detectable by the difference between multilocus and single-locus selfing rates), was not significant in either population. The correlation of paternity (i.e., probability that sibs shared the same father), rp, was significant at Site A, indicating a small "paternal mating pool" of ~23 plants in this population (Table 1). The "paternal mating pool" is calculated as the inverse of rp and is the number of pollen donors that gives rise to this correlation, assuming all paternal trees have equal mating probabilities and consecutive matings are independent (Ritland, 1989).
|
|
|
|
|
| DISCUSSION |
|---|
|
|
|---|
Mating system has traditionally been characterized by selfing or outcrossing rates, which influence patterns of population genetic structure. The pattern of paternity, measured by the parameters of the "correlated mating model" of Ritland (1989), is often missing in plant mating system studies. These parameters describe the genetic relatedness among progeny of a plant. One parameter is the correlation of outcrossed paternity within progeny arrays, rp (defined as the proportion of pairs of progeny that are full-sibs, as opposed to half-sibs; see Ritland, 1989). This correlation arises from repeated matings to the same near-neighbors or from the tendency of pollinator pollen loads to be derived from single plants. In partially selfing populations, correlations among progeny also arise from variation of selfing rates among plants. This is measured by a second parameter, the correlation of selfing within progeny arrays, rs (Ritland, 1989). In addition to selfing or outcrossing rates, correlations among progeny affect local patterns of genetic variation.
This study revealed a significant level of selfing in K. candel, with a mean multilocus selfing rate of 25.3% in the populations assayed. Based on the results of the bagging experiment, automatic self-pollination is obviously negligible, if any, in the species. However, as the plant produces many flowers concurrently at anthesis, and bees often visit several flowers consecutively on the same inflorescence or different inflorescences on the same plant, the potential for geitonogamous selfing is high. The estimated 25.3% selfing rate could be partly or mainly a result of geitonogamous selfing. As no significant differences were detected between tm and ts, biparental inbreeding is apparently not an important factor in the populations of K. candel. The occurrence of outbred parents, as shown by nonsignificant F and f values detected in these populations (Tables 1, 4), is in agreement with nonsignificant biparental inbreeding. Although we detected a significant heterozygote excess in one population and no significant departure from Hardy-Weinberg expectations in 12 other populations studied, the estimated selfing rates in the populations showed a significant departure from random mating. This suggests the presence of selection against inbred progeny, which can occur if genetic load is high, such as in cases of traditionally outcrossing species.
Based on a marker gene for chlorophyll deficiency, the estimated single-locus selfing rates ranged from 17.7 to 19.7% in two populations of Kandelia candel in Fujan (Chen, Lin, and Lin, 1996), not significantly lower than the average single-locus selfing rate (23.9%, based on Table 2) detected in the present study. In contrast, using the same marker gene for chlorophyll deficiency, the selfing rates were estimated to be 71.2 to 95.2% in populations of Rhizophora mangle (Lowenfeld and Klekowski, 1992; Klekowski, Lowenfeld, and Hepler, 1994), indicating that this mangrove species is a predominant selfer.
Sharing of paternity within a family, as measured by the correlations of outcrossed paternity, rp, could arise from multiple deposits of pollen from a single male parent, or from mating to a small number of near neighbors. In K. candel, the extent of correlated paternity between flowers of the same plant is significant at Site A (0.351), equivalent to random mating to an effective pool of 1/0.351 or about three neighbors. The large standard error of rp at Maipo made the estimate insignificant, but it can still be concluded that, on average, at least six neighboring trees contributed genes to the progeny of the maternal parent in the population. Differences in the degree of relatedness among progeny could contribute to the observed variation in genetic structure between local populations (Table 4).
Population genetic consequences
Nearly all species with true vivipary tend to inhabit shallow marine habitats, which are normally characterized by extraordinarily coarse-grained environments for seedling establishment (Elmqvist and Cox, 1996). Selection in such environments would favor vivipary since the seed germinates while still attached to the mother tree and the young seedling grows to a considerable length before falling. Thus, seedling survival and establishment in the parental sites would be easier to achieve. However, viviparous seedlings are considered ill-adapted for dispersal since they lack the protection and nutritional support from the maternal tissue. However, many of the viviparous seedlings are buoyant, capable of long-distance dispersal by ocean currents, such as in the case of K. candel, whose propagules can grow to 2047 cm long and 1.01.3 cm wide, depending on population locations and environmental conditions (Maxwell, 1995). Propagule size may reach 70 cm in Rhizophora mucronata, another member of the Rhizophoraceae. Colonization through long-distance dispersal by ocean current may frequently occur in these species, but likely involves a small number of founders. In addition, the probability of an offspring being dispersed to a site better than the parental patch is considered very low in the marine environment (Elmqvist and Cox, 1996). Thus, the rare success of a few founders in establishing a new population would lead to a genetic erosion. This probably accounts for the observed amount of genetic variation in K. candel in Hong Kong, which is much lower than the mean levels of genetic variation in species with mixed-mating systems (see review in Hamrick and Godt, 1996). Since very little information is available on other viviparous mangrove species regarding their mating systems and population genetic diversity, no conclusion can be made as to whether the pattern of genetic diversity in K. candel is representative of those existing in other viviparous mangrove species. A recent study of a nonviviparous mangrove species, Acanthus ilicifolius Linn. (Acanthaceae), revealed a very low level of within-population genetic variation but substantial genetic polymorphism between populations (Lakshmi et al., 1997). However, no information on its pollination biology or mating system is given in the paper. It is not clear whether the pattern of genetic variation in this species is primarily affected by its mating system or largely due to its being an opportunistic colonizer.
The high genetic identities among the populations of K. candel in Hong Kong suggest their recent coancestry and/or lack of polymorphism in the species. A recent coancestry and/or a high level of gene flow among the local populations could result in the lack of population genetic differentiation (GST = 0.064). However, at the two common marker loci, Skdh-1 and Tpi-2, microgeographic differentiation in allele frequencies is apparent, resulting in the clustering of populations into two major groups separating the west coast from the east coast. Thus, it can be concluded that seedling dispersal in K. candel is somehow restricted between coastlines in Hong Kong. Also, genetic differentiation occurs much faster at the isolated sites. For example, individuals of Site A are located in a shallow impoundment at the Maipo marshes long-established (>50 yr) by farmers for culturing shrimps. Natural dispersal of propagules by tidal waves or ocean current is apparently limited as can be seen by the accumulation of propagules under the canopy at the site. Lee (1989) estimated that <1% of the litter (including propagules) produced by K. candel inside the ponds was exported away from the mangrove stand. The level of migration or gene flow through viviparous seedlings in K. candel apparently played an important role in shaping population genetic structure.
Due to the recent construction and development of an airport, highways, and ports on Lantau Island, many mangroves have been destroyed and those remaining are scattered on all sides of the island (Tam and Wong, 1997). Of the three populations of K. candel sampled from this island, two (Tungchung and Tongfuk) clustered with the West Coast populations, whereas one (Tai O) clustered with the East Coast populations (Fig. 2). Historical factors may account for this pattern of genetic divergence. The sampled site in Tai O was located in an abandoned salt pan that was used to produce salts for several centuries (Tam and Wong, 1997). Some salt pans were later converted to fish ponds and seawalls were built to protect the ponds from strong ocean currents and erosion. Long isolated from the sea, natural gene flow through propagules between this site and others could be very restricted. Genetic drift in allele frequencies in this small population likely accounts for its clustering by chance with the East Coast populations. Due to rapid development in Tai O, very few mangroves remain at the site, and the small patches still present are fragmented and extremely polluted with sewage from neighboring restaurants and households. However, a recent ecological survey suggests that the abandoned salt pans can be made suitable for artificial replanting of mangroves (Tam and Wong, 1997). Enhanced gene flow through transplanting propagules from other populations to this site would also help to maintain the total gene diversity in the species.
Conservation implications in Hong Kong
Compared to the mangroves in the tropics, fewer mangrove species are present in Hong Kong, with smaller stature and limited abundance. Of the eight mangrove species occurring in Hong Kong, four produce viviparous seedlings. Among them, K. candel, Avicennia marina, and Aegiceras corniculatum are the main components in the seaward mangroves. In recent years, urban development in Hong Kong has exerted great pressure on mangroves. Mangroves have been cut and reclaimed at an alarming rate (Morton and Morton, 1983). The remaining mangroves can be found at ~45 locations along the coastlines of the mainland territories and the two major islands, Hong Kong Island and Lantau Island. K. candel is one of the three seaward species and constitutes the major "pioneer" species responsible for seaward extension of the mangrove communities. Research on mangroves on the Chinese coast has been limited to the documentation of species composition and their ecology (see Li and Lee, 1997 for a review). Research on population genetic diversity of mangrove species is still in its infancy worldwide. Such knowledge is important to the management of the mangrove resources. Knowledge of genetic composition of populations and the degree of interpopulational divergence can help managers to make informed decisions regarding how to protect the existing genetic resources. As the difference among populations of K. candel in Hong Kong is mostly in allele frequencies rather than in gene composition, loss of populations at certain locations may not cause immediate loss in genetic diversity but more damage may occur in terms of long-term genetic consequences due to reduced number of populations and smaller population size. Replacement of the lost populations at certain locations can be achieved by transplanting a large number of propagules from long-established, large populations.
| FOOTNOTES |
|---|
| REFERENCES |
|---|
|
|
|---|
Chen, X-Y., P. Lin, and Y-M. Lin.1996Mating systems and spontaneous mutation rates for chlorophyll-deficiency in populations of the mangrove Kandelia candel. Hereditas 125: 4752. [CrossRef][ISI]
Duke, N. C.1992Mangrove floristics and biogeography. In A. Robertson and D. Alongi [eds.], Tropical mangrove ecosystem, 63100. American Geographical Union, Washington, DC.
Elmqvist, T., and P. A. Cox.1996The evolution of vivipary in flowering plants. Oikos 77: 39.
Hamrick, J. L., and M. J. W. Godt.1996Effects of life history traits on genetic diversity in plant species. Philosophical Transactions of the Royal Society of London B 351: 12921298.
Hatcher, B. G., R. E. Johannes, and A. I. Robertson.1989Review of research relevant to conservation of shallow tropical marine ecosystems. Oceanography and Marine Biology: An Annual Review 27: 337414.
Kephart, S. R.1990Starch gel electrophoresis of plant isozymes: a comparative analysis of techniques. American Journal of Botany 77: 693712. [CrossRef][ISI]
Klekowski, E. J., JR., R. Lowenfeld, and P. K. Hepler.1994Mangrove genetics. II. Outcrossing and lower spontaneous mutation rates in Puerto Rican Rhizophora. International Journal of Plant Sciences 155: 373381. [CrossRef]
Lakshmi, M., S. Rajalakshmi, M. Parani, C. S. Anuratha, and A. Parida.1997Molecular phylogeny of mangroves. I. Use of molecular markers in assessing the intraspecific genetic variability in the mangrove species Acanthus ilicifolius Linn. (Acanthaceae). Theoretical and Applied Genetics 94: 11211127. [CrossRef][ISI]
Lee, S. Y.1989Litter production and turnover of the mangrove Kandelia candel (L.) Druce in a Hong Kong tidal shrimp pond. Estuarine, Coastal and Shelf Science 29: 7587. [CrossRef][ISI]
Li, M. S., and S. Y. Lee.1997Mangroves of China: a brief review. Forest Ecology and Management 96: 241259. [CrossRef][ISI]
Lowenfeld, R., and E. J. Klekowski, Jr.1992Mangrove genetics. I. Mating system and mutation rates of Rhizophora mangle in Florida and San Salvador Island, Bahamas. International Journal of Plant Sciences 153: 394399. [CrossRef]
Maxwell, G. S.1995Ecogeographic variation in Kandelia candel from Brunei, Hong Kong and Thailand. Hydrobiologia 295: 5965. [CrossRef][ISI]
Morton, B., and J. Morton.1983The sea shore ecology of Hong Kong. Hong Kong University Press, Hong Kong.
Nei, M.1972Genetic distance between populations. American Naturalist 106: 283292. [CrossRef][ISI]
.1973Analysis of gene diversity in subdivided populations. Proceedings of the National Academy of Sciences, USA 70: 33213323.
Ritland, K.1989Correlated matings in the partial selfer, Mimulus guttatus. Evolution 43: 848859.
, and S. K. Jain.1981A model for the estimation of outcrossing rate and gene frequencies using n independent loci. Heredity 47: 3552. [ISI]
Shields, C. R., T. J. Orton, and C. W. Stuber.1983An outline of general resource needs and procedures for the electrophoretic separation of active enzymes for plant tissue. In S. D. Tanksley and J. J. Orton [eds.], Isozymes in plant genetics and breeding, part A, 443468. Elsevier Science, Amsterdam.
Sun, M., and H. Corke.1992Population genetics of colonizing success of weedy rye in Northern California. Theoretical and Applied Genetics 83: 321329. [ISI]
, and F. R. Ganders.1990Outcrossing rates and allozyme variation in rayed and rayless morphs of Bidens pilosa. Heredity 64: 139143. [ISI]
, and K. Ritland.1998Mating system of yellow starthistle (Centaurea solstitialis), a successful colonizer in North America. Heredity 80: 225232. [CrossRef][ISI]
Tam, N. F. Y., and Y-S. Wong.1997Ecological study on mangrove stands in Hong Kong, vol. 1. Agriculture and Fisheries Department, Hong Kong Special Administrative Region Government.
Tomlinson, P. B.1986The botany of mangroves. Cambridge University Press, Cambridge.
, R. B. Primack, and J. S. Bunt. 1979Preliminary observations on floral biology in mangrove Rhizophoraceae. Biotropica 11: 256277. [CrossRef][ISI]
Umali, R. M., P. M. Zamora, R. R. Gatera, R. S. Jara, A. S. Camacho, and M. Vannucci.1987Mangrove of Asia and the Pacific: status and management. Technical Report of the UNDP/UNESCO Research and Training Pilot Programme on mangrove ecosystems in Asia and the Pacific (RAS/79/002), Quezon City, The Philippines.
Weeden, N. F., and J. F. Wendel.1989Genetics of plant isozymes. In D. E. Soltis and P. S. Soltis [eds.], Isozymes in plant biology, 4672. Dioscorides Press, Portland, OR.
Wendel, J.F., and N. F. Weeden.1989Visualization and interpretation of plant isozymes. In D. E. Soltis and P. S. Soltis [eds.], Isozymes in plant biology, 545. Dioscorides Press, Portland, OR.
Wright, S.1951The genetical structure of populations. Annals of Eugenics 15: 323354. [ISI]
Yipp, M. W., C. H. Hau, and G. Walthew.1995Conservation evaluation of nine Hong Kong mangals. Hydrobiologia 295: 323333. [CrossRef][ISI]
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |