Effects of contaminants on genetic patterns in aquatic organisms: a review.

There is increasing awareness of the need to evaluate the effects of contaminants at the population level. Genetic techniques offer a powerful approach to assess contaminant-induced changes in populations. Yet studies to date are relatively few and not always carefully designed to maximize the utility inherent in this approach. We present a summary of contemporary genetic assessment methods and a review of published studies of genetic effects in field-exposed aquatic organisms. We discuss evaluations of genetic patterns that use genetic adaptation, allozyme variation, and molecular genetic (DNA) variation. Direct tests of genetic adaptation are very effective in establishing a concrete, and potentially deleterious population-level effect of contaminant exposure, but they are difficult to accomplish with most field-exposed organisms. Allozyme surveys are relatively simple and common, and may provide data that are suggestive of contaminant effects. However, these are rarely conclusive, primarily because few allozyme loci are variable and these few loci represent extremely small portions of the genome. Molecular genetic techniques have the potential to be very effective. But, there is a tendency to emphasize the power of the techniques, rather than the underlying causes of the molecular genetic patterns observed. The strength of the conclusions of each study varies widely, partially derived from variation in the strength of the techniques. We caution that all these approaches are greatly improved by careful experimental design that includes adequate numbers of reference and contaminated sites and sample size. In addition, careful exposure assessment is required, including site and tissue chemistry, biomarker responses, and measures of potentially deleterious effects, such as DNA damage, or reduced reproductive output or survival.     

Evidence for shared ancestral polymorphism rather than recurrent gene flow at microsatellite loci differentiating two hybridizing oaks (Quercus spp.)

Quercus petraea and Quercus robur are two closely related oak species, considered to hybridize. Genetic markers, however, indicate that despite sharing most alleles, the two species remain separate genetic units. Analysis of 20 microsatellite loci in multiple populations from both species suggested a genome-wide differentiation. Thus, the allele sharing between both species could be explained either by low rates of gene flow or shared ancestral variation. We performed further analyses of population differentiation in a biogeographical setting and an admixture analysis in mixed oak stands to distinguish between both hypotheses. Based on our results we propose that the low genetic differentiation among these species results from shared ancestry rather than high rates of gene flow.     

Variation in mitochondrial DNA reveals multiple distant glacial refugia in black spruce (Picea mariana), a transcontinental North American conifer

Range-wide genetic variation of black spruce (Picea mariana) was studied using polymerase chain reaction-random fragment length polymorphism markers of the mitochondrial genome. Four polymorphic mitochondrial DNA (mtDNA) loci were surveyed and two or three alleles were detected at each locus, resulting in 10 multilocus mtDNA types or mitotypes. A significant subdivision of population genetic diversity was detected (GST = 0.671; NST = 0.726), suggesting low levels of gene flow among populations. The distribution of mitotypes was not random (NST > GST; P < 0.05) and revealed four partially overlapping zones, presumably representative of different glacial populations. Comparison of the genetic structure derived from mtDNA markers and the colonization paths previously deduced from the fossil and pollen records allow us to infer at least three southern and one northeastern glacial populations for black spruce. The patterns revealed in this study suggest that black spruce shares its biogeographical history with other forest-associated North American species.     

Insights into the biogeographical history of the Lower Guinea Forest Domain: evidence for the role of refugia in the intraspecific differentiation of Aucoumea klaineana

Determining the biogeographical histories of rainforests is central to our understanding of the present distribution of tropical biodiversity. Ice age fragmentation of central African rainforests strongly influenced species distributions. Elevated areas characterized by higher species richness and endemism have been postulated to be Pleistocene forest refugia. However, it is often difficult to separate the effects of history and of present-day ecological conditions on diversity patterns at the interspecific level. Intraspecific genetic variation could yield new insights into history, because refugia hypotheses predict patterns not expected on the basis of contemporary environmental dynamics. Here, we test geographically explicit hypotheses of vicariance associated with the presence of putative refugia and provide clues about their location. We intensively sampled populations of Aucoumea klaineana, a forest tree sensitive to forest fragmentation, throughout its geographical range. Characterizing variation at 10 nuclear microsatellite loci, we were able to obtain phylogeographic data of unprecedented detail for this region. Using Bayesian clustering approaches, we demonstrated the presence of four differentiated genetic units. Their distribution matched that of forest refugia postulated from patterns of species richness and endemism. Our data also show differences in diversity dynamics at leading and trailing edges of the species' shifting distribution. Our results confirm predictions based on refugia hypotheses and cannot be explained on the basis of present-day ecological conditions.     

Apportionment of global human genetic diversity based on craniometrics and skin color

A number of analyses of classical genetic markers and DNA polymorphisms have shown that the majority of human genetic diversity exists within local populations (approximately 85%), with much less among local populations (approximately 5%) or between major geographic regions or "races" (approximately 10%). Previous analysis of craniometric variation (Relethford [1994] Am J Phys Anthropol 95:53-62) found that between 11-14% of global diversity exists among geographic regions, with the remaining diversity existing within regions. The methods used in this earlier paper are extended to a hierarchical partitioning of genetic diversity in quantitative traits, allowing for assessment of diversity among regions, among local populations within regions, and within local populations. These methods are applied to global data on craniometric variation (57 traits) and skin color. Multivariate analysis of craniometric variation shows results similar to those obtained from genetic markers and DNA polymorphisms: roughly 13% of the total diversity is among regions, 6% among local populations within regions, and 81% within local populations. This distribution is concordant with neutral genetic markers. Skin color shows the opposite pattern, with 88% of total variation among regions, 3% among local populations within regions, and 9% within local populations, a pattern shaped by natural selection. The apportionment of genetic diversity in skin color is atypical, and cannot be used for purposes of classification. If racial groups are based on skin color, it appears unlikely that other genetic and quantitative traits will show the same patterns of variation.     

The burgeoning field of statistical phylogeography

In the newly emerging field of statistical phylogeography, consideration of the stochastic nature of genetic processes and explicit reference to theoretical expectations under various models has dramatically transformed how historical processes are studied. Rather than being restricted to ad hoc explanations for observed patterns of genetic variation, assessments about the underlying evolutionary processes are now based on statistical tests of various hypotheses, as well as estimates of the parameters specified by the models. A wide range of demographical and biogeographical processes can be accommodated by these new analytical approaches, providing biologically more realistic models. Because of these advances, statistical phylogeography can provide unprecedented insights about a species' history, including decisive information about the factors that shape patterns of genetic variation, species distributions, and speciation. However, to improve our understanding of such processes, a critical examination and appreciation of the inherent difficulties of historical inference and challenges specific to testing phylogeographical hypotheses are essential. As the field of statistical phylogeography continues to take shape many difficulties have been resolved. Nonetheless, careful attention to the complexities of testing historical hypotheses and further theoretical developments are essential to improving the accuracy of our conclusions about a species' history.     

Identification, distribution and molecular evolution of the pacifastin gene family in Metazoa

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.     

Genetic structure among and within peripheral and central populations of three endangered floodplain violets

Understanding the partitioning of genetic variance in peripheral and central populations may shed more light on the effects of genetic drift and gene flow on population genetic structure and, thereby, improve attempts to conserve genetic diversity. We analysed genetic structure of peripheral and central populations of three insect-pollinated violets (Viola elatior, Viola pumila, Viola stagnina) to evaluate to what extent these patterns can be explained by gene flow and genetic drift. Amplified fragment length polymorphism was used to analyse 930 individuals of 50 populations. Consistent with theoretical predictions, peripheral populations were smaller and more isolated, differentiation was stronger, and genetic diversity and gene flow lower in peripheral populations of V. pumila and V. stagnina. In V. elatior, probably historic fragmentation effects linked to its specific habitat type were superimposed on the plant geographic (peripheral-central) patterns, resulting in lower relative importance of gene flow in central populations. Genetic variation between regions (3-6%), among (30-37%) and within populations (60-64%) was significant. Peripheral populations lacked markers that were rare and localized in central populations. Loss of widespread markers in peripheral V. stagnina populations indicated genetic erosion. Autocorrelation within populations was statistically significant up to a distance of 10-20 m. Higher average genetic similarity in peripheral populations than in central ones indicated higher local gene flow, probably owing to management practices. Peripheral populations contributed significantly to genetic variation and contained unique markers, which made them valuable for the conservation of genetic diversity.     

TESTING MODELS OF MIGRATION AND ISOLATION AMONG POPULATIONS OF CHINOOK SALMON (ONCORHYNCHUS TSCHAWYTSCHA)

The chinook salmon (Oncorhynchus tschawytscha) is a behaviorally, morphologically, and ecologically variable species distributed over a large geographic range. Although previous genetic surveys have revealed considerable genetic differences among populations with different life history types and from different major river drainages, it is not clear to what degree these genetically distinct populations are connected by low levels of gene flow. The work described in this paper addresses this question by surveying DNA restriction site variation at six nuclear genes from nine populations encompassing most of the species's North American range, and then attempting to fit the patterns of variation observed at these genes to five different evolutionary models using computer simulations of the coalescent process. Two commonly used constant population size models, one hypothesizing no migration among populations and one hypothesizing equal rates of migration among populations, were found to be statistically inconsistent with the observed patterns of variation. The other three models, which involved either recent divergence among populations coupled with large changes in populations size, unequal migration rates among populations, or selection, were all found to be consistent with the observed patterns of variation. Assuming selective neutrality, these results suggest that either the populations have all descended from a common ancestral population within the last ~50,000 years and have all suffered large declines in effective population size since that time, or that they have a more ancient divergence time but are connected by low levels of gene flow. These conclusions rest on the assumption of selective neutrality. With the methods employed, it was not possible to simultaneously test hypotheses of both selective neutrality and population structure.     

Relative influences of historical and contemporary forces shaping the distribution of genetic variation in the Atlantic killifish, Fundulus heteroclitus

A major goal of population genetics research is to identify the relative influences of historical and contemporary processes that serve to structure genetic variation. Most population genetic models assume that populations exist in a state of migration-drift equilibrium. However, in the past this assumption has rarely been verified, and is likely rarely achieved in natural populations. We assessed the equilibrium status at both local and regional scales of the Atlantic killifish, Fundulus heteroclitus . This species is a model organism for the study of adaptive clinal variation, but has also experienced a complicated history of range expansion and secondary contact following allopatric divergence, potentially obscuring the influence of contemporary evolutionary processes. Presumptively neutral genetic markers (microsatellites) demonstrated zones of secondary intergradation among coastal populations centred around northern New Jersey and the Chesapeake Bay region. Analysis of genetic variation indicated isolation by distance among some populations and provided supporting evidence that the Delaware Bay, but not the Chesapeake Bay, has acted as a barrier to dispersal among coastal populations. Bayesian estimates indicated large effective population sizes and low migration rates, and were in good agreement with empirically derived estimates of population and neighbourhood size from mark–recapture studies. These data indicate that populations are not in migration-drift equilibrium at a regional scale, and suggest that contributing factors include large population size combined with relatively low migration rates. These conditions should be considered when interpreting the evolutionary significance of the distribution of genetic variation among F. heteroclitus populations.     

Concordance of genetic divergence among sockeye salmon populations at allozyme, nuclear DNA, and mitochondrial DNA markers

Abstract.— We examined genetic variation at 21 polymorphic allozyme loci, 15 nuclear DNA loci, and mitochondrial DNA in four spawning populations of sockeye salmon ( Oncorhynchus nerka ) from Cook Inlet, Alaska, to test for differences in the patterns of divergence among different types of markers. We were specifically interested in testing the suggestion that natural selection at allozyme loci compromises the effectiveness of these markers for describing the amount and patterns of gene flow among populations. We found concordance among markers in the amount of genetic variation within and among populations, with the striking exception of one allozyme locus ( sAH ), which exhibited more than three times the amount of among-population differentiation as other loci. A consideration of reports of discordance between allozymes and other loci indicates that these differences usually result from one or two exceptional loci. We conclude that it is important to examine many loci when estimating genetic differentiation to infer historical amounts of gene flow and patterns of genetic exchange among populations. It is less important whether those loci are allozymes or nuclear DNA markers.     

Historical isolation,range expansion,and secondary contact of two highly divergent mitochondrial lineages in spotted salamanders (Ambystoma maculatum)

Abstract.— The high species diversity of aquatic and terrestrial faunas in eastern North America has been attributed to range reductions and allopatric diversification resulting from historical climate change. The role these processes may have played in speciation is still a matter of considerable debate; however, their impacts on intraspecific genetic structure have been well documented. We use mitochondrial DNA sequences to reconstruct an intraspecific phylogeny of the widespread North American spotted salamander, Ambystoma maculatum , and test whether phylogenetic patterns conform to regional biogeographical hypotheses about the origins of diversity in eastern North America. Specifically, we address the number and locations of historical refugia, the extent and patterns of postglacial colonization by divergent lineages, and the origin and affinities of populations in the Interior Highland region. Despite apparent morphological uniformity, genetic discontinuities throughout the range of this species suggest that populations were historically fragmented in at least two refugia in the southern Appalachian Mountains. The ranges of these two highly divergent clades expanded northward, resulting in two widely distributed lineages that are sympatric in regions previously proposed as suture zones for other taxa. The evolutionary history of spotted salamander populations underscores the generality of biogeographical processes in eastern North America: despite differences in population size, glacial refugia, and vagility, similar signatures of differentiation are evident among and within widespread taxa.     

Isolation of microsatellite loci from the scleractinian corals,Montastraea cavernosa and Porites astreoides

The ability to assess genetic variation is critical for determining genetic diversity and population structure. In corals, slow evolutionary rates in mitochondrial genomes have left allozymes as the only markers presently available to investigate patterns of intraspecific genetic variation. Characteristics of microsatellites render them more informative than allozymes for such analyses; however, few coral microsatellites are available. This study describes polymorphic microsatellite loci isolated from two scleractinian coral species. Most loci exhibit significant heterozygote deficiencies, likely due to nonrandom mating or Wahlund effects. These markers are being used to investigate gene flow among populations, providing insight into reef connectivity.     

Working at the interface of phylogenetics and population genetics: a biogeographical analysis of Triaenops spp. (Chiroptera: Hipposideridae)

New applications of genetic data to questions of historical biogeography have revolutionized our understanding of how organisms have come to occupy their present distributions. Phylogenetic methods in combination with divergence time estimation can reveal biogeographical centres of origin, differentiate between hypotheses of vicariance and dispersal, and reveal the directionality of dispersal events. Despite their power, however, phylogenetic methods can sometimes yield patterns that are compatible with multiple, equally well-supported biogeographical hypotheses. In such cases, additional approaches must be integrated to differentiate among conflicting dispersal hypotheses. Here, we use a synthetic approach that draws upon the analytical strengths of coalescent and population genetic methods to augment phylogenetic analyses in order to assess the biogeographical history of Madagascar's Triaenops bats (Chiroptera: Hipposideridae). Phylogenetic analyses of mitochondrial DNA sequence data for Malagasy and east African Triaenops reveal a pattern that equally supports two competing hypotheses. While the phylogeny cannot determine whether Africa or Madagascar was the centre of origin for the species investigated, it serves as the essential backbone for the application of coalescent and population genetic methods. From the application of these methods, we conclude that a hypothesis of two independent but unidirectional dispersal events from Africa to Madagascar is best supported by the data.     

Local selection modifies phenotypic divergence among Rana temporaria populations in the presence of gene flow

In ectotherms, variation in life history traits among populations is common and suggests local adaptation. However, geographic variation itself is not a proof for local adaptation, as genetic drift and gene flow may also shape patterns of quantitative variation. We studied local and regional variation in means and phenotypic plasticity of larval life history traits in the common frog Rana temporaria using six populations from central Sweden, breeding in either open‐canopy or partially closed‐canopy ponds. To separate local adaptation from genetic drift, we compared differentiation in quantitative genetic traits (Q_ST) obtained from a common garden experiment with differentiation in presumably neutral microsatellite markers (F_ST). We found that R. temporaria populations differ in means and plasticities of life history traits in different temperatures at local, and in F_ST at regional scale. Comparisons of differentiation in quantitative traits and in molecular markers suggested that natural selection was responsible for the divergence in growth and development rates as well as in temperature‐induced plasticity, indicating local adaptation. However, at low temperature, the role of genetic drift could not be separated from selection. Phenotypes were correlated with forest canopy closure, but not with geographical or genetic distance. These results indicate that local adaptation can evolve in the presence of ongoing gene flow among the populations, and that natural selection is strong in this system.     

Effects of geographic distance,sea barriers and habitat on the genetic structure and diversity of all-hybrid water frog populations

The history of population size and migration patterns leaves its mark in the genetics of populations. We investigate the genetic structure of the edible frog, Pelophylax esculentus in the Danish archipelago and adjacent countries. This frog is of particular interest because it is a hybrid that, in this area, forms all-hybrid populations of diploid (LR) and triploid (LLR and LRR) genomotypes with no (or very few) adults of the parental species (LL and RR). This study is the first to cover the entire geographic range of Danish, Swedish and German all-hybrid populations, documenting their extent and providing a broad picture of their diversity of neutral genetic markers and genomotype proportions. With 18 microsatellite markers, we found that genetic diversity declines northwards in agreement with the glacial refuge and central-marginal hypotheses; however, populations on small and medium-sized islands are no less diverse than those on large islands and continental peninsulas. Isolation by distance exists across the archipelago with limited influence of fragmentation by brackish seawater. The extremely low genetic diversity in all-hybrid populations, compared with adjacent populations, may be responsible for the maintenance of their special breeding system. We also show large variation among ponds in proportions of LLR, LR and LRR genomotypes, but little geographic pattern in their distribution. Instead, we found relationships between the genomotype proportions and some of 15 habitat parameters monitored. Body size differences among LLR, LR and LRR further suggest ecological differences.     

Contemporary and historical evolutionary processes interact to shape patterns of within-lake phenotypic divergences in polyphenic pumpkinseed sunfish, Lepomis gibbosus

Historical and contemporary evolutionary processes can both contribute to patterns of phenotypic variation among populations of a species. Recent studies are revealing how interactions between historical and contemporary processes better explain observed patterns of phenotypic divergence than either process alone. Here, we investigate the roles of evolutionary history and adaptation to current environmental conditions in structuring phenotypic variation among polyphenic populations of sunfish inhabiting 12 postglacial lakes in eastern North America. The pumpkinseed sunfish polyphenism includes sympatric ecomorphs specialized for littoral or pelagic lake habitats. First, we use population genetic methods to test the evolutionary independence of within-lake phenotypic divergences of ecomorphs and to describe patterns of genetic structure among lake populations that clustered into three geographical groupings. We then used multivariate analysis of covariance (MANCOVA) to partition body shape variation (quantified with geometric morphometrics) among the effects of evolutionary history (reflecting phenotypic variation among genetic clusters), the shared phenotypic response of all populations to alternate habitats within lakes (reflecting adaptation to contemporary conditions), and unique phenotypic responses to habitats within lakes nested within genetic clusters. All effects had a significant influence on body form, but the effects of history and the interaction between history and contemporary habitat were larger than contemporary processes in structuring phenotypic variation. This highlights how divergence can be better understood against a known backdrop of evolutionary history.     

The roles of pre‐ and post‐hatching growth rates in generating a latitudinal cline of body size in the eastern fence lizard (Sceloporus undulatus)

Countergradient variation in norms of reaction can dampen the direct effects of environmental influences on phenotypic traits, allowing phenotypic similarity among populations despite exposure to different environmental conditions. Such norms of reaction may occur at any phase of the life‐history (e.g. growth rates during both embryonic and postembryonic stages may influence geographical variation in adult body size). We collected gravid female lizards (Sceloporus undulatus) from northern (Indiana), central (Mississippi), and southern (Florida) populations, spanning almost the full latitudinal range of the species. Adult females from the southern population were smaller. Intrinsic growth rates of hatchlings were higher for the central population than for the other two populations. This pattern does not parallel the countergradient variation previously found in embryonic developmental rates among these populations. Earlier hatching enhanced survival rates of juveniles to a similar degree among populations, although juvenile survival rates in the field generally increase with latitude in this species. Our data reveal geographical variation in the ways in which intrinsic developmental/growth rates and survival shift during ontogeny, and suggest that latitudinal patterns in adult body size (such as Bergmann's rule) can result from both faster growth, and longer periods of growth. © 2012 The Linnean Society of London, Biological Journal of the Linnean Society, 2012, 106 , 202–209.     

Sex-biased evolutionary forces shape genomic patterns of human diversity

Comparisons of levels of variability on the autosomes and X chromosome can be used to test hypotheses about factors influencing patterns of genomic variation. While a tremendous amount of nucleotide sequence data from across the genome is now available for multiple human populations, there has been no systematic effort to examine relative levels of neutral polymorphism on the X chromosome versus autosomes. We analyzed ~210 kb of DNA sequencing data representing 40 independent noncoding regions on the autosomes and X chromosome from each of 90 humans from six geographically diverse populations. We correct for differences in mutation rates between males and females by considering the ratio of within-human diversity to human-orangutan divergence. We find that relative levels of genetic variation are higher than expected on the X chromosome in all six human populations. We test a number of alternative hypotheses to explain the excess polymorphism on the X chromosome, including models of background selection, changes in population size, and sex-specific migration in a structured population. While each of these processes may have a small effect on the relative ratio of X-linked to autosomal diversity, our results point to a systematic difference between the sexes in the variance in reproductive success; namely, the widespread effects of polygyny in human populations. We conclude that factors leading to a lower male versus female effective population size must be considered as important demographic variables in efforts to construct models of human demographic history and for understanding the forces shaping patterns of human genomic variability.     

Genetic diversity and differentiation at MHC genes in island populations of tuatara (Sphenodon spp.)

Neutral genetic markers are commonly used to understand the effects of fragmentation and population bottlenecks on genetic variation in threatened species. Although neutral markers are useful for inferring population history, the analysis of functional genes is required to determine the significance of any observed geographical differences in variation. The genes of the major histocompatibility complex (MHC) are well‐known examples of genes of adaptive significance and are particularly relevant to conservation because of their role in pathogen resistance. In this study, we survey diversity at MHC class I loci across a range of tuatara populations. We compare the levels of MHC variation with that observed at neutral microsatellite markers to determine the relative roles of balancing selection, diversifying selection and genetic drift in shaping patterns of MHC variation in isolated populations. In general, levels of MHC variation within tuatara populations are concordant with microsatellite variation. Tuatara populations are highly differentiated at MHC genes, particularly between the northern and Cook Strait regions, and a trend towards diversifying selection across populations was observed. However, overall our results indicate that population bottlenecks and isolation have a larger influence on patterns of MHC variation in tuatara populations than selection.