The study examined soybean (Glycine max) responses to simultaneous drought and Asian soybean rust infection using combined transcriptomic and metabolomic analyses. Weighted Gene Co-expression Network Analysis identified stress-specific gene modules linked to metabolites, while Copula Graphical Models uncovered sparse, condition‑specific networks, revealing distinct molecular signatures for each stress without overlapping genes or metabolites. The integrative approach underscores a hierarchical, modular defense architecture and suggests targets for breeding multi‑stress resilient soybeans.