Estimating Genome-wide Phylogenies Using Probabilistic Topic Modeling.

Methods for rapidly inferring the evolutionary history of species or populations with genome-wide data are progressing, but computational constraints still limit our abilities in this area. We developed an alignment-free method to infer genome-wide phylogenies and implemented it in the Python package T opic C ontml . The method uses probabilistic topic modeling (specifically, Latent Dirichlet Allocation or LDA) to extract 'topic' frequencies from -mers, which are derived from multilocus DNA sequences.

The population genetics of speciation by cascade reinforcement.

Species interactions drive diverse evolutionary outcomes. Speciation by cascade reinforcement represents one example of how species interactions can contribute to the proliferation of species. This process occurs when the divergence of mating traits in response to selection against interspecific hybridization incidentally leads to reproductive isolation among populations of the same species.

Genetic diversity and population structure of Leishmania (Viannia) braziliensis in the Peruvian jungle.

Background: Human cutaneous leishmaniasis caused by Leishmania (Viannia) braziliensis is highly prevalent in the Peruvian jungle, where it affects military forces deployed to fight against drug trafficking and civilian people that migrate from the highland to the lowland jungle for economic activities such as mining, agriculture, construction, and chestnut harvest. We explored the genetic diversity and population structure of 124 L. (V.

Fractional coalescent.

An approach to the coalescent, the fractional coalescent (f-coalescent), is introduced. The derivation is based on the discrete-time Cannings population model in which the variance of the number of offspring depends on the parameter α. This additional parameter α affects the variability of the patterns of the waiting times; values of [Formula: see text] lead to an increase of short time intervals, but occasionally allow for very long time intervals.

A synthesis of genetic connectivity in deep-sea fauna and implications for marine reserve design.

With anthropogenic impacts rapidly advancing into deeper waters, there is growing interest in establishing deep-sea marine protected areas (MPAs) or reserves. Reserve design depends on estimates of connectivity and scales of dispersal for the taxa of interest. Deep-sea taxa are hypothesized to disperse greater distances than shallow-water taxa, which implies that reserves would need to be larger in size and networks could be more widely spaced; however, this paradigm has not been tested.

Strong and stable geographic differentiation of swamp buffalo maternal and paternal lineages indicates domestication in the China/Indochina border region.

The swamp type of the Asian water buffalo is assumed to have been domesticated by about 4000 years BP, following the introduction of rice cultivation. Previous localizations of the domestication site were based on mitochondrial DNA (mtDNA) variation within China, accounting only for the maternal lineage. We carried out a comprehensive sampling of China, Taiwan, Vietnam, Laos, Thailand, Nepal and Bangladesh and sequenced the mtDNA Cytochrome b gene and control region and the Y-chromosomal ZFY, SRY and DBY sequences.

Assortative mating drives linkage disequilibrium between sperm and egg recognition protein loci in the sea urchin Strongylocentrotus purpuratus.

Sperm and eggs have interacting proteins on their surfaces that influence their compatibility during fertilization. These proteins are often polymorphic within species, producing variation in gamete affinities. We first demonstrate the fitness consequences of various sperm bindin protein (Bindin) variants in the sea urchin Strongylocentrotus purpuratus, and assortative mating between males and females based on their sperm Bindin genotype.

A continuous method for gene flow.

Most modern population genetics inference methods are based on the coalescence framework. Methods that allow estimating parameters of structured populations commonly insert migration events into the genealogies. For these methods the calculation of the coalescence probability density of a genealogy requires a product over all time periods between events.

Biogeographic models of gene flow in two waterfowl of the Australo-Papuan tropics.

There are many large, easy-to-observe anseriform birds (ducks, geese, and swans) in northern Australia and New Guinea and they often gather in large numbers. Yet, the structure of their populations and their regional movements are poorly understood. Lack of understanding of population structure limits our capacity to understand source-sink dynamics relevant to their conservation or assess risks associated with avian-borne pathogens, in particular, avian influenza for which waterfowl are the main reservoir species.

Balancing a cline by influx of migrants: a genetic transition in water frogs of eastern Greece.

Variation patterns of allozymes and of ND3 haplotypes of mitochondrial DNA reveal a zone of genetic transition among western Palearctic water frogs extending across northeastern Greece and European Turkey. At the western end of the zone, allozymes characteristic of Central European frogs known as Pelophylax ridibundus predominate, whereas at the eastern end, alleles characteristic of western Anatolian water frogs (P. cf.

A possible explanation for the population size discrepancy in tuna (genus Thunnus) estimated from mitochondrial DNA and microsatellite data.

A recent study using both mitochondrial DNA (mtDNA) and microsatellite data reported on a population size discrepancy in the eastern tiger salamander where the effective population size (N(e)) estimate of the former exceeded that of the latter. That study suggested, among other hypotheses, that homoplasy of microsatellite alleles is responsible for the discrepancy. In this investigation, we report 10 new cases of a similar discrepancy in five species of tuna.

Gene flow and population structure in the Mexican blind cavefish complex (Astyanax mexicanus).

Background: Cave animals converge evolutionarily on a suite of troglomorphic traits, the best known of which are eyelessness and depigmentation. We studied 11 cave and 10 surface populations of Astyanax mexicanus in order to better understand the evolutionary origins of the cave forms, the basic genetic structuring of both cave and surface populations, and the degree to which present day migration among them affects their genetic divergence.

Results: To assess the genetic structure within populations and the relationships among them we genotyped individuals at 26 microsatellite loci.

BEAGLE: an application programming interface and high-performance computing library for statistical phylogenetics.

Phylogenetic inference is fundamental to our understanding of most aspects of the origin and evolution of life, and in recent years, there has been a concentration of interest in statistical approaches such as Bayesian inference and maximum likelihood estimation. Yet, for large data sets and realistic or interesting models of evolution, these approaches remain computationally demanding. High-throughput sequencing can yield data for thousands of taxa, but scaling to such problems using serial computing often necessitates the use of nonstatistical or approximate approaches.

Genetic variation of Anastrepha suspensa (Diptera: Tephritidae) in Florida and the Caribbean using microsatellite DNA markers.

Anastrepha suspensa (Loew) (Diptera: Tephritidae), the Caribbean fruit fly, is indigenous to Florida and the Greater Antilles where it causes economic losses in fruit crops, including citrus. Because of the geographic separation of many of its native locations and anecdotal descriptions of regional differences in host preferences, there have been questions about the population structure of A. suspensa.

Phylogeographic patterns of genetic diversity in eastern Mediterranean water frogs have been determined by geological processes and climate change in the Late Cenozoic.

AIM: Our aims were to assess the phylogeographic patterns of genetic diversity in eastern Mediterranean water frogs and to estimate divergence times using different geological scenarios. We related divergence times to past geological events and discuss the relevance of our data for the systematics of eastern Mediterranean water frogs. LOCATION: The eastern Mediterranean region.

Global migration dynamics underlie evolution and persistence of human influenza A (H3N2).

The global migration patterns of influenza viruses have profound implications for the evolutionary and epidemiological dynamics of the disease. We developed a novel approach to reconstruct the genetic history of human influenza A (H3N2) collected worldwide over 1998 to 2009 and used it to infer the global network of influenza transmission. Consistent with previous models, we find that China and Southeast Asia lie at the center of this global network.

Unified framework to evaluate panmixia and migration direction among multiple sampling locations.

For many biological investigations, groups of individuals are genetically sampled from several geographic locations. These sampling locations often do not reflect the genetic population structure. We describe a framework using marginal likelihoods to compare and order structured population models, such as testing whether the sampling locations belong to the same randomly mating population or comparing unidirectional and multidirectional gene flow models.

In defence of model-based inference in phylogeography.

Recent papers have promoted the view that model-based methods in general, and those based on Approximate Bayesian Computation (ABC) in particular, are flawed in a number of ways, and are therefore inappropriate for the analysis of phylogeographic data. These papers further argue that Nested Clade Phylogeographic Analysis (NCPA) offers the best approach in statistical phylogeography. In order to remove the confusion and misconceptions introduced by these papers, we justify and explain the reasoning behind model-based inference.

Evolution of serum albumin intron-1 is shaped by a 5' truncated non-long terminal repeat retrotransposon in western Palearctic water frogs (Neobatrachia).

A 5' truncated non-LTR CR1-like retrotransposon, named RanaCR1, was identified in the serum albumin intron-1 (SAI-1) of at least seven species of western Palearctic water frogs (WPWF). Based on sequence similarity of the carboxy-terminal region (CTR) of ORF2 and/or the highly conserved 3' untranslated region (3' UTR), RanaCR1-like elements occur also in the genome of Xenopus tropicalis and Rana temporaria. Unlike other CR1 elements, RanaCR1 contains a CA microsatellite in its 3' UTR.

Comparison of Bayesian and maximum-likelihood inference of population genetic parameters.

Unlabelled: Comparison of the performance and accuracy of different inference methods, such as maximum likelihood (ML) and Bayesian inference, is difficult because the inference methods are implemented in different programs, often written by different authors. Both methods were implemented in the program MIGRATE, that estimates population genetic parameters, such as population sizes and migration rates, using coalescence theory. Both inference methods use the same Markov chain Monte Carlo algorithm and differ from each other in only two aspects: parameter proposal distribution and maximization of the likelihood function.