The utility of single nucleotide polymorphisms in inferences of population history

Power for detecting genetic divergence: differences between statistical methods and marker loci

Information on statistical power is critical when planning investigations and evaluating empirical data, but actual power estimates are rarely presented in population genetic studies. We used computer simulations to assess and evaluate power when testing for genetic differentiation at multiple loci through combining test statistics or P values obtained by four different statistical approaches, viz. Pearson's chi-square, the log-likelihood ratio G-test, Fisher's exact test, and an F(ST)-based permutation test. Factors considered in the comparisons include the number of samples, their size, and the number and type of genetic marker loci. It is shown that power for detecting divergence may be substantial for frequently used sample sizes and sets of markers, also at quite low levels of differentiation. The choice of statistical method may be critical, though. For multi-allelic loci such as microsatellites, combining exact P values using Fisher's method is robust and generally provides a high resolving power. In contrast, for few-allele loci (e.g. allozymes and single nucleotide polymorphisms) and when making pairwise sample comparisons, this approach may yield a remarkably low power. In such situations chi-square typically represents a better alternative. The G-test without Williams's correction frequently tends to provide an unduly high proportion of false significances, and results from this test should be interpreted with great care. Our results are not confined to population genetic analyses but applicable to contingency testing in general.

Single nucleotide polymorphism (SNP) discovery in duplicated genomes: intron-primed exon-crossing (IPEC) as a strategy for avoiding amplification of duplicated loci in Atlantic salmon (Salmo salar) and other salmonid fishes

Background

Single nucleotide polymorphisms (SNPs) represent the most abundant type of DNA variation in the vertebrate genome, and their applications as genetic markers in numerous studies of molecular ecology and conservation of natural populations are emerging. Recent large-scale sequencing projects in several fish species have provided a vast amount of data in public databases, which can be utilized in novel SNP discovery in salmonids. However, the suggested duplicated nature of the salmonid genome may hamper SNP characterization if the primers designed in conserved gene regions amplify multiple loci.

Results

Here we introduce a new intron-primed exon-crossing (IPEC) method in an attempt to overcome this duplication problem, and also evaluate different priming methods for SNP discovery in Atlantic salmon (Salmo salar) and other salmonids. A total of 69 loci with differing priming strategies were screened in S. salar, and 27 of these produced ~13 kb of high-quality sequence data consisting of 19 SNPs or indels (one per 680 bp). The SNP frequency and the overall nucleotide diversity (3.99 × 10^-4) in S. salar was lower than reported in a majority of other organisms, which may suggest a relative young population history for Atlantic salmon. A subset of primers used in cross-species analyses revealed considerable variation in the SNP frequencies and nucleotide diversities in other salmonids.

Conclusion

Sequencing success was significantly higher with the new IPEC primers; thus the total number of loci to screen in order to identify one potential polymorphic site was six times less with this new strategy. Given that duplication may hamper SNP discovery in some species, the IPEC method reported here is an alternative way of identifying novel polymorphisms in such cases.

Apparent ‘sympatric’ speciation in ecologically similar herbivorous beetles facilitated by multiple colonizations of an island

Coexistence of recently diverged and ecologically similar sister species in complete sympatry represents a particularly compelling case for sympatric speciation. This study investigates the possible sympatric origin of two coexisting bark beetle taxa that utilize the same host plant on the island of La Palma in the Canary Islands. Aphanarthrum subglabrum and Aphanarthrum glabrum ssp. nudum breed inside dead twigs of Euphorbia lamarckii plants and are closely related to the allopatric A. glabrum ssp. glabrum in Tenerife, El Hierro and La Gomera. We tested the various speciation hypotheses in a genealogical context, using mitochondrial gene fragments from Cytochrome Oxidase I and 16S, and nuclear gene fragments from Enolase, Elongation Factor 1alpha and Histone H3. Phylogenetic analyses of the combined nuclear DNA data strongly supported a sister relationship between two sympatric and reproductively isolated taxa in La Palma. However, network analyses of subdivided nonrecombinant segments of the Enolase locus indicated a closer relationship between the two allopatric A. glabrum subspecies, suggesting multiple colonizations of this island. A bimodal distribution of mtDNA haplotypes in La Palma further documented the independent colonization of this island, with asymmetric introgression of mtDNA between two lineages. Consequently, the sympatric origin of the La Palma species is concluded to have involved allopatric phases before the parallel colonization of this island and subsequent introgression at some loci. The clear genetic and morphological evidence for reproductive isolation between these species suggests that the sympatric completion of divergence was either due to initial genetic incompatibility, morphological character displacement in male genitalia, or a combination of these factors.

Introgression from Lepus europaeus to L. timidus in Russia revealed by mitochondrial single nucleotide polymorphisms and nuclear microsatellites

Hybridisation among wild mammal populations may lead to introgression of genes and genomes over the species barrier. In Sweden, in northern Europe, and on the Iberian Peninsula in southern Europe, mitochondrial DNA from L. timidus occurs among L. europaeus specimens, presumably as a result of interspecific hybridisation. In Russia, the species are believed to hybridise as well, but no investigations have confirmed introgression. Here we develop species diagnostic single nucleotide polymorphisms in the mitochondrial genomes and combine them with analysis of nuclear microsatellite markers to investigate hybridisation and introgression in 71 Lepus specimens from Russia. A total of 58 specimens are typical representatives of either species. An additional nine specimens have slightly intermediate genotypes, potentially as a result of introgression of nuclear genes. Finally, we find three specimens with L. europaeus mitochondrial genome and apparent L. timidus nuclear genome. This indicates that the reciprocal transfer of mtDNA occur among Russian populations of these species. Our observation differs from previous observations of mtDNA introgression in Sweden and Iberia, and provides further support for a reticulated mode of introgression within the genus Lepus.

Use of sequence data from rainbow trout and Atlantic salmon for SNP detection in Pacific salmon

Single nucleotide polymorphisms (SNPs) are a class of genetic markers that are well suited to a broad range of research and management applications. Although advances in genotyping chemistries and analysis methods continue to increase the potential advantages of using SNPs to address molecular ecological questions, the scarcity of available DNA sequence data for most species has limited marker development. As the number and diversity of species being targeted for large-scale sequencing has increased, so has the potential for using sequence from sister taxa for marker development in species of interest. We evaluated the use of Oncorhynchus mykiss and Salmo salar sequence data to identify SNPs in three other species (Oncorhynchus tshawytscha, Oncorhynchus nerka and Oncorhynchus keta). Primers designed based on O. mykiss and S. salar alignments were more successful than primers designed based on Oncorhynchus-only alignments for sequencing target species, presumably due to the much larger number of potential targets available from the former alignments and possibly greater sequence conservation in those targets. In sequencing approximately 89 kb we observed a frequency of 4.30 x 10(-3) SNPs per base pair. Approximately half (53/101) of the subsequently designed validation assays resulted in high-throughput SNP genotyping markers. We speculate that this relatively low conversion rate may reflect the duplicated nature of the salmon genome. Our results suggest that a large number of SNPs could be developed for Pacific salmon using sequence data from other species. While the costs of DNA sequencing are still significant, these must be compared to the costs of using other marker classes for a given application.

Evidence for a large-scale population structure of Arabidopsis thaliana from genome-wide single nucleotide polymorphism markers

Population-based methods for the genetic mapping of adaptive traits and the analysis of natural selection require that the population structure and demographic history of a species are taken into account. We characterized geographic patterns of genetic variation in the model plant Arabidopsis thaliana by genotyping 115 genome-wide single nucleotide polymorphism (SNP) markers in 351 accessions from the whole species range using a matrix-assisted laser desorption/ionization time-of-flight assay, and by sequencing of nine unlinked short genomic regions in a subset of 64 accessions. The observed frequency distribution of SNPs is not consistent with a constant-size neutral model of sequence polymorphism due to an excess of rare polymorphisms. There is evidence for a significant population structure as indicated by differences in genetic diversity between geographic regions. Accessions from Central Asia have a low level of polymorphism and an increased level of genome-wide linkage disequilibrium (LD) relative to accessions from the Iberian Peninsula and Central Europe. Cluster analysis with the structure program grouped Eurasian accessions into K = 6 clusters. Accessions from the Iberian Peninsula and from Central Asia constitute distinct populations, whereas Central and Eastern European accessions represent admixed populations in which genomes were reshuffled by historical recombination events. These patterns likely result from a rapid postglacial recolonization of Eurasia from glacial refugial populations. Our analyses suggest that mapping populations for association or LD mapping should be chosen from regional rather than a species-wide sample or identified genetically as sets of individuals with similar average genetic distances.

Power and Sample Size for Testing Associations of Haplotypes with Complex Traits

Evaluation of the association of haplotypes with either quantitative traits or disease status is common practice, and under some situations provides greater power than the evaluation of individual marker loci. The focus on haplotype analyses will increase as more single nucleotide polymorphisms (SNPs) are discovered, either because of interest in candidate gene regions, or because of interest in genome-wide association studies. However, there is little guidance on the determination of the sample size needed to achieve the desired power for a study, particularly when linkage phase of the haplotypes is unknown, and when a subset of tag-SNP markers is measured. There is a growing wealth of information on the distribution of haplotypes in different populations, and it is not unusual for investigators to measure genetic markers in pilot studies in order to gain knowledge of the distribution of haplotypes in the target population. Starting with this basic information on the distribution of haplotypes, we derive analytic methods to determine sample size or power to test the association of haplotypes with either a quantitative trait or disease status (e.g., a case-control study design), assuming that all subjects are unrelated. Our derivations cover both phase-known and phase-unknown haplotypes, allowing evaluation of the loss of efficiency due to unknown phase. We also extend our methods to when a subset of tag-SNPs is chosen, allowing investigators to explore the impact of tag-SNPs on power. Simulations illustrate that the theoretical power predictions are quite accurate over a broad range of conditions. Our theoretical formulae should provide useful guidance when planning haplotype association studies.

The power of single-nucleotide polymorphisms for large-scale parentage inference

Likelihood-based parentage inference depends on the distribution of a likelihood-ratio statistic, which, in most cases of interest, cannot be exactly determined, but only approximated by Monte Carlo simulation. We provide importance-sampling algorithms for efficiently approximating very small tail probabilities in the distribution of the likelihood-ratio statistic. These importance-sampling methods allow the estimation of small false-positive rates and hence permit likelihood-based inference of parentage in large studies involving a great number of potential parents and many potential offspring. We investigate the performance of these importance-sampling algorithms in the context of parentage inference using single-nucleotide polymorphism (SNP) data and find that they may accelerate the computation of tail probabilities >1 millionfold. We subsequently use the importance-sampling algorithms to calculate the power available with SNPs for large-scale parentage studies, paying particular attention to the effect of genotyping errors and the occurrence of related individuals among the members of the putative mother-father-offspring trios. These simulations show that 60-100 SNPs may allow accurate pedigree reconstruction, even in situations involving thousands of potential mothers, fathers, and offspring. In addition, we compare the power of exclusion-based parentage inference to that of the likelihood-based method. Likelihood-based inference is much more powerful under many conditions; exclusion-based inference would require 40% more SNP loci to achieve the same accuracy as the likelihood-based approach in one common scenario. Our results demonstrate that SNPs are a powerful tool for parentage inference in large managed and/or natural populations.

The importance of immune gene variability (MHC) in evolutionary ecology and conservation

Genetic studies have typically inferred the effects of human impact by documenting patterns of genetic differentiation and levels of genetic diversity among potentially isolated populations using selective neutral markers such as mitochondrial control region sequences, microsatellites or single nucleotide polymorphism (SNPs). However, evolutionary relevant and adaptive processes within and between populations can only be reflected by coding genes. In vertebrates, growing evidence suggests that genetic diversity is particularly important at the level of the major histocompatibility complex (MHC). MHC variants influence many important biological traits, including immune recognition, susceptibility to infectious and autoimmune diseases, individual odours, mating preferences, kin recognition, cooperation and pregnancy outcome. These diverse functions and characteristics place genes of the MHC among the best candidates for studies of mechanisms and significance of molecular adaptation in vertebrates. MHC variability is believed to be maintained by pathogen-driven selection, mediated either through heterozygote advantage or frequency-dependent selection. Up to now, most of our knowledge has derived from studies in humans or from model organisms under experimental, laboratory conditions. Empirical support for selective mechanisms in free-ranging animal populations in their natural environment is rare. In this review, I first introduce general information about the structure and function of MHC genes, as well as current hypotheses and concepts concerning the role of selection in the maintenance of MHC polymorphism. The evolutionary forces acting on the genetic diversity in coding and non-coding markers are compared. Then, I summarise empirical support for the functional importance of MHC variability in parasite resistance with emphasis on the evidence derived from free-ranging animal populations investigated in their natural habitat. Finally, I discuss the importance of adaptive genetic variability with respect to human impact and conservation, and implications for future studies.

A comparison of variation between a MHC pseudogene and microsatellite loci of the little greenbul (Andropadus virens)

Background

We investigated genetic variation of a major histcompatibility complex (MHC) pseudogene (Anvi-DAB1) in the little greenbul (Andropadus virens) from four localities in Cameroon and one in Ivory Coast, West Africa. Previous microsatellite and mitochondrial DNA analyses had revealed little or no genetic differentiation among Cameroon localities but significant differentiation between localities in Cameroon and Ivory Coast.

Results

Levels of genetic variation, heterozygosity, and allelic diversity were high for the MHC pseudogene in Cameroon. Nucleotide diversity of the MHC pseudogene in Cameroon and Ivory Coast was comparable to levels observed in other avian species that have been studied for variation in nuclear genes. An excess of rare variants for the MHC pseudogene was found in the Cameroon population, but this excess was not statistically significant. Pairwise measures of population differentiation revealed high divergence between Cameroon and Ivory Coast for microsatellites and the MHC locus, although for the latter distance measures were much higher than the comparable microsatellite distances.

Conclusion

We provide the first ever comparison of variation in a putative MHC pseudogene to variation in neutral loci in a passerine bird. Our results are consistence with the action of neutral processes on the pseudogene and suggest they can provide an independent perspective on demographic history and population substructure.

The Isthmus of Panama: a major physical barrier to gene flow in a highly mobile pantropical seabird

To further test the hypothesis that the Isthmus of Panama is a major barrier to gene flow in pantropical seabirds, we applied phylogeographic methods to mitochondrial control sequence variation in masked booby (Sula dactylatra) populations on either side of the Isthmus of Panama and the southern tip of Africa. In accord with Steeves et al. (2003), we found that all Caribbean masked boobies with the 'secondary contact' cytochrome b haplotype (m-B) shared a control region haplotype (Sd_100), which grouped with Indian-Pacific haplotypes and not Caribbean-Atlantic haplotypes. In addition, Sd_100 was more closely related to control region haplotypes in the Indian Ocean than in the Pacific. We also found that the 'secondary contact' birds diverged more recently from extant populations in the Indian Ocean than in the Pacific. Thus, it appears that these masked boobies did not breach the Isthmus of Panama. Rather, birds likely dispersed around the southern tip of Africa during favourable oceanographic conditions in the Pleistocene.

Phylogenetics of modern birds in the era of genomics

In the 14 years since the first higher-level bird phylogenies based on DNA sequence data, avian phylogenetics has witnessed the advent and maturation of the genomics era, the completion of the chicken genome and a suite of technologies that promise to add considerably to the agenda of avian phylogenetics. In this review, we summarize current approaches and data characteristics of recent higher-level bird studies and suggest a number of as yet untested molecular and analytical approaches for the unfolding tree of life for birds. A variety of comparative genomics strategies, including adoption of objective quality scores for sequence data, analysis of contiguous DNA sequences provided by large-insert genomic libraries, and the systematic use of retroposon insertions and other rare genomic changes all promise an integrated phylogenetics that is solidly grounded in genome evolution. The avian genome is an excellent testing ground for such approaches because of the more balanced representation of single-copy and repetitive DNA regions than in mammals. Although comparative genomics has a number of obvious uses in avian phylogenetics, its application to large numbers of taxa poses a number of methodological and infrastructural challenges, and can be greatly facilitated by a 'community genomics' approach in which the modest sequencing throughputs of single PI laboratories are pooled to produce larger, complementary datasets. Although the polymerase chain reaction era of avian phylogenetics is far from complete, the comparative genomics era-with its ability to vastly increase the number and type of molecular characters and to provide a genomic context for these characters-will usher in a host of new perspectives and opportunities for integrating genome evolution and avian phylogenetics.

Range-wide molecular analysis of the western pond turtle (Emys marmorata): cryptic variation, isolation by distance, and their conservation implications

We analysed phylogeography and population genetic variation across the range of the western pond turtle (Emys marmorata) using rapidly evolving mitochondrial and nuclear DNA sequence data. Nuclear DNA sequences from two unlinked introns displayed extremely low levels of variation, but phylogenetic analyses based on mtDNA recovered four well-supported and geographically coherent clades. These included a large Northern clade composed of populations from Washington south to San Luis Obispo County, California, west of the Coast Ranges; a San Joaquin Valley clade from the southern Great Central Valley; a geographically restricted Santa Barbara clade from a limited region in Santa Barbara and Ventura counties; and a Southern clade that occurs south of the Tehachapi Mountains and west of the Transverse Range south to Baja California, Mexico. An analysis of molecular variance (amova) based on regional hydrographic units revealed that populations from the Sacramento Valley north to Washington were virtually invariant, with no evidence of population substructure among northern river drainage basins. In other areas, E. marmorata contains considerable unrecognized variation, particularly in central and southern California and in northern Baja California, Mexico. Our northern clade is congruent with the distribution of the subspecies Emys marmorata marmorata (Washington-central California). However, no clade is congruent with the distribution of the southern subspecies Emys marmorata pallida from central California-Baja. Thus, recognition of the current subspecies split is not warranted, based on the available genetic evidence. Our amova and phylogenetic results, in conjunction with a growing comparative database for other codistributed aquatic taxa, confirm the occurrence of genetic breaks across the Tehachapi Mountains and Transverse Range bounding the southern end of the Great Central Valley, and point to southern California as a rich source of cryptic genetic variation.

Sequence polymorphism in polyploid wheat and their d-genome diploid ancestor

Sequencing was used to investigate the origin of the D genome of the allopolyploid species Triticum aestivum and Aegilops cylindrica. A 247-bp region of the wheat D-genome Xwye838 locus, encoding ADP-glucopyrophosphorylase, and a 326-bp region of the wheat D-genome Gss locus, encoding granule-bound starch synthase, were sequenced in a total 564 lines of hexaploid wheat (T. aestivum, genome AABBDD) involving all its subspecies and 203 lines of Aegilops tauschii, the diploid source of the wheat D genome. In Ae. tauschii, two SNP variants were detected at the Xwye838 locus and 11 haplotypes at the Gss locus. Two haplotypes with contrasting frequencies were found at each locus in wheat. Both wheat Xwye838 variants, but only one of the Gss haplotypes seen in wheat, were found among the Ae. tauschii lines. The other wheat Gss haplotype was not found in either Ae. tauschii or 70 lines of tetraploid Ae. cylindrica (genomes CCDD), which is known to hybridize with wheat. It is concluded that both T. aestivum and Ae. cylindrica originated recurrently, with at least two genetically distinct progenitors contributing to the formation of the D genome in both species.

A multidimensional approach for detecting species patterns in Malagasy vertebrates

The biodiversity of Madagascar is extraordinarily distinctive, diverse, and endangered. It is therefore urgent that steps be taken to document, describe, interpret, and protect this exceptional biota. As a collaborative group of field and laboratory biologists, we employ a suite of methodological and analytical tools to investigate the vertebrate portion of Madagascar's fauna. Given that species are the fundamental unit of evolution, where micro- and macroevolutionary forces converge to generate biological diversity, a thorough understanding of species distribution and abundance is critical for understanding the evolutionary, ecological, and biogeographic forces that have shaped Malagasy vertebrate diversity. We illustrate the means by which we apply Mayr's "three basic tasks" of the systematist [Mayr, E. (1942) Systematics and the Origin of Species from the Viewpoint of a Zoologist (Harvard Univ. Press, Cambridge, MA)] to identify, classify, and study the organisms that together constitute Madagascar's vertebrate community. Using field inventory methods, specimen-based studies, and morphological and molecular analyses, we formulate hypotheses of species identity that then serve as the foundation for subsequent studies of biology and history. Our experience, as well as that of other investigators, has shown that much of the vertebrate species diversity in Madagascar is "cryptic" for both biological and practical reasons. Beyond issues of cryptic biological diversity, the resolution of species identity in Madagascar has been hampered because of a lack of vouchered comparative material at the population level. Through our activities, we are attempting to remedy these limitations while simultaneously enhancing research capacity in Madagascar.

SNPs in ecological and conservation studies: a test in the Scandinavian wolf population

Single nucleotide polymorphisms (SNPs) have the potential to become the genetic marker of choice in studies of the ecology and conservation of natural populations because of their capacity to access variability across the genome. In this study, we provide one of the first demonstrations of SNP discovery in a wild population in order to address typical issues of importance in ecology and conservation in the recolonized Scandinavian and neighbouring Finnish wolf Canis lupus populations. Using end sequence from BAC (bacterial artificial chromosome) clones specific for dogs, we designed assays for 24 SNP loci, 20 sites of which had previously been shown to be polymorphic in domestic dogs and four sites were newly identified as polymorphic in wolves. Of the 24 assayed loci, 22 SNPs were found to be variable within the Scandinavian population and, importantly, these were able to distinguish individual wolves from one another (unbiased probability of identity of 4.33 x 10(-8)), providing equivalent results to that derived from 12 variable microsatellites genotyped in the same population. An assignment test shows differentiation between the Scandinavian and neighbouring Finnish wolf populations, although not all known immigrants are accurately identified. An exploration of the misclassification rates in the identification of relationships shows that neither 22 SNP nor 20 microsatellite loci are able to discriminate across single order relationships. Despite the remaining obstacle of SNP discovery in nonmodel organisms, the use of SNPs in ecological and conservation studies is encouraged by the advent of large scale screening methods. Furthermore, the ability to amplify extremely small fragments makes SNPs of particular use for population monitoring, where faecal and other noninvasive samples are routinely used.

Testing nested phylogenetic and phylogeographic hypotheses in the Plethodon vandykei species group

Mesic forests in the North American Pacific Northwest occur in two disjunct areas: along the coastal and Cascade ranges of Oregon, Washington, and British Columbia as well as the Northern Rocky Mountains of Idaho, Montana, and British Columbia. Over 150 species or species complexes have disjunct populations in each area, and a priori hypotheses based on phytogeography and geology potentially explain the disjunction via either dispersal or vicariance. Here, we test these hypotheses in the disjunct salamander complex Plethodon vandykei and P. idahoensisby collecting genetic data (669 bp of Cyt b) from 262 individuals. Maximum likelihood analysis indicated reciprocal monophyly of these species, supporting the ancient vicariance hypothesis, whereas parametric bootstrap and Bayesian hypothesis testing allow rejection of the dispersal hypothesis. The coalescent estimate of the time since population divergence (estimated using MDIV) is 3.75 x 106 years, and the 95%credibility interval of this value overlaps with the geological estimate of vicariance, but not the hypothesized dispersal. These results are congruent with the pattern seen in other mesic forest amphibian lineages and suggest disjunction in amphibians may be a concerted response to a geological/climatological event. WithinP. idahoensis, we tested the corollary hypothesis of an inland Pleistocene refugium in the Clearwater drainage with nested clade analysis and coalescent estimates of population growth rate (g). Both analyses support post-Pleistocene expansion from the Clearwater refugium. We corroborated this result by calculating Tajima's Dand mismatch distribution within each drainage, showing strong evidence for recent population expansion within most drainages. This work demonstrates the utility of statistical phylogeography and contributes two novel analytical tools: tests of stationarity with respect to topology in the Bayesian estimation, and the use of coalescent simulations to test the significance of the population growth-rate parameter.

Distribution of genetic variation in the growth hormone 1 gene in Atlantic salmon (Salmo salar) populations from Europe and North America

The level and hierarchical distribution of genetic variation in complete sequences of the Atlantic salmon (Salmo salar) growth hormone (GH1) gene were investigated in populations from Europe and North America with a view to inferring the major evolutionary forces affecting genetic variation at this locus. Seventeen polymorphic sites were identified in complete sequences from nine populations, with levels of noncoding (intron and untranslated region sequences) nucleotide diversity being similar to those observed in other species. No variation, however, was observed in exonic sequences, indicating that nucleotide diversity in the Atlantic salmon GH1 gene is three and 25 times less than that estimated for human and Drosophila coding sequences, respectively. This suggests that purifying selection is the predominant contemporary force controlling the molecular evolution of GH1 coding sequences. Comparison of haplotype relationships within and between populations indicated that differentiation between populations from Europe and North America was greater than within-continent comparisons. However, several haplotypes observed in the northernmost European populations were more similar to those observed in North American than to any other haplotypes observed in Europe. This is most likely to be a result of historical, rather than contemporary, gene flow. Neutrality test statistics, such as Tajima's D, were significantly positive in the European populations in which North American-like haplotypes were observed. Although a positive Tajima's D is commonly interpreted as the signal of balancing selection, a more likely explanation in this case is that either historical migration or ascertainment bias, rather than within population local adaptation, has given rise to an excess of intermediate frequency alleles.

Using cDNA and genomic sequences as tools to develop SNP strategies in cassava (Manihot esculenta Crantz)

Single nucleotide polymorphisms (SNP) are the most abundant type of DNA polymorphism found in animal and plant genomes. They provide an important new source of molecular markers that are useful in genetic mapping, map-based positional cloning, quantitative trait locus mapping and the assessment of genetic distances between individuals. Very little is known on the frequency of SNPs in cassava. We have exploited the recently-developed collection of cassava expressed sequence tags (ESTs) to detect SNPs in the five cultivars of cassava used to generate the sequences. The frequency of intra-cultivar and inter-cultivar SNPs after analysis of 111 contigs was one polymorphism per 905 and one per 1,032 bp, respectively; totaling 1 each 509 bp. We have obtained further information on the frequency of SNPs in six cassava cultivars by analysis of 33 amplicons obtained from 3' EST and BAC end sequences. Overall, about 11 kb of DNA sequence was obtained for each cultivar. A total of 186 SNPs (136 and 50 from ESTs and BAC ends, respectively) were identified. Among these, 146 were intra-cultivar polymorphisms, while 80 were inter-cultivar polymorphisms. Thus the total frequency of SNPs was one per 62 bp. This information will help to develop new strategies for EST mapping as well as their association with phenotypic characteristics.

FAST-TRACK: Integrating QTL mapping and genome scans towards the characterization of candidate loci under parallel selection in the lake whitefish (Coregonus clupeaformis)

As natural selection must act on underlying genetic variation, discovering the number and location of loci under the influence of selection is imperative towards understanding adaptive divergence in evolving populations. Studies employing genome scans have hypothesized that the action of divergent selection should reduce gene flow at the genomic locations implicated in adaptation and speciation among natural populations, yet once 'outlier' patterns of variation have been identified the function and role of such loci needs to be confirmed. We integrated adaptive QTL mapping and genomic scans among diverging sympatric pairs of the lake whitefish (Coregonus clupeaformis) species complex in order to test the hypothesis that differentiation between dwarf and normal ecotypes at growth-associated QTL was maintained by directional selection. We found evidence of significantly high levels of molecular divergence among eight growth QTL where two of the strongest candidate loci under the influence of directional selection exhibited parallel reductions of gene flow over multiple populations.

Cytochrome b polymorphisms and population structure of two species of Alouatta (Primates)

We carried out a phylogenetic and population study in Alouatta caraya and Alouatta belzebul based on cytochrome b DNA sequence data. Maximum Parsimony and Median-Joining analyses grouped A. caraya from different localities showing a population structure in accordance with geographic distribution. The relation between A. caraya haplotypes could be explained with respect to the species range in the Cerrado, one of the most ancient morphoclimatic domains of South America, and the Chaco. Conversely, A. belzebul from the Amazonas and Atlantic forests grouped in a paraphyletic arrangement without an evident geographic pattern. Recent geologic events resulting in the separation of A. belzebul might explain why these geographically distant groups shared similar haplotypes and why ancestral polymorphisms might have been maintained in this species. Time of divergence estimates indicated that the splitting of the Alouatta lineage leading to A. caraya occurred some 4.58 MYA while the lineage leading to A. belzebul emerged 4.14 MYA.

The expansion of conservation genetics

The 'crisis discipline' of conservation biology has voraciously incorporated many technologies to speed up and increase the accuracy of conservation decision-making. Genetic approaches to characterizing endangered species or areas that contain endangered species are prime examples of this. Technical advances in areas such as high-throughput sequencing, microsatellite analysis and non-invasive DNA sampling have led to a much-expanded role for genetics in conservation. Such expansion will allow for more precise conservation decisions to be made and, more importantly, will allow conservation genetics to contribute to area- and landscape-based decision-making processes.

Origin and diversification of the endemic Hawaiian tree snails (Achatinellidae: Achatinellinae) based on molecular evidence

Tree snails of the endemic subfamily Achatinellinae comprise a diverse and important component of the Hawaiian fauna. In recent decades anthropogenic impacts have resulted in devastating extinction rates in Hawaiian tree snails. To address long-standing biogeographic, systematic, and evolutionary questions we used cytochrome c oxidase subunit I (COI) gene sequences to reconstruct the phylogeny of 23 extant species spanning the range of the subfamily from five Hawaiian Islands. To investigate family-level relationships, data were analyzed from 11 terrestrial pulmonate families. Although nodal support for monophyly of the endemic Pacific family Achatinellidae and endemic Hawaiian subfamily Achatinellinae was strong, bifurcation order among deeper ingroup nodes was not well-supported by bootstrap resampling. We hypothesize that lineage extinction and rapidity of lineage formation may have rendered evolutionary reconstruction difficult using a standard phylogenetic approach. Use of an optimized evolutionary model, however, improved resolution and recovered three main clades. The diversification pattern inferred contradicts the traditional biogeographic hypothesis of a Maui origin of the achatinelline lineage. Taxa comprising the basal ingroup clade (Achatinella spp.) and seeding lineages for subsequent clades originated on O'ahu. Therefore it appears that the ancestral colonizing species of achatinellines arrived first on O'ahu from an unknown source, and that O'ahu is the Hawaiian origin of the subfamily. Species previously defined by morphological criteria were generally found to be phylogenetically distinct, and the overall colonization pattern follows the island-age progression rule with several instances of generic polyphyly and back-colonization.

Single nucleotide polymorphism (SNP) discovery in mammals: a targeted-gene approach

Single nucleotide polymorphisms (SNPs) have rarely been exploited in nonhuman and nonmodel organism genetic studies. This is due partly to difficulties in finding SNPs in species where little DNA sequence data exist, as well as to a lack of robust and inexpensive genotyping methods. We have explored one SNP discovery method for molecular ecology, evolution, and conservation studies to evaluate the method and its limitations for population genetics in mammals. We made use of 'CATS' (or 'EPIC') primers to screen for novel SNPs in mammals. Most of these primer sets were designed from primates and/or rodents, for amplifying intron regions from conserved genes. We have screened 202 loci in 16 representatives of the major mammalian clades. Polymerase chain reaction (PCR) success correlated with phylogenetic distance from the human and mouse sequences used to design most primers; for example, specific PCR products from primates and the mouse amplified the most consistently and the marsupial and armadillo amplifications were least successful. Approximately 24% (opossum) to 65% (chimpanzee) of primers produced usable PCR product(s) in the mammals tested. Products produced generally high but variable levels of readable sequence and similarity to the expected genes. In a preliminary screen of chimpanzee DNA, 12 SNPs were identified from six (of 11) sequenced regions, yielding a SNP on average every 400 base pairs (bp). Given the progress in genome sequencing, and the large numbers of CATS-like primers published to date, this approach may yield sufficient SNPs per species for population and conservation genetic studies in nonmodel mammals and other organisms.

NEW MITOCHONDRIAL DNA DATA AFFIRM THE IMPORTANCE OF PLEISTOCENE SPECIATION IN NORTH AMERICAN BIRDS

The timing of origin of modern North American bird species in relation to Pleistocene glaciations has long been the topic of significant discussion and disagreement. Recently, Klicka and Zink (1997) and Avise and Walker (1998) enlivened this debate by using calibrated molecular distance values to estimate timing of speciations. Here we use new molecular studies to test their conclusions. Molecular distance values for 39 pairs of proven sister species, 27 of which are based on new data, alter the currently perceived pattern that avian species splits occurred mainly in the Pliocene and early-mid-Pleistocene. Mitochondrial DNA divergence values for this set of taxa showed a skewed distribution pointing toward relatively young speciation times, in contrast to the pattern presented by Klicka and Zink (1997) for 35 sister plus non-sister species pairs. Our pattern was not significantly different from that of Avise and Walker (1998) for "intraspecific phylogroups," some of which are species. We conclude that the entire Pleistocene, including the last two glacial cycles (<250,000 years ago), was important in speciations of modern North American birds. A substantial number of speciations were both initiated and completed in the last 250,000 years. Simultaneously, many taxa began to diverge in the Pleistocene but their speciations are not yet complete (per Avise and Walker 1998). The suggestion that durations of speciations average two million years is probably a substantial overestimate.

Evaluating the performance of likelihood methods for detecting population structure and migration

A plethora of statistical models have recently been developed to estimate components of population genetic history. Very few of these methods, however, have been adequately evaluated for their performance in accurately estimating population genetic parameters of interest. In this paper, we continue a research program of evaluation of population genetic methods through computer simulation. Specifically, we examine the software MIGRATEE-N 1.6.8 and test the accuracy of this software to estimate genetic diversity (Theta), migration rates, and confidence intervals. We simulated nucleotide sequence data under a neutral coalescent model with lengths of 500 bp and 1000 bp, and with three different per site Theta values of (0.00025, 0.0025, 0.025) crossed with four different migration rates (0.0000025, 0.025, 0.25, 2.5) to construct 1000 evolutionary trees per-combination per-sequence-length. We found that while MIGRATEE-N 1.6.8 performs reasonably well in estimating genetic diversity (Theta), it does poorly at estimating migration rates and the confidence intervals associated with them. We recommend researchers use this software with caution under conditions similar to those used in this evaluation.

Use of unphased multilocus genotype data in indirect association studies

It is usually assumed that detection of a disease susceptability gene via marker polymorphisms in linkage disequilibrium with it is facilitated by consideration of marker haplotypes. However, capture of the marker haplotype information requires resolution of gametic phase, and this must usually be inferred statistically. Recently, we questioned the value of the marker haplotype information, and suggested that certain analyses of multivariate marker data, not based on haplotypes explicitly and not requiring resolution of gametic phase, are often more powerful than analyses based on haplotypes. Here, we review this work and assess more carefully the situations in which our conclusions might apply. We also relate these analyses to alternative approaches to haplotype analysis, namely those based on haplotype similarity and those inspired by cladistics.

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.

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.

Establishment of a high-efficiency SNP-based framework marker set for Arabidopsis

The major goal of this project was the establishment of a tool for rapid mapping of new mutations and genotyping in Arabidopsis consisting of at least 100 evenly spaced framework markers. We assembled a single nucleotide polymorphism (SNP)-based marker set consisting of 112 polymorphic sites with average spacing of 1.15 Mbp derived from an SNP database that we recently developed. This information was used to set up efficient SNP detection reactions based on multiplexed primer extension assays. The 112 Columbia (Col-0)/C24 framework markers were used to assemble 18 multiplexed SNaPshot assays with which up to eight separate loci can be genotyped in a single-tube/single-capillary format. In addition, for 110 framework markers matrix-assisted laser desorption/ionization time of flight (MALDI-ToF) assays have been established for high throughput analyses. We demonstrated the usefulness and the robustness of both procedures of this tool by genotyping 48 BC3F1 individuals created between the accessions Col-0 and C24. Subsets of 10-62 of the established markers discriminate between various combinations of the accessions Col-0, C24, Landsberg erecta (Ler), Cape Verdi Islands (Cvi) and Niederzenz (Nd). Using a subset of 17 evenly distributed and established SNP markers that are also polymorphic between Ler and Col-0, we were able to rapidly map a mutant gene (tbr1) to an interval of 2.3 Mbp in an Ler (tbr1) x Col-0 cross.