Adaptation at the Extremes of Life: Experimental Evolution with the Extremophile Archaeon Sulfolobus acidocaldarius.

The archaeon Sulfolobus acidocaldarius has emerged as a promising thermophilic model system. Investigating how thermophiles adapt to changing temperatures is a key requirement, not only for understanding fundamental evolutionary processes but also for developing S. acidocaldarius as a chassis for bioengineering. One major obstacle to conducting experimental evolution with thermophiles is the expense of equipment maintenance and energy usage of traditional incubators for high-temperature growth. To address this challenge, a comprehensive experimental protocol for conducting experimental evolution in S. acidocaldarius is presented, utilizing low-cost and energy-efficient bench-top thermomixers. The protocol involves a batch culture technique with relatively small volumes (1.5 mL), enabling tracking of adaptation in multiple independent lineages. This method is easily scalable through the use of additional thermomixers. Such an approach increases the accessibility of S. acidocaldarius as a model system by reducing both initial investment and ongoing costs associated with experimental investigations. Moreover, the technique is transferable to other microbial systems for exploring adaptation to diverse environmental conditions.

Comparative Genomic and Functional Analyses for Insights into Pantoea agglomerans Strains Adaptability in Diverse Ecological Niches.

Pantoea agglomerans inhabit diverse ecological niches, ranging from epiphytes and endophytes in plants, body of animals, and occasionally in the human system. This multifaceted bacterium contributes substantially to plant growth promotion, stress resilience, and biocontrol but can also act as a pathogen to its host. The genetic determinants underlying these diverse functions remain largely unfathomed and to uncover this phenomenon, nineteen strains of Pantoea agglomerans were selected and analyzed. Genome-to-Genome Distance Calculator (GGDC) which uses the Genome Blast Distance Phylogeny (GBDP) technique to calculate digital DDH values. Phylogenetic analysis via Genome-to-Genome distance, Average Nucleotide Identity, and Amino Acid Identity calculation revealed that all strains belonged to the genus Pantoea. However, strain 33.1 had a lower value than the threshold for the same species delineation. Bacterial Pan Genome Analysis (BPGA) Pipeline and MinPath analysis revealed genetic traits associated with environmental resilience, such as oxidative stress, UV radiation, temperature extremes, and metabolism of distinct host-specific carbohydrates. Protein-protein interactome analysis illustrated osmotic stress proteins closely linked with core proteins, while heavy metal tolerance, nitrogen metabolism, and Type III and VI secretion systems proteins generally associated with pathogenicity formed a separate network, indicating strain-specific characteristics. These findings shed new light on the intricate genetic architecture of Pantoea agglomerans, revealing its adaptability to inhabit diverse niches and thrive in varied environments.

Improving adjustment to daylight saving time transitions with light.

Daylight saving time (DST) is currently utilized in many countries with the rationale that it enhances the alignment between daylight hours and activity peaks in the population. The act of transitioning into and out of DST introduces disruptions to the circadian rhythm, thereby impacting sleep and overall health. Despite the substantial number of individuals affected, the consequences of this circadian disruption have often been overlooked. Here, we employ a mathematical model of the human circadian pacemaker to elucidate how the biological clock interacts with daytime and evening exposures to both natural and electrical light. This interaction plays a crucial role in determining the adaptation to the 1 hour time zone shift imposed by the transition to or from DST. In global discussions about DST, there is a prevailing assumption that individuals easily adjust to DST transitions despite a few studies indicating that the human circadian system requires several days to fully adjust to a DST transition. Our study highlights that evening light exposure changes can be the main driving force for re-entrainment, with chronobiological models predicting that people with longer intrinsic period (i.e. later chronotype) entrain more slowly to transitions to or from DST as compared to people with a shorter intrinsic period (earlier chronotype). Moreover, the model forecasts large inter-individual differences in the adaptation speed, in particular during the spring transition. The predictions derived from our model offer circadian biology-based recommendations for light exposure strategies that facilitate a more rapid adaptation to DST-related transitions or travel across a single time zone. As such, our study contributes valuable insights to the ongoing discourse on DST and its implications for human circadian rhythms.

Design Principles for Perfect Adaptation in Biological Networks with Nonlinear Dynamics.

Establishing a mapping between the emergent biological properties and the repository of network structures has been of great relevance in systems and synthetic biology. Adaptation is one such biological property of paramount importance that promotes regulation in the presence of environmental disturbances. This paper presents a nonlinear systems theory-driven framework to identify the design principles for perfect adaptation with respect to external disturbances of arbitrary magnitude. Based on the prior information about the network, we frame precise mathematical conditions for adaptation using nonlinear systems theory. We first deduce the mathematical conditions for perfect adaptation for constant input disturbances. Subsequently, we translate these conditions to specific necessary structural requirements for adaptation in networks of small size and then extend to argue that there exist only two classes of architectures for a network of any size that can provide local adaptation in the entire state space, namely, incoherent feed-forward (IFF) structure and negative feedback loop with buffer node (NFB). The additional positiveness constraints further narrow the admissible set of network structures. This also aids in establishing the global asymptotic stability for the steady state given a constant input disturbance. The proposed method does not assume any explicit knowledge of the underlying rate kinetics, barring some minimal assumptions. Finally, we also discuss the infeasibility of certain IFF networks in providing adaptation in the presence of downstream connections. Moreover, we propose a generic and novel algorithm based on non-linear systems theory to unravel the design principles for global adaptation. Detailed and extensive simulation studies corroborate the theoretical findings.

Multigenerational Adaptation Can Enhance the Pathogen Resistance of Plants via Changes in Rhizosphere Microbial Community Assembly.

Plants withstand pathogen attacks by recruiting beneficial bacteria to the rhizosphere and passing their legacy on to the next generation. However, the underlying mechanisms involved in this process remain unclear. In our study, we combined microbiomic and transcriptomic analyses to reveal how the rhizosphere microbiome assembled through multiple generations and defense-related genes expressed in Arabidopsis thaliana under pathogen attack stress. Our results showed that continuous exposure to the pathogen Pseudomonas syringae pv tomato DC3000 led to improved growth and increased disease resistance in a third generation of rps2 mutant Arabidopsis thaliana. It could be attributed to the enrichment of specific rhizosphere bacteria, such as Bacillus and Bacteroides. Pathways associated with plant immunity and growth in A. thaliana, such as MAPK signaling pathways, phytohormone signal transduction, ABC transporter proteins, and flavonoid biosynthesis, were activated under the influence of rhizosphere bacterial communities. Our findings provide a scientific basis for explaining the relationship between beneficial microbes and defense-related gene expression. Understanding microbial communities and the mechanisms involved in plant responses to disease can contribute to better plant management and reduction of pesticide use.

RNA-binding protein PCBP2 regulates pancreatic ß cell function and adaptation to glucose.

Glucose plays a key role in shaping pancreatic ß cell function. Thus, deciphering the mechanisms by which this nutrient stimulates ß cells holds therapeutic promise for combating ß cell failure in type 2 diabetes (T2D). ß Cells respond to hyperglycemia in part by rewiring their mRNA metabolism, yet the mechanisms governing these changes remain poorly understood. Here, we identify a requirement for the RNA-binding protein PCBP2 in maintaining ß cell function basally and during sustained hyperglycemic challenge. PCBP2 was induced in primary mouse islets incubated with elevated glucose and was required to adapt insulin secretion. Transcriptomic analysis of primary Pcbp2-deficient ß cells revealed impacts on basal and glucose-regulated mRNAs encoding core components of the insulin secretory pathway. Accordingly, Pcbp2-deficient ß cells exhibited defects in calcium flux, insulin granule ultrastructure and exocytosis, and the amplification pathway of insulin secretion. Further, PCBP2 was induced by glucose in primary human islets, was downregulated in islets from T2D donors, and impacted genes commonly altered in islets from donors with T2D and linked to single-nucleotide polymorphisms associated with T2D. Thus, these findings establish a paradigm for PCBP2 in governing basal and glucose-adaptive gene programs critical for shaping the functional state of ß cells.

Perceptual error based on Bayesian cue combination drives implicit motor adaptation.

The sensorimotor system can recalibrate itself without our conscious awareness, a type of procedural learning whose computational mechanism remains undefined. Recent findings on implicit motor adaptation, such as over-learning from small perturbations and fast saturation for increasing perturbation size, challenge existing theories based on sensory errors. We argue that perceptual error, arising from the optimal combination of movement-related cues, is the primary driver of implicit adaptation. Central to our theory is the increasing sensory uncertainty of visual cues with increasing perturbations, which was validated through perceptual psychophysics (Experiment 1). Our theory predicts the learning dynamics of implicit adaptation across a spectrum of perturbation sizes on a trial-by-trial basis (Experiment 2). It explains proprioception changes and their relation to visual perturbation (Experiment 3). By modulating visual uncertainty in perturbation, we induced unique adaptation responses in line with our model predictions (Experiment 4). Overall, our perceptual error framework outperforms existing models based on sensory errors, suggesting that perceptual error in locating one's effector, supported by Bayesian cue integration, underpins the sensorimotor system's implicit adaptation.

Evolution and host-specific adaptation of Pseudomonas aeruginosa.

The major human bacterial pathogen Pseudomonas aeruginosa causes multidrug-resistant infections in people with underlying immunodeficiencies or structural lung diseases such as cystic fibrosis (CF). We show that a few environmental isolates, driven by horizontal gene acquisition, have become dominant epidemic clones that have sequentially emerged and spread through global transmission networks over the past 200 years. These clones demonstrate varying intrinsic propensities for infecting CF or non-CF individuals (linked to specific transcriptional changes enabling survival within macrophages); have undergone multiple rounds of convergent, host-specific adaptation; and have eventually lost their ability to transmit between different patient groups. Our findings thus explain the pathogenic evolution of P. aeruginosa and highlight the importance of global surveillance and cross-infection prevention in averting the emergence of future epidemic clones.

Flexible adaptation of task-positive brain networks predicts efficiency of evidence accumulation.

Efficiency of evidence accumulation (EEA), an individual's ability to selectively gather goal-relevant information to make adaptive choices, is thought to be a key neurocomputational mechanism associated with cognitive functioning and transdiagnostic risk for psychopathology. However, the neural basis of individual differences in EEA is poorly understood, especially regarding the role of largescale brain network dynamics. We leverage data from 5198 participants from the Human Connectome Project and Adolescent Brain Cognitive Development Study to demonstrate a strong association between EEA and flexible adaptation to cognitive demand in the "task-positive" frontoparietal and dorsal attention networks. Notably, individuals with higher EEA displayed divergent task-positive network activation across n-back task conditions: higher activation under high cognitive demand (2-back) and lower activation under low demand (0-back). These findings suggest that brain networks' flexible adaptation to cognitive demands is a key neural underpinning of EEA.

Transcription factor roles in the local adaptation to temperature in the Andean Spiny Toad Rhinella spinulosa.

Environmental temperature strongly influences the adaptation dynamics of amphibians, whose limited regulation capabilities render them susceptible to thermal oscillations. A central element of the adaptive strategies is the transcription factors (TFs), which act as master regulators that orchestrate stress responses, enabling species to navigate the fluctuations of their environment skillfully. Our study delves into the intricate relationship between TF expression and thermal adaptation mechanisms in the Rhinella spinulosa populations. We sought to elucidate the dynamic modulations of TF expression in prometamorphic and metamorphic tadpoles that inhabit two thermally contrasting environments (Catarpe and El Tatio Geyser, Chile) and which were exposed to two thermal treatments (25 °C vs. 20 °C). Our findings unravel an intriguing dichotomy in response strategies between these populations. First, results evidence the expression of 1374 transcription factors. Regarding the temperature shift, the Catarpe tadpoles show a multifaceted approach by up-regulating crucial TFs, including fosB, atf7, and the androgen receptor. These dynamic regulatory responses likely underpin the population's ability to navigate thermal fluctuations effectively. In stark contrast, the El Tatio tadpoles exhibit a more targeted response, primarily up-regulating foxc1. This differential expression suggests a distinct focus on specific TFs to mitigate the effects of temperature variations. Our study contributes to understanding the molecular mechanisms governing thermal adaptation responses and highlights the resilience and adaptability of amphibians in the face of ever-changing environmental conditions.

Genomics-based identification of a cold adapted clade in Deinococcus.

BACKGROUND: Microbes in the cold polar and alpine environments play a critical role in feedbacks that amplify the effects of climate change. Defining the cold adapted ecotype is one of the prerequisites for understanding the response of polar and alpine microbes to climate change. RESULTS: Here, we analysed 85 high-quality, de-duplicated genomes of Deinococcus, which can survive in a variety of harsh environments. By leveraging genomic and phenotypic traits with reverse ecology, we defined a cold adapted clade from eight Deinococcus strains isolated from Arctic, Antarctic and high alpine environments. Genome-wide optimization in amino acid composition and regulation and signalling enable the cold adapted clade to produce CO2 from organic matter and boost the bioavailability of mineral nitrogen. CONCLUSIONS: Based primarily on in silico genomic analysis, we defined a potential cold adapted clade in Deinococcus and provided an updated view of the genomic traits and metabolic potential of Deinococcus. Our study would facilitate the understanding of microbial processes in the cold polar and alpine environments.

Impact of steroid biosynthesis on the aerobic adaptation of eukaryotes.

Steroids are indispensable components of the eukaryotic cellular membrane and the acquisition of steroid biosynthesis was a key factor that enabled the evolution of eukaryotes. The polycyclic carbon structures of steroids can be preserved in sedimentary rocks as chemical fossils for billions of years and thus provide invaluable clues to trace eukaryotic evolution from the distant past. Steroid biosynthesis consists of (1) the production of protosteroids and (2) the subsequent modifications toward "modern-type" steroids such as cholesterol and stigmasterol. While protosteroid biosynthesis requires only two genes for the cyclization of squalene, complete modification of protosteroids involves ~10 additional genes. Eukaryotes universally possess at least some of those additional genes and thus produce modern-type steroids as major final products. The geological biomarker records suggest a prolonged period of solely protosteroid production in the mid-Proterozoic before the advent of modern-type steroids in the Neoproterozoic. It has been proposed that mid-Proterozoic protosteroids were produced by hypothetical stem-group eukaryotes that presumably possessed genes only for protosteroid production, even though in modern environments protosteroid production as a final product is found exclusively in bacteria. The host identity of mid-Proterozoic steroid producers is crucial for understanding the early evolution of eukaryotes. In this perspective, we discuss how geological biomarker data and genetic data complement each other and potentially provide a more coherent scenario for the evolution of steroids and associated early eukaryotes. We further discuss the potential impacts that steroids had on the evolution of aerobic metabolism in eukaryotes, which may have been an important factor for the eventual ecological dominance of eukaryotes in many modern environments.

Purifying selection drove the adaptation of mitochondrial genes along with correlation of gene rearrangements and evolutionary rates in two subfamilies of Whitefly (Insecta: Hemiptera).

Insect mitochondrial genomes (mitogenomes) are usually represented by a conserved gene order. Whiteflies exhibit gene rearrangement in their mitogenomes; however, understanding how nucleotide substitution rates shape gene rearrangement in whiteflies is unclear due to the limited number of mitogenomes. Additionally, the mechanisms by which selection pressure drives adaptations in mitochondrial genes in the two subfamilies of whiteflies are not yet known. Here, we analyzed 18 whitefly mitogenomes, including one newly generated mitogenome, to compare nucleotide substitution rates, selection pressure, and gene arrangements. The newly generated mitogenome is reported along with reannotation of Pealius mori and comparisons to other whitefly mitogenomes. Comparative studies on nucleotide composition of 18 whiteflies revealed the positive GC skewness, confirming the reversal of strand asymmetry. We found 11 rearranged gene orders within two subfamilies of whiteflies with 8-18 breakpoints of gene rearrangements. Members of the subfamily Aleyrodinae exhibit more complex pathways in the evolution of gene order as compared to the subfamily Aleurodicinae. Our findings also revealed that the increase or reduction of nucleotide substitution rates does not have an impact on any of the gene rearrangement scenarios depicting neutral correlation. Selection pressure analysis revealed that the mitogenomes from members of both the subfamilies Aleurodicinae and Aleyrodinae are characterized by intense purifying selection pressure.

Relative numerosity is constructed from size and density information: Evidence from adaptation.

To dissociate aftereffects of size and density in the perception of relative numerosity, large or small adapter sizes were crossed with high or low adapter densities. A total of 48 participants were included in this preregistered design. To adapt the same retinotopic region as the large adapters, the small adapters were flashed in a sequence so as to "paint" the adapting density across the large region. Perceived numerosities and sizes in the adapted region were then compared to those in an unadapted region in separate blocks of trials, so that changes in density could be inferred. These density changes were found to be bidirectional and roughly symmetric, whereas the aftereffects of size and number were not symmetric. A simple account of these findings is that local adaptations to retinotopic density as well as global adaptations to size combine in producing numerosity aftereffects measured by assessing perceived relative number. Accounts based on number adaptation are contraindicated, in particular, by the result of adapting to a large, sparse adapter and testing with a stimulus with a double the density but half number of dots.

The loci of environmental adaptation in a model eukaryote.

While the underlying genetic changes have been uncovered in some cases of adaptive evolution, the lack of a systematic study prevents a general understanding of the genomic basis of adaptation. For example, it is unclear whether protein-coding or noncoding mutations are more important to adaptive evolution and whether adaptations to different environments are brought by genetic changes distributed in diverse genes and biological processes or concentrated in a core set. We here perform laboratory evolution of 3360 Saccharomyces cerevisiae populations in 252 environments of varying levels of stress. We find the yeast adaptations to be primarily fueled by large-effect coding mutations overrepresented in a relatively small gene set, despite prevalent antagonistic pleiotropy across environments. Populations generally adapt faster in more stressful environments, partly because of greater benefits of the same mutations in more stressful environments. These and other findings from this model eukaryote help unravel the genomic principles of environmental adaptation.

Examining the impact of differing caffeine dosages in conjunction with plyometric training on physiological adaptations in basketball players.

The aim of the current study was to investigate the effects of ingesting different dosages of caffeine (CAF) prior to plyometric jump training (PJT) on sport-related performance and physiological parameters in male basketball players. Twenty-four young athletes were randomly divided into 3 groups and performed 6 weeks of PJT while consuming 3 mg·kg-1 of body mass caffeine (CAF3, n = 8), 6 mg·kg-1 body mass caffeine (CAF6, n = 8) or placebo (PL; n = 8) one hour prior to each training session. Before and after the 6-week PJT, the players were evaluated for field-based basketball-specific performance measures (vertical jump, 20-m sprint, Illinois change of direction speed [CODS], and maximal strength) and lab-based physiological (aerobic capacity and anaerobic power) parameters. CAF3, CAF6, and PL groups demonstrated significant improvements in vertical jump (ES = 1.07, 1.45, and 1.1, respectively), 20-m sprint (ES = - 0.50, - 0.61, and - 0.36), change of direction performance (ES = - 1.22, - 1.26, and - 1.09), maximal strength (ES = 1.68, 2.29, and 1.17), maximum oxygen uptake (V̇O2max) (ES = 1.09, 1.59, and 0.92), and peak (ES = 1.82, 1.85, and 0.82) and average power output (ES = 1.39, 1.32, and 1.07) after 6 weeks of training. Comparative analysis of individual adaptive responses to training indicated that the CAF6 led to insignificantly greater effects in vertical jump (ES = 1.45), maximal strength (ES = 2.29), and V̇O2max (ES = 1.59) with lower residuals in individual changes and lower coefficient of variations (CV) in mean group changes. Regarding sprint and CODS performance, both experimental groups indicated similar changes, residuals in individual changes, and CVs in mean group changes. Overall, consuming 6 mg·kg-1 body mass caffeine induces superior adaptations in aerobic fitness, anaerobic power, and sport-specific performance measures, with lower inter-individual variability in the adaptations and more homogenized changes over the training period.

Genomic sequencing reveals convergent adaptation during experimental evolution in two budding yeast species.

Convergent evolution is central in the origins of multicellularity. Identifying the basis for convergent multicellular evolution is challenging because of the diverse evolutionary origins and environments involved. Haploid Kluyveromyces lactis populations evolve multicellularity during selection for increased settling in liquid media. Strong genomic and phenotypic convergence is observed between K. lactis and previously selected S. cerevisiae populations under similar selection, despite their >100-million-year divergence. We find K. lactis multicellularity is conferred by mutations in genes ACE2 or AIM44, with ACE2 being predominant. They are a subset of the six genes involved in the S. cerevisiae multicellularity. Both ACE2 and AIM44 regulate cell division, indicating that the genetic convergence is likely due to conserved cellular replication mechanisms. Complex population dynamics involving multiple ACE2/AIM44 genotypes are found in most K. lactis lineages. The results show common ancestry and natural selection shape convergence while chance and contingency determine the degree of divergence.

Ecological stress memory in wood architecture of two Neotropical hickory species from central-eastern Mexico.

BACKGROUND: Drought periods are major evolutionary triggers of wood anatomical adaptive variation in Lower Tropical Montane Cloud Forests tree species. We tested the influence of historical drought events on the effects of ecological stress memory on latewood width and xylem vessel traits in two relict hickory species (Carya palmeri and Carya myristiciformis) from central-eastern Mexico. We hypothesized that latewood width would decrease during historical drought years, establishing correlations between growth and water stress conditions, and that moisture deficit during past tree growth between successive drought events, would impact on wood anatomical features. We analyzed latewood anatomical traits that developed during historical drought and pre- and post-drought years in both species. RESULTS: We found that repeated periods of hydric stress left climatic signatures for annual latewood growth and xylem vessel traits that are essential for hydric adaptation in tropical montane hickory species. CONCLUSIONS: Our results demonstrate the existence of cause‒effect relationships in wood anatomical architecture and highlight the ecological stress memory linked with historical drought events. Thus, combined time-series analysis of latewood width and xylem vessel traits is a powerful tool for understanding the ecological behavior of hickory species.

Genome evolution and transcriptome plasticity is associated with adaptation to monocot and dicot plants in Colletotrichum fungi.

BACKGROUND: Colletotrichum fungi infect a wide diversity of monocot and dicot hosts, causing diseases on almost all economically important plants worldwide. Colletotrichum is also a suitable model for studying gene family evolution on a fine scale to uncover events in the genome associated with biological changes. RESULTS: Here we present the genome sequences of 30 Colletotrichum species covering the diversity within the genus. Evolutionary analyses revealed that the Colletotrichum ancestor diverged in the late Cretaceous in parallel with the diversification of flowering plants. We provide evidence of independent host jumps from dicots to monocots during the evolution of Colletotrichum, coinciding with a progressive shrinking of the plant cell wall degradative arsenal and expansions in lineage-specific gene families. Comparative transcriptomics of 4 species adapted to different hosts revealed similarity in gene content but high diversity in the modulation of their transcription profiles on different plant substrates. Combining genomics and transcriptomics, we identified a set of core genes such as specific transcription factors, putatively involved in plant cell wall degradation. CONCLUSIONS: These results indicate that the ancestral Colletotrichum were associated with dicot plants and certain branches progressively adapted to different monocot hosts, reshaping the gene content and its regulation.

Identifying critical growth stage and resilient genotypes in cowpea under drought stress contributes to enhancing crop tolerance for improvement and adaptation in Cameroon.

Drought stress following climate change is likely a scenario that will have to face crop growers in tropical regions. In mitigating this constraint, the best option should be the selection and use of resilient varieties that can withstand drought threats. Therefore, a pot experiment was conducted under greenhouse conditions at the Research and Teaching Farm of the Faculty of Agronomy and Agricultural Sciences of the University of Dschang. The objectives are to identify sensitive growth stage, to identify drought-tolerant genotypes with the help of yield-based selection indices and to identify suitable selection indices that are associated with yield under non-stress and stress circumstances. Eighty-eight cowpea genotypes from the sahelian and western regions of Cameroon were subjected to drought stress at vegetative (VDS) and flowering (FDS) stages by withholding water for 28 days, using a split plot design with two factors and three replications. Seed yields under stress (Ys) and non-stress (Yp) conditions were recorded. Fifteen drought indices were calculated for the two drought stress levels against the yield from non-stress plants. Drought Intensity Index (DII) under VDS and FDS were 0.71 and 0.84 respectively, indicating severe drought stress for both stages. However, flowering stage was significantly more sensitive to drought stress compared to vegetative stage. Based on PCA and correlation analysis, Stress Tolerance Index (STI), Relative Efficiency Index (REI), Geometric Mean Productivity (GMP), Mean Productivity (MP), Yield Index (YI) and Harmonic Mean (HM) correlated strongly with yield under stress and non-stress conditions and are therefore suitable to discriminate high-yielding and tolerant genotypes under both stress and non-stress conditions. Either under VDS and FDS, CP-016 exhibited an outstanding performance under drought stress and was revealed as the most drought tolerant genotype as shown by ranking, PCA and cluster analysis. Taking into account all indices, the top five genotypes namely CP-016, CP-021, MTA-22, CP-056 and CP-060 were identified as the most drought-tolerant genotypes under VDS. For stress activated at flowering stage (FDS), CP-016, CP-056, CP-021, CP-028 and MTA-22 were the top five most drought-tolerant genotypes. Several genotypes with insignificant Ys and irrelevant rank among which CP-037, NDT-001, CP-036, CP-034, NDT-002, CP-031, NDT-011 were identified as highly drought sensitive with low yield stability. This study identified the most sensitive stage and drought tolerant genotypes that are proposed for genetic improvement of cowpea.