Natural Genotypic Variation Underpins Root System Response to Drought Stress in Bambara Groundnut [Vigna subterranea (L.) Verdc.]

Bambara groundnut [Vigna subterranea (L.) Verdc.] is grown in rainfed production systems and suffers from periodic drought stress (DS), leading to yield reductions. Natural genotypic variation for root traits is essential for adaptation to water deficit conditions. However, root traits have not been fully utilised as selection criteria to improve DS in bambara groundnut. The present study explored the natural genotypic variation found in single genotypes of bambara groundnut derived from landraces to identify adaptive differences in tap root length (TRL) and root length density (RLD) in response to DS. A diverse core collection of eight bambara groundnut genotypes from various locations (namely, Gresik, LunT, IITA-686, DodR, S19-3, Tiga nicuru, and Ankpa-4, DipC1), were grown for two seasons (2018 and 2019) in polyvinyl chloride (PVC) columns with well-watered (WW) and 30-day DS treatments. Plant samples were collected at 55 days after emergence (DAE) (30 days of DS) and at 105 DAE (30 days of DS plus 50 days of recovery). Under DS, differential TRL among genotypes at 55 DAE was observed, with DodR recording the longest among genotypes with an increase (1% in 2018) in TRL under DS compared to WW, whereas LunT and IITA-686 showed significant (p < 0.001) decrease in TRL (27 and 25%, respectively, in 2018). Average RLD was observed to have the highest reduction under DS in the 90–110 cm layer (42 and 58%, respectively, in 2018 and 2019). Rainy habitat LunT had limited roots in 2018 and recorded the least (0.06 ± 0.013 cm(–3)) RLD in 2019. However, dry-habitat DodR showed an increase in the RLD (60–90 cm) under DS compared to WW, while dry-habitat S19-3 densely occupied all depths with RLD of 0.16 ± 0.05 and 0.18 ± 0.01 cm cm(–3) in the deepest layer in both seasons, respectively. Reduced RLD under DS showed recovery when the plants were re-watered. These plants were additionally observed to have RLD that surpasses the density in WW at all soil depths at 105 DAE. Also, recovery was shown in Tiga nicuru and DodR (0–30 cm) and IITA-686 (90–110 cm) in 2019. Average RLD under DS treatment was associated with substantial grain yield advantage (R(2) = 0.27 and R(2) = 0.49, respectively) in 2018 and 2019. An increase in TRL allowed DodR to quickly explore water at a deeper soil depth in response to gradually declining soil water availability. High RLD in genotypes such as DodR, DipC1 and S19-3 also offered adaptive advantage over other genotypes under DS. Variation in intrinsic RLD in deeper soil depths in the studied genotypes determines root foraging capacity when facing DS. This suggests that different agroecological environments to which bambara groundnut is subjected in its natural habitat have promoted a phenotypic differentiation in root systems to adapt to ecotypic conditions, which may help offset the impact of DS. The natural genotypic variation exhibited, especially by DodR, could be exploited to identify potential quantitative trait loci (QTLs) that control deep rooting and root length density.