Contributions to Zoology, 71 (4) (2002)Daniel F.R. Cleary; Henri Descimon; Steph B.J. Menken: Genetic and ecological differentiation between the butterfly sisterspecies Colias alfacariensis and Colias hyale

To refer to this article use this url: http://dpc.uba.uva.nl/ctz/vol71/nr04/art03

Results

Habitat preference

next section

Along the synthetic environmental axis (PCA-1) C. australis has a highly significant (R2 = 0.599, F1,33 = 49.20, P < 0.001) positive relation to PCA-1 whereas C. hyale has a highly significant (R2 = 0.599, F1,33 = 49.24, P < 0.001) negative relation (Fig. 1). Of the six environmental variables, vegetation and altitude were positively correlated with PCA-1, aspect and humidity negatively correlated, and date and weather uncorrelated. Humidity is thus probably the most important environmental determinant of both species since there is an obvious relation (e.g., exposed ridges tend to be drier than valley bottoms) between humidity and the other three variables that had high loadings for PCA-1 (vegetation, altitude, and aspect).

FIG2

Fig. 1. Regression of C. alfacariensis and C. hyale onto the first Principal component axis (PCA-1). The slope of both lines is highly significant (P < 0.001).

Population analysis

Four of the loci (ald, mdh1, ck2, and got1) are monomorphically fixed for the same allele in the populations of both species. The remaining set of 12 allozyme loci did not contain loci that were diagnostic between C. alfacariensis and C. hyale (Appendix I). Four loci (got2, mpi, pgi, and pgm) show a very high degree of genetic polymorphism, but got2 is monomorphic in all populations of C. hyale.An intermediate level of genetic polymorphism is found among the loci hk, idh1, and me. Low levels of genetic variability are present in the remaining loci (Appendix 1).

In general, C. alfacariensis tends to have both a higher mean number of alleles and percentage of polymorphic loci (Table 2). Consequently, Colias alfacariensis also has a significantly higher heterozygosity (range: 0.148-0.194) than C. hyale (range: 0.084-0.094) (F1,11 = 72.70, P < 0.001). Heterozygosity is highest for C. alfacariensis populations in the geographically complex Alps and drops off towards the lower-lying regions to the South (Provence and Italy) and the North (Metz and Seine), with the lowest values occurring in Northern Europe. Within C. alfacariensis, the Alpine/Italian on the one hand and the Provence/Northern populations on the other hand consistently differ at the PGI locus (Fig. 2).

 

Table 2. Population genetics summary statistics at 12 loci in all populations (standard errors in parentheses). * A locus is considered polymorphic if the frequency of the most common allele does not exceed 0.990

Species

Population

Sample size per Locus

Mean no.of alleles per locus

Percentage of loci polymorphic*

Mean heterozygosity Direct-count

C. australis

A

19

2.2 ± 0.4

56.3

0.17 ± 0.062

C. australis

B

8

1.8 ± 0.3

43.8

0.19 ± 0.073

C. australis

C

19

2.1 ± 0.3

62.5

0.19 ± 0.061

C. australis

D

9

1.9 ± 0.3

56.3

0.19 ± 0.059

C. australis

E

22

2.3 ± 0.4

50.0

0.17 ± 0.057

C. australis

F

21

1.9 ± 0.3

37.5

0.16 ± 0.066

C. australis

G

14

2.2 ± 0.3

62.5

0.17 ± 0.052

C. australis

H

7

1.8 ± 0.3

43.8

0.2 ± 0.077

C. australis

I

32

1.9 ± 0.3

56.3

0.15 ± 0.049

C. australis

J

12

1.8 ± 0.2

50.0

0.16 ± 0.058

C. australis

Mean

2.0

51.9

0.175

C. hyale

K

23

1.8 ± 0.3

50.0

0.08 ± 0.04

C. hyale

L

12

1.6 ± 0.2

31.3

0.09 ± 0.046

C. hyale

M

20

1.9 ± 0.3

43.8

0.08 ± 0.039

C. hyale

Mean

1.8

41.7

0.087

FIG2

Fig. 2. Map of the localities where the populations of C. alfacariensis were sampled. Five main localities are shown, representing a single population with the exception of the Alps and Provence where three and four populations, respectively, were sampled. Within each circle average heterozygosity of 12 loci and the frequency of the pgi-b allele are given. For the Provence and Alps the values are means of the three populations. H = heterozygosity of a single populations, Hm = mean heterozygosity from more than one population.

In C. alfacariensis, estimated single locus FIS significantly differs from zero at me (Table 3). In C. hyale the PGM locus significantly differs from zero. FIS values were generally slightly positive in both species indicating a small deficiency of heterozygotes. In both species total genetic diversity (FIT ) was mainly the result of inbreeding (FIS ), among population variation (FST ) being less important. Within populations of both species the most common allele at a locus is usually the same (Appendix 1). Within C. hyale none of the loci have a FST value significantly different from zero indicating that, at least until recently, there was a panmictic population structure with high levels of gene flow. Over the same geographic range, C. alfacariensis has an estimated Ne M of less than 10. The difference here though is solely due to the PGI locus, which effectively divides the species into a Northern-Provence and Alpen-Italian component. Excluding this locus gives a panmictic population structure, since none of the other loci differ significantly from zero, quite similar to its sister species.

 

Table 3. Estimates of Wrights F-statistics (FIS , FST , and FIT ) at 12 allozyme loci for C. alfacariensis and C. hyale. Significant departures from zero: **: P < 0.01; *** P < 0.001. Average and (standard deviation) obtained by a bootstrap procedure over samples (1000 repeats).

Colias australis

Colias hyale

Locus

FST

FIT

FIS

FST

FIT

FIS

ck1

0.003

-0.006

-0.009

-0.009

-0.003

0.006

got2

0.018

0.015

-0.003

-

-

-

gpd

0.002

-0.006

-0.008

-0.009

-0.003

0.006

hk

0.035

0.175

0.145

-

-

-

idh1

0.010

-0.049

-0.059

-0.013

-0.009

0.004

idh2

0.004

-0.010

-0.014

0.004

-0.017

-0.022

mdh2

0.007

0.140

0.146

-0.013

-0.009

0.004

me

0.007

0.219

0.225***

-0.004

-0.011

-0.006

mpi

0.027

0.176

0.153

-0.014

0.265

0.275

6pgdh

0.003

-0.003

0.000

-0.009

-0.003

0.006

pgi

0.120***

0.122

0.003

-0.022

0.007

0.029

pgm

0.000

0.088

0.088

-0.016

0.350

0.360**

Global

0.037

0.107

0.073

-0.016

0.167

0.180

Average

0.029 (0.011)

0.098 (0.022)

0.071

-0.023 (0.008)

0.158 (0.050)

0.177 (3.862)