Concordant mitochondrial and nuclear DNA phylogenies for populations of the teleost fish Fundulus heteroclitus

Mitochondrial Ecophysiology: Assessing the Evolutionary Forces That Shape Mitochondrial Variation

The mitonuclear species concept hypothesizes that incompatibilities between interacting gene products of the nuclear and mitochondrial genomes are a major factor establishing and maintaining species boundaries. However, most of the data available to test this concept come from studies of genetic variation in mitochondrial DNA, and clines in the mitochondrial genome across contact zones can be produced by a variety of forces. Here, we show that using a combination of population genomic analyses of the nuclear and mitochondrial genomes and studies of mitochondrial function can provide insight into the relative roles of neutral processes, adaptive evolution, and mitonuclear incompatibility in establishing and maintaining mitochondrial clines, using Atlantic killifish (Fundulus heteroclitus) as a case study. There is strong evidence for a role of secondary contact following the last glaciation in shaping a steep mitochondrial cline across a contact zone between northern and southern subspecies of killifish, but there is also evidence for a role of adaptive evolution in driving differentiation between the subspecies in a variety of traits from the level of the whole organism to the level of mitochondrial function. In addition, studies are beginning to address the potential for mitonuclear incompatibilities in admixed populations. However, population genomic studies have failed to detect evidence for a strong and pervasive influence of mitonuclear incompatibilities, and we suggest that polygenic selection may be responsible for the complex patterns observed. This case study demonstrates that multiple forces can act together in shaping mitochondrial clines, and illustrates the challenge of disentangling their relative roles.

Invasion genetics of the mummichog (Fundulus heteroclitus): recent anthropogenic introduction in Iberia

Human activities such as trade and transport have increased considerably in the last decades, greatly facilitating the introduction and spread of non-native species at a global level. In the Iberian Peninsula, Fundulus heteroclitus, a small euryhaline coastal fish with short dispersal, was found for the first time in the mid-1970s. Since then, F. heteroclitus has undergone range expansions, colonizing the southern region of Portugal, southwestern coast of Spain and the Ebro Delta in the Mediterranean Sea. Cytochrome b sequences were used to elucidate the species invasion pathway in Iberia. Three Iberian locations (Faro, Cádiz and Ebro Delta) and 13 other locations along the native range of F. heteroclitus in North America were sampled. Results revealed a single haplotype, common to all invasive populations, which can be traced to the northern region of the species’ native range. We posit that the origin of the founder individuals is between New York and Nova Scotia. Additionally, the lack of genetic structure within Iberia is consistent with a recent invasion scenario and a strong founder effect. We suggest the most probable introduction vector is associated with the aquarium trade. We further discuss the hypothesis of a second human-mediated introduction responsible for the establishment of individuals in the Ebro Delta supported by the absence of adequate muddy habitats linking Cádiz and the Ebro Delta. Although the species has a high tolerance to salinity and temperature, ecological niche modelling indicates that benthic habitat constraints prevent along-shore colonisation suggesting that such expansions would need to be aided by human release.

Genomic islands of divergence in the Yellow Tang and the Brushtail Tang Surgeonfishes

The current ease of obtaining thousands of molecular markers challenges the notion that full phylogenetic concordance, as proposed by phylogenetic species concepts, is a requirement for defining species delimitations. Indeed, the presence of genomic islands of divergence, which may be the cause, or in some cases the consequence, of speciation, precludes concordance. Here, we explore this issue using thousands of RAD markers on two sister species of surgeonfishes (Teleostei: Acanthuridae), Zebrasoma flavescens and Z. scopas, and several populations within each species. Species are readily distinguished based on their colors (solid yellow and solid brown, respectively), yet populations and species are neither distinguishable using mitochondrial markers (cytochrome c oxidase 1), nor using 5193 SNPs (pairwise Φst = 0.034). In contrast, when using outlier loci, some of them presumably under selection, species delimitations, and strong population structure follow recognized taxonomic positions (pairwise Φst = 0.326). Species and population delimitation differences based on neutral and selected markers are likely due to local adaptation, thus being consistent with the idea that these genomic islands of divergence arose as a consequence of isolation. These findings, which are not unique, raise the question of a potentially important pathway of divergence based on local adaptation that is only evident when looking at thousands of loci.

Steep, coincident, and concordant clines in mitochondrial and nuclear-encoded genes in a hybrid zone between subspecies of Atlantic killifish, Fundulus heteroclitus

Steep genetic clines resulting from recent secondary contact between previously isolated taxa can either gradually erode over time or be stabilized by factors such as ecological selection or selection against hybrids. We used patterns of variation in 30 nuclear and two mitochondrial SNPs to examine the factors that could be involved in stabilizing clines across a hybrid zone between two subspecies of the Atlantic killifish, Fundulus heteroclitus. Increased heterozygote deficit and cytonuclear disequilibrium in populations near the center of the mtDNA cline suggest that some form of reproductive isolation such as assortative mating or selection against hybrids may be acting in this hybrid zone. However, only a small number of loci exhibited these signatures, suggesting locus-specific, rather than genomewide, factors. Fourteen of the 32 loci surveyed had cline widths inconsistent with neutral expectations, with two SNPs in the mitochondrial genome exhibiting the steepest clines. Seven of the 12 putatively non-neutral nuclear clines were for SNPs in genes related to oxidative metabolism. Among these putatively non-neutral nuclear clines, SNPs in two nuclear-encoded mitochondrial genes (SLC25A3 and HDDC2), as well as SNPs in the myoglobin, 40S ribosomal protein S17, and actin-binding LIM protein genes, had clines that were coincident and concordant with the mitochondrial clines. When hybrid index was calculated using this subset of loci, the frequency distribution of hybrid indices for a population located at the mtDNA cline center was non-unimodal, suggesting selection against advanced-generation hybrids, possibly due to effects on processes involved in oxidative metabolism.

Molecular systematics of the freshwater stingrays (myliobatiformes: potamotrygonidae) of the Amazon, Orinoco, Magdalena, Esequibo, Caribbean, and Maracaibo basins (Colombia - Venezuela): evidence from three mitochondrial genes

Lack of adequate information about the taxonomic and evolutionary relationships, ecology, biology, and distribution of several species belonging to the family Potamotrygonidae makes these species vulnerable to anthropic activities, including commercial overexploitation for the ornamental fish market. The aim of this study was to investigate the systematic relationships among genera and species belonging to this family by analyses of three mitochondrial gene regions. Samples were collected from the main river basins in Colombia and Venezuela for four genera and seven species of the family, as well as for what appear to be unidentified species. Three mitochondrial molecular markers COI, Cytb, and ATP6 were amplified and sequenced. Maximum likelihood and Bayesian inference analysis were performed to obtain topologies for each marker and for a concatenated dataset including the three genes. Small dataset may compromise some methods estimations of sequence divergence in the ATP6 marker. Monophyly of the four genera in Potamotrygonidae was confirmed and phylogenetic relationships among members of the Potamotrygon genus were not clearly resolved. However, results obtained with the molecular marker Cytb appear to offer a good starting point to differentiate among genera and species as a tool that could be used for barcoding. The application of this gene as a barcode could be applied for management and regulation of extraction practices for these genera. Sequencing complete mitochondrial genomes would be the next step for testing evolutionary hypothesis among these genera. Population structure analyses should be undertaken for Paratrygon, Potamotrygon magdalenae and motoro.

Evidence for a bimodal distribution of hybrid indices in a hybrid zone with high admixture

The genetic structure of a hybrid zone can provide insights into the relative roles of the various factors that maintain the zone. Here, we use a multilocus approach to characterize a hybrid zone between two subspecies of killifish (Fundulus heteroclitus, Walbaum 1792) found along the Atlantic coast of North America. We first analysed clinal variation along the Atlantic coast using a single-nucleotide polymorphism in the mitochondrial DNA (mtDNA) displacement loop (D-loop) and a panel of nine nuclear microsatellite markers. A model constraining all clines to the same width and centre was not significantly different from a model in which the clines were allowed to vary independently. Locus-by-locus analysis indicated that the majority of nuclear clines shared the same centre as the mtDNA cline, and the widths of these clines were also narrower than that predicted by a neutral model, suggesting that selection is operating to maintain the hybrid zone. However, two of the nuclear clines had widths greater than the neutral prediction and had centres that were displaced relative to the mtDNA cline centre. We also found that a marsh located near the centre of the mtDNA cline demonstrated a bimodal distribution of nuclear hybrid index values, suggesting a deficit of first-generation hybrids and backcrossed genotypes. Thus, selection against hybrid genotypes may be playing a role in maintaining this hybrid zone and the associated steep nuclear and mtDNA clines.

Nucleotide variation and balancing selection at the Ckma gene in Atlantic cod: analysis with multiple merger coalescent models

High-fecundity organisms, such as Atlantic cod, can withstand substantial natural selection and the entailing genetic load of replacing alleles at a number of loci due to their excess reproductive capacity. High-fecundity organisms may reproduce by sweepstakes leading to highly skewed heavy-tailed offspring distribution. Under such reproduction the Kingman coalescent of binary mergers breaks down and models of multiple merger coalescent are more appropriate. Here we study nucleotide variation at the Ckma (Creatine Kinase Muscle type A) gene in Atlantic cod. The gene shows extreme differentiation between the North (Canada, Greenland, Iceland, Norway, Barents Sea) and the South (Faroe Islands, North-, Baltic-, Celtic-, and Irish Seas) with F_ST > 0.8 between regions whereas neutral loci show no differentiation. This is evidence of natural selection. The protein sequence is conserved by purifying selection whereas silent and non-coding sites show extreme differentiation. The unfolded site-frequency spectrum has three modes, a mode at singleton sites and two high frequency modes at opposite frequencies representing divergent branches of the gene genealogy that is evidence for balancing selection. Analysis with multiple-merger coalescent models can account for the high frequency of singleton sites and indicate reproductive sweepstakes. Coalescent time scales vary with population size and with the inverse of variance in offspring number. Parameter estimates using multiple-merger coalescent models show that times scales are faster than under the Kingman coalescent.

The Elizabeth River Story: A Case Study in Evolutionary Toxicology

The Elizabeth River system is an estuary in southeastern Virginia, surrounded by the towns of Chesapeake, Norfolk, Portsmouth, and Virginia Beach. The river has played important roles in U.S. history and has been the location of various military and industrial activities. These activities have been the source of chemical contamination in this aquatic system. Important industries, until the 1990s, included wood treatment plants that used creosote, an oil-derived product that is rich in polycyclic aromatic hydrocarbons (PAH). These plants left a legacy of PAH pollution in the river, and in particular Atlantic Wood Industries is a designated Superfund site now undergoing remediation. Numerous studies examined the distribution of PAH in the river and impacts on resident fauna. This review focuses on how a small estuarine fish with a limited home range, Fundulus heteroclitus (Atlantic killifish or mummichog), has responded to this pollution. While in certain areas of the river this species has clearly been impacted, as evidenced by elevated rates of liver cancer, some subpopulations, notably the one associated with the Atlantic Wood Industries site, displayed a remarkable ability to resist the marked effects PAH have on the embryonic development of fish. This review provides evidence of how pollutants have acted as evolutionary agents, causing changes in ecosystems potentially lasting longer than the pollutants themselves. Mechanisms underlying this evolved resistance, as well as mechanisms underlying the effects of PAH on embryonic development, are also described. The review concludes with a description of ongoing and promising efforts to restore this historic American river.

Interindividual variation in Complex I activity in Fundulus heteroclitus along a steep thermocline

The first enzyme in the oxidative phosphorylation pathway is Complex I (EC 1.6.5.3). Complex I is a large heteromeric enzyme complex with 45 protein subunits that translocates H(+) ions across the mitochondrial inner membrane. Among northern and southern populations of the teleost fish Fundulus heteroclitus, Complex I subunits have fixed amino acid substitutions. Additionally, there are differences in oxidative phosphorylation activity among populations of F. heteroclitus. To investigate whether these differences are related to Complex I, enzyme activity was measured in 121 individuals from five populations of F. heteroclitus and its sister species Fundulus grandis acclimated to a constant 20°C temperature. Within each population, Complex I activity is highly variable among individuals of F. heteroclitus (coefficient of variation percentage among individuals has a mean of 90% in the five F. heteroclitus populations), and the mean Complex I activity among populations is significantly different at the latitudinal extremes of the range. Importantly, Complex I activity is more similar between F. heteroclitus from the southernmost population and its sister species F. grandis than to the northern populations of F. heteroclitus, suggesting important evolutionary differences. Unexpectedly, the activity is nearly fourfold higher in southern populations than northern populations. Mitochondrial density appears to compensate partially for decreased activity in northern individuals; activity per wet weight is only twofold higher in southern populations. We suggest that some of the variation in Complex I activity is genetically based and thus is being influenced by directional selection. However, this conclusion presents a conundrum: there should not be so much variation in Complex I activity within a population if this variation is biologically important.

Evolutionary factors affecting Lactate dehydrogenase A and B variation in the Daphnia pulex species complex

BACKGROUND

Evidence for historical, demographic and selective factors affecting enzyme evolution can be obtained by examining nucleotide sequence variation in candidate genes such as Lactate dehydrogenase (Ldh). Two closely related Daphnia species can be distinguished by their electrophoretic Ldh genotype and habitat. Daphnia pulex populations are fixed for the S allele and inhabit temporary ponds, while D. pulicaria populations are fixed for the F allele and inhabit large stratified lakes. One locus is detected in most allozyme surveys, but genome sequencing has revealed two genes, LdhA and LdhB.

RESULTS

We sequenced both Ldh genes from 70 isolates of these two species from North America to determine if the association between Ldh genotype and habitat shows evidence for selection, and to elucidate the evolutionary history of the two genes. We found that alleles in the pond-dwelling D. pulex and in the lake-dwelling D. pulicaria form distinct groups at both loci, and the substitution of Glutamine (S) for Glutamic acid (F) at amino acid 229 likely causes the electrophoretic mobility shift in the LDHA protein. Nucleotide diversity in both Ldh genes is much lower in D. pulicaria than in D. pulex. Moreover, the lack of spatial structuring of the variation in both genes over a wide geographic area is consistent with a recent demographic expansion of lake populations. Neutrality tests indicate that both genes are under purifying selection, but the intensity is much stronger on LdhA.

CONCLUSIONS

Although lake-dwelling D. pulicaria hybridizes with the other lineages in the pulex species complex, it remains distinct ecologically and genetically. This ecological divergence, coupled with the intensity of purifying selection on LdhA and the strong association between its genotype and habitat, suggests that experimental studies would be useful to determine if variation in molecular function provides evidence that LDHA variants are adaptive.

Ecologically and evolutionarily important SNPs identified in natural populations

Evolution by natural selection acts on natural populations amidst migration, gene-by-environmental interactions, constraints, and tradeoffs, which affect the rate and frequency of adaptive change. We asked how many and how rapidly loci change in populations subject to severe, recent environmental changes. To address these questions, we used genomic approaches to identify randomly selected single nucleotide polymorphisms (SNPs) with evolutionarily significant patterns in three natural populations of Fundulus heteroclitus that inhabit and have adapted to highly polluted Superfund sites. Three statistical tests identified 1.4-2.5% of SNPs that were significantly different from the neutral model in each polluted population. These nonneutral patterns in populations adapted to highly polluted environments suggest that these loci or closely linked loci are evolving by natural selection. One SNP identified in all polluted populations using all tests is in the gene for the xenobiotic metabolizing enzyme, cytochrome P4501A (CYP1A), which has been identified previously as being refractory to induction in the three highly polluted populations. Extrapolating across the genome, these data suggest that rapid evolutionary change in natural populations can involve hundreds of loci, a few of which will be shared in independent events.

Niche dimensions in fishes: an integrative view

Current shifts in ecosystem composition and function emphasize the need for an understanding of the links between environmental factors and organism fitness and tolerance. The examples discussed here illustrate how recent progress in the field of comparative physiology may provide a better mechanistic understanding of the ecological concepts of the fundamental and realized niches and thus provide insights into the impacts of anthropogenic disturbance. Here we argue that, as a link between physiological and ecological indicators of organismal performance, the mechanisms shaping aerobic scope and passive tolerance set the dimensions of an animal's niche, here defined as its capacity to survive, grow, behave, and interact with other species. We demonstrate how comparative studies of cod or killifish populations in a latitudinal cline have unraveled mitochondrial mechanisms involved in establishing a species' niche, performance, and energy budget. Riverine fish exemplify how the performance windows of various developmental stages follow the dynamic regimes of both seasonal temperatures and river hydrodynamics, as synergistic challenges. Finally, studies of species in extreme environments, such as the tilapia of Lake Magadi, illustrate how on evolutionary timescales functional and morphological shifts can occur, associated with new specializations. We conclude that research on the processes and time course of adaptations suitable to overcome current niche limits is urgently needed to assess the resilience of species and ecosystems to human impact, including the challenges of global climate change.

Regional differentiation and post-glacial expansion of the Atlantic silverside, Menidia menidia, an annual fish with high dispersal potential

The coastal marine environment of the Northwest Atlantic contains strong environmental gradients that create distinct marine biogeographic provinces by limiting dispersal, recruitment, and survival. This region has also been subjected to numerous Pleistocene glacial cycles, resulting in repeated extirpations and recolonizations in northern populations of marine organisms. In this study, we examined patterns of genetic structure and historical demography in the Atlantic silverside, Menidia menidia, an annual marine fish with high dispersal potential but with well-documented patterns of clinal phenotypic adaptation along the environmental gradients of the Northwest Atlantic. Contrary to previous studies indicating genetic homogeneity that should preclude regional adaptation, results demonstrate subtle but significant (FST = 0.07; P < 0.0001) genetic structure among three phylogeographic regions that partially correspond with biogeographic provinces, suggesting regional limits to gene flow. Tests for non-equilibrium population dynamics and latitudinal patterns in genetic diversity indicate northward population expansion from a single southern refugium following the last glacial maximum, suggesting that phylogeographic and phenotypic patterns have relatively recent origins. The recovery of phylogeographic structure and the partial correspondence of these regions to recognized biogeographic provinces suggest that the environmental gradients that shape biogeographic patterns in the Northwest Atlantic may also limit gene flow in M. menidia, creating phylogeographic structure and contributing to the creation of latitudinal phenotypic clines in this species.

SNP identification, verification, and utility for population genetics in a non-model genus

Background

By targeting SNPs contained in both coding and non-coding areas of the genome, we are able to identify genetic differences and characterize genome-wide patterns of variation among individuals, populations and species. We investigated the utility of 454 sequencing and MassARRAY genotyping for population genetics in natural populations of the teleost, Fundulus heteroclitus as well as closely related Fundulus species (F. grandis, F. majalis and F. similis).

Results

We used 454 pyrosequencing and MassARRAY genotyping technology to identify and type 458 genome-wide SNPs and determine genetic differentiation within and between populations and species of Fundulus. Specifically, pyrosequencing identified 96 putative SNPs across coding and non-coding regions of the F. heteroclitus genome: 88.8% were verified as true SNPs with MassARRAY. Additionally, putative SNPs identified in F. heteroclitus EST sequences were verified in most (86.5%) F. heteroclitus individuals; fewer were genotyped in F. grandis (74.4%), F. majalis (72.9%), and F. similis (60.7%) individuals. SNPs were polymorphic and showed latitudinal clinal variation separating northern and southern populations and established isolation by distance in F. heteroclitus populations. In F. grandis, SNPs were less polymorphic but still established isolation by distance. Markers differentiated species and populations.

Conclusions

In total, these approaches were used to quickly determine differences within the Fundulus genome and provide markers for population genetic studies.

Genomic approaches with natural fish populations

Natural populations v. inbred stocks provide a much richer resource for identifying the effects of nucleotide substitutions because natural populations have greater polymorphism. Additionally, natural populations offer an advantage over most common research organisms because they are subject to natural selection, and analyses of these adaptations can be used to identify biologically important changes. Among fishes, these analyses are enhanced by having a wide diversity of species (>28 000 species, more than any other group of vertebrates) living in a huge range of environments (from below freezing to > 46 degrees C, in fresh water to salinities >40 ppt.). Moreover, fishes exhibit many different life-history and reproductive strategies and have many different phenotypes and social structures. Although fishes provide numerous advantages over other vertebrate models, there is still a dearth of available genomic tools for fishes. Fishes make up approximately half of all known vertebrate species, yet <0.2% of fish species have significant genomic resources. Nonetheless, genomic approaches with fishes have provided some of the first measures of individual variation in gene expression and insights into environmental and ecological adaptations. Thus, genomic approaches with natural fish populations have the potential to revolutionize fundamental studies of diverse fish species that offer myriad ecological and evolutionary questions.

Sharp phylogeographic breaks and patterns of genealogical concordance in the brine shrimp Artemia franciscana

Genealogical concordance is a critical overlay of all phylogenetic analyses, irrespective of taxonomic level. To assess such patterns of congruence we have compiled and derived sequence data for two mitochondrial (16S rRNA, COI) and two nuclear (ITS1, p26) markers in 14 American populations of the hypersaline branchiopod Artemia franciscana. Cladistic analysis revealed three reciprocally monophyletic mitochondrial clades. For nuclear DNA, incomplete lineage sorting was evident presumably as a result of slower coalescence or male-mediated dispersal. Our findings capture the genealogical interval between gene splitting and population divergence. In this sense, strong indications are provided in favour of a superspecies status and ongoing speciation in A. franciscana.

Classic approach revitalizes genomics: Complete characterization of a candidate gene for thermal adaptation in two coral reef fishes

Lactate dehydrogenase-B (ldh-b) encodes a metabolic enzyme (LDH-B) which plays an important role in maintaining aerobic performance and in thermal acclimation and/or adaptation of fish. As the first step in understanding the effect this enzyme has on the ability of tropical coral reef fishes to cope with thermal stress, we characterized both coding and non-coding regions of ldh-b in two congeneric perciformes, Plectropomus leopardus and Plectropomus laevis. Ldh-b was 4666 and 4539bp in length in P. leopardus and P. laevis, respectively, with coding regions comprising 1005bp in both species. We report a high level of sequence homology between the coding regions of ldh-b in these two species, with 98.1% identity of nucleotides corresponding to 100% amino acid identity between the deduced protein sequences. Comparison between non-coding (intron) regions of both species revealed the presence of several indels, despite the high level of homology observed (95.9% identity of intron nucleotides). Potential regulatory motifs and elements, including twenty-six simple sequence repeat motifs (mono-, di-, tri- and tetranucleotide) and twenty-three putative microRNA elements are identified within the introns of both species, further supporting recent demonstrations that such short motifs and elements exhibit widespread positioning throughout non-coding regions of the genome. This novel characterization of ldh-b in these two coral reef fishes allows for a wide range of future studies (e.g. analytical comparisons of ldh-b and LDH-B among different fish genera from different thermal environments and habitats).

Comparative characterization of a temperature responsive gene (lactate dehydrogenase-B, ldh-b) in two congeneric tropical fish, Lates calcarifer and Lates niloticus

The characterization of candidate loci is a critical step in obtaining insight into adaptation and acclimation of organisms. In this study of two non-model tropical (to sub-tropical) congeneric perciformes (Lates calcarifer and Lates niloticus) we characterized both coding and non-coding regions of lactate dehydrogenase-B (ldh-b), a locus which exhibits temperature-adaptive differences among temperate and sub-tropical populations of the North American killifish Fundulus heteroclitus. Ldh-b was 5,004 and 3,527 bp in length in L. calcarifer and L. niloticus, respectively, with coding regions comprising 1,005 bp in both species. A high level of sequence homology existed between species for both coding and non-coding regions of ldh-b (> 97% homology), corresponding to a 98.5% amino acid sequence homology. All six known functional sites within the encoded protein sequence (LDH-B) were conserved between the two Lates species. Ten simple sequence repeat (SSR) motifs (mono-, di-, tri- and tetranucleotide) and thirty putative microRNA elements (miRNAs) were identified within introns 1, 2, 5 and 6 of both Lates species. Five single nucleotide polymorphisms (SNPs) were also identified within miRNA containing intron regions. Such SNPs are implicated in several complex human conditions and/or diseases (as demonstrated by extensive genome-wide association studies). This novel characterization serves as a platform to further examine how non-model species may respond to changes in their native temperatures, which are expected to increase by up to 6 degrees C over the next century.

Ecological genetics in the North Atlantic: environmental gradients and adaptation at specific loci

The North Atlantic intertidal community provides a rich set of organismal and environmental material for the study of ecological genetics. Clearly defined environmental gradients exist at multiple spatial scales: there are broad latitudinal trends in temperature, meso-scale changes in salinity along estuaries, and smaller scale gradients in desiccation and temperature spanning the intertidal range. The geology and geography of the American and European coasts provide natural replication of these gradients, allowing for population genetic analyses of parallel adaptation to environmental stress and heterogeneity. Statistical methods have been developed that provide genomic neutrality tests of population differentiation and aid in the process of candidate gene identification. In this paper, we review studies of marine organisms that illustrate associations between an environmental gradient and specific genetic markers. Such highly differentiated markers become candidate genes for adaptation to the environmental factors in question, but the functional significance of genetic variants must be comprehensively evaluated. We present a set of predictions about locus-specific selection across latitudinal, estuarine, and intertidal gradients that are likely to exist in the North Atlantic. We further present new data and analyses that support and contradict these simple selection models. Some taxa show pronounced clinal variation at certain loci against a background of mild clinal variation at many loci. These cases illustrate the procedures necessary for distinguishing selection driven by internal genomic vs. external environmental factors. We suggest that the North Atlantic intertidal community provides a model system for identifying genes that matter in ecology due to the clarity of the environmental stresses and an extensive experimental literature on ecological function. While these organisms are typically poor genetic and genomic models, advances in comparative genomics have provided access to molecular tools that can now be applied to taxa with well-defined ecologies. As many of the organisms we discuss have tight physiological limits driven by climatic factors, this synthesis of molecular population genetics with marine ecology could provide a sensitive means of assessing evolutionary responses to climate change.

Comparative mitochondrial genomics within and among species of killifish

Background

This study was motivated by the observation of unusual mitochondrial haplotype distributions and associated physiological differences between populations of the killifish Fundulus heteroclitus distributed along the Atlantic coast of North America. A distinct "northern" haplotype is fixed in all populations north of New Jersey, and does not appear south of New Jersey except in extreme upper-estuary fresh water habitats, and northern individuals are known to be more tolerant of hyposmotic conditions than southern individuals. Complete mitochondrial genomes were sequenced from individuals from northern coastal, southern coastal, and fresh water populations (and from out-groups). Comparative genomics approaches were used to test multiple evolutionary hypotheses proposed to explain among-population genome variation including directional selection and hybridization.

Results

Structure and organization of the Fundulus mitochondrial genome is typical of animals, yet subtle differences in substitution patterns exist among populations. No signals of directional selection or hybridization were detected. Mitochondrial genes evolve at variable rates, but all genes exhibit very low dN/dS ratios across all lineages, and the southern population harbors more synonymous polymorphism than other populations.

Conclusion

Evolution of mitochondrial genomes within Fundulus is primarily governed by interaction between strong purifying selection and demographic influences, including larger historical population size in the south. Though directional selection and hybridization hypotheses were not supported, adaptive processes may indirectly contribute to partitioning of variation between populations.

Cryptic ecological diversification of a planktonic estuarine copepod, Acartia tonsa

The recent discovery of cryptic species in marine holoplankton, organisms that 'drift' in oceanic currents throughout their life cycle, contrasts with their potential for long-distance passive dispersal and presumably high gene flow. These observations suggest that holoplankton species are adapting to surprisingly small-scale oceanographic features and imply either limited dispersal or strong selection gradients. Acartia tonsa is a widespread and numerically dominant estuarine copepod containing deep mitochondrial lineages within and among populations along the northwestern Atlantic coast. In this study, we intensively investigated A. tonsa populations in Chesapeake Bay with the goals of testing species status for the deep lineages and testing for their association with environmental features over space and time. Phylogenetic analyses of DNA sequences from mitochondrial cytochrome c oxidase I (mtCOI) and the nuclear ribosomal internal transcribed spacer (nITS) resolved two concordant monophyletic clades. Deep divergence between the two clades (13.7% uncorrected sequence divergence for mtCOI and 32.2% for nITS) and genealogical concordance within sympatric populations strongly suggest that the two clades represent reproductively isolated cryptic species. Based on restriction fragment length polymorphisms of mtCOI, representatives from the two clades were found consistently associated with contrasting salinity regimes (oligohaline vs. meso-polyhaline) with an overlap between 2 and 12 PSU in samples from 1995 to 2005. Finding these patterns in one of the best-known estuarine copepods reinforces the conclusion that marine biodiversity is underestimated, not only in terms of species numbers, but also with respect to niche partitioning and the potential importance of ecological divergence in marine holoplankton.

Fundulus as the premier teleost model in environmental biology: opportunities for new insights using genomics

A strong foundation of basic and applied research documents that the estuarine fish Fundulus heteroclitus and related species are unique laboratory and field models for understanding how individuals and populations interact with their environment. In this paper we summarize an extensive body of work examining the adaptive responses of Fundulus species to environmental conditions, and describe how this research has contributed importantly to our understanding of physiology, gene regulation, toxicology, and ecological and evolutionary genetics of teleosts and other vertebrates. These explorations have reached a critical juncture at which advancement is hindered by the lack of genomic resources for these species. We suggest that a more complete genomics toolbox for F. heteroclitus and related species will permit researchers to exploit the power of this model organism to rapidly advance our understanding of fundamental biological and pathological mechanisms among vertebrates, as well as ecological strategies and evolutionary processes common to all living organisms.

Hippocampus queenslandicus Horne, 2001 - a new seahorse species or yet another synonym?

During the past six years, 15 new seahorse species (Syngnathidae: Hippocampus) have been described on the basis of morphological characters. This approach is known to be problematic, and most species names in Hippocampus are now considered to be synonyms. Genetic methods have great potential to resolve the confused taxonomy of the genus, but none have yet been incorporated into species descriptions. In the present study, mitochondrial control region and cytochrome b DNA sequences, as well as morphological data from the recently described Queensland seahorse, Hippocampus queenslandicus Horne, 2001, were compared with corresponding data from closely related seahorse species to determine whether there is strong support for distinction of this taxon. The haplotypes of H. queenslandicus were nested among haplotypes belonging to two of the three major Southeast Asian lineages of H. spinosissimus Weber, 1913. Although incomplete lineage sorting characteristic of very recently diverged species cannot be ruled out, the genetic results suggest that H. queenslandicus is paraphyletic. Morphometric analysis further fails to provide strong support for the species status of H. queenslandicus. We conclude that support for the distinctness of H. queenslandicus is weak, and indicate that it is a synonym of H. spinosissimus. The taxonomic validity of other recently described seahorse species should be similarly scrutinised using combined genetic and detailed morphological methods.

The Pleistocene history of the sheepshead minnow (Cyprinodon variegatus): Non-equilibrium evolutionary dynamics within a diversifying species complex

The sheepshead minnow, Cyprinodon variegatus, is a widespread fish species that typically inhabits coastal tidal marsh and mangrove swamp environments, ranging from Cape Cod, Massaschusetts to northern Mexico and into the Caribbean. This wide range crosses several biogeographic boundaries which are coincident with genetic structuring within numerous species originating in the Pleistocene. In addition, the more northerly reaches of this species range have been further subject to the evolutionary consequences of Pleistocene glaciation due to local extinction and recolonization of formerly glaciated sites. C. variegatus thus provides an excellent vertebrate model system within which to test the extent of genetic differentiation among populations in a dominant coastal ecosystem and examine patterns of historical demography in populations distributed along a latitudinal gradient. Using mitochondrial control region and ND2 sequence data, we discovered monophyletic clades within C. variegatus with divergence times within the Pleistocene, and very low gene flow between most sites. Intraspecific genetic breaks appear to correspond broadly to biogeographic or oceanic boundaries. Pleistocene climate change appears to have had dramatic impacts on the size and distribution of populations within and near the glacial margins, but has also affected populations far from formerly glaciated regions.

Intraspecific variation in thermal tolerance and heat shock protein gene expression in common killifish,Fundulus heteroclitus

Populations of common killifish, Fundulus heteroclitus, are distributed along the Atlantic coast of North America through a steep latitudinal thermal gradient. We examined intraspecific variation in whole-animal thermal tolerance and its relationship to the heat shock response in killifish from the northern and southern extremes of the species range. Critical thermal maxima were significantly higher in southern than in northern fish by ∼1.5°C at a wide range of acclimation temperatures (from 2-34°C), and critical thermal minima differed by ∼1.5°C at acclimation temperatures above 22°C, converging on the freezing point of brackish water at lower acclimation temperatures. To determine whether these differences in whole-organism thermal tolerance were reflected in differences in either the sequence or regulation of the heat shock protein genes(hsps) we obtained complete cDNA sequences for hsc70, hsp70-1 and hsp70-2, and partial sequences of hsp90α and hsp90β. There were no fixed differences in amino acid sequence between populations in either hsp70-1 or hsp70-2, and only a single conservative substitution between populations in hsc70. By contrast, there were significant differences between populations in the expression of many, but not all, of these genes. Both northern and southern killifish significantly increased hsp70-2 levels above control values(Ton) at a heat shock temperature of 33°C, but the magnitude of this induction was greater in northern fish, suggesting that northern fish may be more susceptible to thermal damage than are southern fish. In contrast, hsp70-1 mRNA levels increased gradually and to the same extent in response to heat shock in both populations. Hsc70 mRNA levels were significantly elevated by heat shock in southern fish, but not in northern fish. Similarly, the more thermotolerant southern killifish had a Ton for hsp90α of 30°C, 2°C lower than that of northern fish. This observation combined with the ability of southern killifish to upregulate hsc70 in response to heat shock suggests a possible role for these hsps in whole-organism differences in thermal tolerance. These data highlight the importance of considering the complexity of the heat shock response across multiple isoforms when attempting to make linkages to whole-organism traits such as thermal tolerance.

Microsatellite analysis of the phylogeography, Pleistocene history and secondary contact hypotheses for the killifish, Fundulus heteroclitus

The mummichog, Fundulus heteroclitus, exhibits extensive latitudinal clinal variation in a number of physiological and biochemical traits, coupled with phylogeographical patterns at mitochondrial and nuclear DNA loci that suggest a complicated history of spatially variable selection and secondary intergradation. This species continues to serve as a model for understanding local and regional adaptation to variable environments. Resolving the influences of historical processes on the distribution of genetic variation within and among extant populations of F. heteroclitus is crucial to a better understanding of how populations evolve in the context of contemporary environments. In this study, we analysed geographical patterns of genetic variation at eight microsatellite loci among 15 populations of F. heteroclitus distributed throughout the North American range of the species from Nova Scotia to Georgia. Genetic variation in Northern populations was lower than in Southern populations and was strongly correlated with latitude throughout the species range. The most common Northern alleles at all eight loci exhibited concordant latitudinal clinal patterns, and the existence of an abrupt transition zone in allele frequencies between Northern and Southern populations was similar to that observed for mitochondrial DNA and allozyme loci. A significant pattern of isolation by distance was observed both within and between northern and southern regions. This pattern was unexpected, particularly for northern populations, given the recent colonization history of post-Pleistocene habitats, and was inconsistent with either a recent northward population expansion or a geographically restricted northern Pleistocene refugium. The data provided no evidence for recent population bottlenecks, and estimates of historical effective population sizes suggest that post-Pleistocene populations have been large throughout the species distribution. These results suggest that F. heteroclitus was broadly distributed throughout most of its current range during the last glacial event and that the abrupt transition in allele frequencies that separate Northern and Southern populations may reflect regional disequilibrium conditions associated with the post-Pleistocene colonization history of habitats in that region.

Neutral and adaptive variation in gene expression

Variation among populations in gene expression should be related to the accumulation of random-neutral changes and evolution by natural selection. The following evolutionary analysis has general applicability to biological and medical science because it accounts for genetic relatedness and identifies patterns of expression variation that are affected by natural selection. To identify genes evolving by natural selection, we allocate the maximum among-population variation to genetic distance and then examine the remaining variation relative to a hypothesized important ecological parameter (temperature). These analyses measure the expression of metabolic genes in common-gardened populations of the fish Fundulus heteroclitus whose habitat is distributed along a steep thermal gradient. Although much of the variation in gene expression fits a null model of neutral drift, the variation in expression for 22% of the genes that regress with habitat temperature was far greater than could be accounted for by genetic distance alone. The most parsimonious explanation for among-population variation for these genes is evolution by natural selection. In addition, many metabolic genes have patterns of variation incongruent with neutral evolution: They have too much or too little variation. These patterns of biological variation in expression may reflect important physiological or ecological functions.

Molecular systematics of Hispaniolan pupfishes (Cyprinodontidae: Cyprinodon): implications for the biogeography of insular Caribbean fishes

We used sequence variation in the mtDNA control-region and ND2 and cyt b genes to assess the systematics and biogeography of the five species of pupfish (Cyprinodon) on Hispaniola. These include four endemics, the relatively large-bodied Cyprinodon bondi, Cyprinodon nichollsi, and Cyprinodon sp., each from a separate lake in southwestern Hispaniola, and Cyprinodon higuey from a coastal lake in eastern Hispaniola. The fifth species consists of coastal populations referable to Cyprinodon variegatus riverendi. The results indicate that Hispaniola has been invaded by at least two forms, first by a late Pliocene progenitor of Cyprinodon variegatus ovinus and the large-bodied Hispaniolan species, and, more recently, by one or more ancestral forms allied with Cyprinodon variegatus variegatus and C. v. riverendi. Levels of divergence indicate that large expanses of open sea have not acted as long-term barriers to inter-island dispersal of cyprinodontiform fishes. This study, together with the molecular systematics of other insular Caribbean fishes, indicates that most insular groups originated from late Neogene dispersal from the mainland. The patterns of mtDNA variation in Cyprinodon showed little congruence with the species/subspecies taxonomy.

Stress and interpopulation variation in glycolytic enzyme activity and expression in a teleost fish Fundulus heteroclitus

Northern populations of Fundulus heteroclitus have twofold greater activity of lactate dehydrogenase-B (LDH-B) than southern populations, but exposure to stress increases LDH-B in southern populations, abolishing this difference. To test whether differences in the activity of other hepatic glycolytic enzymes between populations are sensitive to stress, we injected fish with a pharmacological dose of cortisol in coconut oil (400 microg g(-1)) or exposed them to handling stress and measured the activities of all the glycolytic enzymes. At rest, liver phosphofructokinase (PFK) and aldolase (ALD) activities were greater in southern fish, whereas LDH-B activity was greater in northern fish. No other glycolytic enzymes differed in activity between populations in control fish. Cortisol injection and handling stress decreased PFK and ALD and increased LDH activities in the southern but not the northern population, such that the populations no longer differed in the activity of any enzyme following treatment. Unlike Ldh-B mRNA, Pfk and Ald mRNA levels did not parallel enzyme activity, suggesting complex kinetics or regulation at multiple levels. Plasma cortisol did not differ between populations at rest but was significantly different between populations in treated fish. These data suggest that differences in liver enzyme activity may be related to differences in stress hormone physiology between populations.

Differential gene expression in anthracene-exposed mummichogs (Fundulus heteroclitus)

The polycyclic aromatic hydrocarbon (PAH) anthracene is present in many estuarine systems at concentrations believed to cause sublethal adverse effects, although its exact mode of toxicity remains unclear. Knowledge of the induction or suppression of specific genes as a result of exposure may be useful in explaining these effects. We have generated a fingerprint of anthracene exposure using the mummichog (Fundulus heteroclitus), a non-migratory estuarine fish species. The fish were exposed in 7-day static renewal tests to environmentally relevant concentrations of 0, 27, 50, and 80 microg/l of anthracene. Total RNA was extracted from the livers and differential display reverse transcription polymerase chain reaction (DD RT-PCR) was used to recover 26 differentially expressed cDNA fragments. These cDNAs were isolated, sequenced, and compared to sequences of known genes in order to identify possible physiological consequences of exposure to anthracene. We then constructed macroarrays using these fragments and probed them with RNA from both anthracene-exposed fish and fish from a known PAH-impacted site. Three genes appear to be good indicators of exposure to anthracene in the range of concentrations tested, which included CYP2N2 and two expressed sequence tags (ESTs) termed 15C1 and 18C2. The expression of nine genes was altered in fish collected from a site with multiple PAHs. Band 15C1 and CYP2N2 again showed statistically significant upregulation in the field-caught fish, while a trypsin precursor and fatty acid-binding protein (FABP) all showed similar trends in induction as the laboratory-exposed fish. Further insight into the mechanism of toxicity of contaminants will be gained by the ability to identify and use differentially expressed genes as markers of exposure and effects.

Genetic variation in a cline in a living intertidal snail arose in the Neogastropoda over 100 million years ago

To understand how species adapt and evolve it is necessary to appreciate the relationship between genetic variation and the environment. Here, the fossil record and molecular data from different lineages of the marine Gastropoda are used to understand the evolution of genetic variations found in the nuclear gene for mitochondrial malate dehydrogenase (mMDH; EC 1.1.1.37) in the living intertidal Muricid snail, Nucella lapillus. Assuming a molecular clock, DNA sequences of mMDH indicate that two variants found in N. lapillus, mMDH(9) and mMDH(10), may have arisen as long as 144 MY (million years) ago and at least prior to the evolution of the Muricidae approximately 112-90 MY ago. The Muricidae contain by far the greatest majority of the Neogastropoda specialized for life in the intertidal habitat. In N. lapillus the mMDH(9) and mMDH(10) variants covary with variations in other biochemically defined loci, inherited phenotypic traits (shell shape and physiology) and karyotype frequencies to differentiate two distinct nuclear haplotypes that are associated with different temperature environments. The variations in shell shape that are associated with the haplotypes of N. lapillus represent adaptations to temperature stress and similar variations occur in other related intertidal molluscs whose lineages are much older than Nucella, which arose around 25 million years ago. It is suggested that the divergence of the mMDH variants found in N. lapillus may reflect an ancient genetic event, such as a chromosomal mutation, perhaps involving variation in other linked traits that together became important in the subsequent evolution of the marine Gastropoda.

The study of structured populations--new hope for a difficult and divided science

Natural populations, including those of humans, have complex geographies and histories. Studying how they evolve is difficult, but it is possible with population-based DNA sequence data. However, the study of structured populations is divided by two distinct schools of thought and analysis. The phylogeographic approach is fundamentally graphical and begins with a gene-tree estimate. By contrast, the more traditional approach of using summary statistics is fundamentally mathematical. Both approaches have limitations, but there is promise in newer probabilistic methods that offer the flexibility and data exploitation of the phylogeographic approach in an explicitly model-based mathematical framework.

Phylogeographic congruence between mtDNA and rDNA ITS markers in brown trout

Variation in the internal transcribed spacer (ITS) of rDNA was examined throughout the range of the brown trout (Salmo trutta) to analyze the usefulness of this molecular marker for phylogeographic analysis. The results were compared with those previously obtained with mtDNA, a region exhaustively analyzed along the brown trout distribution. ITS2 was essentially conserved at all populations sampled, no informative characters being detected across the main lineages described in this species. Conversely, ITS1 showed a greater homogenization than other genetic markers at a microgeographic scale, with variation partitioning into several major phylogenetic groups. Phylogeographic patterns were partially congruent between both ITS1 and mtDNA. The main discrepancies were the detection of intra-individual variation and putative recombinant ITS1 sequences in hybridization areas between genetically different, yet historically overlapping, assemblages. Also, the existence of an ancient ITS1 sequence in the Mediterranean-southeastern area (rMEDA), not revealed by mtDNA analysis, was evidenced after rDNA ITS1 analysis.

Strong positive selection and habitat-specific amino acid substitution patterns in MHC from an estuarine fish under intense pollution stress

Population-level studies using the major histocompatibility complex (Mhc) have linked specific alleles with specific diseases, but data requirements are high and the power to detect disease association is low. A novel use of Mhc population surveys involves mapping allelic substitutions onto the inferred structural molecular model to show functional differentiation related to local selective pressures. In the estuarine fish Fundulus heteroclitus, populations experiencing strong differences in antigenic challenges show significant differences in amino acid substitution patterns that are reflected as variation in the structural location of changes between populations. Fish from a population genetically adapted to severe chemical pollution also show novel patterns of DNA substitution at a highly variable Mhc class II B locus including strong signals of positive selection at inferred antigen-binding sites and population-specific signatures of amino acid substitution. Heavily parasitized fish from an extreme PCB-contaminated (U.S. Environmental Protection Agency Superfund) site show enhanced population-specific substitutions in the a-helix portion of the inferred antigen-binding region. In contrast, fish from an unpolluted site show a significantly different pattern focused on the first strand of the B-pleated sheet. Whether Mhc population profile differences represent the direct effects of chemical toxicants or indirect parasite-mediated selection, the result is a composite habitat-specific signature of strong selection and evolution affecting the genetic repertoire of the major histocompatibility complex.

Inferring the history of speciation from multilocus DNA sequence data: the case of Drosophila pseudoobscura and close relatives

The divergence of Drosophila pseudoobscura from its close relatives, D. persimilis and D. pseudoobscura bogotana, was examined using the pattern of DNA sequence variation in a common set of 50 inbred lines at 11 loci from diverse locations in the genome. Drosophila pseudoobscura and D. persimilis show a marked excess of low-frequency variation across loci, consistent with a model of recent population expansion in both species. The different loci vary considerably, both in polymorphism levels and in the levels of polymorphisms that are shared by different species pairs. A major question we address is whether these patterns of shared variation are best explained by gene flow or by persistence since common ancestry. A new test of gene flow, based on patterns of linkage disequilibrium, is developed. The results from these, and other tests, support a model in which D. pseudoobscura and D. persimilis have exchanged genes at some loci. However, the pattern of variation suggests that most gene flow, although occurring after speciation began, was not recent. There is less evidence of gene flow between D. pseudoobscura and D. p. bogotana. The results are compared with recent work on the genomic locations of genes that contribute to reproductive isolation between D. pseudoobscura and D. persimilis. We show that there is a good correspondence between the genomic regions associated with reproductive isolation and the regions that show little or no evidence of gene flow.

Structure and sequence conservation of a putative hypoxia response element in the lactate dehydrogenase-B gene of Fundulus

Many aquatic habitats are characterized by periodic or sustained episodes of low oxygen concentration, or hypoxia, and organisms that survive in these habitats do so by utilizing a suite of behavioral, physiological and biochemical adjustments to low oxygen (1-3). In the killifish Fundulus heteroclitus, one response to prolonged exposure to hypoxia is an increase in the activity of lactate dehydrogenase-B (LDH-B), the terminal enzyme of anaerobic glycolysis, in liver tissue (4). An increase in glycolytic enzyme activity also occurs in mammalian cells during hypoxia, a process due, in part, to increased rates of gene transcription mediated by the hypoxia-inducible transcription factor, HIF-1 (5). Given that a homolog of HIF-1 has been identified in fish (6), we hypothesized that HIF might be involved in the observed up-regulation of LDH-B in F. heteroclitus. Herein, we describe the presence of DNA elements in intron 2 of the Ldh-B gene from F. heteroclitus that resemble hypoxia response elements (HRE) describedfor mammalian genes (7-10). Specifically, over a region of approximately 50 base pairs we identified two consensus HIF-1 binding sites, as well as DNA elements that may bind other transcription factors (e.g., cyclic AMP response elements; CRE). We found that these sites were perfectly conserved among geographically diverse populations of F. heteroclitus, as well as being highly conserved among multiple species in the genus Fundulus. The spacing, orientation, and sequence conservation of these putative regulatory elements suggest that they may be functionally involved in the hypoxic regulation of Ldh-B in these fish.

Environmental adaptations as windows on molecular evolution

Changes in gene regulation may play an important role in adaptive evolution, particularly during adaptation to a changing environment. However, little is known about the molecular mechanisms underlying adaptively significant variation in gene regulation. To address this question, we are using environmental adaptations in populations of a fish, Fundulus heteroclitus as a window into the molecular evolution of gene regulation. F. heteroclitus are found along the East Coast of North America, with populations distributed along a steep thermal gradient. At the extremes of the species range, populations have undergone local adaptation to their habitat temperatures. A variety of genes differ in their regulation between these populations. We have determined the mechanism responsible for changes in lactate dehydrogenase-B (Ldh-B) gene regulation. A limited number of mutations in the regulatory sequence of this gene result in changes in its expression. Both the phenotypic (increased LDH activity) and genotypic (changes in Ldh-B regulatory sequences) differences between populations have been shown to be affected by natural selection, rather than genetic drift. Therefore, even a small number of mutations within important regulatory sequences can provide evolutionarily significant variation and have an impact on environmental adaptation.

Genetic evidence for limited dispersal in the coastal California killifish, Fundulus parvipinnis

The California killifish, Fundulus parvipinnis, is a marine species that lives in salt marshes, estuaries and wetlands along the California and Baja California coasts. In order to estimate levels of dispersal between different coastal habitats over its range, we have studied six populations using morphological and genetic markers. Lateral line scale and vertebrae counts showed significant differences between individuals collected north of Punta Eugenia and south of Punta Eugenia. Morphological differences across Punta Eugenia were accompanied by large genetic differences at the mitochondrial control region (5.8%). Gene flow was in general very reduced over the range of the species (pairwise average F(st)=0.70, Nm=0.30), with a strong break at Punta Eugenia (F(st)=0.95, Nm=0.03). Such limited interchanges between coastal habitats have important theoretical and conservation implications.

An ancient transpecific polymorphism shows extreme divergence in a multitrait cline in an intertidal snail (Nucella lapillus (L.))

Clines in intraspecific genetic variation are frequently associated with an environmental transition. Here, divergence among nucleotide sequences of two nuclear loci, cytosolic and mitochondrial malate dehydrogenase (cMDH and mMDH, respectively), is described, in a multitrait cline over a distance of ca. 3 km where shell phenotype, allozyme, mitochondrial DNA haplotype, and centric fusion (Robertsonian translocations) frequencies covary with temperature and humidity and change abruptly in a continuous population of the dog-whelk (Nucella lapillus), a common intertidal snail of the north temperate Atlantic. Protein electrophoresis has already shown two alleles of mMDH varying from fixation of one allele to near fixation of the other, whereas cMDH appears to be monomorphic. The results of this study show a striking disparity in nucleotide sequence divergence among alleles at the two loci, with extreme molecular differentiation in one of them. Four alleles of cMDH were found to have nucleotide and amino acid sequence divergences of 0.4% and 0.3%, respectively. In contrast, the two mMDH cDNA alleles differed by 23% and 20% at the nucleotide and amino acid levels, respectively. Analysis of a 91-bp partial nucleotide sequence of mMDH from Nucella freycineti, the closest relative of N. lapillus, revealed two similar alleles and indicated that the divergence in mMDH in N. lapillus represents an ancient transpecific polymorphism in these Nucella. Together with earlier studies on variation in N. lapillus, it is argued that the polymorphism in mMDH and the clines in N. lapillus represent the presence of two persistent coadapted gene complexes, multitrait coevolving genetic solutions to environmental variation, which may presently enable this snail to exploit a diverse environment successfully.

Cloning and primary structure of putative cytosolic and mitochondrial malate dehydrogenase from the mollusc Nucella lapillus (L.)

The evolutionary history of the malate dehydrogenase (MDH) gene family [NAD-dependent MDH; EC 1.1.1.37 and NAD(P)-dependent MDH; EC 1.1.1.82] has received much attention. MDHs have also featured extensively as electrophoretic markers in population genetics and evolutionary ecology, and in many cases, intraspecific variation in MDH has been correlated with environmental variables. However, while the amino acid residues essential for MDH function are known, no studies have examined intraspecific nucleotide variation despite evidence indicating that natural selection may be operating on this locus. This study presents two sets of degenerate oligonucleotide PCR primers to facilitate the cloning of cytosolic MDH (cMDH) and mitochondrial MDH (mMDH) from a broad range of animals (cMDH) and animals and plants (mMDH). These primers were used to obtain putative cMDH and mMDH cDNAs from the mollusc Nucella lapillus. The N. lapillus cMDH cDNA was found to encode a putative cMDH protein of 334aa and 36kDa, while the mMDH cDNA encoded a putative mature mMDH protein of 315aa and 33kDa. The putative amino acid sequences of the two compartmentalised N. lapillus MDHs are presented and compared to other known MDH sequences.

Barriers to gene flow in Embiotoca jacksoni, a marine fish lacking a pelagic larval stage

Marine species generally show high dispersal capabilities, which should be accompanied by high levels of gene flow and low speciation rates. However, studies that focused on the relationship between dispersal and gene flow in marine fishes have been inconclusive. This study focuses on the black surfperch, Embiotoca jacksoni, a temperate reef fish that lacks a pelagic larval stage and lives on almost continuous reefs along the California and Baja California coasts. Mitochondrial control-region sequences from 240 individuals were obtained, and phylogeographic patterns were analyzed. A major phylogeographic break was found at Santa Monica Bay, a sandy expanse that prevents adult dispersal. Deep water separating the southern California Channel Islands was also found to be a major barrier to gene flow. Minor phylogeographic breaks were also detected in the Big Sur/Morro Bay and in the Punta Eugenia/Guerrero Negro regions, but none in the Point Conception region. Gene flow levels in E. jacksoni were found to be almost identical to those of another species with limited dispersal, Acanthochromis polyacanthus, thus indicating that the lack of a pelagic larval stage combined with barriers to adult dispersal may have had similar effects on these two species.

Evolution of a functionally related lactate dehydrogenase and pyruvate decarboxylase pseudogene complex in maize

A large proportion of the maize genome is repetitive DNA (60-80%) with retrotransposons contributing significantly to the repetitive DNA component. The majority of retrotransposon DNA is located in intergenic regions and is organized in a nested fashion. Analysis of an 8.2-kb segment of maize genomic DNA demonstrated the presence of three retrotransposons of different reiteration classes in addition to lactate dehydrogenase and pyruvate decarboxylase pseudogenes. Both of the pseudogenes were located within a defective retrotransposon element (LP-like element) which possessed identical long terminal repeats (LTRs) with inverted repeats at each end, a primer binding site, a polypurine tract, and generated a 5-bp target site duplication. A model describing the events leading to the formation of the LP-like element is proposed.Key words: lactate dehydrogenase, LP-like element, pseudogene, pyruvate decarboxylase, retrotransposon.