Integrated dual transcriptome sequencing and experimental validation reveal potential mechanisms of baicalin against pneumocystis pneumonia in immunosuppressed rats

Abstract

Background Pneumocystis pneumonia (PCP) remains a major cause of morbidity and mortality in immunocompromised individuals. Although baicalin (Ba), a natural bioactive flavonoid, has demonstrated protective and therapeutic effects against PCP, its molecular mechanisms remain undefined. We employed dual RNA sequencing (dual RNA-seq) to characterize host and pathogen transcriptional responses to Ba treatment in an immunosuppressed rat model of PCP. Methods Comparative transcriptomic analyses identified differentially expressed genes in both the host and Pneumocystis, followed by Gene Ontology, Kyoto Encyclopedia of Genes and Genomes, and gene set enrichment analyses. Candidate targets were further investigated using network pharmacology, protein–protein interaction analysis, molecular docking, and molecular dynamics simulations. Key findings were validated by immunohistochemistry, enzyme-linked immunosorbent assay, and quantitative PCR. Results Ba markedly remodeled host and pathogen transcriptomes. Host transcriptomic analyses showed that Ba attenuated inflammatory and oxidative stress responses by modulating immune-related pathways, including Toll-like receptor, NF-κB, cytokine-cytokine receptor interaction, chemokine signaling, Th17 cell differentiation, and antigen processing and presentation. Experimental validation demonstrated that Ba reduced pulmonary expression of indoleamine 2,3-dioxygenase 1 (IDO1), Toll-like receptor 2 (TLR2), and TLR4 while increasing nuclear factor erythroid 2-related factor 2 (Nrf2) and its downstream antioxidant enzyme heme oxygenase-1 (HO-1). Pathogen transcriptomic analysis identified Pneumocystis Rtt109 (PcRtt109), a fungal histone acetyltransferase, as a potential pathogen-specific target that was significantly downregulated after Ba treatment. Molecular docking and molecular dynamics simulations supported stable interactions between Ba and IDO1, Nrf2, TLR2, TLR4, and PcRtt109, with the strongest predicted binding observed for PcRtt109. Conclusion Dual RNA-seq revealed that Ba exerts anti-PCP activity through coordinated modulation of host and pathogen molecular networks. Its therapeutic effects are associated with suppression of inflammatory signaling, enhancement of antioxidant defenses, and inhibition of a fungal virulence-associated target. These findings provide mechanistic insights into host-pathogen interactions during PCP and support Ba as a potential therapeutic candidate for PCP.