Locus-specific genetic differentiation at Rw among warfarin-resistant rat (Rattus norvegicus) populations

Accelerated evolution of Vkorc1 in desert rodent species reveals genetic preadaptation to anticoagulant rodenticides

BACKGROUND

Some rodent species living in arid areas show elevated physiological tolerance to anti-vitamin K rodenticides (AVKs), which seems to be due to some unknown selective pressures that rodents may experience in desert habitats. Genes involved in the ϒ-carboxylation of blood coagulation, including vitamin K epoxide reductase complex, subunit 1 (Vkorc1), ϒ-glutamyl-carboxylase (Ggcx) and NAD(P)H quinone one dehydrogenase (Nqo1) are associated with anticoagulant resistance, or some levels of elevated tolerance, in rodents. To detect whether the DNA sequences of the three genes are also under natural selection in the desert rodent species, we analyzed the Vkorc1, Ggcx and Nqo1 genes of the desert rodents and compared them with other rodent species.

RESULTS

We found an accelerated evolutionary rate in Vkorc1 of desert rodents, especially in Mus spretus, Nannospalax galili and Psammomys obesus. By contrast, signals of positive selection were absent for Ggcx and Nqo1 in all species. Mapping the amino acid variations on the VKORC1 protein three-dimensional model suggested most interspecific amino acid variations occur on the outer surface of the VKORC1 pocket, whereas most intraspecific amino acid changes and known AVK resistance mutations occurred on the inner surface and endoplasmic reticulum luminal loop regions. Some desert-species-specific amino acid variations were found on the positions where known resistance mutations occurred, indicating these variations might be related to the elevated physical tolerance to AVKs in desert rodents.

CONCLUSION

The evolution of Vkorc1 has been accelerated in some desert rodent species, indicating genetic preadaptation to anticoagulant rodenticides. Positive selection and relaxed selection have been detected in Psammomys obesus and Nannospalax galili, indicating the two rodent species might also show tolerance to AVKs, which needs further verification. © 2022 Society of Chemical Industry.

A VKORC1-based SNP survey of anticoagulant rodenticide resistance in the house mouse, Norway rat and roof rat in the USA

BACKGROUND

We conducted a vitamin K epoxide reductase subcomponent 1 (Vkorc1)-based nonsynonymous Single Nucleotide Polymorphism (nsSNP) screen with focus on the house mouse (Mus musculus domesticus), but that also considered the Norway rat (Rattus norvegicus) and roof rat (R. rattus) in the USA.

RESULTS

We detected six Vkorc1 nsSNPs underlying the amino-acid polymorphisms Ala21Thr, Trp59Leu, Ile104Val, Val118Leu, Leu128Ser and Tyr139Cys in house mice (average coverage/SNP; n = 182 individuals), two nsSNPs underlying Arg35Pro and Gly46Ser in the Norway rat (n = 93), with the notable absence of Tyr139Cys (n = 179), and one nsSNP underlying Tyr25Phe in the roof rat (n = 27). Inferred resistance frequency is 29.1% for mice (variability of states 0-98.8%), 6.5% (0-33.3%) for the Norway rat, and 39.3% (0-52.6%) for the roof rat based on Tyr25Phe frequencies.

CONCLUSIONS

Resistance detected in the USA in the 1980s likely was the consequence of Vkorc1 mutations in mice (Leu128Ser and Tyr139Cys), Norway rats (Arg35Pro) and roof rats (Tyr25Phe). Patterns of variant sharing between the USA and Europe indicate the importance of convergent evolution and gene flow in spreading resistance. The spread of nsSNPs in mice between continents appears to have been more effective than in Norway rats. We hypothesize that Arg35Pro may have originated in Norway rats in the USA, whereas Tyr139Cys variants originated in Europe. Tyr25Phe is the likely cause for resistance in roof rats. Further genetic testing in the USA is required to close sampling gaps, and population genomic data are needed to study the origin and spread of this adaptive trait.

The Effects of Epistasis and Pleiotropy on Genome-Wide Scans for Adaptive Outlier Loci

With the advent of next-generation sequencing approaches, the search for individual loci underlying local adaptation has become a major enterprise in evolutionary biology. One promising method to identify such loci is to examine genome-wide patterns of differentiation, using an FST-outlier approach. The effects of pleiotropy and epistasis on this approach are not yet known. Here, we model 2 populations of a sexually reproducing, diploid organism with 2 quantitative traits, one of which is involved in local adaptation. We consider genetic architectures with and without pleiotropy and epistasis. We also model neutral marker loci on an explicit genetic map as the 2 populations diverge and apply FST outlier approaches to determine the extent to which quantitative trait loci (QTL) are detectable. Our results show, under a wide range of conditions, that only a small number of QTL are typically responsible for most of the trait divergence between populations, even when inheritance is highly polygenic. We find that the loci making the largest contributions to trait divergence tend to be detectable outliers. These loci also make the largest contributions to within-population genetic variance. The addition of pleiotropy reduces the extent to which quantitative traits can evolve independently but does not reduce the efficacy of outlier scans. The addition of epistasis, however, reduces the mean FST values for causative QTL, making these loci more difficult, but not impossible, to detect in outlier scans.

Low warfarin resistance frequency in Norway rats in two cities in China after 30 years of usage of anticoagulant rodenticides

BACKGROUND

Anticoagulant rodenticides have been widely used in rodent control in China for >30 years and resistant Norway rats have been reported. Mutations in the vitamin K epoxide reductase complex, subunit 1 (Vkorc1) gene can cause anticoagulant resistance in rodents. In this study, we analyzed the Vkorc1 polymorphisms of 681 Norway rats collected in Zhanjiang and Harbin Cities in China from 2008 to 2015 and evaluated the warfarin resistance frequency.

RESULTS

Analysis revealed four mutations, including three not previously reported. Two new synonymous mutations, His68His and Leu105Leu, are not associated with warfarin resistance. One new nonsynonymous mutation, Ala140Thr, was found in Zhanjiang rat samples collected in 3 years with low frequencies (3.3-4.0%) and is probably associated with warfarin resistance. Laboratory resistance tests suggested low warfarin resistance frequencies in rats from Zhanjiang (4.9-17.1%) and Harbin (0-2.5%) Cities.

CONCLUSIONS

Both genetic analysis and laboratory resistance tests suggested low warfarin resistance frequencies in rats from Zhanjiang and Harbin Cities. The alternate usage of first-generation anticoagulant rodenticides (FGARs) and second-generation anticoagulant rodenticides (SGARs) might represent an effective strategy against the development of warfarin resistance in Norway rats in China. © 2018 Society of Chemical Industry.

Towards the identification of the loci of adaptive evolution

1. Establishing the genetic and molecular basis underlying adaptive traits is one of the major goals of evolutionary geneticists in order to understand the connection between genotype and phenotype and elucidate the mechanisms of evolutionary change. Despite considerable effort to address this question, there remain relatively few systems in which the genes shaping adaptations have been identified. 2. Here, we review the experimental tools that have been applied to document the molecular basis underlying evolution in several natural systems, in order to highlight their benefits, limitations and suitability. In most cases, a combination of DNA, RNA and functional methodologies with field experiments will be needed to uncover the genes and mechanisms shaping adaptation in nature.

Testing for beneficial reversal of dominance during salinity shifts in the invasive copepod Eurytemora affinis, and implications for the maintenance of genetic variation

Maintenance of genetic variation at loci under selection has profound implications for adaptation under environmental change. In temporally and spatially varying habitats, non-neutral polymorphism could be maintained by heterozygote advantage across environments (marginal overdominance), which could be greatly increased by beneficial reversal of dominance across conditions. We tested for reversal of dominance and marginal overdominance in salinity tolerance in the saltwater-to-freshwater invading copepod Eurytemora affinis. We compared survival of F1 offspring generated by crossing saline and freshwater inbred lines (between-salinity F1 crosses) relative to within-salinity F1 crosses, across three salinities. We found evidence for both beneficial reversal of dominance and marginal overdominance in salinity tolerance. In support of reversal of dominance, survival of between-salinity F1 crosses was not different from that of freshwater F1 crosses under freshwater conditions and saltwater F1 crosses under saltwater conditions. In support of marginal overdominance, between-salinity F1 crosses exhibited significantly higher survival across salinities relative to both freshwater and saltwater F1 crosses. Our study provides a rare empirical example of complete beneficial reversal of dominance associated with environmental change. This mechanism might be crucial for maintaining genetic variation in salinity tolerance in E. affinis populations, allowing rapid adaptation to salinity changes during habitat invasions.

Effects of freshwater pollution on the genetics of zebra mussels (Dreissena polymorpha) at the molecular and population level

Revealing long-term effects of contaminants on the genetic structure of organisms inhabiting polluted environments should encompass analyses at the population, molecular, and cellular level. Following this concept, we studied the genetic constitution of zebra mussel populations from a polluted (Dp) and reference sites (Cl) at the river Drava, Croatia, and applied microsatellite and DNA damage analyses (Comet assay, micronucleus test (MNT)). Additionally, mussels from both populations were exposed to polluted wastewater in the laboratory for three days, and DNA damage was analyzed to evaluate acclimatization and genetic adaptation of the investigated populations to the polluted environment. The two populations differed in their genetic constitution. Microsatellite analysis suggested that Dp had undergone a genetic bottleneck. Comet assay did not indicate any difference in DNA damage between the two populations, but MNT revealed that Dp had an increased percentage of micronuclei in hemocytes in comparison to Cl. The laboratory experiment revealed that Dp had a lower percentage of tail DNA and a higher percentage of micronuclei than Cl. These differences between populations were possibly caused by an overall decreased fitness of Dp due to genetic drift and by an enhanced DNA repair mechanism due to acclimatization to pollution in the source habitat.

Linkage disequilibrium and genome-wide association mapping in tetraploid wheat (Triticum turgidum L.)

Association mapping is a powerful tool for the identification of quantitative trait loci through the exploitation of the differential decay of linkage disequilibrium (LD) between marker loci and genes of interest in natural and domesticated populations. Using a sample of 230 tetraploid wheat lines (Triticum turgidum ssp), which included naked and hulled accessions, we analysed the pattern of LD considering 26 simple sequence repeats and 970 mostly mapped diversity array technology loci. In addition, to validate the potential for association mapping in durum wheat, we evaluated the same genotypes for plant height, heading date, protein content, and thousand-kernel weight. Molecular and phenotypic data were used to: (i) investigate the genetic and phenotypic diversity; (ii) study the dynamics of LD across the durum wheat genome, by investigating the patterns of LD decay; and (iii) test the potential of our panel to identify marker–trait associations through the analysis of four quantitative traits of major agronomic importance. Moreover, we compared and validated the association mapping results with outlier detection analysis based on population divergence. Overall, in tetraploid wheat, the pattern of LD is extremely population dependent and is related to the domestication and breeding history of durum wheat. Comparing our data with several other studies in wheat, we confirm the position of many major genes and quantitative trait loci for the traits considered. Finally, the analysis of the selection signature represents a very useful complement to validate marker–trait associations.

Mitochondrial DNA phylogeography of the Norway rat

Central Eastern Asia, foremost the area bordering northern China and Mongolia, has been thought to be the geographic region where Norway rats (Rattus norvegicus) have originated. However recent fossil analyses pointed to their origin in southern China. Moreover, whereas analyses of fossils dated the species' origin as ∼1.2–1.6 million years ago (Mya), molecular analyses yielded ∼0.5–2.9 Mya. Here, to study the geographic origin of the Norway rat and its spread across the globe we analyzed new and all published mitochondrial DNA cytochrome-b (cyt-b; N = 156) and D-loop (N = 212) sequences representing wild rats from four continents and select inbred strains. Our results are consistent with an origin of the Norway rat in southern China ∼1.3 Mya, subsequent prehistoric differentiation and spread in China and Asia from an initially weakly structured ancestral population, followed by further spread and differentiation across the globe during historic times. The recent spreading occurred mostly from derived European populations rather than from archaic Asian populations. We trace laboratory strains to wild lineages from Europe and North America and these represent a subset of the diversity of the rat; leaving Asian lineages largely untapped as a resource for biomedical models. By studying rats from Europe we made the observation that mtDNA diversity cannot be interpreted without consideration of pest control and, possibly, the evolution of rodenticide resistance. However, demographic models explored by forward-time simulations cannot fully explain the low mtDNA diversity of European rats and lack of haplotype sharing with their source from Asia. Comprehensive nuclear marker analyses of a larger sample of Norway rats representing the world are needed to better resolve the evolutionary history of wild rats and of laboratory rats, as well as to better understand the evolution of anticoagulant resistance.

On the stability of the Bayenv method in assessing human SNP-environment associations

Background

Phenotypic variation along environmental gradients has been documented among and within many species, and in some cases, genetic variation has been shown to be associated with these gradients. Bayenv is a relatively new method developed to detect patterns of polymorphisms associated with environmental gradients. Using a Bayesian Markov Chain Monte Carlo (MCMC) approach, Bayenv evaluates whether a linear model relating population allele frequencies to environmental variables is more probable than a null model based on observed frequencies of neutral markers. Although this method has been used to detect environmental adaptation in a number of species, including humans, plants, fish, and mosquitoes, stability between independent runs of this MCMC algorithm has not been characterized. In this paper, we explore the variability of results between runs and the factors contributing to it.

Results

Independent runs of the Bayenv program were carried out using genome-wide single-nucleotide polymorphism (SNP) data from samples from 60 worldwide human populations following previous applications of the Bayenv method. To assess factors contributing to the method's stability, we used varying numbers of MCMC iterations and also analyzed a second modified data set that excluded two Siberian populations with extreme climate variables. Between any two runs, correlations between Bayes factors and the overlap of SNPs in the empirical p value tails were surprisingly low. Enrichments of genic versus non-genic SNPs in the empirical tails were more robust than the empirical p values; however, the significance of the enrichments for some environmental variables still varied among runs, contradicting previously published conclusions. Runs with a greater number of MCMC iterations slightly reduced run-to-run variability, and excluding the Siberian populations did not have a large effect on the stability of the runs.

Conclusions

Because of high run-to-run variability, we advise against making conclusions about genome-wide patterns of adaptation based on only one run of the Bayenv algorithm and recommend caution in interpreting previous studies that have used only one run. Moving forward, we suggest carrying out multiple independent runs of Bayenv and averaging Bayes factors between runs to produce more stable and reliable results. With these modifications, future discoveries of environmental adaptation within species using the Bayenv method will be more accurate, interpretable, and easily compared between studies.

Positive selection in the chromosome 16 VKORC1 genomic region has contributed to the variability of anticoagulant response in humans

VKORC1 (vitamin K epoxide reductase complex subunit 1, 16p11.2) is the main genetic determinant of human response to oral anticoagulants of antivitamin K type (AVK). This gene was recently suggested to be a putative target of positive selection in East Asian populations. In this study, we genotyped the HGDP-CEPH Panel for six VKORC1 SNPs and downloaded chromosome 16 genotypes from the HGDP-CEPH database in order to characterize the geographic distribution of footprints of positive selection within and around this locus. A unique VKORC1 haplotype carrying the promoter mutation associated with AVK sensitivity showed especially high frequencies in all the 17 HGDP-CEPH East Asian population samples. VKORC1 and 24 neighboring genes were found to lie in a 505 kb region of strong linkage disequilibrium in these populations. Patterns of allele frequency differentiation and haplotype structure suggest that this genomic region has been submitted to a near complete selective sweep in all East Asian populations and only in this geographic area. The most extreme scores of the different selection tests are found within a smaller 45 kb region that contains VKORC1 and three other genes (BCKDK, MYST1 (KAT8), and PRSS8) with different functions. Because of the strong linkage disequilibrium, it is not possible to determine if VKORC1 or one of the three other genes is the target of this strong positive selection that could explain present-day differences among human populations in AVK dose requirement. Our results show that the extended region surrounding a presumable single target of positive selection should be analyzed for genetic variation in a wide range of genetically diverse populations in order to account for other neighboring and confounding selective events and the hitchhiking effect.

Phenotypic and genetic differentiation in young-of-the-year common sole (Solea solea) at differentially contaminated nursery grounds

Growth-related characters, condition factor, and genetic differentiation were investigated for a single cohort of young-of-the-year (YOY) sole within and among nurseries with differing levels of heavy metals (Cd, Cu and Zn) contamination in the two Charentais Straits, Bay of Biscay, France. Analyses were performed when individuals recruited (May), then after a full summer spent in each nursery (October). Levels of phenotypic and genetic diversity were compared, together with genetic differentiation at a candidate metallothionein (MT) locus and three putatively neutral microsatellite loci. No phenotypic or genetic differentiation was detected among nurseries in May, but significant variation at each phenotypic trait and at the multilocus level in October (P < 0.001). Single locus analysis demonstrated that only the MT locus was significantly differentiated among nurseries, whether corrected for null alleles or not (θ=0.0401 and θ(corr.FreeNA)=0.0326, respectively; P < 0.001). Results indicate that phenotypic differences among YOY sole nurseries present a molecular correlate acting at identical spatio-temporal scales among nurseries, potentially reflecting differential selective pressure among nurseries in response to contamination.

Analysis of vkorc1 polymorphisms in Norway rats using the roof rat as outgroup

Background

Certain mutations in the vitamin K epoxide reductase subcomponent 1 gene (vkorc1) mediate rodent resistance to warfarin and other anticoagulants. Testing for resistance often involves analysis of the vkorc1. However, a genetic test for the roof rat (Rattus rattus) has yet to be developed. Moreover, an available roof rat vkorc1 sequence would enable species identification based on vkorc1 sequence and the evaluation of natural selection on particular vkorc1 polymorphisms in the Norway rat (R. norvegicus).

Results

We report the coding sequence, introns and 5' and 3' termini for the vkorc1 gene of roof rats (R. r. alexandrinus and R. r. frugivorus) from Uganda, Africa. Newly designed PCR primers now enable genetic testing of the roof rat and Norway rat. Only synonymous and noncoding polymorphisms were found in roof rats from Uganda. Both nominal subspecies of roof rats were indistinguishable from each other but were distinct from R. losea and R. flavipectus; however, the roof rat also shares at least three coding sequence polymorphisms with R. losea and R. flavipectus. Many of recently published vkorc1 synonymous and non-synonymous single nucleotide polymorphisms (SNPs) in Norway rats are likely SNPs from roof rats and/or other Rattus species. Tests applied to presumably genuine Norway rat vkorc1 SNPs are consistent with a role for selection in two populations carrying the derived Phe63Cys and Tyr139Cys mutations.

Conclusion

Geographic mapping of vkorc1 SNPs in roof rats should be facilitated by our report. Our assay should be applicable to most species of Rattus, which are intermediate in genetic distance from roof and Norway rats. Vkorc1-mediated resistance due to non-synonymous coding SNPs is not segregating in roof rats from Uganda. By using the roof rat sequence as a reference vkorc1, SNPs now can be assigned to the correct rat species with more confidence. Sampling designs and genotyping strategies employed so far have helped detect candidate mutations underlying vkorc1-mediated resistance, but generally provided unsuitable data to test for selection. We propose that our understanding of vkorc1-mediated evolution of resistance in rodents would benefit from the adoption of sampling and genotyping designs that enable tests for selection on vkorc1.

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.

Distribution of VKORC1 single nucleotide polymorphism in wild Rattus norvegicus in France

BACKGROUND

Anticoagulant rodenticides are commonly used to control rodent pests all over the world. These pesticides inhibit one enzyme of the vitamin K cycle, Vkorc1, and thus prevent blood clotting and cause death by haemorrhage. Resistance to anticoagulants was first observed in Scotland in 1958, and more potent anticoagulants have been developed to overcome this obstacle. Unfortunately, these chemicals are very toxic and cannot be used everywhere. Some authors have shown that resistance to anticoagulants seems closely linked with single nucleotide polymorphism (SNP) in the Vkorc1 gene.

RESULTS

This study draws a map of SNP and haplotypes found in Vkorc1 in rats from different areas of France. Some of them had never been described before. Moreover, the Y139F mutation, described previously in France and Belgium, is the most frequent in France. This mutation is known to be associated with a strong resistance to anticoagulants, and it was found in 28% of the samples.

CONCLUSION

This biomolecular approach to studying and detecting resistance is easier to carry out than the phenotypic approach measuring blood coagulation time because it can be conducted on biological samples from dead animals, and it is less dangerous for the operator.

Introgression from modern hybrid varieties into landrace populations of maize (Zea mays ssp. mays L.) in central Italy

Landraces are domesticated local plant varieties that did not experience a deliberate and intensive selection during a formal breeding programme. In Europe, maize landraces are still cultivated, particularly in marginal areas where traditional farming is often practiced. Here, we have studied the evolution of flint maize landraces from central Italy over 50 years of on-farm cultivation, when dent hybrid varieties were introduced and their use was widespread. We have compared an 'old' collection, obtained during the 1950s, before the introduction of hybrids, and a recent collection of maize landraces. For comparison, a sample of maize landraces from north Italy, and of improved germplasm, including hybrids and inbred lines were also used. A total of 296 genotypes were analysed using 21 microsatellites. Our results show that the maize landraces collected in the last 5-10 years have evolved directly from the flint landrace gene pool cultivated in central Italy before the introduction of modern hybrids. The population structure, diversity and linkage disequilibrium analyses indicate a significant amount of introgression from hybrid varieties into the recent landrace populations. No evidence of genetic erosion of the maize landraces was seen, suggesting that in situ conservation of landraces is an efficient strategy for preserving genetic diversity. Finally, the level of introgression detected was very variable among recent landraces, with most of them showing a low level of introgression; this suggests that coexistence between different types of agriculture is possible, with the adoption of correct practices that are aimed at avoiding introgression from undesired genetic sources.

Signatures of selection in natural populations adapted to chronic pollution

Background

Populations of the teleost fish Fundulus heteroclitus appear to flourish in heavily polluted and geographically separated Superfund sites. Populations from three Superfund sites (New Bedford Harbor, MA, Newark Bay, NJ, and Elizabeth River, VA) have independently evolved adaptive resistance to chemical pollutants. In these polluted populations, natural selection likely has altered allele frequencies of loci that affect fitness or that are linked to these loci. The aim of this study was to identify loci that exhibit non-neutral behavior in the F. heteroclitus genome in polluted populations versus clean reference populations.

Results

To detect signatures of natural selection and thus identify genetic bases for adaptation to anthropogenic stressors, we examined allele frequencies for many hundreds of amplified fragment length polymorphism markers among populations of F. heteroclitus. Specifically, we contrasted populations from three Superfund sites (New Bedford Harbor, MA, Newark Bay, NJ, and Elizabeth River, VA) to clean reference populations flanking the polluted sites. When empirical F_ST values were compared to a simulated distribution of F_ST values, 24 distinct outlier loci were identified among pairwise comparisons of pollutant impacted F. heteroclitus populations and both surrounding reference populations. Upon removal of all outlier loci, there was a strong correlation (R² = 0.79, p < 0.0001) between genetic and geographical distance. This apparently neutral evolutionary pattern was not evident when outlier loci were included (R² = 0.092, p = 0.0721). Two outlier loci were shared between New Bedford Harbor and Elizabeth River populations, and two different loci were shared between Newark Bay and Elizabeth River populations.

Conclusion

In total, 1% to 6% of loci are implicated as being under selection or linked to areas of the genome under selection in three F. heteroclitus populations that reside in polluted estuaries. Shared loci among polluted sites indicate that selection may be acting on multiple loci involved in adaptation, and loci shared between polluted sites potentially are involved in a generalized adaptive response.

Evolutionary ecotoxicology of wild yellow perch (Perca flavescens) populations chronically exposed to a polymetallic gradient

Depending on such factors as the intensity and duration of the exposure, and the genetic diversity and connectedness of the starting population, exposure to elevated metal concentrations can result in population level alterations such as demographic bottlenecks or metal-induced selection. These processes can be revealed using a population genetic approach, and have important implications with respect to population persistence. The main objective of this study was to examine the role of metal contamination in driving evolutionary changes by documenting patterns of genetic diversity within and among populations of wild yellow perch (Perca flavescens) in two major mining regions that have been subjected to metal emissions from smelters for at least 80 years; Rouyn-Noranda, Québec and Sudbury, Ontario. Yellow perch populations from ten lakes representing a gradient of metal contamination in each of the two lake systems were evaluated concurrently to reveal relationships between metal contamination and genetic diversity. These replicated sympatric observations allowed us to evaluate correlations and infer causal relationships between metal exposure and evolutionary responses in this species. Within-population gene diversity over all loci was negatively correlated with liver cadmium contamination (P<0.001; r(2)=0.47). Similarly, a negative correlation between gene diversity and liver copper contamination was observed at a single locus (Pfla L1, P=0.005; r(2)=0.33), suggesting a local effect of copper contamination. Internal relatedness, an index of individual diversity, presented the opposite tendency as the more contaminated individuals were more diverse than were the less contaminated ones in contaminated and reference populations. Our results thus suggest that the selective response to contamination has been large enough to substantially reduce the within-population genetic diversity, despite the fact that the less inbred individuals may be favoured by selection within any given population. Overall, our results reveal that >50 years of metal contamination have significantly impacted patterns of genetic diversity observed among populations of wild yellow perch in mining areas and as such, may have affected the capacity of populations to respond to future environmental changes.

A cardiovascular phenotype in warfarin-resistant Vkorc1 mutant rats

BACKGROUND: The inhibition of the vitamin K cycle by warfarin promotes arterial calcification in the rat. Conceivably, genetically determined vitamin K-deficiency owing to a mutant epoxide reductase subcomponent 1 (Vkorc1) gene, a key component of the vitamin K cycle, might also promote arterial calcification. In the absence of an available Vkorc1 gene knockout model we used a wild-derived Vkorc1 mutant rat strain (Rattus norvegicus) to explore the validity of this hypothesis. METHODS: We provide histopathological descriptions of a naturally occurring Vkorc1 gene knockdown: wild-derived lab-reared rats that are resistant to the anticoagulant warfarin owing to a non-synonymous mutation in the Vkorc1 gene (Vkorc1(Y->C)), which, in vitro, reduces the basal activity of the vitamin K epoxide reductase enzyme complex by ~52%. H&E stained sections of heart and kidney were compared between homozygous Vkorc1(Y->C/ Y->C), heterozygous Vkorc1(Y->C/+) and wildtype Vkorc1(+/+) rats of both sexes. RESULTS: We observed that the aorta of the heart was mineralized in the Vkorc1(Y->C/ Y->C) male rats but lesions were virtually absent from Vkorc1(Y->C/+) and Vkorc1(+/+) male and all female rats. The renal arteries were mineralized in Vkorc1(Y->C/ Y->C) and Vkorc1(Y->C/+) mutant rats, regardless of sex. CONCLUSIONS: Results support a hypothesis that posits that Vkorc1 genetic polymorphisms reducing basal enzyme activity could affect cardiovascular health, with dependencies on genotype, sex, and tissue. The undercarboxylation of the vitamin K-dependent Matrix Gla protein may be the crucial component of the pathway promoting this mineralization.

Combining population genomics and quantitative genetics: finding the genes underlying ecologically important traits

A central challenge in evolutionary biology is to identify genes underlying ecologically important traits and describe the fitness consequences of naturally occurring variation at these loci. To address this goal, several novel approaches have been developed, including 'population genomics,' where a large number of molecular markers are scored in individuals from different environments with the goal of identifying markers showing unusual patterns of variation, potentially due to selection at linked sites. Such approaches are appealing because of (1) the increasing ease of generating large numbers of genetic markers, (2) the ability to scan the genome without measuring phenotypes and (3) the simplicity of sampling individuals without knowledge of their breeding history. Although such approaches are inherently applicable to non-model systems, to date these studies have been limited in their ability to uncover functionally relevant genes. By contrast, quantitative genetics has a rich history, and more recently, quantitative trait locus (QTL) mapping has had some success in identifying genes underlying ecologically relevant variation even in novel systems. QTL mapping, however, requires (1) genetic markers that specifically differentiate parental forms, (2) a focus on a particular measurable phenotype and (3) controlled breeding and maintenance of large numbers of progeny. Here we present current advances and suggest future directions that take advantage of population genomics and quantitative genetic approaches - in both model and non-model systems. Specifically, we discuss advantages and limitations of each method and argue that a combination of the two provides a powerful approach to uncovering the molecular mechanisms responsible for adaptation.

Selective sweeps reveal candidate genes for adaptation to drought and salt tolerance in common sunflower, Helianthus annuus

Here we report the results of an analysis of variation at 128 EST-based microsatellites in wild Helianthus annuus, using populations from the species' typical plains habitat in Kansas and Colorado, as well as two arid desert and two distinct brackish marsh areas in Utah. The test statistics lnRV and lnRH were used to find regions of the genome that were significantly less variable in one population relative to the others and thus are likely to contain genes under selection. A small but detectable percentage (1.5-6%) of genes showed evidence for selection from both statistics in any particular environment, and a total of 17 loci showed evidence of selection in at least one environment. Distance-based measures provided additional evidence of selection for 15 of the 17 loci. Global F(ST)-values were significantly higher for candidate loci, as expected under divergent selection. However, pairwise F(ST)-values were lower for populations that shared a selective sweep. Moreover, while spatially separated populations undergoing similar selective pressures showed evidence of divergence at some loci, they evolved in concert at other loci. Thus, this study illustrates how selective sweeps might contribute both to the integration of conspecific populations and to the differentiation of races or species.

Origin and colonization history of newly virulent strains of the phytopathogenic fungus Venturia inaequalis

Plant resistance based on a gene-for-gene relationship tends not to be durable because virulent strains can emerge through mutation and colonize formerly resistant cultivars over large areas. Our objective was to determine the origin and colonization history of newly virulent strains of the fungus Venturia inaequalis. Four hundred and eighty-one strains from seven multi-cultivar orchards within a 300km wide area were collected from a cultivar carrying the overcome Vf-gene (Vfcv populations) and from non-Vf cultivars (nVfcv populations). Using 9 microsatellite loci, we showed that Vfcv populations were closely related though highly differentiated, indicating they have been recently founded following rare long distance dispersal events from a common origin. Assignment tests revealed that the Vfcv populations had not emerged from any of the nVfcv populations. Vfcv populations had a low diversity and were strongly differentiated from nVfcv populations several years after the first reported breakdown, suggesting reproductive isolation between the two populations.

Warfarin resistance in a French strain of rats

A warfarin-resistant strain and a warfarin-susceptible strain of wild rats (Rattus norvegicus) maintained in enclosures of the National Veterinary School of Lyon (France) were studied to determine the mechanism of the resistance to anticoagulant rodenticides. A low vitamin K epoxide reductase (VKOR) activity has been reported for many resistant rat strains. As recently suggested, mutations in the vitamin K epoxide reductase subunit 1 (VKORC1) gene are the genetic basis of anticoagulant resistance in wild populations of rats from various locations in Europe. Here we report, for our strain, one of the seven described mutations (Tyr139Phe) for VKORC1 in rats. In addition, a low expression of mRNA encoding VKORC1 gene is observed in resistant rats, which could explain their low VKOR activity. We calculated kinetic parameters of VKOR in the warfarin-resistant and warfarin-susceptible rats. The V(max) and the K(m) of the VKOR obtained in resistant rats were lowered by 57 and 77%, respectively, compared to those obtained in susceptible rats. As a consequence, the enzymatic efficiency (V(m)/K(m)) of the VKOR was similar between resistant and susceptible rats. This result could be a good explanation to the observation that no clinical signs of vitamin K deficiency was observed in the warfarin-resistant strain, while a low VKOR activity was found. VKOR activity in warfarin-resistant rats was poorly inhibited by warfarin (K(i) for warfarin is 29 microM and 0.72 microM for resistant and susceptible rats, respectively).

Selective sweep mapping of genes with large phenotypic effects

Many domestic dog breeds have originated through fixation of discrete mutations by intense artificial selection. As a result of this process, markers in the proximity of genes influencing breed-defining traits will have reduced variation (a selective sweep) and will show divergence in allele frequency. Consequently, low-resolution genomic scans can potentially be used to identify regions containing genes that have a major influence on breed-defining traits. We model the process of breed formation and show that the probability of two or three adjacent marker loci showing a spurious signal of selection within at least one breed (i.e., Type I error or false-positive rate) is low if highly variable and moderately spaced markers are utilized. We also use simulations with selection to demonstrate that even a moderately spaced set of highly polymorphic markers (e.g., one every 0.8 cM) has high power to detect regions targeted by strong artificial selection in dogs. Further, we show that a gene responsible for black coat color in the Large Munsterlander has a 40-Mb region surrounding the gene that is very low in heterozygosity for microsatellite markers. Similarly, we survey 302 microsatellite markers in the Dachshund and find three linked monomorphic microsatellite markers all within a 10-Mb region on chromosome 3. This region contains the FGFR3 gene, which is responsible for achondroplasia in humans, but not in dogs. Consequently, our results suggest that the causative mutation is a gene or regulatory region closely linked to FGFR3.

Geographical distribution of selected and putatively neutral SNPs in Southeast Asian malaria parasites

Loci targeted by directional selection are expected to show elevated geographical population structure relative to neutral loci, and a flurry of recent papers have used this rationale to search for genome regions involved in adaptation. Studies of functional mutations that are known to be under selection are particularly useful for assessing the utility of this approach. Antimalarial drug treatment regimes vary considerably between countries in Southeast Asia selecting for local adaptation at parasite loci underlying resistance. We compared the population structure revealed by 10 nonsynonymous mutations (nonsynonymous single-nucleotide polymorphisms [nsSNPs]) in four loci that are known to be involved in antimalarial drug resistance, with patterns revealed by 10 synonymous mutations (synonymous single-nucleotide polymorphisms [sSNPs]) in housekeeping genes or genes of unknown function in 755 Plasmodium falciparum infections collected from 13 populations in six Southeast Asian countries. Allele frequencies at known nsSNPs underlying resistance varied markedly between locations (F(ST) = 0.18-0.66), with the highest frequencies on the Thailand-Burma border and the lowest frequencies in neighboring Lao PDR. In contrast, we found weak but significant geographic structure (F(ST) = 0-0.14) for 8 of 10 sSNPs. Importantly, all 10 nsSNPs showed significantly higher F(ST) (P < 8 x 10(-5)) than simulated neutral expectations based on observed F(ST) values in the putatively neutral sSNPs. This result was unaffected by the methods used to estimate allele frequencies or the number of populations used in the simulations. Given that dense single-nucleotide polymorphism (SNP) maps and rapid SNP assay methods are now available for P. falciparum, comparing genetic differentiation across the genome may provide a valuable aid to identifying parasite loci underlying local adaptation to drug treatment regimes or other selective forces. However, the high proportion of polymorphic sites that appear to be under balancing selection (or linked to selected sites) in the P. falciparum genome violates the central assumption that selected sites are rare, which complicates identification of outlier loci, and suggests that caution is needed when using this approach.

High-diversity genes in the Arabidopsis genome

High-diversity genes represent an important class of loci in organismal genomes. Since elevated levels of nucleotide variation are a key component of the molecular signature for balancing selection or local adaptation, high-diversity genes may represent loci whose alleles are selectively maintained as balanced polymorphisms. Comparison of 4300 random shotgun sequence fragments of the Arabidopsis thaliana Ler ecotype genome with the whole genomic sequence of the Col-0 ecotype identified 60 genes with putatively high levels of intraspecific variability. Eleven of these genes were sequenced in multiple A. thaliana accessions, 3 of which were found to display elevated levels of nucleotide polymorphism. These genes encode the myb-like transcription factor MYB103, a putative soluble starch synthase I, and a homeodomain-leucine zipper transcription factor. Analysis of these genes and 4-7 flanking genes in 14-20 A. thaliana ecotypes revealed that two of these loci show other characteristics of balanced polymorphisms, including broad peaks of nucleotide diversity spanning multiple linked genes and an excess of intermediate-frequency polymorphisms. Scanning genomes for high-diversity genomic regions may be useful in approaches to adaptive trait locus mapping for uncovering candidate balanced polymorphisms.

Using genome scans of DNA polymorphism to infer adaptive population divergence

Elucidating the genetic basis of adaptive population divergence is a goal of central importance in evolutionary biology. In principle, it should be possible to identify chromosomal regions involved in adaptive divergence by screening genome-wide patterns of DNA polymorphism to detect the locus-specific signature of positive directional selection. In the case of spatially separated populations that inhabit different environments or sympatric populations that exploit different ecological niches, it is possible to identify loci that underlie divergently selected traits by comparing relative levels of differentiation among large numbers of unlinked markers. In this review I first address the question of whether diversifying selection on polygenic traits can be expected to produce predictable patterns of allelic variation at the underlying quantitative trait loci (QTL), and whether the locus-specific effects of selection can be reliably detected against the genome-wide backdrop of stochastic variability. I then review different approaches that have been developed to identify loci involved in adaptive population divergence and I discuss the relative merits of model-based approaches that rely on assumptions about population structure vs. model-free approaches that are based on empirical distributions of summary statistics. Finally, I consider the evolutionary and functional insights that might be gained by conducting genome scans for loci involved in adaptive population divergence.

Natural selection drives altitudinal divergence at the albumin locus in deer mice, Peromyscus maniculatus

In populations that are distributed across steep environmental gradients, the potential for local adaptation is largely determined by the spatial scale of fitness variation relative to dispersal distance. Since altitudinal gradients are generally characterized by dramatic ecological transitions over relatively short linear distances, adaptive divergence across such gradients will typically require especially strong selection to counterbalance the homogenizing effect of gene flow. Here we report the results of a study that was designed to test for evidence of adaptive divergence across an altitudinal gradient in a natural population of deer mice, Peromyscus maniculatus. We conducted a multilocus survey of allozyme variation across a steep altitudinal gradient in the southern Rocky Mountains that spanned several distinct biomes, from prairie grassland to alpine tundra. As a control for the effects of altitude, we also surveyed the same loci in mice sampled along a latitudinal transect through the prairie grassland that ran perpendicular to the east-west altitudinal transect. We used a coalescent-based simulation model to identify loci that deviated from neutral expectations, and we then assessed whether locus-specific patterns of variation were nonrandom with respect to altitude. Results indicated that the albumin locus (Alb) reflects a history of diversifying selection across the altitudinal gradient. This conclusion is supported by two main lines of evidence: (1) Alb was characterized by levels of divergence across the altitudinal transect that exceeded neutral expectations in two consecutive years of sampling (in contrast to the spatial pattern of variation across the latitudinal transect), and (2) levels of divergence at the Alb locus exhibited a positive association with altitudinal distance in both years (in contrast to the pattern observed at unlinked loci). We conclude that clinal variation at the Alb locus reflects a balance between gene flow and diversifying selection that results from elevational changes in fitness rankings among alternative genotypes.

Assessing the effect of natural selection in malaria parasites

There are few concepts that have been used across disciplines; one of them is natural selection. The impact that this process has on parasite genetic diversity is reviewed here by discussing examples on drug resistance and vaccine antigens. Emphasis is made on how mechanisms need to be addressed rather than associations, and how such investigations were out of reach of biomedical researchers only a decade ago.

Warfarin and the Vitamin K-Dependent γ-Carboxylation System

Insight into the molecular basis for genetic warfarin resistance has recently been accomplished by the identification of an 18-kDa protein of the endoplasmic reticulum that is targeted by the drug. When expressed in eukaryotic and insect cells, the protein reduces vitamin K1 2,3-epoxide in a warfarin-sensitive reaction. This finding strongly suggests that the protein is part of the vitamin K cycle, which is essential for the production of vitamin K-dependent proteins. Identification of the 18-kDa protein has aided the understanding of the vitamin K-dependent gamma-carboxylation system at the molecular level.

The power and promise of population genomics: from genotyping to genome typing

Population genomics has the potential to improve studies of evolutionary genetics, molecular ecology and conservation biology, by facilitating the identification of adaptive molecular variation and by improving the estimation of important parameters such as population size, migration rates and phylogenetic relationships. There has been much excitement in the recent literature about the identification of adaptive molecular variation using the population-genomic approach. However, the most useful contribution of the genomics model to population genetics will be improving inferences about population demography and evolutionary history.