Patterns of linkage disequilibrium, homoplasy, and incompatibility are difficult to interpret because they depend on several factors, including the recombination process and the population structure. ...Here we introduce a novel model-based framework to infer recombination properties from such summary statistics in bacterial genomes. The underlying model is sequentially Markovian so that data can be simulated very efficiently, and we use approximate Bayesian computation techniques to infer parameters. As this does not require us to calculate the likelihood function, the model can be easily extended to investigate less probed aspects of recombination. In particular, we extend our model to account for the bias in the recombination process whereby closely related bacteria recombine more often with one another. We show that this model provides a good fit to a data set of Bacillus cereus genomes and estimate several recombination properties, including the rate of bias in recombination. All the methods described in this article are implemented in a software package that is freely available for download at http://code.google.com/p/clonalorigin/.
High genetic diversity is a hallmark of the gastric pathogen Helicobacter pylori. We used 454 sequencing technology to perform whole-genome comparisons for five sets of H. pylori strains that had ...been sequentially cultured from four chronically infected Colombians (isolation intervals = 3-16 y) and one human volunteer experimentally infected with H. pylori as part of a vaccine trial. The four sets of genomes from Colombian H. pylori differed by 27-232 isolated SNPs and 16-441 imported clusters of polymorphisms resulting from recombination. Imports (mean length = 394 bp) were distributed nonrandomly over the chromosome and frequently occurred in groups, suggesting that H. pylori first takes up long DNA fragments, which subsequently become partially integrated in multiple shorter pieces. Imports were present at significantly increased frequency in members of the hop family of outer membrane gene paralogues, some of which are involved in bacterial adhesion, suggesting diversifying selection. No evidence of recombination and few other differences were identified in the strain pair from an infected volunteer, indicating that the H. pylori genome is stable in the absence of mixed infection. Among these few differences was an OFF/ON switch in the phase-variable adhesin gene hopZ, suggesting strong in vivo selection for this putative adhesin during early colonization.
Full text
Available for:
BFBNIB, NMLJ, NUK, PNG, SAZU, UL, UM, UPUK
The tuberculosis (TB) epidemic is fueled by a parallel Human Immunodeficiency Virus (HIV) epidemic, but it remains unclear to what extent the HIV epidemic has been a driver for drug resistance in ...Mycobacterium tuberculosis (Mtb). Here we assess the impact of HIV co-infection on the emergence of resistance and transmission of Mtb in the largest outbreak of multidrug-resistant TB in South America to date. By combining Bayesian evolutionary analyses and the reconstruction of transmission networks utilizing a new model optimized for TB, we find that HIV co-infection does not significantly affect the transmissibility or the mutation rate of Mtb within patients and was not associated with increased emergence of resistance within patients. Our results indicate that the HIV epidemic serves as an amplifier of TB outbreaks by providing a reservoir of susceptible hosts, but that HIV co-infection is not a direct driver for the emergence and transmission of resistant strains.
Epidemics and pandemics of cholera, a severe diarrheal disease, have occurred since the early 19th century and waves of epidemic disease continue today. Cholera epidemics are caused by individual, ...genetically monomorphic lineages of Vibrio cholerae: the ongoing seventh pandemic, which has spread globally since 1961, is associated with lineage L2 of biotype El Tor. Previous genomic studies of the epidemiology of the seventh pandemic identified three successive sub-lineages within L2, designated waves 1 to 3, which spread globally from the Bay of Bengal on multiple occasions. However, these studies did not include samples from China, which also experienced multiple epidemics of cholera in recent decades. We sequenced the genomes of 71 strains isolated in China between 1961 and 2010, as well as eight from other sources, and compared them with 181 published genomes. The results indicated that outbreaks in China between 1960 and 1990 were associated with wave 1 whereas later outbreaks were associated with wave 2. However, the previously defined waves overlapped temporally, and are an inadequate representation of the shape of the global genealogy. We therefore suggest replacing them by a series of tightly delineated clades. Between 1960 and 1990 multiple such clades were imported into China, underwent further microevolution there and then spread to other countries. China was thus both a sink and source during the pandemic spread of V. cholerae, and needs to be included in reconstructions of the global patterns of spread of cholera.
Full text
Available for:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
Many of the most virulent bacterial pathogens show low genetic diversity and sexual isolation. Accordingly, Mycobacterium tuberculosis, the deadliest human pathogen, is thought to be clonal and ...evolve by genetic drift. Yet, its genome shows few of the concomitant signs of genome degradation. We analyzed 24 genomes and found an excess of genetic diversity in regions encoding key adaptive functions including the type VII secretion system and the ancient horizontally transferred virulence-related regions. Four different approaches showed evident signs of recombination in M. tuberculosis. Recombination tracts add a high density of polymorphisms, and many are thus predicted to arise from outside the clade. Some of these tracts match Mycobacterium canettii sequences. Recombination introduced an excess of non-synonymous diversity in general and even more in genes expected to be under positive or diversifying selection, e.g., cell wall component genes. Mutations leading to non-synonymous SNPs are effectively purged in MTBC, which shows dominance of purifying selection. MTBC mutation bias toward AT nucleotides is not compensated by biased gene conversion, suggesting the action of natural selection also on synonymous changes. Together, all of these observations point to a strong imprint of recombination and selection in the genome affecting both non-synonymous and synonymous positions. Hence, contrary to some other pathogens and previous proposals concerning M. tuberculosis, this lineage may have come out of its ancestral bottleneck as a very successful pathogen that is rapidly diversifying by the action of mutation, recombination, and natural selection.
miRNAs play important roles in the regulation of gene functions. Maternal dietary modifications during pregnancy and gestation have long-term effects on the offspring, but it is not known whether a ...maternal high fat (HF) diet during pregnancy and lactation alters expression of key miRNAs in the offspring.
We studied the effects of maternal HF diet on the adult offspring by feeding mice with either a HF or a chow diet prior to conception, during pregnancy and lactation, and all offspring were weaned onto the same chow diet until adulthood. Maternal HF fed offspring had markedly increased hepatic mRNA levels of peroxisome proliferator activated receptor-alpha (ppar-alpha) and carnitine palmitoyl transferase-1a (cpt-1a) as well as insulin like growth factor-2 (Igf2). A HF diet induced up-regulation of ppar-alpha and cpt-1a expression in the wild type but not in Igf2 knock out mice. Furthermore, hepatic expression of let-7c was also reduced in maternal HF fed offspring. Among 579 miRNAs measured with microarray, ~23 miRNA levels were reduced by ~1.5-4.9-fold. Reduced expression of miR-709 (a highly expressed miRNA), miR-122, miR-192, miR-194, miR-26a, let-7a, let7b and let-7c, miR-494 and miR-483* (reduced by ~4.9 fold) was validated by qPCR. We found that methyl-CpG binding protein 2 was the common predicted target for miR-709, miR-let7s, miR-122, miR-194 and miR-26a using our own purpose-built computer program.
Maternal HF feeding during pregnancy and lactation induced co-ordinated and long-lasting changes in expression of Igf2, fat metabolic genes and several important miRNAs in the offspring.
Full text
Available for:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
Human networks of sexual contacts are dynamic by nature, with partnerships forming and breaking continuously over time. Sexual behaviours are also highly heterogeneous, so that the number of partners ...reported by individuals over a given period of time is typically distributed as a power-law. Both the dynamism and heterogeneity of sexual partnerships are likely to have an effect in the patterns of spread of sexually transmitted diseases. To represent these two fundamental properties of sexual networks, we developed a stochastic process of dynamic partnership formation and dissolution, which results in power-law numbers of partners over time. Model parameters can be set to produce realistic conditions in terms of the exponent of the power-law distribution, of the number of individuals without relationships and of the average duration of relationships. Using an outbreak of antibiotic resistant gonorrhoea amongst men have sex with men as a case study, we show that our realistic dynamic network exhibits different properties compared to the frequently used static networks or homogeneous mixing models. We also consider an approximation to our dynamic network model in terms of a much simpler branching process. We estimate the parameters of the generation time distribution and offspring distribution which can be used for example in the context of outbreak reconstruction based on genomic data. Finally, we investigate the impact of a range of interventions against gonorrhoea, including increased condom use, more frequent screening and immunisation, concluding that the latter shows great promise to reduce the burden of gonorrhoea, even if the vaccine was only partially effective or applied to only a random subset of the population.
Full text
Available for:
DOBA, IZUM, KILJ, NUK, PILJ, PNG, SAZU, SIK, UILJ, UKNU, UL, UM, UPUK
The development of high-throughput sequencing technology has led to a significant reduction in the time and cost of sequencing whole genomes of bacterial pathogens. Studies can sequence and compare ...hundreds or even thousands of genomes within a given bacterial population. A phylogenetic tree is the most frequently used method of depicting the relationships between these bacterial pathogen genomes. However, the presence of homologous recombination in most bacterial pathogen species can invalidate the application of standard phylogenetic tools. Here we describe a method to produce phylogenetic analyses that accounts for the disruptive effect of recombination. This allows users to investigate the recombination events that have occurred, as well as to produce more meaningful phylogenetic analyses which recover the clonal genealogy representing the clonal relationships between genomes.
Recombination plays a dominant role in the evolution of the bacterial pathogen Helicobacter pylori, but its dynamics remain incompletely understood. Here we use an in vitro transformation system ...combined with genome sequencing to study chromosomal integration patterns after natural transformation. A single transformation cycle results in up to 21 imports, and repeated transformations generate a maximum of 92 imports (8% sequence replacement). Import lengths show a bimodal distribution with averages of 28 and 1,645 bp. Reanalysis of paired H. pylori genomes from chronically infected people demonstrates the same bimodal import pattern in vivo. Restriction endonucleases (REases) of the recipient bacteria fail to inhibit integration of homeologous DNA, independently of methylation. In contrast, REases limit the import of heterologous DNA. We conclude that restriction-modification systems inhibit the genomic integration of novel sequences, while they pose no barrier to homeologous recombination, which reconciles the observed stability of the H. pylori gene content and its highly recombinational population structure.
Corynebacterium diphtheriae, the agent of diphtheria, is a genetically diverse bacterial species. Although antimicrobial resistance has emerged against several drugs including first-line penicillin, ...the genomic determinants and population dynamics of resistance are largely unknown for this neglected human pathogen.
Here, we analyzed the associations of antimicrobial susceptibility phenotypes, diphtheria toxin production, and genomic features in C. diphtheriae. We used 247 strains collected over several decades in multiple world regions, including the 163 clinical isolates collected prospectively from 2008 to 2017 in France mainland and overseas territories.
Phylogenetic analysis revealed multiple deep-branching sublineages, grouped into a Mitis lineage strongly associated with diphtheria toxin production and a largely toxin gene-negative Gravis lineage with few toxin-producing isolates including the 1990s ex-Soviet Union outbreak strain. The distribution of susceptibility phenotypes allowed proposing ecological cutoffs for most of the 19 agents tested, thereby defining acquired antimicrobial resistance. Penicillin resistance was found in 17.2% of prospective isolates. Seventeen (10.4%) prospective isolates were multidrug-resistant (≥ 3 antimicrobial categories), including four isolates resistant to penicillin and macrolides. Homologous recombination was frequent (r/m = 5), and horizontal gene transfer contributed to the emergence of antimicrobial resistance in multiple sublineages. Genome-wide association mapping uncovered genetic factors of resistance, including an accessory penicillin-binding protein (PBP2m) located in diverse genomic contexts. Gene pbp2m is widespread in other Corynebacterium species, and its expression in C. glutamicum demonstrated its effect against several beta-lactams. A novel 73-kb C. diphtheriae multiresistance plasmid was discovered.
This work uncovers the dynamics of antimicrobial resistance in C. diphtheriae in the context of phylogenetic structure, biovar, and diphtheria toxin production and provides a blueprint to analyze re-emerging diphtheria.
Full text
Available for:
IZUM, KILJ, NUK, PILJ, PNG, SAZU, UL, UM, UPUK