Genetic evidence for kin aggregation in the intertidal acorn barnacle (Semibalanus balanoides)

It is generally assumed that larvae of benthic species are thoroughly mixed in the plankton and distributed randomly at settlement. Yet, it has also been hypothesized that a combination of larval gregarious behaviour coupled with particular oceanographic conditions may prevent larvae from mixing completely, and result in nonrandom spatial distributions following settlement. Using microsatellite markers, the main objective of this study was to investigate the occurrence of statistical connections between relatedness and settlement in the intertidal acorn barnacle from the Gulf of St Lawrence, Canada. A second objective was to test the hypothesis that patches of kin-related individuals came from a common parental site. Our results indicated that a significant number of barnacles within a given sample were more closely related than expected by chance despite the enormous potential for admixture during the planktonic phase. Thus, eight out of 37 samples analysed had relatedness values significantly higher than expected from random settlement. Moreover, analyses of sibship network construction and network complexity tests provided evidence for the occurrence of networks within several samples that were characterized by strong connections among individuals. Thus, nonrandom planktonic dispersal associated with relatively stable oceanic currents, as well as additional ecological factors to be rigorously investigated (e.g. behavioural mechanisms), may be more important in determining patterns of genetic structure in marine benthic invertebrates than generally assumed. Therefore, documenting genetic patterns associated with kin aggregation should be a fruitful and an important avenue for future studies in marine invertebrates.     

Stable genetic polymorphism in heterogeneous environments: balance between asymmetrical dispersal and selection in the acorn barnacle

Elucidating the processes responsible for maintaining polymorphism at ecologically relevant genes is intimately related to understanding the interplay between selection imposed by habitat heterogeneity and a species' capacity for dispersal in the face of environmental constraints. In this paper, we used a model-based approach to solve equilibria of balanced polymorphism, given values of fitness and larval dispersal among different habitats in the acorn barnacle Semibalanus balanoides from the Gulf of St Lawrence. Our results showed that allele frequencies observed at both MPI* and GPI* loci represented stable equilibria, given empirical estimates of fitness values, and that considerably more larvae dispersed from one region (north) to the other (south) than vice versa. Dispersal conditions were predicted to be similar for the maintenance of polymorphism at both loci. Moreover, the values of asymmetrical dispersal required by the model to reach stable equilibria were compatible with empirical estimates of larval dispersal and oceanic circulation documented in this system. Overall, this study illustrated the usefulness of a modified and computable version of Bulmer's model (1972) in order to test hypotheses of balanced polymorphism resulting from interactions between spatial selection and asymmetrical dispersal.     

Environmental heterogeneity and balancing selection in the acorn barnacle Semibalanus balanoides

The northern acorn barnacle Semibalans banlanoides occupies several intertidal microhabitats which vary greatly in their degree of physical stress. This environmental heterogeneity creates distinct selection regimes which can maintain genetic variation in natural populations. Despite considerable attention placed on the link between spatial variation in fitness and balancing selection at specific loci, experimental manipulations and fitness estimates for molecular polymorphisms have rarely been conducted in the wild. The aim of this transplant experiment was to manipulate the level of physical stress experienced by a cohort of barnacles in the field and then investigate the spatial variation in fitness for genotypes at three loci: two candidate allozymes and the mitochondrial DNA control region. The viability of mannose-6-phosphate isomerase (Mpi) genotypes was dependent on the level of physical stress experienced in the various treatments; alternative homozygotes were favoured in alternative high stress-low stress environments. In contrast, the fitness of genotypes at other loci was equivalent among treatments and unaffected by the manipulation. Evaluated in the light of balancing selection models, these data indicate that the presence of multiple environmental niches is sufficient to promote a stable Mpi polymorphism in barnacle populations and that allelic variation at this locus reflects the process of adaptation to the heterogeneous intertidal landscape.     

Fluctuating asymmetry in Semibalanus balanoides (Crustacea: Cirripedia) at different heights on the shore

Population level studies of fluctuating asymmetry (FA) need to incorporate clinal variation in FA for proper sampling methodologies to be employed to ensure appropriate comparisons between populations are made. FA of adult acorn barnacles. Semibalanus balanoides , in the intertidal zone offers an excellent model system for looking for such variation. This study looked at the magnitude of FA of S. balanoides at each of three vertical heights in the barnacle zone, in two sites on the Isle of Mull, Scotland. Unsigned FA of width of sutures between tergal and scutal plates, after correcting for size and numbers of others touched, showed significant variation across the upper, mid and lower barnacle zone. Specifically, this was a U-shaped pattern, with individuals in the mid-Zone being more symmetrical than those at the two extremes. It is proposed that high levels of FA in the upper zone are due to exposure stresses, while biotic factors such as competition impose developmental stress in the lower zone.     

Selective settlement of the barnacle Semibalanus balanoides (L.) facilitates its growth and reproduction on mussel beds in the Wadden Sea

On the unstable sedimentary tidal flats of the Wadden Sea, a suitable attachment substrate for sessile organisms is generally lacking. Epibenthic mussel beds (Mytilus edulis L.) provide the only and strongly limited settlement sites available for the barnacle, Semibalanus balanoides (L.). Field investigations showed that barnacles were non-randomly distributed within a mussel bed. They preferentially occurred near the siphonal apertures of living mussels but rarely grew on dead mussels or shell fragments. Field experiments revealed that this was due to selective settlement of barnacle cyprid larvae. Growth of barnacles was significantly higher upon living mussels than on empty mussel shells. Moreover, a higher reproductive output was obtained by individuals on living mussels which produced twice as many nauplii larvae than barnacles attached to empty shells. This study shows that selective settlement of S. balanoides cyprid larvae on living mussels is adaptive with respect to individual fitness. Received in revised form: 15 January 2001 Electronic Publication     

The role of waterholes as 'refugia' in sustaining genetic diversity and variation of two freshwater species in dryland river systems (Western Queensland, Australia)

1. Episodic floods and extended low or no flow periods characterise dryland river systems in Western Queensland, Australia. During protracted intervals between floods, rivers consist of a series of isolated waterholes, which serve as ‘refugia’ for aquatic species and much of the channel is dry. We categorised these waterholes into ‘main waterholes’, which are located in the main part of the river channel and ‘satellite waterholes’, which are located in distributary river channels. 2. We used mitochondrial sequences and allozymes to investigate levels of genetic diversity and patterns of connectivity among waterholes for two obligate freshwater species: Macrobrachium australiense (Decapoda: Palaemonidae) and Notopala sublineata (Gastropoda: Viviparidae). 3. We sampled 31 waterholes for M. australiense and 12 for N. sublineata. Based on a 505‐bp fragment of cytochrome oxidase subunit I, we identified 54 haplotypes in a sample of 232 individuals for M. australiense and based on a 457‐bp fragment of the same gene, 36 haplotypes in a sample of 145 individuals for N. sublineata. 4. Both nuclear and mitochondrial genetic data sets indicated that estimates of genetic diversity were not different in populations inhabiting main and satellite waterholes for either species. Also, there was generally very limited genetic differentiation among populations at any site. 5. We suggest that levels of connectivity among populations inhabiting waterholes at most sites are higher than expected. High levels of connectivity may help to maintain overall high levels of genetic diversity as well as low levels of genetic differentiation among waterholes within sites.     

Spatial pattern of MHC class II variation in the great snipe (Gallinago media)

The genes of the major histocompatibility complex (MHC) code for proteins involved in antigen recognition and triggering of the adaptive immune response, and are therefore likely to be under selection from parasites. These selection regimes may vary in space and time. Here we report a strong geographical structure in MHC class II B genes of a migrating bird, the great snipe (Gallinago media). Genetic differentiation in the MHC between two ecologically distinct distributional regions (Scandinavian mountain populations vs. East European lowland populations) was still present after statistically controlling for the effect of selectively neutral variation (microsatellites) using partial Mantel tests. This suggests a role for selection in generating this spatial structure and that it represents local adaptation to different environments. Differentiation between populations within the two regions was negligible. Overall, we found a high number of MHC alleles (50, from 175 individuals). This, together with a tendency for a higher rate of nonsynonymous than synonymous substitutions in the peptide binding sites, and high Tajima's D in certain regions of the gene, suggests a history of balancing selection. MHC variation is often thought to be maintained by some form of balancing selection, but the nature of this selection remains unclear. Our results support the hypothesis that spatial variation in selection regimes contributes to the high polymorphism.     

Positive Darwinian selection at the pantophysin (Pan I) locus in marine gadid fishes

Maximum-likelihood models of codon substitution were used to test for positive Darwinian selection at the vesicle protein pantophysin in two allelic lineages segregating in the Atlantic cod Gadus morhua and in 18 related species of marine gadid fishes. Positive selection was detected in the two intravesicular loops of the integral membrane protein but not in four membrane-spanning regions or the 3' cytoplasmic tail. The proportion of positively selected sites (24.9%) and the mean nonsynonymous/synonymous rate ratio (omega = d(N)/d(S) = 5.35) were both greater in the first intravesicular (IV1) domain compared with the second intravesicular (IV2) domain (11.0% positively selected sites with mean omega = 3.76). Likelihood ratio tests comparing models that assume identical omega ratios along all branches of the phylogeny to those that allow omega ratios to vary among lineages were not significant for either the IV1 or IV2 domains, indicating that the selective pressures favoring amino acid replacements have operated consistently in both regions during the diversification of the group. Positive selection was observed in the IV1 domain in both G. morhua allelic lineages, and, although three of the four codons that differ between alleles were targets of positive selection in the broader group, no similar polymorphisms were detected in other taxa. The two G. morhua Pan I alleles appeared to have evolved before the speciation event separating it from its sister taxon, Theragra chalcogramma, and on the basis of a standard mtDNA clock are estimated to be at least 2 Myr old. Although the function of pantophysin remains unknown, the strong signal of positive selection at specific sites in the IV1 and IV2 domains may help clarify its role in cellular trafficking pathways.     

Clinal genetic variation and isolation by distance in the European eel Anguilla anguilla (L.)

The genetic variability and structure of the European eel ( Anguilla anguilla L.) in populations throughout Europe was reassessed using 15 allozymic loci, seven of which were polymorphic. Seven sites were sampled on a latitudinal gradient across the natural continental range, extending from southern France to southern Norway. Heterozygosity ( H _e = 0.05) and level of polymorphism (P = 0.43) were comparable to other marine fish. Populations were poorly differentiated ( G _ST = 0.014, F _ST = 0.002), which is not surprising considering the high dispersal capability of the European eel. However, a significant geographical cline was detected at two alleles ( IDH-1 * 100 and GPI-1 * 110 ), and genetic distances ( D _CE) were concordant with geographical coastal distances. Mantel tests, pairwise F _ST's and multidimensional scaling analyses identify three distinct groups: Northern Europe, Western Europe and the Mediterranean Sea. We propose that the clinal genetic structure in the European eel may be due to (1) isolation by distance (as recently detected with microsatellites), (2) temporal reproductive separation, (3) post-larval selective forces, (4) contact between formerly separated groups or (5) some combination thereof.  © 2002 The Linnean Society of London, Biological Journal of the Linnean Society , 2002, 77 , 509–521     

On the relative roles of selection and genetic drift in shaping MHC variation

Genes of the major histocompatibility complex (MHC) have provided some of the clearest examples of how natural selection generates discordances between adaptive and neutral variation in natural populations. The type and intensity of selection as well as the strength of genetic drift are believed to be important in shaping the resulting pattern of MHC diversity. However, evaluating the relative contribution of multiple microevolutionary forces is challenging, and empirical studies have reported contrasting results. For instance, balancing selection has been invoked to explain high levels of MHC diversity and low population differentiation in comparison with other nuclear markers. Other studies have shown that genetic drift can sometimes overcome selection and then patterns of genetic variation at adaptive loci cannot be discerned from those occurring at neutral markers. Both empirical and simulated data also indicate that loss of genetic diversity at adaptive loci can occur faster than at neutral loci when selection and population bottlenecks act simultaneously. Diversifying selection, on the other hand, explains accelerated MHC divergence as the result of spatial variation in pathogen‐mediated selective regimes. Because of all these possible scenarios and outcomes, collecting information from as many study systems as possible, is crucial to enhance our understanding about the evolutionary forces driving MHC polymorphism. In this issue, Miller and co‐workers present an illuminating contribution by combining neutral markers (microsatellites) and adaptive MHC class I loci during the investigation of genetic differentiation across island populations of tuatara Sphenodon punctatus. Their study of geographical variation reveals a major role of genetic drift in shaping MHC variation, yet they also discuss some support for diversifying selection.     

Signals of heterogeneous selection at an MHC locus in geographically proximate ecotypes of sockeye salmon

The genes of the major histocompatibility complex (MHC) are an important component of the vertebrate immune system and can provide insights into the role of pathogen‐mediated selection in wild populations. Here, we examined variation at the MHC class II peptide‐binding region in 27 populations of sockeye salmon (Oncorhynchus nerka), distributed among three distinct spawning ecotypes, from a complex of interconnected rivers and lakes in south‐western Alaska. We also obtained genotypes from 90 putatively neutral single nucleotide polymorphisms for each population to compare the relative roles of demography and selection in shaping the observed MHC variation. We found that MHC divergence was generally partitioned by spawning ecotype (lake beaches, rivers and streams) and was 30 times greater than variation at neutral markers. Additionally, we observed substantial differences in modes of selection and diversity among ecotypes, with beach populations displaying higher levels of directional selection and lower MHC diversity than the other two ecotypes. Finally, the level of MHC differentiation in our study system was comparable to that observed over much larger geographic ranges, suggesting that MHC variation does not necessarily increase with increasing spatial scale and may instead be driven by fine‐scale differences in pathogen communities or pathogen virulence. The low levels of neutral structure and spatial proximity of populations in our study system indicate that MHC differentiation can be maintained through strong selective pressure even when ample opportunities for gene flow exist.     

Effects of natural selection on patterns of DNA sequence variation at the transferrin, somatolactin, and p53 genes within and among chinook salmon (Oncorhynchus tshawytscha) populations

This paper describes DNA sequence variation within and among four populations of chinook salmon (Oncorhynchus tshawytscha) at the transferrin, somatolactin and p53 genes. Patterns of variation among salmon species at the transferrin gene have been hypothesized to be shaped by positive natural selection for new alleles because the rate of nonsynonymous substitution is significantly greater than the rate of synonymous substitution. The twin goals of this study were to determine if the history of selection among salmon species at the transferrin gene is also reflected in patterns of intraspecific variation in chinook salmon, and to look for evidence of local adaptation at the transferrin gene by comparing patterns of nonsynonymous and synonymous variation among chinook salmon populations. The analyses presented here show that unlike patterns of variation between species, there is no evidence of greater differentiation among chinook salmon populations at nonsynonymous compared to synonymous sites. There is also no evidence of a reduction of within-species variation due to the hitchhiking effect at the transferrin gene, although in some populations nonsynonymous and synonymous derived mutations are both at higher frequencies than expected under a simple neutral model. Population size weighted selection coefficients (4Ns) that are consistent with both the inter and intraspecific data range from approximately 10 to approximately 235, and imply that between 1 and 40% of new nonsynonymous mutations at the transferrin gene have been beneficial.     

Allozyme variation in turbot (Psetta maxima) and brill (Scophthalmus rhombus) (Osteichthyes, Pleuronectoformes, Scophthalmidae) throughout their range in Europe

Two species of coastal flatfish (brill and turbot, Scophthalmidae) were analysed electrophoretically at 17 common enzymatic loci in samples taken from 11 sites representing the species ranges in Europe. Brill showed a mean heterozygosity ( H ) of 0·11 while that of turbot was 0·02. The virtual absence of genetic diversity in turbot is probably due to a very low evolutionary rate, and provides little evidence for population substructure even if various signs indicate the possibility of a hidden differentiation (presence of the taxon maoetica in the Black Sea and differentiation of a species-specific cestode parasite on either side of the Strait of Gibraltar). On the other hand, the weak geographic structure in brill seems to result from rapid recolonization following the last ice age.     

Mating pattern and PGI variation in a natural population of the isopod Sphaeroma rugicauda (Leach)

Samples of mating and non-mating animals from a natural population of the estuarine isopod Sphaeroma rugicauda showed no evidence of assortative mating or sexual selection for the three genotypes at the phosphoglucose isomerase locus.     

Strong genetic clines and geographical variation in gene flow in the rocky intertidal barnacle Balanus glandula

A long-standing issue in marine biology is identifying spatial scales at which populations of sessile adults are connected by planktonic offspring. We examined the genetic continuity of the acorn barnacle Balanus glandula, an abundant member of rocky intertidal communities of the northeastern Pacific Ocean, and compared these genetic patterns to the nearshore oceanography described by trajectories of surface drifters. Consistent with its broad dispersal potential, barnacle populations are genetically similar at both mitochondrial (cytochrome oxidase I) and nuclear (elongation factor 1-alpha) loci across broad swaths of the species' range. In central California, however, there is a striking genetic cline across 475 km of coastline between northern and southern populations. These patterns indicate that gene flow within central California is far more restricted spatially than among other populations. Possible reasons for the steep cline include the slow secondary introgression of historically separated populations, a balance between diversifying selection and dispersal, or some mix of both. Geographic trajectories of oceanic drifters closely parallel geographical patterns of gene flow. Drifters placed to the north (Oregon; approximately 44 degrees N) and south (Santa Barbara, California; approximately 34 degrees N) of the cline disperse hundreds of kilometers within 40 days, yet over the long-term their trajectories never overlapped. The lack of communication between waters originating in Oregon and southern California probably helps to maintain strong genetic differentiation between these regions. More broadly, the geographical variation in gene flow implies that focusing on species-level averages of gene flow can mask biologically important variance within species which reflects local environmental conditions and historical events.     

Extensive variation at MHC DRB in the New Zealand sea lion (Phocarctos hookeri) provides evidence for balancing selection

Marine mammals are often reported to possess reduced variation of major histocompatibility complex (MHC) genes compared with their terrestrial counterparts. We evaluated diversity at two MHC class II B genes, DQB and DRB, in the New Zealand sea lion (Phocarctos hookeri, NZSL) a species that has suffered high mortality owing to bacterial epizootics, using Sanger sequencing and haplotype reconstruction, together with next-generation sequencing. Despite this species'' prolonged history of small population size and highly restricted distribution, we demonstrate extensive diversity at MHC DRB with 26 alleles, whereas MHC DQB is dimorphic. We identify four DRB codons, predicted to be involved in antigen binding, that are evolving under adaptive evolution. Our data suggest diversity at DRB may be maintained by balancing selection, consistent with the role of this locus as an antigen-binding region and the species'' recent history of mass mortality during a series of bacterial epizootics. Phylogenetic analyses of DQB and DRB sequences from pinnipeds and other carnivores revealed significant allelic diversity, but little phylogenetic depth or structure among pinniped alleles; thus, we could neither confirm nor refute the possibility of trans-species polymorphism in this group. The phylogenetic pattern observed however, suggests some significant evolutionary constraint on these loci in the recent past, with the pattern consistent with that expected following an epizootic event. These data may help further elucidate some of the genetic factors underlying the unusually high susceptibility to bacterial infection of the threatened NZSL, and help us to better understand the extent and pattern of MHC diversity in pinnipeds.     

Partitioning nuclear and chloroplast variation at multiple spatial scales in the neotropical epiphytic orchid, Laelia rubescens

Insights into processes that lead to the distribution of genetic variation within plant species require recognition of the importance of both pollen and seed movement. Here we investigate the contributions of pollen and seed movement to overall gene flow in the Central American epiphytic orchid, Laelia rubescens. Genetic diversity and structure were examined at multiple spatial scales in the tropical dry forest of Costa Rica using nuclear (allozymes) and chloroplast restriction fragment length polymorphism (RFLP) markers, which were found to be diverse (allozymes, P = 73.3%; HE = 0.174; cpDNA, HE = 0.741). Nuclear genetic structure (FSTn) was low at every spatial scale (0.005-0.091). Chloroplast markers displayed more structure (0.073-0.254) but relatively similar patterns. Neither genome displayed significant isolation-by-distance. Pollen and seed dispersal rates did not differ significantly from one another (mp/ms = 1.40) at the broadest geographical scale, among sites throughout Costa Rica. However, relative contributions of pollen and seeds to gene flow were scale-dependent, with different mechanisms determining the dominant mode of gene flow at different spatial scales. Much seed dispersal is highly localized within the maternal population, while some seeds enter the air column and are dispersed over considerable distances. At the intermediate scale (10s to 100s of metres) pollinators are responsible for substantial pollen flow. This species appears capable of distributing its genes across the anthropogenically altered landscape that now characterizes its Costa Rican dry forest habitat.