Branching plasticity and candidate gene-hormone networks associated with shade responses in soybean under relay strip intercropping.
BACKGROUND: Branching is a key determinant of high-yield plant architecture in soybean, particularly in maize- soybean relay strip intercropping where plants experience an "initially shaded-then fully illuminated" light regime. However, the genetic regulation of branching responses to shading remains poorly understood. METHODS: We evaluated 11 branching-related traits across 202 soybean accessions grown under monoculture (SS) and relay strip intercropping (RI). Branch number (BN), branching incidence (BI), and total branch length (TBL) were assessed together with stress tolerance indices (STI) and relative distance plasticity index (RDPI). Genome-wide association studies (GWAS) using mixed linear model (MLM) and three-variance-component MLM (3VmrMLM) were combined with haplotype and protein structural analyses to refine candidate genes. RESULTS: Based on Pearson correlation analysis of all 11 traits, BN, BI, and TBL measured before maize harvest showed the strongest and most consistent associations with branch seed weight within the corresponding cropping system (BSW_SS under SS and BSW_RI under RI), whereas other traits showed weaker or environment-dependent associations. Higher STI values calculated from these traits during the co-growth phase were negatively associated with BSW_RI, suggesting weaker compensatory recovery after light restoration in genotypes with more stable early branching patterns between SS and RI. In contrast, mediation analysis indicated that RDPI was positively associated with BSW_RI mainly through improved mature branching architecture (MB_index), which accounted for approximately 70% of the total positive effect. GWAS identified 57 and 74 significant QTNs using MLM and 3VmrMLM, respectively, and LD-window genes were filtered for exonic nonsynonymous or premature stop-codon variants, yielding 883 genes with putative functional variants. Two high-confidence genes emerged: Glyma.02G058600 (PP2C55), exhibiting shading-specific haplotype effects likely linked to GA-mediated branch-stem balance, and Glyma.02G059900 (DA1-related protein), showing stable effects across environments and implicated in ABA-mediated suppression of axillary meristems. CONCLUSIONS: These results provide insight into the genetic and physiological basis of soybean branching responses under relay strip intercropping, clarify that branching plasticity and relative shade tolerance represent distinct response dimensions in this system, and identify putative loci that may be useful for breeding soybean cultivars with improved shade adaptation and yield stability.