Predicting peak spectral sensitivities of vertebrate cone visual pigments using atomistic molecular simulations.

Vision is the dominant sensory modality in many organisms for foraging, predator avoidance, and social behaviors including mate selection. Vertebrate visual perception is initiated when light strikes rod and cone photoreceptors within the neural retina of the eye. Sensitivity to individual colors, i.e., peak spectral sensitivities (λmax) of visual pigments, are a function of the type of chromophore and the amino acid sequence of the associated opsin protein in the photoreceptors. Large differences in peak spectral sensitivities can result from minor differences in amino acid sequence of cone opsins. To determine how minor sequence differences could result in large spectral shifts we selected a spectrally-diverse group of 14 teleost Rh2 cone opsins for which sequences and λmax are experimentally known. Classical molecular dynamics simulations were carried out after embedding chromophore-associated homology structures within explicit bilayers and water. These simulations revealed structural features of visual pigments, particularly within the chromophore, that contributed to diverged spectral sensitivities. Statistical tests performed on all the observed structural parameters associated with the chromophore revealed that a two-term, first-order regression model was sufficient to accurately predict λmax over a range of 452-528 nm. The approach was accurate, efficient and simple in that site-by-site molecular modifications or complex quantum mechanics models were not required to predict λmax. These studies identify structural features associated with the chromophore that may explain diverged spectral sensitivities, and provide a platform for future, functionally predictive opsin modeling.

Comparative genomics of Bifidobacterium species isolated from marmosets and humans.

The genus Bifidobacterium is purported to have beneficial consequences for human health and is a major component of many gastrointestinal probiotics. Although species of Bifidobacterium are generally at low relative frequency in the adult human gastrointestinal tract, they can constitute high proportions of the gastrointestinal communities of adult marmosets. To identify genes that might be important for the maintenance of Bifidobacterium in adult marmosets, ten strains of Bifidobacterium were isolated from the feces of seven adult marmosets, and their genomes were sequenced. There were six B. reuteri strains, two B. callitrichos strains, one B. myosotis sp. nov. and one B. tissieri sp. nov. among our isolates. Phylogenetic analysis showed that three of the four species we isolated were most closely related to B. bifidum, B. breve and B. longum, which are species found in high abundance in human infants. There were 1357 genes that were shared by at least one strain of B. reuteri, B. callitrichos, B. breve, and B. longum, and 987 genes that were found in all strains of the four species. There were 106 genes found in B. reuteri and B. callitrichos but not in human bifidobacteria, and several of these genes were involved in nutrient uptake. These pathways for nutrient uptake appeared to be specific to Bifidobacterium from New World monkeys. Additionally, the distribution of Bifidobacterium in fecal samples from captive adult marmosets constituted as much as 80% of the gut microbiome, although this was variable between individuals and colonies. We suggest that nutrient transporters may be important for the maintenance of Bifidobacterium during adulthood in marmosets.

Genetic diversity of potato virus Y (PVY): sequence analyses reveal ten novel PVY recombinant structures.

The potato virus Y (PVY) complex includes five non-recombinant strains and a growing number of recombinants. Up to now, the bulk of these PVY recombinants were found to be composed of genome segments coming from only two "parental" genomes, PVYO and PVYEu-N, with a small minority including segments from other parents, e.g. PVYC, PVYNA-N, and a segment of the recombinant strain PVY-NE11. Here, comprehensive analyses of 396 whole genomes of PVY isolates from 34 countries and a variety of hosts revealed 28 isolates to be rare or novel recombinants. When subjected to a thorough recombination analysis, these 28 PVY isolates were found to represent 25 poorly sampled PVY recombinant structures, ten of which had not been recognized previously. Nine of the ten novel structures carried parental sequences from PVYNA-N or PVY-NE11 strains recombined with PVYO and PVYEu-N sequences, while seven of the new structures contained sequences from three different parents. The number of known PVY recombinant patterns now stands at thirty six. These recombinant structures present a challenge for PVY strain typing, but may shed light on interactions between PVY and host resistance genes.

Annotation of plasmid genes.

Good annotation of plasmid genomes is essential to maximise the value of the rapidly increasing volume of plasmid sequences. This short review highlights some of the current issues and suggests some ways forward. Where a well-studied related plasmid system exists we recommend that new annotation adheres to the convention already established for that system, so long as it is based on sound principles and solid experimental evidence, even if some of the new genes are more similar to homologues in different systems. Where a well-established model does not exist we provide generic gene names that reflect likely biochemical activity rather than overall purpose particularly, for example, where genes clearly belong to a type IV secretion system but it is not known whether they function in conjugative transfer or virulence. We also recommend that annotators use a whole system naming approach to avoid ending up with an illogical mixture of names from other systems based on the highest scoring match from a BLAST search. In addition, where function has not been experimentally established we recommend using just the locus tag, rather than a function-related gene name, while recording possible functions as notes rather than in a provisional name.

Intrinsically disordered regions of p53 family are highly diversified in evolution.

Proteins of the p53 family are expressed in vertebrates and in some invertebrate species. The main function of these proteins is to control and regulate cell cycle in response to various cellular signals, and therefore to control the organism's development. The regulatory functions of the p53 family members originate mostly from their highly-conserved and well-structured DNA-binding domains. Many human diseases (including various types of cancer) are related to the missense mutations within this domain. The ordered DNA-binding domains of the p53 family members are surrounded by functionally important intrinsically disordered regions. In this study, substitution rates and propensities in different regions of p53 were analyzed. The analyses revealed that the ordered DNA-binding domain is conserved, whereas disordered regions are characterized by high sequence diversity. This diversity was reflected both in the number of substitutions and in the types of substitutions to which each amino acid was prone. These results support the existence of a positive correlation between protein intrinsic disorder and sequence divergence during the evolutionary process. This higher sequence divergence provides strong support for the existence of disordered regions in p53 in vivo for if they were structured, they would evolve at similar rates as the rest of the protein.

Inferring the evolutionary history of IncP-1 plasmids despite incongruence among backbone gene trees.

Plasmids of the incompatibility group IncP-1 can transfer and replicate in many genera of the Proteobacteria. They are composed of backbone genes that encode a variety of essential functions and accessory genes that have implications for human health and environmental remediation. Although it is well understood that the accessory genes are transferred horizontally between plasmids, recent studies have also provided examples of recombination in the backbone genes of IncP-1 plasmids. As a consequence, phylogeny estimation based on backbone genes is expected to produce conflicting gene tree topologies. The main goal of this study was therefore to infer the evolutionary history of IncP-1 plasmids in the presence of both vertical and horizontal gene transfer. This was achieved by quantifying the incongruence among gene trees and attributing it to known causes such as 1) phylogenetic uncertainty, 2) coalescent stochasticity, and 3) horizontal inheritance. Topologies of gene trees exhibited more incongruence than could be attributed to phylogenetic uncertainty alone. Species-tree estimation using a Bayesian framework that takes coalescent stochasticity into account was well supported, but it differed slightly from the maximum-likelihood tree estimated by concatenation of backbone genes. After removal of the gene that demonstrated a signal of intergroup recombination, the concatenated tree was congruent with the species-tree estimate, which itself was robust to inclusion/exclusion of the recombinant gene. Thus, in spite of horizontal gene exchange both within and among IncP-1 subgroups, the backbone genome of these IncP-1 plasmids retains a detectable vertical evolutionary history.

HCMV-infected cells maintain efficient nucleotide excision repair of the viral genome while abrogating repair of the host genome.

Many viruses subvert the host cell's ability to mount and complete various DNA damage responses (DDRs) after infection. HCMV infection of permissive fibroblasts activates host DDRs at the time of viral deposition and during replication, but the DDRs remain uncompleted without arrest or apoptosis. We believe this was in part due to partitioning of the damage response and double strand break repair components. After extraction of soluble proteins, the localization of these components fell into three groups: specifically associated with the viral replication centers (RCs), diffused throughout the nucleoplasm and excluded from the RCs. Others have shown that cells are incapable of processing exogenously introduced damage after infection. We hypothesized that the inability of the cells to process damage might be due to the differential association of repair components within the RCs and, in turn, potentially preferential repair of the viral genome and compromised repair of the host genome. To test this hypothesis we used multiple strategies to examine repair of UV-induced DNA damage in mock and virus-infected fibroblasts. Comet assays indicated that repair was initiated, but was not completed in infected cells. Quantitative analysis of immunofluorescent localization of cyclobutane pyrimidine dimers (CPDs) revealed that after 24 h of repair, CPDs were significantly reduced in viral DNA, but not significantly changed in the infected host DNA. To further quantitate CPD repair, we developed a novel dual-color Southern protocol allowing visualization of host and viral DNA simultaneously. Combining this Southern methodology with a CPD-specific T4 endonuclease V alkaline agarose assay to quantitate repair of adducts, we found efficient repair of CPDs from the viral DNA but not host cellular DNA. Our data confirm that NER functions in HCMV-infected cells and almost exclusively repairs the viral genome to the detriment of the host's genome.

Adaptive regulatory substitutions affect multiple stages in the life cycle of the bacteriophage φX174.

BACKGROUND: Previously, we showed that adaptive substitutions in one of the three promoters of the bacteriophage φX174 improved fitness at high-temperature by decreasing transcript levels three- to four-fold. To understand how such an extreme change in gene expression might lead to an almost two-fold increase in fitness at the adaptive temperature, we focused on stages in the life cycle of the phage that occur before and after the initiation of transcription. For both the ancestral strain and two single-substitution strains with down-regulated transcription, we measured seven phenotypic components of fitness (attachment, ejection, eclipse, virion assembly, latent period, lysis rate and burst size) during a single cycle of infection at each of two temperatures. The lower temperature, 37°C, is the optimal temperature at which phages are cultivated in the lab; the higher temperature, 42°C, exerts strong selection and is the condition under which these substitutions arose in evolution experiments. We augmented this study by developing an individual-based stochastic model of this same life cycle to explore potential explanations for our empirical results. RESULTS: Of the seven fitness parameters, three showed significant differences between strains that carried an adaptive substitution and the ancestor, indicating the presence of pleiotropy in regulatory evolution. 1) Eclipse was longer in the adaptive strains at both the optimal and high-temperature environments. 2) Lysis rate was greater in the adaptive strains at the high temperature. 3) Burst size for the mutants was double that of the ancestor at the high temperature, but half that at the lower temperature. Simulation results suggest that eclipse length and latent period variance can explain differences in burst sizes and fitness between the mutant and ancestral strains. CONCLUSIONS: Down-regulating transcription affects several steps in the phage life cycle, and all of these occur after the initiation of transcription. We attribute the apparent tradeoff between delayed progeny production and faster progeny release to improved host resource utilization at high temperature.

Host range diversification within the IncP-1 plasmid group.

Broad-host-range plasmids play a critical role in the spread of antibiotic resistance and other traits. In spite of increasing information about the genomic diversity of closely related plasmids, the relationship between sequence divergence and host range remains unclear. IncP-1 plasmids are currently classified into six subgroups based on the genetic distance of backbone genes. We investigated whether plasmids from two subgroups exhibit a different host range, using two IncP-1γ plasmids, an IncP-1ß plasmid and their minireplicons. Efficiencies of plasmid establishment and maintenance were compared using five species that belong to the Alphaproteobacteria, Betaproteobacteria and Gammaproteobacteria. The IncP-1ß plasmid replicated and persisted in all five hosts in the absence of selection. Of the two IncP-1γ plasmids, both were unable to replicate in alphaproteobacterial host Sphingobium japonicum, and one established itself in Agrobacterium tumefaciens but was very unstable. In contrast, both IncP-1γ minireplicons, which produced higher levels of replication initiation protein than the wild-type plasmids, replicated in all strains, suggesting that poor establishment of the native plasmids is in part due to suboptimal replication initiation gene regulation. The findings suggest that host ranges of distinct IncP-1 plasmids only partially overlap, which may limit plasmid recombination and thus result in further genome divergence.

Diverse broad-host-range plasmids from freshwater carry few accessory genes.

Broad-host-range self-transferable plasmids are known to facilitate bacterial adaptation by spreading genes between phylogenetically distinct hosts. These plasmids typically have a conserved backbone region and a variable accessory region that encodes host-beneficial traits. We do not know, however, how well plasmids that do not encode accessory functions can survive in nature. The goal of this study was to characterize the backbone and accessory gene content of plasmids that were captured from freshwater sources without selecting for a particular phenotype or cultivating their host. To do this, triparental matings were used such that the only required phenotype was the plasmid's ability to mobilize a nonconjugative plasmid. Based on complete genome sequences of 10 plasmids, only 5 carried identifiable accessory gene regions, and none carried antibiotic resistance genes. The plasmids belong to four known incompatibility groups (IncN, IncP-1, IncU, and IncW) and two potentially new groups. Eight of the plasmids were shown to have a broad host range, being able to transfer into alpha-, beta-, and gammaproteobacteria. Because of the absence of antibiotic resistance genes, we resampled one of the sites and compared the proportion of captured plasmids that conferred antibiotic resistance to their hosts with the proportion of such plasmids captured from the effluent of a local wastewater treatment plant. Few of the captured plasmids from either site encoded antibiotic resistance. A high diversity of plasmids that encode no or unknown accessory functions is thus readily found in freshwater habitats. The question remains how the plasmids persist in these microbial communities.

Comparative genomics of IncP-1ε plasmids from water environments reveals diverse and unique accessory genetic elements.

The goal of this study was to determine and compare the complete genome sequences of three new broad-host-range conjugative plasmids. Plasmids pMLUA1, pMLUA3 and pMLUA4 were previously recovered from estuarine water by exogenous plasmid isolation and ranged in size from ∼55 to 59 kb. Comparative genomics showed that their backbone region was identical to the prototype pKJK5 and other IncP1-ε plasmids captured from soils. The accessory region was inserted between the tra region and parA, and presented the typical IncP-1ε ISPa17 and Tn402-like transposon modules. Nevertheless, new class 1 integrons were identified (In794, carrying aadA5 and In795, carrying qacF5-aadA5), as well as a composite transposon IS26-msr(E)-mph(E)-IS26 carrying genes that confer resistance to macrolides. A new insertion sequence, termed ISUnCu17, was also identified on pMLUA3. The architecture of the accessory regions implies the occurrence of multiple insertions and deletions. These data support the notion that IncP-1 plasmids from the ε subgroup are proficient in the capture of diverse genetic elements, including antibiotic resistance genes, and thus may contribute to the co-selection of several resistance determinants. This study constitutes the first report of completely sequenced IncP-1ε plasmids from water environments, and enhances our understanding of the geographic distribution and genetic diversity of these replicons.

The complete genome sequences of four new IncN plasmids from wastewater treatment plant effluent provide new insights into IncN plasmid diversity and evolution.

The dissemination of antibiotic resistance genes among bacteria often occurs by means of plasmids. Wastewater treatment plants (WWTP) were previously recognized as hot spots for the horizontal transfer of genetic material. One of the plasmid groups that is often associated with drug resistance is the incompatibility group IncN. The aim of this study was to gain insights into the diversity and evolutionary history of IncN plasmids by determining and comparing the complete genome sequences of the four novel multi-drug resistance plasmids pRSB201, pRSB203, pRSB205 and pRSB206 that were exogenously isolated from the final effluent of a municipal WWTP. Their sizes range between 42,875 bp and 56,488 bp and they share a common set of backbone modules that encode plasmid replication initiation, conjugative transfer, and plasmid maintenance and control. All plasmids are transferable at high rates between Escherichia coli strains, but did not show a broad host range. Different genes conferring resistances to ampicillin, streptomycin, spectinomycin, sulfonamides, tetracycline and trimethoprim were identified in accessory modules inserted in these plasmids. Comparative analysis of the four WWTP IncN plasmids and IncN plasmids deposited in the NCBI database enabled the definition of a core set of backbone genes for this group. Moreover, this approach revealed a close phylogenetic relationship between the IncN plasmids isolated from environmental and clinical samples. Phylogenetic analysis also suggests the existence of host-specific IncN plasmid subgroups. In conclusion, IncN plasmids likely contribute to the dissemination of resistance determinants between environmental bacteria and clinical strains. This is of particular importance since multi-drug resistance IncN plasmids have been previously identified in members of the Enterobacteriaceae that cause severe infections in humans.

Comparing models of evolution for ordered and disordered proteins.

Most models of protein evolution are based upon proteins that form relatively rigid 3D structures. A significant fraction of proteins, the so-called disordered proteins, do not form rigid 3D structures and sample a broad conformational ensemble. Disordered proteins do not typically maintain long-range interactions, so the constraints on their evolution should be different than ordered proteins. To test this hypothesis, we developed and compared models of evolution for disordered and ordered proteins. Substitution matrices were constructed using the sequences of putative homologs for sets of experimentally characterized disordered and ordered proteins. Separate matrices, at three levels of sequence similarity (>85%, 85-60%, and 60-40%), were inferred for each type of protein structure. The substitution matrices for disordered and ordered proteins differed significantly at each level of sequence similarity. The disordered matrices reflected a greater likelihood of evolutionary changes, relative to the ordered matrices, and these changes involved nonconservative substitutions. Glutamic acid and asparagine were interesting exceptions to this result. Important differences between the substitutions that are accepted in disordered proteins relative to ordered proteins were also identified. In general, disordered proteins have fewer evolutionary constraints than ordered proteins. However, some residues like tryptophan and tyrosine are highly conserved in disordered proteins. This is due to their important role in forming protein-protein interfaces. Finally, the amino acid frequencies for disordered proteins, computed during the development of the matrices, were compared with amino acid frequencies for different categories of secondary structure in ordered proteins. The highest correlations were observed between the amino acid frequencies in disordered proteins and the solvent-exposed loops and turns of ordered proteins, supporting an emerging structural model for disordered proteins.

Broad-host-range plasmids from agricultural soils have IncP-1 backbones with diverse accessory genes.

Broad-host-range plasmids are known to spread genes between distinct phylogenetic groups of bacteria. These genes often code for resistances to antibiotics and heavy metals or degradation of pollutants. Although some broad-host-range plasmids have been extensively studied, their evolutionary history and genetic diversity remain largely unknown. The goal of this study was to analyze and compare the genomes of 12 broad-host-range plasmids that were previously isolated from Norwegian soils by exogenous plasmid isolation and that encode mercury resistance. Complete nucleotide sequencing followed by phylogenetic analyses based on the relaxase gene traI showed that all the plasmids belong to one of two subgroups (ß and ε) of the well-studied incompatibility group IncP-1. A diverse array of accessory genes was found to be involved in resistance to antimicrobials (streptomycin, spectinomycin, and sulfonamides), degradation of herbicides (2,4-dichlorophenoxyacetic acid and 2,4-dichlorophenoxypropionic acid), and a putative new catabolic pathway. Intramolecular transposition of insertion sequences followed by deletion was found to contribute to the diversity of some of these plasmids. The previous observation that the insertion sites of a Tn501-related element are identical in four IncP-1ß plasmids (pJP4, pB10, R906, and R772) was further extended to three more IncP-1ß plasmids (pAKD15, pAKD18, and pAKD29). We proposed a hypothesis for the evolution of these Tn501-bearing IncP-1ß plasmids that predicts recent diversification followed by worldwide spread. Our study increases the available collection of complete IncP-1 plasmid genome sequences by 50% and will aid future studies to enhance our understanding of the evolution and function of this important plasmid family.

Genetic diversity of the ordinary strain of Potato virus Y (PVY) and origin of recombinant PVY strains.

The ordinary strain of Potato virus Y (PVY), PVY(O), causes mild mosaic in tobacco and induces necrosis and severe stunting in potato cultivars carrying the Ny gene. A novel substrain of PVY(O) was recently reported, PVY(O)-O5, which is spreading in the United States and is distinguished from other PVY(O) isolates serologically (i.e., reacting to the otherwise PVY(N)-specific monoclonal antibody 1F5). To characterize this new PVY(O)-O5 subgroup and address possible reasons for its continued spread, we conducted a molecular study of PVY(O) and PVY(O)-O5 isolates from a North American collection of PVY through whole-genome sequencing and phylogenetic analysis. In all, 44 PVY(O) isolates were sequenced, including 31 from the previously defined PVY(O)-O5 group, and subjected to whole-genome analysis. PVY(O)-O5 isolates formed a separate lineage within the PVY(O) genome cluster in the whole-genome phylogenetic tree and represented a novel evolutionary lineage of PVY from potato. On the other hand, the PVY(O) sequences separated into at least two distinct lineages on the whole-genome phylogenetic tree. To shed light on the origin of the three most common PVY recombinants, a more detailed phylogenetic analysis of a sequence fragment, nucleotides 2,406 to 5,821, that is present in all recombinant and nonrecombinant PVY(O) genomes was conducted. The analysis revealed that PVY(N:O) and PVY(N-Wi) recombinants acquired their PVY(O) segments from two separate PVY(O) lineages, whereas the PVY(NTN) recombinant acquired its PVY(O) segment from the same lineage as PVY(N:O). These data suggest that PVY(N:O) and PVY(N-Wi) recombinants originated from two separate recombination events involving two different PVY(O) parental genomes, whereas the PVY(NTN) recombinants likely originated from the PVY(N:O) genome via additional recombination events.

Biosolids as a Source of Antibiotic Resistance Plasmids for Commensal and Pathogenic Bacteria.

Antibiotic resistance (AR) is a threat to modern medicine, and plasmids are driving the global spread of AR by horizontal gene transfer across microbiomes and environments. Determining the mobile resistome responsible for this spread of AR among environments is essential in our efforts to attenuate the current crisis. Biosolids are a wastewater treatment plant (WWTP) byproduct used globally as fertilizer in agriculture. Here, we investigated the mobile resistome of biosolids that are used as fertilizer. This was done by capturing resistance plasmids that can transfer to human pathogens and commensal bacteria. We used a higher-throughput version of the exogenous plasmid isolation approach by mixing several ESKAPE pathogens and a commensal Escherichia coli with biosolids and screening for newly acquired resistance to about 10 antibiotics in these strains. Six unique resistance plasmids transferred to Salmonella typhimurium, Klebsiella aerogenes, and E. coli. All the plasmids were self-transferable and carried 3-6 antibiotic resistance genes (ARG) conferring resistance to 2-4 antibiotic classes. These plasmids-borne resistance genes were further embedded in genetic elements promoting intracellular recombination (i.e., transposons or class 1 integrons). The plasmids belonged to the broad-host-range plasmid (BHR) groups IncP-1 or PromA. Several of them were persistent in their new hosts when grown in the absence of antibiotics, suggesting that the newly acquired drug resistance traits would be sustained over time. This study highlights the role of BHRs in the spread of ARG between environmental bacteria and human pathogens and commensals, where they may persist. The work further emphasizes biosolids as potential vehicles of highly mobile plasmid-borne antibiotic resistance genes.

Predicting plasmid promiscuity based on genomic signature.

Despite the important contribution of self-transmissible plasmids to bacterial evolution, little is understood about the range of hosts in which these plasmids have evolved. Our goal was to infer this so-called evolutionary host range. The nucleotide composition, or genomic signature, of plasmids is often similar to that of the chromosome of their current host, suggesting that plasmids acquire their hosts' signature over time. Therefore, we examined whether the evolutionary host range of plasmids could be inferred by comparing their trinucleotide composition to that of all completely sequenced bacterial chromosomes. The diversity of candidate hosts was determined using taxonomic classification and genetic distance. The method was first tested using plasmids from six incompatibility (Inc) groups whose host ranges are generally thought to be narrow (IncF, IncH, and IncI) or broad (IncN, IncP, and IncW) and then applied to other plasmid groups. The evolutionary host range was found to be broad for IncP plasmids, narrow for IncF and IncI plasmids, and intermediate for IncH and IncN plasmids, which corresponds with their known host range. The IncW plasmids as well as several plasmids from the IncA/C, IncP, IncQ, IncU, and PromA groups have signatures that were not similar to any of the chromosomal signatures, raising the hypothesis that these plasmids have not been ameliorated in any host due to their promiscuous nature. The inferred evolutionary host range of IncA/C, IncP-9, and IncL/M plasmids requires further investigation. In this era of high-throughput sequencing, this genomic signature method is a useful tool for predicting the host range of novel mobile elements.

Positive selection at high temperature reduces gene transcription in the bacteriophage ϕX174.

BACKGROUND: Gene regulation plays a central role in the adaptation of organisms to their environments. There are many molecular components to gene regulation, and it is often difficult to determine both the genetic basis of adaptation and the evolutionary forces that influence regulation. In multiple evolution experiments with the bacteriophage ϕX174, adaptive substitutions in cis-acting regulatory sequences sweep through the phage population as the result of strong positive selection at high temperatures that are non-permissive for laboratory-adapted phage. For one cis-regulatory region, we investigate the individual effects of four adaptive substitutions on transcript levels and fitness for phage growing on three hosts at two temperatures. RESULTS: The effect of the four individual substitutions on transcript levels is to down-regulate gene expression, regardless of temperature or host. To ascertain the conditions under which these substitutions are adaptive, fitness was measured by a variety of methods for several bacterial hosts growing at two temperatures, the control temperature of 37°C and the selective temperature of 42°C. Time to lysis and doublings per hour indicate that the four substitutions individually improve fitness over the ancestral strain at high temperature independent of the bacterial host in which the fitness was measured. Competition assays between the ancestral strain and either of two mutant strains indicate that both mutants out-compete the ancestor at high temperature, but the relative frequencies of each phage remain the same at the control temperature. CONCLUSIONS: Our results strongly suggest that gene transcription plays an important role in influencing fitness in the bacteriophage ϕX174, and different point mutations in a single cis-regulatory region provided the genetic basis for this role in adaptation to high temperature. We speculate that the adaptive nature of these substitutions is due to the physiology of the host at high temperature or the need to maintain particular ratios of phage proteins during capsid assembly. Our investigation of regulatory evolution contributes to interpreting genome-level assessments of regulatory variation, as well as to understanding the molecular basis of adaptation.

Comparative genomics of pAKD4, the prototype IncP-1delta plasmid with a complete backbone.

Plasmids of the incompatibility group IncP-1 are important agents of horizontal gene transfer and contribute to the spread of antibiotic resistance and xenobiotic degradation within bacterial communities. Even though some prototype plasmids have been studied in much detail, the diversity of this plasmid group was still greatly underestimated until recently, as only two of the five currently known divergent sub-groups had been described. To further improve our insight into the diversity and evolutionary history of this family of broad-host-range plasmids, we compared the complete nucleotide sequence of a new IncP-1delta plasmid pAKD4 to the genomes of other IncP-1 plasmids. Plasmid pAKD4 was previously isolated by exogenous plasmid isolation from an agricultural soil in Norway. Its 56,803bp nucleotide sequence shows high similarity in gene sequence and gene order to both plasmids pEST4011 and pIJB1, the only other IncP-1delta plasmids sequenced so far. While all three plasmids have a typical IncP-1 backbone comprising replication, transfer, and stable inheritance/control genes, the low sequence similarity in some regions and presence/absence of some backbone genes compared to other IncP-1 plasmids cluster them in a divergent sub-group. Therefore this study validates the presence of a real IncP-1delta clade with multiple plasmids. Moreover, since both pEST4011 and pIJB1 are missing a portion of their transfer genes, pAKD4 represents the first completely sequenced self-transferable plasmid with a complete IncP-1delta backbone. We therefore propose it to be the prototype IncP-1delta plasmid.

Using Mahalanobis distance to compare genomic signatures between bacterial plasmids and chromosomes.

Plasmids are ubiquitous mobile elements that serve as a pool of many host beneficial traits such as antibiotic resistance in bacterial communities. To understand the importance of plasmids in horizontal gene transfer, we need to gain insight into the 'evolutionary history' of these plasmids, i.e. the range of hosts in which they have evolved. Since extensive data support the proposal that foreign DNA acquires the host's nucleotide composition during long-term residence, comparison of nucleotide composition of plasmids and chromosomes could shed light on a plasmid's evolutionary history. The average absolute dinucleotide relative abundance difference, termed delta-distance, has been commonly used to measure differences in dinucleotide composition, or 'genomic signature', between bacterial chromosomes and plasmids. Here, we introduce the Mahalanobis distance, which takes into account the variance-covariance structure of the chromosome signatures. We demonstrate that the Mahalanobis distance is better than the delta-distance at measuring genomic signature differences between plasmids and chromosomes of potential hosts. We illustrate the usefulness of this metric for proposing candidate long-term hosts for plasmids, focusing on the virulence plasmids pXO1 from Bacillus anthracis, and pO157 from Escherichia coli O157:H7, as well as the broad host range multi-drug resistance plasmid pB10 from an unknown host.