The “milky disease” in Chinese mitten crabs (Eriocheir sinensis), caused by Metschnikowia bicuspidata, poses significant threats to aquaculture, though its pathogenic mechanisms remain poorly understood. This study employs transcriptomic sequencing to analyze gene expression changes in Metschnikowia bicuspidata under hemocyte challenge, iron overload (1 mmol/mL), and combined stress, with functional validation through Common in Fungal Extracellular Membrane (CFEMgene) overexpression strains. Key findings reveal that (1) hemocyte challenge activated base excision repair (−log₁₀[P] = 7.58) and ribosome biogenesis pathways, indicating fungal adaptation through DNA repair and enhanced protein synthesis to counter host immune attacks (e.g., ROS-mediated damage). (2) Iron overload induced glutathione metabolism and pentose phosphate pathway enrichment, demonstrating mitigation of ferroptosis through NADPH/GSH antioxidant systems and autophagy/proteasome coordination. (3) Under combined stress, ribosome biogenesis (−log₁₀[P] = 1.3) and non-homologous end-joining pathways coordinated DNA repair with stress protein synthesis, complemented by vacuolar V-ATPase-mediated iron compartmentalization. (4) CFEM genes showed significant upregulation under hemocyte stress, with overexpression strains exhibiting enhanced biofilm formation (35% increased MTT cytotoxicity) and infectivity (40% higher infection rate), confirming CFEM domains mediate pathogenesis through iron homeostasis and virulence factor production. This work elucidates how M. bicuspidata employs metabolic reprogramming, oxidative stress responses, and CFEM-mediated iron regulation to establish infection, providing critical insights for developing targeted control strategies against milky disease.

Metschnikowia bicuspidata is the specific pathogen for “milky disease” in the Chinese mitten crab (Eriocheir sinensis), accounting for huge losses to the industry. And yet, there is no precise study describing the pathogenesis of M. bicuspidata, largely hindering the development of novel control methods against its causing diseases. Here, we compared the transcriptomes of M. bicuspidata cells collected from a control group (cultured without E. sinensis hemocytes) and a treatment group (cultured with E. sinensis hemocytes), using RNA sequencing. Through comprehensively analyzing the differentially expressed genes (DEGs), both the most regulated ones and the ones involved in crucial enriched KEGG pathways, we found that certain processes might be required for M. bicuspidata’s survival under hemocyte stress. Key genes involved in oxidative phosphorylation, fatty acid metabolism, upper glycolysis, and gluconeogenesis were upregulated, and those for β-glucan unmasking, autophagy, and cell polarity were downregulated, in the treatment group. Our results suggest that M. bicuspidata colonizes and therefore establishes an infection in E. sinensis via enhancing aerobic respiration, glucose-6-phosphate accumulation, and cell-wall masking. In addition, we applied multiple means to evaluate a series of candidate reference genes and found that PMA1 in combination with ACT1 is the most suitable choice for accurate normalization in quantitative real-time PCR (qRT-PCR) assays. Thus, we used this combination as the reference and performed qRT-PCR verification of several DEGs. It is shown that the expression trends of these tested DEGs in qRT-PCR assays are the same as those in RNA-Seq assays. This study not only provides insights into strategies facilitating M. bicuspidata’s survival within E. sinensis, initially elucidating the pathogenesis of this yeast, but also recommends a useful molecular tool regarding qRT-PCR assays in this pathogen.