Integrated comparative transcriptome and weighted gene co-expression network analysis provide valuable insights into the response mechanisms of crayfish (Procambarus clarkii) to copper stress.

One of the significant limitations of aquaculture worldwide is the prevalence of divalent copper (Cu). Crayfish (Procambarus clarkii) are economically important freshwater species adapted to a variety of environmental stimuli, including heavy metal stresses; however, large-scale transcriptomic data of the hepatopancreas of crayfish in response to Cu stress are still scarce. Here, integrated comparative transcriptome and weighted gene co-expression network analyses were initially applied to investigate gene expression profiles of the hepatopancreas of crayfish subjected to Cu stress for different periods. As a result, 4662 significant differentially expressed genes (DEGs) were identified following Cu stress. Bioinformatics analyses revealed that the "focal adhesion" pathway was one of the most significantly upregulated response pathways following Cu stress, and seven DEGs mapped to this pathway were identified as hub genes. Furthermore, the seven hub genes were examined by quantitative PCR, and each was found to have a substantial increase in transcript abundance, suggesting a critical role of the "focal adhesion" pathway in the response of crayfish to Cu stress. Our transcriptomic data can be a good resource for the functional transcriptomics of crayfish, and these results may provide valuable insights into the molecular response mechanisms underlying crayfish to Cu stress.