Since 2022, global outbreaks of monkeypox virus (MPXV) have been repeatedly designated by the World Health Organization (WHO) as a public health emergency of international concern (PHEIC), underscoring the urgent need to elucidate the multidimensional mechanisms underlying viral evolution and transmission. Current understanding remains largely focused on genomic variation, while the critical role of epigenetic regulation has been considerably overlooked. To address this gap, this study integrates high-throughput evolutionary genomic analysis with whole-genome DNA methylation profiling. Using parallel Illumina and Nanopore sequencing platforms, we comprehensively characterized two clinically derived MPXV isolates collected locally. The results revealed that both isolates belonged to the C.1.1 ancestral lineage, diverging into distinct clades (E.3 and E.4, respectively, supporting the presence of at least two independent viral introduction events into the region, each followed by limited local transmission. They had accrued a considerable number of single-nucleotide polymorphisms (SNPs), with APOBEC3-associated substitutions constituting 84.8% and 77.6% of all observed mutations. Furthermore, both 5-hydroxymethylcytosine (5hmC) and N6-methyladenine (6mA) modifications were identified and found to be preferentially enriched within the inverted terminal repeats (ITRs) regions of MPXV genome in both viral strains; moreover, the E.4 lineage viral strain exhibits a markedly more intricate and compositionally diversified modification landscape, a pattern that indicates appreciable epigenetic heterogeneity among MPXV lineages. Our study furnishes a multi-omics framework that presents a systematic evolutionary feature of two clinical MPXV isolates and their genomic DNA 5hmC and 6mA modification topologies, and enhances our understanding of MPXV viral adaptation and diversification.

Severe fever with thrombocytopenia syndrome (SFTS) is a highly fatal tick-borne disease caused by SFTS virus (SFTSV). The Gn/Gc glycoprotein plays a critical role in viral entry, but whether its amino acid polymorphisms directly determine neurological symptoms remains unclear. We sequenced the M segment from 99 SFTS patients (31 with neurological symptoms, 68 without) in Shandong Province, China. After removing recombinant sequences, 91 non-recombinant sequences were analyzed and genotyped based on M segment phylogeny. Systematic signature screening, phylogenetic reconstruction, genotype-stratified analysis, public database validation, cross-genotype comparison, and selection pressure inference were applied to distinguish functional determinants from genotype-specific markers. M segment-based phylogeny classified the sequences into four genotypes (D, E, F, A), with D accounting for 79.3% of NS cases. Although 10 candidate sites initially distinguished NS from Non-NS, their associations disappeared upon genotype-stratified analysis; all were low-frequency variants (<15%) in public databases. Cross-genotype comparison revealed D-specific (“GERHVGIIFL”) and E/F-specific (“DGKQMRVTSF”) patterns, none of which overlapped with positively selected sites. These findings demonstrate that the apparent NS association reflected genotype confounding rather than functional determinants, identifying the 10 sites as potential genotype D-specific markers with high discriminatory potential, warranting further validation in independent datasets. This study provides a cautionary example for virus–phenotype association studies, emphasizing that genotype stratification is essential before interpreting amino acid differences as functional.