Comparative transcriptomic analysis provides insights into antibacterial mechanisms of Branchiostoma belcheri under Vibrio parahaemolyticus infection

Evaluation of Tunisian wheat endophytes as plant growth promoting bacteria and biological control agents against Fusarium culmorum

Plant growth-promoting rhizobacteria (PGPR) applications have emerged as an ideal substitute for synthetic chemicals by their ability to improve plant nutrition and resistance against pathogens. In this study, we isolated fourteen root endophytes from healthy wheat roots cultivated in Tunisia. The isolates were identified based from their 16S rRNA gene sequences. They belonged to Bacillota and Pseudomonadota taxa. Fourteen strains were tested for their growth-promoting and defense-eliciting potentials on durum wheat under greenhouse conditions, and for their in vitro biocontrol power against Fusarium culmorum, an ascomycete responsible for seedling blight, foot and root rot, and head blight diseases of wheat. We found that all the strains improved shoot and/or root biomass accumulation, with Bacillus mojavensis, Paenibacillus peoriae and Variovorax paradoxus showing the strongest promoting effects. These physiological effects were correlated with the plant growth-promoting traits of the bacterial endophytes, which produced indole-related compounds, ammonia, and hydrogen cyanide (HCN), and solubilized phosphate and zinc. Likewise, plant defense accumulations were modulated lastingly and systematically in roots and leaves by all the strains. Testing in vitro antagonism against F. culmorum revealed an inhibition activity exceeding 40% for five strains: Bacillus cereus, Paenibacillus peoriae, Paenibacillus polymyxa, Pantoae agglomerans, and Pseudomonas aeruginosa. These strains exhibited significant inhibitory effects on F. culmorum mycelia growth, sporulation, and/or macroconidia germination. P. peoriae performed best, with total inhibition of sporulation and macroconidia germination. These finding highlight the effectiveness of root bacterial endophytes in promoting plant growth and resistance, and in controlling phytopathogens such as F. culmorum. This is the first report identifying 14 bacterial candidates as potential agents for the control of F. culmorum, of which Paenibacillus peoriae and/or its intracellular metabolites have potential for development as biopesticides.

Integrative analyses of mRNA and microRNA expression profiles reveal the innate immune mechanism for the resistance to Vibrio parahaemolyticus infection in Epinephelus coioides

Vibrio parahaemolyticus, as one of the main pathogens of marine vibriosis, has brought huge losses to aquaculture. However, the interaction mechanism between V. parahaemolyticus and Epinephelus coioides remains unclear. Moreover, there is a lack of comprehensive multi-omics analysis of the immune response of grouper spleen to V. parahaemolyticus. Herein, E. coioides was artificially injected with V. parahaemolyticus, and it was found that the mortality was 16.7% in the early stage of infection, and accompanied by obvious histopathological lesions in the spleen. Furthermore, 1586 differentially expressed genes were screened by mRNA-seq. KEGG analysis showed that genes were significantly enriched in immune-related pathways, Acute-phase immune response, Apoptosis, Complement system and Cytokine-cytokine receptor interaction. As for miRNA-seq analysis, a total of 55 significantly different miRNAs were identified. Further functional annotation analysis indicated that the target genes of differentially expressed miRNAs were enriched in three important pathways (Phosphatidylinositol signaling system, Lysosome and Focal adhesions). Through mRNA-miRNA integrated analysis, 1427 significant miRNA–mRNA pairs were obtained and “p53 signaling pathway”, “Intestinal immune network for IgA production” were considered as two crucial pathways. Finally, miR-144-y, miR-497-x, novel-m0459-5p, miR-7133-y, miR-378-y, novel-m0440-5p and novel-m0084-3p may be as key miRNAs to regulate immune signaling pathways via the miRNA-mRNA interaction network. The above results suggest that the mRNA-miRNA integrated analysis not only sheds new light on the molecular mechanisms underlying the interaction between host and V. parahaemolyticus but also provides valuable and new insights into resistance to vibrio infection.

Biocontrol and Action Mechanism of Bacillus subtilis Lipopeptides' Fengycins Against Alternaria solani in Potato as Assessed by a Transcriptome Analysis

Alternaria solani is an airborne fungus and the primary causal agent of potato early blight worldwide. No available fungicides that are both effective and environmentally friendly are usable to control this fungus. Therefore, biological control is a potential approach for its suppression. In this study, Bacillus subtilis strain ZD01's fermentation broth strongly reduced A. solani pathogenicity under greenhouse conditions. The effects of strain ZD01's secondary metabolites on A. solani were investigated. The exposure of A. solani hyphae to the supernatant resulted in swelling and swollen sacs, and the ZD01 supernatant reduced A. solani conidial germination significantly. Matrix-assisted laser desorption/ionization time of flight mass spectrometry and pure product tests revealed that fengycins were the main antifungal lipopeptide substances. To elucidate the molecular mechanism of the fengycins' biological control, RNA sequencing analyses were performed. A transcriptome analysis revealed that 304 and 522 genes in A. solani were differentially expressed after 2-h and 6-h fengycin treatments, respectively. These genes were respectively mapped to 53 and 57 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. In addition, the most enriched KEGG pathway analysis indicated that the inhibitory mechanisms of fengycins against A. solani regulated the expression of genes related to cell wall, cell membrane, transport, energy process, protein synthesis and genetic information. In particular, cell wall and cell membrane metabolism were the main processes affected by fengycin stress. Scanning and transmission electron microscope results revealed hyphal enlargement and a wide range of abnormalities in A. solani cells after exposure to fengycins. Furthermore, fengycins induced chitin synthesis in treated cells, and also caused the capture of cellular fluorescent green labeling and the release of adenosine triphosphate (ATP) from outer membranes of A. solani cells, which may enhance the fengycins ability to alter cell membrane permeability. Thus, this study increases the transcriptome data resources available and supplies a molecular framework for B. subtilis ZD01 inhibition of A. solani HWC-168 through various mechanisms, especially damaging A. solani cell walls and membranes. The transcriptomic insights may lead to an effective control strategy for potato early blight.

Gene expression kinetics of Exaiptasia pallida innate immune response to Vibrio parahaemolyticus infection

BACKGROUND

Recent sequencing projects on early-diverging metazoans such as cnidarians, have unveiled a rich innate immunity gene repertoire; however, little is known about immunity gene regulation in the host's early response against marine bacterial pathogens over time. Here, we used RNA-seq on the sea anemone Exaiptasia pallida (Ep) strain CC7 as a model to depict the innate immune response during the onset of infection with the marine pathogenic bacteria Vibrio parahaemolyticus (Vp) clinical strain O3:K6, and lipopolysaccharides (LPS) exposure. Pairwise and time series analyses identified the genes responsive to infection as well as the kinetics of innate immune genes over time. Comparisons between the responses to live Vp and purified LPS was then performed.

RESULTS

Gene expression and functional analyses detected hundreds to thousands of genes responsive to the Vp infection after 1, 3, 6 and 12 h, including a few shared with the response to LPS. Our results bring to light the first indications that non-canonical cytoplasmic pattern recognition receptors (PRRs) such as NOD-like and RIG-I-like receptor homologs take part in the immune response of Ep. Over-expression of several members of the lectin-complement pathways in parallel with novel transmembrane and Ig containing ficolins (CniFLs) suggest an active defense against the pathogen. Although lacking typical Toll-like receptors (TLRs), Ep activates a TLR-like pathway including the up-regulation of MyD88, TRAF6, NF-κB and AP-1 genes, which are not induced under LPS treatment and therefore suggest an alternative ligand-to-PRR trigger. Two cytokine-dependent pathways involving Tumor necrosis factor receptors (TNFRs) and several other potential downstream signaling genes likely lead to inflammation and/or apoptosis. Finally, both the extrinsic and intrinsic apoptotic pathways were strongly supported by over-expression of effector and executioner genes.

CONCLUSIONS

To our knowledge, this pioneering study is first to follow the kinetics of the innate immune response in a cnidarian during the onset of infection with a bacterial pathogen. Overall, our findings reveal the involvement of both novel immune gene candidates such as NLRs, RLRs and CniFLs, and previously identified TLR-like and apoptotic pathways in anthozoan innate immunity with a large amount of transcript-level evidence.

Evaluation of IL-17D in Host Immunity to Group A Streptococcus Infection

IL-17D is a cytokine that belongs to the IL-17 family and is conserved in vertebrates and invertebrates. In contrast to IL-17A and IL-17F, which are expressed in Th17 cells, IL-17D is expressed broadly in nonimmune cells. IL-17D can promote immune responses to cancer and viruses in part by inducing chemokines and recruiting innate immune cells such as NK cells. Although bacterial infection can induce IL-17D in fish and invertebrates, the role of mammalian IL-17D in antibacterial immunity has not been established. To determine whether IL-17D has a role in mediating host defense against bacterial infections, we studied i.p. infection by group A Streptococcus (GAS) in wild-type (WT) and Il17d ^-/- mice. Compared with WT animals, mice deficient in IL-17D experienced decreased survival, had greater weight loss, and showed increased bacterial burden in the kidney and peritoneal cavity following GAS challenge. In WT animals, IL-17D transcript was induced by GAS infection and correlated to increased levels of chemokine CCL2 and greater neutrophil recruitment. Of note, GAS-mediated IL-17D induction in nonimmune cells required live bacteria, suggesting that processes beyond recognition of pathogen-associated molecular patterns were required for IL-17D induction. Based on our results, we propose a model in which nonimmune cells can discriminate between nonviable and viable GAS cells, responding only to the latter by inducing IL-17D.

Transcriptome Analysis of Bursaphelenchus xylophilus Uncovers the Impact of Stenotrophomonas maltophilia on Nematode and Pine Wilt Disease

Stenotrophomonas maltophilia influences the reproduction, pathogenicity, and gene expression of aseptic Bursaphelenchus xylophilus after inoculation of aseptic Pinus massoniana. Pine wilt disease is a destructive pine forest disease caused by B. xylophilus, and its pathogenesis is unclear. The role of bacteria associated with B. xylophilus in pine wilt disease has attracted widespread attention. S. maltophilia is one of the most dominant bacteria in B. xylophilus, and its effect is ambiguous. This study aims to explore the role of S. maltophilia in pine wilt disease. The reproduction and virulence of aseptic B. xylophilus and B. xylophilus containing S. maltophilia were examined by inoculating aseptic P. massoniana seedlings. The gene expressions of two nematode treatments were identified by transcriptome sequencing. The reproduction and virulence of B. xylophilus containing S. maltophilia were stronger than that of aseptic nematodes. There were 4240 differentially expressed genes between aseptic B. xylophilus and B. xylophilus containing S. maltophilia after inoculation of aseptic P. massoniana, including 1147 upregulated genes and 2763 downregulated genes. These differentially expressed genes were significantly enriched in some immune-related gene ontology (GO) categories, such as membrane, transporter activity, metabolic processes, and many immune-related pathways, such as the wnt, rap1, PI3K-Akt, cAMP, cGMP-PKG, MAPK, ECM-receptor interaction, and calcium signaling pathways. The polyubiquitin-rich gene, leucine-rich repeat serine/threonine-protein kinase gene, glyceraldehyde-3-phosphate dehydrogenase (GAPDH) gene, acetyl-CoA carboxylase gene, and heat shock protein genes were the key genes associated with immune resistance. Moreover, there were four cell wall hydrolase genes, thirty-six detoxification- and pathogenesis-related protein genes, one effector gene and ten cathepsin L-like cysteine proteinase genes that were differentially expressed. After inoculation of the host pine, S. maltophilia could affect the virulence and reproduction of B. xylophilus by regulating the expression of parasitic, immune, and pathogenicity genes of B. xylophilus.

Bacillus subtilis inhibits Aspergillus carbonarius by producing iturin A, which disturbs the transport, energy metabolism, and osmotic pressure of fungal cells as revealed by transcriptomics analysis

The contamination of Aspergillus carbonarius causes decreases and great decay of agricultural products, and threatens the human and animal health by producing mycotoxins, especially ochratoxin A. Bacillus subtilis has been proved to efficiently inhibit the growth of A. carbonarius. Revealing the major active compound and the mechanisms for the antifungal of B. subtilis are essential to enhance its antifungal activity and control the quality of antifungal products made of it. In this study, we determined that iturin A is the major compound that inhibits Aspergillus carbonarius, a widespread fungal pathogen of grape and other fruits. Iturin A significantly inhibited growth and ochratoxin A production of A. carbonarius with minimal inhibitory concentrations (MICs) of 10 μg/mL and 0.312 μg/mL, respectively. Morphological observations revealed that iturin A caused swelling of the fungal cells and thinning of the cell wall and membrane at 1/2 MIC, whereas it inhibited fungal spore germination and caused mitochondrial swelling at higher concentrations. A differential transcriptomic analysis indicated that the mechanisms used by iturin A to inhibit A. carbonarius were to downregulate the expression of genes related to cell membrane, transport, osmotic pressure, oxidation-reduction processes, and energy metabolism. Among the down-regulated genes, those related to the transport capacity were most significantly influenced, including the increase of energy-related transport pathways and decrease of other pathways. Notably, the genes related to taurine and hypotaurine metabolism were also decreased, indicating iturin A potentially cause the occurrence of osmotic imbalance in A. carbonarius, which may be the intrinsic cause for the swelling of fungal cells and mitochondria. Overall, iturin A produced by B. subtilis played important roles to inhibit A. carbonarius via changing the fungal cell structure and causing perturbations to energy, transport and osmotic pressure metabolisms in fungi. The results indicated a new direction for researches on the mechanisms for lipopeptides and provided useful information to develop more efficient antifungal agents, which are important to agriculture and biomedicine.

Functional Variation of IL-1R-Associated Kinases in the Conserved MyD88-TRAF6 Pathway during Evolution

IL-1R-associated kinases (IRAK) are important regulators in the TLR/IL-1R pathways, but their function appears inconsistent between Drosophila, bony fishes, and vertebrates. This causes a difficulty to understand the IRAK functions. As a step to reveal the evolution of IRAKs, in this study, we performed comparative and functional analysis of IRAKs by exploiting the amphioxus, a pivotal taxon connecting invertebrates and vertebrates. Sequence and phylogenetic analysis indicated three major IRAK lineages: IRAK1/2/3 is a vertebrate-specific lineage, IRAK4 is an ancient lineage conserved between invertebrate and vertebrates, and Pelle is another ancient lineage that is preserved in protostomes and invertebrate deuterostomes but lost in vertebrate deuterostomes. Pelle is closer neither to IRAK4 nor to IRAK1/2/3, hence suggesting no clear functional analogs to IRAK1/2/3 in nonvertebrates. Functional analysis showed that both amphioxus IRAK4 and Pelle could suppress NF-κB activation induced by MyD88 and TRAF6, which are unlike mammalian and Drosophila IRAKs, but, surprisingly, similar to bony fish IRAK4. Also unlike Drosophila IRAKs, no interaction was detected between amphioxus IRAK4 and Pelle, although both of them were shown capable of binding MyD88. These findings, together with previous reports, show that unlike other signal transducers in the TLR/IL-1R pathways, such as MyD88 and TRAF6, the functions of IRAKs are highly variable during evolution and very specialized in different major animal taxa. Indeed, we suggest that the functional variability of IRAKs might confer plasticity to the signal transduction of the TLR/IL-1R pathways, which in return helps the species to evolve against the pathogens.

Comparative transcriptomic analysis provides insights into immune responses of lamprey larvae under three pathogens infections

Based on the large number of deaths of lamprey larvae infected by pathogens in the culture process, the infection response of lamprey larvae to different pathogens was studied. Thus, high-throughput sequencing is firstly used to analyze the lamprey larvae infected by Aeromonas hydrophila, Vibrosplendidus, and Saprolegnia parasitica. A high quality of 338,628,426 Reads were obtained and 203,829 transcripts were assembled. A total of 158,001 Unigenes were annotated by Blast comparison, and a relatively comprehensive transcriptome database of lamprey larvae was established. Based on the transcriptome of lamprey larvae infected by three aquatic pathogens, classification results of COG and GO show that the genes classified into signal transduction mechanism accounted for the largest proportion, and 2316 differentially expressed genes were identified and screened out. After infection by the three pathogens, 16 genes and 19 genes expression were up-regulated and down-regulated, respectively. KEGG pathway analysis of differentially expressed genes showed that lamprey larvae were mainly affected their metabolism and oxidative phosphorylation pathways after infection by two gram-negative bacteria, Vibrosplendidus and Aeromonas hydrophila. As a fungus, Saprolegnia parasitica mainly affected ribosome synthesis in lamprey larvae. The results were validated by detecting the expression levels of 13 DEGs at 24 h after pathogens treatment using Q-PCR. The results provide valuable information for further research into the development of immunity and the innate immune response against pathogen invasion during the artificial propagation of lamprey larvae.

Transcriptomic profiling of cotton Gossypium hirsutum challenged with low-temperature gradients stress

Gossypium hirsutum, a cotton species widely cultivated around the world, is a typical cold-sensitive crop. Low-temperature (LT) stress is one of the main environmental stressors that can affect growth and the quality of cotton fibers. LT is also a major challenge for cotton survival, growth maturity and geographical distribution. However, few genome-wide transcriptional response and profiling datasets are available to explore the LT-tolerant mechanism of cotton. This study treated G. hirsutum with four LT gradients (control at 25 °C and cold temperatures at 4 °C, 10 °C and 15 °C) for 24 hour to generate 12 RNA-Seq datasets (three biological replicates per treatment) with approximately 280 million clean reads per dataset. The quality of the datasets obtained in the current study was validated through a series of quality checks including verification of RNA sample quality and RNA-Seq read quality. Data analyses included novel gene discovery, global gene expression profiling and quantitative real-time PCR. This is the first study to report genome-wide transcriptomic datasets for cotton in response to LT exposure.

Genome-wide gene expression analysis reveals novel insights into the response to nitrite stress in gills of Branchiostoma belcheri

Amphioxus has been widely used as a model for the comparative immunology of vertebrates. Studies have reported that gene expression changes in the amphioxus gill in response to biotic stress, such as microbial and their mimic challenge, but little is known about how gene expression is affected by abiotic stress in the marine environment, such as nitrite. A lack of information regarding gene expression response to abiotic stress hinders a comprehensive understanding of gill defense response in amphioxus. Here, RNA sequencing was used to carry out gene expression profiling analyses of Branchiostoma belcheri gills under nitrite stress. Six libraries were created for the control and treatment groups, including three biological replicates. In total, 2416 differently expressed genes (DEGs) were detected in response to nitrite stress, of which 1522 DEGs were up-regulated in the treatment group in comparison to the control, while the remaining 894 DEGs were down-regulated genes. Functional enrichment revealed that these DEGs are primarily involved in disease, innate immunity, xenobiotic biodegradation and metabolism, and biomolecular processes and apoptosis. We screened 11 key nitrite-responsive DEGs to detect their expression responses to nitrite stress at different time points, and validate the sequencing data using real time quantitative PCR. The results indicated that the expression of gene encoding CYP3A, POD, CASPR1, GST, MAO, DDH, and XDH/XO were induced, while those encoding MRC, GT, DNASE1L, and RIPK5 were reduced, to participate in the anti-nitrite response. This study provides a useful resource for research of molecular toxicology in amphioxus under environmental stress.

Identification of circular RNAs and their altered expression under poly(I:C) challenge in key antiviral immune pathways in amphioxus

Amphioxus is a key model for studying comparative immunity of vertebrates. Circular RNA (circRNA), as RNAs with a circular structure, has received little attention until recently, where several studies have reported that circRNA expression changes are involved in the immune response in animals. However, circRNA and its immune role in amphioxus have not been previously studied. Here, circRNAs in Chinese amphioxus (Branchiostoma belcheri) were sequenced, and 1859 circRNAs were identified using two algorithms (find_circ and CIRI). The analysis of miRNA target sites on circRNAs showed that 332 circRNAs may function as miRNA sponges. Furthermore, we identified circRNAs that were conserved between B. belcheri and vertebrates, tracing the origin of these circRNAs within chordates. Additionally, in combination with several key antiviral immune (poly(I:C), pIC) pathways identified in our previous B. belcheri studies, nine circRNAs potentially involved in these pathways were identified using bioinformatic predictions. Among these nine circRNAs, eight were selected to examine their expression response in B. belcheri challenged by pIC in comparison to control using real-time quantitative PCR. The results showed that four circRNAs were induced as part of the antiviral response against pIC, while expression of two circRNAs was decreased, and the expression levels of the remaining two were not significantly altered after pIC challenge. This work is the first to identify circRNAs and reveal their antiviral role in amphioxus. Therefore, it opens a new window to explore the comparative immunology of circRNAs in chordates and the regulatory roles of circRNAs in antiviral immunity in amphioxus.

Genome-wide transcriptional response of microRNAs to the benzo(a)pyrene stress in amphioxus Branchiostoma belcheri

Amphioxus, a cephalochordate found in sand habitats in shallow in-shore seawaters, has been widely used as a model in comparative immunology of chordates. However, the role of microRNAs (miRNAs) in amphioxus under abiotic stress, particularly xenobiotics with strong toxicity, remains largely unknown. Here, a widespread marine contaminant, benzo(a)pyrene (BaP) is used to evaluate its toxic effects on miRNA expression of amphioxus. Six small RNA libraries were sequenced from Branchiostoma belcheri. A total of 144 known and 157 novel miRNAs were identified using deep sequencing and bioinformatics approaches. A total of 58 differentially expressed miRNAs (DEMs) were screened, including 25 up- and 33 down-regulated DEMs under BaP stress. Target genes possibly regulated by DEMs were predicted, and their functional enrichment analyses were performed. Targets of DEMs are primarily involved in xenobiotic and cellular homeostasis, catabolic and transport process. They could be largely linked to nine immune- and toxin detoxification-related pathways, including metabolism of xenobiotics by cytochrome P450, drug metabolism-other enzymes, and drug metabolism-cytochrome P450, etc. Furthermore, quantitative real-time PCR (qRT-PCR) analysis for 12 key BaP-responsive DEMs validates the accuracy of deep sequencing. Experiments were then conducted to investigate their expression responses to BaP stress at different time intervals in detail to further determine their expression dynamic in responses of B. belcheri towards BaP exposure. This study, to the best of our knowledge, investigates the regulatory roles of miRNAs in the toxicological response of amphioxus for the first time, providing valuable information for the protection of lone existing cephalochordate amphioxus.

Liver transcriptome analysis of the Sparus macrocephlus in response to Vibrio parahaemolyticus infection

The black seabream (Sparus macrocephlus) is an economically pivotal aquaculture species cultured in China and Southeast Asian countries. To understand the molecular immune mechanisms underlying the response to Vibrio parahaemolyticus, a comparative gene transcription analysis were performed with utilized fresh livers of V. parahaemolyticus-immunized Sparus macrocephlus with a control group through RNA-Seq technology. A total of 256663 contigs were obtained after excluded the low-quality sequences and assembly. The average length of contigs collected from this research is 1066.93 bp. Furthermore, blast analysis indicates 30747 contigs were annotated based on homology with matches in the NT, NR, gene, and string databases. A gene ontology analysis was employed to classify 21598 genes according to three major functional categories: molecular function, cellular component, and biological process. A total of 14470 genes were discovered in 303 KEGG pathways. RSEM and EdgeR were introduced to estimate 3841 genes significantly different expressed (False Discovery Rate<0.001) which includes 4072 up-regulated genes and 3771 down-regulated genes. A significant enrichment analysis of these differentially expressed genes and isogenes were conducted to reveal the major immune-related pathways which refer to the toll-like receptor, complement, coagulation cascades, and chemokine signaling pathways. In addition, 92175 potential simple sequence repeats (SSRs) and 121912 candidate single nucleotide polymorphisms (SNPs) were detected and identified sequencely in the Sparus macrocephlus liver transcriptome. This research characterized a gene expression pattern for normal and the V. parahaemolyticus -immunized Sparus macrocephlus for the first time and not only sheds new light on the molecular mechanisms underlying the host-V. parahaemolyticus interaction but contribute to facilitate future studies on Sparus macrocephlus gene expression and functional genomics.

Genome-wide identification of the interactions between key genes and pathways provide new insights into the toxicity of bisphenol F and S during early development in zebrafish

Bisphenol F (BPF) and bisphenol S (BPS) have been widely used as alternatives to bisphenol A (BPA). With their increasing use, BPF and BPS have also been released into the environment; thus, their potential risks to aquatic organisms and humans are drawing attention. The objective of this study was to identify the interactions between key pathways and hub genes in zebrafish following BPF and BPS exposure, and to evaluate the potential risks to human health. We identified three key pathways using KEGG over-representation test and Gene Set Enrichment Analysis (GSEA): 'Necroptosis,' 'Adipocytokine signaling pathway,' and 'C-type lectin receptor signaling pathway.' Moreover, three hub genes (mst1ra, prkcdb, and pik3cb) and detailed interactions among the pathways were examined by the analyses of PPI network, subcellular location, and shortest-pathway. Surprisingly, all three pathways were strongly associated with a potential risk of cancer, as reported previously. In addition, the results of KOBAS shown in 'Pathways in Cancer' and 'Cancers' belong to the top 10 terms in pathway enrichment analyses using genes related to BPF or BPS in human, as was found using GenCLiP. Moreover, the Kaplan-Meier survival analysis was performed using homologenes (MST1R, PIK3CB and PRKCD) of hub genes in human to evaluate whether exposure to bisphenols may adversely affect breast cancer. Taken together, these studies demonstrate the potential carcinogenicity of BPF and BPS. To our knowledge, this is the first study on three overlapping key pathways and three hub genes to investigate BPF and BPS exposure-related mechanisms and subsequent interactions in zebrafish.