DnaSP v5: A software for comprehensive analysis of DNA polymorphism data

Podeu consultar el programari a: http://hdl.handle.net/2445/53451DnaSP is a software package for a comprehensive analysis of DNA polymorphism data. Version 5 implements a number of new features and analytical methods allowing extensive DNA polymorphism analyses on large datasets. Among other features, the newly implemented methods allow for: (i) analyses on multiple data files; (ii) haplotype phasing; (iii) analyses on insertion/deletion polymorphism data; (iv) visualizing sliding window results integrated with available genome annotations in the UCSC browser

DnaSP 6 : DNA sequence polymorphism analysis of large data sets

We present version 6 of the DnaSP (DNA Sequence Polymorphism) software, a new version of the popular tool for performing exhaustive population genetic analyses on multiple sequence alignments. This major upgrade incorporates novel functionalities to analyse large datasets, such as those generated by high-throughput sequencing (HTS) technologies. Among other features, DnaSP 6 implements: i) modules for reading and analysing data from genomic partitioning methods, such as RADseq or hybrid enrichment approaches, ii) faster methods scalable for HTS data, and iii) summary statistics for the analysis of multi-locus population genetics data. Furthermore, DnaSP 6 includes novel modules to perform single- and multi-locus coalescent simulations under a wide range of demographic scenarios. The DnaSP 6 program, with extensive documentation, is freely available at http://www.ub.edu/dnasp

High genetic diversity at the extreme range edge: nucleotide variation at nuclear loci in Scots pine (Pinus sylvestris L.) in Scotland

Nucleotide polymorphism at 12 nuclear loci was studied in Scots pine populations across an environmental gradient in Scotland, to evaluate the impacts of demographic history and selection on genetic diversity. At eight loci, diversity patterns were compared between Scottish and continental European populations. At these loci, a similar level of diversity (θsil=~0.01) was found in Scottish vs mainland European populations, contrary to expectations for recent colonization, however, less rapid decay of linkage disequilibrium was observed in the former (ρ=0.0086±0.0009, ρ=0.0245±0.0022, respectively). Scottish populations also showed a deficit of rare nucleotide variants (multi-locus Tajima's D=0.316 vs D=−0.379) and differed significantly from mainland populations in allelic frequency and/or haplotype structure at several loci. Within Scotland, western populations showed slightly reduced nucleotide diversity (πtot=0.0068) compared with those from the south and east (0.0079 and 0.0083, respectively) and about three times higher recombination to diversity ratio (ρ/θ=0.71 vs 0.15 and 0.18, respectively). By comparison with results from coalescent simulations, the observed allelic frequency spectrum in the western populations was compatible with a relatively recent bottleneck (0.00175 × 4Ne generations) that reduced the population to about 2% of the present size. However, heterogeneity in the allelic frequency distribution among geographical regions in Scotland suggests that subsequent admixture of populations with different demographic histories may also have played a role

Maternal Expression Relaxes Constraint on Innovation of the Anterior Determinant, bicoid

The origin of evolutionary novelty is believed to involve both positive selection and relaxed developmental constraint. In flies, the redesign of anterior patterning during embryogenesis is a major developmental innovation and the rapidly evolving Hox gene, bicoid (bcd), plays a critical role. We report evidence for relaxation of selective constraint acting on bicoid as a result of its maternal pattern of gene expression. Evolutionary theory predicts 2-fold greater sequence diversity for maternal effect genes than for zygotically expressed genes, because natural selection is only half as effective acting on autosomal genes expressed in one sex as it is on genes expressed in both sexes. We sample an individual from ten populations of Drosophila melanogaster and nine populations of D. simulans for polymorphism in the tandem gene duplicates bcd, which is maternally expressed, and zerknüllt (zen), which is zygotically expressed. In both species, we find the ratio of bcd to zen nucleotide diversity to be two or more in the coding regions but one in the noncoding regions, providing the first quantitative support for the theoretical prediction of relaxed selective constraint on maternal-effect genes resulting from sex-limited expression. Our results suggest that the accelerated rate of evolution observed for bcd is owing, at least partly, to variation generated by relaxed selective constraint

VariScan Software

Linux / Mac OS X : The package includes executables for linux (variscan) and Mac OS X (variscan). For other Unix-based platforms you will have to compile it from the source files included in the VariScan package. Windows: The package includes (src directory), the source code, the project (variscan.dev) and makefile (variscan.win) files to be used, for instance, for the Dev-C++ (a free Integrated Development Environment for the C/C++ programming language)Podeu consultar l'article relacionat a: http://hdl.handle.net/2445/7384Podeu consultar la pàgina de desenvolupament del programari: http://www.ub.edu/softevol/VariScan is a software package for the analysis of DNA sequence polymorphisms at the whole genome scale. Among other features, the software: (1) can conduct many population genetic analyses; (2) incorporates a multiresolution wavelet transform-based method that allows capturing relevant information from DNA polymorphism data; and (3) it facilitates the visualization of the results in the most commonly used genome browsers

Genome-wide DNA polymorphism analyses using VariScan

BACKGROUND: DNA sequence polymorphisms analysis can provide valuable information on the evolutionary forces shaping nucleotide variation, and provides an insight into the functional significance of genomic regions. The recent ongoing genome projects will radically improve our capabilities to detect specific genomic regions shaped by natural selection. Current available methods and software, however, are unsatisfactory for such genome-wide analysis. RESULTS: We have developed methods for the analysis of DNA sequence polymorphisms at the genome-wide scale. These methods, which have been tested on a coalescent-simulated and actual data files from mouse and human, have been implemented in the VariScan software package version 2.0. Additionally, we have also incorporated a graphical-user interface. The main features of this software are: i) exhaustive population-genetic analyses including those based on the coalescent theory; ii) analysis adapted to the shallow data generated by the high-throughput genome projects; iii) use of genome annotations to conduct a comprehensive analyses separately for different functional regions; iv) identification of relevant genomic regions by the sliding-window and wavelet-multiresolution approaches; v) visualization of the results integrated with current genome annotations in commonly available genome browsers. CONCLUSION: VariScan is a powerful and flexible suite of software for the analysis of DNA polymorphisms. The current version implements new algorithms, methods, and capabilities, providing an important tool for an exhaustive exploratory analysis of genome-wide DNA polymorphism data

Lineage specific recombination rates and microevolution in Listeria monocytogenes

Background: The bacterium Listeria monocytogenes is a saprotroph as well as an opportunistic human foodborne pathogen, which has previously been shown to consist of at least two widespread lineages (termed lineages I and II) and an uncommon lineage (lineage III). While some L. monocytogenes strains show evidence for considerable diversification by homologous recombination, our understanding of the contribution of recombination to L. monocytogenes evolution is still limited. We therefore used STRUCTURE and ClonalFrame, two programs that model the effect of recombination, to make inferences about the population structure and different aspects of the recombination process in L. monocytogenes. Analyses were performed using sequences for seven loci (including the house-keeping genes gap, prs, purM and ribC, the stress response gene sigB, and the virulence genes actA and inlA) for 195 L. monocytogenes isolates. Results: Sequence analyses with ClonalFrame and the Sawyer's test showed that recombination is more prevalent in lineage II than lineage I and is most frequent in two house-keeping genes (ribC and purM) and the two virulence genes (actA and inlA). The relative occurrence of recombination versus point mutation is about six times higher in lineage II than in lineage I, which causes a higher genetic variability in lineage II. Unlike lineage I, lineage II represents a genetically heterogeneous population with a relatively high proportion (30% average) of genetic material imported from external sources. Phylograms, constructed with correcting for recombination, as well as Tajima's D data suggest that both lineages I and II have suffered a population bottleneck. Conclusion: Our study shows that evolutionary lineages within a single bacterial species can differ considerably in the relative contributions of recombination to genetic diversification. Accounting for recombination in phylogenetic studies is critical, and new evolutionary models that account for the possibility of changes in the rate of recombination would be required. While previous studies suggested that only L. monocytogenes lineage I has experienced a recent bottleneck, our analyses clearly show that lineage II experienced a bottleneck at about the same time, which was subsequently obscured by abundant homologous recombination after the lineage II bottleneck. While lineage I and lineage II should be considered separate species from an evolutionary viewpoint, maintaining single species name may be warranted since both lineages cause the same type of human disease

Balancing Selection at the Tomato RCR3 Guardee Gene Family Maintains Variation in Strength of Pathogen Defense

Coevolution between hosts and pathogens is thought to occur between interacting molecules of both species. This results in the maintenance of genetic diversity at pathogen antigens (or so-called effectors) and host resistance genes such as the major histocompatibility complex (MHC) in mammals or resistance (R) genes in plants. In plant-pathogen interactions, the current paradigm posits that a specific defense response is activated upon recognition of pathogen effectors via interaction with their corresponding R proteins. According to the''Guard-Hypothesis,'' R proteins (the ``guards'') can sense modification of target molecules in the host (the ``guardees'') by pathogen effectors and subsequently trigger the defense response. Multiple studies have reported high genetic diversity at R genes maintained by balancing selection. In contrast, little is known about the evolutionary mechanisms shaping the guardee, which may be subject to contrasting evolutionary forces. Here we show that the evolution of the guardee RCR3 is characterized by gene duplication, frequent gene conversion, and balancing selection in the wild tomato species Solanum peruvianum. Investigating the functional characteristics of 54 natural variants through in vitro and in planta assays, we detected differences in recognition of the pathogen effector through interaction with the guardee, as well as substantial variation in the strength of the defense response. This variation is maintained by balancing selection at each copy of the RCR3 gene. Our analyses pinpoint three amino acid polymorphisms with key functional consequences for the coevolution between the guardee (RCR3) and its guard (Cf-2). We conclude that, in addition to coevolution at the ``guardee-effector'' interface for pathogen recognition, natural selection acts on the ``guard-guardee'' interface. Guardee evolution may be governed by a counterbalance between improved activation in the presence and prevention of auto-immune responses in the absence of the corresponding pathogen

Population Genetic Structure and Species Delimitation of a Widespread, Neotropical Dwarf Gecko

Amazonia harbors the greatest biological diversity on Earth. One trend that spans Amazonian taxa is that most taxonomic groups either exhibit broad geographic ranges or small restricted ranges. This is likely because many traits that determine a species range size, such as dispersal ability or body size, are autocorrelated. As such, it is rare to find groups that exhibit both large and small ranges. Once identified, however, these groups provide a powerful system for isolating specific traits that influence species distributions. One group of terrestrial vertebrates, gecko lizards, tends to exhibit small geographic ranges. Despite one exception, this applies to the Neotropical dwarf geckos of the genus Gonatodes. This exception, Gonatodes humeralis, has a geographic distribution almost 1,000,000 km2 larger than the combined ranges of its 30 congeners. As the smallest member of its genus and a gecko lizard more generally, G. humeralis is an unlikely candidate to be a wide-ranged Amazonian taxon. To test whether or not G. humeralis is one or more species, we generated molecular genetic data using restriction-site associated sequencing (RADseq) and traditional Sanger methods for samples from across its range and conducted a phylogeographic study. We conclude that G. humeralis is, in fact, a single species across its contiguous range in South America. Thus, Gonatodes is a unique clade among Neotropical taxa, containing both wide-ranged and range-restricted taxa, which provides empiricists with a powerful model system to correlate complex species traits and distributions. Additionally, we provide evidence to support species-level divergence of the allopatric population from Trinidad and we resurrect the name Gonatodes ferrugineus from synonymy for this population