Mobilome-driven partitions of the resistome in Salmonella.

IMPORTANCE: Antimicrobial resistance (AMR) has become a significant global challenge, with an estimated 10 million deaths annually by 2050. The emergence of AMR is mainly attributed to mobile genetic elements (MGEs or mobilomes), which accelerate wide dissemination among pathogens. The interaction between mobilomes and AMR genes (or resistomes) in Salmonella, a primary cause of diarrheal diseases that results in over 90 million cases annually, remains poorly understood. The available fragmented or incomplete genomes remain a significant limitation in investigating the relationship between AMR and MGEs. Here, we collected the most extensive closed Salmonella genomes (n = 1,817) from various sources across 58 countries. Notably, our results demonstrate that resistome transmission between Salmonella lineages follows a specific pattern of MGEs and is influenced by external drivers, including certain socioeconomic factors. Therefore, targeted interventions are urgently needed to mitigate the catastrophic consequences of Salmonella AMR.

Compensatory mutations reducing the fitness cost of plasmid carriage occur in plant rhizosphere communities.

Plasmids drive bacterial evolutionary innovation by transferring ecologically important functions between lineages, but acquiring a plasmid often comes at a fitness cost to the host cell. Compensatory mutations, which ameliorate the cost of plasmid carriage, promote plasmid maintenance in simplified laboratory media across diverse plasmid-host associations. Whether such compensatory evolution can occur in more complex communities inhabiting natural environmental niches where evolutionary paths may be more constrained is, however, unclear. Here, we show a substantial fitness cost of carrying the large conjugative plasmid pQBR103 in Pseudomonas fluorescens SBW25 in the plant rhizosphere. This plasmid fitness cost could be ameliorated by compensatory mutations affecting the chromosomal global regulatory system gacA/gacS, which arose rapidly in plant rhizosphere communities and were exclusive to plasmid carriers. These findings expand our understanding of the importance of compensatory evolution in plasmid dynamics beyond simplified lab media. Compensatory mutations contribute to plasmid survival in bacterial populations living within complex microbial communities in their environmental niche.

Systematic In Silico Assessment of Antimicrobial Resistance Dissemination across the Global Plasmidome.

The emergence of pathogenic strains resistant to multiple antimicrobials is a pressing problem in modern healthcare. Antimicrobial resistance is mediated primarily by dissemination of resistance determinants via horizontal gene transfer. The dissemination of some resistance genes has been well documented, but few studies have analyzed the patterns underpinning the dissemination of antimicrobial resistance genes. Analyzing the %GC content of plasmid-borne antimicrobial resistance genes relative to their host genome %GC content provides a means to efficiently detect and quantify dissemination of antimicrobial resistance genes. In this work we automate %GC content analysis to perform a comprehensive analysis of known antimicrobial resistance genes in publicly available plasmid sequences. We find that the degree to which antimicrobial resistance genes are disseminated depends primarily on the resistance mechanism. Our analysis identifies conjugative plasmids as primary dissemination vectors and indicates that most broadly disseminated genes have spread from single genomic backgrounds. We show that resistance dissemination profiles vary greatly among antimicrobials, oftentimes reflecting stewardship measures. Our findings establish %GC content analysis as a powerful, intuitive and scalable method to monitor the dissemination of resistance determinants using publicly available sequence data.

Long-term adaptation of ParA, RelE/ParE partition system, replication protein and phage proteins encoding low-cost plasmids of Escherichia species isolated from diarrheic children of North East India.

AIMS: The prevalent distribution of plasmid-mediated ß-lactam resistance is the most pressing global problem in enteric diseases. The current work aims to characterize plasmid-carrying ß-lactam resistant Enterobacteriaceae isolates from North East India for horizontal gene transfer (HGT) and plasmid adaptation study. METHODS AND RESULTS: In vitro transconjugation and transformation showed overall high conjugation frequency (4.11 × 10-1-9.2 × 10-1) and moderate transformation efficiency/µg DNA (1.02 × 102 -1 × 103), and the highest conjugation frequency (9.2 × 10-1) and transformation efficiency (1 × 103) for Escherichia species S-10. Intra/intergenus plasmid transformation efficiency was highest for the transformation of Klebsiella pneumoniae S-2 to Shigellaflexneri S-42 (1.3 × 103) and lowest for Escherichia species S-10 to Escherichia fergusonii S-30 (2 × 102). In the plasmid stability test, S-10 was detected with the highest plasmid carrying frequency (83.44%) and insignificant segregational loss rate (0.0004) until the 60th day with low plasmid cost on the host. The above findings were also validated by whole-plasmid sequencing of Escherichia species S-10. The genome was identified with two plasmids constituting multiple phage proteins, relaxosomal protein NikA, replication protein RepA, and the plasmid maintenance proteins (ParA, RelE/ParE), thus assisting stable plasmid maintenance. CONCLUSIONS: The results thus indicate that the high conjugation ability and low plasmid fitness cost might lead to horizontal gene transfer of the plasmid to the environment due to their prolonged adaptation in nonselective conditions, intensifying the infection's severity.

Antibiotic and heavy metal resistance genes in sewage sludge survive during aerobic composting.

Municipal sewage sludge has been generated in increasing amounts with the acceleration of urbanization and economic development. The nutrient rich sewage sludge can be recycled by composting that has a great potential to produce stabilized organic fertilizer and substrate for plant cultivation. However, little is known about the metals, pathogens and antibiotic resistance transfer risks involved in applying the composted sludge in agriculture. We studied changes in and relationships between heavy metal contents, microbial communities, and antibiotic resistance genes (ARGs), heavy metal resistance genes (HMRGs) and mobile genetic elements (MGEs) in aerobic composting of sewage sludge. The contents of most of the analyzed heavy metals were not lower after composting. The bacterial α-diversity was lower, and the community composition was different after composting. Firmicutes were enriched, and Proteobacteria and potential pathogens in the genera Arcobacter and Acinetobacter were depleted in the composted sludge. The differences in bacteria were possibly due to the high temperature phase during the composting which was likely to affect temperature-sensitive bacteria. The number of detected ARGs, HMRGs and MGEs was lower, and the relative abundances of several resistance genes were lower after composting. However, the abundance of seven ARGs and six HMRGs remained on the same level after composting. Co-occurrence analysis of bacterial taxa and the genes suggested that the ARGs may spread via horizontal gene transfer during composting. In summary, even though aerobic composting is effective for managing sewage sludge and to decrease the relative abundance of potential pathogens, ARGs and HMRGs, it might include a potential risk for the dissemination of ARGs in the environment.

Assessment of horizontal gene transfer-mediated destabilization of Synechococcus elongatus PCC 7942 biocontainment system.

Biological containment is a biosafety strategy that prevents the dispersal of genetically modified organisms in natural ecosystems. We previously established a biocontainment system that makes bacterial growth dependent on the availability of phosphite (Pt), an ecologically rare form of phosphorus (P), by introducing Pt metabolic pathway genes and disrupting endogenous phosphate and organic phosphate transporter genes. Although this system proved highly effective, horizontal gene transfer (HGT) mediated recovery of a P transporter gene is considered as a potential pathway to abolish the Pt-dependent growth, resulting in escape from the containment. Here, we assessed the risk of HGT driven escape using the Pt-dependent cyanobacterium Synechococcus elongatus PCC 7942. Transformation experiments revealed that the Pt-dependent strain could regain phosphate transporter genes from the S. elongatus PCC 7942 wild-type genome and from the genome of the closely related strain, S. elongatus UTEX 2973. Transformed S. elongatus PCC 7942 became viable in a phosphate-containing medium. Meanwhile, transformation of the Synechocystis sp. PCC 6803 genome or environmental DNA did not yield escape strains, suggesting that only genetic material derived from phylogenetically-close species confer high risk to generate escape. Eliminating a single gene necessary for natural competence from the Pt-dependent strain reduced the escape occurrence rate. These results demonstrate that natural competence could be a potential risk to destabilize Pt-dependence, and therefore inhibiting exogenous DNA uptake would be effective for enhancing the robustness of the gene disruption-dependent biocontainment.

Assessment of a plasmid conjugation procedure to monitor horizontal transfer of an extended-spectrum ß-lactamase resistance gene under food chain scenarios.

Plasmids are relevant reservoirs of antimicrobial resistance genes (ARGs) which confer adaptive advantages to their host and can be horizontally transferred. The aims of this study were to develop a conjugation procedure to monitor the horizontal transfer of a 193 kb plasmid containing the extended-spectrum ß-lactamase production gene bla CTX-M-14 between two Escherichia coli strains under a range of food chain-related scenarios, including temperature (20-37 °C), pH (5.0-9.0) or the presence of some biocidal agents (benzalkonium chloride, sodium hypochlorite or peracetic acid). The average conjugation rate in LB broth after 18 h at 37 °C was 2.09e-04 and similar rates were observed in a food matrix (cow's milk). The conjugation was reduced at temperatures below 37 °C, at alkaline pH (especially at pH 9.0) or in the presence of benzalkonium chloride. Peracetic acid and sodium hypochlorite slightly increased conjugation rates, which reached 5.59e-04 and 6.77e-03, respectively. The conjugation procedure described can be used to identify risk scenarios leading to an enhanced ARGs transmission via plasmid conjugation, as well as to identify novel intervention strategies impairing plasmid conjugation and tackling antimicrobial resistance.

Fitness costs of CRISPR-Cas systems in bacteria.

CRISPR-Cas systems provide bacteria with both specificity and adaptability in defence against invading genetic elements. From a theoretical perspective, CRISPR-Cas systems confer many benefits. However, they are observed at an unexpectedly low prevalence across the bacterial domain. While these defence systems can be gained horizontally, fitness costs may lead to selection against their carriage. Understanding the source of CRISPR-related fitness costs will help us to understand the evolutionary dynamics of CRISPR-Cas systems and their role in shaping bacterial genome evolution. Here, we review our current understanding of the potential fitness costs associated with CRISPR-Cas systems. In addition to potentially restricting the acquisition of genetic material that could confer fitness benefits, we explore five alternative biological factors that from a theoretical perspective may influence the fitness costs associated with CRISPR-Cas system carriage: (1) the repertoire of defence mechanisms a bacterium has available to it, (2) the potential for a metabolic burden, (3) larger-scale population and environmental factors, (4) the phenomenon of self-targeting spacers, and (5) alternative non-defence roles for CRISPR-Cas.

Performance assessment of genomic island prediction tools with an improved version of Design-Island.

Genomic Islands (GIs) play an important role in the evolution and adaptation of prokaryotes. The origin and extent of ecological diversity of prokaryotes can be analyzed by comparing GIs across closely or distantly related prokaryotes. Understanding the importance of GI and to study the bacterial evolution, several GI prediction tools have been generated. An unsupervised method, Design-Island, was developed to identify GIs using Monte-Carlo statistical test on randomly selected segments of a chromosome. Here, in the present study Design-Island was modified with the incorporation of majority voting, multiple hypothesis testing correction. The performance of the modified version, Design-Island-II was tested and compared with the existing GI prediction tools. The performance assessment and benchmarking of the GI prediction tools require experimentally validated dataset, which is lacking. So, different datasets, generated or taken from literature were utilized to compare the sensitivity (SN), specificity (SP), precision (PPV) and accuracy (AC) of Design-Island-II. It showed substantial enhancement in term of SN, SP, PPV and AC, and significantly reduced the computation time of the algorithm. The performance of Design-Island-II has also been compared with several GI prediction tools using curated dataset of putative horizontally transferred genes. Design-Island-II showed the highest sensitivity and F1 score, comparable specificity, precision and accuracy in comparison to the other available methods. IslandViewer4 and Islander outperformed all the available methods in terms of AC and PPV respectively. Our study suggested Design-Island-II, IslandViewer4 and GIHunter among the top performing GI prediction tools considering both sensitivity and specificity of the methods.

Recombination resolves the cost of horizontal gene transfer in experimental populations of Helicobacter pylori.

Horizontal gene transfer (HGT) is important for microbial evolution, yet we know little about the fitness effects and dynamics of horizontally transferred genetic variants. In this study, we evolve laboratory populations of Helicobacter pylori, which take up DNA from their environment by natural transformation, and measure the fitness effects of thousands of transferred genetic variants. We find that natural transformation increases the rate of adaptation but comes at the cost of significant genetic load. We show that this cost is circumvented by recombination, which increases the efficiency of selection by decoupling deleterious and beneficial genetic variants. Our results show that adaptation with HGT, pervasive in natural microbial populations, is shaped by a combination of selection, recombination, and genetic drift not accounted for in existing models of evolution.

Large-scale characterization of the macrolide resistome reveals high diversity and several new pathogen-associated genes.

Macrolides are broad-spectrum antibiotics used to treat a range of infections. Resistance to macrolides is often conferred by mobile resistance genes encoding Erm methyltransferases or Mph phosphotransferases. New erm and mph genes keep being discovered in clinical settings but their origins remain unknown, as is the type of macrolide resistance genes that will appear in the future. In this study, we used optimized hidden Markov models to characterize the macrolide resistome. Over 16 terabases of genomic and metagenomic data, representing a large taxonomic diversity (11 030 species) and diverse environments (1944 metagenomic samples), were searched for the presence of erm and mph genes. From this data, we predicted 28 340 macrolide resistance genes encoding 2892 unique protein sequences, which were clustered into 663 gene families (<70 % amino acid identity), of which 619 (94 %) were previously uncharacterized. This included six new resistance gene families, which were located on mobile genetic elements in pathogens. The function of ten predicted new resistance genes were experimentally validated in Escherichia coli using a growth assay. Among the ten tested genes, seven conferred increased resistance to erythromycin, with five genes additionally conferring increased resistance to azithromycin, showing that our models can be used to predict new functional resistance genes. Our analysis also showed that macrolide resistance genes have diverse origins and have transferred horizontally over large phylogenetic distances into human pathogens. This study expands the known macrolide resistome more than ten-fold, provides insights into its evolution, and demonstrates how computational screening can identify new resistance genes before they become a significant clinical problem.

Translational demand is not a major source of plasmid-associated fitness costs.

Plasmids are key drivers of bacterial evolution because they are crucial agents for the horizontal transfer of adaptive traits, such as antibiotic resistance. Most plasmids entail a metabolic burden that reduces the fitness of their host if there is no selection for plasmid-encoded genes. It has been hypothesized that the translational demand imposed by plasmid-encoded genes is a major mechanism driving the fitness cost of plasmids. Plasmid-encoded genes typically present a different codon usage from host chromosomal genes. As a consequence, the translation of plasmid-encoded genes might sequestrate ribosomes on plasmid transcripts, overwhelming the translation machinery of the cell. However, the pervasiveness and origins of the translation-derived costs of plasmids are yet to be assessed. Here, we systematically altered translation efficiency in the host cell to disentangle the fitness effects produced by six natural antibiotic resistance plasmids. We show that limiting translation efficiency either by reducing the number of available ribosomes or their processivity does not increase plasmid costs. Overall, our results suggest that ribosomal paucity is not a major contributor to plasmid fitness costs. This article is part of the theme issue 'The secret lives of microbial mobile genetic elements'.

Poorly Conserved P15 Proteins of Cileviruses Retain Elements of Common Ancestry and Putative Functionality: A Theoretical Assessment on the Evolution of Cilevirus Genomes.

The genus Cilevirus groups enveloped single-stranded (+) RNA virus members of the family Kitaviridae, order Martellivirales. Proteins P15, scarcely conserved polypeptides encoded by cileviruses, have no apparent homologs in public databases. Accordingly, the open reading frames (ORFs) p15, located at the 5'-end of the viral RNA2 molecules, are considered orphan genes (ORFans). In this study, we have delved into ORFs p15 and the relatively poorly understood biochemical properties of the proteins P15 to posit their importance for viruses across the genus and theorize on their origin. We detected that the ORFs p15 are under purifying selection and that, in some viral strains, the use of synonymous codons is biased, which might be a sign of adaptation to their plant hosts. Despite the high amino acid sequence divergence, proteins P15 show the conserved motif [FY]-L-x(3)-[FL]-H-x-x-[LIV]-S-C-x-C-x(2)-C-x-G-x-C, which occurs exclusively in members of this protein family. Proteins P15 also show a common predicted 3D structure that resembles the helical scaffold of the protein ORF49 encoded by radinoviruses and the phosphoprotein C-terminal domain of mononegavirids. Based on the 3D structural similarities of P15, we suggest elements of common ancestry, conserved functionality, and relevant amino acid residues. We conclude by postulating a plausible evolutionary trajectory of ORFans p15 and the 5'-end of the RNA2 of cileviruses considering both protein fold superpositions and comparative genomic analyses with the closest kitaviruses, negeviruses, nege/kita-like viruses, and unrelated viruses that share the ecological niches of cileviruses.

Quantitative Microbial Risk Assessment of Antimicrobial Resistant and Susceptible Staphylococcus aureus in Reclaimed Wastewaters.

The annual risks of colonization, skin infection, bloodstream infection (BSI), and disease burden from exposures to antibiotic-resistant and susceptible Staphylococcus aureus (S. aureus) were estimated using quantitative microbial risk assessment (QMRA). We estimated the probability of nasal colonization after immersion in wastewater (WW) or greywater (GW) treated across a range of treatment alternatives and subsequent infection. Horizontal gene transfer was incorporated into the treatment model but had little effect on the predicted risk. The cumulative annual probability of infection (resulting from self-inoculation) was most sensitive to the treatment log10 reduction value (LRV), S. aureus concentration, and the newly calculated morbidity ratios and was below the health benchmark of 10-4 infections per person per year (ppy) given a treatment LRV of roughly 3.0. The predicted annual disability-adjusted life years (DALYs), which were dominated by BSI, were below the health benchmark of 10-6 DALYs ppy for resistant and susceptible S. aureus, given LRVs of 4.5 and 3.5, respectively. Thus, the estimated infection risks and disease burdens resulting from nasal colonization are below the relevant health benchmarks for risk-based, nonpotable, or potable reuse systems but possibly above for immersion in minimally treated GW or WW. Strain-specific data to characterize dose-response and concentration in WW are needed to substantiate the QMRA.

Conjugation dynamics depend on both the plasmid acquisition cost and the fitness cost.

Plasmid conjugation is a major mechanism responsible for the spread of antibiotic resistance. Plasmid fitness costs are known to impact long-term growth dynamics of microbial populations by providing plasmid-carrying cells a relative (dis)advantage compared to plasmid-free counterparts. Separately, plasmid acquisition introduces an immediate, but transient, metabolic perturbation. However, the impact of these short-term effects on subsequent growth dynamics has not previously been established. Here, we observed that de novo transconjugants grew significantly slower and/or with overall prolonged lag times, compared to lineages that had been replicating for several generations, indicating the presence of a plasmid acquisition cost. These effects were general to diverse incompatibility groups, well-characterized and clinically captured plasmids, Gram-negative recipient strains and species, and experimental conditions. Modeling revealed that both fitness and acquisition costs modulate overall conjugation dynamics, validated with previously published data. These results suggest that the hours immediately following conjugation may play a critical role in both short- and long-term plasmid prevalence. This time frame is particularly relevant to microbiomes with high plasmid/strain diversity considered to be hot spots for conjugation.

A Novel Small RNA on the Pseudomonas putida KT2440 Chromosome Is Involved in the Fitness Cost Imposed by IncP-1 Plasmid RP4.

Plasmids can provide advantageous traits to host bacteria, although they may impose a fitness cost. Chromosome-encoded factors are important for regulating the expression of genes on plasmids, and host chromosomes may differ in terms of their interactions with a given plasmid. Accordingly, differences in fitness cost loading and compensatory co-evolution may occur for various host chromosome/plasmid combinations. However, the mechanisms of compensatory evolution are highly divergent and require further insights. Here, we reveal novel evolutionally mechanisms of Pseudomonas putida KT2440 to improve the fitness cost imposed by the incompatibility P-1 (IncP-1) multidrug resistance plasmid RP4. A mixed culture of RP4-harboring and -free KT2440 cells was serially transferred every 24 h under non-selective conditions. Initially, the proportion of RP4-harboring cells decreased rapidly, but it immediately recovered, suggesting that the fitness of RP4-harboring strains improved during cultivation. Larger-sized colonies appeared during 144-h mixed culture, and evolved strains isolated from larger-sized colonies showed higher growth rates and fitness than those of the ancestral strain. Whole-genome sequencing revealed that evolved strains had one of two mutations in the same intergenic region of the chromosome. Based on the research of another group, this region is predicted to contain a stress-inducible small RNA (sRNA). Identification of the transcriptional start site in this sRNA indicated that one mutation occurred within the sRNA region, whereas the other was in its promoter region. Quantitative reverse-transcription PCR showed that the expression of this sRNA was strongly induced by RP4 carriage in the ancestral strain but repressed in the evolved strains. When the sRNA region was overexpressed in the RP4-free strain, the fitness decreased, and the colony size became smaller. Using transcriptome analysis, we also showed that the genes involved in amino acid metabolism and stress responses were differentially transcribed by overexpression of the sRNA region. These results indicate that the RP4-inducible chromosomal sRNA was responsible for the fitness cost of RP4 on KT2440 cells, where this sRNA is of key importance in host evolution toward rapid amelioration of the cost.

Extremely fast amelioration of plasmid fitness costs by multiple functionally diverse pathways.

The acquisition of plasmids is often accompanied by fitness costs such that compensatory evolution is required to allow plasmid survival, but it is unclear whether compensatory evolution can be extensive or rapid enough to maintain plasmids when they are very costly. The mercury-resistance plasmid pQBR55 drastically reduced the growth of its host, Pseudomonas fluorescens SBW25, immediately after acquisition, causing a small colony phenotype. However, within 48 h of growth on agar plates we observed restoration of the ancestral large colony morphology, suggesting that compensatory mutations had occurred. Relative fitness of these evolved strains, in lab media and in soil microcosms, varied between replicates, indicating different mutational mechanisms. Using genome sequencing we identified that restoration was associated with chromosomal mutations in either a hypothetical DNA-binding protein PFLU4242, RNA polymerase or the GacA/S two-component system. Targeted deletions in PFLU4242, gacA or gacS recapitulated the ameliorated phenotype upon plasmid acquisition, indicating three distinct mutational pathways to compensation. Our data shows that plasmid compensatory evolution is fast enough to allow survival of a plasmid despite it imposing very high fitness costs upon its host, and indeed may regularly occur during the process of isolating and selecting individual plasmid-containing clones.

Assessment of Antibiotic Levels, Multi-Drug Resistant Bacteria and Genetic Biomarkers in the Waters of the Rio Grande River Between the United States-Mexico Border.

BACKGROUND: The worldwide emergence of multi-drug resistant bacteria has become a health crisis, as fewer or sometimes no antimicrobial agents are effective against these bacteria. The Rio Grande River is the natural boundary between the United States (US) and Mexico. It spans a border region between Texas, New Mexico and Mexico. Underserved populations on the Mexican side use the river for recreational purposes, while on the US side, the river is used for irrigation and as a source of drinking water. OBJECTIVES: The purpose of the present study was to evaluate the concentration of antibiotic residues, to determine the presence of genetic elements conferring antibiotic resistance and to characterize multi-drug resistant bacteria in the waters of the Rio Grande River. METHODS: Water samples were obtained from the Rio Grande River. Deoxyribonucleic acid (DNA) was extracted from both isolated bacteria and directly from the water. Amplification of selected genetic elements was accomplished by polymerase chain reaction. Identification and isolation of bacteria was performed through MicroScan autoSCAN-4. Fecal contamination was assessed by IDEXX Colilert. Antibiotic residues were determined by liquid chromatography and mass spectrometry. RESULTS: Antibiotics were found in 92% of both water and sediment samples. Antibiotic concentrations ranged from 0.38 ng/L - 742.73 ng/L and 0.39 ng/l - 66.3 ng/g dry weight in water and sediment samples, respectively. Genetic elements conferring resistance were recovered from all collection sites. Of the isolated bacteria, 91 (64.08%) were resistant to at least two synergistic antibiotic combinations and 11 (14.79%) were found to be resistant to 20 or more individual antibiotics. Fecal contamination was higher during the months of April and July. CONCLUSIONS: The 26 km segment of the Rio Grande River from Sunland Park NM to El Paso, TX and Juarez, Mexico is an area of concern due to poor water quality. The presence of multidrug resistant bacteria, antibiotics and mobile genetic elements may be a health hazard for the surrounding populations of this binational border region. Policies need to be developed for the appropriate management of the environmental natural resources in this border region. COMPETING INTERESTS: The authors declare no competing financial interests.

IslandCafe: Compositional Anomaly and Feature Enrichment Assessment for Delineation of Genomic Islands.

One of the evolutionary forces driving bacterial genome evolution is the acquisition of clusters of genes through horizontal gene transfer (HGT). These genomic islands may confer adaptive advantages to the recipient bacteria, such as, the ability to thwart antibiotics, become virulent or hypervirulent, or acquire novel metabolic traits. Methods for detecting genomic islands either search for markers or features typical of islands or examine anomaly in oligonucleotide composition against the genome background. The former tends to underestimate, missing islands that have the markers either lost or degraded, while the latter tends to overestimate, due to their inability to discriminate compositional atypicality arising because of HGT from those that are a consequence of other biological factors. We propose here a framework that exploits the strengths of both these approaches while bypassing the pitfalls of either. Genomic islands lacking markers are identified by their association with genomic islands with markers. This was made possible by performing marker enrichment and phyletic pattern analyses within an integrated framework of recursive segmentation and clustering. The proposed method, IslandCafe, compared favorably with frequently used methods for genomic island detection on synthetic test datasets and on a test-set of known islands from 15 well-characterized bacterial species. Furthermore, IslandCafe identified novel islands with imprints of likely horizontal acquisition.

Protein Structure-Guided Hidden Markov Models (HMMs) as A Powerful Method in the Detection of Ancestral Endogenous Viral Elements.

It has been believed for a long time that the transfer and fixation of genetic material from RNA viruses to eukaryote genomes is very unlikely. However, during the last decade, there have been several cases in which "virus-to-host" gene transfer from various viral families into various eukaryotic phyla have been described. These transfers have been identified by sequence similarity, which may disappear very quickly, especially in the case of RNA viruses. However, compared to sequences, protein structure is known to be more conserved. Applying protein structure-guided protein domain-specific Hidden Markov Models, we detected homologues of the Virgaviridae capsid protein in Schizophora flies. Further data analysis supported "virus-to-host" transfer into Schizophora ancestors as a single transfer event. This transfer was not identifiable by BLAST or by other methods we applied. Our data show that structure-guided Hidden Markov Models should be used to detect ancestral virus-to-host transfers.