A. Clinical History
This was a case of a 31-year-old Filipino, male, who consulted for pustules on the left ankle, the left gluteal area, the left upper lip, and one each on both knees. He had a history of recurrent folliculitis (one to two episodes per year). He had no known comorbidities, no known allergies to food or medications, and no previous Varicella or Measles infection. Vaccination against Varicella, Measles, and Smallpox was unrecalled. He was fully vaccinated against COVID-19 with two booster shots. He is a nonsmoker, an occasional alcoholic beverage drinker, and denied illicit drug use. He had a male sexual partner. He took a once-daily emtricitabine plus lamivudine combination tablet as human immunodeficiency virus pre-exposure prophylaxis (HIV PreP) for two months with the latest negative HIV test done just before PrEP initiation.
Four weeks before consultation, the patient began a 3-week-long leisure trip involving several European cities (Figure 1). He attended a concert and several social gatherings during his stay in Europe. He also interacted with stray animals at parks. Thirteen days before consultation (16 days from arrival in Europe), he had one episode of undocumented fever with chills, which resolved after taking a combined anti-inflammatory tablet. He returned to the Philippines seven days before consultation (21 days from arrival in Europe) with no symptoms. Six days before consultation, pruritic vesicular rashes appeared on the left gluteal area, bilateral knees, and left anterior ankle. Five days before consultation, he had bloody stools prompting consultation at a private hospital, and was diagnosed with external hemorrhoids. He was sent home with prescribed diosmin + hesperidin tablets. Four days prior to consult, he developed an ulcer on the right upper lip and itchiness on the surface of his tongue. He also noted the increasing size of the lesions on his knees which urged him to consult.
On physical examination, the patient was seen awake, conscious, ambulatory, and with stable vital signs. A focused physical examination of the skin revealed five well-defined, non-tender pustules with umbilication and erythematous borders on the left gluteal, bilateral knees, left ankle, and right upper lip (Figure 2). No lymphadenopathies were noted. The patient was found to have satisfied the criteria for Monkeypox Suspect hence admitted to an isolation room for further evaluation. Initial laboratory tests showed a WBC count of 8.6 x 109/L, neutrophils 53%, lymphocytes 35%, and the rest of the tests within normal limits including serum creatinine, SGPT, and SGOT. Several other possible causes of rash were ruled out: serologic measles Ag and Ab tests and herpes simplex virus (HSV-1 & 2) PCR of skin lesions were negative. Interestingly, Varicella IgM was positive but Varicella IgG and PCR of the skin lesions were negative for Varicella. The patient was also referred to the Dermatology service which performed skin punch biopsy, with results describing a viral infection. Specimens for monkeypox qPCR test and metagenomic sequencing were obtained and yielded positive results for monkeypox infection. The detailed qPCR and sequencing processes and results are described in the succeeding sections. The patient was discharged on the second hospital day. He was advised to keep the lesions clean and dry, using mild soap and moisturizing lotion daily as supportive management. Home isolation and daily monitoring of symptoms via online consultations were done until all crusts and scabs had completely disappeared. The patient’s isolation period lasted for 23 days without any serious complications (Figure 3). No identified close contact became symptomatic nor developed any rash.
B. Dermatopathologic Evaluation
A skin punch biopsy of the umbilicated pustule on the left gluteal area was done by the Dermatology service. The epidermis revealed scale crusts and ulceration. Some of the keratinocytes in the epidermis and the upper dermis were large with convoluted steel-gray nuclei and the dermis revealed red blood cell extravasation and a moderately dense, perivascular, and interstitial inflammatory infiltrate of lymphocytes, histiocytes, and plasma cells (Figure 4). These findings supported a diagnosis of a viral infection.
C. Specimen Collection and Processing
A total of nine tissue/lesion specimens and nine swab specimens obtained from three sites (right and left knees, and left ankle area) were sent to the Special Pathogens Laboratory for confirmatory monkeypox real-time polymerase chain reaction (RT-PCR) test.
C.1 Nucleic Acid Extraction from Dry Swab Sample
At the laboratory, single dry swab of cotton was separated from the swab base and was added to the lysis buffer solution (400 μL phosphate-buffered saline (1x PBS), 20 μL proteinase K and 400 μL buffer AL). The solution was homogenized by vigorous vortex-mixing and incubated at 56°C for 10 minutes. Total viral DNA was extracted from the sample using the QiaAmp DNA Mini Kit (Qiagen, Hilden, Germany, Cat No: 51306) according to the manufacturer’s instructions (7). The extracted DNA was eluted in 150 μL of nuclease-free water (NFW) and stored at -80°C until use.
C.2 Nucleic Acid Extraction from Tissue Sample
At the laboratory, 2.5mg (3mm in diameter) of the tissue sample was added to 80 μL phosphate-buffered saline (1x PBS) and was homogenized using a micropestle. The solution was briefly vortex-mixed after adding 100 μL buffer ATL and 20 μL proteinase K. Incubation at 56°C for 1 hour (vortex mixed every 20 minutes) and at 70°C for 10 minutes followed. Total viral DNA was extracted from the sample using the QiaAmp DNA Mini Kit (Qiagen, Hilden, Germany, Cat No: 51306) according to the manufacturer’s instructions (7). The extracted DNA was eluted in 200 μL of nuclease-free water (NFW) and stored at -80°C until use.
D. Real-Time Monkeypox PCR Testing and Results
The PCR primers and probes were developed from the sequences described by Li, et al (8). Probe-based real-time PCR assay was performed using Applied Biosystem’s AgPath-ID One Step PCR kit (4387424) (9) and Bio-Rad CFX96 Touch real-time PCR machine as PCR platform. RNase P was the assays’ internal target control. A total of 23 µL of master mix (6.5 µL nuclease-free water, 12.5 µL 2x buffer, 1 µL 25 RT-PCR enzyme mix, 1 µL each of forward and reverse primers, 1 µL probe) and 2 µL of the extracted DNA template was used. The following thermocycling protocols were programmed for monkeypox screening assay: 95 ºC for 10 minutes and 45 cycles of 95 ºC for 15 seconds, 60ºC for 20 seconds; for monkeypox differentiation assay: 95 ºC for 10 minutes and 45 cycles of 95 ºC for 15 seconds, 62ºC for 20 seconds. The optimized primer and probe concentrations were 10µM and 5µM, respectively for both screening and differentiation assays (Table 1).
Table 1. Primers and probes used in Monkeypox quantitative PCR test (8)
Primer/ Probe Name
|
Sequence 5’ to 3’
|
Assay Name
|
G2R_G Forward Primer
|
5'- GGAAAATGTAAAGACAACGAATACAG-3'
|
G2R_G
(Monkeypox Screening)
|
G2R_G Reverse Primer
|
5'- GCTATCACATAATCTGGAAGCGTA-3'
|
G2R_G Probe
|
5' FAM-AAGCCGTAATCTATGTTGTCTATCGTGTCC-3' BHQ1
|
G2R_WA Forward Primer
|
5'- CACACCGTCTCTTCCACAGA-3'
|
G2R_WA
(Monkeypox Differentiation)
|
G2R_WA Reverse Primer
|
5'- GATACAGGTTAATTTCCACATCG-3'
|
G2R_WA Probe
|
FAM 5'-AACCCGTCGTAACCAGCAATACATTT-3' BHQ1
|
C3L Forward Primer
|
5'- TGTCTACCTGGATACAGAAAGCAA-3'
|
C3L
(Monkeypox Differentiation)
|
C3L Reverse Primer
|
5'- GGCATCTCCGTTTAATACATTGAT-3'
|
C3L Probe
|
5' FAM-CCCATATATGCTAAATGTACCGGTACCGGA-3' BHQ1
|
Two sets of the collected lesion dry swab and lesion crust were confirmed to be positive for monkeypox using real-time PCR screening assay with a mean cycle threshold (Ct) value of 19.54. The monkeypox real-time PCR differentiation assay revealed that the same samples were positive for the Western African clade (Clade II) only with a mean Ct value of 19.85 (Table 2).
Table 2. Monkeypox RT-PCR result with mean cycle threshold (Ct) values
Target
|
Sample
|
CT Value
|
Mean CT Value
|
Result
|
G2R_G
|
Lesion Base Dry Swab Right Knee
|
20.68
|
19.54
|
Monkeypox Viral DNA Detected
|
Lesion Base Dry Swab Left Knee
|
19.09
|
Lesion Crust Right knee
|
17.63
|
Lesion Crust Left Ankle
|
20.76
|
G2R_WA
|
Lesion Base Dry Swab Right Knee
|
20.88
|
19.85
|
Monkeypox Clade II** Viral DNA Detected
|
Lesion Base Dry Swab Left Knee
|
19.32
|
Lesion Crust Right knee
|
18.18
|
Lesion Crust Left Ankle
|
21.03
|
C3L
|
Lesion Base Dry Swab Right Knee
|
-
|
-
|
Monkeypox Clade I* Viral DNA Not Detected
|
Lesion Base Dry Swab Left Knee
|
-
|
Lesion Crust Right knee
|
-
|
Lesion Crust Left Ankle
|
-
|
*formerly classified as Congo Basin clade
**formerly classified as West African clade
E. Metagenomic Sequencing of Monkeypox Specimens
To further characterize the etiologic agent of the first laboratory-confirmed Monkeypox case for the country, and properly classify its phylogenetic lineage, four specimens that tested positive for monkeypox RT-PCR were processed by the Molecular Biology Laboratory for metagenomic sequencing. Additionally, two confirmed target-negative samples were included in the sequencing run to be used as MPXV-negative specimen control. The standard Illumina DNA Prep protocol was followed (10). Fifty microliters (50 ul) of DNA extracts were purified using Agencourt AMPure XP (Beckman Coulter Genomics, IN, USA, Cat no. A63881), and quantified using Qubit HS dsDNA reagent kit (ThermoFisher Scientific, MA, USA Cat no. Q32854). The purified DNA extracts were normalized to 148 ng starting DNA input and were used for the Illumina DNA prep kit (Illumina, Inc., CA, USA) library preparation. Constructed libraries were quantified using a Qubit HS dsDNA reagent kit and Agilent D1000 ScreenTape system (Waldbronn, Germany). The samples that passed the QC criteria were pooled and subjected to shotgun metagenomic sequencing using the Illumina Miseq sequencing instrument. Out of the four specimens that tested positive for monkeypox real-time PCR, only three samples qualified for shotgun metagenomic sequencing.
The publicly available ‘monkeypox-nf’ workflow developed by the Public Health Agency of Canada's National Microbiology Laboratory (11) was adapted for generating the MPXV consensus sequence. 10x and 5x depth thresholds and the MT903343.1 sequence from the B.1 hMPXV (human MPXV) lineage was used as the reference sequence for assembly. The script weeSAM was used to generate coverage depth plots (12). Nextclade was used for clade and lineage assignment, identification of single nucleotide variant (SNV) mutations, insertions, deletions, and for the phylogenetic placement of sequences on a reference tree. Tablet was used to inspect aligned reads supporting the identified SNV mutations.
Table 3 and Figure S1 of the supplementary file show the number of reads, sequencing depth, and genome coverage generated by the 3 MPXV positive samples. The total number of reads for each specimen ranges from 0.8M to ~1.2M reads wherein around 2.8k to 7.7k were mapped to the MPXV reference sequence (MT903343.1). Since all three samples were collected from the same patient, sequences generated by the three specimens were pooled to increase the number of reads and genome coverage. A total of ~5M reads with 14.5k reads mapped to the MPXV reference sequence and only 15% genome coverage using 10x default Illumina depth threshold. The genome coverage was increased to 80% when the depth threshold was lowered to 5x however this increases the likelihood of misclassifying mutations.
Table 3. Sequencing reads and MPXV genome coverage generated from the performed shotgun metagenomic sequencing
Sample
|
No. of paired-end reads from sample
|
No. of reads mapped to MPXV reference
|
Mean sequencing depth
|
Depth threshold (default for illumina=10x)
|
# Ns in consensus
|
Coverage over MPXV genome (0.0-1.0)
|
MPOX22-0034CH
|
1,210,779
|
3987
|
2.82958x
|
10x
|
197,134
|
0
|
5x
|
189,399
|
0.0392373
|
MPOX22-0034DSA
|
829,239
|
2824
|
2.10975x
|
10x
|
197,134
|
0
|
5x
|
189,831
|
0.037041
|
MPOX22-0034DSE
|
1,247,310
|
7721
|
5.1794x
|
10x
|
196,808
|
0.0016537
|
5x
|
158,724
|
0.194801
|
Pooled MPOX22-0034
|
5,402,912
|
14,532
|
10.1189x
|
10x
|
167,425
|
0.150705
|
5x
|
38,363
|
0.805382
|
Assembly method: reference-based assembly with monkeypox-nf workflow (bwa, samtools, ivar)
The Nextclade analysis of consensus sequences generated from all three separate MPOX22-0034 samples, the pooled MPOX22-0034 sample, and using 10× and 5× depth thresholds with respect to a reference sequence from the hMPXV outbreak clade (i.e., MPXV_USA_2022_MA001 in NC_063383 coordinates or pseudo_ON563414) is shown in Figure S2 of the supplementary file. Only up to three single-nucleotide variant (SNV) mutations were identified when comparing the consensus sequences to this reference sequence, which was from a sample collected in May 2022.
Figure 5 shows the phylogenetic placement by Nextclade of the MPXV consensus sequences on a reference tree representing the different lineages under the hMPXV clade. The consensus sequences are placed in the B.1/B.1.3 lineage, showing that the first detected monkeypox case in the Philippines belongs to the B.1 or more specifically the B.1.3 lineage. Sequences with lower % coverage (MPOX22-0034DSA with 5× depth threshold and pooled_MPOX22 with 10× depth threshold) are placed at the base of the B.1.3 lineage while the sequence with the highest % coverage (pooled_MPOX22 with 5× depth threshold) is placed in a subtree within the B.1.3 lineage wherein the consensus sequence clusters with sequences from multiple European countries (including France, Germany, Finland, Switzerland, Spain, Belgium, and Slovenia) and the United States.
The three detected SNV mutations from the pooled_MPOX22-0034 consensus sequence are supported by 100% of reads covering their respective genome positions, which are the coordinates 55133, 64426, and 190660 in the pseudo_ON563414 reference sequence (for more details, see the supplementary information, Figure S3). False positive mutations among these three SNVs are thus unlikely. Two of these mutations (G55133A/OPG074:R665C and C64426T/no amino acid change) are unique to B.1.3 (Figure S3A and B), supporting the validity of the lineage assignment of the consensus sequence. C64426T in particular is unique to the cluster of B.1.3 sequences originating from multiple European countries. One of the three mutations (G190660A/NBT03_gp174:R84K) appears in both B.1.3 and in B.1 sequences. Because of the positive varicella IgM result, metagenomic sequence data were reviewed to verify if there is a co-infection. No sequence data indicating varicella was present in the sample which supports the finding that there is no co-infection.