The integrated genomic surveillance system of Andalusia (SIEGA) provides a One Health regional resource connected with the clinic.

The One Health approach, recognizing the interconnectedness of human, animal, and environmental health, has gained significance amid emerging zoonotic diseases and antibiotic resistance concerns. This paper aims to demonstrate the utility of a collaborative tool, the SIEGA, for monitoring infectious diseases across domains, fostering a comprehensive understanding of disease dynamics and risk factors, highlighting the pivotal role of One Health surveillance systems. Raw whole-genome sequencing is processed through different species-specific open software that additionally reports the presence of genes associated to anti-microbial resistances and virulence. The SIEGA application is a Laboratory Information Management System, that allows customizing reports, detect transmission chains, and promptly alert on alarming genetic similarities. The SIEGA initiative has successfully accumulated a comprehensive collection of more than 1900 bacterial genomes, including Salmonella enterica, Listeria monocytogenes, Campylobacter jejuni, Escherichia coli, Yersinia enterocolitica and Legionella pneumophila, showcasing its potential in monitoring pathogen transmission, resistance patterns, and virulence factors. SIEGA enables customizable reports and prompt detection of transmission chains, highlighting its contribution to enhancing vigilance and response capabilities. Here we show the potential of genomics in One Health surveillance when supported by an appropriate bioinformatic tool. By facilitating precise disease control strategies and antimicrobial resistance management, SIEGA enhances global health security and reduces the burden of infectious diseases. The integration of health data from humans, animals, and the environment, coupled with advanced genomics, underscores the importance of a holistic One Health approach in mitigating health threats.

Genotypic characterization and antimicrobial susceptibility of human Campylobacter jejuni isolates in Southern Spain.

Campylobacter jejuni is the main cause of bacterial gastroenteritis and a public health problem worldwide. Little information is available on the genotypic characteristics of human C. jejuni in Spain. This study is based on an analysis of the resistome, virulome, and phylogenetic relationship, antibiogram prediction, and antimicrobial susceptibility of 114 human isolates of C. jejuni from a tertiary hospital in southern Spain from October 2020 to June 2023. The isolates were sequenced using Illumina technology, and a bioinformatic analysis was subsequently performed. The susceptibility of C. jejuni isolates to ciprofloxacin, tetracycline, and erythromycin was also tested. The resistance rates for each antibiotic were 90.3% for ciprofloxacin, 66.7% for tetracycline, and 0.88% for erythromycin. The fluoroquinolone resistance rate obtained is well above the European average (69.1%). CC-21 (n = 23), ST-572 (n = 13), and ST-6532 (n = 13) were the most prevalent clonal complexes (CCs) and sequence types (STs). In the virulome, the cadF, ciaB, and cdtABC genes were detected in all the isolates. A prevalence of 20.1% was obtained for the genes wlaN and cstIII, which are related to the pathogenesis of Guillain-Barré syndrome (GBS). The prevalence of the main antimicrobial resistance markers detected were CmeABC (92.1%), RE-cmeABC (7.9%), the T86I substitution in gyrA (88.9%), blaOXA-61 (72.6%), tet(O) (65.8%), and ant (6)-Ia (17.1%). High antibiogram prediction rates (>97%) were obtained, except for in the case of the erythromycin-resistant phenotype. This study contributes significantly to the knowledge of C. jejuni genomics for the prevention, treatment, and control of infections caused by this pathogen.IMPORTANCEDespite being the pathogen with the greatest number of gastroenteritis cases worldwide, Campylobacter jejuni remains a poorly studied microorganism. A sustained increase in fluoroquinolone resistance in human isolates is a problem when treating Campylobacter infections. The development of whole genome sequencing (WGS) techniques has allowed us to better understand the genotypic characteristics of this pathogen and relate them to antibiotic resistance phenotypes. These techniques complement the data obtained from the phenotypic analysis of C. jejuni isolates. The zoonotic transmission of C. jejuni through the consumption of contaminated poultry supports approaching the study of this pathogen through "One Health" approach. In addition, due to the limited information on the genomic characteristics of C. jejuni in Spain, this study provides important data and allows us to compare the results with those obtained in other countries.

Molecular and phylogenetic characterization of the monkeypox outbreak in the South of Spain.

Background and Aim: Until the May 2022 Monkeypox (MPXV) outbreak, which spread rapidly to many non-endemic countries, the virus was considered a viral zoonosis limited to some African countries. The Andalusian circuit of genomic surveillance was rapidly applied to characterize the MPXV outbreak in the South of Spain. Methods: Whole genome sequencing was used to obtain the genomic profiles of samples collected across the south of Spain, representative of all the provinces of Andalusia. Phylogenetic analysis was used to study the relationship of the isolates and the available sequences of the 2022 outbreak. Results: Whole genome sequencing of a total of 160 MPXV viruses from the different provinces that reported cases were obtained. Interestingly, we report the sequences of MPXV viruses obtained from two patients who died. While one of the isolates bore no noteworthy mutations that explain a potential heightened virulence, in another patient the second consecutive genome sequence, performed after the administration of tecovirimat, uncovered a mutation within the A0A7H0DN30 gene, known to be a prime target for tecovirimat in its Vaccinia counterpart. In general, a low number of mutations were observed in the sequences reported, which were very similar to the reference of the 2022 outbreak (OX044336), as expected from a DNA virus. The samples likely correspond to several introductions of the circulating MPXV viruses from the last outbreak. The virus sequenced from one of the two patients that died presented a mutation in a gene that bears potential connections to drug resistance. This mutation was absent in the initial sequencing before treatment.

C-terminal deletion of RelA protein is suggested as a possible cause of infective endocarditis recurrence with Enterococcus faecium.

Infective endocarditis (IE) caused by Enterococcus spp. represents the third most common cause of IE, with high rates of relapse compared with other bacteria. Interestingly, late relapses (>6 months) have only been described in Enterococcus faecalis, but here we describe the first reported IE relapse with Enterococcus faecium more than a year (17 months) after the initial endocarditis episode. Firstly, by multi locus sequence typing (MLST), we demonstrated that both isolates (EF646 and EF641) belong to the same sequence type (ST117). Considering that EF641 was able to overcome starvation and antibiotic treatment conditions surviving for a long period of time, we performed bioinformatic analysis in identifying potential genes involved in virulence and stringent response. Our results showed a 13-nucleotide duplication (positions 1638-1650) in the gene relA, resulting in a premature stop codon, with a loss of 167 amino acids from the C-terminal domains of the RelA enzyme. RelA mediates the stringent response in bacteria, modulating levels of the alarmone guanosine tetraphosphate (ppGpp). The relA mutant (EF641) was associated with lower growth capacity, the presence of small colony variants, and higher capacity to produce biofilms (compared with the strain EF646), but without differences in antimicrobial susceptibility patterns according to standard procedures during planktonic growth. Instead, EF641 demonstrated tolerance to high doses of teicoplanin when growing in a biofilm. We conclude that all these events would be closely related to the long-term survival of the E. faecium and the late relapse of the IE. These data represent the first clinical evidence of mutations in the stringent response (relA gene) related with E. faecium IE relapse.

A case of fatal monkeypox infection: necropsy and molecular findings, with some considerations related to clinical management.

OBJECTIVE: Human monkeypox (mpox) is usually self-limited infection; however, rising data show a worse outcome in patients with impaired immune status, particularly those co-infected with HIV [Mitjà O, Alemany A, Marks M, Lezama Mora JI, Rodríguez-Aldama JC, Torres Silva MS et al. Mpox in people with advanced HIV infection: A global case series. Lancet. 2023; 401:939-49. DOI:https://doi.org/10.1016/S0140-6736(23)00273-8] [Govind A, Lazarte SM, Kitchell E, Chow JY, Estelle CD, Fixsen E et al. Severe mpox infections in people with uncontrolled human immunodeficiency virus (HIV). Clin Infect Dis. 2023; 76:1843-6. DOI:https://doi.org/10.1093/cid/ciad052]. METHODS: We report the clinical, pathological, and molecular study of a patient with mpox infection and a late HIV diagnosis, with fatal outcome. RESULTS: Necropsy revealed visceral spread of mpox. Mpox virus was sequenced twice during the admission, uncovering an emerging mutation near a genomic region where mutations associated with tecovirimat resistance have been documented. DISCUSSION: Monkeypox can manifest as an opportunistic infection in individuals with advanced HIV-associated immunosuppression.

Role of blaTEM and OmpC in the piperacillin-tazobactam resistance evolution by E. coli in patients with complicated intra-abdominal infection.

Piperacillin-tazobactam resistance (P/T-R) is increasingly reported among Escherichia coli isolates. Although in vitro experiments have suggested that blaTEM gene plays a key role in the P/T-R acquisition, no clinical in vivo study has yet confirmed the role of blaTEM or other genes. Therefore, we aimed to identify the mechanisms underlying P/T-R by following up patients with E. coli complicated intra-abdominal infections (cIAI) who experienced P/T treatment failure. Four pairs of strains, clonally related from four patients, were isolated both before and after treatment with P/T dosed at 4 g/0.5 g intravenously. The P/T MIC was tested using broth microdilution, and ß-lactamase activity was determined in these isolates. Whole-genome sequencing (WGS) was performed to decipher the role of blaTEM and other genes associated with P/T-R. Changes in the outer membrane protein (OMP) profile were analyzed using SDS-PAGE, and blaTEM and ompC transcription levels were measured by RT-qPCR. In addition, in vitro competition fitness was performed between each pairs of strains (P/T-susceptible vs. P/T-resistant). We found a higher copy number of blaTEM gene in P/T-R isolates, generated by three different genetic events: (1) IS26-mediated duplication of the blaTEM gene, (2) generation of a small multicopy plasmid (ColE-like) carrying blaTEM, and (3) adaptive evolution via reduction of plasmid size, leading to a higher plasmid copy number. Moreover, two P/T-R strains showed reduced expression of OmpC. This study describes the mechanisms involved in the acquisition of P/T-R by E. coli in patients with cIAI. The understanding of P/T-R evolution is crucial for effectively treating infected patients and preventing the spread of resistant microorganisms.

Case report: Analysis of phage therapy failure in a patient with a Pseudomonas aeruginosa prosthetic vascular graft infection.

Clinical case of a patient with a Pseudomonas aeruginosa multidrug-resistant prosthetic vascular graft infection which was treated with a cocktail of phages (PT07, 14/01, and PNM) in combination with ceftazidime-avibactam (CZA). After the application of the phage treatment and in absence of antimicrobial therapy, a new P. aeruginosa bloodstream infection (BSI) with a septic residual limb metastasis occurred, now involving a wild-type strain being susceptible to ß-lactams and quinolones. Clinical strains were analyzed by microbiology and whole genome sequencing techniques. In relation with phage administration, the clinical isolates of P. aeruginosa before phage therapy (HE2011471) and post phage therapy (HE2105886) showed a clonal relationship but with important genomic changes which could be involved in the resistance to this therapy. Finally, phenotypic studies showed a decrease in Minimum Inhibitory Concentration (MIC) to ß-lactams and quinolones as well as an increase of the biofilm production and phage resistant mutants in the clinical isolate of P. aeruginosa post phage therapy.

High-level delafloxacin resistance through the combination of two different mechanisms in Staphylococcus aureus.

Delafloxacin is a new fluoroquinolone indicated for the treatment of complicated bacterial skin infections caused by Staphylococcus aureus. Despite its recent approval by the US Food and Drug Administration, the emergence of S. aureus-resistant strains has been reported. As such, this study aimed to investigate the activity of delafloxacin against a collection of S. aureus, and to determine the mechanisms of resistance. The activity of delafloxacin was measured in 59 S. aureus clinical isolates [40 methicillin-resistant S. aureus (MRSA) and 19 methicillin-susceptible S. aureus (MSSA)]. Whole-genome sequencing (WGS) was performed in the isolates resistant to delafloxacin. The minimum inhibitory concentrations required to inactivate 50% and 90% of the isolates (MIC50 and MIC90, respectively) were higher in MRSA (0.19 mg/L and 0.75 mg/L, respectively) than in MSSA (0.008 mg/L and 0.25 mg/L, respectively). Furthermore, 10 S. aureus clinical isolates (16.9%) were categorized as resistant to delafloxacin. Regarding the WGS data, several mutations were found in the quinolone resistance-determining region. Nevertheless, a mutation in the same position (E84K and E84V) of topoisomerase IV (ParC) was found exclusively in the four high-level delafloxacin-resistant isolates. Interestingly, a plasmid-encoded qacC gene (efflux pump) was found to be harboured by the isolate with the highest delafloxacin MIC value (32 mg/L). The use of a wide-spectrum efflux pump inhibitor revealed an important contribution of this system to the acquisition of delafloxacin resistance. In conclusion, delafloxacin has activity against S. aureus, including MRSA. However, this study showed that mutations in position 84 of ParC and the acquisition of a QacC efflux pump are key factors for the development of delafloxacin resistance in S. aureus.

Detection of High Level of Co-Infection and the Emergence of Novel SARS CoV-2 Delta-Omicron and Omicron-Omicron Recombinants in the Epidemiological Surveillance of Andalusia.

Recombination is an evolutionary strategy to quickly acquire new viral properties inherited from the parental lineages. The systematic survey of the SARS-CoV-2 genome sequences of the Andalusian genomic surveillance strategy has allowed the detection of an unexpectedly high number of co-infections, which constitute the ideal scenario for the emergence of new recombinants. Whole genome sequence of SARS-CoV-2 has been carried out as part of the genomic surveillance programme. Sample sources included the main hospitals in the Andalusia region. In addition to the increase of co-infections and known recombinants, three novel SARS-CoV-2 delta-omicron and omicron-omicron recombinant variants with two break points have been detected. Our observations document an epidemiological scenario in which co-infection and recombination are detected more frequently. Finally, we describe a family case in which co-infection is followed by the detection of a recombinant made from the two co-infecting variants. This increased number of recombinants raises the risk of emergence of recombinant variants with increased transmissibility and pathogenicity.

Assessing the Impact of SARS-CoV-2 Lineages and Mutations on Patient Survival.

OBJECTIVES: More than two years into the COVID-19 pandemic, SARS-CoV-2 still remains a global public health problem. Successive waves of infection have produced new SARS-CoV-2 variants with new mutations for which the impact on COVID-19 severity and patient survival is uncertain. METHODS: A total of 764 SARS-CoV-2 genomes, sequenced from COVID-19 patients, hospitalized from 19th February 2020 to 30 April 2021, along with their clinical data, were used for survival analysis. RESULTS: A significant association of B.1.1.7, the alpha lineage, with patient mortality (log hazard ratio (LHR) = 0.51, C.I. = [0.14,0.88]) was found upon adjustment by all the covariates known to affect COVID-19 prognosis. Moreover, survival analysis of mutations in the SARS-CoV-2 genome revealed 27 of them were significantly associated with higher mortality of patients. Most of these mutations were located in the genes coding for the S, ORF8, and N proteins. CONCLUSIONS: This study illustrates how a combination of genomic and clinical data can provide solid evidence for the impact of viral lineage on patient survival.

MDR Shigella sonnei in Spain: an ever-evolving emerging threat?

Background: Seven CTX-M-27-producing Shigella sonnei strains were isolated at the University Hospital Virgen del Rocío (Seville, Spain) microbiology service from October to November 2021. Objectives: To offer extensive information on the microbiological and molecular epidemiology results of the seven S. sonnei isolates and compare them with other previously documented CTX-M-27-producing S. sonnei associated with MSM transmission. Methods: S. sonnei isolated from stool samples of patients with acute diarrhoea were identified through biochemical and serological typing. Whole characterization of the seven isolates was performed by sequencing with MinION Mk1C followed by genomic and molecular analysis. Results: All the isolates were resistant to penicillins, cephalosporins, fluoroquinolones, cotrimoxazole and azithromycin. Sequencing showed the presence of several resistance determinants, outstanding bla CTX-M-27, azithromycin resistance genes [ermB and mph(A)], qnrB19 and mutations in the QRDRs. All isolates belonged to the same hierarchical clustering of cgMLST (HierCC) with five allele distance (HC5) scheme v1 from EnteroBase. However, they presented differences in plasmid composition, with all seven isolates harbouring IncFII, IncB/O/K/Z and ColE1-like while SH2, SH6 and SH7 had IncFIB only. Our isolates were closely related to others from Spain (HC5; 98748), Australia (HC5; 98748) and the UK (HC5; 98748), which were also associated with MSM transmission. Nevertheless, the structure of the non-chromosomal genetic elements and the genetic context of bla CTX-M-27 presented a certain variability compared with isolates from other countries and among them. Conclusions: This study confirms the emergence of CTX-M-27-producing S. sonnei (ST152) associated with MSM transmission in Spain, adding it to the Europe outbreak list and reinforcing the necessity of active surveillance and control of this high-risk clone.

Towards a metagenomics machine learning interpretable model for understanding the transition from adenoma to colorectal cancer.

Gut microbiome is gaining interest because of its links with several diseases, including colorectal cancer (CRC), as well as the possibility of being used to obtain non-intrusive predictive disease biomarkers. Here we performed a meta-analysis of 1042 fecal metagenomic samples from seven publicly available studies. We used an interpretable machine learning approach based on functional profiles, instead of the conventional taxonomic profiles, to produce a highly accurate predictor of CRC with better precision than those of previous proposals. Moreover, this approach is also able to discriminate samples with adenoma, which makes this approach very promising for CRC prevention by detecting early stages in which intervention is easier and more effective. In addition, interpretable machine learning methods allow extracting features relevant for the classification, which reveals basic molecular mechanisms accounting for the changes undergone by the microbiome functional landscape in the transition from healthy gut to adenoma and CRC conditions. Functional profiles have demonstrated superior accuracy in predicting CRC and adenoma conditions than taxonomic profiles and additionally, in a context of explainable machine learning, provide useful hints on the molecular mechanisms operating in the microbiota behind these conditions.

Highly accurate whole-genome imputation of SARS-CoV-2 from partial or low-quality sequences.

BACKGROUND: The current SARS-CoV-2 pandemic has emphasized the utility of viral whole-genome sequencing in the surveillance and control of the pathogen. An unprecedented ongoing global initiative is producing hundreds of thousands of sequences worldwide. However, the complex circumstances in which viruses are sequenced, along with the demand of urgent results, causes a high rate of incomplete and, therefore, useless sequences. Viral sequences evolve in the context of a complex phylogeny and different positions along the genome are in linkage disequilibrium. Therefore, an imputation method would be able to predict missing positions from the available sequencing data. RESULTS: We have developed the impuSARS application, which takes advantage of the enormous number of SARS-CoV-2 genomes available, using a reference panel containing 239,301 sequences, to produce missing data imputation in viral genomes. ImpuSARS was tested in a wide range of conditions (continuous fragments, amplicons or sparse individual positions missing), showing great fidelity when reconstructing the original sequences, recovering the lineage with a 100% precision for almost all the lineages, even in very poorly covered genomes (<20%). CONCLUSIONS: Imputation can improve the pace of SARS-CoV-2 sequencing production by recovering many incomplete or low-quality sequences that would be otherwise discarded. ImpuSARS can be incorporated in any primary data processing pipeline for SARS-CoV-2 whole-genome sequencing.

Phylogenetic Analysis of the 2020 West Nile Virus (WNV) Outbreak in Andalusia (Spain).

During recent decades West Nile Virus (WNV) outbreaks have continuously occurred in the Mediterranean area. In August 2020 a new WNV outbreak affected 71 people with meningoencephalitis in Andalusia and six more cases were detected in Extremadura (south-west of Spain), causing a total of eight deaths. The whole genomes of four viruses were obtained and phylogenetically analyzed in the context of recent outbreaks. The Andalusian viral samples belonged to lineage 1 and were relatively similar to those of previous outbreaks which occurred in the Mediterranean region. Here we present a detailed analysis of the outbreak, including an extensive phylogenetic study. As part on this effort, we implemented a local Nextstrain server, which has become a constituent piece of regional epidemiological surveillance, wherein forthcoming genomes of environmental samples or, eventually, future outbreaks, will be included.

Taxonomic variations in the gut microbiome of gout patients with and without tophi might have a functional impact on urate metabolism.

OBJECTIVE: To evaluate the taxonomic composition of the gut microbiome in gout patients with and without tophi formation, and predict bacterial functions that might have an impact on urate metabolism. METHODS: Hypervariable V3-V4 regions of the bacterial 16S rRNA gene from fecal samples of gout patients with and without tophi (n = 33 and n = 25, respectively) were sequenced and compared to fecal samples from 53 healthy controls. We explored predictive functional profiles using bioinformatics in order to identify differences in taxonomy and metabolic pathways. RESULTS: We identified a microbiome characterized by the lowest richness and a higher abundance of Phascolarctobacterium, Bacteroides, Akkermansia, and Ruminococcus_gnavus_group genera in patients with gout without tophi when compared to controls. The Proteobacteria phylum and the Escherichia-Shigella genus were more abundant in patients with tophaceous gout than in controls. Fold change analysis detected nine genera enriched in healthy controls compared to gout groups (Bifidobacterium, Butyricicoccus, Oscillobacter, Ruminococcaceae_UCG_010, Lachnospiraceae_ND2007_group, Haemophilus, Ruminococcus_1, Clostridium_sensu_stricto_1, and Ruminococcaceae_UGC_013). We found that the core microbiota of both gout groups shared Bacteroides caccae, Bacteroides stercoris ATCC 43183, and Bacteroides coprocola DSM 17136. These bacteria might perform functions linked to one-carbon metabolism, nucleotide binding, amino acid biosynthesis, and purine biosynthesis. Finally, we observed differences in key bacterial enzymes involved in urate synthesis, degradation, and elimination. CONCLUSION: Our findings revealed that taxonomic variations in the gut microbiome of gout patients with and without tophi might have a functional impact on urate metabolism.

Genome-wide analysis of DNA methylation in Hirschsprung enteric precursor cells: unraveling the epigenetic landscape of enteric nervous system development.

BACKGROUND: Hirschsprung disease (HSCR, OMIM 142623) is a rare congenital disorder that results from a failure to fully colonize the gut by enteric precursor cells (EPCs) derived from the neural crest. Such incomplete gut colonization is due to alterations in EPCs proliferation, survival, migration and/or differentiation during enteric nervous system (ENS) development. This complex process is regulated by a network of signaling pathways that is orchestrated by genetic and epigenetic factors, and therefore alterations at these levels can lead to the onset of neurocristopathies such as HSCR. The goal of this study is to broaden our knowledge of the role of epigenetic mechanisms in the disease context, specifically in DNA methylation. Therefore, with this aim, a Whole-Genome Bisulfite Sequencing assay has been performed using EPCs from HSCR patients and human controls. RESULTS: This is the first study to present a whole genome DNA methylation profile in HSCR and reveal a decrease of global DNA methylation in CpG context in HSCR patients compared with controls, which correlates with a greater hypomethylation of the differentially methylated regions (DMRs) identified. These results agree with the de novo Methyltransferase 3b downregulation in EPCs from HSCR patients compared to controls, and with the decrease in the global DNA methylation level previously described by our group. Through the comparative analysis of DMRs between HSCR patients and controls, a set of new genes has been identified as potential susceptibility genes for HSCR at an epigenetic level. Moreover, previous differentially methylated genes related to HSCR have been found, which validates our approach. CONCLUSIONS: This study highlights the relevance of an adequate methylation pattern for a proper ENS development. This is a research area that provides a novel approach to deepen our understanding of the etiopathogenesis of HSCR.

Virulence role of the outer membrane protein CarO in carbapenem-resistant Acinetobacter baumannii.

Novel approaches to treat carbapenem-resistant Acinetobacter baumannii (CRAB) infections are urgently needed and anti-virulence drugs represent promising alternatives, but our knowledge on potential targets is scarce. We searched for potential A. baumannii virulence factors by whole-genome sequencing-based comparisons of CRAB clinical isolates causing bloodstream infections secondary to ventilator-associated pneumonia from demographics and clinically homogeneous patients, who received optimal treatment but with different clinical outcomes. Thus, the carO gene was interrupted in CRAB isolates from surviving patients, while it was intact in isolates from non-surviving patients, and proteomic/immunoblot techniques corroborated it. When the virulence role of A. baumannii CarO was analyzed in model systems, isogenic ΔcarO mutants and a CRAB clinical isolate with truncated CarO, showed lower ability to adhere and invade A549 cells and in vivo virulence. This unnoticed virulence role for CarO postulate this A. baumannii outer membrane protein as a potential target for new therapies against CRAB infections.

Ancient evolutionary signals of protein-coding sequences allow the discovery of new genes in the Drosophila melanogaster genome.

BACKGROUND: The current growth in DNA sequencing techniques makes of genome annotation a crucial task in the genomic era. Traditional gene finders focus on protein-coding sequences, but they are far from being exhaustive. The number of this kind of genes continuously increases due to new experimental data and development of improved bioinformatics algorithms. RESULTS: In this context, AnABlast represents a novel in silico strategy, based on the accumulation of short evolutionary signals identified by protein sequence alignments of low score. This strategy potentially highlights protein-coding regions in genomic sequences regardless of traditional homology or translation signatures. Here, we analyze the evolutionary information that the accumulation of these short signals encloses. Using the Drosophila melanogaster genome, we stablish optimal parameters for the accurate gene prediction with AnABlast and show that this new strategy significantly contributes to add genes, exons and pseudogenes regions, yet to be discovered in both already annotated and new genomes. CONCLUSIONS: AnABlast can be freely used to analyze genomic regions of whole genomes where it contributes to complete the previous annotation.

Antibiotic resistance and metabolic profiles as functional biomarkers that accurately predict the geographic origin of city metagenomics samples.

BACKGROUND: The availability of hundreds of city microbiome profiles allows the development of increasingly accurate predictors of the origin of a sample based on its microbiota composition. Typical microbiome studies involve the analysis of bacterial abundance profiles. RESULTS: Here we use a transformation of the conventional bacterial strain or gene abundance profiles to functional profiles that account for bacterial metabolism and other cell functionalities. These profiles are used as features for city classification in a machine learning algorithm that allows the extraction of the most relevant features for the classification. CONCLUSIONS: We demonstrate here that the use of functional profiles not only predict accurately the most likely origin of a sample but also to provide an interesting functional point of view of the biogeography of the microbiota. Interestingly, we show how cities can be classified based on the observed profile of antibiotic resistances. REVIEWERS: Open peer review: Reviewed by Jin Zhuang Dou, Jing Zhou, Torsten Semmler and Eran Elhaik.

AnABlast: Re-searching for Protein-Coding Sequences in Genomic Regions.

AnABlast is a computational tool that highlights protein-coding regions within intergenic and intronic DNA sequences which escape detection by standard gene prediction algorithms. DNA sequences with small protein-coding genes or exons, complex intron-containing genes, or degenerated DNA fragments are efficiently targeted by AnABlast. Furthermore, this algorithm is particularly useful in detecting protein-coding sequences with nonsignificant homologs to sequences in databases. AnABlast can be executed online at http://www.bioinfocabd.upo.es/anablast/ .