Recurrence, Microevolution, and Spatiotemporal Dynamics of Legionella pneumophila Sequence Type 1905, Portugal, 2014-2022

We investigated molecular evolution and spatiotemporal dynamics of atypical Legionella pneumophila serogroup 1 sequence type 1905 and determined its long-term persistence and linkage to human disease in dispersed locations, far beyond the large 2014 outbreak epicenter in Portugal. Our finding highlights the need for public health interventions to prevent further disease spread.

INSaFLU-TELEVIR: an open web-based bioinformatics suite for viral metagenomic detection and routine genomic surveillance

BACKGROUND

Implementation of clinical metagenomics and pathogen genomic surveillance can be particularly challenging due to the lack of bioinformatics tools and/or expertise. In order to face this challenge, we have previously developed INSaFLU, a free web-based bioinformatics platform for virus next-generation sequencing data analysis. Here, we considerably expanded its genomic surveillance component and developed a new module (TELEVIR) for metagenomic virus identification.

RESULTS

The routine genomic surveillance component was strengthened with new workflows and functionalities, including (i) a reference-based genome assembly pipeline for Oxford Nanopore technologies (ONT) data; (ii) automated SARS-CoV-2 lineage classification; (iii) Nextclade analysis; (iv) Nextstrain phylogeographic and temporal analysis (SARS-CoV-2, human and avian influenza, monkeypox, respiratory syncytial virus (RSV A/B), as well as a "generic" build for other viruses); and (v) algn2pheno for screening mutations of interest. Both INSaFLU pipelines for reference-based consensus generation (Illumina and ONT) were benchmarked against commonly used command line bioinformatics workflows for SARS-CoV-2, and an INSaFLU snakemake version was released. In parallel, a new module (TELEVIR) for virus detection was developed, after extensive benchmarking of state-of-the-art metagenomics software and following up-to-date recommendations and practices in the field. TELEVIR allows running complex workflows, covering several combinations of steps (e.g., with/without viral enrichment or host depletion), classification software (e.g., Kaiju, Kraken2, Centrifuge, FastViromeExplorer), and databases (RefSeq viral genome, Virosaurus, etc.), while culminating in user- and diagnosis-oriented reports. Finally, to potentiate real-time virus detection during ONT runs, we developed findONTime, a tool aimed at reducing costs and the time between sample reception and diagnosis.

CONCLUSIONS

The accessibility, versatility, and functionality of INSaFLU-TELEVIR are expected to supply public and animal health laboratories and researchers with a user-oriented and pan-viral bioinformatics framework that promotes a strengthened and timely viral metagenomic detection and routine genomics surveillance. INSaFLU-TELEVIR is compatible with Illumina, Ion Torrent, and ONT data and is freely available at https://insaflu.insa.pt/ (online tool) and https://github.com/INSaFLU (code).

Raw milk cheeses from Beira Baixa, Portugal-A contributive study for the microbiological hygiene and safety assessment

Due to specific bacterial microbiota, raw milk cheeses have appreciated sensory properties. However, they may pose a threat to consumer safety due to potential pathogens presence. This study evaluated the microbiological contamination of 98 raw milk cheeses from Beira Baixa, Portugal. Presence and enumeration of Coagulase Positive Staphylococci (CPS), Listeria monocytogenes, Salmonella spp., pathogenic Escherichia coli, and indicator microorganisms (non-pathogenic E. coli and Listeria spp.) was attained. E. coli antimicrobial resistance (AMR) was also evaluated. PCR and/or Whole genome sequencing (WGS) was used to characterize E. coli, Salmonella spp. and L. monocytogenes isolates. Sixteen cheeses (16.3%) were classified as Satisfactory, 59 (60.2%) as Borderline and 23 (23.5%) as Unsatisfactory/Potential Injurious to Health. L. monocytogenes, CPS > 104 cfu g-1, Extraintestinal pathogenic E. coli (ExPEC) and Salmonella spp. were detected in 4.1%, 6.1%, 3.1% and 1.0% of the samples, respectively. Listeria innocua (4.1%) and E. coli > 104 cfu g-1 (16.3%) were also detected. AMR E. coli was detected in 23/98 (23.5%) of the cheese samples, of which two were multidrug resistant. WGS identified genotypes already associated to human disease and Listeria spp. cluster analysis indicated that cheese contamination might be related with noncompliance with Good Hygiene Practices during cheese production.

Unveiling a Listeria monocytogenes Outbreak in a Rabbit Farm: Clinical Manifestation, Antimicrobial Resistance, Genomic Insights and Environmental Investigation

Listeria monocytogenes poses a threat to both human and animal health. This work describes an L. monocytogenes outbreak in a Portuguese rabbit farm, detailing the isolates' clinical manifestations, necropsy findings, and phenotypic and genomic profiles. Clinical signs, exclusively observed in does, included lethargy and reproductive signs. Post-mortem examination of does revealed splenomegaly, hepatomegaly with a reticular pattern, pulmonary congestion, and haemorrhagic lesions in the uterus, with thickening of the uterine wall and purulent greyish exudates. Positive L. monocytogenes samples were identified in fattening and maternity units across different samples, encompassing does and environmental samples. Core-genome Multi Locus Sequence Typing (cgMLST) analysis confirmed the outbreak, with the 16 sequenced isolates (lineage II, CC31, and ST325) clustering within a ≤2 allelic difference (AD) threshold. Antimicrobial susceptibility testing for five antibiotics revealed that 15 out of 19 outbreak isolates were resistant to sulfamethoxazole-trimethoprim (SXT). Concordantly, all SXT-resistant sequenced isolates were found to exclusively harbour a plasmid containing a trimethoprim-resistance gene (dfrD), along with loci linked to resistance to lincosamides (lnuG), macrolides (mphB), and polyether ionophores (NarAB operon). All sequenced outbreak isolates carried the antibiotic resistance-related genes tetM, fosX, lin, norB, lmrB, sul, and mprF. The outbreak cluster comprises isolates from does and the environment, which underscores the ubiquitous presence of L. monocytogenes and emphasizes the importance of biosecurity measures. Despite limited data on listeriosis in rabbit farming, this outbreak reveals its significant impact on animal welfare and production.

Strengthening the genomic surveillance of Francisella tularensis by using culture-free whole-genome sequencing from biological samples

Introduction

Francisella tularensis is a highly infectious bacterium that causes the zoonotic disease tularemia. The development of genotyping methods, especially those based on whole-genome sequencing (WGS), has recently increased the knowledge on the epidemiology of this disease. However, due to the difficulties associated with the growth and isolation of this fastidious pathogen in culture, the availability of strains and subsequently WGS data is still limited.

Methods

To surpass these constraints, we aimed to implement a culture-free approach to capture and sequence F. tularensis genomes directly from complex samples. Biological samples obtained from 50 common voles and 13 Iberian hares collected in Spain were confirmed as positive for F. tularensis subsp. holarctica and subjected to a WGS target capture and enrichment protocol, using RNA oligonucleotide baits designed to cover F. tularensis genomic diversity.

Results

We obtained full genome sequences of F. tularensis from 13 animals (20.6%), two of which had mixed infections with distinct genotypes, and achieved a higher success rate when compared with culture-dependent WGS (only successful for two animals). The new genomes belonged to different clades commonly identified in Europe (B.49, B.51 and B.262) and subclades. Despite being phylogenetically closely related to other genomes from Spain, the detected clusters were often found in other countries. A comprehensive phylogenetic analysis, integrating 599 F. tularensis subsp. holarctica genomes, showed that most (sub)clades are found in both humans and animals and that closely related strains are found in different, and often geographically distant, countries.

Discussion

Overall, we show that the implemented culture-free WGS methodology yields timely, complete and high-quality genomic data of F. tularensis, being a highly valuable approach to promote and potentiate the genomic surveillance of F. tularensis and ultimately increase the knowledge on the genomics, ecology and epidemiology of this highly infectious pathogen.

Viral genetic clustering and transmission dynamics of the 2022 mpox outbreak in Portugal

Pathogen genome sequencing during epidemics enhances our ability to identify and understand suspected clusters and investigate their relationships. Here, we combine genomic and epidemiological data of the 2022 mpox outbreak to better understand early viral spread, diversification and transmission dynamics. By sequencing 52% of the confirmed cases in Portugal, we identified the mpox virus sublineages with the highest impact on case numbers and fitted them into a global context, finding evidence that several international sublineages probably emerged or spread early in Portugal. We estimated a 62% infection reporting rate and that 1.3% of the population of men who have sex with men in Portugal were infected. We infer the critical role played by sexual networks and superspreader gatherings, such as sauna attendance, in the dissemination of mpox virus. Overall, our findings highlight genomic epidemiology as a tool for the real-time monitoring and control of mpox epidemics, and can guide future vaccine policy in a highly susceptible population.

Recurrent Campylobacter jejuni Infections with In Vivo Selection of Resistance to Macrolides and Carbapenems: Molecular Characterization of Resistance Determinants

We present two independent cases of recurrent multidrug-resistant Campylobacter jejuni infection in immunocompromised hosts and the clinical challenges encountered due to the development of high-level carbapenem resistance. The mechanisms associated with this unusual resistance for Campylobacters were characterized. Initial macrolide and carbapenem-susceptible strains acquired resistance to erythromycin (MIC > 256mg/L), ertapenem (MIC > 32mg/L), and meropenem (MIC > 32mg/L) during treatment. Carbapenem-resistant isolates developed an in-frame insertion resulting in an extra Asp residue in the major outer membrane protein PorA, within the extracellular loop L3 that connects β-strands 5 and 6 and forms a constriction zone involved in Ca2+ binding. The isolates presenting the highest MIC to ertapenem exhibited an extra nonsynonymous mutation (G167A|Gly56Asp) at PorA's extracellular loop L1. IMPORTANCE Carbapenem susceptibility patterns suggest drug impermeability, related to either insertion and/or single nucleotide polymorphism (SNP) within porA. Similar molecular events occurring in two independent cases support the association of these mechanisms with carbapenem resistance in Campylobacter spp.

ReporTree: a surveillance-oriented tool to strengthen the linkage between pathogen genetic clusters and epidemiological data

BACKGROUND

Genomics-informed pathogen surveillance strengthens public health decision-making, playing an important role in infectious diseases' prevention and control. A pivotal outcome of genomics surveillance is the identification of pathogen genetic clusters and their characterization in terms of geotemporal spread or linkage to clinical and demographic data. This task often consists of the visual exploration of (large) phylogenetic trees and associated metadata, being time-consuming and difficult to reproduce.

RESULTS

We developed ReporTree, a flexible bioinformatics pipeline that allows diving into the complexity of pathogen diversity to rapidly identify genetic clusters at any (or all) distance threshold(s) or cluster stability regions and to generate surveillance-oriented reports based on the available metadata, such as timespan, geography, or vaccination/clinical status. ReporTree is able to maintain cluster nomenclature in subsequent analyses and to generate a nomenclature code combining cluster information at different hierarchical levels, thus facilitating the active surveillance of clusters of interest. By handling several input formats and clustering methods, ReporTree is applicable to multiple pathogens, constituting a flexible resource that can be smoothly deployed in routine surveillance bioinformatics workflows with negligible computational and time costs. This is demonstrated through a comprehensive benchmarking of (i) the cg/wgMLST workflow with large datasets of four foodborne bacterial pathogens and (ii) the alignment-based SNP workflow with a large dataset of Mycobacterium tuberculosis. To further validate this tool, we reproduced a previous large-scale study on Neisseria gonorrhoeae, demonstrating how ReporTree is able to rapidly identify the main species genogroups and characterize them with key surveillance metadata, such as antibiotic resistance data. By providing examples for SARS-CoV-2 and the foodborne bacterial pathogen Listeria monocytogenes, we show how this tool is currently a useful asset in genomics-informed routine surveillance and outbreak detection of a wide variety of species.

CONCLUSIONS

In summary, ReporTree is a pan-pathogen tool for automated and reproducible identification and characterization of genetic clusters that contributes to a sustainable and efficient public health genomics-informed pathogen surveillance. ReporTree is implemented in python 3.8 and is freely available at https://github.com/insapathogenomics/ReporTree .

Development of an amplicon-based sequencing approach in response to the global emergence of mpox

The 2022 multicountry mpox outbreak concurrent with the ongoing Coronavirus Disease 2019 (COVID-19) pandemic further highlighted the need for genomic surveillance and rapid pathogen whole-genome sequencing. While metagenomic sequencing approaches have been used to sequence many of the early mpox infections, these methods are resource intensive and require samples with high viral DNA concentrations. Given the atypical clinical presentation of cases associated with the outbreak and uncertainty regarding viral load across both the course of infection and anatomical body sites, there was an urgent need for a more sensitive and broadly applicable sequencing approach. Highly multiplexed amplicon-based sequencing (PrimalSeq) was initially developed for sequencing of Zika virus, and later adapted as the main sequencing approach for Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2). Here, we used PrimalScheme to develop a primer scheme for human monkeypox virus that can be used with many sequencing and bioinformatics pipelines implemented in public health laboratories during the COVID-19 pandemic. We sequenced clinical specimens that tested presumptively positive for human monkeypox virus with amplicon-based and metagenomic sequencing approaches. We found notably higher genome coverage across the virus genome, with minimal amplicon drop-outs, in using the amplicon-based sequencing approach, particularly in higher PCR cycle threshold (Ct) (lower DNA titer) samples. Further testing demonstrated that Ct value correlated with the number of sequencing reads and influenced the percent genome coverage. To maximize genome coverage when resources are limited, we recommend selecting samples with a PCR Ct below 31 Ct and generating 1 million sequencing reads per sample. To support national and international public health genomic surveillance efforts, we sent out primer pool aliquots to 10 laboratories across the United States, United Kingdom, Brazil, and Portugal. These public health laboratories successfully implemented the human monkeypox virus primer scheme in various amplicon sequencing workflows and with different sample types across a range of Ct values. Thus, we show that amplicon-based sequencing can provide a rapidly deployable, cost-effective, and flexible approach to pathogen whole-genome sequencing in response to newly emerging pathogens. Importantly, through the implementation of our primer scheme into existing SARS-CoV-2 workflows and across a range of sample types and sequencing platforms, we further demonstrate the potential of this approach for rapid outbreak response.

Comparative complete scheme and booster effectiveness of COVID-19 vaccines in preventing SARS-CoV-2 infections with SARS-CoV-2 Omicron (BA.1) and Delta (B.1.617.2) variants: A case-case study based on electronic health records

Background

Information on vaccine effectiveness in a context of novel variants of concern (VOC) emergence is of key importance to inform public health policies. This study aimed to estimate a measure of comparative vaccine effectiveness between Omicron (BA.1) and Delta (B.1.617.2 and sub-lineages) VOC according to vaccination exposure (primary or booster).

Methods

We developed a case-case study using data on RT-PCR SARS-CoV-2-positive cases notified in Portugal during Weeks 49-51, 2021. To obtain measure of comparative vaccine effectiveness, we compared the odds of vaccination in Omicron cases versus Delta using logistic regression adjusted for age group, sex, region, week of diagnosis, and laboratory of origin.

Results

Higher odds of vaccination were observed in cases infected by Omicron VOC compared with Delta VOC cases for both complete primary vaccination (odds ratio [OR] = 2.1; 95% confidence interval [CI]: 1.8 to 2.4) and booster dose (OR = 5.2; 95% CI: 3.1 to 8.8), equivalent to reduction of vaccine effectiveness from 44.7% and 92.8%, observed against infection with Delta, to -6.0% (95% CI: 29.2% to 12.7%) and 62.7% (95% CI: 35.7% to 77.9%), observed against infection with Omicron, for complete primary vaccination and booster dose, respectively.

Conclusion

Consistent reduction in vaccine-induced protection against infection with Omicron was observed. Complete primary vaccination may not be protective against SARS-CoV-2 infection in regions where Omicron variant is dominant.

Comparative Effectiveness of COVID-19 Vaccines in Preventing Infections and Disease Progression from SARS-CoV-2 Omicron BA.5 and BA.2, Portugal

We estimated comparative primary and booster vaccine effectiveness (VE) of SARS-CoV-2 Omicron BA.5 and BA.2 lineages against infection and disease progression. During April-June 2022, we implemented a case-case and cohort study and classified lineages using whole-genome sequencing or spike gene target failure. For the case-case study, we estimated the adjusted odds ratios (aORs) of vaccination using a logistic regression. For the cohort study, we estimated VE against disease progression using a penalized logistic regression. We observed no reduced VE for primary (aOR 1.07 [95% CI 0.93-1.23]) or booster (aOR 0.96 [95% CI 0.84-1.09]) vaccination against BA.5 infection. Among BA.5 case-patients, booster VE against progression to hospitalization was lower than that among BA.2 case-patients (VE 77% [95% CI 49%-90%] vs. VE 93% [95% CI 86%-97%]). Although booster vaccination is less effective against BA.5 than against BA.2, it offers substantial protection against progression from BA.5 infection to severe disease.

Molecular Capture of Mycobacterium tuberculosis Genomes Directly from Clinical Samples: A Potential Backup Approach for Epidemiological and Drug Susceptibility Inferences

The application of whole genome sequencing of Mycobacterium tuberculosis directly on clinical samples has been investigated as a means to avoid the time-consuming need for culture isolation that can lead to a potential prolonged suboptimal antibiotic treatment. We aimed to provide a proof-of-concept regarding the application of the molecular capture of M. tuberculosis genomes directly from positive sputum samples as an approach for epidemiological and drug susceptibility predictions. Smear-positive sputum samples (n = 100) were subjected to the SureSelectXT HS Target Enrichment protocol (Agilent Technologies, Santa Clara, CA, USA) and whole-genome sequencing analysis. A higher number of reads on target were obtained for higher smear grades samples (i.e., 3+ followed by 2+). Moreover, 37 out of 100 samples showed ≥90% of the reference genome covered with at least 10-fold depth of coverage (27, 9, and 1 samples were 3+, 2+, and 1+, respectively). Regarding drug-resistance/susceptibility prediction, for 42 samples, ≥90% of the >9000 hits that are surveyed by TB-profiler were detected. Our results demonstrated that M. tuberculosis genome capture and sequencing directly from clinical samples constitute a potential valid backup approach for phylogenetic inferences and resistance prediction, essentially in settings when culture is not routinely performed or for samples that fail to grow.

Development of an amplicon-based sequencing approach in response to the global emergence of human monkeypox virus

The 2022 multi-country monkeypox (mpox) outbreak concurrent with the ongoing COVID-19 pandemic has further highlighted the need for genomic surveillance and rapid pathogen whole genome sequencing. While metagenomic sequencing approaches have been used to sequence many of the early mpox infections, these methods are resource intensive and require samples with high viral DNA concentrations. Given the atypical clinical presentation of cases associated with the outbreak and uncertainty regarding viral load across both the course of infection and anatomical body sites, there was an urgent need for a more sensitive and broadly applicable sequencing approach. Highly multiplexed amplicon-based sequencing (PrimalSeq) was initially developed for sequencing of Zika virus, and later adapted as the main sequencing approach for SARS-CoV-2. Here, we used PrimalScheme to develop a primer scheme for human monkeypox virus that can be used with many sequencing and bioinformatics pipelines implemented in public health laboratories during the COVID-19 pandemic. We sequenced clinical samples that tested presumptive positive for human monkeypox virus with amplicon-based and metagenomic sequencing approaches. We found notably higher genome coverage across the virus genome, with minimal amplicon drop-outs, in using the amplicon-based sequencing approach, particularly in higher PCR cycle threshold (lower DNA titer) samples. Further testing demonstrated that Ct value correlated with the number of sequencing reads and influenced the percent genome coverage. To maximize genome coverage when resources are limited, we recommend selecting samples with a PCR cycle threshold below 31 Ct and generating 1 million sequencing reads per sample. To support national and international public health genomic surveillance efforts, we sent out primer pool aliquots to 10 laboratories across the United States, United Kingdom, Brazil, and Portugal. These public health laboratories successfully implemented the human monkeypox virus primer scheme in various amplicon sequencing workflows and with different sample types across a range of Ct values. Thus, we show that amplicon based sequencing can provide a rapidly deployable, cost-effective, and flexible approach to pathogen whole genome sequencing in response to newly emerging pathogens. Importantly, through the implementation of our primer scheme into existing SARS-CoV-2 workflows and across a range of sample types and sequencing platforms, we further demonstrate the potential of this approach for rapid outbreak response.

Pathogenic Escherichia coli, Salmonella spp. and Campylobacter spp. in Two Natural Conservation Centers of Wildlife in Portugal: Genotypic and Phenotypic Characterization

Human–wildlife coexistence may increase the potential risk of direct transmission of emergent or re-emergent zoonotic pathogens to humans. Intending to assess the occurrence of three important foodborne pathogens in wild animals of two wildlife conservation centers in Portugal, we investigated 132 fecal samples for the presence of Escherichia coli (Shiga toxin-producing E. coli (STEC) and non-STEC), Salmonella spp. and Campylobacter spp. A genotypic search for genes having virulence and antimicrobial resistance (AMR) was performed by means of PCR and Whole-Genome Sequencing (WGS) and phenotypic (serotyping and AMR profiles) characterization. Overall, 62 samples tested positive for at least one of these species: 27.3% for STEC, 11.4% for non-STEC, 3.0% for Salmonella spp. and 6.8% for Campylobacter spp. AMR was detected in four E. coli isolates and the only Campylobacter coli isolated in this study. WGS analysis revealed that 57.7% (30/52) of pathogenic E. coli integrated genetic clusters of highly closely related isolates (often involving different animal species), supporting the circulation and transmission of different pathogenic E. coli strains in the studied areas. These results support the idea that the health of humans, animals and ecosystems are interconnected, reinforcing the importance of a One Health approach to better monitor and control public health threats.

Antimicrobial Resistance and Biofilms Underlying Catheter-Related Bloodstream Coinfection by Enterobacter cloacae Complex and Candida parapsilosis

Biofilm-associated infections are a public health concern especially in the context of healthcare-associated infections such as catheter-related bloodstream infections (CRBSIs). We evaluated the biofilm formation and antimicrobials resistance (AMR) of Enterobacter cloacae complex and Candida parapsilosis co-isolated from a CRBSI patient. Antimicrobial susceptibility of central venous catheters (CVCs) and hemoculture (HC) isolates was evaluated, including whole genome sequencing (WGS) resistome analysis and evaluation of gene expression to obtain insight into their AMR determinants. Crystal violet assay was used to assess dual biofilm biomass and microscopy was used to elucidate a microorganism’s distribution within biofilms assembled on different materials. Bacteria were multidrug-resistant including resistance to colistin and beta-lactams, likely linked to the mcr-9-like phosphoethanolamine transferase and to an ACT family cephalosporin-hydrolyzing class C beta-lactamase, respectively. The R398I and Y132F mutations in the ERG11 gene and its differential expression might account for C. parapsilosis resistance to fluconazole. The phenotype of dual biofilms assembled on glass, polystyrene and polyurethane depends on the material and how biofilms were initiated by one or both pathogens. Biofilms assembled on polyurethane were denser and richer in the extracellular polymeric matrix, and microorganisms were differently distributed on the inner/outer surface of the CVC.

Bioinformatics study of expression from genomes of epidemiologically related MRSA CC398 isolates from human and wild animal samples

One of the most important livestock-associated methicillin-resistant Staphylococcus aureus (LA-MRSA) genetic lineages is the clonal complex (CC) 398, which can cause typical S. aureus-associated infections in people. In this work, whole-genome sequencing, RNA-sequencing, and gel-based comparative proteomics were applied to study the genetic characteristics of three MRSA CC398 isolates recovered from humans (strains C5621 and C9017), and from an animal (strain OR418). Of the three strains, C9017 presented the broadest resistance genotype, including resistance to fluroquinolone, clindamycin, tiamulin, macrolide and aminoglycoside antimicrobial classes. The scn, sak, and chp genes of the immune evasion cluster system were solely detected in OR418. Pangenome analysis showed a total of 288 strain-specific genes, most of which are hypothetical or phage-related proteins. OR418 had the most pronounced genetic differences. RNAIII (δ-hemolysin) gene was clearly the most expressed gene in OR418 and C5621, but it was not detected in C9017. Significant differences in the proteome profiles were found between strains. For example, the immunoglobulin-binding protein Sbi was more abundant in OR418. Considering that Sbi is a multifunctional immune evasion factor in S. aureus, the results point to OR418 strain having high zoonotic potential. Overall, multiomics biomarker signatures can assume an important role to advance precision medicine in the years to come. SIGNIFICANCE: MRSA is one of the most representative drug-resistant pathogens and its dissemination is increasing due to MRSA capability of establishing new reservoirs. LA-MRSA is considered an emerging problem worldwide and CC398 is one of the most important genetic lineages. In this study, three MRSA CC398 isolates recovered from humans and from a wild animal were analyzed through whole-genome sequencing, RNA-sequencing, and gel-based comparative proteomics in order to gather systems-wide omics data and better understand the genetic characteristics of this lineage to identify distinctive markers and genomic features of relevance to public health. The scn, sak, and chp genes of the immune evasion cluster system were solely detected in OR418. Pangenome analysis showed a total of 288 strain-specific genes, most of which are hypothetical or phage-related proteins. OR418 had the most pronounced genetic differences. RNAIII (δ-hemolysin) gene was clearly the most expressed gene in OR418 and C5621, but it was not detected in C9017. Significant differences in the proteome profiles were found between strains.

Urgent need for a non-discriminatory and non-stigmatizing nomenclature for monkeypox virus

The current nomenclature for monkeypox virus is stigmatising and misleading. This Perspective article proposes a practical and neutral system of nomenclature that will allow efficient communication without the risk of further misconceptions, discrimination and stigmatisation.

Phylogenomic characterization and signs of microevolution in the 2022 multi-country outbreak of monkeypox virus

The largest monkeypox virus (MPXV) outbreak described so far in non-endemic countries was identified in May 2022 (refs. 1-6). In this study, shotgun metagenomics allowed the rapid reconstruction and phylogenomic characterization of the first MPXV outbreak genome sequences, showing that this MPXV belongs to clade 3 and that the outbreak most likely has a single origin. Although 2022 MPXV (lineage B.1) clustered with 2018-2019 cases linked to an endemic country, it segregates in a divergent phylogenetic branch, likely reflecting continuous accelerated evolution. An in-depth mutational analysis suggests the action of host APOBEC3 in viral evolution as well as signs of potential MPXV human adaptation in ongoing microevolution. Our findings also indicate that genome sequencing may provide resolution to track the spread and transmission of this presumably slow-evolving double-stranded DNA virus.

Global gene expression analysis of Streptococcus agalactiae at exponential growth phase

A comparative overview of the global gene expression levels of S. agalactiae reference strain NEM316 at the exponential growth phase was done through RNA-sequencing. The expression levels of 47 genes potentially linked to virulence evidenced that: i) the major nuclease, GBS_RS03720/gbs0661, presented higher mean expression values than the remainder of DNase genes; ii) the genetic pilus island PI-2a genes presented higher mean expression values than PI-1 coding genes; and, iii) three virulence-associated genes ranked among the top-100 most expressed genes (GBS_RS07760, GBS_RS09445 and GBS_RS03485). Among this top-100, genes encoding proteins involved in "Translation, ribosomal structure and biogenesis" represented 46%. Curiously, genes with no assigned function were grouped in the category of highly expressed genes. As very little is known about the molecular mechanisms behind the release of DNases, preliminary assays were developed to understand whether direct DNA exposure would affect gene expression at the exponential growth phase. No differentially expressed genes were detected, indicating that follow-up studies are needed to disclose the complex molecular pathways (and stimuli) triggering the release of DNases. In general, our insights on the global expression levels of NEM316 at exponential growth phase with and without DNA exposure should open novel research lines to decipher S. agalactiae puzzling adaptation and virulence mechanisms, such as DNase production.

The Type III Secretion Effector CteG Mediates Host Cell Lytic Exit of Chlamydia trachomatis

Chlamydia trachomatis is an obligate intracellular bacterium causing ocular and urogenital infections in humans that are a significant burden worldwide. The completion of its characteristic infectious cycle relies on the manipulation of several host cell processes by numerous chlamydial type III secretion effector proteins. We previously identified the C. trachomatis CteG effector and showed it localizes at the host cell plasma membrane at late stages of infection. Here, we showed that, from 48 h post-infection, mammalian cells infected by wild-type C. trachomatis contained more infectious chlamydiae in the culture supernatant than cells infected by a CteG-deficient strain. This phenotype was CteG-dependent as it could be complemented in cells infected by the CteG-deficient strain carrying a plasmid encoding CteG. Furthermore, we detected a CteG-dependent defect on host cell cytotoxicity, indicating that CteG mediates chlamydial lytic exit. Previous studies showed that Pgp4, a global regulator of transcription encoded in the C. trachomatis virulence plasmid, also mediates chlamydial lytic exit. However, by using C. trachomatis strains encoding or lacking Pgp4, we showed that production and localization of CteG are not regulated by Pgp4. A C. trachomatis strain lacking both CteG and Pgp4 was as defective in promoting host cell cytotoxicity as mutant strains lacking only CteG or Pgp4. Furthermore, CteG overproduction in a plasmid suppressed the host cell cytotoxic defect of CteG- and Pgp4-deficient chlamydiae. Overall, we revealed the first chlamydial type III secretion effector involved in host cell lytic exit. Our data indicates that CteG and Pgp4 participate in a single cascade of events, but involving multiple layers of regulation, leading to lysis of host cells and release of the infectious chlamydiae.

A Search for Novel Legionella pneumophila Effector Proteins Reveals a Strain Specific Nucleotropic Effector

Legionella pneumophila is an accidental human pathogen that causes the potentially fatal Legionnaires’ disease, a severe type of pneumonia. The main virulence mechanism of L. pneumophila is a Type 4B Secretion System (T4SS) named Icm/Dot that transports effector proteins into the host cell cytosol. The concerted action of effectors on several host cell processes leads to the formation of an intracellular Legionella-containing vacuole that is replication competent and avoids phagolysosomal degradation. To date over 300 Icm/Dot substrates have been identified. In this study, we searched the genome of a L. pneumophila strain (Pt/VFX2014) responsible for the second largest L. pneumophila outbreak worldwide (in Vila Franca de Xira, Portugal, in 2014) for genes encoding potential novel Icm/Dot substrates. This strain Pt/VFX2014 belongs to serogroup 1 but phylogenetically segregates from all other serogroup 1 strains previously sequenced, displaying a unique mosaic genetic backbone. The ability of the selected putative effectors to be delivered into host cells by the T4SS was confirmed using the TEM-1 β-lactamase reporter assay. Two previously unknown Icm/Dot effectors were identified, VFX05045 and VFX10045, whose homologs Lpp1450 and Lpp3070 in clinical strain L. pneumophila Paris were also confirmed as T4SS substrates. After delivery into the host cell cytosol, homologs VFX05045/Lpp1450 remained diffused in the cell, similarly to Lpp3070. In contrast, VFX10045 localized to the host cell nucleus. To understand how VFX10045 and Lpp3070 (94% of identity at amino acid level) are directed to distinct sites, we carried out a comprehensive site-directed mutagenesis followed by analyses of the subcellular localization of the mutant proteins. This led to the delineation of region in the C-terminal part (residues 380 to 534) of the 583 amino acid-long VFX10045 as necessary and sufficient for nuclear targeting and highlighted the fundamental function of the VFX10045-specific R440 and I441 residues in this process. These studies revealed a strain-specific nucleotropism for new effector VFX10045/Lpp3070, which anticipates distinct functions between these homologs.

Mutation rate of SARS-CoV-2 and emergence of mutators during experimental evolution

Background and objectives

To understand how organisms evolve, it is fundamental to study how mutations emerge and establish. Here, we estimated the rate of mutation accumulation of SARS-CoV-2 in vitro and investigated the repeatability of its evolution when facing a new cell type but no immune or drug pressures.

Methodology

We performed experimental evolution with two strains of SARS-CoV-2, one carrying the originally described spike protein (CoV-2-D) and another carrying the D614G mutation that has spread worldwide (CoV-2-G). After 15 passages in Vero cells and whole genome sequencing, we characterized the spectrum and rate of the emerging mutations and looked for evidences of selection across the genomes of both strains.

Results

From the frequencies of the mutations accumulated, and excluding the genes with signals of selection, we estimate a spontaneous mutation rate of 1.3 × 10 ^-6 ± 0.2 × 10^-6 per-base per-infection cycle (mean across both lineages of SARS-CoV-2 ± 2SEM). We further show that mutation accumulation is larger in the CoV-2-D lineage and heterogeneous along the genome, consistent with the action of positive selection on the spike protein, which accumulated five times more mutations than the corresponding genomic average. We also observe the emergence of mutators in the CoV-2-G background, likely linked to mutations in the RNA-dependent RNA polymerase and/or in the error-correcting exonuclease protein.

Conclusions and implications

These results provide valuable information on how spontaneous mutations emerge in SARS-CoV-2 and on how selection can shape its genome toward adaptation to new environments. Lay Summary: Each time a virus replicates inside a cell, errors (mutations) occur. Here, via laboratory propagation in cells originally isolated from the kidney epithelium of African green monkeys, we estimated the rate at which the SARS-CoV-2 virus mutates-an important parameter for understanding how it can evolve within and across humans. We also confirm the potential of its Spike protein to adapt to a new environment and report the emergence of mutators-viral populations where mutations occur at a significantly faster rate.

Comparative Effectiveness of Coronavirus Vaccine in Preventing Breakthrough Infections among Vaccinated Persons Infected with Delta and Alpha Variants

We developed a case-case study to compare mRNA vaccine effectiveness against Delta versus Alpha coronavirus variants. We used data on 2,097 case-patients with PCR-positive severe acute respiratory syndrome coronavirus 2 infections reported in Portugal during May-July 2021. We estimated the odds of vaccine breakthrough infection in Delta-infected versus Alpha-infected patients by using conditional logistic regression adjusted for age group and sex and matched by the week of diagnosis. We compared reverse-transcription PCR cycle threshold values by vaccination status and variant as an indirect measure of viral load. We found significantly higher odds of vaccine breakthrough infection in Delta-infected patients than in Alpha-infected patients (odds ratio 1.96 [95% CI 1.22-3.14]), suggesting lower effectiveness of the mRNA vaccines in preventing infection with the Delta variant. We estimated lower mean cycle threshold values for the Delta cases (mean difference -2.10 [95% CI -2.74 to -1.47]), suggesting higher infectiousness than the Alpha variant.

Emergence and global spread of Listeria monocytogenes main clinical clonal complex

The bacterial foodborne pathogen Listeria monocytogenes clonal complex 1 (Lm-CC1) is the most prevalent clonal group associated with human listeriosis and is strongly associated with cattle and dairy products. Here, we analyze 2021 isolates collected from 40 countries, covering Lm-CC1 first isolation to present days, to define its evolutionary history and population dynamics. We show that Lm-CC1 spread worldwide from North America following the Industrial Revolution through two waves of expansion, coinciding with the transatlantic livestock trade in the second half of the 19th century and the rapid growth of cattle farming and food industrialization in the 20th century. In sharp contrast to its global spread over the past century, transmission chains are now mostly local, with limited inter- and intra-country spread. This study provides an unprecedented insight into L. monocytogenes phylogeography and population dynamics and highlights the importance of genome analyses for a better control of pathogen transmission.

Potential recurrence of COVID-19 in a healthcare professional: SARS-CoV-2 genome sequencing confirms contagiousness after re-positivity

Re-positivity of SARS-CoV-2 tests is widely reported, raising discussion about guidance for patient discharge and ending isolation. The unsuccessful recovery of replication-competent virus and/or absence of secondary cases has suggested that re-positive patients are not contagious. This study reports SARS-CoV-2 re-positivity in a healthcare professional 16 days after three negative tests, with viral genome sequencing supporting contagiousness leading to secondary cases.

Tracking the international spread of SARS-CoV-2 lineages B.1.1.7 and B.1.351/501Y-V2 with grinch

Late in 2020, two genetically-distinct clusters of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) with mutations of biological concern were reported, one in the United Kingdom and one in South Africa. Using a combination of data from routine surveillance, genomic sequencing and international travel we track the international dispersal of lineages B.1.1.7 and B.1.351 (variant 501Y-V2). We account for potential biases in genomic surveillance efforts by including passenger volumes from location of where the lineage was first reported, London and South Africa respectively. Using the software tool grinch (global report investigating novel coronavirus haplotypes), we track the international spread of lineages of concern with automated daily reports, Further, we have built a custom tracking website (cov-lineages.org/global_report.html) which hosts this daily report and will continue to include novel SARS-CoV-2 lineages of concern as they are detected.

Genomic insights on DNase production in Streptococcus agalactiae ST17 and ST19 strains

Streptococcus agalactiae evasion from the human defense mechanisms has been linked to the production of DNases. These were proposed to contribute to the hypervirulence of S. agalactiae ST17/capsular-type III strains, mostly associated with neonatal meningitis. We performed a comparative genomic analysis between ST17 and ST19 human strains with different cell tropism and distinct DNase production phenotypes. All S. agalactiae ST17 strains, with the exception of 2211-04, were found to display DNase activity, while the opposite scenario was observed for ST19, where 1203-05 was the only DNase(+) strain. The analysis of the genetic variability of the seven genes putatively encoding secreted DNases in S. agalactiae revealed an exclusive amino acid change in the predicted signal peptide of GBS0661 (NucA) of the ST17 DNase(-), and an exclusive amino acid change alteration in GBS0609 of the ST19 DNase(+) strain. Further core-genome analysis identified some specificities (SNVs or indels) differentiating the DNase(-) ST17 2211-04 and the DNase(+) ST19 1203-05 from the remaining strains of each ST. The pan-genomic analysis evidenced an intact phage without homology in S. agalactiae and a transposon homologous to TnGBS2.3 in ST17 DNase(-) 2211-04; the transposon was also found in one ST17 DNase(+) strain, yet with a different site of insertion. A group of nine accessory genes were identified among all ST17 DNase(+) strains, including the Eco47II family restriction endonuclease and the C-5 cytosine-specific DNA methylase. None of these loci was found in any DNase(-) strain, which may suggest that these proteins might contribute to the lack of DNase activity. In summary, we provide novel insights on the genetic diversity between DNase(+) and DNase(-) strains, and identified genetic traits, namely specific mutations affecting predicted DNases (NucA and GBS0609) and differences in the accessory genome, that need further investigation as they may justify distinct DNase-related virulence phenotypes in S. agalactiae.

Tracking the international spread of SARS-CoV-2 lineages B.1.1.7 and B.1.351/501Y-V2

Late in 2020, two genetically-distinct clusters of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) with mutations of biological concern were reported, one in the United Kingdom and one in South Africa. Using a combination of data from routine surveillance, genomic sequencing and international travel we track the international dispersal of lineages B.1.1.7 and B.1.351 (variant 501Y-V2). We account for potential biases in genomic surveillance efforts by including passenger volumes from location of where the lineage was first reported, London and South Africa respectively. Using the software tool grinch (global report investigating novel coronavirus haplotypes), we track the international spread of lineages of concern with automated daily reports, Further, we have built a custom tracking website (cov-lineages.org/global_report.html) which hosts this daily report and will continue to include novel SARS-CoV-2 lineages of concern as they are detected.

Characteristics of SARS-CoV-2 variants of concern B.1.1.7, B.1.351 or P.1: data from seven EU/EEA countries, weeks 38/2020 to 10/2021

We compared 19,207 cases of SARS-CoV-2 variant B.1.1.7/S gene target failure (SGTF), 436 B.1.351 and 352 P.1 to non-variant cases reported by seven European countries. COVID-19 cases with these variants had significantly higher adjusted odds ratios for hospitalisation (B.1.1.7/SGTF: 1.7, 95% confidence interval (CI): 1.0-2.9; B.1.351: 3.6, 95% CI: 2.1-6.2; P.1: 2.6, 95% CI: 1.4-4.8) and B.1.1.7/SGTF and P.1 cases also for intensive care admission (B.1.1.7/SGTF: 2.3, 95% CI: 1.4-3.5; P.1: 2.2, 95% CI: 1.7-2.8).

Nosocomial Outbreak of SARS-CoV-2 in a "Non-COVID-19" Hospital Ward: Virus Genome Sequencing as a Key Tool to Understand Cryptic Transmission

Dissemination of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in healthcare institutions affects both patients and health-care workers (HCW), as well as the institutional capacity to provide essential health services. Here, we investigated an outbreak of SARS-CoV-2 in a "non-COVID-19" hospital ward unveiled by massive testing, which challenged the reconstruction of transmission chains. The contacts network during the 15-day period before the screening was investigated, and positive SARS-CoV-2 RNA samples were subjected to virus genome sequencing. Of the 245 tested individuals, 48 (21 patients and 27 HCWs) tested positive for SARS-CoV-2. HCWs were mostly asymptomatic, but the mortality among patients reached 57.1% (12/21). Phylogenetic reconstruction revealed that all cases were part of the same transmission chain. By combining contact tracing and genomic data, including analysis of emerging minor variants, we unveiled a scenario of silent SARS-CoV-2 dissemination, mostly driven by the close contact within the HCWs group and between HCWs and patients. This investigation triggered enhanced prevention and control measures, leading to more timely detection and containment of novel outbreaks. This study shows the benefit of combining genomic and epidemiological data for disclosing complex nosocomial outbreaks, and provides valuable data to prevent transmission of COVID-19 in healthcare facilities.

Tracking SARS-CoV-2 lineage B.1.1.7 dissemination: insights from nationwide spike gene target failure (SGTF) and spike gene late detection (SGTL) data, Portugal, week 49 2020 to week 3 2021

We show that the SARS-CoV-2 B.1.1.7 lineage is highly disseminated in Portugal, with the odds of B.1.1.7 proportion increasing at an estimated 89% (95% confidence interval: 83-95%) per week until week 3 2021. RT-PCR spike gene target late detection (SGTL) can constitute a useful surrogate to track B.1.1.7 spread, besides the spike gene target failure (SGTF) proxy. SGTL/SGTF samples were associated with statistically significant higher viral loads, but not with substantial shift in age distribution compared to non-SGTF/SGTL cases.

Neisseria gonorrhoeae clustering to reveal major European whole-genome-sequencing-based genogroups in association with antimicrobial resistance

Neisseria gonorrhoeae, the bacterium responsible for the sexually transmitted disease gonorrhoea, has shown an extraordinary ability to develop antimicrobial resistance (AMR) to multiple classes of antimicrobials. With no available vaccine, managing N. gonorrhoeae infections demands effective preventive measures, antibiotic treatment and epidemiological surveillance. The latter two are progressively being supported by the generation of whole-genome sequencing (WGS) data on behalf of national and international surveillance programmes. In this context, this study aims to perform N. gonorrhoeae clustering into genogroups based on WGS data, for enhanced prospective laboratory surveillance. Particularly, it aims to identify the major circulating WGS-genogroups in Europe and to establish a relationship between these and AMR. Ultimately, it enriches public databases by contributing with WGS data from Portuguese isolates spanning 15 years of surveillance. A total of 3791 carefully inspected N. gonorrhoeae genomes from isolates collected across Europe were analysed using a gene-by-gene approach (i.e. using cgMLST). Analysis of cluster composition and stability allowed the classification of isolates into a two-step hierarchical genogroup level determined by two allelic distance thresholds revealing cluster stability. Genogroup clustering in general agreed with available N. gonorrhoeae typing methods [i.e. MLST (multilocus sequence typing), NG-MAST (N. gonorrhoeae multi-antigen sequence typing) and PubMLST core-genome groups], highlighting the predominant genogroups circulating in Europe, and revealed that the vast majority of the genogroups present a dominant AMR profile. Additionally, a non-static gene-by-gene approach combined with a more discriminatory threshold for potential epidemiological linkage enabled us to match data with previous reports on outbreaks or transmission chains. In conclusion, this genogroup assignment allows a comprehensive analysis of N. gonorrhoeae genetic diversity and the identification of the WGS-based genogroups circulating in Europe, while facilitating the assessment (and continuous monitoring) of their frequency, geographical dispersion and potential association with specific AMR signatures. This strategy may benefit public-health actions through the prioritization of genogroups to be controlled, the identification of emerging resistance carriage, and the potential facilitation of data sharing and communication.

Transcontinental Dissemination of the L2b/D-Da Recombinant Chlamydia trachomatis Lymphogranuloma venereum (LGV) Strain: Need of Broad Multi-Country Molecular Surveillance

Previously, we identified a Chlamydia trachomatis lymphogranuloma venereum (LGV) recombinant strain possessing a non-LGV ompA genotype. Here, culture-independent genome sequencing confirms its circulation in Europe, Middle East, and North America, and unveils emergence of antibiotic resistance. Broad surveillance is needed.

Massive dissemination of a SARS-CoV-2 Spike Y839 variant in Portugal

Genomic surveillance of SARS-CoV-2 was rapidly implemented in Portugal by the National Institute of Health in collaboration with a nationwide consortium of >50 hospitals/laboratories. Here, we track the geotemporal spread of a SARS-CoV-2 variant with a mutation (D839Y) in a potential host-interacting region involving the Spike fusion peptide, which is a target motif of anti-viral drugs that plays a key role in SARS-CoV-2 infectivity. The Spike Y839 variant was most likely imported from Italy in mid-late February and massively disseminated in Portugal during the early epidemic, becoming prevalent in the Northern and Central regions of Portugal where it represented 22% and 59% of the sampled genomes, respectively, by 30 April. Based on our high sequencing sampling during the early epidemics [15.5% (1275/8251) and 6.0% (1500/24987) of all confirmed cases until the end of March and April, respectively], we estimate that, between 14 March and 9 April (covering the epidemic exponential phase) the relative frequency of the Spike Y839 variant increased at a rate of 12.1% (6.1%-18.2%, CI 95%) every three days, being potentially associated with 24.8% (20.8-29.7%, CI 95%; 3177-4542 cases, CI 95%) of all COVID-19 cases in Portugal during this period. Our data supports population/epidemiological (founder) effects contributing to the Y839 variant superspread. The potential existence of selective advantage is also discussed, although experimental validation is required. Despite huge differences in genome sampling worldwide, SARS-CoV-2 Spike D839Y has been detected in 13 countries in four continents, supporting the need for close surveillance and functional assays of Spike variants.

Multiomics Substrates of Resistance to Emerging Pathogens? Transcriptome and Proteome Profile of a Vancomycin-Resistant Enterococcus faecalis Clinical Strain

Antibiotic resistance and hospital acquired infections are on the rise worldwide. Vancomycin-resistant enterococci have been reported in clinical settings in recent decades. In this multiomics study, we provide comprehensive proteomic and transcriptomic analyses of a vancomycin-resistant Enterococcus faecalis clinical isolate from a patient with a urinary tract infection. The previous genotypic profile of the strain C2620 indicated the presence of antibiotic resistance genes characteristic of the vanB cluster. To further investigate the transcriptome of this pathogenic strain, we used whole genome sequencing and RNA-sequencing to detect and quantify the genes expressed. In parallel, we used two-dimensional gel electrophoresis followed by MALDI-TOF/MS (Matrix-assisted laser desorption/ionization-Time-of-flight/Mass spectrometry) to identify the proteins in the proteome. We studied the membrane and cytoplasm subproteomes separately. From a total of 207 analysis spots, we identified 118 proteins. The protein list was compared to the results obtained from the full transcriptome assay. Several genes and proteins related to stress and cellular response were identified, as well as some linked to antibiotic and drug responses, which is consistent with the known state of multiresistance. Even though the correlation between transcriptome and proteome data is not yet fully understood, the use of multiomics approaches has proven to be increasingly relevant to achieve deeper insights into the survival ability of pathogenic bacteria found in health care facilities.

Mitogenome diversity of Aedes (Stegomyia) albopictus: Detection of multiple introduction events in Portugal

Aedes albopictus, along with Ae. aegypti, are key arbovirus vectors that have been expanding their geographic range over the last decades. In 2017, Ae. albopictus was detected for the first time at two distinct locations in Portugal. In order to understand how the Ae. albopictus populations recently introduced in Portugal are genetically related and which is their likely route of invasion, we performed an integrative cytochrome C oxidase I gene (COI)- and mitogenome-based phylogeographic analysis of mosquitoes samples collected in Portugal in 2017 and 2018 in the context of the global Ae. albopictus diversity. COI-based analysis (31 partial sequences obtained from 83 mosquitoes) revealed five haplotypes (1 to 5), with haplotype 1 (which is widely distributed in temperate areas worldwide) being detected in both locations. Haplotypes 2 and 3 were exclusively found in Southern region (Algarve), while haplotype 4 and 5 were only detected in the North of Portugal (Penafiel, Oporto region). Subsequent high discriminatory analyses based on Ae. albopictus mitogenome (17 novel sequences) not only confirmed a high degree of genetic variability within and between populations at both geographic locations (compatible with the Ae. albopictus mosquito populations circulating in Europe), but also revealed two mitogenome mutational signatures not previously reported at worldwide level. While our results generally sustain the occurrence of multiple introduction events, fine mitogenome sequence inspection further indicates a possible Ae. albopictus migration within the country, from the Northern introduction locality to the Southern region. In summary, the observed scenario of high Ae. albopictus genetic diversity in Portugal, together with the detection of mosquitoes in successive years since 2017 in Algarve and Penafiel, points that both Ae. albopictus populations seem to be already locally established, as its presence has been reported for three consecutive years, raising the public health awareness for future mosquito-borne diseases outbreaks.

In 2017, Aedes albopictus was reported for the first time in Portugal at two distinct locations, in the premises of a tire company in Penafiel, in the North, and nearby a golf course in Algarve, a tourism destination in the southernmost country region. The geographical spread of this species is promoted by egg and larvae transport in aquatic trade goods, as tires and aquatic plants, and adult anthropophilic behavior that favors passive land transportation. In Portugal, especially in the Southern region, temperate climate conditions are adequate for adult mosquitoes to survive most of the year. In a way to understand the genetic variability of Ae. albopictus populations introduced in Portugal, we analyzed 31 cytochrome C oxidase I gene (COI) partial sequences and 17 mitogenome sequences, integrating them in the context of the global Ae. albopictus phylogeographic diversity (i.e., 183 COI and 26 mitogenome sequences previously reported at worldwide level). Although COI haplotype 1 predominated, four additional haplotypes (2 to 5) were detected in Portugal. Subsequent in-depth mitogenome analysis revealed considerable genetic diversity, including not only sequences relating to mitogenomes reported mainly from Italy, Japan and China, but also two novel mitogenome mutational signatures.

Our study indicates that Ae. albopictus is locally established in Portugal and intra-country dispersal may have already happened, highlighting the challenges for vector surveillance and control programs aiming at restraining arbovirus disease burden in the future.

The anti-apoptotic Coxiella burnetii effector protein AnkG is a strain specific virulence factor

The ability to inhibit host cell apoptosis is important for the intracellular replication of the obligate intracellular pathogen Coxiella burnetii, as it allows the completion of the lengthy bacterial replication cycle. Effector proteins injected into the host cell by the C. burnetii type IVB secretion system (T4BSS) are required for the inhibition of host cell apoptosis. AnkG is one of these anti-apoptotic effector proteins. The inhibitory effect of AnkG requires its nuclear localization, which depends on p32-dependent intracellular trafficking and importin-α1-mediated nuclear entry of AnkG. Here, we compared the sequences of ankG from 37 C. burnetii isolates and classified them in three groups based on the predicted protein size. The comparison of the three different groups allowed us to identify the first 28 amino acids as essential and sufficient for the anti-apoptotic activity of AnkG. Importantly, only the full-length protein from the first group is a bona fide effector protein injected into host cells during infection and has anti-apoptotic activity. Finally, using the Galleria mellonella infection model, we observed that AnkG from the first group has the ability to attenuate pathology during in vivo infection, as it allows survival of the larvae despite bacterial replication.

MRSA CC398 recovered from wild boar harboring new SCCmec type IV J3 variant

A methicillin-resistant Staphylococcus aureus CC398 was recovered from a wild female boar (Sus scrofa) in the north of Portugal, in 2013 (Sousa et al. 2017). Whole genome sequencing (WGS) revealed this strain carries a new variant of a mecA-containing staphylococcal chromosomal gene cassette (SCCmec) type IV with an uncommon J3 region. WGS studies can facilitate surveillance and provide more detailed characterization of bacterial clones circulating in the wild, reinforcing the need for a one health perspective to better understand and control antimicrobial resistance.

Chlamydia trachomatis: when the virulence-associated genome backbone imports a prevalence-associated major antigen signature

Chlamydia trachomatis is the most prevalent sexually transmitted bacterium worldwide and the causative agent of trachoma. Its strains are classified according to their ompA genotypes, which are strongly linked to differential tissue tropism and disease outcomes [ocular disease, urogenital disease and lymphogranuloma venereum (LGV)]. While the genome-based species phylogenetic tree presents four main clades correlating with tropism/prevalence, namely ocular, LGV, urogenital T1 (more prevalent genotypes) and urogenital T2 (less prevalent genotypes), inter-clade exchange of ompA is considered a rare phenomenon probably mediating marked tropism alterations. An LGV epidemic, associated with the clonal expansion of the L2b genotype, has emerged in the last few decades, raising concerns particularly due to its atypical clinical presentation (ulcerative proctitis) and circulation among men who have sex with men (MSM). Here, we report an LGV outbreak, mostly affecting human immunodeficiency virus-positive MSM engaging in high-risk sexual practices, caused by an L2b strain with a rather unique non-LGV ompA signature that precluded the laboratory notification of this outbreak as LGV. C. trachomatis whole-genome capture and sequencing directly from clinical samples was applied to deeply characterize the genomic backbone of this novel LGV outbreak-causing clone. It revealed a chimeric genome structure due to the genetic transfer of ompA and four neighbouring genes from a serovar D/Da strain, likely possessing the genomic backbone associated with the more prevalent urogenital genotypes (T1 clade), to an LGV (L2b) strain. The hybrid L2b/D-Da strain presents the adhesin and immunodominant antigen MOMP (major outer membrane protein) (encoded by ompA) with an epitope repertoire typical of non-invasive genital strains, while keeping the genome-dispersed virulence fingerprint of a classical LGV strain. As previously reported for inter-clade ompA exchange among non-LGV clades, this novel C. trachomatis genomic mosaic involving a contemporary epidemiologically and clinically relevant LGV strain may have implications on its transmission, tissue tropism and pathogenic capabilities. The emergence of variants with epidemic and pathogenic potential highlights the need for more focused surveillance strategies to capture C. trachomatis evolution in action.

Virulence and antibiotic resistance plasticity of Arcobacter butzleri: Insights on the genomic diversity of an emerging human pathogen

Arcobacter butzleri is a foodborne emerging human pathogen, frequently displaying a multidrug resistant character. Still, the lack of comprehensive genome-scale comparative analysis has limited our knowledge on A. butzleri diversification and pathogenicity. Here, we performed a deep genome analysis of A. butzleri focused on decoding its core- and pan-genome diversity and specific genetic traits underlying its pathogenic potential and diverse ecology. A. butzleri (genome size 2.07-2.58 Mbp) revealed a large open pan-genome with 7474 genes (about 50% being singletons) and a small but diverse core-genome with 1165 genes. It presents a plastic virulome (including newly identified determinants), marked by the differential presence of multiple adaptation-related virulence factors, such as the urease cluster ureD(AB)CEFG (phenotypically confirmed), the hypervariable hemagglutinin-encoding hecA, a type I secretion system (T1SS) harboring another agglutinin and a novel VirB/D4 T4SS likely linked to interbacterial competition and cytotoxicity. In addition, A. butzleri harbors a large repertoire of efflux pumps (EPs) and other antibiotic resistant determinants. We unprecedentedly describe a genetic mechanism of A. butzleri macrolides resistance, (inactivation of a TetR repressor likely regulating an EP). Fluoroquinolones resistance correlated with Thr-85-Ile in GyrA and ampicillin resistance was linked to an OXA-15-like β-lactamase. Remarkably, by decoding the polymorphism pattern of the main antigen PorA, we show that A. butzleri is able to exchange porA as a whole and/or hypervariable epitope-encoding regions separately, leading to a multitude of chimeric PorA presentations that can impact pathogen-host interaction during infection. Ultimately, our unprecedented screening of short sequence repeats indicates that phase variation likely modulates A. butzleri key adaptive functions. In summary, this study constitutes a turning point on A. butzleri comparative genomics revealing that this human gastrointestinal pathogen is equipped with vast and diverse virulence and antibiotic resistance arsenals that open a multitude of phenotypic fingerprints for environmental/host adaptation and pathogenicity.

Towards a rapid sequencing-based molecular surveillance and mosaicism investigation of Toxoplasma gondii

Advances in molecular epidemiology of Toxoplasma gondii are hampered by technical and cost-associated hurdles underlying the acquisition of genomic data from parasites. In order to implement an enhanced genotyping approach for molecular surveillance of T. gondii, we applied a multi-locus amplicon-based sequencing strategy to samples associated with human infection. This approach, targeting genome-dispersed polymorphic loci potentially involved in adaptation and virulence, genetically discriminated almost all 68 studied strains and revealed a scenario of marked genomic mosaicism. Two-thirds (n = 43) of all strains were classified as recombinant, although recombination seemed to be linked to the classical archetypal lineage. While 92% of the Sag2 archetype I strains revealed genetic mosaicism, only 45% of Sag2 archetype II strains were identified as recombinant. Contrarily to the virulence-associated archetype I, most type II strains (regardless of their recombination background) were non-virulent in mouse. Besides Sag2, some of the newly studied loci (namely the type I/I-like alleles of Sag1, B17, PK1, and Sag3 and type III/III-like alleles of TgM-A) constitute promising candidates to rapidly infer T. gondii mouse virulence. Our successful attempt to capture microsatellite length variation launches good perspectives for the straightforward transition from the laborious intensive historical method to more informative next-generation sequencing (NGS)/bioinformatics-based methodologies. Overall, while T. gondii whole-genome sequencing will be hardly feasible in most laboratories, this study shows that a discrete loci panel has the potential to improve the molecular epidemiology of T. gondii towards a better monitoring of circulating genotypes with clinical importance.

Proficiency Testing of Virus Diagnostics Based on Bioinformatics Analysis of Simulated In Silico High-Throughput Sequencing Data Sets

Quality management and independent assessment of high-throughput sequencing-based virus diagnostics have not yet been established as a mandatory approach for ensuring comparable results. The sensitivity and specificity of viral high-throughput sequence data analysis are highly affected by bioinformatics processing using publicly available and custom tools and databases and thus differ widely between individuals and institutions. Here we present the results of the COMPARE [Collaborative Management Platform for Detection and Analyses of (Re-)emerging and Foodborne Outbreaks in Europe] in silico virus proficiency test. An artificial, simulated in silico data set of Illumina HiSeq sequences was provided to 13 different European institutes for bioinformatics analysis to identify viral pathogens in high-throughput sequence data. Comparison of the participants' analyses shows that the use of different tools, programs, and databases for bioinformatics analyses can impact the correct identification of viral sequences from a simple data set. The identification of slightly mutated and highly divergent virus genomes has been shown to be most challenging. Furthermore, the interpretation of the results, together with a fictitious case report, by the participants showed that in addition to the bioinformatics analysis, the virological evaluation of the results can be important in clinical settings. External quality assessment and proficiency testing should become an important part of validating high-throughput sequencing-based virus diagnostics and could improve the harmonization, comparability, and reproducibility of results. There is a need for the establishment of international proficiency testing, like that established for conventional laboratory tests such as PCR, for bioinformatics pipelines and the interpretation of such results.