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.This work is supported by grant PT17/0009/0006 from the Spanish Ministry of Economy and Competitiveness, COVID-0012–2020 from Consejería de Salud y Familias, Junta de Andalucía, and postdoctoral contract PAIDI2020- DOC_00350 for C.L., from Junta de Andalucía, co-funded by the European Social Fund (FSE) 2014–2020.Peer reviewe

Insight into potential probiotic markers predicted in Lactobacillus pentosus MP-10 genome sequence

Lactobacillus pentosus MP-10 is a potential probiotic lactic acid bacterium (LAB) originally isolated from naturally fermented Aloreña green table olives. The entire genome sequence was annotated to in-silico analyze the molecular mechanisms involved in the adaptation of L. pentosus MP-10 to the human gastrointestinal tract (GIT), such as carbohydrate metabolism (related with prebiotic utilization) and the proteins involved in bacteria-host interactions. We predicted an arsenal of genes coding for carbohydrate-modifying enzymes to modify oligo- and polysaccharides, such as glycoside hydrolases, glycoside transferases and isomerases, and other enzymes involved in complex carbohydrate metabolism especially starch, raffinose and levan. These enzymes represent key indicators of the bacteria’s adaptation to the GIT environment, since they involve the metabolism and assimilation of complex carbohydrates not digested by human enzymes. We also detected key probiotic ligands (surface proteins, excreted or secreted proteins) involved in the adhesion to host cells such as adhesion to mucus, epithelial cells or extracellular matrix, and plasma components; also, moonlighting proteins or multifunctional proteins were found that could be involved in adhesion to epithelial cells and/or extracellular matrix proteins and also affect host immunomodulation. In-silico analysis of the genome sequence of L. pentosus MP-10 is an important initial step to screen for genes encoding for proteins that may provide probiotic features, and thus provides one new routes for screening and studying this potentially probiotic bacterium

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 bla gene plays a key role in the P/T-R acquisition, no clinical in vivo study has yet confirmed the role of bla 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 bla and other genes associated with P/T-R. Changes in the outer membrane protein (OMP) profile were analyzed using SDS-PAGE, and bla 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 bla gene in P/T-R isolates, generated by three different genetic events: (1) IS26-mediated duplication of the bla gene, (2) generation of a small multicopy plasmid (ColE-like) carrying bla, 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.This work was funded by the Instituto de Salud Carlos III, Subdirección General de Redes y Centros de Investigación Cooperativa, Ministerio de Economía, Industria y Competitividad (grant PI19/01009), by Plan Nacional de I+D+i 2013–2016, and by the Instituto de Salud Carlos III, Subdirección General de Redes y Centros de Investigación Cooperativa, Ministerio de Economía, Industria y Competitividad, Spanish Network for Research in Infectious Diseases (grant RD16/0016/0009), cofinanced by the European Development Regional Fund “A way to make Europe”/”Investing in your future”. A.R.V, R.A-M, J.A.L. and J.M.C. were supported (grant CB21/13/00006) by CIBERINFEC - Consorcio Centro de Investigación Biomédica en Red, Instituto de Salud Carlos III, Ministerio de Ciencia e Innovación and Unión Europea – NextGenerationEU. L.G.B. was partly financed by a grant from the Spanish Society of Infectious Diseases and Clinical Microbiology (SEIMC) and is supported by the Subprograma Rio Hortega, Instituto Carlos III, Subdirección General de Redes y Centros de Investigación Cooperativa, Ministerio de Ciencia, Innovación y Universidades, Spain (CM22/00196). J.M.O.R. is supported by the Subprograma Sara Borrell, Instituto de Salud Carlos III, Subdirección General de Redes y Centros de Investigación Cooperativa, Ministerio de Ciencia, Innovación y Universidades, Spain (CD21/00098). A.R.V. is supported by the Subprograma Juan Rodés, Instituto de Salud Carlos III, Subdirección General de Redes y Centros de Investigación Cooperativa, Ministerio de Ciencia, Innovación y Universidades, Spain (JR20/00023), and Y.S. has received a Miguel Servet Tipo II contract from the Instituto de Salud Carlos III, Subdirección General de Redes y Centros de Investigación Cooperativa, Ministerio de Economía y Competitividad, Spain (grant CPII20/00018). Funding for open access publishing: Universidad Pablo de OlavidePeer reviewe

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.This study was supported by the Instituto Nacional de Rehabilitación Luis Guillermo Ibarra Ibarra” and the Grant INF-2016-01-269675 from the Consejo Nacional de Ciencia y Tecnología (CONACYT)

SARS-CoV-2 genome sequencing in Andalusia, methodology and study of variants

La incorporación de las técnicas de secuenciación genómica mediante secuenciación de nueva generación ha revolucionado la microbiología clínica, innovando y mejorando el diagnóstico clínico de las enfermedades infecciosas. Hoy en día, la secuenciación de genoma completo en enfermedades infecciosas tiene multitud de aplicaciones en virología, bacteriología, resistencia antibiótica, epidemiología y salud pública. Con la aparición del SARS-CoV-2, se ha visto subrayada la importancia del análisis y estudio de las secuencias genéticas. Desde la identificación inicial del SARS-CoV-2, hasta la fecha, se han compartido, a nivel mundial, más de 414.575 secuencias genómicas completas a través de bases de datos de acceso público. La capacidad de monitorizar la evolución viral casi en tiempo real tiene un impacto directo en la respuesta de salud pública a la pandemia de COVID-19. En este trabajo se presenta la importancia de la secuenciación genómica en microbiología, enfermedades infecciosas, epidemiología y salud pública, y se describe cómo se ha implementado la secuenciación de SARS-CoV-2 en Andalucía, y cuales son los principales resultados hasta la fecha.The incorporation of genomic sequencing techniques through next-generation sequencing has revolutionized clinical microbiology, innovating and improving the clinical diagnosis of infectious diseases. Today, whole genome sequencing in infectious diseases has many applications in virology, bacteriology, antibiotic resistance, epidemiology, and public health. With the appearance of SARS-CoV-2, the importance of the analysis and study of genetic sequences has been underlined. Since the initial identification of SARS-CoV-2, to date, more than 414,575 complete genomic sequences have been shared worldwide through public access databases. The ability to monitor viral evolution in near real time has a direct impact on the public health response to the COVID-19 pandemic. This paper presents the importance of genomic sequencing in microbiology, infectious diseases, epidemiology and public health, and describes how SARS-CoV-2 sequencing has been implemented in Andalusia, and what the main results are to date

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

Until the May 2022 Monkeypox 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 Monkeypox outbreak in the South of Spain. 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. Whole genome sequencing of a total of 160 monkeypox viruses from the different provinces that reported cases were obtained. Interestingly, we report the sequences of monkeypox 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 monkeypox 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 prior to treatmentThis work was supported by Spanish Ministry of Science and Innovation (grants PID2020- 117979RB-I00 and FJC2021-046546-I), the Instituto de Salud Carlos III (ISCIII), co-funded with European Regional Development Funds (ERDF) (grant IMP/00019), it has also been funded by Consejería de Salud y Consumo, Junta de Andalucía (grants COVID-0012-2020), and by grant ELIXIR-CONVERGE - Connect and align ELIXIR Nodes to deliver sustainable FAIR lifescience data management services (AMD-871075-16), funded by EU – H2020.N

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.This research was funded by Spanish Ministry of Science and Innovation (grant PID2020-117979RB-I00), the Instituto de Salud Carlos III (ISCIII), co-funded with European Regional Development Funds (ERDF) (grant IMP/00019), and has also been funded by Consejería de Salud y Familias, Junta de Andalucía (grants COVID-0012-2020, PS-2020-342 and IE19_259 FPS).Peer reviewe

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.This work was supported by Spanish Ministry of Science and Innovation (grant PID2020- 117979RB-I00), the Instituto de Salud Carlos III (ISCIII), co-funded with European Regional Development Funds (ERDF) (grant IMP/00019), and has also been funded by Consejería de Salud y Familias, Junta de Andalucía (grants COVID-0012-2020 and PS-2020-342) and the postdoctoral contract of Carlos Loucera (PAIDI2020- DOC_00350), co-funded by the European Social Fund (FSE) 2014-2020. ELIXIR-CONVERGE—H2020 (871075).Peer reviewe

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

PlantFuncSSR: Integrating first and next generation transcriptomics for mining of SSR-functional domains markers

et al.Analysis of repetitive DNA sequence content and divergence among the repetitive functional classes is a well-accepted approach for estimation of inter- and intra-generic differences in plant genomes. Among these elements, microsatellites, or Simple Sequence Repeats (SSRs), have been widely demonstrated as powerful genetic markers for species and varieties discrimination. We present PlantFuncSSRs platform having more than 364 plant species with more than 2 million functional SSRs. They are provided with detailed annotations for easy functional browsing of SSRs and with information on primer pairs and associated functional domains. PlantFuncSSRs can be leveraged to identify functional-based genic variability among the species of interest, which might be of particular interest in developing functional markers in plants. This comprehensive on-line portal unifies mining of SSRs from first and next generation sequencing datasets, corresponding primer pairs and associated in-depth functional annotation such as gene ontology annotation, gene interactions and its identification from reference protein databases.Gaurav Sablok thanks Plant Functional Biology and Climate Change Cluster (C3), University of Technology Sydney, PO Box 123, Broadway, NSW 2007, Australia, for providing the computational facilities. An internal grant number to GS (2226018) supported this work. JAH and TYS were partially supported by High Impact Research Chancellory Grant UM.C/625/1/HIR/MOHE/SCI/19 from the University of Malaya.Peer reviewedPeer Reviewe