Hemophagocytic lymphohistiocytosis (HLH) is a severe, rapidly progressive disease. While viral infection is considered a common etiology of pediatric HLH, specific causative viruses other than the Epstein-Barr virus (EBV) have been rarely identified. This study utilized metagenomic next-generation sequencing (NGS) to identify potential causative pathogens in plasma samples from 17 pediatric patients with suspected HLH. Additionally, one case each of confirmed EBV- and cytomegalovirus (CMV)-associated HLH was analyzed for methodological validation. Plasma cell-free RNA (cfRNA) profiling was performed using NGS data to assess the host transcriptome response. Significant viral reads of human herpesvirus-6B, human herpesvirus-7, and Hubei reo-like virus (HRLV) 14 were detected using metagenomic NGS in one patient each. Plasma cfRNA profiles from five patients with viral infection (including EBV and CMV) were compared to those of 14 patients without viral infection. By comparing the two patient groups, 1053 differentially expressed genes were identified. The gene ontology (GO) term of "adaptive immune response" (GO: 0002250) was significantly enriched among upregulated genes in the virus-positive group. Furthermore, an isolated cluster consisting specifically of mitochondrial RNAs, was identified in the upregulated genes of the virus-positive group. Using metagenomic NGS, several candidate viral pathogens were identified in patients with suspected infection-related HLH. The viral genome of HRLV 14, previously undetected in human clinical samples, was identified in one patient. The results from plasma cfRNA profiling suggest that mitochondrial RNAs may reflect the underlying pathogenesis of virus-associated HLH and have potential utility as disease biomarkers.