Integrated longitudinal multiomics study identifies immune programs associated with acute COVID-19 severity and mortality.

BACKGROUNDPatients hospitalized for COVID-19 exhibit diverse clinical outcomes, with outcomes for some individuals diverging over time even though their initial disease severity appears similar to that of other patients. A systematic evaluation of molecular and cellular profiles over the full disease course can link immune programs and their coordination with progression heterogeneity.METHODSWe performed deep immunophenotyping and conducted longitudinal multiomics modeling, integrating 10 assays for 1,152 Immunophenotyping Assessment in a COVID-19 Cohort (IMPACC) study participants and identifying several immune cascades that were significant drivers of differential clinical outcomes.RESULTSIncreasing disease severity was driven by a temporal pattern that began with the early upregulation of immunosuppressive metabolites and then elevated levels of inflammatory cytokines, signatures of coagulation, formation of neutrophil extracellular traps, and T cell functional dysregulation. A second immune cascade, predictive of 28-day mortality among critically ill patients, was characterized by reduced total plasma Igs and B cells and dysregulated IFN responsiveness. We demonstrated that the balance disruption between IFN-stimulated genes and IFN inhibitors is a crucial biomarker of COVID-19 mortality, potentially contributing to failure of viral clearance in patients with fatal illness.CONCLUSIONOur longitudinal multiomics profiling study revealed temporal coordination across diverse omics that potentially explain the disease progression, providing insights that can inform the targeted development of therapies for patients hospitalized with COVID-19, especially those who are critically ill.TRIAL REGISTRATIONClinicalTrials.gov NCT04378777.FUNDINGNIH (5R01AI135803-03, 5U19AI118608-04, 5U19AI128910-04, 4U19AI090023-11, 4U19AI118610-06, R01AI145835-01A1S1, 5U19AI062629-17, 5U19AI057229-17, 5U19AI125357-05, 5U19AI128913-03, 3U19AI077439-13, 5U54AI142766-03, 5R01AI104870-07, 3U19AI089992-09, 3U19AI128913-03, and 5T32DA018926-18); NIAID, NIH (3U19AI1289130, U19AI128913-04S1, and R01AI122220); and National Science Foundation (DMS2310836).

No Remdesivir Resistance Observed in the Phase 3 Severe and Moderate COVID-19 SIMPLE Trials.

Remdesivir (RDV) is a broad-spectrum nucleotide analog prodrug approved for the treatment of COVID-19 in hospitalized and non-hospitalized patients with clinical benefit demonstrated in multiple Phase 3 trials. Here we present SARS-CoV-2 resistance analyses from the Phase 3 SIMPLE clinical studies evaluating RDV in hospitalized participants with severe or moderate COVID-19 disease. The severe and moderate studies enrolled participants with radiologic evidence of pneumonia and a room-air oxygen saturation of ≤94% or >94%, respectively. Virology sample collection was optional in the study protocols. Sequencing and related viral load data were obtained retrospectively from participants at a subset of study sites with local sequencing capabilities (10 of 183 sites) at timepoints with detectable viral load. Among participants with both baseline and post-baseline sequencing data treated with RDV, emergent Nsp12 substitutions were observed in 4 of 19 (21%) participants in the severe study and none of the 2 participants in the moderate study. The following 5 substitutions emerged: T76I, A526V, A554V, E665K, and C697F. The substitutions T76I, A526V, A554V, and C697F had an EC50 fold change of ≤1.5 relative to the wildtype reference using a SARS-CoV-2 subgenomic replicon system, indicating no significant change in the susceptibility to RDV. The phenotyping of E665K could not be determined due to a lack of replication. These data reveal no evidence of relevant resistance emergence and further confirm the established efficacy profile of RDV with a high resistance barrier in COVID-19 patients.

Prophage-encoded methyltransferase drives adaptation of community-acquired methicillin-resistant Staphylococcus aureus.

We recently described the evolution of a community-acquired methicillin-resistant Staphylococcus aureus (CA-MRSA) USA300 variant responsible for an outbreak of skin and soft tissue infections. Acquisition of a mosaic version of the Φ11 prophage (mΦ11) that increases skin abscess size was an early step in CA-MRSA adaptation that primed the successful spread of the clone. The present report shows how prophage mΦ11 exerts its effect on virulence for skin infection without encoding a known toxin or fitness genes. Abscess size and skin inflammation were associated with DNA methylase activity of an mΦ11-encoded adenine methyltransferase (designated pamA). pamA increased expression of fibronectin-binding protein A (fnbA; FnBPA), and inactivation of fnbA eliminated the effect of pamA on abscess virulence without affecting strains lacking pamA. Thus, fnbA is a pamA-specific virulence factor. Mechanistically, pamA was shown to promote biofilm formation in vivo in skin abscesses, a phenotype linked to FnBPA's role in biofilm formation. Collectively, these data reveal a novel mechanism-epigenetic regulation of staphylococcal gene expression-by which phage can regulate virulence to drive adaptive leaps by S. aureus.

A compendium of multi-omics data illuminating host responses to lethal human virus infections.

Human infections caused by viral pathogens trigger a complex gamut of host responses that limit disease, resolve infection, generate immunity, and contribute to severe disease or death. Here, we present experimental methods and multi-omics data capture approaches representing the global host response to infection generated from 45 individual experiments involving human viruses from the Orthomyxoviridae, Filoviridae, Flaviviridae, and Coronaviridae families. Analogous experimental designs were implemented across human or mouse host model systems, longitudinal samples were collected over defined time courses, and global multi-omics data (transcriptomics, proteomics, metabolomics, and lipidomics) were acquired by microarray, RNA sequencing, or mass spectrometry analyses. For comparison, we have included transcriptomics datasets from cells treated with type I and type II human interferon. Raw multi-omics data and metadata were deposited in public repositories, and we provide a central location linking the raw data with experimental metadata and ready-to-use, quality-controlled, statistically processed multi-omics datasets not previously available in any public repository. This compendium of infection-induced host response data for reuse will be useful for those endeavouring to understand viral disease pathophysiology and network biology.

Transcriptomic landscape of human induced pluripotent stem cell-derived osteogenic differentiation identifies a regulatory role of KLF16.

Osteogenic differentiation is essential for bone development and metabolism, but the underlying gene regulatory networks have not been well investigated. We differentiated mesenchymal stem cells, derived from 20 human induced pluripotent stem cell lines, into preosteoblasts and osteoblasts, and performed systematic RNA-seq analyses of 60 samples for differential gene expression. We noted a highly significant correlation in expression patterns and genomic proximity among transcription factor (TF) and long noncoding RNA (lncRNA) genes. We identified TF-TF regulatory networks, regulatory roles of lncRNAs on their neighboring coding genes for TFs and splicing factors, and differential splicing of TF, lncRNA, and splicing factor genes. TF-TF regulatory and gene co-expression network analyses suggested an inhibitory role of TF KLF16 in osteogenic differentiation. We demonstrate that in vitro overexpression of human KLF16 inhibits osteogenic differentiation and mineralization, and in vivo Klf16+/- mice exhibit increased bone mineral density, trabecular number, and cortical bone area. Thus, our model system highlights the regulatory complexity of osteogenic differentiation and identifies novel osteogenic genes.

SARS-CoV-2-infection- and vaccine-induced antibody responses are long lasting with an initial waning phase followed by a stabilization phase.

It is thought that mRNA-based vaccine-induced immunity to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) wanes quickly, based mostly on short-term studies. Here, we analyzed the kinetics and durability of the humoral responses to SARS-CoV-2 infection and vaccination using >8,000 longitudinal samples collected over a 3-year period in New York City. Upon primary immunization, participants with pre-existing immunity mounted higher antibody responses faster and achieved higher steady-state antibody titers than naive individuals. Antibody kinetics were characterized by two phases: an initial rapid decay, followed by a stabilization phase with very slow decay. Booster vaccination equalized the differences in antibody concentration between participants with and without hybrid immunity, but the peak antibody titers decreased with each successive antigen exposure. Breakthrough infections increased antibodies to similar titers as an additional vaccine dose in naive individuals. Our study provides strong evidence that SARS-CoV-2 antibody responses are long lasting, with initial waning followed by stabilization.

Features of acute COVID-19 associated with post-acute sequelae of SARS-CoV-2 phenotypes: results from the IMPACC study.

Post-acute sequelae of SARS-CoV-2 (PASC) is a significant public health concern. We describe Patient Reported Outcomes (PROs) on 590 participants prospectively assessed from hospital admission for COVID-19 through one year after discharge. Modeling identified 4 PRO clusters based on reported deficits (minimal, physical, mental/cognitive, and multidomain), supporting heterogenous clinical presentations in PASC, with sub-phenotypes associated with female sex and distinctive comorbidities. During the acute phase of disease, a higher respiratory SARS-CoV-2 viral burden and lower Receptor Binding Domain and Spike antibody titers were associated with both the physical predominant and the multidomain deficit clusters. A lower frequency of circulating B lymphocytes by mass cytometry (CyTOF) was observed in the multidomain deficit cluster. Circulating fibroblast growth factor 21 (FGF21) was significantly elevated in the mental/cognitive predominant and the multidomain clusters. Future efforts to link PASC to acute anti-viral host responses may help to better target treatment and prevention of PASC.

IgM N-glycosylation correlates with COVID-19 severity and rate of complement deposition.

The glycosylation of IgG plays a critical role during human severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection, activating immune cells and inducing cytokine production. However, the role of IgM N-glycosylation has not been studied during human acute viral infection. The analysis of IgM N-glycosylation from healthy controls and hospitalized coronavirus disease 2019 (COVID-19) patients reveals increased high-mannose and sialylation that correlates with COVID-19 severity. These trends are confirmed within SARS-CoV-2-specific immunoglobulin N-glycan profiles. Moreover, the degree of total IgM mannosylation and sialylation correlate significantly with markers of disease severity. We link the changes of IgM N-glycosylation with the expression of Golgi glycosyltransferases. Lastly, we observe antigen-specific IgM antibody-dependent complement deposition is elevated in severe COVID-19 patients and modulated by exoglycosidase digestion. Taken together, this work links the IgM N-glycosylation with COVID-19 severity and highlights the need to understand IgM glycosylation and downstream immune function during human disease.

Clinical, microbiological and genomic characterization of Gram-negative bacteria with dual carbapenemases as identified by rapid molecular testing.

Objective: Dual carbapenemase-producing organisms (DCPOs) are an emerging threat that expands the spectrum of antimicrobial resistance. There is limited literature on the clinical and genetic epidemiology of DCPOs. Methods: DCPO isolates were identified by Xpert® Carba-R PCR testing of routine diagnostic cultures performed from 2018 to 2021 at a New York City health system. WGS was performed by Illumina and/or PacBio. Medical records of patients were reviewed for clinical and epidemiological data. Results: Twenty-six DCPO isolates were obtained from 13 patients. Klebsiella pneumoniae (n = 22) was most frequent, followed by Pseudomonas aeruginosa (n = 2), Escherichia coli (n = 1) and Enterobacter cloacae (n = 1). The most common DCPO combination was blaNDM/blaOXA-48-like (n = 16). Notably, 1.05% (24/2290) of carbapenem-resistant Enterobacterales isolates were identified as DCPOs. The susceptibility profiles matched the identified resistance genes, except for a K. pneumoniae (blaKPC/blaOXA-48-like) isolate that was phenotypically susceptible to meropenem. Eleven patients were hospitalized within the year prior to admission, and received antibiotic(s) 1 month prior. Seven patients were originally from outside the USA. Hypertension, kidney disease and diabetes were frequent comorbidities. Death in two cases was attributed to DCPO infection. WGS of eight isolates showed that carbapenemases were located on distinct plasmids, except for one K. pneumoniae isolate where NDM and KPC carbapenemases were located on a single IncC-type plasmid backbone. Conclusions: Here we characterized a series of DCPOs from New York City. Foreign travel, prior hospitalization, antibiotic usage and comorbidities were common among DCPO cases. All carbapenemases were encoded on plasmids, which may facilitate horizontal transfer.

The impact of S2 mutations on Omicron SARS-CoV-2 cell surface expression and fusogenicity.

SARS-CoV-2 Omicron subvariants are still emerging and spreading worldwide. These variants contain a high number of polymorphisms in the spike (S) glycoprotein that could potentially impact their pathogenicity and transmission. We have previously shown that the S:655Y and P681H mutations enhance S protein cleavage and syncytia formation. Interestingly, these polymorphisms are present in Omicron S protein. Here, we characterized the cleavage efficiency and fusogenicity of the S protein of different Omicron sublineages. Our results showed that Omicron BA.1 subvariant is efficiently cleaved but it is poorly fusogenic compared to previous SARS-CoV-2 strains. To understand the basis of this phenotype, we generated chimeric S protein using combinations of the S1 and S2 domains from WA1, Delta and Omicron BA.1 variants. We found that the S2 domain of Omicron BA.1 hindered efficient cell-cell fusion. Interestingly, this domain only contains six unique polymorphisms never detected before in ancestral SARS-CoV-2 variants. WA1614G S proteins containing the six individuals S2 Omicron mutations were assessed for their fusogenicity and S surface expression after transfection in cells. Results showed that the S:N856K and N969K substitutions decreased syncytia formation and impacted S protein cell surface levels. However, we observed that "first-generation" Omicron sublineages that emerged subsequently, had convergently evolved to an enhanced fusogenic activity and S expression on the surface of infected cells while "second-generation" Omicron variants have highly diverged and showed lineage-specific fusogenic properties. Importantly, our findings could have potential implications in the improvement and redesign of COVID-19 vaccines.

Cellular mechanisms associated with sub-optimal immune responses to SARS-CoV-2 bivalent booster vaccination in patients with Multiple Myeloma.

BACKGROUND: The real-world impact of bivalent vaccines for wild type (WA.1) and Omicron variant (BA.5) is largely unknown in immunocompromised patients with Multiple Myeloma (MM). We characterize the humoral and cellular immune responses in patients with MM before and after receiving the bivalent booster, including neutralizing assays to identify patterns associated with continuing vulnerability to current variants (XBB1.16, EG5) in the current post-pandemic era. METHODS: We studied the humoral and cellular immune responses before and after bivalent booster immunization in 48 MM patients. Spike binding IgG antibody levels were measured by SARS-CoV-2 spike binding ELISA and neutralization capacity was assessed by a SARS-CoV-2 multi-cycle microneutralization assays to assess inhibition of live virus. We measured spike specific T-cell function using the QuantiFERON SARS-CoV-2 (Qiagen) assay as well as flow-cytometry based T-cell. In a subset of 38 patients, high-dimensional flow cytometry was performed to identify immune cell subsets associated with lack of humoral antibodies. FINDINGS: We find that bivalent vaccination provides significant boost in protection to the omicron variant in our MM patients, in a treatment specific manner. MM patients remain vulnerable to newer variants with mutations in the spike portion. Anti-CD38 and anti-BCMA therapies affect the immune machinery needed to produce antibodies. INTERPRETATION: Our study highlights varying immune responses observed in MM patients after receiving bivalent COVID-19 vaccination. Specifically, a subgroup of MM patients undergoing anti-CD38 and anti-BCMA therapy experience impairment in immune cells such DCs, B cells, NK cells and TFH cells, leading to an inability to generate adequate humoral and cellular responses to vaccination. FUNDING: National Cancer Institute (National Institutes of Health), National Institute of Allergy and Infectious Diseases (National Institutes of Health), NCI Serological Sciences Network for COVID-19 (SeroNet) and The Icahn School of Medicine at Mount Sinai.

Integrated longitudinal multi-omics study identifies immune programs associated with COVID-19 severity and mortality in 1152 hospitalized participants.

Hospitalized COVID-19 patients exhibit diverse clinical outcomes, with some individuals diverging over time even though their initial disease severity appears similar. A systematic evaluation of molecular and cellular profiles over the full disease course can link immune programs and their coordination with progression heterogeneity. In this study, we carried out deep immunophenotyping and conducted longitudinal multi-omics modeling integrating ten distinct assays on a total of 1,152 IMPACC participants and identified several immune cascades that were significant drivers of differential clinical outcomes. Increasing disease severity was driven by a temporal pattern that began with the early upregulation of immunosuppressive metabolites and then elevated levels of inflammatory cytokines, signatures of coagulation, NETosis, and T-cell functional dysregulation. A second immune cascade, predictive of 28-day mortality among critically ill patients, was characterized by reduced total plasma immunoglobulins and B cells, as well as dysregulated IFN responsiveness. We demonstrated that the balance disruption between IFN-stimulated genes and IFN inhibitors is a crucial biomarker of COVID-19 mortality, potentially contributing to the failure of viral clearance in patients with fatal illness. Our longitudinal multi-omics profiling study revealed novel temporal coordination across diverse omics that potentially explain disease progression, providing insights that inform the targeted development of therapies for hospitalized COVID-19 patients, especially those critically ill.

Mapping SARS-CoV-2 antigenic relationships and serological responses.

During the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic, multiple variants escaping preexisting immunity emerged, causing reinfections of previously exposed individuals. Here, we used antigenic cartography to analyze patterns of cross-reactivity among 21 variants and 15 groups of human sera obtained after primary infection with 10 different variants or after messenger RNA (mRNA)-1273 or mRNA-1273.351 vaccination. We found antigenic differences among pre-Omicron variants caused by substitutions at spike-protein positions 417, 452, 484, and 501. Quantifying changes in response breadth over time and with additional vaccine doses, our results show the largest increase between 4 weeks and >3 months after a second dose. We found changes in immunodominance of different spike regions, depending on the variant an individual was first exposed to, with implications for variant risk assessment and vaccine-strain selection.

Impact of SARS-CoV-2 ORF6 and its variant polymorphisms on host responses and viral pathogenesis.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) encodes several proteins that inhibit host interferon responses. Among these, ORF6 antagonizes interferon signaling by disrupting nucleocytoplasmic trafficking through interactions with the nuclear pore complex components Nup98-Rae1. However, the roles and contributions of ORF6 during physiological infection remain unexplored. We assessed the role of ORF6 during infection using recombinant viruses carrying a deletion or loss-of-function (LoF) mutation in ORF6. ORF6 plays key roles in interferon antagonism and viral pathogenesis by interfering with nuclear import and specifically the translocation of IRF and STAT transcription factors. Additionally, ORF6 inhibits cellular mRNA export, resulting in the remodeling of the host cell proteome, and regulates viral protein expression. Interestingly, the ORF6:D61L mutation that emerged in the Omicron BA.2 and BA.4 variants exhibits reduced interactions with Nup98-Rae1 and consequently impairs immune evasion. Our findings highlight the role of ORF6 in antagonizing innate immunity and emphasize the importance of studying the immune evasion strategies of SARS-CoV-2.

IgM N-glycosylation correlates with COVID-19 severity and rate of complement deposition.

The glycosylation of IgG plays a critical role during human SARS-CoV-2, activating immune cells and inducing cytokine production. However, the role of IgM N-glycosylation has not been studied during acute viral infection in humans. In vitro evidence suggests that the glycosylation of IgM inhibits T cell proliferation and alters complement activation rates. The analysis of IgM N-glycosylation from healthy controls and hospitalized COVID-19 patients reveals that mannosylation and sialyation levels associate with COVID-19 severity. Specifically, we find increased di- and tri-sialylated glycans and altered mannose glycans in total serum IgM in severe COVID-19 patients when compared to moderate COVID-19 patients. This is in direct contrast with the decrease of sialic acid found on the serum IgG from the same cohorts. Moreover, the degree of mannosylation and sialylation correlated significantly with markers of disease severity: D-dimer, BUN, creatinine, potassium, and early anti-COVID-19 amounts of IgG, IgA, and IgM. Further, IL-16 and IL-18 cytokines showed similar trends with the amount of mannose and sialic acid present on IgM, implicating these cytokines' potential to impact glycosyltransferase expression during IgM production. When examining PBMC mRNA transcripts, we observe a decrease in the expression of Golgi mannosidases that correlates with the overall reduction in mannose processing we detect in the IgM N-glycosylation profile. Importantly, we found that IgM contains alpha-2,3 linked sialic acids in addition to the previously reported alpha-2,6 linkage. We also report that antigen-specific IgM antibody-dependent complement deposition is elevated in severe COVID-19 patients. Taken together, this work links the immunoglobulin M N-glycosylation with COVID-19 severity and highlights the need to understand the connection between IgM glycosylation and downstream immune function during human disease.

Continued Circulation of Highly Pathogenic H5 Influenza Viruses in Vietnamese Live Bird Markets in 2018-2021.

We isolated 77 highly pathogenic avian influenza viruses during routine surveillance in live poultry markets in northern provinces of Vietnam from 2018 to 2021. These viruses are of the H5N6 subtype and belong to HA clades 2.3.4.4g and 2.3.4.4h. Interestingly, we did not detect viruses of clade 2.3.4.4b, which in recent years have dominated in different parts of the world. The viruses isolated in this current study do not encode major determinants of mammalian adaptation (e.g., PB2-E627K or PB1-D701N) but possess amino acid substitutions that may affect viral receptor-binding, replication, or the responses to human antiviral factors. Several of the highly pathogenic H5N6 virus samples contained other influenza viruses, providing an opportunity for reassortment. Collectively, our study demonstrates that the highly pathogenic H5 viruses circulating in Vietnam in 2018-2021 were different from those in other parts of the world, and that the Vietnamese H5 viruses continue to evolve through mutations and reassortment.

Sequential intrahost evolution and onward transmission of SARS-CoV-2 variants.

Persistent severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections have been reported in immune-compromised individuals and people undergoing immune-modulatory treatments. Although intrahost evolution has been documented, direct evidence of subsequent transmission and continued stepwise adaptation is lacking. Here we describe sequential persistent SARS-CoV-2 infections in three individuals that led to the emergence, forward transmission, and continued evolution of a new Omicron sublineage, BA.1.23, over an eight-month period. The initially transmitted BA.1.23 variant encoded seven additional amino acid substitutions within the spike protein (E96D, R346T, L455W, K458M, A484V, H681R, A688V), and displayed substantial resistance to neutralization by sera from boosted and/or Omicron BA.1-infected study participants. Subsequent continued BA.1.23 replication resulted in additional substitutions in the spike protein (S254F, N448S, F456L, M458K, F981L, S982L) as well as in five other virus proteins. Our findings demonstrate not only that the Omicron BA.1 lineage can diverge further from its already exceptionally mutated genome but also that patients with persistent infections can transmit these viral variants. Thus, there is, an urgent need to implement strategies to prevent prolonged SARS-CoV-2 replication and to limit the spread of newly emerging, neutralization-resistant variants in vulnerable patients.

Multi-omic longitudinal study reveals immune correlates of clinical course among hospitalized COVID-19 patients.

The IMPACC cohort, composed of >1,000 hospitalized COVID-19 participants, contains five illness trajectory groups (TGs) during acute infection (first 28 days), ranging from milder (TG1-3) to more severe disease course (TG4) and death (TG5). Here, we report deep immunophenotyping, profiling of >15,000 longitudinal blood and nasal samples from 540 participants of the IMPACC cohort, using 14 distinct assays. These unbiased analyses identify cellular and molecular signatures present within 72 h of hospital admission that distinguish moderate from severe and fatal COVID-19 disease. Importantly, cellular and molecular states also distinguish participants with more severe disease that recover or stabilize within 28 days from those that progress to fatal outcomes (TG4 vs. TG5). Furthermore, our longitudinal design reveals that these biologic states display distinct temporal patterns associated with clinical outcomes. Characterizing host immune responses in relation to heterogeneity in disease course may inform clinical prognosis and opportunities for intervention.

Genomic and ultrastructural analysis of monkeypox virus in skin lesions and in human/animal infected cells reveals further morphofunctional insights into viral pathogenicity.

Monkeypox (MPOX) is a zoonotic disease that affects humans and other primates, resulting in a smallpox-like illness. It is caused by monkeypox virus (MPXV), which belongs to the Poxviridae family. Clinically manifested by a range of cutaneous and systemic findings, as well as variable disease severity phenotypes based on the genetic makeup of the virus, the cutaneous niche and respiratory mucosa are the epicenters of MPXV pathogenicity. Herein, we describe the ultrastructural features of MPXV infection in both human cultured cells and cutaneous clinical specimens collected during the 2022-2023 MPOX outbreak in New York City that were revealed through electron microscopy. We observed typical enveloped virions with brick-shaped morphologies that contained surface protrusions, consistent with the classic ultrastructural features of MPXV. In addition, we describe morpho-functional evidence that point to roles of distinct cellular organelles in viral assembly during clinical MPXV infection. Interestingly, in skin lesions, we found abundant melanosomes near viral assembly sites, particularly in the vicinity of mature virions, which provides further insight into virus-host interactions at the subcellular level that contribute to MPXV pathogenesis. These findings not only highlight the importance of electron microscopic studies for further investigation of this emerging pathogen but also in characterizing MPXV pathogenesis during human infection.