Sequencing-based network analysis provides a core set of gene resource for understanding kidney immune response against Edwardsiella tarda infection in Japanese flounder.

Marine organisms are under a frequent threat from various pathogens. Edwardsiella tarda is one of the major fish pathogens infecting both cultured and wild fish species. It can also infect a variety of other vertebrates, including amphibians, reptiles, and mammals, and bacteremia caused by E. tarda can be fatal in humans. The kidney is the largest lymphoid organ in fish, and generating kidney transcriptomic information under different stresses is crucial for understanding molecular mechanisms underlying the immune responses in the kidneys. In this study, we performed transcriptome-wide gene expression profiling of the Japanese flounder (Paralichthys olivaceus) challenged by 8 and 48 h of E. tarda infection. An average of 40 million clean reads per library was obtained, and approximately 81.6% of these reads were successfully mapped to the reference genome. In addition, 1319 and 4439 differentially expressed genes (DEGs) were found at 8 and 48 h post-injection, respectively. Gene Ontology (GO) functional enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were performed to search immune-related DEGs. A protein-protein interaction network was constructed to ascertain the relationship between interacting immune genes during pathogen-induced stress. Based on the KEGG and protein association networks analysis, 24 hub genes were discovered and validated by qRT-PCR. To our knowledge, this study is the first to represent the kidney transcriptome analysis based on protein-protein interaction networks in fish. Our results provide valuable gene resources for further research on kidney immune response in fish, which can significantly improve our understanding of the molecular mechanisms underlying the immune response to E. tarda in humans and other vertebrates.