The utility of single nucleotide polymorphisms in inferences of population history

Whole-genome-based phylogeny and divergence of the genus Brucella

Brucellae are worldwide bacterial pathogens of livestock and wildlife, but phylogenetic reconstructions have been challenging due to limited genetic diversity. We assessed the taxonomic and evolutionary relationships of five Brucella species-Brucella abortus, B. melitensis, B. suis, B. canis, and B. ovis-using whole-genome comparisons. We developed a phylogeny using single nucleotide polymorphisms (SNPs) from 13 genomes and rooted the tree using the closely related soil bacterium and opportunistic human pathogen, Ochrobactrum anthropi. Whole-genome sequencing and a SNP-based approach provided the requisite level of genetic detail to resolve species in the highly conserved brucellae. Comparisons among the Brucella genomes revealed 20,154 orthologous SNPs that were shared in all genomes. Rooting with Ochrobactrum anthropi reveals that the B. ovis lineage is basal to the rest of the Brucella lineage. We found that B. suis is a highly divergent clade with extensive intraspecific genetic diversity. Furthermore, B. suis was determined to be paraphyletic in our analyses, only forming a monophyletic clade when the B. canis genome was included. Using a molecular clock with these data suggests that most Brucella species diverged from their common B. ovis ancestor in the past 86,000 to 296,000 years, which precedes the domestication of their livestock hosts. Detailed knowledge of the Brucella phylogeny will lead to an improved understanding of the ecology, evolutionary history, and host relationships for this genus and can be used for determining appropriate genotyping approaches for rapid detection and diagnostic assays for molecular epidemiological and clinical studies.

EST, COSII, and arbitrary gene markers give similar estimates of nucleotide diversity in cultivated tomato (Solanum lycopersicum L.)

Because cultivated tomato (Solanum lycopersicum L.) is low in genetic diversity, public, verified single nucleotide polymorphism (SNP) markers within the species are in demand. To promote marker development we resequenced approximately 23 kb in a diverse set of 31 tomato lines including TA496. Three classes of markers were sampled: (1) 26 expressed-sequence tag (EST), all of which were predicted to be polymorphic based on TA496, (2) 14 conserved ortholog set II (COSII) or unigene, and (3) ten published sequences, composed of nine fruit quality genes and one anonymous RFLP marker. The latter two types contained mostly noncoding DNA. In total, 154 SNPs and 34 indels were observed. The distributions of nucleotide diversity estimates among marker types were not significantly different from each other. Ascertainment bias of SNPs was evaluated for the EST markers. Despite the fact that the EST markers were developed using SNP prediction within a sample consisting of only one TA496 allele and one additional allele, the majority of polymorphisms in the 26 EST markers were represented among the other 30 tomato lines. Fifteen EST markers with published SNPs were more closely examined for bias. Mean SNP diversity observations were not significantly different between the original discovery sample of two lines (53 SNPs) and the 31 line diversity panel (56 SNPs). Furthermore, TA496 shared its haplotype with at least one other line at 11 of the 15 markers. These data demonstrate that public EST databases and noncoding regions are a valuable source of unbiased SNP markers in tomato.

New technologies for ultra-high throughput genotyping in plants

Molecular genetic markers represent one of the most powerful tools for the analysis of plant genomes and the association of heritable traits with underlying genetic variation. Molecular marker technology has developed rapidly over the last decade, with the development of high-throughput genotyping methods. Two forms of sequence-based marker, simple sequence repeats (SSRs), also known as microsatellites and single nucleotide polymorphisms (SNPs) now predominate applications in modern plant genetic analysis, along the anonymous marker systems such as amplified fragment length polymorphisms (AFLPs) and diversity array technology (DArT). The reducing cost of DNA sequencing and increasing availability of large sequence data sets permits the mining of this data for large numbers of SSRs and SNPs. These may then be used in applications such as genetic linkage analysis and trait mapping, diversity analysis, association studies and marker-assisted selection. Here, we describe automated methods for the discovery of molecular markers and new technologies for high-throughput, low-cost molecular marker genotyping. Genotyping examples include multiplexing of SSRs using Multiplex-Ready marker technology (MRT); DArT genotyping; SNP genotyping using the Invader assay, the single base extension (SBE), oligonucleotide ligation assay (OLA) SNPlex system, and Illumina GoldenGate and Infinium methods.

Mining for SNPs and SSRs using SNPServer, dbSNP and SSR taxonomy tree

Molecular genetic markers represent one of the most powerful tools for the analysis of genomes and the association of heritable traits with underlying genetic variation. The development of high-throughput methods for the detection of single nucleotide polymorphisms (SNPs) and simple sequence repeats (SSRs) has led to a revolution in their use as molecular markers. The availability of large sequence data sets permits mining for these molecular markers, which may then be used for applications such as genetic trait mapping, diversity analysis and marker assisted selection in agriculture. Here we describe web-based automated methods for the discovery of SSRs using SSR taxonomy tree, the discovery of SNPs from sequence data using SNPServer and the identification of validated SNPs from within the dbSNP database. SSR taxonomy tree identifies pre-determined SSR amplification primers for virtually all species represented within the GenBank database. SNPServer uses a redundancy based approach to identify SNPs within DNA sequences. Following submission of a sequence of interest, SNPServer uses BLAST to identify similar sequences, CAP3 to cluster and assemble these sequences and then the SNP discovery software autoSNP to detect SNPs and insertion/deletion (indel) polymorphisms. The NCBI dbSNP database is a catalogue of molecular variation, hosting validated SNPs for several species within a public-domain archive.

Natural genetic variation of Arabidopsis thaliana is geographically structured in the Iberian peninsula

To understand the demographic history of Arabidopsis thaliana within its native geographical range, we have studied its genetic structure in the Iberian Peninsula region. We have analyzed the amount and spatial distribution of A. thaliana genetic variation by genotyping 268 individuals sampled in 100 natural populations from the Iberian Peninsula. Analyses of 175 individuals from 7 of these populations, with 20 chloroplast and nuclear microsatellite loci and 109 common single nucleotide polymorphisms, show significant population differentiation and isolation by distance. In addition, analyses of one genotype from 100 populations detected significant isolation by distance over the entire Iberian Peninsula, as well as among six Iberian subregions. Analyses of these 100 genotypes with different model-based clustering algorithms inferred four genetic clusters, which show a clear-cut geographical differentiation pattern. On the other hand, clustering analysis of a worldwide sample showed a west-east Eurasian longitudinal spatial gradient of the commonest Iberian genetic cluster. These results indicate that A. thaliana genetic variation displays significant regional structure and consistently support the hypothesis that Iberia has been a glacial refugium for A. thaliana. Furthermore, the Iberian geographical structure indicates a complex regional population dynamics, suggesting that this region contained multiple Pleistocene refugia with a different contribution to the postglacial colonization of Europe.

Scanning the genome for gene single nucleotide polymorphisms involved in adaptive population differentiation in white spruce

Conifers are characterized by a large genome size and a rapid decay of linkage disequilibrium, most often within gene limits. Genome scans based on noncoding markers are less likely to detect molecular adaptation linked to genes in these species. In this study, we assessed the effectiveness of a genome-wide single nucleotide polymorphism (SNP) scan focused on expressed genes in detecting local adaptation in a conifer species. Samples were collected from six natural populations of white spruce (Picea glauca) moderately differentiated for several quantitative characters. A total of 534 SNPs representing 345 expressed genes were analysed. Genes potentially under natural selection were identified by estimating the differentiation in SNP frequencies among populations (F(ST)) and identifying outliers, and by estimating local differentiation using a Bayesian approach. Both average expected heterozygosity and population differentiation estimates (H(E) = 0.270 and F(ST) = 0.006) were comparable to those obtained with other genetic markers. Of all genes, 5.5% were identified as outliers with F(ST) at the 95% confidence level, while 14% were identified as candidates for local adaptation with the Bayesian method. There was some overlap between the two gene sets. More than half of the candidate genes for local adaptation were specific to the warmest population, about 20% to the most arid population, and 15% to the coldest and most humid higher altitude population. These adaptive trends were consistent with the genes' putative functions and the divergence in quantitative traits noted among the populations. The results suggest that an approach separating the locus and population effects is useful to identify genes potentially under selection. These candidates are worth exploring in more details at the physiological and ecological levels.

Consistent geographic structure among multiple nuclear sequences and cpDNA polymorphisms of Cardamine nipponica Franch. et Savat. (Brassicaceae)

Molecular phylogeography has inferred the history of differentiation between regions and/or among populations following the Pleistocene climatic oscillations, mostly based on the genetic structure of organelle DNA. However, such genetic structure only reflects the history of a single gene, and studies based on single-copy genes of nuclear DNA (nDNA) are required for phylogeography, although their efficiency remains unclear. To examine the utility of nDNA loci, the genetic structures of three genes from Cardamine nipponica, which is closely related to the model species Arabidopsis thaliana, were elucidated: the nDNA genes DET1, PHYA, PHYE, as well as chloroplast DNA (cpDNA). In 279 individuals collected from throughout the range of the species, strong genetic differentiation between northern and central Japan was found for all loci. This result suggested that populations in central Japan experienced a different history from those in northern Japan during the Pleistocene climatic oscillations. In addition, the evidence of refugia at the edges of the distribution, where the genetic structure was less influenced by colonization following range expansion, was shown for several loci. The specific genetic structure within the southernmost populations of northern Japan suggested that this region was also isolated during range expansion. Hence, the consistent history among loci and a more detailed history from several loci indicated that cpDNA can represent the history of vicariance and demonstrated the efficiency of single-copy nuclear genes in phylogeography.

Ecological and evolutionary consequences of biotic homogenization

Biotic homogenization, the gradual replacement of native biotas by locally expanding non-natives, is a global process that diminishes floral and faunal distinctions among regions. Although patterns of homogenization have been well studied, their specific ecological and evolutionary consequences remain unexplored. We argue that our current perspective on biotic homogenization should be expanded beyond a simple recognition of species diversity loss, towards a synthesis of higher order effects. Here, we explore three distinct forms of homogenization (genetic, taxonomic and functional), and discuss their immediate and future impacts on ecological and evolutionary processes. Our goal is to initiate future research that investigates the broader conservation implications of homogenization and to promote a proactive style of adaptive management that engages the human component of the anthropogenic blender that is currently mixing the biota on Earth.

Recombination and genetic differentiation among natural populations of the ectomycorrhizal mushroom Tricholoma matsutake from southwestern China

Effective conservation and utilization strategies for natural biological resources require a clear understanding of the natural populations of the target organisms. Tricholoma matsutake is an ectomycorrhizal mushroom that forms symbiotic associations with plants and plays an important ecological role in natural forest ecosystems in many parts of the world. It is also an economically very important gourmet mushroom. Because no artificial cultivation is available, natural populations of this species are under increasing threats, primarily from habitat disturbance and destruction. Despite its economical and ecological importance, little is known about its genetics and population biology. Here, using 14 polymerase chain reaction-restriction fragment length polymorphism markers, we analysed 154 strains from 17 geographical locations in southwestern China, a region where over 25% of the global T. matsutake harvest comes from. Our results revealed abundant genetic variation within individual populations. The analyses of gene and genotype frequencies within populations indicated that most loci did not deviate from Hardy-Weinberg equilibrium in most populations and that alleles among loci were in linkage equilibrium in the majority of the local populations. These results are consistent with the hypothesis that sexual reproduction and recombination play an important role in natural populations of this species. Our analyses indicated low but significant genetic differentiation among the geographical populations, with a significant positive correlation between genetic distance and geographical distance. We discuss the implications of our results to the ecology and resource management of this species.

Haplotype reconstruction for scnp DNA: a consensus vote approach with extensive sequence data from populations of the migratory locust (Locusta migratoria)

Single copy nuclear polymorphic (scnp) DNA is potentially a powerful molecular marker for evolutionary studies of populations. However, a practical obstacle to its employment is the general problem of haplotype determination due to the common occurrence of heterozygosity in diploid organisms. We explore here a 'consensus vote' (CV) approach to this question, combining statistical haplotype reconstruction and experimental verification using as an example an indel-free scnp DNA marker from the flanking region of a microsatellite locus of the migratory locust. The raw data comprise 251-bp sequences from 526 locust individuals (1052 chromosomes), with 71 (28.3%) polymorphic nucleotide sites (including seven triallelic sites) and 141 distinct genotypes (with frequencies ranging from 0.2 to 25.5%). Six representative statistical haplotype reconstruction algorithms are employed in our CV approach, including one parsimony method, two expectation-maximization (EM) methods and three Bayesian methods. The phases of 116 ambiguous individuals inferred by this approach are verified by molecular cloning experiments. We demonstrate the effectiveness of the CV approach compared to inferences based on individual statistical algorithms. First, it has the unique power to partition the inferrals into a reliable group and an uncertain group, thereby allowing the identification of the inferrals with greater uncertainty (12.7% of the total sample in this case). This considerably reduces subsequent efforts of experimental verification. Second, this approach is capable of handling genotype data pooled from many geographical populations, thus tolerating heterogeneity of genetic diversity among populations. Third, the performance of the CV approach is not influenced by the number of heterozygous sites in the ambiguous genotypes. Therefore, the CV approach is potentially a reliable strategy for effective haplotype determination of nuclear DNA markers. Our results also show that rare variations and rare inferrals tend to be more vulnerable to inference error, and hence deserve extra surveillance.

Development and application of a phylogenomic toolkit: resolving the evolutionary history of Madagascar's lemurs

Lemurs and the other strepsirrhine primates are of great interest to the primate genomics community due to their phylogenetic placement as the sister lineage to all other primates. Previous attempts to resolve the phylogeny of lemurs employed limited mitochondrial or small nuclear data sets, with many relationships poorly supported or entirely unresolved. We used genomic resources to develop 11 novel markers from nine chromosomes, representing approximately 9 kb of nuclear sequence data. In combination with previously published nuclear and mitochondrial loci, this yields a data set of more than 16 kb and adds approximately 275 kb of DNA sequence to current databases. Our phylogenetic analyses confirm hypotheses of lemuriform monophyly and provide robust resolution of the phylogenetic relationships among the five lemuriform families. We verify that the genus Daubentonia is the sister lineage to all other lemurs. The Cheirogaleidae and Lepilemuridae are sister taxa and together form the sister lineage to the Indriidae; this clade is the sister lineage to the Lemuridae. Divergence time estimates indicate that lemurs are an ancient group, with their initial diversification occurring around the Cretaceous-Tertiary boundary. Given the power of this data set to resolve branches in a notoriously problematic area of primate phylogeny, we anticipate that our phylogenomic toolkit will be of value to other studies of primate phylogeny and diversification. Moreover, the methods applied will be broadly applicable to other taxonomic groups where phylogenetic relationships have been notoriously difficult to resolve.

Insertion-deletion polymorphisms (indels) as genetic markers in natural populations

BACKGROUND

We introduce the use of short insertion-deletion polymorphisms (indels) for genetic analysis of natural populations.

RESULTS

Sequence reads from light shot-gun sequencing efforts of different dog breeds were aligned to the dog genome reference sequence and gaps corresponding to indels were identified. One hundred candidate markers (4-bp indels) were selected and genotyped in unrelated dogs (n = 7) and wolves (n = 18). Eighty-one and 76 out of 94 could be validated as polymorphic loci in the respective sample. Mean indel heterozygosity in a diverse set of wolves was 19%, and 74% of the loci had a minor allele frequency of >10%. Indels found to be polymorphic in wolves were subsequently genotyped in a highly bottlenecked Scandinavian wolf population. Fifty-one loci turned out to be polymorphic, showing their utility even in a population with low genetic diversity. In this population, individual heterozygosity measured at indel and microsatellite loci were highly correlated.

CONCLUSION

With an increasing amount of sequence information gathered from non-model organisms, we suggest that indels will come to form an important source of genetic markers, easy and cheap to genotype, for studies of natural populations.

A smörgåsbord of markers for avian ecology and evolution

The polymerase chain reaction has been a boon to the study of molecular ecology and population genetics of birds. But the nagging truth is that for many bird species, the number of polymerase chain reaction (PCR) primer pairs that one can pick off the shelf and expect to amplify their target loci with ease is frustratingly small. Now, studying DNA sequence variation in natural populations of birds just got a whole lot easier. This issue of Molecular Ecology reports a large-scale bioinformatics search for exonic sequences conserved between the chicken and zebra finch genomes and flanking polymorphic introns that has generated a staggering 242 PCR primer pairs that readily amplify their single-copy target loci in five avian species spanning approximately 100 million years of avian evolution (Backström et al. 2008). As proof of principle, these primers have also been used to survey the genomic landscape in over 110 kb of intronic sequence in the collared flycatcher, a model species in ecology and evolution. These resources pave the way for easy multilocus study of evolving populations and lineages of birds, and bring the goal of quickly turning nonmodel species in to ecological genomic models tantalizingly close.

Inferring patterns of migration

The historical movement of organisms, whether recent or in the distant past, forms a central aspect of evolutionary studies. Inferring patterns of migration can be difficult and requires reliance on a large suite of bioinformatic tools. As it is primarily the movement of groups of related individuals or populations that are of interest, population genetic and phylogeographic methods form the core of tools used to decipher migration patterns. Following a description of these tools, we discuss the most critical--and potentially most difficult--aspect of these studies: the inference process used. Designing a study, determining which data to collect, how to analyze the data, and how to coordinate these results into a coherent narrative are all a part of this inference process. Furthermore, using different types of data (e.g., genotypic and DNA sequence) from different types of sources (direct, or from the organisms of interest; and indirect, from symbiotic organisms) produces a powerful suite of techniques that are used to infer patterns of migration.

Population genetic studies of Alouatta caraya (Alouattinae, Primates): inferences on geographic distribution and ecology

Cytochrome b DNA sequence data (ca. 1,140 bp) of 44 Alouatta caraya, including 42 specimens from three localities of Brazil and two from Bolivia, were used for phylogenetic reconstructions and population studies. Seventeen haplotypes were identified, eight of which were present in more than one individual. Seven of these eight haplotypes were shared by individuals from a same locality and one by individuals from two localities. We found 26 variable sites along the entire gene, consisting of 18 transitions and eight transversions; most replacements occurring at the third codon position (65.39%) in contrast to first and second positions (26.92 and 7.69%, respectively). In the sample collected at Chapada dos Guimarães (Brazil), nucleotide and haplotype diversity estimates were pi=0.002325 and h=0.8772, respectively. Maximum parsimony analysis grouped all haplotypes in two clades, separating Bolivian haplotypes from Brazilian haplotypes, the grouping of which did not show a straightforward correspondence with geographic distribution. Median-joining and TCS network pointed to haplotypes 11 or 12 as the most likely ancestral ones. Mismatch distribution and the goodness-of-fit test (SSD estimate=0.0027; P=0.6999) indicated that the population from Chapada dos Guimarães experienced a demographic expansion, in agreement with the median-joining star-like pattern, although this finding could not be confirmed by Fu's F(s) test.

Single nucleotide polymorphisms in the ectomycorrhizal mushroom Tricholoma matsutake

Single nucleotide polymorphisms (SNPs) are becoming increasingly popular markers for studying a variety of biological phenomena. This paper describes the development and analysis of a set of SNP markers for the basidiomycete fungus Tricholoma matsutake. T. matsutake is a gourmet mycorrhizal mushroom primarily associated with pine forests. However, little is known about its genetics and genomic variation, including SNP variation. To identify and analyse SNPs in T. matsutake, a genomic library was constructed and >72,000 nt were analysed from >200 random clones. Primers from 20 sequenced fragments were then designed to amplify and sequence >10,000 bp sequences from the original strain, from which the genomic library was constructed, as well as another strain from >350 km away; both strains were from south-western China. These two strains had similar intra-strain SNP frequencies (1.104 and 1.278 % per nucleotide, respectively). The combined analysis revealed that 14 of the 20 examined fragments contained SNPs, ranging from two to 47 per fragment, and yielding a total of 178 SNPs out of the 10 428 sequenced nucleotides (an SNP frequency of 1.707 %). Among the 178 SNPs, one site had three alternative nucleotides, while the remaining 177 had two each, with 148 transitions and 29 transversions, resulting in a combined transition to transversion ratio of 5.1. In addition, the haplotype phases of all SNPs within individual fragments for both strains were determined. Phylogenetic analyses of these haplotypes revealed three kinds of haplotype relationship, including haplotype sharing both within and between strains. Furthermore, a subset of the SNPs detectable by restriction enzyme digests was screened for its distribution among 31 additional wild strains from five distinct locations in south-western China. The implications of these SNPs and haplotypes for our understanding of the genetics, population history, ecology and evolution of this important mushroom species are discussed.

Differential gene flow of mitochondrial and nuclear DNA markers among chromosomal races of Australian morabine grasshoppers (Vandiemenella, viatica species group)

Recent theoretical developments have led to a renewed interest in the potential role of chromosomal rearrangements in speciation. Australian morabine grasshoppers (genus Vandiemenella, viatica species group) provide an excellent study system to test this potential role of chromosomal rearrangements because they show extensive chromosomal variation and formed the basis of a classic chromosomal speciation model. There are three chromosomal races, viatica19, viatica17, and P24(XY), on Kangaroo Island, South Australia, forming five parapatric populations with four putative contact zones among them. We investigate the extent to which chromosomal variation among these populations may be associated with barriers to gene flow. Population genetic and phylogeographical analyses using 15 variable allozyme loci and the elongation factor-1alpha (EF-1alpha) gene indicate that the three races represent genetically distinct taxa. In contrast, analyses of the mitochondrial cytochrome c oxidase subunit I (COI) gene show the presence of three distinctive and geographically localized groups that do not correspond with the distribution of the chromosomal races. These discordant population genetic patterns are likely to result from introgressive hybridization between the chromosomal races and range expansions/contractions. Overall, these results suggest that reduction of nuclear gene flow may be associated with chromosomal variation, or underlying genetic variation linked with chromosomal variation, whereas mitochondrial gene flow appears to be independent of this variation in these morabine grasshoppers. The identification of an intact contact zone between P24(XY) and viatica17 offers considerable potential for further investigation of molecular mechanisms that maintain distinct nuclear genomes among the chromosomal races.

Transcriptional regulation of NADP-dependent malate dehydrogenase: comparative genetics and identification of DNA-binding proteins

The transcriptional regulation of NADP-malate dehydrogenase (NADP-MDH) was analyzed in Arabidopsis ecotypes and other Brassicaceae. The amount of transcript increased twofold after transfer into low temperature (12 degrees C) or high light (750 microE) in all species. Analysis of the genomic DNA reveals that the NADP-MDH gene (At5g58330 in A. thaliana) in Brassicaceae is located between two other genes (At5g58320 and At5g58340 in Arabidopsis), both encoded on the opposite DNA strand. No promoter elements were identified in 5' direction of the NADP-MDH gene, and the expression of NADP-MDH was not affected in knock-out plants carrying a DNA insert in the 5' region. A yeast-one hybrid approach yielded only three DNA-binding proteins for the 500-bp fragment located upstream of the ATG sequence, but 34 proteins for its coding region. However, in Chlamydomonas and in some Poaceae, which do not possess any genes within the 1200 bp upstream region, typical promoter elements were identified. Alignments of genomic DNA reveal that, in contrast to Poaceae, the introns are highly conserved within Brassicaceae. We conclude that in Brassicaceae the majority of regulatory elements are located within the coding region. The NADP-MDH gene of both families evolved from a common precursor, similar to the gene in Chlamydomonas. Changes in the selection pressure allowed the insertion of At5g58340 into the promoter region of a common ancestor. When the demand for transcriptional regulation increased, At5g58340 disappeared in Poaceae, and a promoter developed in the 5' region. In contrast, Brassicaceae maintained At5g58340 and shifted all regulatory elements into the coding region of NADP-MDH.

Multi-locus DNA sequence data reveal a history of deep cryptic vicariance and habitat-driven convergence in the desert night lizardXantusia vigilisspecies complex (Squamata: Xantusiidae)

The lizard genus Xantusia of southwestern North America has received recent attention in relation to delimiting species. Using more than 500 lizards from 156 localities, we further test hypothesized species boundaries and clarify phylogeographical patterns, particularly in regions of potential secondary contact. We sequenced the entire mitochondrial cytochrome b gene for every lizard in the study, plus a second mitochondrial DNA (mtDNA) region and two nuclear introns for subsets of the total sample. Phylogenetic analyses of the mtDNA recover a well-resolved, novel hypothesis for species in the Xantusia vigilis complex. The nuclear DNA (nDNA) data provide independent support for the recognition of X. arizonae, X. bezyi and X. wigginsi. Differences between the respective mtDNA and nDNA topologies result from either the effects of lineage sorting or ancient introgression. Nuclear data confirm the inference that some populations of X. vigilis in northwestern Arizona converged on rock-crevice-dwelling morphology and are not X. arizonae with an introgressed X. vigilis mtDNA genome. The historical independence of ancient cryptic lineages of Xantusia in southern California is also corroborated, though limited introgression is detected. Our proposed biogeographical scenario indicates that diversification of this group was driven by vicariance beginning in the late Miocene. Additionally, Pleistocene climatical changes influenced Xantusia distribution, and the now inhospitable Colorado Desert previously supported night lizard presence. The current taxonomy of the group likely underestimates species diversity within the group, and our results collectively show that while convergence on the rock-crevice-dwelling morphology is one hallmark of Xantusia evolution, morphological stasis is paradoxically another.

Leishmania and the leishmaniases: a parasite genetic update and advances in taxonomy, epidemiology and pathogenicity in humans

Leishmaniases remain a major public health problem today despite the vast amount of research conducted on Leishmania pathogens. The biological model is genetically and ecologically complex. This paper explores the advances in Leishmania genetics and reviews population structure, taxonomy, epidemiology and pathogenicity. Current knowledge of Leishmania genetics is placed in the context of natural populations. Various studies have described a clonal structure for Leishmania but recombination, pseudo-recombination and other genetic processes have also been reported. The impact of these different models on epidemiology and the medical aspects of leishmaniases is considered from an evolutionary point of view. The role of these parasites in the expression of pathogenicity in humans is also explored. It is important to ascertain whether genetic variability of the parasites is related to the different clinical expressions of leishmaniasis. The review aims to put current knowledge of Leishmania and the leishmaniases in perspective and to underline priority questions which 'leishmaniacs' must answer in various domains: epidemiology, population genetics, taxonomy and pathogenicity. It concludes by presenting a number of feasible ways of responding to these questions.

Molecular evolutionary genomics of birds

Insight into the molecular evolution of birds has been offered by the steady accumulation of avian DNA sequence data, recently culminating in the first draft sequence of an avian genome, that of chicken. By studying avian molecular evolution we can learn about adaptations and phenotypic evolution in birds, and also gain an understanding of the similarities and differences between mammalian and avian genomes. In both these lineages, there is pronounced isochore structure with highly variable GC content. However, while mammalian isochores are decaying, they are maintained in the chicken lineage, which is consistent with a biased gene conversion model where the high and variable recombination rate of birds reinforces heterogeneity in GC. In Galliformes, GC is positively correlated with the rate of nucleotide substitution; the mean neutral mutation rate is 0.12-0.15% at each site per million years but this estimate comes with significant local variation in the rate of mutation. Comparative genomics reveals lower d(N)/d(S) ratios on micro- compared to macrochromosomes, possibly due to population genetic effects or a non-random distribution of genes with respect to chromosome size. A non-random genomic distribution is shown by genes with sex-biased expression, with male-biased genes over-represented and female-biased genes under-represented on the Z chromosome. A strong effect of selection is evident on the non-recombining W chromosome with high d(N)/d(S) ratios and limited polymorphism. Nucleotide diversity in chicken is estimated at 4-5 x 10(-3) which might be seen as surprisingly high given presumed bottlenecks during domestication, but is lower than that recently observed in several natural populations of other species. Several important aspects of the molecular evolutionary process of birds remain to be understood and it can be anticipated that the upcoming genome sequence of a second bird species, the zebra finch, as well as the integration of data on gene expression, shall further advance our knowledge of avian evolution.

Molecular evolution and genome divergence at RPB2 gene of the St and H genome in Elymus species

Molecular evolution of the second largest subunit of low copy nuclear RNA polymerase II (RPB2) in allotetrploid StH genomic species of Elymus is characterized here. Our study first reported a 39-bp MITE stowaway element insertion in the genic region of RPB2 gene for all tetraploid Elymus St genome and diploid Pseudoroegneria spicata and P. stipifolia St genome. The sequences on 3'-end are highly conserved, with AGTA in all sequences but H10339 (E. fibrosis), in which the AGTA was replaced with AGAA. All 12 Stowaway-containing sequences encompassed a 9 bp conserved TIRs (GAGGGAGTA). Interestingly, the 5'-end sequence of GGTA which was changed to AGTA or deleted resulted in Stowaway excision in the H genome of Elymus sepcies, in which Stowaway excision did not leave footprint. Another two large insertions in all St genome sequences are also transposable-like elements detected in the genic region of RPB2 gene. Our results indicated that these three transposable element indels have occurred prior to polyploidization, and shaped the homoeologous RPB2 loci in St and H genome of Eymus species. Nucleotide diversity analysis suggested that the RPB2 sequence may evolve faster in the polyploid species than in the diploids. Higher level of polymorphism and genome-specific amplicons generated by this gene indicated that RPB2 is an excellent tool for investigating the phylogeny and evolutionary dynamics of speciation, and the mode of polyploidy formation in Elymus species.

Historical demography of brown trout (Salmo trutta) in the Adriatic drainage including the putative S. letnica endemic to Lake Ohrid

We explore the historical demography of the Adriatic lineage of brown trout and more explicitly the colonization and phylogenetic placement of Ohrid trout, based on variation at 12 microsatellite loci and the mtDNA control region. All Adriatic basin haplotypes reside in derived positions in a network that represents the entire lineage. The central presumably most ancestral haplotype in this network is restricted to the Iberian Peninsula, where it is very common, supporting a Western Mediterranean origin for the lineage. The expansion statistic R2, Bayesian based estimates of demographic parameters, and star-like genealogies support expansions on several geographic scales, whereas application of pairwise mismatch analysis was somewhat ambiguous. The estimated time since expansion (155,000 years ago) for the Adriatic lineage was supported by a narrow confidence interval compared to previous studies. Based on microsatellite and mtDNA sequence variation, the endemic Ohrid trout represents a monophyletic lineage isolated from other Adriatic basin populations, but nonetheless most likely evolving from within the Adriatic lineage of brown trout. Our results do not support the existence of population structuring within Lake Ohrid, even though samples included two putative intra-lacustrine forms. In the interests of protecting the unique biodiversity of this ancient ecosystem, we recommend retaining the taxonomic epithet Salmo letnica for the endemic Ohrid trout.

AFLP and single-strand conformation polymorphism studies of recombination in the entomopathogenic fungus Nomuraea rileyi

In most putative asexual fungi analysed through population genetic studies, recombination has been detected. However, the mechanism by which it is achieved is still not known. A parasexual cycle is known to occur in asexual fungi but there is no evidence, as yet, of its prevalence in natural populations. This study was undertaken to investigate the possibility of a parasexual cycle mediating recombination in the mitosporic fungus Nomuraea rileyi. The genotypic diversity in isolates sampled from an epizootic population from South India was studied through AFLP. The AFLP data were subjected to analysis of molecular variance (AMOVA) and cluster analysis. Great genetic variation was observed in the population including the isolates from a single insect. To assess the occurrence of recombination in the population, single-strand conformation polymorphism (SSCP) of partial regions of two mitochondrial (mt) genes (rRNA genes of LSU and SSU) and a nuclear gene (beta tubulin) was performed. The SSCP data were analysed using MP, the tree length permutation test, and multilocus analysis. Recombination was inferred from the SSCP analysis. The occurrence of isolates with diverse genotypes in a single insect; the fact that fungi multiply as hyphal bodies (cell wall-less) in the insect haemolymph; and the inference of recombination in mitochondrial genes (suggesting cytomixis), all indicate that recombination is accomplished by fusion of hyphal bodies of different isolates infecting the insect.

DNA barcoding: how it complements taxonomy, molecular phylogenetics and population genetics

DNA barcoding aims to provide an efficient method for species-level identifications and, as such, will contribute powerfully to taxonomic and biodiversity research. As the number of DNA barcode sequences accumulates, however, these data will also provide a unique 'horizontal' genomics perspective with broad implications. For example, here we compare the goals and methods of DNA barcoding with those of molecular phylogenetics and population genetics, and suggest that DNA barcoding can complement current research in these areas by providing background information that will be helpful in the selection of taxa for further analyses.

Reliable allele detection using SNP-based PCR primers containing Locked Nucleic Acid: application in genetic mapping

BACKGROUND

The diploid, Solanum caripense, a wild relative of potato and tomato, possesses valuable resistance to potato late blight and we are interested in the genetic base of this resistance. Due to extremely low levels of genetic variation within the S. caripense genome it proved impossible to generate a dense genetic map and to assign individual Solanum chromosomes through the use of conventional chromosome-specific SSR, RFLP, AFLP, as well as gene- or locus-specific markers. The ease of detection of DNA polymorphisms depends on both frequency and form of sequence variation. The narrow genetic background of close relatives and inbreds complicates the detection of persisting, reduced polymorphism and is a challenge to the development of reliable molecular markers. Nonetheless, monomorphic DNA fragments representing not directly usable conventional markers can contain considerable variation at the level of single nucleotide polymorphisms (SNPs). This can be used for the design of allele-specific molecular markers. The reproducible detection of allele-specific markers based on SNPs has been a technical challenge.

RESULTS

We present a fast and cost-effective protocol for the detection of allele-specific SNPs by applying Sequence Polymorphism-Derived (SPD) markers. These markers proved highly efficient for fingerprinting of individuals possessing a homogeneous genetic background. SPD markers are obtained from within non-informative, conventional molecular marker fragments that are screened for SNPs to design allele-specific PCR primers. The method makes use of primers containing a single, 3'-terminal Locked Nucleic Acid (LNA) base. We demonstrate the applicability of the technique by successful genetic mapping of allele-specific SNP markers derived from monomorphic Conserved Ortholog Set II (COSII) markers mapped to Solanum chromosomes, in S. caripense. By using SPD markers it was possible for the first time to map the S. caripense alleles of 16 chromosome-specific COSII markers and to assign eight of the twelve linkage groups to consensus Solanum chromosomes.

CONCLUSION

The method based on individual allelic variants allows for a level-of-magnitude higher resolution of genetic variation than conventional marker techniques. We show that the majority of monomorphic molecular marker fragments from organisms with reduced heterozygosity levels still contain SNPs that are sufficient to trace individual alleles.

Automated discovery of single nucleotide polymorphism and simple sequence repeat molecular genetic markers

Molecular genetic markers represent one of the most powerful tools for the analysis of genomes. Molecular marker technology has developed rapidly over the last decade, and two forms of sequence-based markers, simple sequence repeats (SSRs), also known as microsatellites, and single nucleotide polymorphisms (SNPs), now predominate applications in modern genetic analysis. The availability of large sequence data sets permits mining for SSRs and SNPs, which may then be applied to genetic trait mapping and marker-assisted selection. Here, we describe Web-based automated methods for the discovery of these SSRs and SNPs from sequence data. SSRPrimer enables the real-time discovery of SSRs within submitted DNA sequences, with the concomitant design of PCR primers for SSR amplification. Alternatively, users may browse the SSR Taxonomy Tree to identify predetermined SSR amplification primers for any species represented within the GenBank database. SNPServer uses a redundancy-based approach to identify SNPs within DNA sequence data. Following submission of a sequence of interest, SNPServer uses BLAST to identify similar sequences, CAP3 to cluster and assemble these sequences, and then the SNP discovery software autoSNP to detect SNPs and insertion/deletion (indel) polymorphisms.

Types and rates of sequence evolution at the high-molecular-weight glutenin locus in hexaploid wheat and its ancestral genomes

The Glu-1 locus, encoding the high-molecular-weight glutenin protein subunits, controls bread-making quality in hexaploid wheat (Triticum aestivum) and represents a recently evolved region unique to Triticeae genomes. To understand the molecular evolution of this locus region, three orthologous Glu-1 regions from the three subgenomes of a single hexaploid wheat species were sequenced, totaling 729 kb of sequence. Comparing each Glu-1 region with its corresponding homologous region from the D genome of diploid wheat, Aegilops tauschii, and the A and B genomes of tetraploid wheat, Triticum turgidum, revealed that, in addition to the conservation of microsynteny in the genic regions, sequences in the intergenic regions, composed of blocks of nested retroelements, are also generally conserved, although a few nonshared retroelements that differentiate the homologous Glu-1 regions were detected in each pair of the A and D genomes. Analysis of the indel frequency and the rate of nucleotide substitution, which represent the most frequent types of sequence changes in the Glu-1 regions, demonstrated that the two A genomes are significantly more divergent than the two B genomes, further supporting the hypothesis that hexaploid wheat may have more than one tetraploid ancestor.

Development of a lab-made microarray for analyzing the genetic diversity of nitrogen fixing symbionts Sinorhizobium meliloti and Sinorhizobium medicae

Some bacterial species, like nitrogen-fixing Sinorhizobium that interact with Medicago plants, are prone to frequent horizontal gene transfers. Investigation of their genetic structure requires to study polymorphism patterns at many loci. Although DNA microarrays represent a method of choice for high throughput analysis of polymorphisms, this technology yet remains an expensive and heavy approach, thus depriving most of research groups from this powerful tool. In an attempt to overcome this limitation, we have developed a simple genotyping procedure by DNA microarrays, and have evaluated its ability to characterize a Sinorhizobium population. Thirty 18- to 24-mer oligonucleotide probes were designed to target the most frequent mutations in three polymorphic loci of Sinorhizobium meliloti and S. medicae. Probe hybridization efficiency was compared on two spotting surfaces: nylon membranes and epoxy-coated glass slides. Epoxy-coated glass slides revealed more sensitive than nylon membranes and allowed discrimination of single mismatches. Using this procedure, an uncharacterized population consisting of 33 S. meliloti/S. medicae isolates was successfully genotyped.

Multiple nuclear gene sequences identify phylogenetic species boundaries in the rapidly radiating clade of Mexican ambystomatid salamanders

Delimiting the boundaries of species involved in radiations is critical to understanding the tempo and mode of lineage formation. Single locus gene trees may or may not reflect the underlying pattern of population divergence and lineage formation, yet they constitute the vast majority of the empirical data in species radiations. In this study we make use of an expressed sequence tag (EST) database to perform nuclear (nDNA) and mitochondrial (mtDNA) genealogical tests of species boundaries in Ambystoma ordinarium, a member of an adaptive radiation of metamorphic and paedomorphic salamanders (the Ambystoma tigrinum complex) that have diversified across terrestrial and aquatic environments. Gene tree comparisons demonstrate extensive nonmonophyly in the mtDNA genealogy of A. ordinarium, while seven of eight independent nuclear loci resolve the species as monophyletic or nearly so, and diagnose it as a well-resolved genealogical species. A differential introgression hypothesis is supported by the observation that western A. ordinarium localities contain mtDNA haplotypes that are identical or minimally diverged from haplotypes sampled from a nearby paedomorphic species, Ambystoma dumerilii, while most nDNA trees place these species in distant phylogenetic positions. These results provide a strong example of how historical introgression can lead to radical differences between gene trees and species histories, even among currently allopatric species with divergent life history adaptations and morphologies. They also demonstrate how EST-based nuclear resources can be used to more fully resolve the phylogenetic history of species radiations.