Molecular Clock Model Research Articles

The substitution rate of HIV-1 subtypes: a genomic approach.

HIV-1M causes most infections in the AIDS pandemic. Its genetic diversity is defined by nine pure subtypes and more than sixty recombinant forms. We have performed a comparative analysis of the evolutionary rate of five pure subtypes (A1, B, C, D, and G) and two circulating recombinant forms (CRF01_AE and CRF02 AG) using data obtained from nearly complete genome coding sequences. Times to the most recent common ancestor (tMRCA) and substitution rates of these HIV genomes, and their genomic partitions, were estimated by Bayesian coalescent analyses. Genomic substitution rate estimates were compared between the HIV-1 datasets analyzed by means of randomization tests. Significant differences in the rate of evolution were found between subtypes, with subtypes C and A1 and CRF01_AE displaying the highest rates. On the other hand, CRF02_AG and subtype D were the slowest evolving types. Using a different molecular clock model for each genomic partition led to more precise tMRCA estimates than when linking the same clock along the HIV genome. Overall, the earliest tMRCA corresponded to subtype A1 (median = 1941, 95% HPD = 1943–55), whereas the most recent tMRCA corresponded to subtype G and CRF01_AE subset 3 (median = 1971, 95% HPD = 1967–75 and median = 1972, 95% HPD = 1970–75, respectively). These results suggest that both biological and epidemiological differences among HIV-1M subtypes are reflected in their evolutionary dynamics. The estimates obtained for tMRCAs and substitution rates provide information that can be used as prior distributions in future Bayesian coalescent analyses of specific HIV-1 subtypes/CRFs and genes.

Molecular phylogeography ofTroglophiluscave crickets (Orthoptera, Rhaphidophoridae): A combination of vicariance and dispersal drove diversification in the East Mediterranean region

In this study, we investigated the molecular phylogenetic divergence and historical biogeography of cave crickets belonging to the genus Troglophilus (Orthoptera, Rhaphidophoridae) from caves in eastern Mediterranean and Anatolia regions. Three mitochondrial DNA genes (COI, 12S rDNA, and 16S rDNA) and two nuclear ones (18S rDNA and 28S rDNA) were amplified and partially sequenced to reconstruct phylogenetic relationships among most of the known Troglophilus species. Results showed a well-resolved phylogeny with three main clades representing the Balkan, the Anatolian, and the Cycladian–Cretan lineages. Based on Bayesian analyses, we applied a relaxed molecular clock model to estimate the divergence times between these three lineages. Dating estimates indicate that radiation of the ingroup might have been triggered by the opening of the Mid-Aegean trench, while the uplift of the Anatolian Plateau in Turkey and the changes of relief, emergence, and disappearance of orographic and hydrographical barriers in the Balkan Peninsula are potential paleogeographic events responsible for the initial diversification of the genus Troglophilus. A possible biogeographic scenario, reconstructed using S-DIVA with RASP software, suggested that the current distribution of Troglophilus species can be explained by a combination of both dispersal and vicariance events that occurred in particular in the ancestral populations. The radiation of Troglophilus species likely started from the Aegean and proceeded eastward to Anatolia and westward to the Balkan region. Results are additionally compared to those available for Dolichopoda, the only other representative genus of Rhaphidophoridae present in the Mediterranean area.

Evolutionary analysis of whole-genome sequences confirms inter-farm transmission of Aleutian mink disease virus.

Aleutian mink disease virus (AMDV) is a frequently encountered pathogen associated with mink farming. Previous phylogenetic analyses of AMDV have been based on shorter and more conserved parts of the genome, e.g. the partial NS1 gene. Such fragments are suitable for detection but are less useful for elucidating transmission pathways while sequencing entire viral genomes provides additional informative sites and often results in better-resolved phylogenies. We explore how whole-genome sequencing can benefit investigations of AMDV transmission by reconstructing the relationships between AMDV field samples from a Danish outbreak. We show that whole-genome phylogenies are much better resolved than those based on the partial NS1 gene sequences extracted from the same alignment. Well-resolved phylogenies contain more information about the underlying transmission trees and are useful for understanding the spread of a pathogen. In the main case investigated here, the transmission path suggested by the tree structure was supported by epidemiological data. The use of molecular clock models further improved tree resolution and provided time estimates for the viral ancestors consistent with the proposed direction of spread. It was however impossible to infer transmission pathways from the partial NS1 gene tree, since all samples from the case farms branched out from a single internal node. A sliding window analysis showed that there were no shorter genomic regions providing the same phylogenetic resolution as the entire genome. Altogether, these results suggest that phylogenetic analyses based on whole-genome sequencing taking into account sampling dates and epidemiological data is a promising set of tools for clarifying AMDV transmission.

RelTime Rates Collapse to a Strict Clock When Estimating the Timeline of Animal Diversification.

Establishing an accurate timescale for the history of life is crucial to understand evolutionary processes. For this purpose, relaxed molecular clock models implemented in a Bayesian MCMC framework are generally used. However, these methods are time consuming. RelTime, a non-Bayesian method implementing a fast, ad hoc, algorithm for relative dating, was developed to overcome the computational inefficiencies of Bayesian software. RelTime was recently used to investigate the timing of origin of animals, yielding results consistent with early strict clock studies from the 1980s and 1990s, estimating metazoans to have a Mesoproterozoic origin—over a billion years ago. RelTime results are unexpected and disagree with the largest majority of modern, relaxed, Bayesian molecular clock analyses, which suggest animals originated in the Tonian-Cryogenian (less that 850 million years ago). Here, we demonstrate that RelTime-inferred divergence times for the origin of animals are spurious, a consequence of the inability of RelTime to relax the clock along the internal branches of the animal phylogeny. RelTime-inferred divergence times are comparable to strict-clock estimates because they are essentially inferred under a strict clock. Our results warn us of the danger of using ad hoc algorithms making implicit assumptions about rate changes along a tree. Our study roundly rejects a Mesoproterozoic origin of animals; metazoans emerged in the Tonian-Cryogenian, and diversified in the Ediacaran, in the immediate prelude to the routine fossilization of animals in the Cambrian associated with the emergence of readily preserved skeletons.

Bayesian inference reveals ancient origin of simian foamy virus in orangutans

Simian foamy viruses (SFVs) infect most nonhuman primate species and appears to co-evolve with its hosts. This co-evolutionary signal is particularly strong among great apes, including orangutans (genus Pongo). Previous studies have identified three distinct orangutan SFV clades. The first of these three clades is composed of SFV from P. abelii from Sumatra, the second consists of SFV from P. pygmaeus from Borneo, while the third clade is mixed, comprising an SFV strain found in both species of orangutan. The existence of the mixed clade has been attributed to an expansion of P. pygmaeus into Sumatra following the Mount Toba super-volcanic eruption about 73,000years ago. Divergence dating, however, has yet to be performed to establish a temporal association with the Toba eruption. Here, we use a Bayesian framework and a relaxed molecular clock model with fossil calibrations to test the Toba hypothesis and to gain a more complete understanding of the evolutionary history of orangutan SFV. As with previous studies, our results show a similar three-clade orangutan SFV phylogeny, along with strong statistical support for SFV-host co-evolution in orangutans. Using Bayesian inference, we date the origin of orangutan SFV to >4.7 million years ago (mya), while the mixed species clade dates to approximately 1.7mya, >1.6 million years older than the Toba super-eruption. These results, combined with fossil and paleogeographic evidence, suggest that the origin of SFV in Sumatran and Bornean orangutans, including the mixed species clade, likely occurred on the mainland of Indo-China during the Late Pliocene and Calabrian stage of the Pleistocene, respectively.

Historical biogeography and mitogenomics of two endemic Mediterranean gorgonians (Holaxonia, Plexauridae)

Among the Mediterranean plexaurids, four species are endemic and despite their ecological importance, comprehensive studies on the evolution and biogeography of these organisms are lacking. Here, we explore the mitogenomic variability of two endemic, ecologically important Mediterranean Paramuricea species. We assess their phylogenetic relationships and provide first insights into their evolution and biogeography. Complete mitogenome sequences of Paramuricea clavata and Paramuricea macrospina were obtained using long-range PCR, primer-walking and Sanger sequencing. For an enlarged sample of Paramuricea species, maximum likelihood and Bayesian phylogenetic trees of the mitochondrial gene mtMutS were obtained and used to study the biogeographic history of Paramuricea through a statistical Dispersal-Vicariance (S-DIVA) method and a Dispersal Extinction Cladogenesis (DEC) model. Divergence time was estimated under strict and relaxed molecular clock models in BEAST using published octocoral mutation rates. Our results revealed high nucleotide diversity (2.6%) among the two Mediterranean endemics; the highest mutation rates were found in the mtMutS, Nad4 and Nad5. In addition, we found length polymorphisms in several intergenic regions and differences in mitochondrial genome size. The red gorgonian P. clavata was closely related to the Eastern Atlantic Paramuricea grayi rather than its Mediterranean congener, P. macrospina. Our biogeographic results provide evidence for the independent speciation of the Mediterranean species and point to a Miocene origin of the two endemics, highlighting the role played by the Messinian Salinity Crisis in the evolutionary history of Mediterranean organisms.

Comparative genomics provides a timeframe for Wolbachia evolution and exposes a recent biotin synthesis operon transfer.

The genus Wolbachia (Alphaproteobacteria) comprises the most abundant inherited intracellular bacteria1. Despite their relevance as manipulators of human pathogen transmission2 and arthropod reproduction3, many aspects of their evolutionary history are not well understood4. In arthropods, Wolbachia infections are typically transient on evolutionary timescales5,6 and co-divergence between hosts and Wolbachia is supposedly rare. Consequently, much of our knowledge of Wolbachia genome evolution derives from very recently diverged strains, and a timescale for Wolbachia is lacking. Here, we investigated the genomes of four Wolbachia strains that have persisted within and co-diverged with their host lineage for ∼2 million years. Although the genomes showed very little evolutionary change on a nucleotide level, we found evidence for a recent lateral transfer of a complete biotin synthesis operon that has the potential to transform Wolbachia-host relationships7. Furthermore, this evolutionary snapshot enabled us to calibrate the divergence times of the supergroup A and B Wolbachia lineages using genome-wide data sets and relaxed molecular clock models. We estimated the origin of Wolbachia supergroups A and B to be ∼200 million years ago (Ma), which is considerably older than previously appreciated. This age coincides with the diversification of many insect lineages8 that represent most of Wolbachia's host spectrum.

A Statistical Framework to Identify Deviation from Time Linearity in Epigenetic Aging.

In multiple studies DNA methylation has proven to be an accurate biomarker of age. To develop these biomarkers, the methylation of multiple CpG sites is typically linearly combined to predict chronological age. By contrast, in this study we apply the Universal PaceMaker (UPM) model to investigate changes in DNA methylation during aging. The UPM was initially developed to study rate acceleration/deceleration in sequence evolution. Rather than identifying which linear combinations of sites predicts age, the UPM models the rates of change of multiple CpG sites, as well as their starting methylation levels, and estimates the age of each individual to optimize the model fit. We refer to the estimated age as the “epigenetic age”, which is in contrast to the known chronological age of each individual. We construct a statistical framework and devise an algorithm to determine whether a genomic pacemaker is in effect (i.e rates of change vary with age). The decision is made by comparing two competing likelihood based models, the molecular clock (MC) and UPM. For the molecular clock model, we use the known chronological age of each individual and fit the methylation rates at multiple sites, and express the problem as a linear least squares and solve it in polynomial time. For the UPM case, the search space is larger as we are fitting both the epigenetic age of each individual as well as the rates for each site, yet we succeed to reduce the problem to the space of individuals and polynomial in the more significant space—the methylated sites. We first tested our algorithm on simulated data to elucidate the factors affecting the identification of the pacemaker model. We find that, provided with enough data, our algorithm is capable of identifying a pacemaker even when a weak signal is present in the data. Based on these results, we applied our method to DNA methylation data from human blood from individuals of various ages. Although the improvement in variance across sites between the UPM and MC was small, the results suggest that the existence of a pacemaker is highly significant. The PaceMaker results also suggest a decay in the rate of change in DNA methylation with age.

Concordant patterns of mtDNA and nuclear phylogeographic structure reveal Pleistocene vicariant event in the green crab Carcinus aestuarii across the Siculo-Tunisian Strait

This study focuses on the population genetic structure of the green crab Carcinus aestuarii along part of the African Mediterranean coast, with the main target to confirm genetic subdivision across the well documented genetic boundary of the Siculo-Tunisian Strait. For this purpose, the mitochondrial COI (cytochrome oxidase I) gene and five polymorphic microsatellite loci were analysed in 144 and 120 specimens, respectively. Our results show the existence of two distinct haplogroups, separated by 16 mutational steps and revealed a non random distribution of the genetic variation along the African Mediterranean coast. Dating analyses, based on the use of different molecular clock models and rates, placed the divergence among both haplogroups at 1.91 Myr (95% HPD: 1.11–2.68 Myr) to 0.69 Myr (95% HPD: 0.44–0.98 Myr). This range of divergence time estimation corresponds to the Early Pleistocene. The particular pattern of genetic divergence among Eastern and Western African Mediterranean populations of C. aestuarii, detected by 2-level AMOVA at the mitochondrial level, was consistent with that inferred from microsatellite analysis and suggests a vicariant event in C. aestuarii. Demographic reconstruction, inferred from mismatch distribution and BSP analyses, yielded different patterns of demographic history between both African Mediterranean groups. The distribution pattern of the two haplogroups across the African Mediterranean coast, along with results of Bayesian analysis of genetic structure revealing an intermediate geographic group between the two divergent groups of the African coast, support the hypothesis of a secondary contact between two historically isolated groups. Although this hypothetical contact zone, thought to be located near the Siculo-Tunisian Strait, still needs to be verified, the asymmetric gene flow from Western to Eastern African Mediterranean, as inferred by the results of a MIGRATE analysis, reinforces the previously mentioned results.

Cross-validation to select Bayesian hierarchical models in phylogenetics

BackgroundRecent developments in Bayesian phylogenetic models have increased the range of inferences that can be drawn from molecular sequence data. Accordingly, model selection has become an important component of phylogenetic analysis. Methods of model selection generally consider the likelihood of the data under the model in question. In the context of Bayesian phylogenetics, the most common approach involves estimating the marginal likelihood, which is typically done by integrating the likelihood across model parameters, weighted by the prior. Although this method is accurate, it is sensitive to the presence of improper priors. We explored an alternative approach based on cross-validation that is widely used in evolutionary analysis. This involves comparing models according to their predictive performance.ResultsWe analysed simulated data and a range of viral and bacterial data sets using a cross-validation approach to compare a variety of molecular clock and demographic models. Our results show that cross-validation can be effective in distinguishing between strict- and relaxed-clock models and in identifying demographic models that allow growth in population size over time. In most of our empirical data analyses, the model selected using cross-validation was able to match that selected using marginal-likelihood estimation. The accuracy of cross-validation appears to improve with longer sequence data, particularly when distinguishing between relaxed-clock models.ConclusionsCross-validation is a useful method for Bayesian phylogenetic model selection. This method can be readily implemented even when considering complex models where selecting an appropriate prior for all parameters may be difficult.Electronic supplementary materialThe online version of this article (doi:10.1186/s12862-016-0688-y) contains supplementary material, which is available to authorized users.

Evolution of a Heavy Metal Homeostasis/Resistance Island Reflects Increasing Copper Stress in Enterobacteria.

Copper homeostasis in bacteria is challenged by periodic elevation of copper levels in the environment, arising from both natural sources and human inputs. Several mechanisms have evolved to efflux copper from bacterial cells, including the cus (copper sensing copper efflux system), and pco (plasmid-borne copper resistance system) systems. The genes belonging to these two systems can be physically clustered in a Copper Homeostasis and Silver Resistance Island (CHASRI) on both plasmids and chromosomes in Enterobacteria. Increasing use of copper in agricultural and industrial applications raises questions about the role of human activity in the evolution of novel copper resistance mechanisms. Here we present evidence that CHASRI emerged and diversified in response to copper deposition across aerobic and anaerobic environments. An analysis of diversification rates and a molecular clock model suggest that CHASRI experienced repeated episodes of elevated diversification that could correspond to peaks in human copper production. Phylogenetic analyses suggest that CHASRI originated in a relative of Enterobacter cloacae as the ultimate product of sequential assembly of several pre-existing two-gene modules. Once assembled, CHASRI dispersed via horizontal gene transfer within Enterobacteriaceae and also to certain members of Shewanellaceae, where the original pco module was replaced by a divergent pco homolog. Analyses of copper stress mitigation suggest that CHASRI confers increased resistance aerobically, anaerobically, and during shifts between aerobic and anaerobic environments, which could explain its persistence in facultative anaerobes and emergent enteric pathogens.

Exploring the temporal structure of heterochronous sequences using TempEst (formerly Path-O-Gen).

Gene sequences sampled at different points in time can be used to infer molecular phylogenies on a natural timescale of months or years, provided that the sequences in question undergo measurable amounts of evolutionary change between sampling times. Data sets with this property are termed heterochronous and have become increasingly common in several fields of biology, most notably the molecular epidemiology of rapidly evolving viruses. Here we introduce the cross-platform software tool, TempEst (formerly known as Path-O-Gen), for the visualization and analysis of temporally sampled sequence data. Given a molecular phylogeny and the dates of sampling for each sequence, TempEst uses an interactive regression approach to explore the association between genetic divergence through time and sampling dates. TempEst can be used to (1) assess whether there is sufficient temporal signal in the data to proceed with phylogenetic molecular clock analysis, and (2) identify sequences whose genetic divergence and sampling date are incongruent. Examination of the latter can help identify data quality problems, including errors in data annotation, sample contamination, sequence recombination, or alignment error. We recommend that all users of the molecular clock models implemented in BEAST first check their data using TempEst prior to analysis.

Fast Dating Using Least-Squares Criteria and Algorithms.

Phylogenies provide a useful way to understand the evolutionary history of genetic samples, and data sets with more than a thousand taxa are becoming increasingly common, notably with viruses (e.g., human immunodeficiency virus (HIV)). Dating ancestral events is one of the first, essential goals with such data. However, current sophisticated probabilistic approaches struggle to handle data sets of this size. Here, we present very fast dating algorithms, based on a Gaussian model closely related to the Langley–Fitch molecular-clock model. We show that this model is robust to uncorrelated violations of the molecular clock. Our algorithms apply to serial data, where the tips of the tree have been sampled through times. They estimate the substitution rate and the dates of all ancestral nodes. When the input tree is unrooted, they can provide an estimate for the root position, thus representing a new, practical alternative to the standard rooting methods (e.g., midpoint). Our algorithms exploit the tree (recursive) structure of the problem at hand, and the close relationships between least-squares and linear algebra. We distinguish between an unconstrained setting and the case where the temporal precedence constraint (i.e., an ancestral node must be older that its daughter nodes) is accounted for. With rooted trees, the former is solved using linear algebra in linear computing time (i.e., proportional to the number of taxa), while the resolution of the latter, constrained setting, is based on an active-set method that runs in nearly linear time. With unrooted trees the computing time becomes (nearly) quadratic (i.e., proportional to the square of the number of taxa). In all cases, very large input trees (>10,000 taxa) can easily be processed and transformed into time-scaled trees. We compare these algorithms to standard methods (root-to-tip, r8s version of Langley–Fitch method, and BEAST). Using simulated data, we show that their estimation accuracy is similar to that of the most sophisticated methods, while their computing time is much faster. We apply these algorithms on a large data set comprising 1194 strains of Influenza virus from the pdm09 H1N1 Human pandemic. Again the results show that these algorithms provide a very fast alternative with results similar to those of other computer programs. These algorithms are implemented in the LSD software (least-squares dating), which can be downloaded from http://www.atgc-montpellier.fr/LSD/, along with all our data sets and detailed results. An Online Appendix, providing additional algorithm descriptions, tables, and figures can be found in the Supplementary Material available on Dryad at http://dx.doi.org/10.5061/dryad.968t3.

Molecular epidemiology of coxsackievirus type B1.

Coxsackievirus type B1 (CVB1) has emerged globally as the predominant enterovirus serotype and is associated with epidemics of meningitis and chronic diseases. In this report, the phylogeny of CVB1 was studied based on the VP1 sequences of 11 North African isolates and 81 published sequences. All CVB1 isolates segregated into four distinct genogroups and 10 genotypes. Most of the identified genotypes of circulating CVB1 strains appear to have a strict geographical specificity. The North African strains were of a single genotype and probably evolved distinctly. Using a relaxed molecular clock model and three different population models (constant population, exponential growth and Bayesian skyline demographic models) in coalescent analysis using the BEAST program, the substitution rate in CVB1 varied between 6.95 × 10(-3) and 7.37 × 10(-3) substitutions/site/year in the VP1 region. This study permits better identification of circulating CVB1, which has become one of the most predominant enterovirus serotypes in humans.

Time scale evolution of avipoxviruses

Avipoxviruses are divided into three clades: canarypox-like viruses, fowlpox-like viruses, and psittacinepox-like viruses. Several molecular clock and demographic models available in the BEAST package were compared on three avipoxvirus genes (P4b, cnpv186 and DNA polymerase genes), which enabled to determine that avipoxviruses evolved at a rate of 2–8×10−5substitution/site/year, in the range of poxviruses previously reported evolution rates. In addition, the date of mean time of divergence of avipoxviruses from a common ancestor was extrapolated to be about 10,000–30,000years ago, at the same period as modern poxvirus species. Our findings will facilitate epidemiological investigations on avipoxviruses’ spread, origin and circulation.

Characterisation of the relationship between Spirometra erinaceieuropaei and Diphyllobothrium species using complete cytb and cox1 genes

Spirometra erinaceieuropaei is an important diphyllobothriid cestode whose larval plerocercoid (sparganum) can parasitise humans, causing sparganosis. In China, sparganosis is an important foodborne parasitic zoonosis and poses a serious threat to human health. However, our knowledge of the taxonomic classification of S. erinaceieuropaei and its relationship with other diphyllobothriid cestodes, such as Diphyllobothrium, remains fragmentary. In this study, we determined the complete sequences of cytochrome b (cytb) and cytochrome c oxidase subunit I (cox1) genes from 61 isolates of S. erinaceieuropaei collected from 11 geographical locations in eastern China and compared them to the sequences of Diphyllobothrium species from the GenBank database. Then, we evaluated the genetic variation of S. erinaceieuropaei and the phylogenetic relationship between S. erinaceieuropaei and Diphyllobothrium. Both negative Tajima’s D and Fu’s FS values of neutrality tests and unimodal curve of mismatch distribution analyses supported demographic population expansion. Bayesian skyline plot (BSP) analyses indicated that a population expansion occurred in the Pleistocene. Phylogenetic analysis based on cytb, cox1, and cytb+cox1 under all three methods tested (neighbour-joining, maximum parsimony, and maximum likelihood) suggested that the S. erinaceieuropaei isolates from eastern China and the genus Diphyllobothrium are separate monophyletic groups. Furthermore, the uncorrelated log-normal relaxed molecular-clock model suggested that divergence between Spirometra and Diphyllobothrium began in the late Miocene.

Evaluating the Adequacy of Molecular Clock Models Using Posterior Predictive Simulations.

Molecular clock models are commonly used to estimate evolutionary rates and timescales from nucleotide sequences. The goal of these models is to account for rate variation among lineages, such that they are assumed to be adequate descriptions of the processes that generated the data. A common approach for selecting a clock model for a data set of interest is to examine a set of candidates and to select the model that provides the best statistical fit. However, this can lead to unreliable estimates if all the candidate models are actually inadequate. For this reason, a method of evaluating absolute model performance is critical. We describe a method that uses posterior predictive simulations to assess the adequacy of clock models. We test the power of this approach using simulated data and find that the method is sensitive to bias in the estimates of branch lengths, which tends to occur when using underparameterized clock models. We also compare the performance of the multinomial test statistic, originally developed to assess the adequacy of substitution models, but find that it has low power in identifying the adequacy of clock models. We illustrate the performance of our method using empirical data sets from coronaviruses, simian immunodeficiency virus, killer whales, and marine turtles. Our results indicate that methods of investigating model adequacy, including the one proposed here, should be routinely used in combination with traditional model selection in evolutionary studies. This will reveal whether a broader range of clock models to be considered in phylogenetic analysis.

Microbial DNA records historical delivery of anthropogenic mercury.

Mercury (Hg) is an anthropogenic pollutant that is toxic to wildlife and humans, but the response of remote ecosystems to globally distributed Hg is elusive. Here, we use DNA extracted from a dated sediment core to infer the response of microbes to historical Hg delivery. We observe a significant association between the mercuric reductase gene (merA) phylogeny and the timing of Hg deposition. Using relaxed molecular clock models, we show a significant increase in the scaled effective population size of the merA gene beginning ~200 years ago, coinciding with the Industrial Revolution and a coincident strong signal for positive selection acting on residues in the terminal region of the mercuric reductase. This rapid evolutionary response of microbes to changes in the delivery of anthropogenic Hg indicates that microbial genomes record ecosystem response to pollutant deposition in remote regions.

Bridging epidemiology with population genetics in a low incidence MSM-driven HIV-1 subtype B epidemic in Central Europe.

BackgroundThe HIV-1 epidemic in Slovenia, a small Central European country, has some characteristics that make it an ideal model to study HIV-1 transmission. The epidemic is predominantly affecting men who have sex with men infected with subtype B (89% of all patients), has a low prevalence (less than 1/1000) and is growing slowly. The aim of the present study was to analyze in detail the evolutionary history and the determinants of transmission.MethodsA total of 223 partial pol gene sequences from therapy naïve individuals were included, representing 52% of all patients newly diagnosed in 13 years (2000–2012) and analyzed together with genetically similar worldwide sequences, selected in a BLAST search.ResultsCombined analysis (maximum likelihood and Bayesian) of HIV-1 transmission chains revealed 8 major clusters (n ≥ 10 patients), 1 group of 4 patients, 2 trios and 12 transmission pairs, thus leaving only 43 (19.3%) Slovenian patients infected with subtype B without a local epidemiological link, indicating a predominance of local transmission of HIV-1 infection. Bayesian analysis performed on a full set of sequences estimated several introductions of HIV-1 into Slovenia, with the most recent common ancestor (tMRCA) of the earliest Slovenian cluster dated to the late 1980s, although tMRCAs obtained from separate independent analysis of each cluster showed considerably more recent estimates. These findings indicate inconsistencies in molecular clock estimation, which we further explored.We hypothesize that these inconsistent dating estimates across the tree could be caused by an evolutionary rate acceleration of HIV-1 after entering the Slovenia epidemic that is not taken into account by the molecular clock model. It could be caused by a lower transmission rate in this setting, as demonstrated by the low epidemic growth rate estimated by Bayesian skyline plot analysis.ConclusionsHIV-1 subtype B was introduced into Slovenia at several time points from the late 80s onward. The majority of patients had a local transmission link, indicating a closed HIV community. The observed slower epidemic rate suggests that individuals with a long-lasting infection are the driving force of the epidemic in this region.

Jerzy Jurka – 1950–2014

The passing of Dr. Jerzy Jurka (Jurek to his friends) represents a heartfelt loss to his many friends in the Mobile DNA field. The whole field of Mobile DNA has lost the scientific efforts of a brilliant colleague and for many of us an esteemed friend. The accompanying biography provides an outstanding outline of his life and career, so I will focus primarily on his impact to the field of Mobile DNA. Jurek’s first publication related to Mobile DNA was the discovery of an Alu element in an alpha globin gene in the late 1980s while working with Temple Smith. Jurek was one of the pioneers of well-trained bioinformatics experts to focus his expertise on mobile elements, with one of the early and important findings of subfamilies in the Alu elements. Through the years, his ability to handle large genomic datasets led him to a number of important discoveries, including a bioinformatic prediction of the target site for SINEs and LINEs, the discovery of several new families and clades of mobile elements, as well as overall distribution, evolution and structure of mobile elements in whole genomes. Among his more interesting new mobile elements were the discovery of helitrons and polintrons, identified from whole genome sequence analyses. In the process of his analyses, Jurek developed several new bioinformatic tools, including MASE, an early program for multiple alignment editing useful for aligning repetitive sequences, and CENSOR, a program for identifying various classes of repetitive DNA sequences in genomic sequences. However, I believe that his real strength in this research was his ability to couple his understanding of bioinformatics with a genuine understanding of biological processes. This led him to make some amazing findings in terms of genome evolution, evolution of mobile elements, and even functional conclusions such as the likelihood of an endonuclease sequence preference leading to the preferred Alu insertion sites. These studies resulted in over 120 publications, making many significant and even groundbreaking discoveries in the field of mobile elements. As outlined in the accompanying biography, Jurek’s career did not always take a traditional track. After training in several laboratories in Europe, he went to the University of Michigan for postdoctoral training. However, it was as a fellow with Temple Smith at the Dana Farber Cancer Institute that he began his love affair with mobile elements. He then briefly joined BIONET in the San Francisco Bay area which was associated with the early days of GENBANK. He was greatly influenced at that time by his interactions with two great minds who originally proposed the molecular clock model that was so critical to Jurek’s future work, Emile Zuckerkandl and Linus Pauling. Jurek served as the Assistant Director of Research for the Linus Pauling Institute for several years. In 1994, Jurek founded the Genetic Information Research Institute (GIRI) in the Bay Area and served as its president until his death. Jurek often stated that the strong support of his wife, Elzbieta, was critical to his ability to work through an independent research institute. Through the Institute, he supported his research on mobile elements and mentored many other outstanding young bioinformatics scientists. However, the Institute also served as a tremendous resource to our entire Mobile DNA community, largely funded by grants from the National Library of Medicine, to database information on mobile elements in a database called REPBASE. This database, and its associated online publication REPBASE Updates, represented the ultimate collection of mobile element sequence data for members of the Mobile DNA community. Critically, this database serves as the reference data for sequence analysis programs like REPEATMASKER to identify mobile elements in genomes. To further facilitate development of our field, over the last decade, Jurek organized and held three outstanding meetings on the “Genomic Impact of Eukaryotic Transposable Elements.” These meetings brought a wide variety of scientists together to discuss the structure, evolution and function of mobile DNA, including workshops for those wanting to utilize bioinformatics tools more effectively. These extremely influential meetings are often referred to as the “Asilomar” meetings by investigators in the field. Dr. Jerzy ‘Jurek’ Jurka was not only an outstanding scientist, he was a tremendous leader for our field. He was an investigator who loved a passionate scientific discussion, and equally loved to discuss politics and the family of which he was so proud. He has left an indelible impression on our field with his seminal contributions. He will be greatly missed.