Genetic variation, its assessment and implications to the conservation of seagrasses

In a study of the widespread Australian endemic seagrass Posidonia australis , allozyme analysis identified a wide range in population genetic structure assessed using the multilocus genotype diversity statistic ( D _G). Values of D _G between zero and one were obtained; however, RAPD analysis generally detected higher levels of diversity, where D _G values were all greater than 0.5 ( D _G = 0.67 – 1). Some populations were allozymically monomorphic using allozyme analysis yet were highly polymorphic using RAPD analysis. The differences observed between methods, particularly among allozymically uniform populations, demonstrate the importance of choosing an appropriate method when assessing genotypic diversity. Different methods may reflect different historical aspects of population processes where allozymes reflect broader-scale gene flow and population establishment and DNA fingerprinting methods such as RAPDs may reflect fine-scale local recruitment events and shorter-term population processes. Using either method alone, particularly in genotypically depauperate organisms such as seagrasses and other clonal organisms, will be problematic in assessing their population genetic potential, a parameter being used by conservation managers to decide upon management strategies in rare and endangered organisms. It is recommended that the impact of disturbance assessed using genotypic diversity measures requires more than one technique to provide the most appropriate information for designing subsequent conservation strategies.     

Allozyme, mitochondrial-DNA, and morphometrie variability indicate cryptic species of anchovy (Engraulis encrasicolus)

Previous surveys of population structure in the Atlantic-Mediterranean anchovy Engraulis encrasicolus L. have reported heterogeneity in morphology, allozyme frequencies, and mitochondrial DNA haplotype frequencies at a regional scale. In particular, two stocks of anchovy have been detected in the Adriatic Sea. In this paper, the available data is reviewed with the aim to relate genetic variation to geography at the widest possible geographical scale, for investigating the evolutionary mechanisms underlying stock structure in anchovy. Correspondence analysis of allozyme frequencies (24 samples, three polymorphic loci) compiled from the literature indicates three distinct entities in the Mediterranean Sea. Open-sea or oceanic anchovy populations are genetically different from inshore-water populations within a region (Nei's ^ G _ST = 0.035–0.067), while broadscale geographical variation is weak for each of these two habitat-specific forms (^ G _ST = 0.005–0.006). Mitochondrial-DNA haplotype frequencies support the distinction between an inshore form and an oceanic form (^ G _ST = 0.067–0.107), with virtually no genetic differences among oceanic populations across the Gulf of Biscay, the western Mediterranean and the Ionian Sea (^ G _ST = −0.001). If natural selection on marker loci is unimportant, these results indicate the occurrence of two parapatric, genetically distinct, habitat-specific forms that are widely distributed throughout the Mediterranean Sea. Persistent allele and haplotype-frequency differences between these forms indicate reproductive isolation and the presence of an E. encrasicolus species complex in the Mediterranean. © 2002 The Linnean Society of London, Biological Journal of the Linnean Society 75 : 261–269.     

Population structure, gene flow and evolutionary relationships in four species of the genera Tomocerus and Pogonognathellus (Collembola, Tomoceridae)

Genetic diversity and gene flow at 10 allozyme loci were investigated in 13 natural populations of four species of the collembolan genera Tomocerus and Pogonognathellus. Levels of observed heterozygosity were found to vary from 0.033 in P. longicornis to 0.120 in T. vulgaris. Average Nei's genetic distances (D) ranged from D = 0.222 between populations of P. ftavescens to D = 0.647 between populations of T. vulgaris. Genetic distances between species were always 1.000 and the highest value (D = 4.321) was between P. longicomis and one population of T. vulgaris. Values of F_sr were very high in all species (from 0.567 to 0.696) and levels of gene flow (Nm) derived from F_ST and the private allele method were low (Nm <1). Gene flow was significandy higher only in two subsets of populations of T. vulgaris. The Plio-Pleistocenic geological rearrangements and the effect of stochastic events, such as genetic drift, are invoked to explain the different levels of genetic divergence between and within species. Distance- and character-based approaches were used to reconstruct evolutionary relationships between and within species. While the monophyly of all species was confirmed, the results did not univocally support the monophyly of the two genera, leaving the question of their generic or subgeneric status unresolved.     

Karyotype structure, supernumerary chromosomes and heterochromatin distribution suggest a pathway of karyotype evolution in Dactylorhiza (Orchidaceae)

The relationships between Dactylorhiza romana and D. saccifera from southern Italy were analysed. These two species, both with 2 n = 2 x = 40 chromosomes and belonging to different sections of the genus, were distinguishable on the basis of karyotype structure and heterochromatin amounts and distribution. Their C-banded karyotypes differed considerably. D. saccifera showed most chromosomes with banded regions in the short arms, whereas in D. romana the bands were located mostly at telomeric regions of longer arms. Several individuals of D. romana had one or two very large heterochromatic supernumerary chromosomes. Based on evidence resulting from karyotype structure and heterochromatin distribution in the two species and on the genetic distances derived from the comparison of ITS sequences, it is suggested that D. romana represents a primitive form with respect to D. saccifera and is a possible intermediate step in the evolution of the genus Dactylorhiza from the 42-chromosome Orchis group. © 2002 The Linnean Society of London, Botanical Journal of the Linnean Society , 138 , 85–91.     

Population genetic structure of wild rice Oryza glumaepatula distributed in the Amazon flood area influenced by its life-history traits

Wild plant species develop their own way of living to adapt to the specific environment of their habitats. Their life-history traits strongly affect the genetic structure of the population. The wild species Oryza glumaepatula Steud. growing in the Amazon basin seems to have characteristic life-history traits suited for the flood condition. At the vegetative growth stage, the culms frequently break at internodes. With no roots anchoring on the ground, plant bodies floating in the water move downriver by water current and wind. To examine the association between the life-history traits and genetic population structure of Amazonian O. glumaepatula , we analysed allozyme variability at 29 loci of 16 enzymes using 37 populations from five regions. Allozymes were not so variable (total gene diversity H _E = 0.044) compared with Asian wild rice, O. rufipogon Griff. The bottleneck effect and rare opportunity of interspecies gene flow may prevent the development of allozyme variability. Population genotypes tended to be differentiated among geographically isolated regions. Observed heterozygosities were much lower than expected heterozygosities, or gene diversity ( H _O = 0.003 for whole population) and F _IS over polymorphic loci was 0.931, indicating that O. glumaepatula has developed an inbreeding system. But, the intrapopulation gene diversity ( H _S) was higher than interpopulation gene diversity ( D _ST), as generally observed in outbreeding populations. The migration ability of O. glumaepatula makes long-distance seed dispersal possible. This might have led to frequent gene flow among populations.     

Population genetics of a cyclostome species pair, river lamprey (Lampetra fluviatilis L.) and brook lamprey (Lampetra planeri Bloch)

Horizontal agarose gel electrophoresis of 24 allozyme loci in four species of Central European lampreys (321 Lampetra planeri , 83 L. fluviatilis , 11 Eudontomyzon mariae and nine Petromyzon marinus ) was used to study the 'paired species' L. fluviatilis and L. planeri . The genetic differentiation of the anadromous river lamprey ( L. fluviatilis ) from the stationary brook lamprey ( L. planeri ) was within the range of ingroup differentiation of the latter, but L. fluviatilis exhibited much greater population cohesion over a more extended geographic range: G _ST = 0.0537 versus G _ST = 0.3398, N _em = 4.402 versus N _em = 0.4856, mean genetic among-stock distances D = 0.0047 versus D = 0.0257. L. planeri populations coexisting geographically with L. fluviatilis in the Rhine and Elbe river systems were genetically more cohesive than L. planeri stocks from the Danubian basin where L. fluviatilis is absent. Danubian L. planeri populations exhibit a lower degree of heterozygosity than brook lampreys from the Rhine river system, but comprise deeper genetic lineages ( G _ST = 0.4629 versus G _ST = 0.2434), despite being sampled from a much more restricted area. Isolation-by-distance is observed for L. planeri from the Danubian but not from the Atlantic drainage basins. Transspecific gene flow between L. planeri from Atlantic drainage basins and the long-distant migrating L. fluviatilis is inferred, raising doubt on the validity of two separate biospecies. E. mariae and P. marinus are clearly differentiated from Lampetra spp. at several allozyme loci.     

Genotypic diversity revealed by allozymes and oligonucleotide DNA fingerprinting in French populations of the aquatic macrophyte Sparganium erectum

The vast majority of perennial plants reproduce sexually and vegetatively at the same time. This may lead to important variation among clonal plant populations in their amount of genotypic diversity. In order to verify this assumption, we compare the clonal diversity of 10 natural populations of the aquatic clonal macrophyte Sparganium erectum in France. Diversity was quantified by DNA fingerprinting and allozyme electrophoresis for a sample of 10 shoots per population. Two DNA probes (CA)₈ and (TAA)₆ TA, were selected among 10 synthetic oligonucleotide probes containing simple repeat motifs. Both allozymes and DNA fingerprints revealed different amounts of diversity among populations. Five populations consist of a single genotype, whereas two populations were genetically highly diverse. In four of the monomorphic populations, absence of fingerprints diversity was combined with uniformly heterozygous allozyme loci, suggesting that each population was composed of a single clone. In the highly diverse populations, the level of clonal diversity combined with the allele segregation of the two allozyme loci Lap and Est suggests frequent seedling recruitment. The origin of new genotypes remains unclear but the absence of widespread genotypes together with the discrete distribution of Sparganium erectum populations implies that new genotypes are locally produced through sexual reproduction.     

Genetic variability in Irish populations of the invasive zebra mussel, Dreissena polymorpha: discordant estimates of population differentiation from allozymes and microsatellites

1. The recent arrival and explosive spread of the zebra mussel, Dreissena polymorpha (Pallas), in Ireland provided a rare opportunity to study the population genetics of an invasive species.
2. Eight polymorphic allozyme loci ( ACO-1, ACO-2 , EST-D, GPI, IDH-2, MDH, OPDH and PGM ) were used to investigate genetic diversity and population structure in five Irish populations, and the results were compared with those from a previous microsatellite study on the same samples.
3. The mean number of alleles per locus (2.7 ± 0.1) was similar to the mean for the same loci in European populations, suggesting that Irish founder populations were large and/or multiple colonization events took place after foundation. A deficiency of heterozygotes was observed in all populations, but was uneven across loci.
4. Pairwise comparisons, using Fisher's exact tests and F _ST values, revealed significant genetic differentiation among populations. The overall multilocus F _ST estimate was 0.118 ± 0.045, which contrasted with an estimate of 0.015 ± 0.007 from five microsatellite loci on the same samples in a previous study.
5. Assuming that microsatellites can be used as a neutral baseline, the discordant results from allozymes and microsatellites suggest that selection may be acting on some allozyme loci, specifically ACO-1, ACO-2 , IDH-2 and MDH, which contributed most to the significant differentiation between samples.     

Extensive clonality in the endangered shrub Haloragodendron lucasii (Haloragaceae) revealed by allozymes and RAPDs

The occurrence of clonality in threatened plants can have important implications for their conservation. In this study, allozymes and RAPDs were used to determine the extent of clonality in the endangered shrub Haloragodendron lucasii (Haloragaceae), which is known from only four sites within an 8 km range. Allozyme markers identified only six multilocus genotypes among the 53 ramets sampled across the four sites, although a total of 54 different genotypes were possible with the three polymorphic allozyme loci detected. The polymorphic bands detected in the RAPD analysis were capable of producing 2⁴⁶ genotypes, but again only six multilocus genotypes were delineated. The allozyme and RAPD data were congruent at three of the four sites. At the fourth site two genotypes were detected by each marker; however, once combined, three multilocus genotypes were observed. The probabilities that the observed number of replicates of each combined allozyme and RAPD genotype could be generated by sexual reproduction were less than 10^–18, leaving little doubt that clonality is the explanation for the observed patterns of genotypes. The genetic conclusions are supported by root excavations which show potential for vegetative reproduction and the observation of no sexual reproduction in the species. The recognition of extensive clonality in H. lucasii has had immediate implications for the conservation management of the species and resulted in changes to the management priorities for the species. Thus it is clear that appropriate genetic studies can play an important role in the management of threatened species.     

Allozyme Diversity in Populations of Cymbidium goeringii (Orchidaceae)

Abstract: Using 14 allozyme loci, we investigated levels of genetic diversity within populations, and degree of genetic divergence among 24 populations of Cymbidium goeringii (Orchidaceae) in Korea and Japan. Cymbidium goeringii maintains high levels of genetic diversity both at population (mean expected heterozygosity, H _e = 0.238) and species levels (0.260). Means of H _e found in 24 populations were not significantly different from each other. About 90 % of the total variation in the species is common to all populations (mean G _ST = 0.108). No unique allele was found in any population. The indirect estimate of gene flow based on the mean G _ST was high ( Nm = 2.06). Nei's genetic identities for pairs of populations had high values (mean = 0.974 [SD = 0.013]). The Mantel-Z test showed a significant correlation between genetic distance and geographic distance. However, the mean G _ST value between 17 populations in Korea and seven Japanese populations was relatively low (0.029), even though the land connection between the southern Korean peninsula and southern Japanese archipelagos has not existed since the middle Pleistocene. Large numbers of small seeds of C. goeringii might travel long distances by wind from populations to populations both in Korea and Japan, increasing genetic diversity within populations and maintaining low genetic differentiation among populations.     

Chloroplast DNA phylogeography of the argan tree of Morocco

Polymorphisms in the chloroplast genome of the argan tree (Sapotaceae), an endemic species of south-western Morocco, have been detected by restriction site studies of PCR-amplified fragments. A total of 12 chloroplast DNA (cpDNA) and two mitochondrial DNA (mtDNA) fragments were amplified and digested with a single restriction enzyme ( Hin fI). Polymorphisms were identified in six of the cpDNA fragments, whereas no mtDNA polymorphisms were detected in a survey of 95 individuals from 19 populations encompassing most of the natural range of the species. The cpDNA polymorphisms allowed the identification of 11 haplotypes. Two lineages, one in the south-east and the other in the north-west, divide the range of the argan tree into two distinct areas. The level of genetic differentiation measured at the haplotype level ( G _STc= 0.60) (i.e. with unordered haplotypes) was smaller than when phylogenetic relationships were taken into account ( N _STc= 0.71–0.74) (ordered haplotypes), indicating that population history must be considered in the study of the geographical distribution of cpDNA lineages in this species. If contrasted with the level of nuclear genetic differentiation measured in a previous study with isozymes ( G _STn= 0.25), the results indicate a relatively high level of gene flow by seeds, or conversely a relatively low level of gene flow by pollen, as compared with other tree species. Goats and camels could have played an important role in disseminating the fruits of this tree.     

Morphometric and genetic characterization of sympatric populations of Clarias gariepinus and C. anguillaris from Senegal

A sample of African Clarias catfishes from the Senegal River was studied using morphometry, allozyme variation, microsatellites and RFLPs of mitochondrial DNA. They all confirmed the presence of two species , C. gariepinus and C. anguillaris . The two species were closely related genetically and no diagnostic loci were found in allozymes and microsatellites studies. Two of the 11 haplotypes of mtDNA observed were shared by both species. Three of the four assays (morphometry, allozymes and microsatellites) allowed a precise characterization of both. One specimen occupied an intermediate position in the analysis of the data; it was considered an F_l hybrid whose possible origin is discussed.     

PATTERNS OF GENETIC ARCHITECTURE FOR LIFE-HISTORY TRAITS AND MOLECULAR MARKERS IN A SUBDIVIDED SPECIES

Abstract Understanding the utility and limitations of molecular markers for predicting the evolutionary potential of natural populations is important for both evolutionary and conservation genetics. To address this issue, the distribution of genetic variation for quantitative traits and molecular markers is estimated within and among 14 permanent lake populations of Daphnia pulicaria representing two regional groups from Oregon. Estimates of population subdivision for molecular and quantitative traits are concordant, with Q _ST generally exceeding G _ST. There is no evidence that microsatellites loci are less informative about subdivision for quantitative traits than are allozyme loci. Character-specific comparison of Q _ST and G _ST support divergent selection pressures among populations for the majority of life-history traits in both coast and mountain regions. The level of within-population variation for molecular markers is uninformative as to the genetic variation maintained for quantitative traits. In D. pulicaria , regional differences in the frequency of sex may contribute to variation in the maintenance of expressed within-population quantitative-genetic variation without substantially impacting diversity at the genic level. These data are compared to an identical dataset for 17 populations of the temporary-pond species, D. pulex .     

Allozyme Variation and the Genetic Structure of Populations of the Rare Sedge Carex misera (Cyperaceae)

Abstract Carex misera is a rare sedge, endemic to rocky outcrops and mountain summits within the southern Appalachian Mountains from northern Georgia to northern North Carolina. We assessed allozyme diversity for 406 individuals from nine populations over most of the geographic range. Twenty-seven putative loci were assayed and eight (30%) were found to be polymorphic. Nei's gene diversity statistics ( H _T =0.043, H _S =0.019, G _ST =0.551) indicated low levels of variation but relatively highly differentiated populations, suggesting little gene flow. Significant deviations from genotypic expectations under Hardy-Weinberg equilibrium, high positive fixation indices, and the existence of small genetic neighborhoods within populations suggest that at least some inbreeding occurs. Cluster analysis of Nei's genetic identity statistics and principal component analysis of allele frequency data showed high similarity among the six southern populations with the two northern populations more differentiated from them and from each other. These results suggest that preservation of the northern populations is necessary to conserve the already low levels of genetic diversity within the species.     

Local differentiation and gene dispersal distances in a Danish population of Armaria maritima

In a continuously distributed population of the self-incompatible species Armeria maritima , nine artificially delimited subpopulations were sampled in a transect with a length of 1300 m. Genetic variation was quantified with enzyme electrophoresis, which revealed five polymorphic loci. The genotypic distribution within subpopulations accorded to Hardy-Weinberg proportions, but revealed a Wahlund effect at one of the loci when the total population was considered. The Wahlund effect was supported by significant differences in allele frequencies among the nine subpopulations. The differentiation was low with a F_ST value of 0.028, which is much smaller than the previously reported value of 0.18 for 17 geographically separated Danish populations of this species. The low differentiation of the subpopulations was also revealed by an assignment test. In this test, the multilocus genotype of an individual is used to assign it to the subpopulation it has the highest likelihood to appear in. With this test, a lower percentage of individuals were assigned to the subpopulation of origin (24) than was seen with assignment of individuals to populations in Denmark (34). No abrupt changes in allele frequencies were observed. A spatial autocorrelation analysis suggested a breeding patch size of 500–700 m in diameter along the transect.     

Genetic characterization of Dirofilaria repens and D. immitis by electrophoretic analysis of gene-enzyme systems

A multilocus genetic analysis based on protein electrophoresis was carried out on 39 adult specimens of Dirofilaria repens and 31 of D. immitis. Seventeen enzymatic loci provided reliable electrophoretic patterns which were compared in the two species. All loci except one were monomorphic. Fixed alternative allozymes were found at 13 loci while only 3 loci shared apparently the same allozyme in the two species. The polymorphic locus, Pgm, also showed alternative allozymes. This remarkably high genetic divergence, which presumably reflects a very old speciation process, allows an easy characterization of D. repens and D. immitis and supports their classification in the subgenera Nochtiella and Dirofilaria, respectively.     

Natural hybridization between endangered and introduced species of Pseudorasbora, with their genetic relationships and characteristics inferred from allozyme analyses

From 25 populations of Pseudorasbora in Japan, polymorphism at 22 allozyme loci indicated that the level of genetic differentiation between Pseudorasbora parva and Pseudorasbora pumila ( D  = 0·421–0·517) was greater than that between the two subspecies of P. pumila , P. pumila pumila and P. pumila subsp. ( D  = 0·164), consistent with morphological differences. While P. parva displayed genetic variation ( H  = 0·003–0·100) similar to other freshwater fishes, P. pumila pumila and P. pumila subsp. populations showed no genetic variation. In five of the 15 populations collected from the contact zone between P. parva and P. pumila pumila , hybrids were detected by allozyme analyses. All hybrids were presumed to represent the F₁ generation, because they were heterozygous at all 12 loci diagnostic between P. parva and P. pumila pumila . Although four populations were characterized by high frequencies of F₁ hybrids ( c . 40%), only one of the parental species was observed in each case. The results indicated that the two species have hybridized easily under natural conditions, but cannot coexist in the long-term. It is suggested that continued invasion of P. parva would hasten the extinction of P. pumila pumila .     

Evolutionary relationships of the Lake Malawi Oreochromis species: evidence from allozymes

Three tilapiine species belonging to the endemic Lake Malawi species flock known as' chambo' Oreochromis ( Nyasalapia ) karongae , O. (N.) lidole and O. (N.) squamipinnis and the species Oreochromis (Oreochromis) shiranus were collected from the wild. These four species were analysed at 43 enzyme loci using starch gel electrophoresis. No significant deviations from Hardy-Weinberg equilibrium were observed after correcting for multiple simultaneous testing. The expected heterozygosity ( H _e) was lowest in O. (O.) shiranus ( H _e=0.082). The three chambo species had higher levels: H _e=0.103 to H _e=0.116. The three chambo species as a whole were polymorphic (99%) at 16 different loci sharing 10 of these in common. No fixed differences between the three chambo species were observed but highly significant allele frequency differences existed between all the chambo species. The F _ST calculated for the three chambo species was 0051, closer to intra-specific than to inter-specific levels. O. (O.) shiranus could be clearly separated from all chambo species at five fixed loci. Comparison of the allozyme data from these species with five species from the same sub genera support the hypothesis that the chambo form a monophyletic group and that O. (N.) macrochir or a related species represents the sister taxon.     

Systematics and evolution of the Dactylorhiza romana/sambucina polyploid complex (Orchidaceae)

The European–Mediterranean–Oriental Dactylorhiza romana/sambucina polyploid complex was studied with regard to genetic and morphological variation patterns. Allozyme and morphometric data were collected from 24 and 19 populations, respectively, initially identified as D. flavescens, D. insularis, D. markusii, D. romana, D. sambucina, and an indeterminate taxon. Genetic distances were calculated and illustrated by an unweighted pair‐group method using arithmetic averages (UPGMA) dendrogram, and principal components analyses (PCAs) were used to summarize morphological variation patterns. Another PCA was performed on combined allozyme and morphometric data. On the basis of the dendrogram and the PCA plots, main groups of populations were delimited, and the probability that each morphological character would distinguish correctly between these groups was estimated. After combining morphometric interpretations with studies of herbarium material and information from the literature, the following taxa were confidently accepted: D. romana ssp. romana, D. romana ssp. guimaraesii (comb. et stat. nov.) , D. romana ssp. georgica, D. sambucina, D. cantabrica (sp. nov.) , and D. insularis. Levels of genetic diversity suggest that D. romana s.s. is the least derived member of the complex. The evolutionary divergence of the diploid species, D. romana and D. sambucina, was probably the outcome of vicariant speciation, whereas D. romana ssp. georgica and D. romana ssp. guimaraesii appear to have evolved from D. romana s.s. through incomplete vicariant and peripheral isolate speciation events, respectively. In some populations of the diploid taxa, a significant deficiency in heterozygotes was found at one to three loci. It is proposed that this pattern may indicate a Wahlund effect, hypothesizing that local populations are subdivided into demes determined by the commonly sympatric occurrence of two distinct colour morphs combined with partial morph constancy of individual pollinators (bumblebees). Several pathways are possible for the origin of the allotriploid D. insularis and the apparently allotetraploid D. cantabrica. A taxonomic revision is provided. © 2006 The Linnean Society of London, Botanical Journal of the Linnean Society, 2006, 152 , 405–434.     

Genetic differentiation between alpine and lowland populations of a butterfly is related to PGI enzyme genotype

Understanding the ecological process of population differentiation and identifying the molecular changes that contribute to adaptation is a central aspect of evolutionary biology. In this study we analyzed the geographic variation in allozyme allele frequencies (based on 15 enzyme systems representing 18 loci) across 18 populations of the butterfly Lycaena tityrus from different altitudes. Population genetic analyses showed that within population genetic diversity (e.g. mean number of alleles per loci: 1.8; expected heterozygosity: 12%) was within the typical value range for the Lepidoptera. The populations of L. tityrus investigated showed a remarkable genetic differentiation (F_ST: 0.065), being clearly divided into an alpine (high-altitude) and a non-alpine (low-altitude) cluster. This differentiation was almost entirely caused by variation at a single enzyme locus, PGI. Although the involvement of historical events cannot be ruled out, several lines of evidence suggest that the specific pattern of allozyme (and in this case particularly PGI) variation found is caused by thermal selection. For example, genetic variation was highly locus-specific rather than relatively uniform, as should be expected for effects of natural selection. Further, the PGI 2–2 genotype dominating in alpine (in contrast to low-altitude) populations is known to exhibit increased cold stress resistance and relatively long development times typical of alpine populations. Thus, PGI (or possibly a closely linked gene) is an obvious target for thermal selection in L. tityrus . This study exemplifies how allozyme analyses can be used to detect candidate loci likely to be under selection.