Prevalence of cefixime-resistant Neisseria gonorrhoeae in Melbourne, Australia, 2021-2022.

While ceftriaxone remains the first-line treatment for gonorrhoea, the US CDC recommended cefixime as a second-line treatment in 2021. We tested 1176 Neisseria gonorrhoeae isolates among clients attending the Melbourne Sexual Health Centre in 2021-2022. The prevalence of cefixime resistance was 6.3% (74/1176), azithromycin resistance was 4.9% (58/1176) and ceftriaxone resistance was 0% (0/1176). Cefixime resistance was the highest among women (16.4%, 10/61), followed by men-who-have-sex-with-women (6.4%, 7/109), and men-who-have-sex-with-men (5.8%, 57/982). The prevalence of cefixime-resistant N. gonorrhoeae exceeds the threshold of the 5% resistance level recommended by the World Health Organization; and thus, cefixime treatment would have limited benefits in Australia.

Emergence and clonal expansion of a qacA-harbouring sequence type 45 lineage of methicillin-resistant Staphylococcus aureus.

The past decade has seen an increase in the prevalence of sequence type (ST) 45 methicillin-resistant Staphylococcus aureus (MRSA), yet the underlying drivers for its emergence and spread remain unclear. To better understand the worldwide dissemination of ST45 S. aureus, we performed phylogenetic analyses of Australian isolates, supplemented with a global population of ST45 S. aureus genomes. Our analyses revealed a distinct lineage of multidrug-resistant ST45 MRSA harbouring qacA, predominantly found in Australia and Singapore. Bayesian inference predicted that the acquisition of qacA occurred in the late 1990s. qacA was integrated into a structurally variable region of the chromosome containing Tn552 (carrying blaZ) and Tn4001 (carrying aac(6')-aph(2")) transposable elements. Using mutagenesis and in vitro assays, we provide phenotypic evidence that qacA confers tolerance to chlorhexidine. These findings collectively suggest both antimicrobial resistance and the carriage of qacA may play a role in the successful establishment of ST45 MRSA.

Phage resistance in Klebsiella pneumoniae and bidirectional effects impacting antibiotic susceptibility.

OBJECTIVES: Bacteriophage (phage) therapy is a promising anti-infective option to combat antimicrobial resistance. However, the clinical utilization of phage therapy has been severely compromised by the potential emergence of phage resistance. Although certain phage resistance mechanisms can restore bacterial susceptibility to certain antibiotics, a lack of knowledge of phage resistance mechanisms hinders optimal use of phages and their combination with antibiotics. METHODS: Genome-wide transposon screening was performed with a mutant library of Klebsiella pneumoniae MKP103 to identify phage pKMKP103_1-resistant mutants. Phage-resistant phenotypes were evaluated by time-kill kinetics and efficiency of plating assays. Phage resistance mechanisms were investigated with adsorption, one-step growth, and mutation frequency assays. Antibiotic susceptibility was determined with broth microdilution and population analysis profiles. RESULTS: We observed a repertoire of phage resistance mechanisms in K pneumoniae, such as disruption of phage binding (fhuA::Tn and tonB::Tn), extension of the phage latent period (mnmE::Tn and rpoN::Tn), and increased mutation frequency (mutS::Tn and mutL::Tn). Notably, in contrast to the prevailing view that phage resistance re-sensitizes antibiotic-resistant bacteria, we observed a bidirectional steering effect on bacterial antibiotic susceptibility. Specifically, rpoN::Tn increased susceptibility to colistin while mutS::Tn and mutL::Tn increased resistance to rifampicin and colistin. DISCUSSION: Our findings demonstrate that K pneumoniae employs multiple strategies to overcome phage infection, which may result in enhanced or reduced antibiotic susceptibility. Mechanism-guided phage steering should be incorporated into phage therapy to better inform clinical decisions on phage-antibiotic combinations.

Using whole genome sequencing to characterize Clostridioides difficile isolates at a tertiary center in Melbourne, Australia.

Objective: Clostridioides difficile infection (CDI) is the commonest cause of healthcare-associated diarrhea and undergoes standardized surveillance and mandatory reporting in most Australian states and territories. Historically attributed to nosocomial spread, local and international whole genome sequencing (WGS) data suggest varied sources of acquisition. This study describes C. difficile genotypes isolated at a tertiary center in Melbourne, Australia, their likely source of acquisition, and common risk factors. Design: Retrospective observational study. Setting: The Royal Melbourne Hospital (RMH), a 570-bed tertiary center in Victoria, Australia. Methods: Short-read whole genome sequencing was performed on 75 out of 137 C. difficile isolates obtained from 1/5/2021 to 28/2/2022 and compared to previous data from 8/11/2015 to 1/11/2016. Existing data from infection control surveillance and electronic medical records were used for epidemiological and risk factor analysis. Results: Eighty-five (62.1%) of the 137 cases were defined as healthcare-associated from epidemiological data. On genome sequencing, 33 different multi-locus sequence type (MLST) subtypes were identified, with changes in population structure compared to the 2015-16 period. Risk factors for CDI were present in 130 (94.9%) cases, including 108 (78.8%) on antibiotics, 86 (62.8%) on acid suppression therapy, and 25 (18.2) on chemotherapy. Conclusion: In both study periods, most C. difficile isolates were not closely related, suggesting varied sources of acquisition and that spread of C. difficile within the hospital was unlikely. Current infection control precautions may therefore warrant review. Underlying risk factors for CDI were common and may contribute to the proportion of healthcare-associated infections in the absence of proven hospital transmission.

Intratumoral presence of the genotoxic gut bacteria pks+ E. coli, Enterotoxigenic Bacteroides fragilis, and Fusobacterium nucleatum and their association with clinicopathological and molecular features of colorectal cancer.

BACKGROUND: This study aimed to investigate clinicopathological and molecular tumour features associated with intratumoral pks+ Escherichia coli (pks+E.coli+), pks+E.coli- (non-E.coli bacteria harbouring the pks island), Enterotoxigenic Bacteroides fragilis (ETBF) and Fusobacterium nucleatum (F. nucleatum). METHODS: We screened 1697 tumour-derived DNA samples from the Australasian Colorectal Cancer Family Registry, Melbourne Collaborative Cohort Study and the ANGELS study using targeted PCR. RESULTS: Pks+E.coli+ was associated with male sex (P < 0.01) and APC:c.835-8 A > G somatic mutation (P = 0.03). The association between pks+E.coli+ and APC:c.835-8 A > G was specific to early-onset CRCs (diagnosed<45years, P = 0.02). The APC:c.835-A > G was not associated with pks+E.coli- (P = 0.36). F. nucleatum was associated with DNA mismatch repair deficiency (MMRd), BRAF:c.1799T>A p.V600E mutation, CpG island methylator phenotype, proximal tumour location, and high levels of tumour infiltrating lymphocytes (Ps < 0.01). In the stratified analysis by MMRd subgroups, F. nucleatum was associated with Lynch syndrome, MLH1 methylated and double MMR somatic mutated MMRd subgroups (Ps < 0.01). CONCLUSION: Intratumoral pks+E.coli+ but not pks+E.coli- are associated with CRCs harbouring the APC:c.835-8 A > G somatic mutation, suggesting that this mutation is specifically related to DNA damage from colibactin-producing E.coli exposures. F. nucleatum was associated with both hereditary and sporadic MMRd subtypes, suggesting the MMRd tumour microenvironment is important for F. nucleatum colonisation irrespective of its cause.

Whole genome sequencing of drug-resistant Mycobacterium tuberculosis isolates in Victoria, Australia.

OBJECTIVES: Whole genome sequencing (WGS) can identify clusters, transmission patterns, and drug resistance mutations. This is important in low-burden settings such as Australia, as it can assist in efficient contact tracing and surveillance. METHODS: We conducted a retrospective cohort study using WGS from 155 genomically defined drug-resistant Mycobacterium tuberculosis (DR-TB) isolates collected between 2018-2021 in Victoria, Australia. Bioinformatic analysis was performed to identify resistance-conferring mutations, lineages, clusters and understand how local sequences compared with international context. RESULTS: Of the 155 sequences, 42% were identified as lineage 2 and 35% as lineage 1; 65.8% (102/155) were isoniazid mono-resistant, 8.4% were multi-drug resistant TB and 5.8% were pre-extensively drug-resistant / extensively drug-resistant TB. The most common mutations were observed in katG and fabG1 genes, especially at Ser315Thr and fabG1 -15 C>T for first-line drugs. Ser450Leu was the most frequent mutation in rpoB gene. Phylogenetic analysis confirmed that Victorian DR-TB were associated with importation events. There was little evidence of local transmission with only five isolate pairs. CONCLUSION: Isoniazid-resistant TB is the commonest DR-TB in Victoria, and the mutation profile is similar to global circulating DR-TB. Most cases are diagnosed among migrants with limited transmission. This study highlights the value of WGS in identification of clusters and resistance-conferring mutations. This information is crucial in supporting disease mitigation and treatment strategies.

Comparison of contemporary invasive and non-invasive Streptococcus pneumoniae isolates reveals new insights into circulating anti-microbial resistance determinants.

Streptococcus pneumoniae is a major human pathogen with a high burden of disease. Non-invasive isolates (those found in non-sterile sites) are thought to be a key source of invasive isolates (those found in sterile sites) and a reservoir of anti-microbial resistance (AMR) determinants. Despite this, pneumococcal surveillance has almost exclusively focused on invasive isolates. We aimed to compare contemporaneous invasive and non-invasive isolate populations to understand how they interact and identify differences in AMR gene distribution. We used a combination of whole-genome sequencing and phenotypic anti-microbial susceptibility testing and a data set of invasive (n = 1,288) and non-invasive (n = 186) pneumococcal isolates, collected in Victoria, Australia, between 2018 and 2022. The non-invasive population had increased levels of antibiotic resistance to multiple classes of antibiotics including beta-lactam antibiotics penicillin and ceftriaxone. We identified genomic intersections between the invasive and non-invasive populations and no distinct phylogenetic clustering of the two populations. However, this analysis revealed sub-populations overrepresented in each population. The sub-populations that had high levels of AMR were overrepresented in the non-invasive population. We determined that WamR-Pneumo was the most accurate in silico tool for predicting resistance to the antibiotics tested. This tool was then used to assess the allelic diversity of the penicillin-binding protein genes, which acquire mutations leading to beta-lactam antibiotic resistance, and found that they were highly conserved (≥80% shared) between the two populations. These findings show the potential of non-invasive isolates to serve as reservoirs of AMR determinants.

The importance of utilizing travel history metadata for informative phylogeographical inferences: a case study of early SARS-CoV-2 introductions into Australia.

Inferring the spatiotemporal spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) via Bayesian phylogeography has been complicated by the overwhelming sampling bias present in the global genomic dataset. Previous work has demonstrated the utility of metadata in addressing this bias. Specifically, the inclusion of recent travel history of SARS-CoV-2-positive individuals into extended phylogeographical models has demonstrated increased accuracy of estimates, along with proposing alternative hypotheses that were not apparent using only genomic and geographical data. However, as the availability of comprehensive epidemiological metadata is limited, many of the current estimates rely on sequence data and basic metadata (i.e. sample date and location). As the bias within the SARS-CoV-2 sequence dataset is extensive, the degree to which we can rely on results drawn from standard phylogeographical models (i.e. discrete trait analysis) that lack integrated metadata is of great concern. This is particularly important when estimates influence and inform public health policy. We compared results generated from the same dataset, using two discrete phylogeographical models: one including travel history metadata and one without. We utilized sequences from Victoria, Australia, in this case study for two unique properties. Firstly, the high proportion of cases sequenced throughout 2020 within Victoria and the rest of Australia. Secondly, individual travel history was collected from returning travellers in Victoria during the first wave (January to May) of the coronavirus disease 2019 (COVID-19) pandemic. We found that the implementation of individual travel history was essential for the estimation of SARS-CoV-2 movement via discrete phylogeography models. Without the additional information provided by the travel history metadata, the discrete trait analysis could not be fit to the data due to numerical instability. We also suggest that during the first wave of the COVID-19 pandemic in Australia, the primary driving force behind the spread of SARS-CoV-2 was viral importation from international locations. This case study demonstrates the necessity of robust genomic datasets supplemented with epidemiological metadata for generating accurate estimates from phylogeographical models in datasets that have significant sampling bias. For future work, we recommend the collection of metadata in conjunction with genomic data. Furthermore, we highlight the risk of applying phylogeographical models to biased datasets without incorporating appropriate metadata, especially when estimates influence public health policy decision making.

Public health implementation of pathogen genomics: the role for accreditation and application of ISO standards.

Pathogen genomics has transitioned rapidly from the research setting into a powerful tool now routinely used in public health microbiology, for surveillance, outbreak investigations and disease control. As these investigations can have significant public health, treatment and legal impacts, we must ensure the accuracy of these results through validation of testing processes. For laboratories working in this space, it is important to approach this work with a quality and accreditation framework in mind, working towards implementation of quality systems and test validation that meet international regulatory standards. Here we outline the key international standards and processes that lead toward accreditation for pathogen genomics.

Genomic epidemiology offers high resolution estimates of serial intervals for COVID-19.

Serial intervals - the time between symptom onset in infector and infectee - are a fundamental quantity in infectious disease control. However, their estimation requires knowledge of individuals' exposures, typically obtained through resource-intensive contact tracing efforts. We introduce an alternate framework using virus sequences to inform who infected whom and thereby estimate serial intervals. We apply our technique to SARS-CoV-2 sequences from case clusters in the first two COVID-19 waves in Victoria, Australia. We find that our approach offers high resolution, cluster-specific serial interval estimates that are comparable with those obtained from contact data, despite requiring no knowledge of who infected whom and relying on incompletely-sampled data. Compared to a published serial interval, cluster-specific serial intervals can vary estimates of the effective reproduction number by a factor of 2-3. We find that serial interval estimates in settings such as schools and meat processing/packing plants are shorter than those in healthcare facilities.

Population structure, serotype distribution and antibiotic resistance of Streptococcus pneumoniae causing invasive disease in Victoria, Australia.

Streptococcus pneumoniae is a major human pathogen and can cause a range of conditions from asymptomatic colonization to invasive pneumococcal disease (IPD). The epidemiology and distribution of IPD-causing serotypes in Australia has undergone large changes following the introduction of the 7-valent pneumococcal conjugate vaccine (PCV) in 2005 and the 13-valent PCV in 2011. In this study, to provide a contemporary understanding of the IPD causing population in Victoria, Australia, we aimed to examine the population structure and prevalence of antimicrobial resistance using whole-genome sequencing and comprehensive antimicrobial susceptibility data of 1288 isolates collected between 2018 and 2022. We observed high diversity among the isolates with 52 serotypes, 203 sequence types (STs) and 70 Global Pneumococcal Sequencing Project Clusters (GPSCs) identified. Serotypes contained in the 13v-PCV represented 35.3 % (n=405) of isolates. Antimicrobial resistance (AMR) to at least one antibiotic was identified in 23.8 % (n=358) of isolates with penicillin resistance the most prevalent (20.3 %, n=261 using meningitis breakpoints and 5.1 % n=65 using oral breakpoints). Of the AMR isolates, 28 % (n=101) were multidrug resistant (MDR) (resistant to three or more drug classes). Vaccination status of cases was determined for a subset of isolates with 34 cases classified as vaccine failure events (fully vaccinated IPD cases of vaccine serotype). However, no phylogenetic association with failure events was observed. Within the highly diverse IPD population, we identified six high-risk sub-populations of public health concern characterized by high prevalence, high rates of AMR and MDR, or serotype inclusion in vaccines. High-risk serotypes included serotypes 3, 19F, 19A, 14, 11A, 15A and serofamily 23. In addition, we present our data validating seroBA for in silico serotyping to facilitate ISO-accreditation of this test in routine use in a public health reference laboratory and have made this data set available. This study provides insights into the population dynamics, highlights non-vaccine serotypes of concern that are highly resistant, and provides a genomic framework for the ongoing surveillance of IPD in Australia which can inform next-generation IPD prevention strategies.

An ISO-certified genomics workflow for identification and surveillance of antimicrobial resistance.

Realising the promise of genomics to revolutionise identification and surveillance of antimicrobial resistance (AMR) has been a long-standing challenge in clinical and public health microbiology. Here, we report the creation and validation of abritAMR, an ISO-certified bioinformatics platform for genomics-based bacterial AMR gene detection. The abritAMR platform utilises NCBI's AMRFinderPlus, as well as additional features that classify AMR determinants into antibiotic classes and provide customised reports. We validate abritAMR by comparing with PCR or reference genomes, representing 1500 different bacteria and 415 resistance alleles. In these analyses, abritAMR displays 99.9% accuracy, 97.9% sensitivity and 100% specificity. We also compared genomic predictions of phenotype for 864 Salmonella spp. against agar dilution results, showing 98.9% accuracy. The implementation of abritAMR in our institution has resulted in streamlined bioinformatics and reporting pathways, and has been readily updated and re-verified. The abritAMR tool and validation datasets are publicly available to assist laboratories everywhere harness the power of AMR genomics in professional practice.

Hidden Resistances: How Routine Whole-Genome Sequencing Uncovered an Otherwise Undetected blaNDM-1 Gene in Vibrio alginolyticus from Imported Seafood.

Vibrio alginolyticus causes vibriosis of marine vertebrates, invertebrates, and humans, and while there have been several reports of multidrug resistance in V. alginolyticus, carbapenem resistance is rare. V. alginolyticus strain AUSMDU00064140 was isolated in Melbourne, Australia, from imported prawns. Routine genomic surveillance detected the presence of a full-length blaNDM-1 gene, subsequently shown to be collocated with additional acquired antimicrobial resistance genes on a resistance cassette on the largest chromosome, flanked by mobilization gene annotations. Comparisons to a previously described V. alginolyticus plasmid, pC1349, revealed differing gene content and arrangements between the resistance cassettes. Phylogenetic analysis was performed against a local and global data set (n = 109), demonstrating that AUSMDU00064140 was distinct and did not cluster with any other strains. Despite the presence of the complete blaNDM-1 gene and positive phenotypic assays for carbapenemase production, carbapenem MICs were low (meropenem MIC ≤0.5 mg/liter). However, it is still possible that this gene may be transferred to another species in the environment or a host, causing phenotypic carbapenem resistance and presenting a risk of great public health concern. IMPORTANCE Carbapenems are last-line antimicrobials, vital for use in human medicine. Antimicrobial resistance determinants such as blaNDM (New Delhi metallo-ß-lactamase producing) genes conferring resistance to the carbapenem class of antimicrobials, are typically found in Enterobacterales (first described in 2009 from a Klebsiella pneumoniae isolate). Our study shows that Vibrio alginolyticus isolated from cooked prawn is able to harbor antimicrobial resistance (AMR) genes of public health concern, specifically a chromosomally located blaNDM-1 gene, and there is the potential for transmission of resistance genes. This may be linked with antimicrobial use in low- and middle-income settings, which has typically been high, unregulated, or not reported. Many countries, including Thailand, have implemented national strategic plans to incorporate the World Health Organization (WHO)'s Global Action Plan (2015) recommendations of a global One Health approach, including increased resources for surveillance of antimicrobial usage and AMR; however, efficient antimicrobial surveillance systems incorporating genomic and phenotypic testing of isolates are still lacking in many jurisdictions.

Whole genome sequencing for tuberculosis in Victoria, Australia: A genomic implementation study from 2017 to 2020.

Background: Whole genome sequencing (WGS) is increasingly used by tuberculosis (TB) programs to monitor Mycobacterium tuberculosis (Mtb) transmission. We aimed to characterise the molecular epidemiology of TB and Mtb transmission in the low-incidence setting of Victoria, Australia, and assess the utility of WGS. Methods: WGS was performed on all first Mtb isolates from TB cases from 2017 to 2020. Potential clusters (≤12 single nucleotide polymorphisms [SNPs]) were investigated for epidemiological links. Transmission events in highly-related (≤5 SNPs) clusters were classified as likely or possible, based on the presence or absence of an epidemiological link, respectively. Case characteristics and transmission settings (as defined by case relationship) were summarised. Poisson regression was used to examine associations with secondary case number. Findings: Of 1844 TB cases, 1276 (69.2%) had sequenced isolates, with 182 (14.2%) in 54 highly-related clusters, 2-40 cases in size. Following investigation, 140 cases (11.0% of sequenced) were classified as resulting from likely/possible local-transmission, including 82 (6.4%) for which transmission was likely. Common identified transmission settings were social/religious (26.4%), household (22.9%) and family living in different households (7.1%), but many were uncertain (41.4%). While household transmission featured in many clusters (n = 24), clusters were generally smaller (median = 3 cases) than the fewer that included transmission in social/religious settings (n = 12, median = 7.5 cases). Sputum-smear-positivity was associated with higher secondary case numbers. Interpretation: WGS results suggest Mtb transmission commonly occurs outside the household in our low-incidence setting. Further work is required to optimise the use of WGS in public health management of TB. Funding: The Victorian Tuberculosis Program receives block funding for activities including case management and contact tracing from the Victorian Department of Health. No specific funding for this report was received by manuscript authors or the Victorian Tuberculosis Program, and the funders had no role in the study design, data collection, data analysis, interpretation or report writing.

Genomic characterisation reveals a dominant lineage of SARS-CoV-2 in Papua New Guinea.

The coronavirus disease pandemic has highlighted the utility of pathogen genomics as a key part of comprehensive public health response to emerging infectious diseases threats, however, the ability to generate, analyse, and respond to pathogen genomic data varies around the world. Papua New Guinea (PNG), which has limited in-country capacity for genomics, has experienced significant outbreaks of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) with initial genomics data indicating a large proportion of cases were from lineages that are not well defined within the current nomenclature. Through a partnership between in-country public health agencies and academic organisations, industry, and a public health genomics reference laboratory in Australia a system for routine SARS-CoV-2 genomics from PNG was established. Here we aim to characterise and describe the genomics of PNG's second wave and examine the sudden expansion of a lineage that is not well defined but very prevalent in the Western Pacific region. We generated 1797 sequences from cases in PNG and performed phylogenetic and phylodynamic analyses to examine the outbreak and characterise the circulating lineages and clusters present. Our results reveal the rapid expansion of the B.1.466.2 and related lineages within PNG, from multiple introductions into the country. We also highlight the difficulties that unstable lineage assignment causes when using genomics to assist with rapid cluster definitions.

State-wide genomic epidemiology investigations of COVID-19 in healthcare workers in 2020 Victoria, Australia: Qualitative thematic analysis to provide insights for future pandemic preparedness.

Background: COVID-19 has affected many healthcare workers (HCWs) globally. We performed state-wide SARS-CoV-2 genomic epidemiological investigations to identify HCW transmission dynamics and provide recommendations to optimise healthcare system preparedness for future outbreaks. Methods: Genome sequencing was attempted on all COVID-19 cases in Victoria, Australia. We combined genomic and epidemiologic data to investigate the source of HCW infections across multiple healthcare facilities (HCFs) in the state. Phylogenetic analysis and fine-scale hierarchical clustering were performed for the entire dataset including community and healthcare cases. Facilities provided standardised epidemiological data and putative transmission links. Findings: Between March-October 2020, approximately 1,240 HCW COVID-19 infection cases were identified; 765 are included here, requested for hospital investigations. Genomic sequencing was successful for 612 (80%) cases. Thirty-six investigations were undertaken across 12 HCFs. Genomic analysis revealed that multiple introductions of COVID-19 into facilities (31/36) were more common than single introductions (5/36). Major contributors to HCW acquisitions included mobility of staff and patients between wards and facilities, and characteristics and behaviours of patients that generated numerous secondary infections. Key limitations at the HCF level were identified. Interpretation: Genomic epidemiological analyses enhanced understanding of HCW infections, revealing unsuspected clusters and transmission networks. Combined analysis of all HCWs and patients in a HCF should be conducted, supported by high rates of sequencing coverage for all cases in the population. Established systems for integrated genomic epidemiological investigations in healthcare settings will improve HCW safety in future pandemics. Funding: The Victorian Government, the National Health and Medical Research Council Australia, and the Medical Research Future Fund.

SARS-CoV-2 infection results in immune responses in the respiratory tract and peripheral blood that suggest mechanisms of disease severity.

Respiratory tract infection with SARS-CoV-2 results in varying immunopathology underlying COVID-19. We examine cellular, humoral and cytokine responses covering 382 immune components in longitudinal blood and respiratory samples from hospitalized COVID-19 patients. SARS-CoV-2-specific IgM, IgG, IgA are detected in respiratory tract and blood, however, receptor-binding domain (RBD)-specific IgM and IgG seroconversion is enhanced in respiratory specimens. SARS-CoV-2 neutralization activity in respiratory samples correlates with RBD-specific IgM and IgG levels. Cytokines/chemokines vary between respiratory samples and plasma, indicating that inflammation should be assessed in respiratory specimens to understand immunopathology. IFN-α2 and IL-12p70 in endotracheal aspirate and neutralization in sputum negatively correlate with duration of hospital stay. Diverse immune subsets are detected in respiratory samples, dominated by neutrophils. Importantly, dexamethasone treatment does not affect humoral responses in blood of COVID-19 patients. Our study unveils differential immune responses between respiratory samples and blood, and shows how drug therapy affects immune responses during COVID-19.