Molecular genetic analysis of two native desert palm genera,Washingtonia and Brahea,from the Baja California Peninsula and Guadalupe Island

The complex geological and ecological processes that have generated high levels of biodiversity and endemism in the Baja California Peninsula have been the subject of intensive study. However, relatively little is known about phylogeography of the iconic endemic palm species of this region. We therefore analyzed a total of 2,294 bp of chloroplast and 738 bp of nuclear sequence data in 169 samples of five native palm species from Baja California, Sonora and Guadalupe Island. We found that Washingtonia and Brahea palms had low levels of genetic diversity and were highly structured, with the majority of species and major geographic regions being characterized by distinct haplotypes. We also found strong support for currently recognized species in Washingtonia, but our results were less clear cut for Brahea due to haplotype sharing. Furthermore, patterns of population structure were broadly consistent with historical vicariant events such as the inundation of the Isthmus of La Paz, the formation of the Sea of Cortez, and the more recent colonization and isolation of Guadalupe Island's palms. Our findings contribute toward a growing appreciation of the complexity of plant responses to past geological changes and also provide valuable baseline genetic data on relict American palm species.     

Genomewide patterns of variation in genetic diversity are shared among populations,species and higher‐order taxa

Genomewide screens of genetic variation within and between populations can reveal signatures of selection implicated in adaptation and speciation. Genomic regions with low genetic diversity and elevated differentiation reflective of locally reduced effective population sizes (N_e) are candidates for barrier loci contributing to population divergence. Yet, such candidate genomic regions need not arise as a result of selection promoting adaptation or advancing reproductive isolation. Linked selection unrelated to lineage‐specific adaptation or population divergence can generate comparable signatures. It is challenging to distinguish between these processes, particularly when diverging populations share ancestral genetic variation. In this study, we took a comparative approach using population assemblages from distant clades assessing genomic parallelism of variation in N_e. Utilizing population‐level polymorphism data from 444 resequenced genomes of three avian clades spanning 50 million years of evolution, we tested whether population genetic summary statistics reflecting genomewide variation in N_e would covary among populations within clades, and importantly, also among clades where lineage sorting has been completed. All statistics including population‐scaled recombination rate (ρ), nucleotide diversity (π) and measures of genetic differentiation between populations (F_ST, PBS, d_xy) were significantly correlated across all phylogenetic distances. Moreover, genomic regions with elevated levels of genetic differentiation were associated with inferred pericentromeric and subtelomeric regions. The phylogenetic stability of diversity landscapes and stable association with genomic features support a role of linked selection not necessarily associated with adaptation and speciation in shaping patterns of genomewide heterogeneity in genetic diversity.     

QTL mapping and GWAS reveal candidate genes controlling capsaicinoid content in Capsicum

Capsaicinoids are unique compounds produced only in peppers (Capsicum spp.). Several studies using classical quantitative trait loci (QTLs) mapping and genomewide association studies (GWAS) have identified QTLs controlling capsaicinoid content in peppers; however, neither the QTLs common to each population nor the candidate genes underlying them have been identified due to the limitations of each approach used. Here, we performed QTL mapping and GWAS for capsaicinoid content in peppers using two recombinant inbred line (RIL) populations and one GWAS population. Whole‐genome resequencing and genotyping by sequencing (GBS) were used to construct high‐density single nucleotide polymorphism (SNP) maps. Five QTL regions on chromosomes 1, 2, 3, 4 and 10 were commonly identified in both RIL populations over multiple locations and years. Furthermore, a total of 109 610 SNPs derived from two GBS libraries were used to analyse the GWAS population consisting of 208 C. annuum‐clade accessions. A total of 69 QTL regions were identified from the GWAS, 10 of which were co‐located with the QTLs identified from the two biparental populations. Within these regions, we were able to identify five candidate genes known to be involved in capsaicinoid biosynthesis. Our results demonstrate that QTL mapping and GBS‐GWAS represent a powerful combined approach for the identification of loci controlling complex traits.     

Genetic divergence with ongoing gene flow is maintained by the use of different hosts in phytophagous ladybird beetles genus Henosepilachna

Adaptation to different environments can promote population divergence via natural selection even in the presence of gene flow – a phenomenon that typically occurs during ecological speciation. To elucidate how natural selection promotes and maintains population divergence during speciation, we investigated the population genetic structure, degree of gene flow and heterogeneous genomic divergence in three closely related Japanese phytophagous ladybird beetles: Henosepilachna pustulosa, H. niponica and H. yasutomii. These species act as a generalist, a wild thistle (Cirsium spp.) specialist and a blue cohosh (Caulophyllum robustum) specialist, respectively, and their ranges differ accordingly. The two specialist species widely co‐occur but are reproductively isolated solely due to their high specialization to a particular host plant. Genomewide amplified fragment‐length polymorphism (AFLP) markers and mitochondrial cytochrome c oxidase subunit I (COI) gene sequences demonstrated obvious genomewide divergence associated with both geographic distance and ecological divergence. However, a hybridization assessment for both AFLP loci and the mitochondrial sequences revealed a certain degree of unidirectional gene flow between the two sympatric specialist species. Principal coordinates analysis (PCoA) based on all of the variable AFLP loci demonstrated that there are genetic similarities between populations from adjacent localities irrespective of the species (i.e. host range). However, a further comparative genome scan identified a few fractions of loci representing approximately 1% of all loci as different host‐associated outliers. These results suggest that these three species had a complex origin, which could be obscured by current gene flow, and that ecological divergence can be maintained with only a small fraction of the genome is related to different host use even when there is a certain degree of gene flow between sympatric species pairs.     

Population structure and local selection yield high genomic variation in Mimulus guttatus

Across western North America, Mimulus guttatus exists as many local populations adapted to site‐specific environmental challenges. Gene flow between locally adapted populations will affect genetic diversity both within demes and across the larger metapopulation. Here, we analyse 34 whole‐genome sequences from the intensively studied Iron Mountain population (IM) in conjunction with sequences from 22 Mimulus individuals sampled from across western North America. Three striking features of these data address hypotheses about migration and selection in a locally adapted population. First, we find very high levels of intrapopulation polymorphism (synonymous π = 0.033). Variation outside of genes is likely even higher but difficult to estimate because excessive divergence reduces the efficiency of read mapping. Second, IM exhibits a significantly positive genomewide average for Tajima's D. This indicates allele frequencies are typically more intermediate than expected from neutrality, opposite the pattern observed in many other species. Third, IM exhibits a distinctive haplotype structure with a genomewide excess of positive associations between rarer alleles at linked loci. This suggests an important effect of gene flow from other Mimulus populations, although a residual effect of population founding might also contribute. The combination of multiple analyses, including a novel tree‐based analytic method, illustrates how the balance of local selection, limited dispersal and metapopulation dynamics manifests across the genome. The overall genomic pattern of sequence diversity suggests successful gene flow of divergent immigrant genotypes into IM. However, many loci show patterns indicative of local adaptation, particularly at SNPs associated with chromosomal inversions.     

SNP discovery in wild and domesticated populations of blue catfish,Ictalurus furcatus,using genotyping‐by‐sequencing and subsequent SNP validation

Blue catfish, Ictalurus furcatus, are valued in the United States as a trophy fishery for their capacity to reach large sizes, sometimes exceeding 45 kg. Additionally, blue catfish × channel catfish (I. punctatus) hybrid food fish production has recently increased the demand for blue catfish broodstock. However, there has been little study of the genetic impacts and interaction of farmed, introduced and stocked populations of blue catfish. We utilized genotyping‐by‐sequencing (GBS) to capture and genotype SNP markers on 190 individuals from five wild and domesticated populations (Mississippi River, Missouri, D&B, Rio Grande and Texas). Stringent filtering of SNP‐calling parameters resulted in 4275 SNP loci represented across all five populations. Population genetics and structure analyses revealed potential shared ancestry and admixture between populations. We utilized the Sequenom MassARRAY to validate two multiplex panels of SNPs selected from the GBS data. Selection criteria included SNPs shared between populations, SNPs specific to populations, number of reads per individual and number of individuals genotyped by GBS. Putative SNPs were validated in the discovery population and in two additional populations not used in the GBS analysis. A total of 64 SNPs were genotyped successfully in 191 individuals from nine populations. Our results should guide the development of highly informative, flexible genotyping multiplexes for blue catfish from the larger GBS SNP set as well as provide an example of a rapid, low‐cost approach to generate and genotype informative marker loci in aquatic species with minimal previous genetic information.     

Population genetic structure of Grauer's Swamp Warbler Bradypterus graueri,an Albertine Rift endemic

The endangered warbler Bradypterus graueri is endemic to the Albertine Rift, where it is restricted to montane swamps above 1900 m across the region. We studied genetic structure among six populations sampled across the species' distribution in northern Burundi, Rwanda, Uganda and the eastern Democratic Republic of Congo. A total of 2117 base pairs of mitochondrial data were sequenced. Phylogenetic analyses and network reconstruction of B. graueri haplotypes recovered three clades with a defined geographical pattern: clade 1, Virunga Volcanoes and Kigezi Highlands; clade 2, Rugege Highlands; and clade 3, Kahuzi‐Biega Highlands; clades 2 and 3 are sisters to each other. Both landscape dynamics and historical climate are likely to have played a role in the diversification of this species. The divergence between clade 1 and clades 2 and 3 (168.5 ka, 95% HPD 108.5, 244.4) coincides with a prolonged period of aridity in tropical Africa between 130 and 270 ka. Similarly, the divergence between clades 2 and 3 (99.4 ka, 95% HPD 55.4, 153.8) corresponds with a period of aridity just prior to 94 ka. Populations sampled from the eastern arm of the central Albertine Rift (Kigezi and Rugege Highlands) show a coincident increase in effective population size after the Last Glacial Maximum at c. 15 ka, whereas those sampled from Kahuzi‐Biega on the western arm of the rift do not. Despite the perceived higher vagility of bird species relative to other vertebrates, the degree of phylogeographical structure among populations of B. graueri is similar to that reported for small mammals (Hylomyscus vulcanorum, Lophuromys woosnami, Sylvisorex vulcanorum) and a frog Hyperolius castaneus sampled across the central Albertine Rift. Collectively our results suggest that climate dynamics associated with late Pleistocene cycles had a significant influence on driving the population genetic structure and associated levels of genetic diversity in B. graueri and other small terrestrial vertebrates. Our results have implications for the conservation of B. graueri and other endemics to the Albertine Rift, particularly in the context of other phylogegeographical studies centred on this biodiversity hotspot.     

Comparative phylogeography of mainland and insular species of Neotropical molossid bats (Molossus)

Historical events, habitat preferences, and geographic barriers might result in distinct genetic patterns in insular versus mainland populations. Comparison between these two biogeographic systems provides an opportunity to investigate the relative role of isolation in phylogeographic patterns and to elucidate the importance of evolution and demographic history in population structure. Herein, we use a genotype‐by‐sequencing approach (GBS) to explore population structure within three species of mastiff bats (Molossus molossus, M. coibensis, and M. milleri), which represent different ecological histories and geographical distributions in the genus. We tested the hypotheses that oceanic straits serve as barriers to dispersal in Caribbean bats and that isolated island populations are more likely to experience genetic drift and bottlenecks in comparison with highly connected ones, thus leading to different phylogeographic patterns. We show that population structures vary according to general habitat preferences, levels of population isolation, and historical fluctuations in climate. In our dataset, mainland geographic barriers played only a small role in isolation of lineages. However, oceanic straits posed a partial barrier to the dispersal for some populations within some species (M. milleri), but do not seem to disrupt gene flow in others (M. molossus). Lineages on distant islands undergo genetic bottlenecks more frequently than island lineages closer to the mainland, which have a greater exchange of haplotypes.     

History,geography and host use shape genomewide patterns of genetic variation in the redheaded pine sawfly (Neodiprion lecontei)

Divergent host use has long been suspected to drive population differentiation and speciation in plant‐feeding insects. Evaluating the contribution of divergent host use to genetic differentiation can be difficult, however, as dispersal limitation and population structure may also influence patterns of genetic variation. In this study, we use double‐digest restriction‐associated DNA (ddRAD) sequencing to test the hypothesis that divergent host use contributes to genetic differentiation among populations of the redheaded pine sawfly (Neodiprion lecontei), a widespread pest that uses multiple Pinus hosts throughout its range in eastern North America. Because this species has a broad range and specializes on host plants known to have migrated extensively during the Pleistocene, we first assess overall genetic structure using model‐based and model‐free clustering methods and identify three geographically distinct genetic clusters. Next, using a composite‐likelihood approach based on the site frequency spectrum and a novel strategy for maximizing the utility of linked RAD markers, we infer the population topology and date divergence to the Pleistocene. Based on existing knowledge of Pinus refugia, estimated demographic parameters and patterns of diversity among sawfly populations, we propose a Pleistocene divergence scenario for N. lecontei. Finally, using Mantel and partial Mantel tests, we identify a significant relationship between genetic distance and geography in all clusters, and between genetic distance and host use in two of three clusters. Overall, our results indicate that Pleistocene isolation, dispersal limitation and ecological divergence all contribute to genomewide differentiation in this species and support the hypothesis that host use is a common driver of population divergence in host‐specialized insects.     

Outlier SNP markers reveal fine‐scale genetic structuring across European hake populations (Merluccius merluccius)

Shallow population structure is generally reported for most marine fish and explained as a consequence of high dispersal, connectivity and large population size. Targeted gene analyses and more recently genome‐wide studies have challenged such view, suggesting that adaptive divergence might occur even when neutral markers provide genetic homogeneity across populations. Here, 381 SNPs located in transcribed regions were used to assess large‐ and fine‐scale population structure in the European hake (Merluccius merluccius), a widely distributed demersal species of high priority for the European fishery. Analysis of 850 individuals from 19 locations across the entire distribution range showed evidence for several outlier loci, with significantly higher resolving power. While 299 putatively neutral SNPs confirmed the genetic break between basins (F_CT = 0.016) and weak differentiation within basins, outlier loci revealed a dramatic divergence between Atlantic and Mediterranean populations (F_CT range 0.275–0.705) and fine‐scale significant population structure. Outlier loci separated North Sea and Northern Portugal populations from all other Atlantic samples and revealed a strong differentiation among Western, Central and Eastern Mediterranean geographical samples. Significant correlation of allele frequencies at outlier loci with seawater surface temperature and salinity supported the hypothesis that populations might be adapted to local conditions. Such evidence highlights the importance of integrating information from neutral and adaptive evolutionary patterns towards a better assessment of genetic diversity. Accordingly, the generated outlier SNP data could be used for tackling illegal practices in hake fishing and commercialization as well as to develop explicit spatial models for defining management units and stock boundaries.     

Are habitat fragmentation,local adaptation and isolation‐by‐distance driving population divergence in wild rice Oryza rufipogon?

Habitat fragmentation weakens the connection between populations and is accompanied with isolation by distance (IBD) and local adaptation (isolation by adaptation, IBA), both leading to genetic divergence between populations. To understand the evolutionary potential of a population and to formulate proper conservation strategies, information on the roles of IBD and IBA in driving population divergence is critical. The putative ancestor of Asian cultivated rice (Oryza sativa) is endangered in China due to habitat loss and fragmentation. We investigated the genetic variation in 11 Chinese Oryza rufipogon populations using 79 microsatellite loci to infer the effects of habitat fragmentation, IBD and IBA on genetic structure. Historical and current gene flows were found to be rare (m_h = 0.0002–0.0013, m_c = 0.007–0.029), indicating IBD and resulting in a high level of population divergence (F_ST = 0.343). High within‐population genetic variation (H_E = 0.377–0.515), relatively large effective population sizes (N_e = 96–158), absence of bottlenecks and limited gene flow were found, demonstrating little impact of recent habitat fragmentation on these populations. Eleven gene‐linked microsatellite loci were identified as outliers, indicating local adaptation. Hierarchical AMOVA and partial Mantel tests indicated that population divergence of Chinese O. rufipogon was significantly correlated with environmental factors, especially habitat temperature. Common garden trials detected a significant adaptive population divergence associated with latitude. Collectively, these findings imply that IBD due to historical rather than recent fragmentation, followed by local adaptation, has driven population divergence in O. rufipogon.     

Isolation and characterization of eight polymorphic microsatellite loci in the rainforest canopy tree,Syzygium sayeri (Myrtaceae)

Eight polymorphic microsatellite loci were isolated from the rainforest canopy tree, Syzygium sayeri, in order to study parentage and subsequently pollen dispersal among individuals in wild populations. Screening of one natural population (n = 64) mapped in a 500 × 500 m area at Cape Tribulation, north Queensland, Australia, yielded two to 11 alleles per locus with observed levels of heterozygosity ranging from 0.07 to 0.70. One locus was significantly out of Hardy–Weinberg equilibrium after Bonferroni correction for multiple tests. These loci should provide a useful tool in further understanding the dispersal patterns of this species.     

Spatial genetic structure in Beta vulgaris subsp. maritima and Beta macrocarpa reveals the effect of contrasting mating system,influence of marine currents,and footprints of postglacial recolonization routes

Understanding the factors that contribute to population genetic divergence across a species' range is a long‐standing goal in evolutionary biology and ecological genetics. We examined the relative importance of historical and ecological features in shaping the present‐day spatial patterns of genetic structure in two related plant species, Beta vulgaris subsp. maritima and Beta macrocarpa. Using nuclear and mitochondrial markers, we surveyed 93 populations from Brittany (France) to Morocco – the southern limit of their species' range distribution. Whereas B. macrocarpa showed a genotypic structure and a high level of genetic differentiation indicative of selfing, the population genetic structure of B. vulgaris subsp. maritima was consistent with an outcrossing mating system. We further showed (1) a strong geographic clustering in coastal B. vulgaris subsp. maritima populations that highlighted the influence of marine currents in shaping different lineages and (2) a peculiar genetic structure of inland B. vulgaris subsp. maritima populations that could indicate the admixture of distinct evolutionary lineages and recent expansions associated with anthropogenic disturbances. Spatial patterns of nuclear diversity and differentiation also supported a stepwise recolonization of Europe from Atlantic‐Mediterranean refugia after the last glacial period, with leading‐edge expansions. However, cytoplasmic diversity was not impacted by postglacial recolonization: stochastic long‐distance seed dispersal mediated by major oceanic currents may mitigate the common patterns of reduced cytoplasmic diversity observed for edge populations. Overall, the patterns we documented here challenge the general view of reduced genetic diversity at the edge of a species' range distribution and provide clues for understanding how life‐history and major geographic features interact to shape the distribution of genetic diversity.     

Complex longitudinal diversification across South China and Vietnam in Stejneger's pit viper,Viridovipera stejnegeri (Schmidt, 1925) (Reptilia: Serpentes: Viperidae)

Viridovipera stejnegeri is one of the most common pit vipers in Asia, with a wide distribution in southern China and Vietnam. We investigated historical demography and explored how the environment and climatic factors have shaped genetic diversity and the evolutionary history of this venomous snake. A total of 171 samples from 47 localities were sequenced and analysed for two mitochondrial gene fragments and three nuclear genes. Gene trees reveal the existence of two well‐supported clades (Southwest China and Southeast China) with seven distinct and strongly supported, geographically structured subclades within V. stejnegeri. Estimation of divergence time and ancestral area suggests that V. stejnegeri originated at ~6.0 Ma in the late Miocene on the Yunnan–Guizhou Plateau. The estimated date of origin and divergence of the island populations of Taiwan and Hainan closely matches the geological origin of the both islands. The mtDNA gene tree reveals the presence of west–east diversification in V. stejnegeri populations. Complex orogenesis and heterogeneous habitats, as well as climate‐mediated habitat differentiation including glacial cycles, all have influenced population structure and the distribution of this taxon. The validity of V. stejnegeri chenbihuii is questionable, and this subspecies most probably represents an invalid taxon.     

Genetic variation within and between populations of hermaphroditic <Emphasis Type="Italic">Bulinus truncatus</Emphasis> tetraploid freshwater snails of the Albertine Rift,East Africa

Genetic variability within and among Bulinus truncatus of the Albertine Rift freshwater bodies were assessed to investigate the degree of inbreeding and gene flow in the snail populations. The effect of ploidy on the genetic structuring of B. truncatus is also described. We characterized the genetic structure of seven B. truncatus populations from Lake Albert, Lake Kivu, and Katosho swamp in Tanzania using five polymorphic microsatellite loci. Genetic differentiation was quantified using pairwise FST values and Nei’s standard genetic distances. Different alleles were observed across all loci and genetic diversity was low although it varied greatly across populations; observed heterozygosity was, however, higher than the expected heterozygosity in three of the populations studied. Significant heterozygote deficiencies were observed coupled with significant linkage disequilibria in five populations for all the five loci examined in this study. We found significant genetic differentiation among the seven freshwater bodies; private alleles were observed across all loci indicating restricted or absence of gene flow between populations. Limited snail dispersal and the reproductive biology of B. truncatus are the major forces shaping the genetic variation observed. Low genetic variation within B. truncatus populations exposes them to a high parasite infection risk as predicted in the Red Queen hypothesis.     

Microsatellite multiplexes for high‐throughput genotyping of French grunts (Haemulon flavolineatum,Pisces: Haemulidae) and their utility in other grunt species

Grunts of the genus Haemulon are abundant in the western Atlantic and are an important component of regional fisheries. Management efforts for these and other Caribbean reef fish require an understanding of population structure and the extent of dispersal among populations. We characterized 11 polymorphic microsatellite loci for French grunts (Haemulon flavolineatum) and grouped them into three PCR multiplexes to facilitate high‐throughput genotyping. These loci will be useful for population genetic studies of French grunts. Cross‐species amplifications suggest these loci are also suitable for population genetic studies of H. aurolineatum, H. plumieri, and H. sciurus.     

Phylogeography and the origins of range disjunctions in a north temperate fish,the pygmy whitefish (Prosopium coulterii), inferred from mitochondrial and nuclear DNA sequence analysis

Aim To investigate the degree of phylogeographical divergence within pygmy whitefish (Prosopium coulterii) and to test hypotheses concerning the origin of disjunct populations within North America. Location North America from western Alaska to Lake Superior. Methods Mitochondrial (ATPase subunit VI) and nuclear (ITS‐1, ITS‐2) DNA sequence variation was assessed across the species’ North American range to test for the existence of distinct phylogeographical groupings of pygmy whitefish associated with known glacial refugia. Coalescent simulations of the mitochondrial DNA (mtDNA) data were used to test hypotheses of population structure. Results This species is composed of two monophyletic mitochondrial clades across its North American range. The two mtDNA clades differed by an average 3.3% nucleotide sequence divergence. These clades were also distinguished by ITS‐2, but the relationships among lineages were not resolved by the ITS‐1 analysis. Coalescent analyses rejected the null hypothesis of the current disjunct distributions being a result of fragmentation of a single widespread ancestral lineage across a variety of effective population sizes and divergence times. Main conclusions The current range disjunctions of pygmy whitefish in North America probably resulted from isolation, genetic divergence, and selective dispersal from at least two major Pleistocene glacial refugia: Beringia and Cascadia. More recent isolation and dispersal from an upper Mississippi refugium is suggested by relationships within one of the clades and by distributional evidence from co‐distributed species. The Beringian and Cascadian refugia have played major roles in the zoogeography of Nearctic temperate aquatics, but the roles of smaller refugia appear more variable among other species.     

Elevation as a barrier: genetic structure for an Atlantic rain forest tree (Bathysa australis) in the Serra do Mar mountain range,SE Brazil

Distance and discrete geographic barriers play a role in isolating populations, as seed and pollen dispersal become limited. Nearby populations without any geographic barrier between them may also suffer from ecological isolation driven by habitat heterogeneity, which may promote divergence by local adaptation and drift. Likewise, elevation gradients may influence the genetic structure and diversity of populations, particularly those marginally distributed. Bathysa australis (Rubiaceae) is a widespread tree along the elevation gradient of the Serra do Mar, SE Brazil. This self‐compatible species is pollinated by bees and wasps and has autochoric seeds, suggesting restricted gene dispersal. We investigated the distribution of genetic diversity in six B. australis populations at two extreme sites along an elevation gradient: a lowland site (80–216 m) and an upland site (1010–1100 m.a.s.l.). Nine microsatellite loci were used to test for genetic structure and to verify differences in genetic diversity between sites. We found a marked genetic structure on a scale as small as 6 km (F_ST = 0.21), and two distinct clusters were identified, each corresponding to a site. Although B. australis is continuously distributed along the elevation gradient, we have not observed a gene flow between the extreme populations. This might be related to B. australis biological features and creates a potential scenario for adaptation to the different conditions imposed by the elevation gradient. We failed to find an isolation‐by‐distance pattern; although on the fine scale, all populations showed spatial autocorrelation until ~10‐20 m. Elevation difference was a relevant factor though, but we need further sampling effort to check its correlation with genetic distance. The lowland populations had a higher allelic richness and showed higher rare allele counts than the upland ones. The upland site may be more selective, eliminating rare alleles, as we did not find any evidence for bottleneck.     

Fin whale MDH‐1 and MPI allozyme variation is not reflected in the corresponding DNA sequences

The appeal of genetic inference methods to assess population genetic structure and guide management efforts is grounded in the correlation between the genetic similarity and gene flow among populations. Effects of such gene flow are typically genomewide; however, some loci may appear as outliers, displaying above or below average genetic divergence relative to the genomewide level. Above average population, genetic divergence may be due to divergent selection as a result of local adaptation. Consequently, substantial efforts have been directed toward such outlying loci in order to identify traits subject to local adaptation. Here, we report the results of an investigation into the molecular basis of the substantial degree of genetic divergence previously reported at allozyme loci among North Atlantic fin whale (Balaenoptera physalus) populations. We sequenced the exons encoding for the two most divergent allozyme loci (MDH‐1 and MPI) and failed to detect any nonsynonymous substitutions. Following extensive error checking and analysis of additional bioinformatic and morphological data, we hypothesize that the observed allozyme polymorphisms may reflect phenotypic plasticity at the cellular level, perhaps as a response to nutritional stress. While such plasticity is intriguing in itself, and of fundamental evolutionary interest, our key finding is that the observed allozyme variation does not appear to be a result of genetic drift, migration, or selection on the MDH‐1 and MPI exons themselves, stressing the importance of interpreting allozyme data with caution. As for North Atlantic fin whale population structure, our findings support the low levels of differentiation found in previous analyses of DNA nucleotide loci.