Phylogenetic Inference and Neanderthal Mitochondrial DNA: Comparison of Parsimony and Distance Models

Recently, mtDNA was successfully extracted and sequenced from the Neanderthal type specimen (Krings et al, 1997, 1999). Researches attempted to determine the genetic relationship between the Neanderthal specimen and modem human populations using phylogenetic analysis and concluded that the variation existing between the Neanderthal specimen and the modem lineages falls outside the range of variation of modem human populations. Using molecular mutation rate assumptions, it has been concluded that the Neanderthal line diverged from the line leading to modem humans hundreds of thousands of years previous to earlier estimates. This suggests that Neanderthals went extinct without contributing genes to the lineage of modem humans. There are many techniques that can be used in the phylogenetic analysis of molecular data. There is much discussion over the merits of individual techniques and which techniques are best suited for different analysis. I will examine these debates within the framework of the Krings et al. studies and late hominid evolution. Similar analysis was done on the Neanderthal sequences using distance and parsimony methods. A unique database of contemporary human sequences was used. The goals are to test the validity of the results published by Krings et al., and to gain a clearer understanding of the processes of phylogenetic analysis and a greater appreciation of the significance and impact of its results on the field of hominid evolution

Deep-Sea Origin and In-Situ Diversification of Chrysogorgiid Octocorals

The diversity, ubiquity and prevalence in deep waters of the octocoral family Chrysogorgiidae Verrill, 1883 make it noteworthy as a model system to study radiation and diversification in the deep sea. Here we provide the first comprehensive phylogenetic analysis of the Chrysogorgiidae, and compare phylogeny and depth distribution. Phylogenetic relationships among 10 of 14 currently-described Chrysogorgiidae genera were inferred based on mitochondrial (mtMutS, cox1) and nuclear (18S) markers. Bathymetric distribution was estimated from multiple sources, including museum records, a literature review, and our own sampling records (985 stations, 2345 specimens). Genetic analyses suggest that the Chrysogorgiidae as currently described is a polyphyletic family. Shallow-water genera, and two of eight deep-water genera, appear more closely related to other octocoral families than to the remainder of the monophyletic, deep-water chrysogorgiid genera. Monophyletic chrysogorgiids are composed of strictly (Iridogorgia Verrill, 1883, Metallogorgia Versluys, 1902, Radicipes Stearns, 1883, Pseudochrysogorgia Pante & France, 2010) and predominantly (Chrysogorgia Duchassaing & Michelotti, 1864) deep-sea genera that diversified in situ. This group is sister to gold corals (Primnoidae Milne Edwards, 1857) and deep-sea bamboo corals (Keratoisidinae Gray, 1870), whose diversity also peaks in the deep sea. Nine species of Chrysogorgia that were described from depths shallower than 200 m, and mtMutS haplotypes sequenced from specimens sampled as shallow as 101 m, suggest a shallow-water emergence of some Chrysogorgia species

Evolutionary Relationships Among Duiker Antelope (Bovidae: Cephalophinae)

Duikers are a species rich subfamily of threatened African antelope whose recent origin poses a challenge to the molecular identification of taxa and estimation of their phylogeny. I test the ability of DNA barcodes to identify all taxa within this group. I then use mitochondrial and nuclear genes to estimate a multi-locus species tree and to date divergence times. DNA barcodes are unable to distinguish many sister taxa, calling into question the utility of barcodes for the regulation of duiker trade or in identification of field-collected feces. The multi-locus phylogeny provides support for the relationships among major duiker lineages and placement of two problematic taxa, but challenges the validity of the savanna genus and identifies hybridization between taxa. This study reveals that most duikers diverged during the Pleistocene, meriting further inquiry into the role that Pleistocene glacial cycling played in the diversification and population structuring of duikers

Lineage-specific positive selection at the merozoite surface protein 1 (msp1) locus of Plasmodium vivax and related simian malaria parasites

<p>Abstract</p> <p>Background</p> <p>The 200 kDa merozoite surface protein 1 (MSP-1) of malaria parasites, a strong vaccine candidate, plays a key role during erythrocyte invasion and is a target of host protective immune response. <it>Plasmodium vivax</it>, the most widespread human malaria parasite, is closely related to parasites that infect Asian Old World monkeys, and has been considered to have become a parasite of man by host switch from a macaque malaria parasite. Several Asian monkey parasites have a range of natural hosts. The same parasite species shows different disease manifestations among host species. This suggests that host immune responses to <it>P. vivax</it>-related malaria parasites greatly differ among host species (albeit other factors). It is thus tempting to invoke that a major immune target parasite protein such as MSP-1 underwent unique evolution, depending on parasite species that exhibit difference in host range and host specificity.</p> <p>Results</p> <p>We performed comparative phylogenetic and population genetic analyses of the gene encoding MSP-1 (<it>msp1</it>) from <it>P. vivax </it>and nine <it>P. vivax</it>-related simian malaria parasites. The inferred phylogenetic tree of <it>msp1 </it>significantly differed from that of the mitochondrial genome, with a striking displacement of <it>P. vivax </it>from a position close to <it>P. cynomolgi </it>in the mitochondrial genome tree to an outlier of Asian monkey parasites. Importantly, positive selection was inferred for two ancestral branches, one leading to <it>P. inui </it>and <it>P. hylobati </it>and the other leading to <it>P. vivax</it>, <it>P. fieldi </it>and <it>P. cynomolgi</it>. This ancestral positive selection was estimated to have occurred three to six million years ago, coinciding with the period of radiation of Asian macaques. Comparisons of <it>msp1 </it>polymorphisms between <it>P. vivax</it>, <it>P. inui </it>and <it>P. cynomolgi </it>revealed that while some positively selected amino acid sites or regions are shared by these parasites, amino acid changes greatly differ, suggesting that diversifying selection is acting species-specifically on <it>msp1</it>.</p> <p>Conclusions</p> <p>The present results indicate that the <it>msp1 </it>locus of <it>P. vivax </it>and related parasite species has lineage-specific unique evolutionary history with positive selection. <it>P. vivax </it>and related simian malaria parasites offer an interesting system toward understanding host species-dependent adaptive evolution of immune-target surface antigen genes such as <it>msp1</it>.</p

Lutzomyia Sand Fly Diversity and Rates of Infection by Wolbachia and an Exotic Leishmania Species on Barro Colorado Island, Panama

Certain sand fly species living inside or on the edge of tropical forests are well known to transmit a protozoan to humans, which in lowland Panama develops into a cutaneous form of leishmaniasis; open, itching sores on the face and extremities requiring aggressive treatment with antimonial compounds. Morphological characters and DNA sequence from mitochondrial and nuclear gene fragments permitted us to identify and then establish historical relationships among 20 common sand fly species occurring in the understory of Barro Colorado Island, a forested preserve in the middle of the Panama Canal. Individuals in three of these sand fly species were found to be 26–43% infected by Leishmania naiffi, a species hitherto known only from the Amazonian region and the Caribbean. We then screened the same 20 sand fly species for the cytoplasmically transmitted bacteria Wolbachia pipientis, finding three infected at high rates, each by a distinct strain. Lutzomyia trapidoi, the most likely transmitter of Leishmania to humans in Panama, was among the Wolbachia-infected species, thus marking it as a possible high-value target for future biocontrol studies using the bacteria either to induce mating incompatabilities or to drive selected genes into the population

DNA Sequence Classification: It’s Easier Than You Think: An open-source k-mer based machine learning tool for fast and accurate classification of a variety of genomic datasets

Supervised classification of genomic sequences is a challenging, well-studied problem with a variety of important applications. We propose an open-source, supervised, alignment-free, highly general method for sequence classification that operates on k-mer proportions of DNA sequences. This method was implemented in a fully standalone general-purpose software package called Kameris, publicly available under a permissive open-source license. Compared to competing software, ours provides key advantages in terms of data security and privacy, transparency, and reproducibility. We perform a detailed study of its accuracy and performance on a wide variety of classification tasks, including virus subtyping, taxonomic classification, and human haplogroup assignment. We demonstrate the success of our method on whole mitochondrial, nuclear, plastid, plasmid, and viral genomes, as well as randomly sampled eukaryote genomes and transcriptomes. Further, we perform head-to-head evaluations on the tasks of HIV-1 virus subtyping and bacterial taxonomic classification with a number of competing state-of-the-art software solutions, and show that we match or exceed all other tested software in terms of accuracy and speed

Phylogenetic Reconstruction and DNA Barcoding for Closely Related Pine Moth Species (Dendrolimus) in China with Multiple Gene Markers

Unlike distinct species, closely related species offer a great challenge for phylogeny reconstruction and species identification with DNA barcoding due to their often overlapping genetic variation. We tested a sibling species group of pine moth pests in China with a standard cytochrome c oxidase subunit I (COI) gene and two alternative internal transcribed spacer (ITS) genes (ITS1 and ITS2). Five different phylogenetic/DNA barcoding analysis methods (Maximum likelihood (ML)/Neighbor-joining (NJ), “best close match” (BCM), Minimum distance (MD), and BP-based method (BP)), representing commonly used methodology (tree-based and non-tree based) in the field, were applied to both single-gene and multiple-gene analyses. Our results demonstrated clear reciprocal species monophyly for three relatively distant related species, Dendrolimus superans, D. houi, D. kikuchii, as recovered by both single and multiple genes while the phylogenetic relationship of three closely related species, D. punctatus, D. tabulaeformis, D. spectabilis, could not be resolved with the traditional tree-building methods. Additionally, we find the standard COI barcode outperforms two nuclear ITS genes, whatever the methods used. On average, the COI barcode achieved a success rate of 94.10–97.40%, while ITS1 and ITS2 obtained a success rate of 64.70–81.60%, indicating ITS genes are less suitable for species identification in this case. We propose the use of an overall success rate of species identification that takes both sequencing success and assignation success into account, since species identification success rates with multiple-gene barcoding system were generally overestimated, especially by tree-based methods, where only successfully sequenced DNA sequences were used to construct a phylogenetic tree. Non-tree based methods, such as MD, BCM, and BP approaches, presented advantages over tree-based methods by reporting the overall success rates with statistical significance. In addition, our results indicate that the most closely related species D. punctatus, D. tabulaeformis, and D. spectabilis, may be still in the process of incomplete lineage sorting, with occasional hybridizations occurring among them

Core set approach to reduce uncertainty of gene trees

BACKGROUND: A genealogy based on gene sequences within a species plays an essential role in the estimation of the character, structure, and evolutionary history of that species. Because intraspecific sequences are more closely related than interspecific ones, detailed information on the evolutionary process may be available by determining all the node sequences of trees and provide insight into functional constraints and adaptations. However, strong evolutionary correlations on a few lineages make this determination difficult as a whole, and the maximum parsimony (MP) method frequently allows a number of topologies with a same total branching length. RESULTS: Kitazoe et al. developed multidimensional vector-space representation of phylogeny. It converts additivity of evolutionary distances to orthogonality among the vectors expressing branches, and provides a unified index to measure deviations from the orthogoality. In this paper, this index is used to detect and exclude sequences with large deviations from orthogonality, and then selects a maximum subset ("core set") of sequences for which MP generates a single solution. Once the core set tree is formed whose all the node sequences are given, the excluded sequences are found to have basically two phylogenetic positions on this tree, respectively. Fortunately, since multiple substitutions are rare in intra-species sequences, the variance of nucleotide transitions is confined to a small range. By applying the core set approach to 38 partial env sequences of HIV-1 in a single patient and also 198 mitochondrial COI and COII DNA sequences of Anopheles dirus, we demonstrate how consistently this approach constructs the tree. CONCLUSION: In the HIV dataset, we confirmed that the obtained core set tree is the unique maximum set for which MP proposes a single tree. In the mosquito data set, the fluctuation of nucleotide transitions caused by the sequences excluded from the core set was very small. We reproduced this core-set tree by simulation based on random process, and applied our approach to many sets of the obtained endpoint sequences. Consequently, the ninety percent of the endpoint sequences was identified as the core sets and the obtained node sequences were perfectly identical to the true ones