Mapping QTLs for drought tolerance in a SEA 5 x AND 277 common bean cross with SSRs and SNP markers

Genetic mapping of loci associated with yield and their components in black common bean (Phaseolus vulgaris L.)

The increase in world population linked to climate change leads to the need to develop more productive and more adapted cultivars of food species. Quantitative trait loci (QTLs) mapping is a useful tool although, interaction between genotype and the environment is still a challenge. In this study, we sought to identify QTL related to grain yield and the production components in common beans (Phaseolus vulgaris L.) supported by QTL × $\times $ environment interaction. Two hundred eight recombinant inbred lines obtained from the Awauna UEM × $ \times $ IPR88 Uirapuru common bean cross were evaluated in 2017, 2018, and 2019 in field conditions under a 15 × 15 triple lattice experimental design. QTL mapping was estimated using genotypic means and a genetic linkage map with 288 single nucleotide polymorphism markers. Five QTLs associated with plant height (PH), number of pods per plant (NPP), first pod height (FPH), 100-seed weight (SW), and grain yield per plant (GYP) were identified on chromosomes Pv01, Pv04, Pv08, and Pv10. Interestingly, three of these QTLs were co-localized for more than one trait, where the QTL for PH, NPP, and GYP co-locate on Pv01, the QTL for PH and FPH co-locate on Pv04, and the QTL for NPP and SW co-locate on Pv08. In turn, on Pv10, two distinct QTLs were found for SW. The identification of these QTLs stands out in Brazil since relatively little research is directed at this economically important commercial group. It is noteworthy that the molecular markers found linked to the QTLs must later be validated to be used in a multi-trait marker-assisted selection.

Whole plant transpiration responses of common bean (Phaseolus vulgaris L.) to drying soil: Water channels and transcription factors

Common bean (Phaseolus vulgaris L.) is the main legume crop for direct human consumption worldwide. Among abiotic factors affecting common bean, drought is the most limiting. This study aimed at characterizing genetic variability and architecture of transpiration, stomatal regulation and whole plant water use within the Mesoamerican germplasm. A critical fraction of transpirable soil water (FTSWc) was estimated as the inflection point at which NTR starts decreasing linearly. Genome-wide association (GWA) analyses for mean NTR and FTSWc were performed. High variation on mean NTR and FTSWc was found among genotypes. Unreported genomic signals controlling the variation of these traits were identified on Pv01 and Pv07 some located in intergenic, intronic and exonic regions. A set of novel candidate genes and putative regulatory elements located in these QTL were identified. Some of the genes have been previously reported to be involved in abiotic tolerance in model species, including some of the five transcription factors (TF) identified. Four candidate genes, one with potential water transportation activity and three TFs were validated. The gene Phvul.001G108800, an aquaporin SIP2-1 related gene, showed water channel activity through oocyte water assays. Mutant Arabidopsis thaliana (Ath) lines for the homologous genes of common bean were evaluated in transpiration experiments. Two of the three evaluated TFs, UPBEAT1 and C2H2-type ZN finger protein, were involved in the control of transpiration responses to drying soil. Our results provide evidence of novel genes to accelerate the drought tolerance improvement in the crop and study the physiological basis of drought response in plants.

A review on strategies for crop improvement against drought stress through molecular insights

The demand for food goods is rising along with the world population growth, which is directly related to the yield of agricultural crops around the world. However, a number of environmental factors, including floods, salinity, moisture, and drought, have a detrimental effect on agricultural production around the world. Among all of these stresses, drought stress (DS) poses a constant threat to agricultural crops and is a significant impediment to global agricultural productivity. Its potency and severity are expected to increase in the future years. A variety of techniques have been used to generate drought-resistant plants in order to get around this restriction. Different crop plants exhibit specific traits that contribute to drought resistance (DR), such as early flowering, drought escape (DE), and leaf traits. We are highlighting numerous methods that can be used to overcome the effects of DS in this review. Agronomic methods, transgenic methods, the use of sufficient fertilizers, and molecular methods such as clustered regularly interspaced short palindromic repeats (CRISPRs)-associated nuclease 9 (Cas9), virus-induced gene silencing (VIGS), quantitative trait loci (QTL) mapping, microRNA (miRNA) technology, and OMICS-based approaches make up the majority of these techniques. CRISPR technology has rapidly become an increasingly popular choice among researchers exploring natural tolerance to abiotic stresses although, only a few plants have been produced so far using this technique. In order to address the difficulties imposed by DS, new plants utilizing the CRISPR technology must be developed.

Genome-wide association mapping for component traits of drought tolerance in dry beans (Phaseolus vulgaris L.)

Understanding the genetic basis of traits of economic importance under drought stressed and well-watered conditions is important in enhancing genetic gains in dry beans (Phaseolus vulgaris L.). This research aims to: (i) identify markers associated with agronomic and physiological traits for drought tolerance and (ii) identify drought-related putative candidate genes within the mapped genomic regions. An andean and middle-american diversity panel (AMDP) comprising of 185 genotypes was screened in the field under drought stressed and well-watered conditions for two successive seasons. Agronomic and physiological traits, viz., days to 50% flowering (DFW), plant height (PH), days to physiological maturity (DPM), grain yield (GYD), 100-seed weight (SW), leaf temperature (LT), leaf chlorophyll content (LCC) and stomatal conductance (SC) were phenotyped. Principal component and association analysis were conducted using the filtered 9370 Diversity Arrays Technology sequencing (DArTseq) markers. The mean PH, GYD, SW, DPM, LCC and SC of the panel was reduced by 12.1, 29.6, 10.3, 12.6, 28.5 and 62.0%, respectively under drought stressed conditions. Population structure analysis revealed two sub-populations, which corresponded to the andean and middle-american gene pools. Markers explained 0.08-0.10, 0.22-0.23, 0.29-0.32, 0.43-0.44, 0.65-0.66 and 0.69-0.70 of the total phenotypic variability (R2) for SC, LT, PH, GYD, SW and DFW, respectively under drought stressed conditions. For well-watered conditions, R2 varied from 0.08 (LT) to 0.70 (DPM). Overall, 68 significant (p < 10-03) marker-trait associations (MTAs) and 22 putative candidate genes were identified across drought stressed and well-watered conditions. Most of the identified genes had known biological functions related to regulating the response to drought stress. The findings provide new insights into the genetic architecture of drought stress tolerance in common bean. The findings also provide potential candidate SNPs and putative genes that can be utilized in gene discovery and marker-assisted breeding for drought tolerance after validation.

The aquaporin gene PvXIP1;2 conferring drought resistance identified by GWAS at seedling stage in common bean

A whole-genome resequencing-derived SNP dataset used for genome-wide association analysis revealed 12 loci significantly associated with drought stress based on survival rate after drought stress at seedling stage. We further confirmed the drought-related function of an aquaporin gene (PvXIP1;2) located at Locus_10. A variety of adverse conditions, including drought stress, severely affect common bean production. Molecular breeding for drought resistance has been proposed as an effective and practical way to improve the drought resistance of common bean. A genome-wide association analysis was conducted to identify drought-related loci based on survival rates at the seedling stage using a natural population consisting of 400 common bean accessions and 3,832,340 SNPs. The coefficient of variation ranged from 40.90 to 56.22% for survival rates in three independent experiments. A total of 12 associated loci containing 89 significant SNPs were identified for survival rates at the seedling stage. Four loci overlapped in the region of the QTLs reported to be associated with drought resistance. According to the expression profiles, gene annotations and references of the functions of homologous genes in Arabidopsis, 39 genes were considered potential candidate genes selected from 199 genes annotated within all associated loci. A stable locus (Locus_10) was identified on chromosome 11, which contained LEA, aquaporin, and proline-rich protein genes. We further confirmed the drought-related function of an aquaporin (PvXIP1;2) located at Locus_10 by expression pattern analysis, phenotypic analysis of PvXIP1;2-overexpressing Arabidopsis and Agrobacterium rhizogenes-mediated hairy root transformation systems, indicating that the association results can facilitate the efficient identification of genes related to drought resistance. These loci and their candidate genes provide a foundation for crop improvement via breeding for drought resistance in common bean.

Genome-Wide Association Study Reveals Genomic Regions Associated with Fusarium Wilt Resistance in Common Bean

Fusarium wilt (Fusarium oxysporum f. sp. phaseoli, Fop) is one of the main fungal soil diseases in common bean. The aim of the present study was to identify genomic regions associated with Fop resistance through genome-wide association studies (GWAS) in a Mesoamerican Diversity Panel (MDP) and to identify potential common bean sources of Fop's resistance. The MDP was genotyped with BARCBean6K_3BeadChip and evaluated for Fop resistance with two different monosporic strains using the root-dip method. Disease severity rating (DSR) and the area under the disease progress curve (AUDPC), at 21 days after inoculation (DAI), were used for GWAS performed with FarmCPU model. The p-value of each SNP was determined by resampling method and Bonferroni test. For UFV01 strain, two significant single nucleotide polymorphisms (SNPs) were mapped on the Pv05 and Pv11 for AUDPC, and the same SNP (ss715648096) on Pv11 was associated with AUDPC and DSR. Another SNP, mapped on Pv03, showed significance for DSR. Regarding IAC18001 strain, significant SNPs on Pv03, Pv04, Pv05, Pv07 and on Pv01, Pv05, and Pv10 were observed. Putative candidate genes related to nucleotide-binding sites and carboxy-terminal leucine-rich repeats were identified. The markers may be important future tools for genomic selection to Fop disease resistance in beans.

Genetic dissection of drought resistance based on root traits at the bud stage in common bean

A whole-genome resequencing-derived SNP dataset used for genome-wide association analysis revealed 196 loci significantly associated with drought stress based on root traits. Candidate genes identified in the regions of these loci include homologs of known drought resistance genes in A. thaliana. Drought is the main abiotic constraint of the production of common bean. Improved adaptation to drought environments has become a main goal of crop breeding due to the increasing scarcity of water that will occur in the future. The overall objective of our study was to identify genomic regions associated with drought resistance based on root traits using genome-wide association analysis. A natural population of 438 common bean accessions was evaluated for root traits: root surface area, root average diameter, root volume, total root length, taproot length, lateral root number, root dry weight, lateral root length, special root weight/length, using seed germination pouches under drought conditions and in well-watered environments. The coefficient of variation ranged from 11.24% (root average diameter) to 38.19% (root dry weight) in the well-watered environment and from 9.61% (root average diameter) to 39.05% (lateral root length) under drought stress. A whole-genome resequencing-derived SNP dataset revealed 196 loci containing 230 candidate SNPs associated with drought resistance. Seventeen candidate SNPs were simultaneously associated with more than two traits. Forty-one loci were simultaneously associated with more than two traits, and eleven loci were colocated with loci previously reported to be related to drought resistance. Candidate genes of the associated loci included the ABA-responsive element-binding protein family, MYB, NAC, the protein kinase superfamily, etc. These results revealed promising alleles linked to drought resistance or root traits, providing insights into the genetic basis of drought resistance and roots, which will be useful for common bean improvement.

Common bean SNP alleles and candidate genes affecting photosynthesis under contrasting water regimes

Water deficit is a major worldwide constraint to common bean (Phaseolus vulgaris L.) production, being photosynthesis one of the most affected physiological processes. To gain insights into the genetic basis of the photosynthetic response of common bean under water-limited conditions, a collection of 158 Portuguese accessions was grown under both well-watered and water-deficit regimes. Leaf gas-exchange parameters were measured and photosynthetic pigments quantified. The same collection was genotyped using SNP arrays, and SNP-trait associations tested considering a linear mixed model accounting for the genetic relatedness among accessions. A total of 133 SNP-trait associations were identified for net CO2 assimilation rate, transpiration rate, stomatal conductance, and chlorophylls a and b, carotenes, and xanthophyll contents. Ninety of these associations were detected under water-deficit and 43 under well-watered conditions, with only two associations common to both treatments. Identified candidate genes revealed that stomatal regulation, protein translocation across membranes, redox mechanisms, hormone, and osmotic stress signaling were the most relevant processes involved in common bean response to water-limited conditions. These candidates are now preferential targets for common bean water-deficit-tolerance breeding. Additionally, new sources of water-deficit tolerance of Andean, Mesoamerican, and admixed origin were detected as accessions valuable for breeding, and not yet explored.

Genome-Wide Association Studies Detect Multiple QTLs for Productivity in Mesoamerican Diversity Panel of Common Bean Under Drought Stress

Drought stress is an important abiotic factor limiting common bean yield, with great impact on the production worldwide. Understanding the genetic basis regulating beans' yield and seed weight (SW) is a fundamental prerequisite for the development of superior cultivars. The main objectives of this work were to conduct genome-wide marker discovery by genotyping a Mesoamerican panel of common bean germplasm, containing cultivated and landrace accessions of broad origin, followed by the identification of genomic regions associated with productivity under two water regimes using different genome-wide association study (GWAS) approaches. A total of 11,870 markers were genotyped for the 339 genotypes, of which 3,213 were SilicoDArT and 8,657 SNPs derived from DArT and CaptureSeq. The estimated linkage disequilibrium extension, corrected for structure and relatedness (r ² _sv ), was 98.63 and 124.18 kb for landraces and breeding lines, respectively. Germplasm was structured into landraces and lines/cultivars. We carried out GWASs for 100-SW and yield in field environments with and without water stress for 3 consecutive years, using single-, segment-, and gene-based models. Higher number of associations at high stringency was identified for the SW trait under irrigation, totaling ∼185 QTLs for both single- and segment-based, whereas gene-based GWASs showed ∼220 genomic regions containing ∼650 genes. For SW under drought, 18 QTLs were identified for single- and segment-based and 35 genes by gene-based GWASs. For yield, under irrigation, 25 associations were identified, whereas under drought the total was 10 using both approaches. In addition to the consistent associations detected across experiments, these GWAS approaches provided important complementary QTL information (∼221 QTLs; 650 genes; r ² from 0.01% to 32%). Several QTLs were mined within or near candidate genes playing significant role in productivity, providing better understanding of the genetic mechanisms underlying these traits and making available molecular tools to be used in marker-assisted breeding. The findings also allowed the identification of genetic material (germplasm) with better yield performance under drought, promising to a common bean breeding program. Finally, the availability of this highly diverse Mesoamerican panel is of great scientific value for the analysis of any relevant traits in common bean.

Genomic Prediction of Agronomic Traits in Common Bean (Phaseolus vulgaris L.) Under Environmental Stress

In plant and animal breeding, genomic prediction models are established to select new lines based on genomic data, without the need for laborious phenotyping. Prediction models can be trained on recent or historic phenotypic data and increasingly available genotypic data. This enables the adoption of genomic selection also in under-used legume crops such as common bean. Beans are an important staple food in the tropics and mainly grown by smallholders under limiting environmental conditions such as drought or low soil fertility. Therefore, genotype-by-environment interactions (G × E) are an important consideration when developing new bean varieties. However, G × E are often not considered in genomic prediction models nor are these models implemented in current bean breeding programs. Here we show the prediction abilities of four agronomic traits in common bean under various environmental stresses based on twelve field trials. The dataset includes 481 elite breeding lines characterized by 5,820 SNP markers. Prediction abilities over all twelve trials ranged between 0.6 and 0.8 for yield and days to maturity, respectively, predicting new lines into new seasons. In all four evaluated traits, the prediction abilities reached about 50-80% of the maximum accuracies given by phenotypic correlations and heritability. Predictions under drought and low phosphorus stress were up to 10 and 20% improved when G × E were included in the model, respectively. Our results demonstrate the potential of genomic selection to increase the genetic gain in common bean breeding. Prediction abilities improved when more phenotypic data was available and G × E could be accounted for. Furthermore, the developed models allowed us to predict genotypic performance under different environmental stresses. This will be a key factor in the development of common bean varieties adapted to future challenging conditions.

Exploration of the Yield Potential of Mesoamerican Wild Common Beans From Contrasting Eco-Geographic Regions by Nested Recombinant Inbred Populations

Genetic analyses and utilization of wild genetic variation for crop improvement in common bean (Phaseolus vulgaris L.) have been hampered by yield evaluation difficulties, identification of advantageous variation, and linkage drag. The lack of adaptation to cultivation conditions and the existence of highly structured populations make association mapping of diversity panels not optimal. Joint linkage mapping of nested populations avoids the later constraint, while populations crossed with a common domesticated parent allow the evaluation of wild variation within a more adapted background. Three domesticated by wild backcrossed-inbred-line populations (BC₁S₄) were developed using three wild accessions representing the full range of rainfall of the Mesoamerican wild bean distribution crossed to the elite drought tolerant domesticated parent SEA 5. These populations were evaluated under field conditions in three environments, two fully irrigated trials in two seasons and a simulated terminal drought in the second season. The goal was to test if these populations responded differently to drought stress and contained progenies with higher yield than SEA 5, not only under drought but also under water-watered conditions. Results revealed that the two populations derived from wild parents of the lower rainfall regions produced lines with higher yield compared to the domesticated parent in the three environments, i.e., both in the drought-stressed environment and in the well-watered treatments. Several progeny lines produced yields, which on average over the three environments were 20% higher than the SEA 5 yield. Twenty QTLs for yield were identified in 13 unique regions on eight of the 11 chromosomes of common bean. Five of these regions showed at least one wild allele that increased yield over the domesticated parent. The variation explained by these QTLs ranged from 0.6 to 5.4% of the total variation and the additive effects ranged from -164 to 277 kg ha^-1, with evidence suggesting allelic series for some QTLs. Our results underscore the potential of wild variation, especially from drought-stressed regions, for bean crop improvement as well the identification of regions for efficient marker-assisted introgression.

Effect of Water Deficit on Morphoagronomic and Physiological Traits of Common Bean Genotypes with Contrasting Drought Tolerance

Water deficit is considered one of the most limiting factors of the common bean. Understanding the adaptation mechanisms of the crop to this stress is fundamental for the development of drought-tolerant cultivars. In this sense, the objective of this study was to analyze the influence of water deficit on physiological and morphoagronomic traits of common bean genotypes with contrasting drought tolerance, aiming to identify mechanisms associated with tolerance to water deficit. The experiment was carried out in a greenhouse, arranged in a randomized complete block 4 × 2 factorial design, consisting of four common bean genotypes under two water regimes (with and without water stress), with six replications. The morphoagronomic and physiological traits of four cultivars, two drought-tolerant (IAPAR 81 and BAT 477) and two drought-sensitive (IAC Tybatã and BRS Pontal), were measured for 0, 4, 8, and 12 days, under water deficit, initiated in the phenological stage R5. Water-deficit induced physiological changes in the plants, altering the evaluated morphoagronomic traits. The drought tolerance of cultivar BAT 477 is not only a direct result of the low influence of water deficit on its yield components, but also a consequence of the participation of multiple adaptive physiological mechanisms, such as higher intrinsic water use efficiency, net photosynthesis rate, transpiration, carboxylation efficiency, stomatal conductance, and intracellular concentration of CO2 under water deficit conditions. On the other hand, cultivar IAPAR 81 can be considered drought-tolerant for short water-deficit periods only, since after the eighth day of water deficit, the physiological activities decline drastically.

Construction, characteristics and high throughput molecular screening methodologies in some special breeding populations: a horticultural perspective

Advanced marker technologies are widely used for evaluation of genetic diversity in cultivated crops, wild ancestors, landraces or any special plant genotypes. Developing agricultural cultivars requires the following steps: (i) determining desired characteristics to be improved, (ii) screening genetic resources to help find a superior cultivar, (iii) intercrossing selected individuals, (iv) generating genetically hybrid populations and screening them for agro-morphological or molecular traits, (v) evaluating the superior cultivar candidates, (vi) testing field performance at different locations, and (vii) certifying. In the cultivar development process valuable genes can be identified by creating special biparental or multiparental populations and analysing their association using suitable markers in given populations. These special populations and advanced marker technologies give us a deeper knowledge about the inherited agronomic characteristics. Unaffected by the changing environmental conditions, these provide a higher understanding of genome dynamics in plants. The last decade witnessed new applications for advanced molecular techniques in the area of breeding,with low costs per sample. These, especially, include next-generation sequencing technologies like reduced representation genome sequencing (genotyping by sequencing, restriction site-associated DNA). These enabled researchers to develop new markers, such as simple sequence repeat and single- nucleotide polymorphism, for expanding the qualitative and quantitative information onpopulation dynamics. Thus, the knowledge acquired from novel technologies is a valuable asset for the breeding process and to better understand the population dynamics, their properties, and analysis methods.

Research Progress and Perspective on Drought Stress in Legumes: A Review

Climate change, food shortage, water scarcity, and population growth are some of the threatening challenges being faced in today's world. Drought stress (DS) poses a constant challenge for agricultural crops and has been considered a severe constraint for global agricultural productivity; its intensity and severity are predicted to increase in the near future. Legumes demonstrate high sensitivity to DS, especially at vegetative and reproductive stages. They are mostly grown in the dry areas and are moderately drought tolerant, but severe DS leads to remarkable production losses. The most prominent effects of DS are reduced germination, stunted growth, serious damage to the photosynthetic apparatus, decrease in net photosynthesis, and a reduction in nutrient uptake. To curb the catastrophic effect of DS in legumes, it is imperative to understand its effects, mechanisms, and the agronomic and genetic basis of drought for sustainable management. This review highlights the impact of DS on legumes, mechanisms, and proposes appropriate management approaches to alleviate the severity of water stress. In our discussion, we outline the influence of water stress on physiological aspects (such as germination, photosynthesis, water and nutrient uptake), growth parameters and yield. Additionally, mechanisms, various management strategies, for instance, agronomic practices (planting time and geometry, nutrient management), plant growth-promoting Rhizobacteria and arbuscular mycorrhizal fungal inoculation, quantitative trait loci (QTLs), functional genomics and advanced strategies (CRISPR-Cas9) are also critically discussed. We propose that the integration of several approaches such as agronomic and biotechnological strategies as well as advanced genome editing tools is needed to develop drought-tolerant legume cultivars.

Quantitative trait loci (QTL) mapping for intermittent drought tolerance in BRB 191 × SEQ 1027 Andean Intragene cross recombinant inbred line population of common bean (Phaseolus vulgaris L.)

Drought is a major constraint of common bean (Phaseolus vulgaris L.) production in Uganda where irrigation for the crop is very uncommon. This study aimed to identify quantitative trait loci (QTLs) underlying drought tolerance in 128 F5 RILs derived from an Andean intra-gene cross between drought-tolerant SEQ 1027 and BRB 191. Eighteen traits were evaluated under drought stress and non-stress conditions in the field for 2 years and in the greenhouse for 1 year, respectively. A linkage map spanning 486.29 cM was constructed using 53 single nucleotide polymorphic markers (SNP) markers obtained from the KASP genotyping assay. Eleven consistent QTLs were detected on five linkage groups at a threshold of Logarithm of Odds (LOD) ≥ 3.0. Four QTLs were constitutive, seven were adaptive and were associated with 100 seed weight, grain yield, chlorophyll content, harvest index, dry weight of leaf and stem biomass and yield production efficiency. The QTL associated with a 100 seed weight (sw3.1^BS) was the most consistent with the highest percentage of variation explained (21%). Co-localization of five drought-related factors QTLs was detected on pv10 suggesting pleiotropic effects on this chromosome. Identification of molecular markers closely linked to the QTLs identified in this study will facilitate marker assisted breeding for drought tolerance.

Physiological, biochemical and morphoagronomic characterization of drought-tolerant and drought-sensitive bean genotypes under water stress

Water deficit is one of the main factors that reduce grain yield. A better understanding of the mechanisms related to this abiotic stress is a key aspect to design and act upon drought tolerance improvement in crop plants. Therefore, the major objective of this study was to investigate four common bean genotypes for drought tolerance and to establish their tolerance mechanisms. The experiment was carried out in a greenhouse, using the completely randomized design in a factorial arrangement (2 × 4), composed by 2 water conditions (well-watered and water deficit) and 4 cultivars, with six replicates per treatment. The four cultivars, two drought-sensitive (IAC Tybatã and BRS Pontal) and two drought-tolerant (IAPAR 81 and BAT 477), were evaluated for some physiological, biochemical and morphoagronomic traits. Drought promoted physiological and metabolic changes in the plants, reflecting on the morphoagronomic traits. Under water deficit, the genotype IAPAR 81 stood out from the others in terms of physiological characters, however, it presented a low efficiency concerning biochemical activities and a significant reduction in the morphoagronomic characters. The cultivar BAT 477 demonstrated to be drought-adapted presenting more efficient biochemical and morphoagronomic adaptions and the genotype BRS Pontal obtained morphoagronomic values similar to BAT 477, thus it may be classified as moderately tolerant to drought.