To establish a whole-genome sequencing (WGS) workflow and molecular tracing method for monkeypox virus (MpxV) based on amplicon sequencing technology combined with a high-throughput sequencing platform and bioinformatics analysis techniques, providing technical support for molecular epidemiological research of monkeypox cases in Dalian city.
By taking total nucleic acid extracted from vesicle fluid samples of monkeypox cases as templates, a cDNA library was constructed on the Ion Torrent sequencing platform via WGS technology. Template preparation and sequencing were completed. Using the built-in plug-ins of Ion GeneStudio S5 Plus and NC063383.1 as the reference sequence, data quality control, base identification, sequence alignment, variant analysis, and whole-genome sequence assembly were performed to obtain consensus sequences (Fasta files). Online analysis tools were used to analyze the whole-genome characteristics, genotyping, and genetic evolutionary relationships of MpxV.
The cycle threshold (Ct) values of MpxV nucleic acid detection in vesicle fluid samples from the two patients were 28.03 and 29.60, respectively. The consensus sequence obtained from sample 1 was a 197 084 bp MpxV whole-genome sequence with a total of 3 470 722 reads, an average sequencing depth of 3 205, an average read length of 182 bp, and approximately 98.8% coverage. The consensus sequence obtained from sample 2 was a 197 099 bp MpxV whole-genome sequence with a total of 1 440 296 reads, an average sequencing depth of 746.7, an average read length of 105 bp, and approximately 97.8% coverage. Sample 1 had 84 single nucleotide mutations and 40 amino acid mutations, while sample 2 had 85 single nucleotide mutations and 42 amino acid mutations. Both samples belonged to the C.1.1 sub-branch of the MpxV IIb clade (West African clade). hMpxV/China/ DLCDC-01 was closely related to the sequences EPI_ISL18213375, EPI_ISL19058867, EPI_ISL18213374, EPI_ISL19058866, and EPI_ISL19058873 submitted from Guangdong province, China, and EPIISL_18059184 submitted from Yunnan province, China. hMpxV/China/DLCDC-02 was not homologous to hMpxV/China/DLCDC-01 or the 20 valid sequences from China uploaded to GISAID.
Consensus sequences of MpxV were obtained from vesicle fluid samples of two patients, successfully establishing WGS technology, bioinformatics analysis, and molecular tracing methods for MpxV in Dalian city.
| 科 Family | 属数 Number of genus | 种数 Number of species | 占总种数比例 Percentage of total species (%) | 属 Genus | 种数 Number of species | 占总种数比例 Percentage of total species (%) |
|---|---|---|---|---|---|---|
| 鹅膏菌科Amanitaceae | 2 | 11 | 5.26 | 鹅膏菌属 Amanita | 10 | 4.78 |
| 小菇科 Mycenaceae | 2 | 12 | 5.74 | 丝盖伞属 Inocybe | 5 | 2.39 |
| 多孔菌科 Polyporaceae | 8 | 14 | 6.70 | 蜡蘑属 Laccaria | 5 | 2.39 |
| 红菇科 Russulaceae | 3 | 23 | 11.00 | 小皮伞属 Marasmius | 6 | 2.87 |
| 小菇属 Mycena | 11 | 5.26 | ||||
| 光柄菇属 Pluteus | 5 | 2.39 | ||||
| 红菇属 Russula | 17 | 8.13 | ||||
| 栓菌属 Trametes | 5 | 2.39 |