Abstract

Background

In May 2022, mpox cases were reported in nonendemic countries, including the United States. We examined mpox infections in the Veterans Health Administration (VHA).

Methods

Mpox diagnostic and whole genome sequencing (WGS) results, demographics, risk factors, hospitalizations, exposures, deaths, and pharmacy and immunization data were obtained from VHA data sources (23 May 2022–31 May 2023).

Results

Of 1144 Veterans tested, 251 (21.9%) were presumptive positive for nonvariola orthopoxvirus (NVO) or confirmed positive for NVO and Monkeypox virus (MPXV). Incidence rate was 7.5 per 100 000 Veterans in care, with the highest rate observed in Veterans aged 25–34 years (13.83 cases per 100 000). Higher odds of NVO or NVO/MPXV positivity was associated with male sex; non-Hispanic Black race/ethnicity; syphilis or human immunodeficiency virus (HIV) positivity; or genital/rectal sample site, whereas older age and vaccination with JYNNEOS or vaccinia (smallpox) had lower odds. Among 209 with confirmatory testing, 90.4% reported intimate contact and/or an epidemiological link, 84.5% were men who have sex with men (MSM), 24.2% received tecovirimat, and 8.1% were hospitalized with 1 death. Eighty-six sequenced samples had evaluable WGS results. All were clade IIb, representing 10 different lineages from 20 states and the District of Columbia.

Conclusions

Mpox affected younger, MSM, non-Hispanic Black, and HIV/syphilis-positive men among US Veterans. Viral diversity was noted across geographic regions. At-risk Veterans would benefit from vaccination and risk reduction strategies for mpox and other sexually transmitted infections.

Mpox is a viral zoonotic disease caused by Monkeypox virus (MPVX), an orthopoxvirus in the Poxviridae family and endemic to Western and Central Africa. A global outbreak was first recognized in May 2022 with cases reported concurrently in endemic and nonendemic countries, including the United States (US). The World Health Organization declared a Public Health Emergency of International Concern on 23 July 2022, followed by the US Department of Health and Human Services declaring a US Public Health Emergency on 4 August 2022 [1, 2].

There are 2 recognized genetic clades of the virus, with clade I originating in Central Africa and clade II associated with West Africa. Clade I is typically more transmissible and deadlier, while clade II tends to be more self-limited. Endemic spread is generally via animal–human transmission, while nonendemic spread is primarily via human-to-human contact [3, 4]. With human-to-human transmission, virus enters through broken skin, mucosal surfaces, or possibly via the respiratory tract [3]. MPXV can also be transmitted via contaminated objects (such as clothing) or sharps injuries in healthcare settings or community settings (such as tattoo parlors). Mpox is most often characterized by a self-limited illness; however, some cases may lead to severe illness and death. Patients typically present 1–21 days following exposure with symptoms that may include rash, enlarged lymph nodes, fever, sore throat, headache, muscle aches, back pain, and fatigue. Rash is typically seen on the face, mouth, throat, hands, feet, groin, and genital and/or anal region and may be painful and/or itchy, often evolving from macules and papules to vesicles and pustules that eventually crust over. Most infections are mild and resolve within 2–4 weeks after symptom onset, most often without requiring antiviral treatment [5]. However, symptoms may persist longer or progress to disseminated infections in individuals with weakened immune systems, especially persons with advanced human immunodeficiency virus (HIV) disease [6]. Concomitant sexually transmitted infections (STIs), particularly gonorrhea, chlamydia, and syphilis, have been reported in 17%–31% of those tested [7–10].

The global outbreak beginning in 2022 was caused mostly by clade IIb virus and with somewhat atypical clinical presentation in that the lesions often started in the groin, anus, or mouth, appearing before or concurrently with other symptoms and not always associated with spread to other parts of the body. Transmission dynamics were also different compared to historical reports, with spread most often via sexual contact and occasionally via close household contact or sharps injuries [11–14]. Initial cases in the US were associated with international travel, but within a short time period, domestic transmission was reported in all 50 US states, the District of Columbia (DC), and Puerto Rico. Most cases reported to the Centers for Disease Control and Prevention (CDC) were among men who have sex with men (MSM), aged 21–55 years, and disproportionately affecting persons with HIV and those from racial and ethnic minority groups, including Black and Hispanic/Latino persons. Additionally, most fatal cases occurred in immunocompromised Black persons [13, 15]. In the first year of the outbreak, >30 000 cases occurred, 1.2 million JYNNEOS (Bavarian Nordic, Hellerup, Denmark) third-generation smallpox vaccines were administered, and >6900 individuals were treated with the antiviral medication tecovirimat (TPOXX, SIGA Technologies, New York, New York) in the US via an expanded-access, Investigational New Drug protocol [15].

US Veterans represent a potentially high-risk population for MPXV infections, including severe infections, given most enrollees in the Veterans Health Administration (VHA) are adult males, with higher rates of underlying comorbidities [16, 17] and because VHA is the largest single provider of HIV care in the US [18]. As the largest integrated healthcare system in the US, VHA provides care to approximately 6.75 million patients annually across 1321 care sites nationwide and to >31 000 Veterans living with HIV [19, 20].

The objectives of this study were to identify risk factors for nonvariola orthopoxvirus (NVO)/MPXV test positivity by comparing characteristics of NVO/MPXV-positive individuals with those who tested negative and to characterize the epidemiology of MPXV confirmed infections, including infection severity, hospitalization, and death as well as risk factors for infection exposure. In addition, we sought to evaluate VHA's immunization efforts since the start of the outbreak and summarize genetic characterization of samples submitted for whole genome sequencing (WGS) to VHA's Public Health Reference Laboratory (PHRL). Analyzing the impact of the mpox outbreak in VHA, particularly among racial and ethnic minority groups and individuals living with HIV, may help guide prevention, testing, treatment, and control strategies, including vaccination efforts for mpox.

METHODS

Data Collection and Epidemiologic Analysis

VHA inpatient and outpatient medical encounters are maintained in a nationwide database of electronic health records (EHRs), known as Veterans Information Systems and Technology Architecture (VistA). All data for this study were obtained from the Corporate Data Warehouse (a repository and health data warehouse comprised of VistA clinical data and other data systems) and VHA's Praedico Public Health Surveillance System (Bitscopic, Los Angeles, California), which compiles public health data from VistA across multiple EHR domains, including diagnostic testing. The time frame for laboratory samples included in this study was 23 May 2022 through 24 May 2023.

For comparison of patients presumptive positive/confirmed positive for NVO/MPXV real-time polymerase chain reaction (qPCR) to those who were NVO qPCR negative, demographic variables extracted from the EHR included age, birth sex, race, ethnicity, facility state location, rurality of residence (based on Rural-Urban Commuting Area [RUCA] codes) [21], sample site, JYNNEOS or vaccinia (smallpox) vaccination status, HIV or STI status, and most recent CD4 count and percentage for those who were HIV positive prior to mpox diagnosis (obtained starting from 1 January 2021). HIV-positive status was defined as presence of an International Classification of Diseases, Clinical Modification, Ninth or Tenth Revision (ICD-9 or ICD-10 diagnosis code for HIV [042, V08, 079.53; B20, B97.35, O98.7*, Z21]) during any outpatient or inpatient encounter from 1 January 2000, through the study time frame. STI-positive status was defined as presence of an ICD-10 diagnosis code for syphilis, gonorrhea, chlamydia, or anogenital herpes (A51*, A52*, A53*, A54*, A55*, A60*) during any outpatient or inpatient encounter from 1 March 2022 through the study time frame. Individuals were determined to have been vaccinated for smallpox/mpox if they fell into any of the following categories: (1) documented receipt of 1 or more JYNNEOS vaccine doses prior to NVO/MPXV test date; (2) documented receipt of vaccinia (smallpox) vaccine prior to NVO/MPXV test date; or (3) birth date prior to 1972 and therefore assumed to have been vaccinated with vaccinia vaccine. Individuals who either had no documented receipt of vaccination or received vaccination after NVO/MPXV test date (and were not born before 1972) were considered unvaccinated. Incidence rates for NVO/MPXV positivity were calculated by age group per 100 000 Veterans in care, based on the VHA Support Service Center Capital Assets Unique Patients data cube.

The distributions of characteristics at the time of testing were compared between patients who tested positive for NVO/MPXV and those who tested NVO negative using Pearson χ2 test for categorical variables and Mann–Whitney–Wilcoxon test for continuous variables. We then performed binary logistic regression analysis using R version 4.3.0 (R Foundation for Statistical Computing, Vienna, Austria) to estimate the association of potentially relevant covariates with test positivity, adjusting for age, birth sex, race/ethnicity, rurality of residence, sample anatomic site, HIV, STI, and vaccination status, and expressed results as adjusted odds ratios (AORs) with 95% confidence intervals (CIs).

For the subset of patients with confirmatory MPXV testing, we additionally collected outpatient and inpatient encounters with monkeypox (B04)/orthopoxvirus (B08.09) diagnosis codes and performed chart reviews for exposure history; sexual contact among gay, bisexual, and MSM; confirmation of HIV status; TPOXX receipt; and deaths.

Sample Testing and Whole Genome Sequencing

For patients who presented to VHA facilities, suspect lesions collected for testing were swabbed by clinical providers according to CDC instructions [22]. Samples were submitted for NVO qPCR screen testing to state or local public health laboratories according to local criteria or VHA PHRL using the CDC NVO Screen Real-Time PCR Assay as described. Samples could also be sent to commercial laboratories that used their in-house NVO or NVO/MPXV qPCR assays. We performed manual chart review to determine whether specimens ordered by VHA providers but processed outside VHA had confirmatory testing for MPXV either via scanned laboratory reports or provider documentation of confirmatory testing. Samples sent to PHRL were shipped frozen in blue or dry ice as dry swabs. At PHRL, viral nucleic acid was extracted using the MagnaPure 96 (Roche Diagnostics, Indianapolis, Indiana) automated nucleic acid extraction platform and a QuantStudio Dx (ThermoFisher Scientific, Waltham, Massachusetts) was used for qPCR as described in the CDC NVO instructions for use [23]. Samples positive with the CDC NVO qPCR assay were confirmed at PHRL using a laboratory-developed MPXV confirmatory qPCR assay analyzing clade II-West African Strain (WAS) and Generic Strain (GEN) tumor necrosis factor receptor (TNFR) gene as previously described [24]. Additional VHA patients who were cared for at non-VHA healthcare facilities were identified based on mpox diagnosis codes and/or facility issue briefs submitted to Veterans Affairs (VA) Central Office, and testing results were obtained via EHR review.

Samples confirmed at PHRL for MPXV that had cycle threshold (Ct) results ≤29 underwent MPXV WGS using previously extracted MPXV DNA used for the NVO screening assay and Nextera XT DNA Sample Preparation and Index Kits (Illumina, San Diego, California) to generate libraries for WGS on an MiSeq (Illumina) according to manufacturer instructions. Generated FASTQ files were subjected to a previously published analytic pipeline [25], and generated FASTA files were further analyzed using Nextclade to determine viral clade and lineage type [26, 27]. MPXV sequences were analyzed for single-nucleotide polymorphisms (SNPs) that defined the B.1 parent lineage and SNPs that defined B.1 sublineages that were distinct from the parent B.1 lineage [28]. If the nucleotide substitution(s) were present, the lineage was reported. Phylogenetic analysis was performed using Nextclade [27]. VA lineages were compared to type strains for clade IIb. VA lineage distribution and frequency was compared to other US lineages from sequences submitted to the Global Initiative on Sharing All Influenza Data (GISAID) [29]. Ct results were collected from samples tested at PHRL for orthopoxvirus DNA polymerase gene (E9L, VAC1) and human RNase P gene (RNP) using the NVO screening assay; and WAS and GEN TNFR genes from the MPXV confirmatory assay. Ct results from VAC1 were then compared to WAS and GEN among MPXV-positive samples, and RNP Ct results were compared between positive and negative samples using Z-test.

The data utilized in this study were obtained for the purpose of public health operations in VHA. No additional analyses outside of public health operational activities were performed. Therefore, this study was deemed to meet the requirements of public health surveillance as defined in 45 Code of Federal Regulations 46.102(I)(2). This project was approved by the Stanford University Institutional Review Board (Protocol ID 47191, “Public Health Surveillance in the Department of Veterans Affairs”) and written informed consent was waived.

RESULTS

We identified 1144 unique patients tested for NVO (and if presumptive positive using MPXV confirmatory testing) from 48 states/territories/District of Columbia, with 1055 (92.2%) male and median age of 46 (interquartile range [IQR], 35–60) years. Two hundred fifty-one (21.9%) were presumptive positive/positive, 802 (70.1%) were negative, and 91 (8%) equivocal or inconclusive by the NVO screening assay and not further tested, with positive tests peaking in August 2022 (Supplementary Figure 1). The age group 25–34 years had the highest incidence rate (13.83 cases per 100 000 unique patients in VHA care) (Supplementary Table 1).

The median age of those testing positive for NVO/MPXV was 41 (IQR, 33–49) years, and 51 (IQR, 36–63) years for those testing NVO negative (P < .001) (Table 1). Risk factors associated with higher odds of NVO/MPXV test positivity (compared to those testing negative) included male birth sex (AOR, 6.77 [95% CI, 2.03–42.06]), non-Hispanic Black race/ethnicity (AOR, 1.98 [95% CI, 1.12–3.63]), genital or rectal sample site (AOR, 1.80 [95% CI, 1.26–2.58]), HIV-positive status (AOR, 2.43 [95% CI, 1.68–3.52]), and syphilis-positive status (AOR, 1.83 [95% CI, 1.25–2.69]). Residence in urban areas and gonorrhea/chlamydia STI diagnoses were not significantly associated when adjusted for age, birth sex, race/ethnicity, sample site, vaccination status, and other diagnoses in the logistic regression model. Among individuals with positive HIV status, median CD4 count did not differ between those testing positive for NVO/MPXV (586.5 [IQR, 406.75–778.25] cells/µL) compared to those testing negative (590 [IQR, 404–768] cells/µL) (P = .905). Vaccination with JYNNEOS or vaccinia (smallpox) vaccine had lower odds of NVO/MPXV test positivity (AOR, 0.33 [95% CI, .22–.48]). Breakdown of vaccination status by NVO negative or NVO/MPXV presumptive/confirmed positive test result is detailed in Supplementary Table 2. Overall, from July 2022 through May 2023, >8000 JYNNEOS vaccination doses were administered in VHA to >4400 unique Veterans.

Table 1.

Characteristics of Patients Testing Positive or Negative for Nonvariola Orthopoxvirus or Monkeypox virus, and Risk Factors Associated With Test Positivity, Veterans Health Administration, May 2022–May 2023a

VariableNVO/MPXV PositiveNVO/MPXV NegativeP ValuebAdjusted Odds Ratio (95% CI)cP Valuec
Total251 (23.8)802 (76.2)
Age at diagnosis, y, median (IQR)41 (33–49)51 (36–63)<.0010.98 (.96–.99).004
Birth sex
 Male249 (99.2)726 (90.5)<.0016.77 (2.03–42.06).009
Race/ethnicity
 Hispanic/Latino27 (10.8)67 (8.4).2441.94 (.94–4.07).075
 Non-Hispanic Black139 (55.4)302 (37.7)<.0011.98 (1.12–3.63).022
 Non-Hispanic White65 (25.9)336 (41.9)<.0011.24 (.68–2.32).493
 Other or unknownd20 (8.0)93 (11.6).071
Rurality of residencee
 Urban residence231 (92.0)681 (84.9).0041.17 (.67–2.13).584
Sample site
 Genital or rectal sample site89 (35.5)140 (17.5)<.0011.80 (1.26–2.58).001
Vaccination statusf
 Vaccinated65 (25.9)520 (64.8)<.0010.33 (.22–.48)<.001
HIV or STI diagnosesg
 Anogenital herpes25 (10.0)53 (6.6).077
 Chlamydia24 (9.6)27 (3.4)<.0011.02 (.51–2.03).962
 Gonorrhea30 (12.0)34 (4.2)<.0011.17 (.62–2.22).626
 Syphilis103 (41.0)122 (15.2)<.0011.83 (1.25–2.69).002
 HIV145 (57.8)182 (22.7)<.0012.43 (1.68–3.52)<.001
CD4 count among HIV positive, cells/μLh, median (IQR)586.5 (406.75–778.25)590 (404–768).905
VariableNVO/MPXV PositiveNVO/MPXV NegativeP ValuebAdjusted Odds Ratio (95% CI)cP Valuec
Total251 (23.8)802 (76.2)
Age at diagnosis, y, median (IQR)41 (33–49)51 (36–63)<.0010.98 (.96–.99).004
Birth sex
 Male249 (99.2)726 (90.5)<.0016.77 (2.03–42.06).009
Race/ethnicity
 Hispanic/Latino27 (10.8)67 (8.4).2441.94 (.94–4.07).075
 Non-Hispanic Black139 (55.4)302 (37.7)<.0011.98 (1.12–3.63).022
 Non-Hispanic White65 (25.9)336 (41.9)<.0011.24 (.68–2.32).493
 Other or unknownd20 (8.0)93 (11.6).071
Rurality of residencee
 Urban residence231 (92.0)681 (84.9).0041.17 (.67–2.13).584
Sample site
 Genital or rectal sample site89 (35.5)140 (17.5)<.0011.80 (1.26–2.58).001
Vaccination statusf
 Vaccinated65 (25.9)520 (64.8)<.0010.33 (.22–.48)<.001
HIV or STI diagnosesg
 Anogenital herpes25 (10.0)53 (6.6).077
 Chlamydia24 (9.6)27 (3.4)<.0011.02 (.51–2.03).962
 Gonorrhea30 (12.0)34 (4.2)<.0011.17 (.62–2.22).626
 Syphilis103 (41.0)122 (15.2)<.0011.83 (1.25–2.69).002
 HIV145 (57.8)182 (22.7)<.0012.43 (1.68–3.52)<.001
CD4 count among HIV positive, cells/μLh, median (IQR)586.5 (406.75–778.25)590 (404–768).905

Data are presented as No. (%) unless otherwise indicated. Bold values denote statistical significance at the P < .05 level.

Abbreviations: CI, confidence interval; HIV, human immunodeficiency virus; IQR, interquartile range; MPXV, Monkeypox virus; NVO, nonvariola orthopoxvirus; STI, sexually transmitted infection.

aPatients testing equivocal or inconclusive for NVO/MPXV were excluded from this analysis (n = 91).

bPearson χ2 test for categorical variables and Mann–Whitney Wilcoxon test for continuous variables.

cBinary logistic regression analysis.

dIncludes non-Hispanic American Indian or Alaska Native, non-Hispanic Asian, non-Hispanic Native Hawaiian or other Pacific Islander, and those with unknown race/ethnicity.

eRurality of residence based on Rural-Urban Commuting Area (RUCA) codes.

fIncludes vaccination with JYNNEOS 1 or 2 doses, or vaccinia vaccine, any time prior to mpox test date. Vaccinations given after mpox test date, or individuals with no recorded vaccination, were excluded.

gHIV and STI diagnosis status based on International Classification of Diseases diagnosis codes.

hMost recent CD4 count prior to mpox diagnosis.

Table 1.

Characteristics of Patients Testing Positive or Negative for Nonvariola Orthopoxvirus or Monkeypox virus, and Risk Factors Associated With Test Positivity, Veterans Health Administration, May 2022–May 2023a

VariableNVO/MPXV PositiveNVO/MPXV NegativeP ValuebAdjusted Odds Ratio (95% CI)cP Valuec
Total251 (23.8)802 (76.2)
Age at diagnosis, y, median (IQR)41 (33–49)51 (36–63)<.0010.98 (.96–.99).004
Birth sex
 Male249 (99.2)726 (90.5)<.0016.77 (2.03–42.06).009
Race/ethnicity
 Hispanic/Latino27 (10.8)67 (8.4).2441.94 (.94–4.07).075
 Non-Hispanic Black139 (55.4)302 (37.7)<.0011.98 (1.12–3.63).022
 Non-Hispanic White65 (25.9)336 (41.9)<.0011.24 (.68–2.32).493
 Other or unknownd20 (8.0)93 (11.6).071
Rurality of residencee
 Urban residence231 (92.0)681 (84.9).0041.17 (.67–2.13).584
Sample site
 Genital or rectal sample site89 (35.5)140 (17.5)<.0011.80 (1.26–2.58).001
Vaccination statusf
 Vaccinated65 (25.9)520 (64.8)<.0010.33 (.22–.48)<.001
HIV or STI diagnosesg
 Anogenital herpes25 (10.0)53 (6.6).077
 Chlamydia24 (9.6)27 (3.4)<.0011.02 (.51–2.03).962
 Gonorrhea30 (12.0)34 (4.2)<.0011.17 (.62–2.22).626
 Syphilis103 (41.0)122 (15.2)<.0011.83 (1.25–2.69).002
 HIV145 (57.8)182 (22.7)<.0012.43 (1.68–3.52)<.001
CD4 count among HIV positive, cells/μLh, median (IQR)586.5 (406.75–778.25)590 (404–768).905
VariableNVO/MPXV PositiveNVO/MPXV NegativeP ValuebAdjusted Odds Ratio (95% CI)cP Valuec
Total251 (23.8)802 (76.2)
Age at diagnosis, y, median (IQR)41 (33–49)51 (36–63)<.0010.98 (.96–.99).004
Birth sex
 Male249 (99.2)726 (90.5)<.0016.77 (2.03–42.06).009
Race/ethnicity
 Hispanic/Latino27 (10.8)67 (8.4).2441.94 (.94–4.07).075
 Non-Hispanic Black139 (55.4)302 (37.7)<.0011.98 (1.12–3.63).022
 Non-Hispanic White65 (25.9)336 (41.9)<.0011.24 (.68–2.32).493
 Other or unknownd20 (8.0)93 (11.6).071
Rurality of residencee
 Urban residence231 (92.0)681 (84.9).0041.17 (.67–2.13).584
Sample site
 Genital or rectal sample site89 (35.5)140 (17.5)<.0011.80 (1.26–2.58).001
Vaccination statusf
 Vaccinated65 (25.9)520 (64.8)<.0010.33 (.22–.48)<.001
HIV or STI diagnosesg
 Anogenital herpes25 (10.0)53 (6.6).077
 Chlamydia24 (9.6)27 (3.4)<.0011.02 (.51–2.03).962
 Gonorrhea30 (12.0)34 (4.2)<.0011.17 (.62–2.22).626
 Syphilis103 (41.0)122 (15.2)<.0011.83 (1.25–2.69).002
 HIV145 (57.8)182 (22.7)<.0012.43 (1.68–3.52)<.001
CD4 count among HIV positive, cells/μLh, median (IQR)586.5 (406.75–778.25)590 (404–768).905

Data are presented as No. (%) unless otherwise indicated. Bold values denote statistical significance at the P < .05 level.

Abbreviations: CI, confidence interval; HIV, human immunodeficiency virus; IQR, interquartile range; MPXV, Monkeypox virus; NVO, nonvariola orthopoxvirus; STI, sexually transmitted infection.

aPatients testing equivocal or inconclusive for NVO/MPXV were excluded from this analysis (n = 91).

bPearson χ2 test for categorical variables and Mann–Whitney Wilcoxon test for continuous variables.

cBinary logistic regression analysis.

dIncludes non-Hispanic American Indian or Alaska Native, non-Hispanic Asian, non-Hispanic Native Hawaiian or other Pacific Islander, and those with unknown race/ethnicity.

eRurality of residence based on Rural-Urban Commuting Area (RUCA) codes.

fIncludes vaccination with JYNNEOS 1 or 2 doses, or vaccinia vaccine, any time prior to mpox test date. Vaccinations given after mpox test date, or individuals with no recorded vaccination, were excluded.

gHIV and STI diagnosis status based on International Classification of Diseases diagnosis codes.

hMost recent CD4 count prior to mpox diagnosis.

Among 209 with positive confirmatory MPXV testing (from PHRL and other testing laboratories) from 29 states plus DC, 207 (99.0%) were birth sex male (includes 2 transgender women), primarily in the age group 25–44 years (131 [62.7%]). The majority were non-Hispanic Black race/ethnicity (108 [51.7%]) and 24 (11.5%) were Hispanic/Latino ethnicity. California, Texas, and DC had the highest number of confirmed positives (45, 26, and 22, respectively). VHA outpatient mpox ICD-10–coded encounter(s) were present in 141 (67.5%), and 189 (90.4%) reported intimate contact exposure and/or epidemiological link. Additionally, 175 (84.5%) were MSM and 113 (54.1%) were HIV positive. Besides HIV, 109 (52.2%) had 1 or more ICD-10–coded STIs, including syphilis (40.7%), herpesvirus (10.1%), gonococcal (12.0%), and chlamydial (9.1%) infections.

Treatment of confirmed MPXV-positive individuals with TPOXX in VHA was identified in 51 (24.4%) and 17 (8.1%) were hospitalized with an assigned B04/B08.09 diagnosis code. Among those hospitalized, only 2 had history of vaccination (in 2005 and 2007, respectively); 13 received TPOXX. There was 1 death in a patient with end-stage HIV/AIDS with multiple comorbidities and opportunistic infections, including mpox. The patient had a single documented vaccinia (smallpox) vaccine administered while serving in the military more than 16 years prior to mpox infection. There was no additional documentation of more recent receipt of JYNNEOS vaccine prior to death.

Among samples from 824 patients tested at PHRL (thereby providing Ct values for anatomic collection site comparisons and sequencing that were otherwise not accessible or available from other testing sites) using the NVO screening assay, 178 (21.6%) were NVO presumptive positive. Of these, 173 were confirmed MPXV positive, 4 were MPXV negative, and 1 was not tested for MPXV. Of the 4 patient samples that were NVO presumptive positive and MPXV negative, the confirmatory MPXV GEN and WAS qPCR Ct values were negative. When Ct for samples using the NVO screening assay was compared to the Ct using the MPXV confirmatory assays, there was no significant difference comparing VAC1 (mean, 24.4 [range, 11.9–37]) to GEN (mean, 24.0 [range, 13.5–39.7]; P = .37) or WAS (mean, 25.0 [range, 15.2–40]; P = .31), respectively. There was a small but significant difference in RNP Ct comparing NVO presumptive positive (mean, 26.9 [range, 19.5–39.7]) and negative (mean, 29.8 [range, 0–39.31]) samples (P < .01).

No significant difference was noted comparing the NVO Ct values of genital or rectal samples (mean, 24.4 [range, 14.7–37.0]) to those of nongenital/rectal samples (mean, 24.8 [range, 15.3–38.33]) (P = .71). No significant difference was also noted comparing the Ct values of MPXV GEN and WAS between the anatomic groups described.

Of 173 MPXV confirmed samples, 118 (68%) samples based on Ct threshold <29 were sequenced. We observed that 19 of 27 (70%) and 5 of 91 (5%) with whole genome coverage of <70% or >70%, respectively, did not have all the required SNPs for lineage designation. Therefore, using a conservative approach, we considered that 86 of 91 patient samples with >70% whole genome coverage had evaluable WGS results based on specific SNPs that defined lineages. All MPXV variants were clade IIb, representing 10 different lineages from 20 states and DC (Table 2 and Supplementary Figure 2). The B.1.2 variant was most common (33) detected in 13 states, followed by B.1 (28) in 14 states. In our cohort, all B.1 sublineages were contemporaneously represented with B.1.

Table 2.

Mpox Virus Variant Type and Number Among Veteran Patients by Testing Facility State, Veterans Health Administration, May 2022–May 2023

StateTotal PatientsVariants
B.1B.1.1B.1.2B.1.3B.1.7B.1.8B.1.11B.1.12B.1.13B.1.17
Alabama10010000000
Arizona11000000000
Arkansas10001000000
California166321111001
Colorado10010000000
Florida61020111000
Georgia101090000000
Illinois43000000010
Indiana10010000000
Michigan21010000000
Mississippi10001000000
Missouri22000000000
New Jersey11000000000
New York21000000100
North Carolina60022200000
Ohio11000000000
Pennsylvania31010100000
South Carolina20020000000
Texas126060000000
Virginia41012000000
District of Columbia92042100000
Totals86283339622111
StateTotal PatientsVariants
B.1B.1.1B.1.2B.1.3B.1.7B.1.8B.1.11B.1.12B.1.13B.1.17
Alabama10010000000
Arizona11000000000
Arkansas10001000000
California166321111001
Colorado10010000000
Florida61020111000
Georgia101090000000
Illinois43000000010
Indiana10010000000
Michigan21010000000
Mississippi10001000000
Missouri22000000000
New Jersey11000000000
New York21000000100
North Carolina60022200000
Ohio11000000000
Pennsylvania31010100000
South Carolina20020000000
Texas126060000000
Virginia41012000000
District of Columbia92042100000
Totals86283339622111

Bold values indicates sample number of lineage or sublineage present for US state or District of Columbia.

Table 2.

Mpox Virus Variant Type and Number Among Veteran Patients by Testing Facility State, Veterans Health Administration, May 2022–May 2023

StateTotal PatientsVariants
B.1B.1.1B.1.2B.1.3B.1.7B.1.8B.1.11B.1.12B.1.13B.1.17
Alabama10010000000
Arizona11000000000
Arkansas10001000000
California166321111001
Colorado10010000000
Florida61020111000
Georgia101090000000
Illinois43000000010
Indiana10010000000
Michigan21010000000
Mississippi10001000000
Missouri22000000000
New Jersey11000000000
New York21000000100
North Carolina60022200000
Ohio11000000000
Pennsylvania31010100000
South Carolina20020000000
Texas126060000000
Virginia41012000000
District of Columbia92042100000
Totals86283339622111
StateTotal PatientsVariants
B.1B.1.1B.1.2B.1.3B.1.7B.1.8B.1.11B.1.12B.1.13B.1.17
Alabama10010000000
Arizona11000000000
Arkansas10001000000
California166321111001
Colorado10010000000
Florida61020111000
Georgia101090000000
Illinois43000000010
Indiana10010000000
Michigan21010000000
Mississippi10001000000
Missouri22000000000
New Jersey11000000000
New York21000000100
North Carolina60022200000
Ohio11000000000
Pennsylvania31010100000
South Carolina20020000000
Texas126060000000
Virginia41012000000
District of Columbia92042100000
Totals86283339622111

Bold values indicates sample number of lineage or sublineage present for US state or District of Columbia.

DISCUSSION

Among a national cohort of VHA patients who were tested for NVO/MPXV during the 2022–2023 US outbreak, risk factors having the highest odds for NVO/MPXV positivity included male birth sex, HIV positivity, non-Hispanic Black race/ethnicity, syphilis diagnosis, and genital/rectal specimen collection. JYNNEOS/vaccinia vaccination and older age were protective factors. CDC reports that the highest-risk race/ethnicity groups in the US were non-Hispanic Black and Hispanic persons [15], while our study found only non-Hispanic Blacks at significant risk. Like other studies [7, 30], we found that NVO/MPXV positivity was associated with HIV and syphilis positivity. As a large national healthcare system, VHA surveillance capability was enhanced by the availability of data on patients who presented with mpox signs/symptoms but tested negative, which is an advantage compared to other prior evaluations and reporting systems.

Our finding of sample sites taken from genital or rectal lesions being associated with NVO/MPXV positivity was not surprising as Ct values from these sites have been shown to be lower than from other body sites [31], indicating that proper sampling from high-yield body sites may decrease false-negative or equivocal/inconclusive results. However, our results are different from controlled studies documenting the appearance, duration, and swabbing of rash/lesions and corresponding Ct values. These controlled studies demonstrated lower Ct values from skin followed by rectal swab sites, with Ct levels related to lesion duration [31, 32]. Ultimately, the Ct values collected at the ideal time and anatomic sites could be used to determine the MPXV infectivity. A study from Israel described a correlation between Cq (quantification cycle like Ct) values and in vitro infectivity of African green monkey kidney cell line (BSC-1). Samples with Cq > 35 resulted in either no or minimal infection [33].

Although JYNNEOS vaccination occurred significantly less among those with an mpox diagnosis in a recent study [34], documentation of JYNNEOS vaccination in our study was very low in both NVO/MPXV positive and negative patients. Prior vaccination with vaccinia (smallpox) vaccine is thought to provide little protection against mpox if given >3 years before exposure. Our data pose an interesting question regarding the possibility that in our large cohort of individuals vaccinated with vaccinia (during their military service and/or due to being born prior to 1972), protection against mpox may possibly have been conferred by vaccinia vaccination rather than JYNNEOS (Supplementary Table 3). This phenomenon was reported in a recent systematic review [35], and an additional recent study confirmed the existence of long-term cross-neutralizing antibodies elicited by prior vaccinia vaccination and seen in higher levels among individuals born before 1980 [36].

Among our confirmed MPXV cases, we found a high percentage with intimate contact exposure histories, MSM, HIV, and STIs documented in medical records. These findings align with previous studies that have reported high proportion of MPXV-infected individuals among adult MSM, those living with HIV, and from specific racial and ethnic minority groups [13, 15]. Other studies have also documented high rates of concomitant STIs in patients with mpox [7–10]. High rates of HIV, MSM, and STIs among mpox cases suggest the importance of testing for HIV and STIs, as well as discussing vaccination, HIV preexposure prophylaxis, and STI prevention strategies, where applicable, for patients who present with possible mpox.

Although VHA had relatively few mpox cases, testing in VHA peaked in August 2022 with the highest incidence rate observed among Veterans aged 25–34 years (13.83 per 100 000 Veterans in care). While the VHA peak was consistent with that reported by CDC, the incidence rate for Veterans aged 18–64 years (7.5 per 100 000 Veterans in care) was lower than published US population rate of 13.5 per 100 000 for persons aged 15–64 years [15, 37]. TPOXX administration was documented in 24.4% of VHA-confirmed mpox cases compared with a 22.8% administration rate for the US overall through 11 May 2023 [15]. Our hospitalization rate of 8.1% was similar to albeit slightly higher than a recent meta-analysis of global mpox reports, which found a hospitalization rate of approximately 7.3% and a reported hospitalization rate of 7.5% among unvaccinated individuals from 29 US jurisdictions [38, 39].

Our study has several limitations. Our sample of NVO/MPXV positives was small and our Veteran population is older and predominantly male, which may not represent the general population. Race/ethnicity data were missing in approximately 8% of our patients, which could have impacted our results. Some NVO presumptive positive individuals were not confirmed using a specific MPXV qPCR assay (or documentation of confirmatory testing was unavailable in the EHR). Some patients with mild, self-limited illness may not have presented for care or may have presented with epidemiological link/exposure history but before development of rash, so they were not tested. Poor sampling/swabbing technique may have led to equivocal or inconclusive test results, which was seen in 8% of VHA samples. Documentation of sample anatomic site was not available in 38% of samples (Supplementary Table 3). Early in the epidemic, some facilities were not aware that testing could be performed within VHA, and patients were referred to health departments or non-VHA facilities for testing. Our data sources do not contain clinical records for patients who received care in non-VHA healthcare settings unless these services were ordered and paid for by VHA. For these reasons, our data likely undercapture mpox in VHA. Similarly, some VHA patients were referred to their local health department or other facilities to obtain TPOXX and JYNNEOS vaccine until VHA facilities were able to obtain it locally. Therefore, MPXV testing, JYNNEOS vaccination, and ICD-coded mpox encounters occurring under these circumstances outside VHA were not represented in our study.

Only samples received by PHRL had WGS performed. More than 40% of the samples had high Ct values and were not sequenced, limiting our analysis of variants. However, 97.5% of nonsequenced samples came from the same 21 states and DC and generally had the same proportions of patients by state as those that were sequenced (Supplementary Table 4). In addition, we did not include samples that were sequenced and had low coverage (<70%) or in a few with >90% coverage, even though a designation was made by Nextclade, because we were not confident that a variant designation was accurate if lineage-specific SNPs were missing. Using the MPXV WGS metagenomic approach, we encountered occasional difficulty designating some sublineages, even with lower Ct values and high genome-wide coverage [40]. An amplicon-based MPXV WGS approach for improved results is now available [41]. Nevertheless we confirmed seeing many of the sublineages found in the US and represented in GISAID (Supplementary Table 5) and add to literature demonstrating contemporaneously geographic diversity of variants in the US among Veterans.

In conclusion, this report highlights the epidemiology of mpox among VHA enrollees during the 2022–2023 national outbreak. Demographic features of patients testing positive for NVO/MPXV were consistent with previous reports, including among HIV-positive, MSM, and non-Hispanic Black individuals. Timing and incidence rates of mpox among Veterans followed similar trends compared with reported national cases, particularly among younger age groups. As a national healthcare system report, this study provides a supplemental source of data on mpox, including observed MPXV lineages and viral diversity, and among patients who tested negative. Providers should continue to encourage vaccination among patients at highest risk for infection with mpox and test for mpox in patients with compatible exposure histories and documented risk factors, even in the absence of travel history to an endemic region. Additional, targeted public health messages, education, and prevention strategies are needed to ensure at-risk patients are receiving timely testing, treatment, and follow-up.

Supplementary Data

Supplementary materials are available at The Journal of Infectious Diseases online (http://jid.oxfordjournals.org/). Supplementary materials consist of data provided by the author that are published to benefit the reader. The posted materials are not copyedited. The contents of all supplementary data are the sole responsibility of the authors. Questions or messages regarding errors should be addressed to the author.

Notes

Acknowledgments. We thank Carlos Flores, PhD, and colleagues at the Genomics Division, Instituto Tecnológico y de Energías Renovables, Santa Cruz de Tenerife, Spain, and Dr Richard Neher at Nextclade for helpful discussions on mpox WGS procedures and analyses. We gratefully acknowledge all data contributors, that is, the authors and their originating laboratories responsible for obtaining the specimens, and their submitting laboratories for generating the genetic sequence and metadata and sharing via the GISAID Initiative, on which this research is based.

Author contributions. Conceptualization: C. L.-O., G. O., M. H. Methodology: C. L.-O., C. T., C. E., G. O., M. H. Validation: C. L.-O., G. O., P. S., M. H. Formal analysis: C. L.-O., G. O., C. T., T. E., C. E., M. H. Resources: M. H. Data curation: C. L.-O., G. O., C. E., C. T., T. E. Writing–original draft preparation: M. H., C. L.-O. Writing–review and editing: P. S., G. O., C. T., C. E., T. E., M. H. Supervision: M. H., C. L.-O. Project administration: M. H., C. L.-O. All authors have read and agreed to the final version of the manuscript.

Data availability. The data that support the findings of this study are available from the corresponding author upon reasonable request.

Disclaimer. The contents and opinions of this manuscript are those of the authors and do not necessarily represent those of the US government or the Department of Veterans Affairs.

Supplement sponsorship. This article appears as part of the supplement “Mpox: Challenges and Opportunities Following the Global 2022 Outbreak,” sponsored by the Centers for Disease Control and Prevention (Atlanta, GA).

Financial support. This work was supported by Department of Veterans Affairs intramural operational funding.

References

1

World Health Organization.
Second meeting of the International Health Regulations emergency committee regarding the multi-country outbreak of monkeypox. 2005. https://www.who.int/news/item/23-07-2022-second-meeting-of-the-international-health-regulations-(2005)-(ihr)-emergency-committee-regarding-the-multi-country-outbreak-of-monkeypox. Accessed 31 July 2023.

2

Administration for Strategic Preparedness and Response
. Determination that a public health emergency exists. 2022. https://aspr.hhs.gov/legal/PHE/Pages/monkeypox-4Aug22.aspx. Accessed 31 July 2023.

3

Sukhdeo
S
,
Mishra
S
,
Walmsley
S
.
Human monkeypox: a comparison of the characteristics of the new epidemic to the endemic disease
.
BMC Infect Dis
2022
;
22
:
928
.

4

Gessain
A
,
Nakoune
E
,
Yazdanpanah
Y
.
Monkeypox
.
N Engl J Med
2022
;
387
:
1783
93
.

5

Mitjà
O
,
Ogoina
D
,
Titanji
BK
, et al. 
Monkeypox
.
Lancet
2023
;
401
:
60
74
.

6

Miller
MJ
,
Cash-Goldwasser
S
,
Marx
GE
, et al. 
Severe monkeypox in hospitalized patients—United States, August 10–October 10, 2022
.
MMWR Morb Mortal Wkly Rep
2022
;
71
:
1412
7
.

7

Thornhill
JP
,
Barkati
S
,
Walmsley
S
, et al. 
Monkeypox virus infection in humans across 16 countries—April–June 2022
.
N Engl J Med
2022
;
387
:
679
91
.

8

Kim
M
,
Acevedo Martinez
E
,
Suárez Moscoso
NP
, et al. 
A retrospective study on 198 mpox cases in South Florida: clinical characteristics and outcomes with focus on human immunodeficiency virus status
.
Int J STD AIDS
2023
;
34
:
884
9
.

9

Patel
A
,
Bilinska
J
,
Tam
JCH
, et al. 
Clinical features and novel presentations of human monkeypox in a central London centre during the 2022 outbreak: descriptive case series
.
BMJ
2022
;
378
:
e072410
.

10

Tarín-Vicente
EJ
,
Alemany
A
,
Agud-Dios
M
, et al. 
Clinical presentation and virological assessment of confirmed human monkeypox virus cases in Spain: a prospective observational cohort study
.
Lancet
2022
;
400
:
661
9
.

11

Blackburn
D
,
Roth
NM
,
Gold
JAW
, et al. 
Epidemiologic and clinical features of mpox in transgender and gender-diverse adults—United States, May–November 2022
.
MMWR Morb Mortal Wkly Rep
2022
;
71
:
1605
9
.

12

Hennessee
I
,
Shelus
V
,
McArdle
CE
, et al. 
Epidemiologic and clinical features of children and adolescents aged <18 years with monkeypox—United States, May 17–September 24, 2022
.
MMWR Morb Mortal Wkly Rep
2022
;
71
:
1407
11
.

13

Kava
CM
,
Rohraff
DM
,
Wallace
B
, et al. 
Epidemiologic features of the monkeypox outbreak and the public health response—United States, May 17–October 6, 2022
.
MMWR Morb Mortal Wkly Rep
2022
;
71
:
1449
56
.

14

Oakley
LP
,
Hufstetler
K
,
O’Shea
J
, et al. 
Mpox cases among cisgender women and pregnant persons—United States, May 11–November 7, 2022
.
MMWR Morb Mortal Wkly Rep
2023
;
72
:
9
14
.

15

McQuiston
JH
,
Braden
CR
,
Bowen
MD
, et al. 
The CDC domestic mpox response—United States, 2022–2023
.
MMWR Morb Mortal Wkly Rep
2023
;
72
:
547
52
.

16

Rogers
WH
,
Kazis
LE
,
Miller
DR
, et al. 
Comparing the health status of VA and non-VA ambulatory patients: the Veterans’ health and medical outcomes studies
.
J Ambul Care Manage
2004
;
27
:
249
62
.

17

Eibner
C
,
Krull
H
,
Brown
KM
, et al. 
Current and projected characteristics and unique health care needs of the patient population served by the Department of Veterans Affairs
.
Rand Health Q
2016
;
5
:
13
.

18

US Department of Veterans Affairs
. VA HIV care. https://www.hiv.va.gov/about-index.asp. Accessed 31 July 2023.

19

National Center for Veterans Analysis and Statistics
. At-a-glance pocket card, fiscal year 2023, 3rd quarter. 2023. https://www.va.gov/vetdata/docs/pocketcards/fy2023q3.pdf. Accessed 31 July 2023.

20

US Department of Veterans Affairs
. Human immunodeficiency virus (HIV) fact sheet. 2020. https://www.hiv.va.gov/pdf/HIV-program-factsheet.pdf. Accessed 31 July 2023.

21

US Department of Agriculture
. Rural-urban commuting area codes. 2019. https://www.ers.usda.gov/data-products/rural-urban-commuting-area-codes/. Accessed 31 July 2023.

22

Centers for Disease Control and Prevention
. Guidelines for collecting and handling specimens for mpox testing. 2022. https://www.cdc.gov/poxvirus/mpox/clinicians/prep-collection-specimens.html. Accessed 31 July 2023.

23

Centers for Disease Control and Prevention
.
Test procedure: non-variola orthopoxvirus generic real-time PCR test.
2022
. https://www.cdc.gov/poxvirus/mpox/pdf/non-variola-orthopoxvirus-generic-real-time-pcr-test.pdf. Accessed 31 July 2023.

24

Li
Y
,
Zhao
H
,
Wilkins
K
,
Hughes
C
,
Damon
IK
.
Real-time PCR assays for the specific detection of Monkeypox virus West African and Congo Basin strain DNA
.
J Virol Methods
2010
;
169
:
223
7
.

25

Muñoz-Barrera
A
,
Ciuffreda
L
,
Alcoba-Florez
J
, et al. 
Bioinformatic approaches to draft the viral genome sequence of Canary Islands cases related to the multicountry mpox virus 2022—outbreak
.
Comput Struct Biotechnol J
2023
;
21
:
2197
203
.

26

Nextclade
. Clade assignment, mutation calling, and sequence quality checks. https://clades.nextstrain.org/. Accessed 31 July 2023.

27

Aksamentov
I
,
Roemer
C
,
Hodcroft
EB
,
Neher
RA
.
Nextclade: clade assignment, mutation calling and quality control for viral genomes
.
J Open Source Software
2021
;
6
:
3773
.

28

Nextclade
. MPVX-lineage designations. https://github.com/mpxv-lineages/lineage-designation/tree/master/lineages. Accessed 31 July 2023.

29

Khare
S
,
Gurry
C
,
Freitas
L
, et al. 
GISAID’s role in pandemic response
.
China CDC Wkly
2021
;
3
:
1049
51
.

30

Curran
KG
,
Eberly
K
,
Russell
OO
, et al. 
HIV and sexually transmitted infections among persons with monkeypox—eight U.S. jurisdictions, May 17–July 22, 2022
.
MMWR Morb Mortal Wkly Rep
2022
;
71
:
1141
7
.

31

Edman-Wallér
J
,
Jonsson
O
,
Backlund
G
,
Muradrasoli
S
,
Sondén
K
.
Results of PCR analysis of mpox clinical samples, Sweden, 2022
.
Emerg Infect Dis
2023
;
29
:
1220
2
.

32

Palich
R
,
Burrel
S
,
Monsel
G
, et al. 
Viral loads in clinical samples of men with monkeypox virus infection: a French case series
.
Lancet Infect Dis
2023
;
23
:
74
80
.

33

Paran
N
,
Yahalom-Ronen
Y
,
Shifman
O
, et al. 
Monkeypox DNA levels correlate with virus infectivity in clinical samples, Israel, 2022
.
Euro Surveill
2022
;
27
:
2200636
.

34

Deputy
NP
,
Deckert
J
,
Chard
AN
, et al. 
Vaccine effectiveness of JYNNEOS against mpox disease in the United States
.
N Engl J Med
2023
;
388
:
2434
43
.

35

Akter
F
,
Hasan
TB
,
Alam
F
, et al. 
Effect of prior immunisation with smallpox vaccine for protection against human mpox: a systematic review
.
Rev Med Virol
2023
;
33
:
e2444
.

36

Hubert
M
,
Guivel-Benhassine
F
,
Bruel
T
, et al. 
Complement-dependent mpox-virus-neutralizing antibodies in infected and vaccinated individuals
.
Cell Host Microbe
2023
;
31
:
937
948.e4
.

37

Zelaya
CE
,
Smith
BP
,
Riser
AP
, et al. 
Urban and rural mpox incidence among persons aged 15–64 years—United States, May 10–December 31, 2022
.
MMWR Morb Mortal Wkly Rep
2023
;
72
:
574
8
.

38

DeWitt
ME
,
Polk
C
,
Williamson
J
, et al. 
Global monkeypox case hospitalisation rates: a rapid systematic review and meta-analysis
.
EClinicalMedicine
2022
;
54
:
101710
.

39

Farrar
JL
,
Lewis
NM
,
Houck
K
, et al. 
Demographic and clinical characteristics of mpox in persons who had previously received 1 dose of JYNNEOS vaccine and in unvaccinated persons—29 U.S. jurisdictions, May 22–September 3, 2022
.
MMWR Morb Mortal Wkly Rep
2022
;
71
:
1610
5
.

40

Isidro
J
,
Borges
V
,
Pinto
M
, et al. 
Phylogenomic characterization and signs of microevolution in the 2022 multi-country outbreak of monkeypox virus
.
Nat Med
2022
;
28
:
1569
72
.

41

Chen
NFG
,
Chaguza
C
,
Gagne
L
, et al. 
Development of an amplicon-based sequencing approach in response to the global emergence of mpox
.
PLoS Biol
2023
;
21
:
e3002151
.

Author notes

Presented in part: IDWeek, Boston, Massachusetts, 11–15 October 2023. Abstract 899.

Potential conflicts of interest. All authors: No reported conflicts.

All authors have submitted the ICMJE Form for Disclosure of Potential Conflicts of Interest. Conflicts that the editors consider relevant to the content of the manuscript have been disclosed.

This work is written by (a) US Government employee(s) and is in the public domain in the US.

Supplementary data