Double trouble: medical implications of genetic duplication and amplification in bacteria

Streptococcus suis manganese transporter mutant as a live attenuated vaccine: Safety, efficacy, and virulence reversion mechanisms

Streptococcus suis is the leading cause of mortality in piglets and is responsible for severe economic losses in the global pork industry. Severe invasive diseases caused by S. suis include sepsis, meningitis, arthritis, and endocarditis. S. suis disease prevention is hampered by the lack of safe and efficacious vaccines. In this study, we constructed an S. suis live attenuated vaccine candidate lacking the major streptococcal manganese transporter, a known virulence determinant of this organism. The safety and efficacy of this live vaccine were evaluated in swine. Our clinical study results showed that when administered at a dose of 1010 CFU, the vaccine strain was safe and efficacious. However, a lower dose of 109 CFU failed to generate significant immune protection. To investigate if an adjuvant could enhance the efficacy of the vaccine at a lower dose, we spiked the vaccine with a polymeric adjuvant and evaluated its performance. Surprisingly, four pigs receiving the adjuvanted vaccine died during the vaccination phase. Pathology, microbiology, and genetic analyses suggested that the vaccine strain reverted to virulence in these animals. Functional genetic analysis found that the vaccine strain acquired compensatory mutations that upregulated the expression of a secondary manganese transporter, which in turn restored the virulence of the vaccine strain. Our results provide a new understanding of S. suis host adaptation mechanisms and useful information for the design of future live-attenuated vaccines.

Replication stress increases de novo CNVs across the malaria parasite genome

Changes in the copy number of large genomic regions, termed copy number variations (CNVs), contribute to important phenotypes in many organisms. CNVs are readily identified using conventional approaches when present in a large fraction of the cell population. However, CNVs that are present in only a few genomes across a population are often overlooked but important; if beneficial under specific conditions, a de novo CNV that arises in a single genome can expand during selection to create a larger population of cells with novel characteristics. While the reach of single cell methods to study de novo CNVs is increasing, we continue to lack information about CNV dynamics in rapidly evolving microbial populations. Here, we investigated de novo CNVs in the genome of the Plasmodium parasite that causes human malaria. The highly AT-rich P. falciparum genome readily accumulates CNVs that facilitate rapid adaptation to new drugs and host environments. We employed a low-input genomics approach optimized for this unique genome as well as specialized computational tools to evaluate the de novo CNV rate both before and after the application of stress. We observed a significant increase in genomewide de novo CNVs following treatment with a replication inhibitor. These stress-induced de novo CNVs encompassed genes that contribute to various cellular pathways and tended to be altered in clinical parasite genomes. This snapshot of CNV dynamics emphasizes the connection between replication stress, DNA repair, and CNV generation in this important microbial pathogen.

Unraveling the evolutionary dynamics of toxin-antitoxin systems in diverse genetic lineages of Escherichia coli including the high-risk clonal complexes

IMPORTANCE

Large-scale genomic studies of E. coli provide an invaluable opportunity to understand how genomic fine-tuning contributes to the transition of bacterial lifestyle from being commensals to mutualists or pathogens. Within this context, through machine learning-based studies, it appears that TA systems play an important role in the classification of high-risk clonal lineages and could be attributed to their epidemiological success. Due to these profound indications and assumptions, we attempted to provide unique insights into the ordered world of TA systems at the population level by investigating the diversity and evolutionary patterns of TA genes across 19 different STs of E. coli. Further in-depth analysis of ST-specific TA structures and associated genetic coordinates holds the potential to elucidate the functional implications of TA systems in bacterial cell survival and persistence, by and large.

CAMPYLOBACTER HYOINTESTINALIS ISOLATION FROM HOWLER (ALOUATTA CARAYA) AND SPIDER MONKEYS (ATELES FUSCICEPS ROBUSTUS) AT A ZOOLOGIC FACILITY IN CENTRAL ILLINOIS

Campylobacter hyointestinalis was initially isolated from an asymptomatic black howler monkey (Alouatta caraya) in a routine fecal culture examination. Fecal cultures from other individuals in this group and an adjacently housed black-headed spider monkey (Ateles fusciceps robustus) group recovered C. hyointestinalis from all but one of the individuals sampled (1.1 spider monkeys and 2.1 howler monkeys). Concurrently, one spider monkey presented with acute onset severe rectal prolapse and diarrhea. Whole-genome sequencing results of C. hyointestinalis isolates from all individuals were homologous and closely related to Campylobacter hyointestinalis subsp. hyointestinalis TTU_618, a strain typically associated with environmental samples. In addition, two cytolethal distending toxin (CDT) expressing gene clusters, cdt-I and cdt-II, were identified in all isolates. These results suggest C. hyointestinalis is transmissible to both howler monkeys and spider monkeys, though the origin of infection and whether it is transmissible between these species is undetermined.

Gene amplification mutations originate prior to selective stress in Acinetobacter baylyi

The controversial theory of adaptive amplification states gene amplification mutations are induced by selective environments where they are enriched due to the stress caused by growth restriction on unadapted cells. We tested this theory with three independent assays using an Acinetobacter baylyi model system that exclusively selects for cat gene amplification mutants. Our results demonstrate all cat gene amplification mutant colonies arise through a multistep process. While the late steps occur during selection exposure, these mutants derive from low-level amplification mutant cells that form before growth-inhibiting selection is imposed. During selection, these partial mutants undergo multiple secondary steps generating higher amplification over several days to multiple weeks to eventually form visible high-copy amplification colonies. Based on these findings, amplification in this Acinetobacter system can be explained by a natural selection process that does not require a stress response. These findings have fundamental implications to understanding the role of growth-limiting selective environments on cancer development. We suggest duplication mutations encompassing growth factor genes may serve as new genomic biomarkers to facilitate early cancer detection and treatment, before high-copy amplification is attained.

Third Generation Genome Sequencing Reveals That Endobacteria in Nematophagous Fungi Esteya vermicola Contain Multiple Genes Encoding for Nematicidal Proteins

Esteya vermicola is the first recorded endoparasitic nematophagous fungus with high infectivity capacity, attacking the pinewood nematode Bursaphelenchus xylophilus which causes pine wilt disease. Endosymbionts are found in the cytoplasm of E. vermicola from various geographical areas. We sequenced the genome of endobacteria residing in E. vermicola to discover possible biological functions of these widespread endobacteria. Multilocus phylogenetic analyses showed that the endobacteria form a previously unidentified lineage sister to Phyllobacterium myrsinacearum species. The number of genes in the endobacterium was 4542, with 87.8% of the proteins having a known function. It contained a high proportion of repetitive sequences, as well as more Acyl-CoA synthetase genes and genes encoding the electron transport chain, compared with compared with plant-associated P. zundukense Tri 48 and P. myrsinacearum DSM 5893. Thus, this symbiotic bacterium is likely to be more efficient in regulating gene expression and energy release. Furthermore, the endobacteria in nematophagous fungi Esteya vermicola contained multiple nematicidal subtilase/subtilisin encoding genes, so it is likely that endobacteria cooperate with the host to kill nematodes.

Subfunctionalization influences the expansion of bacterial multidrug antibiotic resistance

Background

Antibiotic resistance is a major problem for human health. Multidrug resistance efflux pumps, especially those of the Resistance-Nodulation-Cell Division (RND) family, are major contributors to high-level antibiotic resistance in Gram-negative bacteria. Most bacterial genomes contain several copies of the different classes of multidrug resistance efflux pumps. Gene duplication and gain of function by the duplicate copies of multidrug resistance efflux pump genes plays a key role in the expansion and diversification of drug-resistance mechanisms.

Results

We used two members of the Burkholderia RND superfamily as models to understand how duplication events affect the antibiotic resistance of these strains. First, we analyzed the conservation and distribution of these two RND systems and their regulators across the Burkholderia genus. Through genetic manipulations, we identified both the exact substrate range of these transporters and their eventual interchangeability. We also performed a directed evolution experiment, combined with next generation sequencing, to evaluate the role of antibiotics in the activation of the expression of these systems. Together, our results indicate that the first step to diversify the functions of these pumps arises from changes in their regulation (subfunctionalization) instead of functional mutations. Further, these pumps could rewire their regulation to respond to antibiotics, thus maintaining high genomic plasticity.

Conclusions

Studying the regulatory network that controls the expression of the RND pumps will help understand and eventually control the development and expansion of drug resistance.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4222-4) contains supplementary material, which is available to authorized users.

Copy number change: evolving views on gene amplification

The rapid pace of genomic sequence analysis is increasing the awareness of intrinsically dynamic genetic landscapes. Gene duplication and amplification (GDA) contribute to adaptation and evolution by allowing DNA regions to expand and contract in an accordion-like fashion. This process affects diverse aspects of bacterial infection, including antibiotic resistance and host-pathogen interactions. In this review, microbial GDA is discussed, primarily using recent bacterial examples that demonstrate medical and evolutionary consequences. Interplay between GDA and horizontal gene transfer further impact evolutionary trajectories. Complementing the discovery of gene duplication in clinical and environmental settings, experimental evolution provides a powerful method to document genetic change over time. New methods for GDA detection highlight both its importance and its potential application for genetic engineering, synthetic biology and biotechnology.

Detection of copy number variation and single nucleotide polymorphisms in genes involved in drug resistance and other phenotypic traits in P. falciparum clinical isolates collected from Uganda

There is an increasing interest in mapping the genes of pathogens which underlie important phenotypic traits such as virulence and drug resistance. The Plasmodium falciparum genome exhibits sequence variation that contributes to the pathogenic mechanisms of the parasite. Determining the prevalence of resistance markers could provide a prediction about drug efficacy. Copy number polymorphism (CNP) of genes has been shown to influence important parasite phenotypes. In this work, CNPs within genes involved in drug resistance and other phenotypic traits namely P. falciparum multidrug resistance 1 (pfmdr-1), GTP cyclo hydrolase (gch1), Ring infected erythrocyte surface antigen precursor (resa) and a hypothetical protein coding gene were analyzed by quantitative real time-polymerase reaction (qRT-PCR) among clinical isolates collected from Uganda. The pfmdr-1 codons 86 and 1246 and P. falciparum chloroquine resistance (pfcrt) codon 76 were genotyped for single nucleotide polymorphisms (SNPs) by polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP), and the proportion of resistance associated mutations were determined among mild and severe malaria cases using the chi-square test. Forty and 42 P. falciparum isolates collected from children with mild and severe malaria respectively were analyzed for CNPs. Seventy five and 81 P. falciparum isolates from children with mild or severe malaria were analyzed for SNPs. No pfmdr-1, gch1 or novel gene amplifications were identified among the P. falciparum clinical isolates. Although chloroquine was officially withdrawn from policy use since 7 years, all P. falciparum isolates presented the associated pfcrt K76T mutation, whatever the clinical status and no specific mutation in the pfmdr-1 gene was associated with disease type. In conclusion, this study provides baseline measures for continued surveillance for changes in copy number and SNP types among genes implicated in drug resistance and other important phenotypes that may have a potential role in parasite virulence mechanisms or drug treatment outcomes.

Gene duplication as a mechanism of genomic adaptation to a changing environment

A subject of extensive study in evolutionary theory has been the issue of how neutral, redundant copies can be maintained in the genome for long periods of time. Concurrently, examples of adaptive gene duplications to various environmental conditions in different species have been described. At this point, it is too early to tell whether or not a substantial fraction of gene copies have initially achieved fixation by positive selection for increased dosage. Nevertheless, enough examples have accumulated in the literature that such a possibility should be considered. Here, I review the recent examples of adaptive gene duplications and make an attempt to draw generalizations on what types of genes may be particularly prone to be selected for under certain environmental conditions. The identification of copy-number variation in ecological field studies of species adapting to stressful or novel environmental conditions may improve our understanding of gene duplications as a mechanism of adaptation and its relevance to the long-term persistence of gene duplications.

Analysis of IS1236-mediated gene amplification events in Acinetobacter baylyi ADP1

Recombination between insertion sequence copies can cause genetic deletion, inversion, or duplication. However, it is difficult to assess the fraction of all genomic rearrangements that involve insertion sequences. In previous gene duplication and amplification studies of Acinetobacter baylyi ADP1, an insertion sequence was evident in approximately 2% of the characterized duplication sites. Gene amplification occurs frequently in all organisms and has a significant impact on evolution, adaptation, drug resistance, cancer, and various disorders. To understand the molecular details of this important process, a previously developed system was used to analyze gene amplification in selected mutants. The current study focused on amplification events in two chromosomal regions that are near one of six copies of the only transposable element in ADP1, IS1236 (an IS3 family member). Twenty-one independent mutants were analyzed, and in contrast to previous studies of a different chromosomal region, IS1236 was involved in 86% of these events. IS1236-mediated amplification could occur through homologous recombination between insertion sequences on both sides of a duplicated region. However, this mechanism presupposes that transposition generates an appropriately positioned additional copy of IS1236. To evaluate this possibility, PCR and Southern hybridization were used to determine the chromosomal configurations of amplification mutants involving IS1236. Surprisingly, the genomic patterns were inconsistent with the hypothesis that intramolecular homologous recombination occurred between insertion sequences following an initial transposition event. These results raise a novel possibility that the gene amplification events near the IS1236 elements arise from illegitimate recombination involving transposase-mediated DNA cleavage.

Contribution of gene amplification to evolution of increased antibiotic resistance in Salmonella typhimurium

The use of beta-lactam antibiotics has led to the evolution and global spread of a variety of resistance mechanisms, including beta-lactamases, a group of enzymes that degrade the beta-lactam ring. The evolution of increased beta-lactam resistance was studied by exposing independent lineages of Salmonella typhimurium to progressive increases in cephalosporin concentration. Each lineage carried a beta-lactamase gene (bla(TEM-1)) that provided very low resistance. In most lineages, the initial response to selection was an amplification of the bla(TEM-1) gene copy number. Amplification was followed in some lineages by mutations (envZ, cpxA, or nmpC) that reduced expression of the uptake functions, the OmpC, OmpD, and OmpF porins. The initial resistance provided by bla(TEM-1) amplification allowed the population to expand sufficiently to realize rare secondary point mutations. Mathematical modeling showed that amplification often is likely to be the initial response because events that duplicate or further amplify a gene are much more frequent than point mutations. These models show the importance of the population size to appearance of later point mutations. Transient gene amplification is likely to be a common initial mechanism and an intermediate in stable adaptive improvement. If later point mutations (allowed by amplification) provide sufficient adaptive improvement, the amplification may be lost.

Adaptive copy number evolution in malaria parasites

Copy number polymorphism (CNP) is ubiquitous in eukaryotic genomes, but the degree to which this reflects the action of positive selection is poorly understood. The first gene in the Plasmodium folate biosynthesis pathway, GTP-cyclohydrolase I (gch1), shows extensive CNP. We provide compelling evidence that gch1 CNP is an adaptive consequence of selection by antifolate drugs, which target enzymes downstream in this pathway. (1) We compared gch1 CNP in parasites from Thailand (strong historical antifolate selection) with those from neighboring Laos (weak antifolate selection). Two percent of chromosomes had amplified copy number in Laos, while 72% carried multiple (2–11) copies in Thailand, and differentiation exceeded that observed at 73 synonymous SNPs. (2) We found five amplicon types containing one to greater than six genes and spanning 1 to >11 kb, consistent with parallel evolution and strong selection for this gene amplification. gch1 was the only gene occurring in all amplicons suggesting that this locus is the target of selection. (3) We observed reduced microsatellite variation and increased linkage disequilibrium (LD) in a 900-kb region flanking gch1 in parasites from Thailand, consistent with rapid recent spread of chromosomes carrying multiple copies of gch1. (4) We found that parasites bearing dhfr-164L, which causes high-level resistance to antifolate drugs, carry significantly (p = 0.00003) higher copy numbers of gch1 than parasites bearing 164I, indicating functional association between genes located on different chromosomes but linked in the same biochemical pathway. These results demonstrate that CNP at gch1 is adaptive and the associations with dhfr-164L strongly suggest a compensatory function. More generally, these data demonstrate how selection affects multiple enzymes in a single biochemical pathway, and suggest that investigation of structural variation may provide a fast-track to locating genes underlying adaptation.

Recent comparative genomic hybridization studies have revealed extensive copy number variation in eukaryotic genomes. The first gene in the Plasmodium folate biosynthesis pathway, GTP-cyclohydrolase I (gch1), shows extensive copy number polymorphism (CNP). We provide compelling evidence that gch1 CNP is adaptive and most likely results from selection by antifolate drugs, which target enzymes downstream in this pathway. Gch1 CNP shows extreme geographical differentiation; hitchhiking reduces diversity and increases LD in flanking sequence, indicating recent rapid spread within Thailand, while amplicon structure reveals multiple origins and parallel evolution. Furthermore, strong association between elevated copy number and a critical mutation dhfr-I164L that underlies high-level antifolate resistance indicates functional linkage and fitness epistasis between genes on different chromosomes. These data reveal hidden complexity in the evolutionary response to antifolate treatment and demonstrate that analysis of structural variation can provide a fast-track to locating genes that underlie adaptation.

MEMS and the microbe

In recent years, relatively simple MEMS fabrications have helped accelerate our knowledge of the microbial cell. Current progress and challenges in the application of lab-on-a-chip devices to the viable microbe are reviewed. Furthermore, the degree to which microbiologists are becoming the engineers and are tailoring microbial cells and protocells as potential components for bioMEMS devices is highlighted. We conclude this is a highly productive time for microbiologists and microengineers to unite their shared interest in the micron scale world.

Plastic architecture of bacterial genome revealed by comparative genomics of Photorhabdus variants

BACKGROUND

The phenotypic consequences of large genomic architecture modifications within a clonal bacterial population are rarely evaluated because of the difficulties associated with using molecular approaches in a mixed population. Bacterial variants frequently arise among Photorhabdus luminescens, a nematode-symbiotic and insect-pathogenic bacterium. We therefore studied genome plasticity within Photorhabdus variants.

RESULTS

We used a combination of macrorestriction and DNA microarray experiments to perform a comparative genomic study of different P. luminescens TT01 variants. Prolonged culturing of TT01 strain and a genomic variant, collected from the laboratory-maintained symbiotic nematode, generated bacterial lineages composed of primary and secondary phenotypic variants and colonial variants. The primary phenotypic variants exhibit several characteristics that are absent from the secondary forms. We identify substantial plasticity of the genome architecture of some variants, mediated mainly by deletions in the 'flexible' gene pool of the TT01 reference genome and also by genomic amplification. We show that the primary or secondary phenotypic variant status is independent from global genomic architecture and that the bacterial lineages are genomic lineages. We focused on two unusual genomic changes: a deletion at a new recombination hotspot composed of long approximate repeats; and a 275 kilobase single block duplication belonging to a new class of genomic duplications.

CONCLUSION

Our findings demonstrate that major genomic variations occur in Photorhabdus clonal populations. The phenotypic consequences of these genomic changes are cryptic. This study provides insight into the field of bacterial genome architecture and further elucidates the role played by clonal genomic variation in bacterial genome evolution.

Genomic fluidity and pathogenic bacteria: applications in diagnostics, epidemiology and intervention

The increasing availability of DNA-sequence information for multiple pathogenic and non-pathogenic variants of individual bacterial species has indicated that both DNA acquisition and genome reduction have important roles in genome evolution. Such genomic fluidity, which is found in human pathogens such as Escherichia coli, Helicobacter pylori and Mycobacterium tuberculosis, has important consequences for the clinical management of the diseases that are caused by these pathogens and for the development of diagnostics and new molecular epidemiological methods.

Comparing whole genomes using DNA microarrays

The rapid accumulation of complete genomic sequences offers the opportunity to carry out an analysis of inter- and intra-individual genome variation within a species on a routine basis. Sequencing whole genomes requires resources that are currently beyond those of a single laboratory and therefore it is not a practical approach for resequencing hundreds of individual genomes. DNA microarrays present an alternative way to study differences between closely related genomes. Advances in microarray-based approaches have enabled the main forms of genomic variation (amplifications, deletions, insertions, rearrangements and base-pair changes) to be detected using techniques that are readily performed in individual laboratories using simple experimental approaches.

A process for analysis of microarray comparative genomics hybridisation studies for bacterial genomes

Background

Microarray based comparative genomic hybridisation (CGH) experiments have been used to study numerous biological problems including understanding genome plasticity in pathogenic bacteria. Typically such experiments produce large data sets that are difficult for biologists to handle. Although there are some programmes available for interpretation of bacterial transcriptomics data and CGH microarray data for looking at genetic stability in oncogenes, there are none specifically to understand the mosaic nature of bacterial genomes. Consequently a bottle neck still persists in accurate processing and mathematical analysis of these data. To address this shortfall we have produced a simple and robust CGH microarray data analysis process that may be automated in the future to understand bacterial genomic diversity.

Results

The process involves five steps: cleaning, normalisation, estimating gene presence and absence or divergence, validation, and analysis of data from test against three reference strains simultaneously. Each stage of the process is described and we have compared a number of methods available for characterising bacterial genomic diversity, for calculating the cut-off between gene presence and absence or divergence, and shown that a simple dynamic approach using a kernel density estimator performed better than both established, as well as a more sophisticated mixture modelling technique. We have also shown that current methods commonly used for CGH microarray analysis in tumour and cancer cell lines are not appropriate for analysing our data.

Conclusion

After carrying out the analysis and validation for three sequenced Escherichia coli strains, CGH microarray data from 19 E. coli O157 pathogenic test strains were used to demonstrate the benefits of applying this simple and robust process to CGH microarray studies using bacterial genomes.