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.     

Insights into the evolutionary process from patterns of DNA sequence variability

Studies of DNA sequence diversity, particularly in Drosophila, reveal both the complexity of natural selection and the importance of the interaction of natural selection and variation in rates of recombination within genomes and between species in determining levels of sequence variability in different genes and different species. The impact of both adaptive and deleterious mutations are evident. Extension of these types of studies to other organisms has begun in earnest.     

Historical population size change of bowhead whales inferred from DNA sequence polymorphism data

Nucleotide sequence data from the mitochondrial control region were used from a phylogenetic context to investigate the long-term history of a population of bowhead whales (Balaena mysticetus). In addition, the coalescence time of these sequences was used to estimate the age of the inferred patterns of population size change. The results indicate that mitochondrial genetic polymorphism was not affected by a recent bottleneck that occurred near the turn of the 20th century, thereby preserving the signature of historical population size change in the mitochondrial genome. Further analysis showed that this population underwent an expansion initiated in the Middle to Late Pleistocene. As such, early Holocene changes in Arctic sea ice distribution appear to have had little influence on patterns of genetic variability in this population.     

eBURST: inferring patterns of evolutionary descent among clusters of related bacterial genotypes from multilocus sequence typing data

The introduction of multilocus sequence typing (MLST) for the precise characterization of isolates of bacterial pathogens has had a marked impact on both routine epidemiological surveillance and microbial population biology. In both fields, a key prerequisite for exploiting this resource is the ability to discern the relatedness and patterns of evolutionary descent among isolates with similar genotypes. Traditional clustering techniques, such as dendrograms, provide a very poor representation of recent evolutionary events, as they attempt to reconstruct relationships in the absence of a realistic model of the way in which bacterial clones emerge and diversify to form clonal complexes. An increasingly popular approach, called BURST, has been used as an alternative, but present implementations are unable to cope with very large data sets and offer crude graphical outputs. Here we present a new implementation of this algorithm, eBURST, which divides an MLST data set of any size into groups of related isolates and clonal complexes, predicts the founding (ancestral) genotype of each clonal complex, and computes the bootstrap support for the assignment. The most parsimonious patterns of descent of all isolates in each clonal complex from the predicted founder(s) are then displayed. The advantages of eBURST for exploring patterns of evolutionary descent are demonstrated with a number of examples, including the simple Spain(23F)-1 clonal complex of Streptococcus pneumoniae, "population snapshots" of the entire S. pneumoniae and Staphylococcus aureus MLST databases, and the more complicated clonal complexes observed for Campylobacter jejuni and Neisseria meningitidis.     

proseq: A software for preparation and evolutionary analysis of DNA sequence data sets

proseq is an integrated user‐friendly windows based program for convenient sequence editing and evolutionary analysis. It is designed to simplify preparation and analysis of DNA sequence data sets in population genetic, phylogenetic and molecular ecology studies. Sequence editor features include editing of chromatogram files, contig assembly, sequence alignment, translation and other utilities. Analysis features include calculation of genetic diversity, divergence, population subdivision and gene flow with permutation‐based significance testing and various tests of neutrality. A tool for coalescent simulations implements models with intragenic recombination, population subdivision and population growth.     

Neandertal evolutionary genetics: mitochondrial DNA data from the iberian peninsula

Mitochondrial DNA (mtDNA) was retrieved for the first time from a Neandertal from the Iberian Peninsula, excavated from the El Sidrón Cave (Asturias, North of Spain), and dated to ca. 43,000 years ago. The sequence suggests that Iberian Neandertals were not genetically distinct from those of other regions. An estimate of effective population size indicates that the genetic history of the Neandertals was not shaped by an extreme population bottleneck associated with the glacial maximum of 130,000 years ago. A high level of polymorphism at sequence position 16258 reflects deeply rooted mtDNA lineages, with the time to the most recent common ancestor at ca. 250,000 years ago. This coincides with the full emergence of the "classical" Neandertal morphology and fits chronologically with a proposed speciation event of Homo neanderthalensis.     

A practical method for calculating evolutionary trees from sequence data

In a previous paper (Klotz et a1., 1979) we described a method for determining evolutionary trees from sequence data when rates of evolution of the sequences might differ greatly. It was shown theoretically that the method always gave the correct topology and root when the exact number of mutation differences between sequences and from their common ancestor was known. However, the method is impractical to use in most situations because it requires some knowledge of the ancestor. In this present paper we describe another method, related to the previous one, in which a present-day sequence can serve temporarily as an ancestor for purposes of determining the evolutionary tree regardless of the rates of evolution of the sequences involved. This new method can be carried out with high precision without the aid of a computer, and it does not increase in difficulty rapidly as the number of sequences involved in the study increases, unlike other methods.     

Detecting bottlenecks and selective sweeps from DNA sequence polymorphism

A coalescence-based maximum-likelihood method is presented that aims to (i) detect diversity-reducing events in the recent history of a population and (ii) distinguish between demographic (e.g., bottlenecks) and selective causes (selective sweep) of a recent reduction of genetic variability. The former goal is achieved by taking account of the distortion in the shape of gene genealogies generated by diversity-reducing events: gene trees tend to be more star-like than under the standard coalescent. The latter issue is addressed by comparing patterns between loci: demographic events apply to the whole genome whereas selective events affect distinct regions of the genome to a varying extent. The maximum-likelihood approach allows one to estimate the time and strength of diversity-reducing events and to choose among competing hypotheses. An application to sequence data from an African population of Drosophila melanogaster shows that the bottleneck hypothesis is unlikely and that one or several selective sweeps probably occurred in the recent history of this population.     

Inferring the number of evolutionary events from DNA coding sequence differences

The estimation of the amount of evolutionary divergence that has taken place between two DNA coding sequences depends strongly on the degree of constraint on amino acid replacements. If amino acid replacements are relatively unconstrained, the individual nucleotide is the appropriate unit of analysis and the method of Tajima and Nei can be used. If amino acid replacements are constrained, however, this method is shown to be inapplicable. For sequences with strong amino acid constraints, a method is outlined analogous to the Tajima and Nei method using codons as the unit of analysis. Only synonymous substitutions are used. Codon usage data can be employed to estimate the necessary parameters of the calculation, or a priori models of substitution may be employed. Sequences with significant but intermediate constraints on amino acid replacements are, in principle, unanalyzable.      

CopywriteR: DNA copy number detection from off-target sequence data

Current methods for detection of copy number variants (CNV) and aberrations (CNA) from targeted sequencing data are based on the depth of coverage of captured exons. Accurate CNA determination is complicated by uneven genomic distribution and non-uniform capture efficiency of targeted exons. Here we present CopywriteR, which eludes these problems by exploiting ‘off-target’ sequence reads. CopywriteR allows for extracting uniformly distributed copy number information, can be used without reference, and can be applied to sequencing data obtained from various techniques including chromatin immunoprecipitation and target enrichment on small gene panels. CopywriteR outperforms existing methods and constitutes a widely applicable alternative to available tools.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0617-1) contains supplementary material, which is available to authorized users.     

Analysis of sequence variation in Gnathostoma spinigerum mitochondrial DNA by single-strand conformation polymorphism analysis and DNA sequence

Morphological variations were observed in the advance third stage larvae of Gnathostoma spinigerum collected from swamp eel (Fluta alba), the second intermediate host. Larvae with typical and three atypical types were chosen for partial cytochrome c oxidase subunit I (COI) gene sequence analysis. A 450 bp polymerase chain reaction product of the COI gene was amplified from mitochondrial DNA. The variations were analyzed by single-strand conformation polymorphism and DNA sequencing. The nucleotide variations of the COI gene in the four types of larvae indicated the presence of an intra-specific variation of mitochondrial DNA in the G. spinigerum population.     

MultiSeq: unifying sequence and structure data for evolutionary analysis

Background  

Since the publication of the first draft of the human genome in 2000, bioinformatic data have been accumulating at an overwhelming pace. Currently, more than 3 million sequences and 35 thousand structures of proteins and nucleic acids are available in public databases. Finding correlations in and between these data to answer critical research questions is extremely challenging. This problem needs to be approached from several directions: information science to organize and search the data; information visualization to assist in recognizing correlations; mathematics to formulate statistical inferences; and biology to analyze chemical and physical properties in terms of sequence and structure changes.     

Detecting the form of selection from DNA sequence data

Clues to our evolutionary history lie hidden within DNA sequence data. One of the great challenges facing population geneticists is to identify and accurately interpret these clues. This task is made especially difficult by the fact that many different evolutionary processes can lead to similar observations. For example, low levels of polymorphism within a region can be explained by a low local mutation rate, by selection having eliminated deleterious mutations, or by the recent spread to fixation of a beneficial allele. Theoretical advances improve our ability to distinguish signals left by different evolutionary processes. In particular, a new test might better detect the footprint of selection having favored the spread of a beneficial allele.     

Accuracy of phylogenetic trees estimated from DNA sequence data

The relative merits of four different tree-making methods in obtaining the correct topology were studied by using computer simulation. The methods studied were the unweighted pair-group method with arithmetic mean (UPGMA), Fitch and Margoliash's (FM) method, thd distance Wagner (DW) method, and Tateno et al.'s modified Farris (MF) method. An ancestral DNA sequence was assumed to evolve into eight sequences following a given model tree. Both constant and varying rates of nucleotide substitution were considered. Once the DNA sequences for the eight extant species were obtained, phylogenetic trees were constructed by using corrected (d) and uncorrected (p) nucleotide substitutions per site. The topologies of the trees obtained were then compared with that of the model tree. The results obtained can be summarized as follows: (1) The probability of obtaining the correct rooted or unrooted tree is low unless a large number of nucleotide differences exists between different sequences. (2) When the number of nucleotide substitutions per sequence is small or moderately large, the FM, DW, and MF methods show a better performance than UPGMA in recovering the correct topology. The former group of methods is particularly good for obtaining the correct unrooted tree. (3) When the number of substitutions per sequence is large, UPGMA is at least as good as the other methods, particularly for obtaining the correct rooted tree. (4) When the rate of nucleotide substitution varies with evolutionary lineage, the FM, DW, and MF methods show a better performance in obtaining the correct topology than UPGMA, except when a rooted tree is to be produced from data with a large number of nucleotide substitutions per sequence.(ABSTRACT TRUNCATED AT 250 WORDS)      

Evolutionary and biogeographic patterns of the Badidae (Teleostei: Perciformes) inferred from mitochondrial and nuclear DNA sequence data

We reconstructed phylogenetic relationships of the family Badidae using both mitochondrial and nuclear nucleotide sequence data to address badid systematics and to evaluate the role of vicariant speciation on their evolution and current distribution. Phy-logenetic hypotheses were derived from complete cytochrome b (1,140 base pairs) sequences of 33 individuals representing 13 badid species, and using three species of Nandidae as outgroups. Additionally, we sequenced the nuclear RAG1 (1,473 base pairs) and Tmo-4C4 (511 base pairs) genes from each of the badid species and one representative of the outgroup. Our molecular data provide the first phylogenetic hypothesis of badid intrarelationships. Analysis of the mitochondrial and nuclear nucleotide sequence data sets resulted in well-supported trees, indicating a basal split between the genera Dario and Badis, and further supporting the division of the genus Badis into five species groups as suggested by a previous taxonomic revision of the Badidae. Within the genus Badis, mitochondrial and nuclear phylogenies differed in the relative position of B. kyar. We also used our molecular phylogeny to test a vicariant speciation hypothesis derived from geological evidence of large-scale changes in drainage patterns in the Miocene affecting the Irrawaddy- and Tsangpo-Brahmaputra drainages, in the southeastern Himalaya. Within both genera, Badis and Dario, we observed a divergence into Irrawaddy- and Tsangpo-Brahmaputra clades. Using a cytb substitution rate of 8.2 x 10(-9) (substitutions x base pair(-1) x year(-1), we tentatively date this vicariant event at the Oligocene-Miocene boundary (19-24Myr). It is concordant with a hypothesized paleo connection of the Tsangpo river with the Irrawaddy drainage that was most likely interrupted during Miocene orogenic events through tectonic uplifts in eastern Tibet. Our data, therefore, indicate a substantial role of vicariant-based speciation shaping the current distribution patterns of badids.     

Complex evolutionary history of the Aeromonas veronii group revealed by host interaction and DNA sequence data

Aeromonas veronii biovar sobria, Aeromonas veronii biovar veronii, and Aeromonas allosaccharophila are a closely related group of organisms, the Aeromonas veronii Group, that inhabit a wide range of host animals as a symbiont or pathogen. In this study, the ability of various strains to colonize the medicinal leech as a model for beneficial symbiosis and to kill wax worm larvae as a model for virulence was determined. Isolates cultured from the leech out-competed other strains in the leech model, while most strains were virulent in the wax worms. Three housekeeping genes, recA, dnaJ and gyrB, the gene encoding chitinase, chiA, and four loci associated with the type three secretion system, ascV, ascFG, aexT, and aexU were sequenced. The phylogenetic reconstruction failed to produce one consensus tree that was compatible with most of the individual genes. The Approximately Unbiased test and the Genetic Algorithm for Recombination Detection both provided further support for differing evolutionary histories among this group of genes. Two contrasting tests detected recombination within aexU, ascFG, ascV, dnaJ, and gyrB but not in aexT or chiA. Quartet decomposition analysis indicated a complex recent evolutionary history for these strains with a high frequency of horizontal gene transfer between several but not among all strains. In this study we demonstrate that at least for some strains, horizontal gene transfer occurs at a sufficient frequency to blur the signal from vertically inherited genes, despite strains being adapted to distinct niches. Simply increasing the number of genes included in the analysis is unlikely to overcome this challenge in organisms that occupy multiple niches and can exchange DNA between strains specialized to different niches. Instead, the detection of genes critical in the adaptation to specific niches may help to reveal the physiological specialization of these strains.