Random amplified polymorphic DNA analysis of DNA extracted from Trichuris trichiura (Linnaeus, 1771) eggs and its prospective application to paleoparasitological studies

Random amplified polymorphic DNA analysis was applied to DNAs extracted from Trichuris trichiura eggs recovered from human fecal samples. Four out of 6 primers tested displayed 18 distinct and well defined polymorphic patterns, ranging from 650 to 3200 base pairs. These results, upon retrieval and DNA sequencing of some of these bands from agarose gels, might help in establishing T. trichiura specific genetic markers, not available yet, and an important step to design primers to be used in molecular diagnosis approaches.

Analysis of ancient DNA from coprolites: a perspective with random amplified polymorphic DNA-polymerase chain reaction approach

The aim of this work was to determine approaches that would improve the quality of ancient DNA (aDNA) present in coprolites to enhance the possibility of success in retrieving specific sequence targets. We worked with coprolites from South American archaeological sites in Brazil and Chile dating up to 7,000 years ago. Using established protocols for aDNA extraction we obtained samples showing high degradation as usually happens with this kind of material. The reconstructive polymerization pretreatment was essential to overcome the DNA degradation and the serial dilutions helped with to prevent polymerase chain reaction (PCR) inhibitors. Moreover, the random amplified polymorphic DNA-PCR has been shown to be a reliable technique for further experiments to recover specific aDNA sequences.

Analysis of the first and second internal transcribed spacer sequences of the ribosomal DNA in Biomphalaria tenagophila complex (Mollusca: Planorbidae)

The first and second internal transcribed spacer regions (ITS1 and ITS2) of the ribosomal DNA of Biomphalaria tenagophila complex (B. tenagophila, B. occidentalis, and B. t. guaibensis) were sequenced and compared. The alignment lengths of these regions were about 655 bp and 481 bp, respectively. Phylogenetic relationships among the Biomphalaria species were inferred by Maximum Parsimony and Neighbor-joining methods. The phylogenetic trees produced, in most of the cases, were in accordance with morphological systematics and other molecular data previously obtained by polymerase chain reaction and restriction fragment length polymorphism analysis. The present results provide support for the proposal that B. tenagophila represents a complex comprising B. tenagophila, B. occidentalis and B. t. guaibensis.

Identification of a genetic marker associated with the resistance to Schistosoma mansoni infection using random amplified polymorphic DNA analysis

In schistosomiasis, the host/parasite interaction remains not completely understood. Many questions related to the susceptibility of snails to infection by respective trematode still remain unanswered. The control of schistosomiasis requires a good understanding of the host/parasite association. In this work, the susceptibility/resistance to Schistosoma mansoni infection within Biomphalaria alexandrina snails were studied starting one month post infection and continuing thereafter weekly up to 10 weeks after miracidia exposure. Genetic variations between susceptible and resistant strains to Schistosoma infection within B. alexandrina snails using random amplified polymorphic DNA analysis technique were also carried out. The results showed that 39.8% of the examined field snails were resistant, while 60.2% of these snails showed high infection rates.In the resistant genotype snails, OPA-02 primer produced a major low molecular weight marker 430 bp. Among the two snail strains there were interpopulational variations, while the individual specimens from the same snail strain, either susceptible or resistant, record semi-identical genetic bands. Also, the resistant character was ascendant in contrast to a decline in the susceptibility of snails from one generation to the next.

Sequence analysis of coding DNA fragments of pfcrt and pfmdr-1 genes in Plasmodium falciparum isolates from Odisha, India

The global emergence and spread of malaria parasites resistant to antimalarial drugs is the major problem in malaria control. The genetic basis of the parasite's resistance to the antimalarial drug chloroquine (CQ) is well-documented, allowing for the analysis of field isolates of malaria parasites to address evolutionary questions concerning the origin and spread of CQ-resistance. Here, we present DNA sequence analyses of both the second exon of the Plasmodium falciparum CQ-resistance transporter (pfcrt) gene and the 5' end of the P. falciparum multidrug-resistance 1 (pfmdr-1) gene in 40 P. falciparum field isolates collected from eight different localities of Odisha, India. First, we genotyped the samples for the pfcrt K76T and pfmdr-1 N86Y mutations in these two genes, which are the mutations primarily implicated in CQ-resistance. We further analyzed amino acid changes in codons 72-76 of the pfcrt haplotypes. Interestingly, both the K76T and N86Y mutations were found to co-exist in 32 out of the total 40 isolates, which were of either the CVIET or SVMNT haplotype, while the remaining eight isolates were of the CVMNK haplotype. In total, eight nonsynonymous single nucleotide polymorphisms (SNPs) were observed, six in the pfcrt gene and two in the pfmdr-1 gene. One poorly studied SNP in the pfcrt gene (A97T) was found at a high frequency in many P. falciparum samples. Using population genetics to analyze these two gene fragments, we revealed comparatively higher nucleotide diversity in the pfcrt gene than in the pfmdr-1 gene. Furthermore, linkage disequilibrium was found to be tight between closely spaced SNPs of the pfcrt gene. Finally, both the pfcrt and the pfmdr-1 genes were found to evolve under the standard neutral model of molecular evolution.

Sources of pre-analytical variations in yield of DNA extracted from blood samples: analysis of 50,000 DNA samples in EPIC.

The European Prospective Investigation into Cancer and nutrition (EPIC) is a long-term, multi-centric prospective study in Europe investigating the relationships between cancer and nutrition. This study has served as a basis for a number of Genome-Wide Association Studies (GWAS) and other types of genetic analyses. Over a period of 5 years, 52,256 EPIC DNA samples have been extracted using an automated DNA extraction platform. Here we have evaluated the pre-analytical factors affecting DNA yield, including anthropometric, epidemiological and technical factors such as center of subject recruitment, age, gender, body-mass index, disease case or control status, tobacco consumption, number of aliquots of buffy coat used for DNA extraction, extraction machine or procedure, DNA quantification method, degree of haemolysis and variations in the timing of sample processing. We show that the largest significant variations in DNA yield were observed with degree of haemolysis and with center of subject recruitment. Age, gender, body-mass index, cancer case or control status and tobacco consumption also significantly impacted DNA yield. Feedback from laboratories which have analyzed DNA with different SNP genotyping technologies demonstrate that the vast majority of samples (approximately 88%) performed adequately in different types of assays. To our knowledge this study is the largest to date to evaluate the sources of pre-analytical variations in DNA extracted from peripheral leucocytes. The results provide a strong evidence-based rationale for standardized recommendations on blood collection and processing protocols for large-scale genetic studies.

DNA multigene sequencing of topotypic specimens of the fascioliasis vector Lymnaea diaphana and phylogenetic analysis of the genus Pectinidens (Gastropoda)

Freshwater lymnaeid snails are crucial in defining transmission and epidemiology of fascioliasis. In South America, human endemic areas are related to high altitudes in Andean regions. The species Lymnaea diaphana has, however, been involved in low altitude areas of Chile, Argentina and Peru where human infection also occurs. Complete nuclear ribosomal DNA 18S, internal transcribed spacer (ITS)-2 and ITS-1 and fragments of mitochondrial DNA 16S and cytochrome c oxidase (cox)1 genes of L. diaphana specimens from its type locality offered 1,848, 495, 520, 424 and 672 bp long sequences. Comparisons with New and Old World Galba/Fossaria, Palaearctic stagnicolines, Nearctic stagnicolines, Old World Radix and Pseudosuccinea allowed to conclude that (i) L. diaphana shows sequences very different from all other lymnaeids, (ii) each marker allows its differentiation, except cox1 amino acid sequence, and (iii) L. diaphana is not a fossarine lymnaeid, but rather an archaic relict form derived from the oldest North American stagnicoline ancestors. Phylogeny and large genetic distances support the genus Pectinidens as the first stagnicoline representative in the southern hemisphere, including colonization of extreme world regions, as most southern Patagonia, long time ago. The phylogenetic link of L. diaphana with the stagnicoline group may give light to the aforementioned peculiar low altitude epidemiological scenario of fascioliasis.

Functional analysis of a conserved DNA methyltransferase in caulobacter crescentus

CcrM is a DNA methyltransferase that methylates the adenine in GANTC motifs in the chromo-some of the bacterial model Caulobacter crescentus. The loss of the CcrM homolog is lethal in C. crescentus and in several other species of Alphaproteobacteria. In this research, we used different experimental and bioinformatic approaches to determine why CcrM is so critical to the physiology of C. crescentus. We first showed that CcrM is a resident orphan DNA methyltransferase in non-Rickettsiales Alphaproteobacteria and that its gene is strictly conserved in this clade (with only one ex¬ception among the genomes sequenced so far). In C. crescentus, cells depleted in CcrM in rich medium quickly lose viability and present an elongated phenotype characteristic of an im¬pairment in cell division. Using minimal medium instead of rich medium as selective and main¬tenance substrate, we could generate a AccrM mutant that presents a viability comparable to the wild type strain and only mild morphological defects. On the basis of a transcriptomic ap¬proach, we determined that several genes essential for cell division were downregulated in the AccrM strain in minimal medium. We offered decisive arguments to support that the efficient transcription of two of these genes, ftsZ and mipZ, coding respectively for the Z-ring forming GTPase FtsZ and an inhibitor of FtsZ polymerization needed for the correct positioning of the Z- ring at mid-cell, requires the methylation of an adenine in a conserved GANTC motif located in their core promoter region. We propose a model, according to which the genome of C. crescentus encodes a transcriptional activator that requires a methylated adenine in a GANTC context to bind to DNA and suggest that this transcriptional regulator might be the global cell-cycle regulator GcrA. In addition, combining a classic genetic approach and in vitro evolution experiments, we showed that the mortality and cell division defects of the AccrM strain in rich medium are mainly due to limiting intracellular levels of the FtsZ protein. We also studied the dynamics of GANTC methylation in C. crescentus using the SMRT technol¬ogy developed by Pacific Biosciences. Our findings support the commonly accepted model, accord¬ing to which the methylation state of GANTC motifs varies during the cell cycle of C. crescentus: before the initiation of DNA replication, the GANTC motifs are fully-methylated (methylated on both strands); when the DNA gets replicated, the GANTC motifs become hemi-methylated (methyl¬ated on one strand only) and this occurs at different times during replication for different loci along the chromosome depending on their position relative to the origin of replication; the GANTC mo¬tifs are only remethylated after DNA replication has finished as a consequence of the massive and short-lived expression of CcrM in predivisional cells. About 30 GANTC motifs in the C. crescentus chromosome were found to be undermethylated in most of the bacterial population; these might be protected from CcrM activity by DNA binding proteins and some of them could be involved in methylation-based bistable transcriptional switches. - CcrM est une ADN méthyltransférase qui méthyle les adénines dans le contexte GANTC dans le génome de la bactérie modèle Caulobacter crescentus. La perte de l'homologue de CcrM chez C. crescentus et chez plusieurs autres espèces d'Alphaproteobactéries est létale. Dans le courant de cette recherche, nous tentons de déterminer pourquoi la protéine CcrM est cruciale pour la survie de C. crescentus. Nous démontrons d'abord que CcrM est une adénine méthyltransférase orpheline résidente, dont le gène fait partie du génome minimal partagé par les Alphaprotéobactéries non-Rickettsiales (à une exception près). Lorsqu'une souche de C. crescentus est privée de CcrM, sa viabilité décroît rapi¬dement et ses cellules présentent une morphologie allongée qui suggère que la division cellulaire est inhibée. Nous sommes parvenus à créer une souche AccrM en utilisant un milieu minimum, au lieu du milieu riche classiquement employé, comme milieu de sélection et de maintenance pour la souche. Lorsque nous avons étudié le transcriptome de cette souche de C. crescentus privée de CcrM, nous avons pu constater que plusieurs gènes essentiels pour le bon déroulement de la division cellulaire bactérienne étaient réprimés. En particulier, l'expression adéquate des gènes ftsZ et mipZ - qui codent, respectivement, pour FtsZ, la protéine qui constitue, au milieu de la cellule, un anneau protéique qui initie le processus de division et pour MipZ, un inhibiteur de la polymérisation de FtsZ qui est indispensable pour le bon positionnement de l'anneau FtsZ - est dépendante de la présence d'une adénine méthylée dans un motif GANTC conservé situé dans leur région promotrice. Nous présentons un modèle selon lequel le génome de C. crescentus code pour un facteur de transcription qui exige la présence d'une adénine méthylée dans un contexte GANTC pour s'attacher à l'ADN et nous suggérons qu'il pourrait s'agir du régulateur global du cycle cellulaire GcrA. En outre, nous montrons, en combinant la génétique classique et une approche basée sur l'évolution expérimentale, que la mortalité et l'inhibition de la division cellulaire caractéristiques de la souche àccrMeη milieu riche sont dues à des niveaux excessivement bas de protéine FtsZ. Nous avons aussi étudié la dynamique de la méthylation du chromosome de C. crescentus sur la base de la technologie SMRT développée par Pacific Biosciences. Nous confirmons le modèle communément accepté, qui affirme que l'état de méthylation des motifs GANTC change durant le cycle cellulaire de C. crescentus: les motifs GANTC sont complètement méthylés (méthylés sur les deux brins) avant de début de la réplication de l'ADN; ils deviennent hémi-méthylés (méthylés sur un brin seulement) une fois répliqués, ce qui arrive à différents moments durant la réplication pour différents sites le long du chromosome en fonction de leur position par rapport à l'origine de répli-cation; finalement, les motifs GANTC sont reméthylés après la fin de la réplication du chromosome lorsque la protéine CcrM est massivement, mais très transitoirement, produite. Par ailleurs, nous identifions dans le chromosome de C. crescentus environ 30 motifs GANTC qui restent en perma-nence non-méthylés dans une grande partie de la population bactérienne; ces motifs sont probable-ment protégés de l'action de CcrM par des protéines qui s'attachent à l'ADN et certains d'entre eux pourraient être impliqués dans des mécanismes de régulation générant une transcription bistable.

Functional analysis of human POT1 and RPA : insights into the role of single stranded DNA binding proteins at telomeres

Abstract Telomeres, the natural ends of chromosomes, need to be protected from chromosome end fusions, aberrant homologous recombination and degradation. In humans, chromosome ends are specified through arrays of tandemly repeated 5'-TTAGGG-3' hexamers, ending in a 3' overhang. A complex formed by the six proteins TRF1, TRF2, hRap1, TIN2, TPP1 and POT1 specifically assocìates with and protects telomeres. Telomeres are maintained by semiconservative DNA replication and by a specialized reverse transcriptase, telomerase, that carries an RNA subunit which templates new telomeric repeat synthesis. The telomeric single stranded (ss) DNA binding protein POT1 protects the telomeric 3' overhang and modulates telomerase-mediated telomere elongation. It is possible that POT1 also influences DNA synthesis during semiconservative DNA replication, which is initiated by the DNA polymerase alpha-primase complex. The heterotrimeric ss DNA-binding protein RPA plays essential roles during DNA replication. RPA binds to ss DNA with high affinity in order to stabilize ss DNA and facilitate nascent strand synthesis at the replication fork. Here we investigate how the two proteins RPA and POT1 contribute to telomere maintenance by regulating semi-conservative DNA replication and telomerase. Using chromatin immunoprecipitation experiments, we show that RPA associates with telomeres during S-phase. Analysis of telomere structure in cells shRNA-depleted for RPA and POT1 reveals that loss of RPA and POT1 causes exposure of single-stranded DNA at telomeres, suggestive of incomplete DNA replication. Biochemical experiments using purified recombinant POT1 and RPA show that saturating telomeric oligonucleotides with POT1 or RPA reduces the primase activity of the DNA polymerase alpha-primase complex and the overall activity of telomerase. POT1 and RPA also increase the primer extension by DNA polymerase alpha-primase complex and the processivity of telomerase under certain conditions, although POT1 increases the activities to a greater extent than RPA. We propose that POT1 is required for proper replication of the lagging strand of telomeres and that some phenotypes observed in POT1-depleted cells may stern from incomplete DNA replication rather than de-protection of the single-stranded overhang. Résumé Les télomères, les extrémités normales des chromosomes linéaires, doivent être protégés des fusions chromosomiques, d'événements de recombinaison homologue aberrants et de phénomènes de dégradation. Chez l'Homme, les extrémités des chromosomes sont constitués d'ADN double brin répétitif de séquence 5'-TTAGGG-3', d'une extension simple brin 3' sortante et d'un complexe protéique formé des six facteurs TRF1, TRF2, hRap1, TIN2, TPP1 et POT1 qui, s'associant à cette séquence, protègent l'ADN télomèrique. Les télomères sont maintenus par la télomérase, une transcriptase inverse capable d'allonger l'extension 3' sortante télomérique. POT1 lie l'ADN simple brin télomérique et module l'élongation des télomères par la télomérase. POT1 pourrait en théorie également influencer la réplication semi-conservative de l'ADN. L'ADN-polymérase Pal alpha-primase amorce et initie la synthèse d'ADN. Pendant la réplication, l'ADN simple brin est stabilisé par RPA, un complexe hétérotrimèrique qui lie l'ADN simple brin. RPA facilite la synthèse du brin naissant à la fourche de réplication. Ici nous avons étudié comment ces deux protéines qui lient l'ADN simple brin, RPA et POT1, régulent la réplication des télomères par la télomérase et la machinerie classique de réplication de l'ADN. Par immunoprécipitation de chromatine (ChIP), nous montrons que RPA est localisé aux télomères lors de la phase S du cycle cellulaire. De plus, l'analyse de la structure des télomeres indique que !a perte de RPA ou de POT1 conduit à l'apparition d'ADN simple brin télomérique, suggérant une réplication incomplète de l'ADN télomérique in vivo. Par une approche complémentaire biochimique utilisant les protéines POT1 et RPA recombinantes purifiées, nous montrons également que la liaison de POT1 ou de RPA à des oligonucléotides télomériques bloque l'activité primase du complexe polymérase alpha/primase et réduit l'activité télomérase sur ces substrats. En revanche, leur liaison augmente l'activité ADN-polymérase du complexe polymérase alpha/primase, ainsi que fa processivité de la télomérase dans certaines conditions, POT1 étant le plus efficace des deux facteurs. Nous proposons que POT1 est nécessaire à la réplication du brin retardé au niveau des télomères, ce qui suggère que certains phénotypes des cellules déplétés en POT1 puissent résulter d'une réplication incomplète de l'ADN télémétrique plutôt que d'une déprotection de l'extrémité sortante des télomères.

Meta-analysis of estrogen response in MCF-7 distinguishes early target genes involved in signaling and cell proliferation from later target genes involved in cell cycle and DNA repair.

ABSTRACT: BACKGROUND: Many studies have been published outlining the global effects of 17 beta-estradiol (E2) on gene expression in human epithelial breast cancer derived MCF-7 cells. These studies show large variation in results, reporting between ~100 and ~1500 genes regulated by E2, with poor overlap. RESULTS: We performed a meta-analysis of these expression studies, using the Rank product method to obtain a more accurate and stable list of the differentially expressed genes, and of pathways regulated by E2. We analyzed 9 time-series data sets, concentrating on response at 3-4 hrs (early) and at 24 hrs (late). We found >1000 statistically significant probe sets after correction for multiple testing at 3-4 hrs, and >2000 significant probe sets at 24 hrs. Differentially expressed genes were examined by pathway analysis. This revealed 15 early response pathways, mostly related to cell signaling and proliferation, and 20 late response pathways, mostly related to breast cancer, cell division, DNA repair and recombination. CONCLUSIONS: Our results show that meta-analysis identified more differentially expressed genes than the individual studies, and that these genes act together in networks. These results provide new insight into E2 regulated mechanisms, especially in the context of breast cancer.

Mutational and computational analysis of the alpha(1b)-adrenergic receptor. Involvement of basic and hydrophobic residues in receptor activation and G protein coupling.

To investigate their role in receptor coupling to G(q), we mutated all basic amino acids and some conserved hydrophobic residues of the cytosolic surface of the alpha(1b)-adrenergic receptor (AR). The wild type and mutated receptors were expressed in COS-7 cells and characterized for their ligand binding properties and ability to increase inositol phosphate accumulation. The experimental results have been interpreted in the context of both an ab initio model of the alpha(1b)-AR and of a new homology model built on the recently solved crystal structure of rhodopsin. Among the twenty-three basic amino acids mutated only mutations of three, Arg(254) and Lys(258) in the third intracellular loop and Lys(291) at the cytosolic extension of helix 6, markedly impaired the receptor-mediated inositol phosphate production. Additionally, mutations of two conserved hydrophobic residues, Val(147) and Leu(151) in the second intracellular loop had significant effects on receptor function. The functional analysis of the receptor mutants in conjunction with the predictions of molecular modeling supports the hypothesis that Arg(254), Lys(258), as well as Leu(151) are directly involved in receptor-G protein interaction and/or receptor-mediated activation of the G protein. In contrast, the residues belonging to the cytosolic extensions of helices 3 and 6 play a predominant role in the activation process of the alpha(1b)-AR. These findings contribute to the delineation of the molecular determinants of the alpha(1b)-AR/G(q) interface.

Using graphical probability analysis (Bayes Nets) to evaluate a conditional DNA inclusion

This paper discusses the analysis of cases in which the inclusion or exclusion of a particular suspect, as a possible contributor to a DNA mixture, depends on the value of a variable (the number of contributors) that cannot be determined with certainty. It offers alternative ways to deal with such cases, including sensitivity analysis and object-oriented Bayesian networks, that separate uncertainty about the inclusion of the suspect from uncertainty about other variables. The paper presents a case study in which the value of DNA evidence varies radically depending on the number of contributors to a DNA mixture: if there are two contributors, the suspect is excluded; if there are three or more, the suspect is included; but the number of contributors cannot be determined with certainty. It shows how an object-oriented Bayesian network can accommodate and integrate varying perspectives on the unknown variable and how it can reduce the potential for bias by directing attention to relevant considerations and distinguishing different sources of uncertainty. It also discusses the challenge of presenting such evidence to lay audiences.

Understanding mutational and selective processes that govern gene duplication in mammals

Abstract : Gene duplication is an essential source of material for the origin of genetic novelties. The reverse transcription of source gene mRNA followed by the genomic insertion of the resulting cDNA - retroposition - has provided the human genome with at least ~3600 detectable retrocopies. We find that ~30% of these retrocopies are transcribed, generally in testes. Their transcription often relies on preexisting regulatory elements (or open chromatin) close to their insertion site, which is illustrated by mRNA molecules containing retrocopies fused to their neighboring genes. Retrocopies appear to have been profoundly shaped by selection. Consistently, human retrocopies with an intact open reading (ORF) are more often transcribed than retropseudogenes, which leads to a minimal estimate of 120 functional retrogenes present in our genome. We also performed an analysis of Ka/Ks for human retrocopies. This analysis demonstrates that several intact retrocopies evolved under purifying selection and yields an estimated formation rate of ~1 retrogene per million year in the primate lineage. Using DNA sequencing and evolutionary simulations, we have identified 7 such primate-specific retrogenes that emerged on the lineage leading to humans In therian genomes, we found an excess of retrogenes with X-linked parents. Expression analyses support the idea that this "out of X" movement was driven by natural selection to produce autosomal functional counterparts for X-linked genes, which are silenced during male meiosis. Phylogenetic dating of this "out of X" movement suggests that our sex chromosomes arose about 180 MYA ago and are thus much younger than previously thought. Finally, we have also analyzed young gene duplications (and deletions) that arose by non allelic-homologous recombination and are not fixed in species. Using wild-caught and laboratory animals, we detected thousands of DNA segments that are polymorphic in copy number in mice. These copy number variants were found to profoundly alter the transcriptome of several mouse tissues. Strikingly, their influence on gene expression is not limited to the gene they contain but seems to extend to genes located up to 1.5 million bases away.

Experimental analysis and computer prediction of CTF/NFI transcription factor DNA binding sites.

Accurate prediction of transcription factor binding sites is needed to unravel the function and regulation of genes discovered in genome sequencing projects. To evaluate current computer prediction tools, we have begun a systematic study of the sequence-specific DNA-binding of a transcription factor belonging to the CTF/NFI family. Using a systematic collection of rationally designed oligonucleotides combined with an in vitro DNA binding assay, we found that the sequence specificity of this protein cannot be represented by a simple consensus sequence or weight matrix. For instance, CTF/NFI uses a flexible DNA binding mode that allows for variations of the binding site length. From the experimental data, we derived a novel prediction method using a generalised profile as a binding site predictor. Experimental evaluation of the generalised profile indicated that it accurately predicts the binding affinity of the transcription factor to natural or synthetic DNA sequences. Furthermore, the in vitro measured binding affinities of a subset of oligonucleotides were found to correlate with their transcriptional activities in transfected cells. The combined computational-experimental approach exemplified in this work thus resulted in an accurate prediction method for CTF/NFI binding sites potentially functioning as regulatory regions in vivo.

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

DnaSP 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.