Identification of Crucial Modules and Genes Associated with Bt Gene Expression in Cotton

The expression of Bacillus thuringiensis (Bt) toxins in transgenic cotton confers resistance to insect pests. However, it has been demonstrated that its effectiveness varies among cotton cultivars and different tissues. In this study, we evaluated the expression of Bt protein in 28 cotton cultivars and selected 7 cultivars that differed in Bt protein expression for transcriptome analysis. Based on their Bt protein expression levels, the selected cultivars were categorized into three groups: H (high Bt protein expression), M (moderate expression), and L (low expression). In total, 342, 318, and 965 differentially expressed genes were detected in the H vs. L, M vs. L, and H vs. M comparison groups, respectively. And three modules significantly associated with Bt protein expression were identified by weighted gene co-expression network analysis. Three hub genes were selected to verify their relationships with Bt protein expression using virus-induced gene silencing (VIGS). Silencing GhM_D11G1176, encoding an MYC transcription factor, was confirmed to significantly decrease the expression of Bt protein. The present findings contribute to an improved understanding of the mechanisms that influence Bt protein expression in transgenic cotton.