Cécile Fauvelot
email: cecile.fauvelot@unibo.it.

Daniel F R Cleary

Steph B J Menken

F R Cleary

´cile Ce

Daniel F R Fauvelot

AND Steph B J Cleary

Menken

Short-Term Impact of 1997/1998 ENSO-Induced Disturbance on Abundance and Genetic Variation in a Tropical Butterfly

In order to assess the short-term impact of habitat loss after disturbance, we studied Arhopala epimuta (Lepidoptera: Lycaenidae) populations in 5 landscapes in Borneo that were differentially affected by the 1997/1998 El Nin ˜o Southern Oscillationinduced drought and fire. Sampling was conducted before (1997) and after (1998 and 2000) disturbance. This study combined demographic and genetic data inferred from the analysis of 5 microsatellite loci and mitochondrial DNA (mtDNA) control region sequences. Over all 5 landscapes, a total of 313 A. epimuta were sampled over the 3-year survey. Butterfly abundance varied greatly both spatially and temporally (within disturbed landscapes). After the disturbance, a 4-fold population expansion was observed in a small unburned isolate, whereas population extinction was observed in one of the severely burned areas. The analysis of mtDNA sequences in a subsample of 106 A. epimuta revealed no significant spatial or temporal genetic structure. The analysis of 5 microsatellite loci revealed high frequencies of null alleles. Genetic evidence of recent change in population size was found in all 3 unburned landscapes using microsatellites. Congruent to mtDNA, microsatellites failed to detect significant genetic structure according to sampling year or landscapes. Our results suggest that, for mobile species within recently fragmented habitat, habitat loss after disturbance may lead to local population extinction but may augment genetic diversity in remnant local populations because of increased gene flow.

The effect of habitat disturbance on genetic and demographic viability of species has become an important issue over recent decades [START_REF] Wilcox | Conservation strategy: the effects of fragmentation on extinction[END_REF]Saunders and others 1991). The fragmentation of continuous suitable habitat is usually related to disturbance and leads to patchily distributed populations that can be both demographically and genetically modified (Saccheri and others 1998). The extent to which a particular species will be affected depends on the magnitude of the disturbance, the degree of habitat specialization of the species, and its dispersal potential and behavioral response to habitat fragmentation [START_REF] Weins | Metapopulation dynamics and landscape ecology[END_REF]. For instance, strongly dispersing generalist species may actively move among patchy habitats, increasing local genetic diversity and maintaining gene flow among isolated locations [START_REF] Porter | Refugees from lost habitat and reorganization of genetic population structure[END_REF]. In contrast, decreased habitat patch connectivity will be observed in species with limited vagility, leading to increases in genetic divergence between remnant populations [START_REF] Harrison | Genetic and evolutionary consequences of metapopulation structure[END_REF]. Furthermore, habitat losses are expected to result in population size reductions that may in turn increase inbreeding rates, losses of genetic variation, and fixation of deleterious mutations, further reducing the evolutionary potential of remnant isolated populations and increasing their risk of extinction [START_REF] Chakraborty | Bottleneck effects on average heterozygosity and genetic distance with the stepwise mutation model[END_REF][START_REF] Frankham | Inbreeding depression: a threshold effect[END_REF]England and others 2003). After these bottlenecks, population recoveries (both in size and in genetic diversity) will depend on species turnover, number of remaining breeders, stochasticity of reproductive success, species dispersal abilities, and the geographic scale at which the disturbance occurred.

El Nin ˜o Southern Oscillation (ENSO)-induced droughts and fires represent a major disturbance and cause of habitat loss and fragmentation in tropical rain forests. The Indonesian province of East Kalimantan (Borneo) was severely affected by the 1997/1998 ENSO event. More than 5 million ha burned in East Kalimantan alone, making it by far the most severe and widespread event recorded in Borneo [START_REF] Harrison | Repercussions of El Nin ˜o: drought causes extinction and the breakdown of mutualism in Borneo[END_REF]Siegert and others 2001). Human-induced forest fires associated with ENSO-related droughts had a profound impact on ecosystem dynamics, community composition, and species richness (Holmgren and others 2001;[START_REF] Cleary | An examination of scale of assessment, logging and ENSO-induced fires on butterfly diversity in Borneo[END_REF]Laurance 2003).

Butterflies in Borneo have previously been extensively studied in terms of community responses to both logging and ENSO-induced droughts and fires (Willott and others 2000;[START_REF] Cleary | An examination of scale of assessment, logging and ENSO-induced fires on butterfly diversity in Borneo[END_REF]Hamer and others 2003;[START_REF] Cleary | Changes in rain forest butterfly diversity following major ENSO-induced fires in Borneo[END_REF]Cleary and Mooers 2004;Cleary and others 2004). In East Kalimantan, the fires created a mosaic of unburned forests (some of which were affected by drought), isolating local butterfly habitats and reducing species richness in the affected areas [START_REF] Cleary | An examination of scale of assessment, logging and ENSO-induced fires on butterfly diversity in Borneo[END_REF][START_REF] Cleary | Changes in rain forest butterfly diversity following major ENSO-induced fires in Borneo[END_REF]. Butterflies are a particularly good model for studying the genetic effects of disturbance and habitat loss because landscape composition is an important determinant of butterfly population structure (Brookes and others 1997;Lewis and others 1997;Saccheri and others 1998;Keyghobadi and others 1999;[START_REF] Schmitt | Influence of habitat fragmentation on the genetic structure of Polyommatus coridon (Lepidoptera: Lycaenidae): implications for conservation[END_REF]Harper and others 2003;Williams and others 2003).

Here, we studied the effects of disturbance on abundance and population genetic structure in a tropical butterfly, Arhopala epimuta (Lepidoptera, Lycaenidae). Arhopala epimuta is a common forest-dependent lycaenid found in Borneo. Lycaenids represent more than 47% of butterfly (Papilinoidea) species [START_REF] Robins | How many butterfly species?[END_REF], and the Arhopala group of species is extraordinarily diverse, with more than 200 described species [START_REF] Bridges | Catalogue of Lycaenidae and Riodinidae (Lepidoptera: Rhopalocera)[END_REF][START_REF] D'abrera | Butterflies of the Australian region[END_REF] of which 88 are found in Borneo (Seki and others 1991). Little is known, however, about genetic variation within species, and this study is the first to attempt a detailed analysis of this diversity. Arhopala species are usually found in pristine forest and are presumably moderate to poor dispersers as are most lycaenids [START_REF] Dennis | The ecology of butterflies in Britain[END_REF]. The larvae only feed on tree species, with some species-specific variation in diet (some are narrowly oligophagous, others polyphagous). In addition, lycaenid larvae including Arhopala often display a variety of interactions with ants [START_REF] Fiedler | Systematic, evolutionary, and ecological implications of myrmecophily within the Lycaenidae (Insecta: Lepidoptera: Papilionoidea)[END_REF]Pierce and others 2002). Little is known about host plant and ant association in A. epimuta. However, according to Megens and others (2005), A. epimuta feeds on young foliage of Fagaceae (beech family), as do most closely related species. Thus, A. epimuta can be classified as a specialist in terms of its larval host and obligate myrmecophily (i.e., larvae are associated with ants), even though Arhopala atosia (the closest related species) has been associated with at least 3 ant species (Megens and others 2005).

The main goal of our study was to assess the change in A. epimuta abundance and genetic makeup relative to the 1997/1998 ENSO-induced disturbance. We therefore assessed A. epimuta abundance and genetic composition through 3 years of surveys (1997, 1998, and 2000) in 5 different landscapes differentially affected by the 1997/1998 ENSO-induced drought and fire. A standardized method was used to esti-mate abundance within landscapes, and the genetic composition was assessed using 2 types of markers: sequences of the maternally inherited mitochondrial DNA (mtDNA) and 5 microsatellite loci.

Materials and Methods

Study Area

Butterfly surveys were conducted in 1997 (before the fires), 1998, and 2000 in the Balikpapan-Samarinda region of East Kalimantan, Indonesian Borneo [START_REF] Cleary | An examination of scale of assessment, logging and ENSO-induced fires on butterfly diversity in Borneo[END_REF][START_REF] Cleary | Changes in rain forest butterfly diversity following major ENSO-induced fires in Borneo[END_REF]) (Figure 1). During these surveys, 5 differentially disturbed landscapes from 3 regions were sampled. Landscapes in unburned isolates (I1, I2, and I3) and burned forests (B1 and B2) were located in the 5.2 million ha of East Kalimantan that changed from a habitat mosaic of forest with areas of secondary growth to an area dominated by secondary (burned) forest with only remnant unburned patches (Siegert and others 2001). All landscapes are fully described in terms of exact location and dominant vegetation in [START_REF] Cleary | An examination of scale of assessment, logging and ENSO-induced fires on butterfly diversity in Borneo[END_REF]. Briefly, I1 (unlogged primary) and I3 (logged in 1993/1994) landscapes were located in a 108 000-ha isolate within the socalled Gunung Meratus region, which includes a part of the International Timber Concessions Indonesia and adjacent Balikpapan Forest Industries concessions and the Gunung Meratus Protected Forest Reserve (ca. 30 000 ha) (Figure 1, Table 1). Within the so-called Sungai Wain region, I2 landscape was located in a 3500-ha unburned primary isolate that is all that remains of the Sungai Wain Protection Forest (6500 ha of the original 10 000 ha burned during the 1997/1998 ENSO event); the burned landscape B1 was located in the burned part of the nature reserve. The second burned landscape B2 is located in the Wanariset Samboja Research Forest that partially burned during the 1982/1983 ENSO event and severely burned during the 1997/1998 ENSO event. None of the burned landscapes have been commercially logged, although they may have been subjected to unrecorded illegal logging.

Sampling, Preservation, and DNA Extraction

Specimens of A. epimuta were sampled along 8-16 transects of 300 Â 30 m each within each landscape, as described in [START_REF] Cleary | An examination of scale of assessment, logging and ENSO-induced fires on butterfly diversity in Borneo[END_REF]. Along transects, the number of observed A. epimuta was first recorded prior to sampling. Abundance within transects was estimated as the number of recorded individuals per transect per day in order to correct for variations in the time used to sample transects. Abundance within landscapes was estimated as the average abundance observed within sampled transects. Collected individuals were immediately killed after which 3 legs were removed and stored in 96% ethanol. These voucher specimens were preserved (using silica gel) and have been deposited in the Zoological Museum of the University of Amsterdam, the Netherlands. Genomic DNA was isolated from one ground leg using a cetyltrimethyl ammonium bromide protocol [START_REF] Hillis | Molecular systematics[END_REF].

The mtDNA Sequencing and Data Analysis

Amplifications of the mtDNA control region (CR) plus a fragment of the 12S rDNA gene were conducted using primers SeqLepMet and LepAT2B [START_REF] Vila | The utility of the neglected mitochondrial control region for evolutionary studies in Lepidoptera (Insecta)[END_REF]. Polymerase chain reaction (PCR) was performed in 25 ll containing 1Â PCR buffer (HT Biotechnology, Cambridge, UK), 3 mM MgCl 2 , 0.12 mM of each deoxynucleoside triphosphate, 0.5 mg/ml bovine serum albumin, 0.05 lM of each primer, 0.4 U Taq polymerase (HT Biotechnology), and approximately 20 ng DNA. PCRs were carried out in PTC 100 Thermocycler (MJ Research, Waltham, MA): 2 min showing the location of landscapes in which individuals of Arhopala epimuta were sampled. Black rectangles represent burned landscapes, and white rectangles represent unburned isolates. I1 and I3 landscapes were located in the so-called Gunung Meratus region and were sampled in 2000. B1 and I2 landscapes were located in the Sungai Wain Protection Forest and were sampled in 1998 and in 2000. B2 was located in the Wanariset forest and was sampled in 1997, 1998, 1999, and 2000.

Table 1. Sampling information, number of recorded and analyzed Arhopala epimuta (numbers in parentheses represent numbers of individuals for which the mtDNA CR was sequenced), and estimated abundance (average number of individuals per transect per day) within each sampled landscape. Samples are named according to the sampled landscape (see Figure 1) and the year of sampling. The forest status refers to the degree of disturbance of the landscape at the time of sampling. For example, B1 was originally a pristine forest that burned during the 1997/98 ENSO event CA). Sequences were run on an ABI 3100 Genetic Analyzer (Applied Biosystems). Forward and reverse sequences were compared for each individual in BIOEDIT [START_REF] Hall | BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT[END_REF]. Individual sequences were then aligned by eye in GENEDOC [START_REF] Nicholas | GeneDoc: a tool for editing and annotating multiple sequence alignments[END_REF]. For further analysis, indels were coded as substitutions, and dinucleotide indels situated in the [TA] n -like stretch were coded as single-base substitutions (see Results).

Haplotypes were identified, and their relative frequency within populations was calculated using ARLEQUIN version 2.000 (Schneider and others 2000). A minimum spanning network displaying evolutionary relationships between haplotypes was obtained in ARLEQUIN and drawn by hand in Adobe Illustrator 9.0 (Adobe Systems Incorporated, San Jose, CA). Genetic diversity within the population was measured as the haplotype diversity (h, Nei 1987), the nucleotide diversity (p, [START_REF] Nei | Molecular evolutionary genetics[END_REF], and the mean number of pairwise nucleotide differences between individuals (k, Tajima 1983) using DNASP 4.0 (Rozas and others 2003). Tajima's D-test [START_REF] Tajima | Statistical method for testing the neutral mutation hypothesis by DNA polymorphism[END_REF] as well as the haplotype diversity test (h-test; Depaulis and Veuille 1998) was conducted within each sample using DNASP 4.0. For both tests, computer simulations were based on the coalescent process for a neutral infinite-sites model under the assumption of a large constant population size [START_REF] Hudson | Gene genealogies and the coalescent process[END_REF]. DNASP generates the empirical distribution of D and h, obtained from 10 000 simulations under the neutral coalescent process and given the number of segregating sites observed in each sample (the model therefore considers that the number of segregating sites is fixed and that mutations are uniformly and randomly distributed along lineages). DNASP thus estimates the probability of obtaining lower or higher values of these parameters computed from the simulation than the one observed in each sample.

Data analyses of spatial and temporal genetic structure among populations were performed using ARLEQUIN. Genetic divergences between pairwise samples were estimated using U-statistics (U ST based on haplotype frequencies and molecular divergence using pairwise differences). P values were obtained using a nonparametric permutation procedure with 10 000 permutations. The distribution of genetic variation was examined using an analysis of molecular variance (AMOVA; Excoffier and others 1992) to test the null hypothesis that genetic variation was not associated with temporal variation within landscapes or with spatial structure among regions (i.e., Gunung Meratus, Sungai Wain, and Wanariset).

Microsatellite Genotyping and Data Analysis

Genotypes at 5 microsatellite loci recently isolated for A. epimuta [START_REF] Fauvelot | Isolation and characterization of microsatellites in two tropical butterflies, Drupadia theda and Arhopala epimuta (Lepidoptera: Lycaenidae)[END_REF] were amplified according to the described protocol. PCR products were diluted 20 times, and run on a Li-Cor 4200 automatic sequencer. To ensure scoring accuracy and consistency across gels, a locus-specific size standard from individuals with known genotypes was run at both sides of each gel, as well as a commercial size standard (50-350 bp, Li-Cor, Lincoln, NE). Individuals for which no PCR products were visualized on gel were reamplified during a new PCR using a 2 degrees lower annealing temperature.

Genetic diversity within samples was estimated as the observed heterozygosity (H O ) and the mean number of alleles per locus in GENETIX (Belkhir and others 2002), as well as the allelic richness [START_REF] El Mousadik | High level of genetic differentiation for allelic richness among populations of the argan tree [Argania spinosa (L.) Skeels] endemic to Morocco[END_REF] and the unbiased gene diversity (H N , [START_REF] Nei | Molecular evolutionary genetics[END_REF] in FSTAT version 2.9.3.2 [START_REF] Goudet | FSTAT, a program to estimate and test gene diversities and fixation indices[END_REF]. Deviations from Hardy-Weinberg (HW) equilibrium were estimated using the inbreeding coefficient f of [START_REF] Weir | Estimating F-statistics for the analysis of population structure[END_REF] (equivalent to Wright's F IS ), and departures from HW expectations were assessed using exact tests with a Markov chain method in GENEPOP version 3.4 [START_REF] Raymond | An exact test for population differentiation[END_REF]. Linkage disequilibrium between all pairwise locus combinations was tested using the exact probability test in GENEPOP.

Null alleles were detected with MICRO-CHECKER (van Oosterhout and others 2004). When a significant heterozygosity deficiency was observed, null allele frequencies were estimated following [START_REF] Brookfield | A simple new method for estimating null allele frequency from heterozygote deficiency[END_REF] null allele estimator 2, and the genotypes were adjusted by 1) estimating the number of true homozygotes according to HW expectations and 2) replacing the ''false'' homozygotes by genotypes made of one known allele and one missing allele. Even though adjusted data sets no longer contain multilocus genotypes, allelic frequencies can still be used for population genetic analyses, assuming that loci with null alleles are unlinked to other loci (van Oosterhout and others 2004). This corrected data set was therefore used to reestimate the inbreeding coefficient f of [START_REF] Weir | Estimating F-statistics for the analysis of population structure[END_REF] and assess locus conformance to HW expectations as described previously.

To test for evidence of a recent reduction in effective population size, the adjusted microsatellite data set was analyzed using BOTTLENECK 1.2.02 [START_REF] Cornuet | Description and power analysis of two tests for detecting recent population bottlenecks from allele frequency data[END_REF]Piry and others 1999). This test detects significant differences between the gene diversity (H Exp ) and the expected equilibrium heterozygosity (H Eq ) calculated through simulations from the observed number of alleles at each locus, under various mutation models [START_REF] Luikart | Empirical evaluation of a test for identifying recently bottlenecked populations from allele frequency data[END_REF]. A significant excess or deficiency of H Exp compared with the H Eq can be interpreted as a signature of a recent change in population size [START_REF] Luikart | Empirical evaluation of a test for identifying recently bottlenecked populations from allele frequency data[END_REF]. H Eq was calculated under the stepwise mutation model (SMM) and the 2-phased mutation model (TPM) with 95% single-step mutations. Significant differences between H Exp and H Eq were detected using the Wilcoxon signed-rank test.

Pairwise genetic divergences were estimated from the noncorrected data set using multilocus the estimators of F ST of [START_REF] Weir | Estimating F-statistics for the analysis of population structure[END_REF], and their significance was tested with 10 000 permutations of individuals among samples in GENETIX. As the adjusted data set no longer contained multilocus genotypes, we could only perform single-locus AMOVAs using microsatellite allele frequency including null alleles, as described in Keyghobadi and others (2005).

Results

Arhopala epimuta Abundance

Over all 5 landscapes, a total of 587 A. epimuta were recorded during the 3 years survey (Table 1). Abundance varied greatly between landscapes at a set time and between years in specific landscapes. For example, in the Sungai Wain Protection Forest, abundance were estimated as 0.38 and 2.51 individuals per transect per day in 2000 in proximate B1 and I2, respectively, and in this same forest, within the landscape I2, estimations were 0.6 and 2.51 individuals per transect per day in 1998 and 2000, respectively. Within landscapes, there were large variations in A. epimuta abundance among the sampled transects, with the exception of the unburned landscape I2 from which the number of recorded individuals was similar between all 16 sampled transects (data not shown). No A. epimuta were recorded after the 1997/1998 ENSO event in the severely burned landscape B2 in Wanariset forest, despite relatively high butterfly abundance in 1997 (1.3 individuals per transect per day) and similar surveys conducted in 1998, 1999, and 2000 (including a substantial increase in sampling effort in 2000). A total of 313 A. epimuta were collected over the 5 landscapes and 3 years in order to conduct the genetic analysis.

Genetic Diversity and Neutrality Tests-mtDNA

A subsample of 106 A. epimuta was randomly selected to conduct mtDNA sequence analysis. The amplified mtDNA fragment, from which 463-467 continuous nucleotides could be reliably read for all specimens, contained the CR (327-331 bp) plus a fragment of 136 bp of the 5#-end of the 12S rDNA gene. Sequence length variations observed among individuals within the CR originated from length polymorphism within region D (the [TA] n -like stretch) and the homopolymer situated at the 3#-end of the CR described in [START_REF] Vila | The utility of the neglected mitochondrial control region for evolutionary studies in Lepidoptera (Insecta)[END_REF]. Because length polymorphism observed in region D originated from dinucleotide [TA] indels (bp 68 and 69, Table 2), each of these dinucleotide indels were treated as a single-base substitution. Among the 106 individuals of A. epimuta sequenced, we identified 37 haplotypes distinguished by 29 polymorphic sites, from which 14 were singletons (Table 2).

Among haplotypes, the uncorrected sequence divergence ranged from 0.21% to 1.50%. Overall, haplotype diversity (h) was 0.861, nucleotide diversity (p) was 0.00391, and the mean number of pairwise differences (k) was 1.821.

The minimum spanning network (Figure 2) presents the relationships among the 37 mtDNA haplotypes. Haplotypes H1 and H2 (differentiated by one transition) were the most common, representing 32% and 17.9% of overall samples, respectively (Table 3). A total of 7 haplotypes were shared by at least 2 landscapes. All other 30 haplotypes were shared by a maximum of 2 individuals restricted to a single landscape. Within samples, the genetic diversity varied widely both spatially and temporally (Table 3). Tajima's D-tests did not detect deviations from neutral expectations in any of the landscapes (all P . 0.05), but the overall sample significantly deviated from neutral expectations (P 5 0.006). The observed haplotype diversity within I2-98 was significantly higher than expected under the neutral coalescent process (Table 3; h-test, I2-98, average h exp 5 0.704; P . 0.999).

Genetic Diversity and HW equilibrium-Microsatellites

Among 313 individuals analyzed at 5 microsatellite loci, the number of alleles per locus varied from 13 (in AeF10) to 32 (in AeB9). Over all loci, the expected heterozygosity ranged from 0.690 (in B1-00) to 0.778 (in B1-98), whereas the observed heterozygosity ranged from 0.362 (in I3-00) to 0.441 (in I2-98) (Table 4). Allelic richness, based on a minimum sample size of 18 individuals (i.e., omitting B1 samples), ranged from 7.16 (in I2-98) to 8.51 (in I1-00) and was not correlated with population abundance (P . 0.05).

Across all populations, no linkage disequilibrium was detected between loci (all P . 0.05). Sixty percent of locus-population combinations showed significant deviations from HW expectations, after Bonferroni correction of the significance level. Four loci (AeA12, AeB9, AeF10, and AeG5) deviated from equilibrium conditions in almost all samples, whereas one locus (AeH5) met HW expectations in every single sample (Table 4). In all cases, the deviations were caused by heterozygote deficiencies. The most probable cause of the departures was null alleles identified using MICRO-CHECKER (van Oosterhout and others 2004). This was confirmed by the presence of nonamplifying individuals at these 4 loci despite several attempts. Null allele frequencies were estimated for each locus within each sample, and the allelic frequencies were adjusted to account for this bias (see Appendix). Null allele frequencies averaged across all sites were 0.20 ± 0.075, 0.32 ± 0.113, 0.27 ± 0.066, and 0.38 ± 0.077 for Ae12, AeB9, AeF10, and AeG5, respectively (means ± standard deviation). After the allele frequency adjustments, 7 exact tests conducted (20%) revealed deviations from HW proportions after Bonferroni correction. However, these deviations were not associated with a particular sample or loci.

This corrected data set could be used to test for recent population changes because the BOTTLENECK software allows for missing alleles in the genotypes. However, as a minimum of 10 individuals are needed to conduct this test, we excluded B1 samples (B1-98 and B1-00) from this analysis. Recent changes in effective population size (P , 0.05) were detected through heterozygote deficiencies in I1-00, I3-00, I2-00, and B2-97 under the SMM and TPM (Table 5). Using the uncorrected data set to conduct this test led to identical conclusions (data not shown).

Spatial and Temporal Genetic Structure

Over all pairwise comparisons, estimates of genetic differentiation were low with both genetic markers (Table 6). Within regions at time t (Gunung Meratus in 2000, Sungai Wain in 1998 and[START_REF] Harrison | Repercussions of El Nin ˜o: drought causes extinction and the breakdown of mutualism in Borneo[END_REF], none of the pairwise comparisons were significant. Using the uncorrected microsatellite data set, 2 pairwise samples appeared genetically divergent: B2-97 and I3-00 (F ST 5 0.012, P 5 0.009) and B2-97 and I2-00 (F ST 5 0.008, P 5 0.02). Using the adjusted data set, these pairwise comparisons remained significant after Bonferroni correction of the significance level (B2-97 and I3-00, F ST 5 0.004, P , 0.001; B2-97 and I2-00, F ST 5 0.002, P , 0.002), whereas no other pairwise comparison was significant.

Using the mtDNA marker, 100% of the total variance was attributed to variation within the 7 samples (1-group AMOVA; global U ST 5 0.000). Microsatellite singlelocus AMOVAs conducted on adjusted allele frequencies (including null alleles) gave similar findings. Within landscapes, no significant genetic divergence was found between sample years using either the mtDNA or the microsatellite markers.

Discussion

Microsatellite Markers in Lepidoptera and Null Alleles

Developing microsatellite markers in Lepidoptera is a major challenge (Megle ´cz and others 2004; Zhang 2004). Authors report small numbers of useful loci per species, ranging from 4 to 15 (Harper and others 2000;Caldas and others 2002;[START_REF] Cassel | Characterization of microsatellite loci in Coenonympha hero (Lepidoptera: Nymphalidae)[END_REF]Flanagan and others 2002;Keyghobadi and others 2002;Williams and others 2002;Ji and others 2003; Table 2. Sequence differences among the 37 haplotypes found in Arhopala epimuta samples. Dots indicate identical bases. The entire sequence of haplotype H1 was deposited in GenBank (accession number DQ491574) mtDNA CR 5# / 3# 12S RNA gene 1 1 1 1 1 1 1 1 2 2 2 2 2 2 2 3 3 3 3 3 3 3 4 4 6 6 8 9 9 3 5 6 6 7 7 8 8 1 3 4 5 5 8 8 0 0 2 5 5 8 9 1 7 8 9 4 4 7 1 7 1 8 6 9 0 1 1 9 8 2 3 1 5 7 8 6 5 6 1 5 1 As in the majority of studies using microsatellite markers in Lepidoptera (Megle ´cz and Solignac 1998; Keyghobadi and others 1999Keyghobadi and others , 2002Keyghobadi and others , 2005;;Harper andothers 2000, 2003;Anthony and others 2001;Caldas and others 2002;Flanagan and others 2002;Ji and others 2003;Williams and others 2003;Bezzerides and others 2004;Franck and others 2005;van't Hof and others 2005), we found homozygote excesses when compared with HW expectations. Inbreeding and/or undetected population subdivision (Wahlund effect), both affecting all loci, selection, and/or the presence of null alleles (both locus specific), might create such deficiencies. In our study, 1 of the 5 loci met HW expectations, and a notable set of individuals did not amplify at the 4 loci showing HW disequilibrium despite several attempts. These are indications of the presence of null alleles at the 4 loci in disequilibrium. Null alleles can be generated from poor primer annealing due to point mutations or indels in the flanking primer (Pemberton and others 1995), short allele dominance (i.e., differential amplification of size-variant alleles, Wattier and others 1998), or inconsistent DNA template quality. Because we did not sequence microsatelliteflanking regions, we cannot specify the mutation occurrence in the primer region. However, because 1) the null alleles were found over all classes of alleles, 2) null homozygotes were encountered in all loci/populations where deficiencies were observed, and 3) DNA extractions were reconducted for individuals with PCR failures without further success, it is very likely that null alleles identified in 4 of the 5 scored loci resulted from mishybridization of at least one primer.

H1 A --T A A A T A A T A A T T A G G

A A A --C G A A G
Estimated null allele frequencies in our study (0.20-0.38) are high but in the range of null allele frequencies presented in [START_REF] Dakin | Microsatellite null alleles in parentage analysis[END_REF] based on 74 microsatellite loci from a wide range of organisms including Lepidoptera. They are furthermore similar to null allele frequencies obtained for microsatellite loci in other Lepidoptera, for example, the Adonis blue butterfly Polyommatus bellargus (from 0.03 to 0.2 at 5 loci; Harper and others 2003) and the Rocky Mountain Apollo butterfly Parnassius smintheus (from 0.03 to 0.34 at 7 loci; Keyghobadi and others 2005). Although null alleles lead to underestimated heterozygosity within samples, it is a minor source of error in estimating heterozygosity excess [START_REF] Cornuet | Description and power analysis of two tests for detecting recent population bottlenecks from allele frequency data[END_REF] and in parental assessments [START_REF] Dakin | Microsatellite null alleles in parentage analysis[END_REF]. Recently, O'Reilly and others ( 2004) and Peijnenburg and others (2006) showed that estimates of differentiation and the probability of detecting genetic differences among populations both diminished when locus heterozygosities were high and data corrected for null alleles. We found that adjusting our data set according to the presence of null alleles did not significantly modify our results concerning the detection of past population size changes in the samples. Likewise, we obtained similar population differentiation estimations using the mtDNA sequence data set and the uncorrected and the adjusted microsatellite data set.

Together, these findings suggest that null alleles likely do not have a major impact on our results, though we are aware that this survey should be conducted with more microsatellite loci in HW to conclude so.

Population Size Changes

The highly unpredictable nature of fires makes it impossible to design a perfect sampling scheme (comprising pre-and postdisturbance samples) to study the short-term impact of fires on the genetic diversity of butterfly populations. Among our 5 sampled landscapes, 2 were located in severely burned forests (B1 and B2). By chance, the Wanariset landscape (B2) was sampled just before the 1997/1998 ENSO event. In 1997, the Wanariset forest presented a mosaic of 3 for haplotype distribution). Each circle corresponds to one haplotype, and each line between 2 haplotypes corresponds to one mutation (see Table 2 for reference). Small black circles indicate missing haplotypes.

forest in an advanced stage of regeneration together with primary dipterocarp forest (Slik and others 2002) because it had already been partially burned during the 1982/1983 ENSO event. Despite this past disturbance, A. epimuta was the second-most dominant species within Wanariset in 1997 [START_REF] Cleary | Butterfly response to severe ENSO-induced forest fires in Borneo[END_REF], and its relative abundance was almost 3-fold that of the unburned and unlogged forest from Gunung Meratus sampled in 2000. The Wanariset forest was severely affected by the 1997/1998 ENSO fires (more than 90% of the forest burned, [START_REF] Cleary | Butterfly response to severe ENSO-induced forest fires in Borneo[END_REF], and no A. epimuta were recorded in this forest in 1998[START_REF] Harrison | Repercussions of El Nin ˜o: drought causes extinction and the breakdown of mutualism in Borneo[END_REF], and even 2004 (C Fauvelot and DFR Cleary, personal observations), further preventing the assessment of genetic diversity in this landscape after the fires.

After the 1997/1998 ENSO event, A. epimuta was recorded in proximate sampled landscapes within the Sungai Wain Protection Forest (ca. 20 km from Wanariset); abundance, however, within the burned landscape B1 was almost 5 times lower than in the adjoining unburned isolate I2. Because of the lack of sampling in Sungai Wain in 1997, we cannot be sure that the lower abundance observed in the burned landscape is a consequence of the fires. However, the landscapes B1 and I2 were both sampled in an originally large unburned forest prior to the disturbance (the Sungai Wain Protection Forest), so that it is very likely that similar abundances of A. epimuta were observed in both landscapes prior to the 1997/1998 ENSO disturbance. In fact, abundance estimates within the unburned landscape I2 in 1998 and in 2000 were more or less similar among all sampled transects (data available from the authors on request). Therefore, because 1) A. epimuta abundance was lower in B1 than I2, 2) a population expansion was observed postdisturbance in B1 (the relative abundance almost tripled between 1998 and 2000), and 3) the fires exterminated A. epimuta populations in the proximate Wanariset forest that was also severely burned, we can deduce that the fires caused A. epimuta population bottlenecks within the burned forest of Sungai Wain.

2 5 4 1 1 2 3 3 1 9 H 3 2 3 1 1 7 H4 2 2 H5 1 1 H6 1 1 H7 1 1 H8 1 1 H9 1 1 H10 1 1 H11 1 1 2 H12 1 1 H13 1 1 1 3 H14 1 1 H15 1 1 H16 1 1 H17 1 1 H18 1 1 2 H19 1 1 H20 1 1 H21 1 1 H 2 2 2 2 1 1 6 H23 2 2 H24 1 1 H25 1 1 H26 1 1 H27 1 1 H28 1 1 H29 1 1 H30 1 1 H31 1 1 H32 1 1 H33 1 1 H34 1 1 H35 1 1 H36 2 2
According to theoretical predictions with respect to the genetic consequences of population bottlenecks, reduced genetic diversity should be observed in A. epimuta populations affected by the 1997/1998 ENSO event. However, we were not able to obtain reliable estimates of genetic diversity within the burned landscape sampled after the fires because of the low numbers of recorded A. epimuta within this landscape (4 in B1-98 and six in B1-00).

A relative abundance increase was observed in the small unburned isolate within Sungai Wain (I2) between 1998 and 2000 (increase by a factor of 4 in a 2-year period). This increase in A. epimuta abundance could result from 1) random fluctuations associated either with variation in the reproductive success of this species or with the sampling method, 2) a recovery after a population bottleneck, and/or 3) an influx of individuals from the surrounding burned forests. However, because a population size change was reflected in the microsatellite data set in I2-00, and considering that abundance was estimated using a standardized method (i.e., no sampling bias; Cleary 2003), this expansion probably did not result from random fluctuations. Moreover, a significant departure from neutrality was observed in this landscape in 1998 using the mtDNA sequence data set, suggesting that A. epimuta present in the pristine isolate of Sungai Wain forest were affected by the drought (the unburned isolate was surrounded by burned forest and suffered severe drought; [START_REF] Slik | Fire survival of lowland tropical rain forest trees in relation to stem diameter and topographic position[END_REF]. Also, the occurrence of newly recorded mitochondrial haplotypes in I2 in 2000 suggests immigration.

Nevertheless, because we lack a pre-1997/1998 ENSO sampling, we cannot be sure that the population expansion corresponds to a recovery because we cannot totally exclude the possibility of a sampling bias in explaining observed variations in haplotype distributions between 1998 and 2000 (not all A. epimuta were sequenced). Within the Sungai Wain Protection Forest, a much slower population expansion was observed in the burned landscape B1 than the unburned I2. The explanation is certainly that fires affected most of the larval host plants of A. epimuta within B1 when compared with the proximate I2 [START_REF] Cleary | Vegetation responses to burning in a rain forest in Borneo[END_REF]. This also happened in Wanariset, leading to the local extirpation of A. epimuta there. However, the lack of a pre-ENSO sample within Sungai Wain (because of the unpredictability of the forest fires) weakens our conclusions. Nevertheless, the contrasted situation observed in nearby habitats differentially affected (by fires or by drought) is an indication of the importance of unburned forest for the survival of A. epimuta populations.

Within the unburned landscapes I1 and I3 (both sampled in 2000 in Gunung Meratus region) and the burned landscape B2 (sampled in 1997 in Wanariset forest), genetic evidence of population expansions was detected using the microsatellite data set. The population densities of A. epimuta prior to sampling within these landscapes are unknown. However, I3 and B2 are landscapes that suffered from some form of disturbance; I3 was commercially logged in 1993-1994, and B2 was a primary forest that partially burned during the less severe 1982/1983 ENSO event (Slik and others 2002). Therefore, population expansions may have followed these intermediate disturbances in both landscapes. Nevertheless, the population sampled in the pristine Gunung Meratus forest (I1) was expected to be at mutation-drift equilibrium because no major disturbance was recorded prior to 2000 (Cleary and others 2004). As null alleles may cause an underestimation of the heterozygosity within sample, the observation of heterozygosity deficiencies should be interpreted with caution. In addition, Harper and others (2003) found that 12.5% of populations documented as having remained stable displayed excesses of heterozygosity, indicating recent bottlenecks. Yet, because we lack population size estimations of the A. epimuta population in I1 prior to 2000, we cannot conclude that a recent population expansion resulted from migrants from proximate burned landscapes, as in I2, or if the results of the Wilcoxon signed-rank test are biased by low loci number coupled with the occurrence of null alleles.

Nevertheless, population size changes (both reduction and expansion) were detected within A. epimuta populations sampled from the Balikpapan-Samarinda region of East Kalimantan, and the demographic bottlenecks were likely the consequences of the 1997/1998 ENSO disturbance.

Altogether, our results show that even though A. epimuta local populations may be extirpated by the disturbance, the proximate remnant populations located in patchy habitat may ''gain'' from the subsequent fragmentation by receiving immigrants from disturbed sites. Under this model, known as the refugee model [START_REF] Porter | Refugees from lost habitat and reorganization of genetic population structure[END_REF], genetic reorganization by movements of refugees (immigrants) causes deviations from genetic equilibrium, increasing genetic variation within the remaining populations and decreasing differentiation among them. Further investigations at larger spatial scales are therefore needed to infer the consequences of the 1997/1998 ENSO disturbance as we now need to disentangle the levels of genetic diversity resulting from the refugee model from the original levels of genetic diversity within undisturbed forests.

Figure 1 .

 1 Figure 1. Map of the research area, located in the Balikpapan-Samarinda region of East Kalimantan (Indonesian Borneo), showing the location of landscapes in which individuals of Arhopala epimuta were sampled. Black rectangles represent burned landscapes, and white rectangles represent unburned isolates. I1 and I3 landscapes were located in the so-called Gunung Meratus region and were sampled in 2000. B1 and I2 landscapes were located in the Sungai Wain Protection Forest and were sampled in 1998 and in 2000. B2 was located in the Wanariset forest and was sampled in1997, 1998, 1999, and 2000.

Figure 2 .

 2 Figure 2. Minimum spanning network displaying evolutionary relationships among the 37 Arhopala epimuta mtDNA haplotypes (see Table3for haplotype distribution). Each circle corresponds to one haplotype, and each line between 2 haplotypes corresponds to one mutation (see Table2for reference). Small black circles indicate missing haplotypes.

Table 3 .

 3 Distribution of mtDNA haplotypes and genetic diversity in Arhopala epimuta samples

		I1-00	I3-00	B1-98	B1-00	I2-98	I2-00	B2-97	Overall
	H 1	5	9	1	2	4	6	7	3 4
	H

Table 4 .

 4 Genetic diversity at 5 microsatellites loci in Arhopala epimuta samples 5 number of individuals scored per sample, H N 5 unbiased gene diversity[START_REF] Nei | Molecular evolutionary genetics[END_REF]), H O 5 observed heterozygosity, f 5 inbreeding coefficientof[START_REF] Weir | Estimating F-statistics for the analysis of population structure[END_REF], All/loc 5 mean number of alleles per locus, and Ac(18) 5 allelic richness based on 18 individuals (El Mousadik and Petit 1996). * Non conformance to HW proportions.

		I1-00	I3-00	B1-98	B1-00	I2-98	I2-00	B2-97
	N	48	77	6	6	18	110	48
	AeA12
	H N	0.721	0.741	0.783	0.433	0.686	0.731	0.763
	H O	0.458	0.284	0.333	0.500	0.278	0.477	0.500
	f	0.365*	0.617*	0.574	ÿ0.154	0.595*	0.347*	0.345*
	AeB9
	H N	0.886	0.873	0.950	0.767	0.913	0.894	0.911
	H O	0.447	0.378	0.200	0.167	0.556	0.388	0.447
	f	0.496*	0.566*	0.789	0.783	0.392*	0.567*	0.509*
	AeF10
	H N	0.736	0.689	0.833	0.733	0.746	0.698	0.798
	H O	0.489	0.263	0.167	0.667	0.353	0.426	0.283
	f	0.335*	0.618*	0.800	0.091	0.527*	0.390*	0.646*
	AeG5n
	H N	0.693	0.530	0.625	0.833	0.761	0.664	0.708
	H O	0.239	0.192	0.400	0.000	0.294	0.182	0.316
	f	0.655*	0.638*	0.360	1.000	0.614*	0.726*	0.554*
	AeH5
	H N	0.675	0.677	0.700	0.683	0.706	0.687	0.629
	H O	0.565	0.693	0.833	0.667	0.722	0.578	0.596
	f	0.162	ÿ0.024	ÿ0.190	0.024	ÿ0.023	0.158	0.053
	Overall
	H N	0.742	0.702	0.778	0.690	0.762	0.735	0.762
	H O	0.440	0.362	0.387	0.400	0.441	0.410	0.428
	f	0.407*	0.484*	0.503*	0.420*	0.422*	0.442*	0.438*
	All/loc	12.2	12.4	4	3.8	7.2	13.6	10.2
	Ac(18)	8.51	7.51	-	-	7.16	8.24	7.96

N

Table 5 .

 5 Results of the Wilcoxon signed-rank tests computed using BOTTLENECK(Piry and others 1999) for heterozygosity excess or deficiency under TPM and SMM for 5 Arhopala epimuta samples. Significant P values are represented in bold

		No. loci		P value
		TPM	SMM	TPM	SMM
	I1-00
	Deficiency	4	5	0.031	0.016
	Excess	1	0	-	-
	I3-00
	Deficiency	5	5	0.016	0.016
	Excess	0	0	-	-
	I2-98
	Deficiency	2	2	0.594	0.406
	Excess	3	3	0.500	0.688
	I2-00
	Deficiency	5	5	0.016	0.016
	Excess	0	0	-	-
	B2-97
	Deficiency	4	5	0.031	0.016
	Excess	1	1	-	-

Table 6 .

 6 Pairwise F ST estimates using 5 microsatellite loci (above the diagonal) and U ST estimates using mtDNA CR sequences (below the diagonal) between Arhopala epimuta samples

	Appendix A

Table A1 .

 A1 Within-sample microsatellite allele frequencies observed at 5 microsatellite loci and microsatellite allele frequencies adjusted for null allele occurrence according toBrookflied (1996)

	Table A1. Continued
	Locus	I1-00 (48)	I3-00 (77)	B1-98 (6)	B1-00 (6)	I2-98 (18)	I2-00 (110)	B2-97 (48)
	allele	Obs.	Adj.	Obs.	Adj.	Obs.	Adj.	Obs.	Adj.	Obs.	Adj.	Obs.	Adj.	Obs.	Adj.
	Locus 130 133	I1-00 (48) 0.043 0.035 0.149 0.122	I3-00 (77) 0.105 0.077 0.151 0.110	B1-98 (6) 0.250 0.171	B1-00 (6) 0.083 -	I2-98 (18) 0.265 0.184	I2-00 (110) 0.051 0.041 0.167 0.134	0.033 B2-97 (48) 0.021 0.207 0.136
	allele 136	Obs. 0.032	Adj. 0.026	Obs. 0.033	Adj. 0.024	Obs. 0.167	Adj. 0.114	Obs. 0.083	Adj. -	Obs. 0.029	Adj. 0.021	Obs. 0.060	Adj. 0.048	Obs. 0.054	Adj. 0.036
	139 AeA12 144 142 146 145 150 148 151 163 152 Null 153 AeG5 154 73 156 76 158 77 160 78 162 79 164 80 174 82 176 83 Null 84 85 AeB9 118 87 124 88 126 91 128 94 130 97 132 100 134 103 136 106 138 109 140 112 142 Null 144 AeH5 146 155 148 167 150 170 152 171 154 172 156 173 158 174 160 175 162 176 164 177 166 178 168 179 170 182 172 185 174 188 176 206	0.085 0.032 0.135 0.011 0.043 0.208 0.010 0.188 0.438 0.022 0.021 0.011 0.022 0.011 0.011 0.033 0.022 0.033 0.533 0.152 0.043 0.087 0.011 0.033 0.022 0.021 0.011 0.032 0.021 0.021 0.011 0.021 0.011 0.106 0.064 0.011 0.298 0.064 0.053 0.337 0.096 0.011 0.021 0.011 0.457 0.087 0.011 0.011 0.054 0.011 0.033 0.021	0.070 0.026 0.115 0.009 0.035 0.178 0.183 0.009 0.160 0.373 0.015 0.018 0.007 0.015 0.007 0.007 0.022 0.148 0.015 0.022 0.357 0.102 0.032 0.058 0.008 0.022 0.015 0.016 0.007 0.024 0.016 0.016 0.330 0.008 0.016 -0.079 0.047 -0.220 0.047 0.039 -0.071 -0.016 0.008 --0.008 0.008 --0.008 0.016	0.053 0.014 0.169 0.007 0.013 0.007 0.007 0.237 0.169 0.021 0.392 0.007 0.034 0.007 0.034 0.007 0.007 0.089 0.678 0.062 0.007 0.014 0.021 0.020 0.021 0.014 0.014 0.014 0.014 0.014 0.101 0.007 0.142 0.020 0.277 0.007 0.135 0.061 0.313 0.041 0.013 0.027 0.007 0.020 0.467 0.007 0.087 0.007 0.053 0.020 0.007 0.014	0.038 0.009 0.115 0.005 0.010 0.005 0.005 0.160 0.273 0.115 0.014 0.266 0.005 0.023 0.005 0.023 0.005 0.005 0.322 0.061 0.464 0.042 0.005 0.009 0.014 0.014 0.014 0.009 0.009 0.009 0.316 0.009 0.009 0.069 0.097 -0.189 -0.092 0.042 -0.028 -0.018 -0.014 -0.005 -0.005 ---0.009	0.167 0.333 0.167 0.083 0.417 0.100 0.600 0.300 0.100 0.200 0.300 0.200 0.200 0.500 0.250 0.167 0.083	0.114 --0.317 -----0.048 0.096 0.144 0.096 0.096 ----	0.083 0.083 0.167 0.750 0.250 0.500 0.250 0.500 0.083 0.167 0.500 0.250 0.333 0.083 0.083	-------0.353 0.059 0.118 -0.177 ---	0.029 0.111 0.389 0.083 0.417 0.059 0.177 0.029 0.059 0.471 0.088 0.028 0.059 0.028 0.059 0.056 0.083 0.222 0.056 0.083 0.194 0.056 0.083 0.306 0.056 0.056 0.444 0.056 0.111 0.028 0.056	0.021 0.086 0.300 0.304 0.064 0.322 0.039 0.117 0.228 0.020 0.039 0.312 0.059 0.039 0.023 0.039 0.023 0.046 0.337 0.069 0.184 0.046 -0.161 0.046 0.069 -0.046 0.046 ---0.023 0.046	0.046 0.028 0.225 0.009 0.014 0.211 0.151 0.395 0.025 0.005 0.010 0.010 0.015 0.005 0.091 0.005 0.010 0.061 0.005 0.030 0.005 0.566 0.066 0.040 0.005 0.010 0.015 0.035 0.005 0.010 0.056 0.005 0.026 0.005 0.015 0.015 0.077 0.056 0.117 0.028 0.260 0.087 0.009 0.051 0.326 0.066 0.028 0.026 0.026 0.445 0.026 0.023 0.009 0.101 0.005 0.014 0.018 0.010 0.026	0.037 0.022 0.188 0.011 0.008 0.176 0.198 0.126 0.329 0.015 0.004 0.006 0.006 0.009 0.052 0.004 0.003 0.166 0.006 0.035 0.017 0.003 0.003 0.324 0.038 0.023 0.003 0.006 0.020 0.009 0.006 0.003 0.033 0.003 0.003 0.015 0.009 0.427 0.009 0.045 0.033 0.070 -0.154 0.051 -0.030 -0.039 -0.015 0.015 -0.015 ----0.003 -0.006 0.015	0.087 0.033 0.156 0.044 0.156 0.219 0.385 0.040 0.031 0.021 0.031 0.145 0.158 0.079 0.011 0.500 0.011 0.040 0.026 0.011 0.013 0.011 0.021 0.032 0.202 0.128 0.085 0.011 0.085 0.064 0.298 0.128 0.043 0.075 0.011 0.532 0.021 0.021 0.043 0.053 0.032 0.011	0.057 0.021 0.029 0.134 0.343 0.134 0.187 0.020 0.330 0.027 0.018 0.073 0.027 0.080 0.144 0.040 0.253 0.008 0.020 0.008 0.013 0.008 0.007 0.008 0.495 0.016 0.023 0.148 0.093 -0.062 0.062 -0.047 -0.093 0.054 -0.008 -0.016 -0.031 --
	178			0.047	0.032							0.010	0.006	0.043	0.031
	180	0.043	0.032	0.027	0.018							0.015	0.009	0.011	0.008
	182	0.011	0.008	0.014	0.009
	186													0.011	0.008
	Null		0.261		0.318		0.520		0.294		0.172		0.407		0.270
	AeF10
	121	0.011	0.009	0.007	0.005					0.059	0.041	0.023	0.019	0.054	0.036
	123	0.021	0.017	0.007	0.005
	124	0.479	0.391	0.520	0.378	0.417	0.285	0.500	-	0.412	0.287	0.509	0.409	0.380	0.250
	127	0.096	0.078	0.099	0.072			0.250	-	0.206	0.143	0.097	0.078	0.109	0.071
	128											0.005	0.004

Journal of Heredity 2006:97(4)

Number of individuals sequenced per landscape (N), total number of haplotypes (H), haplotype diversity (h, Nei 1987), nucleotide diversity (p, Nei 1987), mean number of pairwise nucleotide differences (k, Tajima 1983), Tajima's D statistic (D, Tajima 1989; *P 5 0.006), and P value associated with the h-test (Depaulis and Veuille 1998) conducted in DNASP.

* P , 0.05.

Obs., observed; adj., adjusted.

Acknowledgments

All field assistants in Kalimantan are thanked for their help. We further thank W. van Ginkel for sorting all the materials and P. Kuperus and B. Voetdijk for laboratory assistance. We thank M. Vila for her valuable advice concerning PCR amplifications and J. A. J. Breeuwer and K. T. C. A Peijnenburg who contributed to this work through insightful discussions. B. Bowen and 2 anonymous reviewers provided helpful comments on earlier versions of the manuscript. This work was supported by grant 895.100.029 of the Netherlands Foundation for the Advancement of Tropical Research, within the Priority Program ''Biodiversity in disturbed Ecosystems.''