Next Article in Journal
Narrow-Band 311 nm Ultraviolet-B Radiation Evokes Different Antioxidant Responses from Broad-Band Ultraviolet
Next Article in Special Issue
Genome-Environment Association Analysis for Bio-Climatic Variables in Common Bean (Phaseolus vulgaris L.) from Brazil
Previous Article in Journal
Identification and Characterization of Arbutus unedo L. Endophytic Bacteria Isolated from Wild and Cultivated Trees for the Biological Control of Phytophthora cinnamomi
Previous Article in Special Issue
Differential Regulation of Drought Responses in Two Phaseolus vulgaris Genotypes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genetic Mapping for Agronomic Traits in IAPAR 81/LP97-28 Population of Common Bean (Phaseolus vulgaris L.) under Drought Conditions

by
Júlio César Ferreira Elias
1,
Maria Celeste Gonçalves-Vidigal
1,*,
Mariana Vaz Bisneta
1,
Giseli Valentini
2,
Pedro Soares Vidigal Filho
1,
Thiago Alexandre Santana Gilio
3,
Vânia Moda-Cirino
4 and
Qijian Song
2
1
Departamento de Agronomia, Universidade Estadual de Maringá—UEM, Av. Colombo 5790, Maringá 87020-900, PR, Brazil
2
Soybean Genomics and Improvement Laboratory, US Department of Agriculture, Agricultural Research Service (USDA-ARS), Beltsville, MD 20705, USA
3
Programa de Pós-graduação em Genética e Melhoramento de Plantas, Universidade do Estado de Mato Grosso, Cáceres 78217-900, MT, Brazil
4
Instituto Agronômico do Paraná—IAPAR, Rua Celso Garcia Cid, km 375, Londrina 86047-902, PR, Brazil
*
Author to whom correspondence should be addressed.
Plants 2021, 10(8), 1568; https://doi.org/10.3390/plants10081568
Submission received: 3 July 2021 / Revised: 24 July 2021 / Accepted: 28 July 2021 / Published: 30 July 2021
(This article belongs to the Special Issue Drought Tolerance in Common Bean)

Abstract

:
One of the significant challenges of common bean breeding is developing cultivars with high yields under drought conditions. The present study attempted to map quantitative trait loci (QTLs) and identify molecular markers that are linked to drought tolerance in the common bean. We evaluated 160 recombinant inbred lines (RILs), derived from the cross between the carioca cultivars IAPAR 81 (drought tolerant) and LP97-28 (susceptible to drought). In 2014 and 2015, two experiments were conducted (DS—drought stress, and NS—no drought stress). In the DS experiment, water suppression was performed at the flowering stages R5 to R6. The results of our experiments showed that drought conditions play an essential role in reducing most of the traits that were evaluated. RILs under drought conditions reduced the grain yield by 62.03% and 24% in 2014 and 2015, respectively. We identified 15 quantitative trait loci distributed on the chromosomes Pv01, Pv02, Pv03, Pv07, Pv08, Pv09, Pv10, and Pv11, related to grain yield, seed yield per day, 100-seed weight, number of pods per plant, plant height, number of days for flowering, and number of days to maturity. The characteristics of seed yield per day, 100-seed weight, and number of days to maturity showed that QTLs colocalized on Pv07. Identifying QTLs that are linked to drought tolerance in the RIL population IAPAR 81 × LP97-28 is of particular importance for common bean breeding programs seeking to improve carioca beans that are cultivated in regions with drought conditions, such as Brazil.

1. Introduction

The common bean is an essential legume for the human diet, as it represents a vital source of nutrition and income for people, especially those from African and Latin American countries [1,2,3,4]. The common bean is considered to be a subsistence crop in several regions globally, especially regions that are characterized by soil with low-quality, adequate fertilization, and susceptibility to biotic and abiotic stresses. Among these stresses, a lack of rain is the primary cause of the low productivity of the common bean, since the cultivation of the mentioned legume relies mainly on rainfall.
Water deficit in the common bean (Phaseolus vulgaris L.) causes severe restrictions to productivity, which can directly affect the morphophysiology of the crop during any cycle stage. The incidence of short-term drought periods in the main common bean-producing regions, such as Brazil, Mexico, the United States, and African countries, tends to increase due to climate change, which negatively impacts global production and food security [5,6,7]. According to Diaz [2], drought represents a high-risk factor for bean production, and improved varieties that withstand the altered climatic conditions constitute a valuable resource. Therefore, efforts are necessary to improve the crop tolerance to stress conditions, particularly drought, to reduce poverty and malnutrition in developing countries in tropical regions [8].
One of the most significant breeding program challenges is developing cultivars with a substantial production capacity under water stress conditions [9]. In this context, plant breeding that is focused on drought tolerance offers a useful, efficient, and enduring method to ensure sustainable common bean productivity, even under unfavorable water-deficit conditions [8,10].
Brazil is the third largest global producer, and one of the largest consumer markets, of grains from the Phaseolus genus, characterized either by subsistence or mechanized agriculture. In both the scenarios, abiotic factors (e.g., high aluminum toxicity, sandy soil with low nutrients, and drought) can restrict bean cultivation, and consequently reduce bean yield [11].
Common bean productivity in Brazil can vary according to the season, and this is closely related to the rain supply. During rainy-season production, sowing occurs in the months that show a high frequency of rain, with the harvest occurring between November and February. Dry-season production, with irrigation supplementation, is concentrated in some regions of Brazil, mainly for the central region, where highly technical cropping systems and central-pivot irrigation have been adopted. Winter, or third-season, production takes place from May to September. This scenario shows that beans can be grown year round, guaranteeing the permanent supply of this product on the market; however, the average production is highly unstable, and varies according to the growing region [12].
The common beans of Mesoamerican origin are the most cultivated in Brazil [13], because of the consumer-preferred grains of black and carioca commercial types. The vast diversity that is present in Mesoamerican and Andean genotypes provides levels of drought tolerance. Consequently, common bean breeders can profit from these germplasms, to develop drought-tolerant cultivars [8,14].
Mesoamerica and Durango races have been frequently described as drought-tolerance sources to develop cultivars for environments that are affected by abiotic stress [15,16,17]. Despite that, studies have revealed little progress in the transference of drought-tolerant genes to commercial cultivars, through hybridization between the two bean gene pools [18,19].
Biomass accumulation in aerial parts and roots, translocation of photoassimilates to seeds, water stress harvesting, and pod harvest indices are examples of traits that are associated with grain production stability under water stress [8,20,21,22]. An optimal strategy to tolerate drought would be a deep root system with long secondary and tertiary roots that allow for water extraction from the deeper soil layers [23,24,25].
The change in characteristics, such as the number of days for flowering and maturation for water supply, is an essential strategy for adapting common beans to water shortages. Hoyos-Villegas et al. [26] reported that water conservation might limit photosynthesis, and the growth and development of plants. In response to severe drought stress, bean plants may produce smaller or fewer leaves, resulting in suboptimal leaf area and reduced net photosynthesis. Schneider et al. [27] suggested selection for increased biomass under drought stress in order to avoid a reduction in photosynthetic capacity. This practice, however, may result in indirect selection for later maturity.
Beebe et al. [1] noted that mechanisms to escape drought include early maturity, phenotypic plasticity, and rapid partitioning of photosynthates to seeds. Selection for earlier maturity may help avoid terminal drought, but early maturity may reduce the seed yield potential during more favorable growing seasons. An early, defined, and irreversible shift to reproductive development, and a shorter period of pod filling, could reduce the exposure during the sensitive reproductive period of growth, and shorten the growing season, thus increasing the chances of escape [28]. The characteristics that are related to drought tolerance are usually controlled by the genetic factors of polygenic inheritance [7,10,27]. However, further detailed studies are needed in order to better understand all the mechanisms that are involved in drought tolerance, and ease gene transference to new cultivars.
Efficient water use strategies, through root systems and crop development, are crucial factors for common bean adaptation to drought conditions [15,29]. Therefore, drought affects common beans directly, by reducing their growth and causing plant death [29]. Likewise, a plant under a high water stress index becomes more susceptible to plagues and diseases [30]. In contrast, a healthy plant, with genetic mechanisms for drought tolerance, could maintain production stability under drought conditions. The genetic control of the traits that determine the drought tolerance of the common bean may suffer variations according to the stress conditions to which the genotypes are submitted [1,31].
The investigation of the genetic control of this type of trait can be performed by genetic mapping and association analysis, aiming to identify genomic regions controlling phenotypic traits. Until recently, the common bean genotyping ability was limited, which compromised the accuracy of the QTLs that were mapped for water stress, but this situation has changed in recent years, with technological evolution in genotyping, such as the bead chip, genotyping by sequencing, reference genome, and phenotyping advances.
Traub et al. [9] evaluated the metabolic and photosynthetic responses to drought stress in P. vulgaris and P.acutifolius. Mukeshimana et al. [8] identified QTLs that were associated with the yield components under drought tolerance in an intergenic common bean population (SEA 5 × CAL 93), on chromosomes Pv03, Pv07, and Pv09. Asfaw et al. [24] identified QTLs for the characteristic yield and SPAD chlorophyll meter reading on Pv08 and Pv06, respectively, under water stress conditions. Additionally, with the Mesoamerican population from the Buster × Roza cross, significant QTLs were mapped on the chromosomes Pv01 and Pv02, for yield and phenological traits [10]. By nested association mapping (NAM), Hoyos-Villegas et al. [26] reported QTNs on Pv03 and Pv10, with effects that were associated with yield components under drought stress. Dramadri et al. [32], using an RIL population from Ecuador (drought tolerant) × Red Hawk, mapped eighteen QTLs that were associated with the phenology, yield components, and partitioning traits on the chromosomes Pv01, Pv02, Pv03, Pv04, Pv06, and Pv11. The purpose of this study was to recognize essential traits under water stress, define the genetic effects of drought stress in 160 common bean recombinant inbred lines (RILs) from cross IAPAR 81 (tolerant) × LP 97-28 (susceptible), and identify qualitative trait loci (QTLs) that are related to drought-tolerant traits in this population.

2. Results

2.1. Level of Water Stress and Its Effects on Genotype Performance

In the area where we performed the experiments, during 2014 and 2015, the precipitation was insufficient for the optimal development of common bean plants. Consequently, we implemented irrigation water management. The average temperatures were suitable for normal plant development during the same period; the average daily temperature in 2014 was 23.73 °C, whereas in 2015, it was 23.45 °C (Figure 1).
The precipitation in 2014 was 220 mm, of which 132 mm accumulated until the beginning of the flowering season (R2). Furthermore, in 2015, a total of 599 mm of accumulated rainfall was observed during the whole experimental period. Only 386 mm was observed during the maturation phase. Thus, the experiment described as DS (with drought stress), in 2014 and 2015, suffered a sporadic intraseasonal drought at the beginning of flowering, contributing to an adequate evaluation of the drought effect on the bean lines at the flowering stage.
The drought intensity index (DII) evaluates the behavior of genotypes under DS (with drought stress) and NS (no drought stress) conditions, as shown in Table 1. The drought intensity index for the grain yield was moderate, and the average value for 2014 and 2015 was DII = 0.44. The drought effect on grain yield was more severe in 2014, with DII = 0.59, whereas in 2015, this index value was 0.30. The parents exhibited mean index values of DII = 0.70 and DII = 0.38 in 2014 and 2015, respectively. However, contrasting DII values were observed for the parents in 2015; IAPAR 81 (tolerant to drought) showed a low drought index (DII = 0.21), but LP 97-28 (not susceptible to drought) underwent severe drought and exhibited a DII value of 0.50.

2.2. Phenotypic Data

In 2014, all the evaluated traits were significantly affected by the drought stress conditions during the flowering stage, except for the number of seeds per pod (Table 2). The average grain yields of the 160 recombinant inbred lines (RILs), under NS and DS conditions, were 2363.49 and 897.38 kg ha−1, respectively. The outcome was a reduction in 62.03% of the grain yield in the RILs under water drought conditions. The grain yield estimates for the parents were 3381.50 kg ha−1 for IAPAR 81, and 2048.17 kg ha−1 for LP 97-28, under no-drought stress conditions (NS). Under DS conditions, the IAPAR 81 grain yield was 16.39% higher than LP 97-28 (IAPAR 81 GY = 980.23 kg ha−1 and LP 97-28 GY = 819.52 kg ha−1). The weights of 100 seeds of RILs were 17.79 g and 19.39 g for the DS and NS conditions, respectively.
The SW for the parents presented intermediate values of those shown by the RILs. The SW ranged from 16.77 to 19.70 g in the DS condition, and it varied from 17.84 to 25.68 g in the NS condition, demonstrating that water deficit enhanced its reduction. The RILs exhibited an average number of 11.41 pods per plant and 5.60 seeds per pod under the NS condition, while these values were 8.03 pods and 5.13 seeds under the DS condition. There were no observed differences in the number of days for flowering on the parents and RILs under the NS and DS conditions. However, genotypes that were cultivated under drought conditions (DS) showed a reduced interval for grain filling and maturity than those that were cultivated under non-drought conditions (NS).
In 2015, highly significant differences were observed between the study conditions DS and NS, for GY, NP, SW, PH, and NDF (Table 2). Moreover, genotype-by-environment condition interactions (G × E) exhibited significant differences for all the evaluated characteristics. We observed a decline in all the traits for both the experimental treatments during the agricultural year of 2015. GY for the RILs ranged from 744.72 to 1508.50 kg ha−1 under the NS conditions, whereas this trait varied from 521.96 to 1621.87 kg ha−1 under the DS conditions. Exposure to drought stress caused an average reduction in GY and SYD of 24% and 30%, respectively. SW varied from 19.60 to 23.78 in the NS treatment, while the interval for DS treatment ranged from 18.97 to 22.68 g. The mean NP was 8.61 (NS) and 7.23 (DS), whereas the average PH values for the NS and DS experiments were 50.49 and 45.65 cm, respectively. NDF, NDM, and NDE showed the highest mean values throughout the year compared to the previous experimental year (2014).
In both the years (2014 and 2015), the RILs showed higher values than the parents for the grain yield and yield components, which characterized the presence of transgressive segregation in this population (Figure 1). Additionally, the normal distribution of all the evaluated traits shows that genetic control is inherited quantitatively.
The phenotypic correlation coefficients for GY and other variables were higher in the NS experiment than in the DS experiment (Table 3). Positive correlations between SYD and the yield components have previously been reported under drought conditions [8,33,34,35]. The GY was positively correlated with SYD (r = 0.73), NP (r = 0.65), and SW (r = 0.42). The correlation between NP and SYD was positive and significant (r = 0.57). Moreover, the 100-seed weight and SYD revealed a positive and significant correlation (r = 0.35). Under the DS condition, the correlation was positive and significant between the yield components for NP and SYD (r = 0.72), and SW and SYD (r = 0.46). The correlations among NDE and NDM were positive and highly significant in the DS (r = 0.93) and NS (r = 0.84) experimental conditions. Conversely, the grain yield and its primary components and phenology traits showed negative correlations(Figure 2).

2.3. Genetic Mapping

This study genotyped 160 RILs from the IAPAR 81 × LP97-28 cross with 5.398 SNPs. A total of 1.287 SNPs were polymorphic between the parents. After excluding the markers that did not meet the minimum allele frequency of 30%, we selected a total of 773 SNP markers for mapping procedures. Mapping of the IAPAR81 × LP97-28 population covered 815.8 cM, with an average distance between the markers of 1.34 cM (Table 4). The number of markers per chromosome ranged from 23 (Pv04 and Pv11) to 195 in Pv02. The chromosome Pv04 had the shortest length, at 13.1 cM, while Pv01 had the longest length, at 117.9 cM.

2.4. QTLs Mapping

A total of 15 QTLs, associated with GY, SYD, SW, NP, PH, NDF, and NDM, were identified on Pv01, Pv02, Pv03, Pv07, Pv08, Pv09, Pv10, and Pv11 (Table 5, Figure 3). The QTLs that were associated with the GY (NS) conditions were mapped in Pv09, linked to the marker ss715640302, in the interval of 46.7 to 59.9 cM, and accounted for 8.7% of the phenotypic variance. SYD exhibited QTLs mapped on the chromosomes Pv07, Pv08, and Pv09. Additionally, QTLs in Pv07 in 2014, under drought stress, were linked to the ss715640487 marker, within a mapping interval from 49.8 to 57.1 cM. For the combined analysis (2014 and 2015), under drought stress, we identified QTLs that were mapped on Pv08 in 2014 and 2015, and linked to the marker ss715648231, with a peak position of 43.0 cM. The QTL in Pv09 was mapped between 13.5 and 27.3 cM to the ss715649357 marker, under NS conditions, in 2014. SYD suggests the presence of alleles that were originally from IAPAR81, under both the experimental conditions DS and NS. These QTLs demonstrated phenotypic variations of 12.1, 9.0, and 8.3%.
Figure 3 shows the order and location of the SNP marker in the linkage groups that agreed with the ranking and assignment of chromosomes on the common bean map from the Stampede/Red Hawk RIL population [36].

2.4.1. QTLs for 100-Seed Weight (SW)

Interestingly, we observed that the QTLs on Pv07 and Pv08, controlling SW in the population, were derived from the cross IAPAR81 × LP97-28. On Pv07, QTL is linked to the marker ss715647728 at 61 cM, for experiments DS and NS, and the combination of years.
In 2014, this QTL was identified under the DS condition linked to the SNP marker ss715639385. Likewise, the QTLs for SW on Pv08, which were inherited from the parental LP97-28, were found under NS and DS conditions. The most common SNP markers linked to this QTL were ss715648929, ss715648219, and ss715648220. The QTLs for SW showed a more constant behavior across the experimental years in comparison to other agricultural traits, such as SYD, NP, and PH (Table 5).

2.4.2. QTLs for Number of Pods per Plant (NP)

QTL mapping for the number of pods per plant showed stable performance across the experimental years. Beneficial alleles, derived from IAPAR81, were expressed under DS conditions, whereas negative alleles, donated from LP97-28, were expressed under DS conditions. One QTL was found in Pv11 under the DS condition, in both the experimental years. This QTL is linked to the marker ss715647465 at intervals of 31.1 and 41.4 cM, and contributes 13% of the phenotypic variation (Table 5).

2.4.3. QTLs for Plant Height (PH)

Three QTLs for PH were mapped in Pv01, Pv02, and Pv03, under DS conditions for each year and the combined analysis (both the experimental conditions). The QTL in Pv01 linked to ss715646868 marked at 67.0 cM for both the years and combined research under DS conditions. LP97-28 drought tolerance alleles contributed to PH, and accounted for 8.9% of the phenotypic variation. Two QTLs were found on Pv02, the first being in the NS condition, located in the interval of 18.5 to 29.8 cM, with a phenotypic variation of 8.7%, and is linked to the marker ss715640941. The second QTL was found in both the years and the combined analysis, mainly under drought conditions. The markers ss715646144 and ss715646929 were linked to the second QTL, explaining a phenotypic variation of approximately 9%. Only a single QTL for PH (under DS conditions) was identified in Pv03, through combined analysis across the years. Beneficial alleles (donated from IAPAR81) were linked to the marker ss715650580, and explained 9.5% of the phenotypic variation (Table 5).

2.4.4. QTLs for Number of Days to Flowering (NDF)

We detected a total of three QTLs flowering on the chromosomes Pv02, Pv03, and Pv10 for several days (Table 5). The QTL found in Pv02 in 2014, and in the combined analysis under NS conditions, explained more than 10% of the phenotypic variance, and was linked to the markers ss715648819, ss715647526, and ss715647652. In the experiment (NS) that was conducted in 2014, we found a QTL on Pv03 that was linked to the marker ss715647689, accounting for 12.3% of the phenotypic variance. The QTL (controlled by IAPAR81) was mapped in Pv10, linked to the marker ss71539783, under DS conditions in the combined analyses of the years and accounted for 8.6% of the phenotypic variation.

2.4.5. QTLs for Number of Days to Maturity (NDM)

Regarding the number of days to maturity, one QTL was identified in Pv07. In the experiment DS, 2014 and NS, 2015, and for the combined analysis, we identified QTLs that were linked to five markers at intervals starting from 50.6 to 74.5 cM.

3. Discussion

3.1. Genetic Diversity

At the molecular level, recent studies based on sequence data indicate that common beans are of Mesoamerican origin and originate in Mexico [12,14]. In Brazil, there is a predominance of Mesoamerican cultivars, which present a high level of genetic diversity when compared to the Andean cultivars [7].
In the present study, the RIL population derived from the cross IAPAR 81 × LP97-28 was used to map the QTLs that were related to drought tolerance in the common bean. The two parents IAPAR 81 and LP97-28 that were involved in the cross have Mesoamerican origin and belong to the carioca market class, which is the most consumed market class in Brazil. Both the parents contrasted in several traits related to drought tolerance, such as days to flowering and guiding, seed protein content, and grain yield. Most of the QTLs and markers associated with drought in the common bean have been based on crosses between the Mesoamerican and Andean genotypes, to facilitate abundant DNA polymorphisms [33]. According to the same authors, since most genetic breeding is carried out within the gene pools or market classes, there is still a need for markers that will permit selection among genetically similar materials, such as carioca beans.

3.2. Phenotypic Data

The common bean requires at least 100 mm of rainfall per month for appropriate development. In 2014, the precipitation in the area where the experiments were conducted was not sufficient for the optimal development of the bean plants. The rainfall was 200 mm throughout the entire crop cycle in 2014. On the other hand, in 2015, there was sufficient precipitation (599 mm), ensuring good crop development. However, the rainfall during the experiments conducted in 2015 was not regular, providing intervals that were characterized by no rain and high temperatures.
In both 2014 and 2015, the rain pattern during the experiment led to prolonged drought in the initial stages of plant development. Thus, the IAPAR 81 × LP97-28 RIL population was exposed to a range of stresses throughout the development stage (water deficit combined with high temperatures), representing the actual conditions in Brazilian producing areas, and contributing to identifying QTLs for drought tolerance in the common bean. Common bean plants are susceptible to water stress, which, especially when combined with high temperatures, produces effects such as abortion of flowers, pods, and grains that directly affect the grain yield [1].
In general, the IAPAR 81 × LP97-28 population had a maturation period that was longer in 2015 than in 2014. In 2014, the minimum temperature was 14.82 °C, and the maximum temperature was 31.34 °C. In 2015, the minimum and maximum temperatures were more pronounced, at 17.20 and 37 °C, respectively. Several authors have described that the common bean grows faster when the temperatures are elevated, due to flower, pod, and grain abortion, directly reflecting the low grain yield [8,16]. Thus, the average yield of grains under water deficit was higher in 2014 than in 2015.
The performance of the RILs was characterized by a transgressive genetic effect on the grain yield and the grain yield components that were inherited by the drought-tolerant parent used in the cross, that is, IAPAR 81 (Figure 2). We identified transgressive segregation for the grain yield in intragenic Mesoamerican crosses. Similar results were observed in the DOR364 × BAT 477 population [7] and Yuste-Lisbona et al. [37] in the intergenic population PMB0222 × PHA1037. Beebe et al. [1] reported that the transgressive segregation of RILs is due to the complementarity and/or mechanisms of the parents that allow the expression for drought tolerance. Diaz et al. [2] studied a multiparent advanced generation intercross (MAGIC) and noticed significant genetic variability within all the traits, and transgressive segregation was in most of them.
The average 100-seed weight that was observed in the experiment for the RILs was slightly lower than the 100-seed weight from the parent IAPAR 81 under drought stress and no drought stress conditions in both the years. According to other studies, the effect of water stress shows a slight reduction in the 100-seed weight of common beans [8,10,11,16]. Grain filling is inhibited under water stress, so RILs with greater than 100-seed weight are indicated as the most tolerant to drought [20]. When the DII (stress intensity) is light, the 100-seed weight presents itself as a stable characteristic over the years, which would make it worthwhile for the drought tolerance breeding of common beans [7]. Schneider et al. [27] demonstrated that the 100-seed weight in common beans is relatively stable, with the DII ranging from 0.19 to 0.49. The fact that the 100-seed weight correlated with the number of pods per plant under both the NS and DS stress conditions in 2014 and 2015 indicates that the genes that are responsible for the grain yield are linked to the genes that control the 100-seed weight, signifying the presence of pleiotropy. This claim is further supported by identifying QTLs that are related to the 100-seed weight and yield per day, co-locating in the same link group Pv07. This co-segregation of genes for the 100-seed weight and yield per day facilitates the simultaneous selection of lines with high grain yield in this Mesoamerican population.
Nonsignificant negative correlations between phenology and grain yield characteristics in the NS condition, and significant negative correlations in the DS condition, can be explained by the rainfall variation in the experimental locations. That is, days without rain in the initial flowering stage characterize an intermittent drought instead of the terminal drought that reaches the grain ripening phase.
The characteristics of grain yield per day, number of pods per plant, and 100-seed weight showed a correlation with grain yield for the non-drought stress condition in 2014 and 2015 (Table 3). On the other hand, no association between grain yield and yield-related traits was observed under the drought stress conditions. The number of days for flowering and maturation, and the days for filling, were negatively correlated with the grain yield components under drought stress, suggesting that maturation was necessary in order to escape drought in the two agricultural years. These correlation effects indicate that, depending on the intensity of the deficit for each year, several resistance mechanisms contributed to the complete performance of this population. These effects also demonstrate the complexity of breeding for drought tolerance in the common bean, mainly to address the local stress conditions and annual variation [8].

3.3. QTLs Mapping

In this study, the presence of transgressive segregation had an essential effect on the populations of Mesoamerican origin, and it was very noticeable for the grain yield and yield components (Figure 2). Similar results were observed by Diaz et al. [2] for DF, DM, SW, and productivity characteristics. We identified quantitative trait loci for several characteristics that were mapped on the same chromosome region, such as those on Pv07 and Pv09. This fact may be related to the cause–effect structure, which indicates that these characteristics may be conditioned by physically linked genes (pleiotropic) [8,24]. According to the criteria of Frei et al. [38], Miklas et al. [31], Assefa et al. [39], Beebe et al. [1], and Mukeshimana et al. [8], for studies that evaluate the performance of genotypes under water stress, identical genotypes must be assessed under normal irrigation conditions and under water restriction. In the present study, we met these conditions, showing that the genotypes’ performance grew in a stressed field condition, where water supplementation was restricted at the beginning of the flowering stage, and a stress-free field condition, with good water maintenance, as required for normal plant development.
The QTL GY9IL for grain yield was identified on chromosome Pv09, under NS conditions, in 2014. No association was found for grain yield under drought stress in Pv09, since the allele contribution was from the parent LP97-28 in the individual harvests and combined analyses. Several studies have described the presence of QTLs for grain yields in Pv09 [8,40]. An essential challenge for breeders working with certain classes of common beans has been the negative association between seed weight and productive potential in some crosses [8,40]. In the present study, positive correlations were observed in grain yield, concerning production components. In addition, the grain yield alleles per day (SYD7IL) and 100-seed weight (SW7IL) that were identified under the stress condition, were contributed by the parent LP97-28, consequently favoring the simultaneous selection of two characteristics in this cross between the Mesoamerican cultivars.
QTLs for the number of pods per plant trait, mapped in Pv11 (NP11IL), were identified under DS conditions, which indicates that this QTL has an indispensable effect on the number of pods per plant, and can be used to select beans with a higher number of pods per plant under drought-stressed conditions. Quantitative trait loci for pods per plant were previously mapped to Pv07, Pv09, and Pv11 in populations in an intergenic population with advanced backcrossing [41]. In studies conducted by Mukeshimana et al. [8], QTLs in Pv03 and Pv08 were identified. Moreover, Checa and Blair [42] mapped QTLs for the number of pods per plant on Pv04 and Pv10, although Beattie et al. [43] mapped QTLs for the number of pods per plant on Pv03. Although it is difficult to compare these QTLs, due to the different marker technologies that were used, there is strong evidence for the presence of QTLs affecting the number of pods per plant on Pv03.
Seed-weight QTLs were exclusively mapped in two genomic regions on Pv07 and Pv08. A highly expressive QTL for 100-seed weight, SW8IL, explained between 8.5 and 21% of the phenotypic variation. This QTL was found in 2015 and studied for two years, and combined for DS and NS conditions and under DS conditions in 2014. A second QTL, SW7IL, was consistently mapped to Pv07 under DS and NS conditions. This QTL was mapped in the same QTL regions for grain yield per day (SYD7IL) and the number of days for maturation (NDM7IL). The connection of different QTLs in Pv07 provides an option for simultaneous selection for 100-seed weight, and indirect selection for grain yield and maturation. These results agree with several previous studies that mapped QTLs to 100-seed weight in Pv07 [8,40,44]. The location of QTLs for 100-seed weight in Pv07 reinforces the link with the phaseolin protein locus in seeds (Phs), which was previously mentioned in several populations using different protein identification methodologies [42,44,45].
In this work, three (PH1IL, PH2IL, PH3IL) QTLs for plant height were observed on the chromosomes Pv01, Pv02, and Pv03, under DS, showing consistency between 2014 and 2015. Previous studies, conducted by Diaz et al. [46] and Asfaw et al. [24], described PHI QTLs on the chromosomes Pv06 and Pv08, and these chromosomal regions would be considered alternatives for breeding programs. Previous studies identified QTLs that were related to plant height, located on the Pv03 chromosomes, Pv10 and Pv11 [42,43,47], and in Pv6 [48].
An essential mechanism for adapting to water stress is phenological flexibility [8,22,40]. In the present study, QTLs for the number of days for flowering for the stressed condition (DS) were mapped on Pv10, corroborating the studies that were carried out by Kwak et al. [49] and Mukeshimana et al. [8]. In comparison, the QTLs for NDF under NS conditions were found in Pv02 and Pv03. The QTL for the number of maturation days located on Pv07, under DS and NS conditions, explained 8.4% to 12.1% of the phenotypic variation, respectively, constituting only the favorable alleles of the parental IAPAR81. Similar results were observed in several studies involving different populations of common beans under water stress [8,9,10,24,26,31].
DNA markers improve the efficiency and precision of plant breeding via assisted selection in the development of new common bean cultivars that are tolerant to water deficiency. Therefore, there is extensive variability in the parent under field conditions, making the selection of promising lines under field conditions a critical part of the breeding program. Using a high-resolution genetic map and phenotypic field data under different water deficit conditions, QTLs are linked to traits that guarantee the mapping of adaptation of common beans to water stress environments. Additionally, the availability of molecular markers that are linked to drought tolerance traits, which were identified through crosses between genotypes of the carioca commercial group (IAPAR81 × LP97-28), is of particular importance for common bean breeding programs in Brazil, as the commercial group carioca is the most produced and consumed.

4. Material and Methods

4.1. Population Development

In the present study, 160 recombinant inbred lines (RILs) derived from the cross between IAPAR 81 × LP 97-28 were evaluated. The F1 generation of the cross IAPAR 81 × LP 97-28 was obtained at the Instituto Agronômico do Paraná (IAPAR) (Londrina, PR, Brazil). The F2 to F6 generations were conducted through the single-seed descendant (SSD) method in a greenhouse at the Núcleo de Pesquisa Aplicada à Agricultura (Nupagri) of the Universidade Estadual de Maringá (UEM) (Maringá, PR, Brazil). The IAPAR 81 cultivar is known to be tolerant to water deficiency.
According to Moda-Cirino et al. [50], the IAPAR 81 cultivar was developed from multiple crossings (IAPAR 81 = (A248 × EM P 117) × {Araúna Selection × [BAT 93 × (Carioca Selection 99 × Great Northern Nebraska 1 Selection 27)]}) between lines considered resistant to disease and was released in July 1997 in Paraná state. This cultivar has an indeterminate growth habit, which allows mechanical harvesting. The average flowering time of IAPAR 81 was 43 days and 92 days from emergence to harvest. IAPAR 81 exhibits moderate tolerance to drought and moderate resistance to anthracnose, rust, powdery mildew, and mosaic virus. On the other hand, it is susceptible to bacterial blight and angular leaf spot. Line LP 97-28 was created based on multiple crossings (LP97-28 = [bz1977-2//IAPAR 14/Carioca 80 Selection]/[BAT93/2/Selection Carioca 99/Great Northern Nebraska/3//Araúna Selection] also conducted by IAPAR, but it has not yet been released. It belongs to the commercial group carioca, upright bush type, indeterminate growth habit (type II), brown-striped cream grains, and average protein content of 27.65%. It shows low tolerance to water deficit [51].

4.2. Field Trial Design, Evaluation of Productivity and Morpho Agronomic Characteristics

We evaluated the parents (IAPAR 81 and LP 97-28), 160 recombinant inbred lines from this cross, and seven controls (SCS Guará, Flor Diniz, IPR Juriti, BRS Pérola, IPR Tangará, BRS Talismã, and BAT93) in field conditions during two agricultural years (2014 and 2015) and two environmental conditions, drought stress (DS) and no drought stress (NS). Using this methodology, we could evaluate the performance of the genotypes in a situation with water supply according to the crop requirements (NS). After the flowering period, the water supply was interrupted, simulating drought conditions for ten days (DS).
Field experiments to evaluate the 160 recombinant inbred lines from IAPAR 81 × LP 97-28 cross to drought tolerance were conducted at the Centro Técnico de Irrigação (CTI) from Universidade Estadual de Maringá (UEM) (Maringá, PR, Brazil). The experiment was specifically located at latitude 23°25′ S and longitude 51°57′. The climate is characterized as warm temperate. The annual average rainfall for the local area is estimated to be 1276 mm, and the annual average temperature is 22 °C at an elevation of 515 m above sea level. The soil is classified as dystroferric red latosol (oxisol) with a pH value ranging from 5.7 to 6.5 and clayey soil texture. Soil contents included sand (122.6 g kg−1), silt (120.6 g kg−1) and clay (756.8 g kg−1).
In 2014, plots were composed of two 1.5 m rows spaced 0.5 m apart, with 1 m between plots. The design experiment was a triple 13 × 13 lattice experiment with two repetitions for each DS and NS experiment in 2014. In 2015, the plot design was 2 m rows spaced 0.5 m apart and 1 m between plots with three repetitions for each DS and NS experiment. The design experiment was a triple 13 × 13 lattice experiment with three repetitions for each DS and NS experiment in 2015.
For both experiments, a 15 m border of cultivar IAPAR 81 avoided lateral water movements through the soil between the DS and NS experimental conditions. All experiments were carried out at the end of July to increase the chances of extending the drought season during flowering and grain filling periods. Harvesting in both agricultural years occurred at the beginning of November.
Every three days, a water supply was established for the two environmental conditions (DS and NS) from the sowing data until the 20th day after emergence, allowing soil humidity to achieve field capacity. After the 20th day, irrigation was implemented to support soil water tension between 6 KPa and 50 KPa. The 6 kPa tension corresponds to field capacity humidity, while 50 kPa was applied as a low soil water tension limit recommended for common bean crops. Soil humidity was measured with 12 tensiometers installed in each DS and NS experiment.
At the phenological stage of flowering (R5), two irrigation management approaches were used for DS and NS conditions. For NS (no drought), irrigation was continued until the end of the crop cycle, keeping the soil water tension between 6 KPa and 50 KPa. For the DS (with drought) condition, irrigation was discontinued after the flowering period (R5) for 10 days. The period of ten days of stress for drought began when the soil water tension was 50 KPa. After a 10-day drought stress simulation, the water supply was re-established using the same strategy as the NS experiment. Environmental data, such as rainfall (mm), temperatures (°C) and soil water tension (KPa), were collected by the meteorological station located at Centro Técnico de Irrigação (CTI).
The following agronomic traits were collected: (a) number of days to flowering (NDF)—days between emergence until flowering, when 50% of the plants from the plots exhibited at least one flower; (b) number of days to maturity (NDM)—days between emergence until 90% of the pods lost their green coloration and became too dry; (c) number of days to grain filling (NDG)—calculation of the days between flowering and maturation; (d) grain yield (GY)—a measurement based on each plot at 14% humidity and calculated in kg ha−1; (e) seed yield per day (SYD)—refers to the total weight of seeds per plot/number of days to maturity and is estimated in g.m−2; (f) 100-seeds weight (SW)—the weight of a plot sample containing 100 seeds expressed in grams; (g) number of pods per plant (NP)—total number of pods per plot divided by the number of plants per plot; (h) number of seeds per pod (NSP)—number of seeds of 10 pods randomly selected within a plot; (i) plant height (PH)—the distance between soil and tip of leaves measured in cm; (j) drought intensity index (DII)—assessed for each experiment as 1-(Xd/Xp), where Xd corresponds to the mean yield of genotypes under DS condition and Xp is the mean yield of genotypes under ND condition.

4.3. Phenotypic Data Analysis

The phenotypic data of each trial were analyzed by fitting SELEGEN software [52]. The prediction of genotypic values was performed via BLUP (best linear unbiased predictor) using estimates of variance components of random factors obtained by the restricted maximum likelihood (REML) method. Mixed models were used to analyze the phenotypic data for one location and one harvest season, and the following mathematical model for an incomplete block design was adopted (Model 17 Selegen):
y = X r + Z g + W b + e
y = data vector; r = replication effect vector (assumed as fixed) added to the overall mean; g = genotypic effect vector (assumed as random); b = block effect vector (assumed as random); and e = error or residuals value (random). The uppercase letters represent the incidence matrices for the aforementioned effects.
The significance of the effects of the models was estimated using the likelihood ratio test (LRT), thus configuring a deviance analysis (ANADEV). The goodness of fit was obtained based on studies with and without genotype effects and genotypes × environments interaction. Then, the deviance without that effect was subtracted from each deviance of the complete model, comparing it with the chi-square value (x2) with 1 degree of freedom at 1% and 5% probability [53].
Mathematically, as follows:
L R T = 2 l n ( M L   o f   r e d u c e d   m o d e l M L V   o f   c o m p l e t e   m o d e l   )
ln represents the Naperian logarithm, and ML is the maximum likelihood estimation.
The analyses fit the effect of each line on the phenotype as fixed and random to obtain the best linear unbiased estimators (BLUEs) and best linear unbiased predictors (BLUPs), respectively. The best linear unbiased estimates (BLUEs) obtained through the REML/BLUP analysis were used to calculate Pearson’s correlation coefficients. Pearson’s correlation between grain yield (GY), seed yield per day (SYD), number of pods per plant (NP), 100-seed weight (SW), plant height (PH), number of days for flowering (NDF), number of days to mature (NDM) and the number of days to grain filling (NDE) were calculated independently for each of the environmental conditions (NS—no drought stress; DS—with drought stress) evaluated in two agricultural years. Significant differences were assumed at 0.05, 0.01 and 0.001 probability levels.

4.4. Genotyping of the RILs Using BARCBean6K_3 BeadChip

Total genomic DNA was isolated from the 160 RIL families (F8 generation) and parents (IAPAR 81 × LP 97-28) using the DNeasy plant mini kit (Qiagen, CA, USA) following the manufacturer’s instructions. The DNA was quantified using a 1.5% agarose gel (Agarose SFR, Amresco, IL, USA) with TBE buffer (tris-borate-ethylenediamine tetraacetic acid) and stained with 1 μg mL−1 Sybr Green (Sigma-Aldrich, St. Louis, MI, USA).
The DNA samples were screened with 5398 SNP DNA markers on the BARCBean6K_3 Illumina BeadChip by following the Infinium HD assay ultra protocol (Illumina, Inc., San Diego, CA, USA). The BeadChip was imaged using the Illumina beadarray reader to measure fluorescence intensity. Automatic allele calling for each locus was performed with Genome Studio Genotyping Module v1.8.4 software (Illumina, San Diego, CA, USA) and visually inspected. Any errors in allele calling due to improper cluster identification were corrected, resulting in 4633 high-quality SNPs.
Subsequently, after checking for the quality of SNPs with the software GenomeStudio, we used an electronic spreadsheet in Microsoft Excel to perform the filtering process by excluding monomorphic SNPs among the parents IAPAR 81 (tolerant parent) and LP 97-28 (not tolerant parent), with lost data greater than 5%, and with a minimum allele frequency minimum allele frequency less than 30%. In addition, SNP markers were placed on their respective mapping chromosomes and assigned their nomenclature deposited at the National Center for Biotechnology Information advances (NCBI).

4.5. Genetic Mapping

SNP data from the 160 recombinant inbred lines from IAPAR 81 × LP 97-28 genotyped with the BARCBean6K_3 BeadChip and previously filtered were used to construct genetic linkage maps for P. vulgaris. Only SNP markers polymorphic between parents IAPAR 81 and LP 97-28 that segregated in the RIL population and exhibited a minimum allelic frequency of 30% were used for the linkage map construction. The genetic distances (cM) for each of the eleven linkage groups of P. vulgaris were calculated using JoinMap 4.1 software [54]. Default settings for the regression mapping algorithm based on the Kosambi map function were used to define linkage order and distances (cM). The genetic distances (cM) were later used to graphically design the linkage groups using MapChat 2.3 [55].

4.6. QTL Analysis

Quantitative trait loci (QTL) analysis was performed using the new adjusted genotypic values of the evaluated characteristics obtained through the REML/BLUP analysis for the years 2014 and 2015 and the combination of agricultural years for each of the environmental conditions (NS—no drought stress; DS—with drought stress) evaluated. QTL analysis was conducted using Win Cartographer V2.5-011 [56]. For this study, we used the composite interval mapping method (CIM), model 6. We set the parameters at a window size of 10 cM and a walking speed of 1 cM. The QTL was considered significant based on the 1000-permutation test for each trait in Win QTL Cartographer to determine the LOD threshold at p = 0.05, and the peak LOD score was considered the location of the QTL. The proportion of phenotypic variance accounted for by the QTL was determined using the R2 value. The QTLs were named according to the common bean QTL nomenclature guidelines proposed by Miklas and Porch [57]. Finally, using MapChat 2.3 [55], the QTLs were graphically represented in the genetic linkage maps for P. vulgaris developed using the RIL population IAPAR 81 × LP 97-28.

5. Conclusions

The evaluation of the RIL population from the IAPAR81 × LP97-28 cross, under two different conditions of water supplementation, no drought stress (NS) and drought stress (DS), showed a reduction in grain yield, yield per day, and the number of pods per plant when plants were submitted to an environment with a reduction in water supplementation. However, the lines derived from the IAPAR81 × LP97-28 cross have shown potential for overcoming water deficiency during the flowering stages, which is one of the most critical phases for cultivating common beans. Using 773 polymorphic SNP markers in the RIL population, and phenotypic data obtained in the DS and NS conditions, we identified consistent QTLs associated with six traits that ensure adaptation of the RILs to water stress in common beans. Thus, a total of 15 QTLs associated with grain yield, seed yield per day, 100-seed weight, number of pods per plant, plant height, number of days for flowering, and number of days for maturation, were identified on the chromosomes Pv01, Pv02, Pv03, Pv07, Pv08, Pv09, Pv10, and Pv11.

Author Contributions

Design and supervision of the research—J.C.F.E., P.S.V.F., M.C.G.-V., V.M.-C., Q.S.; conduct of the experiments—J.C.F.E., Q.S.; analyzed the data—J.C.F.E., T.A.S.G.; G.V.; M.C.G.-V., M.V.B.; original draft of the manuscript—J.C.F.E., G.V., T.A.S.G.; revision and editing of the manuscript: M.C.G.-V., M.V.B., P.S.V.F. All authors have read and agreed to the published version of the manuscript.

Funding

Federal Funding Institution National Council for Scientific and Technological Development (CNPq), Grant/Award Numbers: 311599/2017-6, 408472/2018-9, 301861/2017-0; and from Coordination for the Improvement of Higher Education Personnel (Capes), Award number 88887.470337/2019-00.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data are presented within the article.

Acknowledgments

J.C.F.E. is grateful for the grant from CNPq (Brasilia, Brazil) and for the support and assistance provided by professors and students of the Núcleo de Pesquisa Aplicada à Agricultura (Nupagri) and Centro Técnico de Irrigação (CTI) of the Universidade Estadual de Maringá (UEM). M.C.G.-V. is thankful for a grant from CNPq. M.V.B. is supported by a scholarship from CAPES (Brasilia, Brazil). We also appreciate the critical reading and writing of the original draft in English of the manuscript by Pedrina Gonçalves Vidigal.

Conflicts of Interest

The authors declare that there are no conflict of interest.

References

  1. Beebe, S.E.; Rao, I.M.; Blair, M.W.; Acosta-Gallegos, J.A. Phenotyping common beans for adaptation to drought. Front. Physiol. 2013, 4, 35. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Diaz, S.; Ariza-Suarez, D.; Izquierdo, P.; Lobaton, J.D.; de la Hoz, J.F.; Acevedo, F.; Duitama, J.; Guerrero, A.F.; Cajiao, C.; Mayor, V.; et al. Genetic mapping for agronomic traits in a MAGIC population of common bean (Phaseolus vulgaris L.) under drought conditions. BMC Genom. 2020, 21, 799. [Google Scholar] [CrossRef]
  3. Gepts, P. Tropical environments, biodiversity, and the origin of crops. In Genomics of Tropical Crop Plants. Plant Genetics and Genomics: Crops and Models; Springer: New York, NY, USA, 2008; pp. 1–20. [Google Scholar] [CrossRef]
  4. Barili, L.D.; Vale, N.M.; Moura, L.M.; Paula, R.G.; Silva, F.F.; Carneiro, J.E. Genetic progress resulting from forty-three years of breeding of the carioca common bean in Brazil. Gen. Mol. Res. 2016, 15, 1–11. [Google Scholar] [CrossRef] [PubMed]
  5. McClean, P.E.; Burridge, J.; Beebe, S.; Rao, I.M.; Porch, T.G. Crop improvement in the era of climate change: An integrated, multi-disciplinary approach for common bean (Phaseolus vulgaris). Funct. Plant Biol. 2011, 38, 927–933. [Google Scholar] [CrossRef] [PubMed]
  6. Beebe, S. Common Bean Breeding in the Tropics. Plant Breed. Rev. 2012, 36, 357–426. [Google Scholar]
  7. Blair, M.W.; Galeano, C.H.; Tovar, E.; Torres, M.C.M.; Castrillón, A.L.; Beebe, S.E.; Rao, I.M. Development of a Mesoamerican intra-gene pool genetic map for quantitative trait loci detection in a drought tolerant × susceptible common bean (Phaseolus vulgaris L.) cross. Mol. Breed. 2012, 29, 71–88. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Mukeshimana, G.; Butare, L.; Cregan, P.B.; Blair, M.W.; Kelly, J.D. Quantitative trait loci associated with drought tolerance in common bean. Crops Sci. 2014, 54, 923–938. [Google Scholar] [CrossRef] [Green Version]
  9. Traub, J.; Kelly, J.D.; Loescher, W. Early Metabolic and Photosynthetic Responses to Drought Stress in Common and Tepary Bean. Crops Sci. 2017, 57, 1670–1686. [Google Scholar] [CrossRef]
  10. Trapp, J.J.; Urrea, C.A.; Zhou, J.; Khotc, L.R.; Sankaranc, S.; Miklas, P.N. Selective phenotyping traits related to multiple stress and drought response in dry bean. Crops Sci. 2016, 56, 1460–1472. [Google Scholar] [CrossRef] [Green Version]
  11. Blair, M.W.; Brondani, R.V.P.; Díaz, L.M.; Del Peloso, M.J. Diversity and population structure of common bean from Brazil. Crops Sci. 2013, 53, 1983–1993. [Google Scholar] [CrossRef]
  12. Bitocchi, E.; Rau, D.; Bellucci, E.; Rodriguez, M.; Murgia, M.L.; Gioia, T.; Santo, D.; Nanni, L.; Attene, G.; Papa, R. Beans (Phaseolus ssp) as a model for understanding crop evolution. Front. Plant Sci. 2017, 8, 722. [Google Scholar] [CrossRef] [Green Version]
  13. Possobom, M.T.D.F.; Ribeiro, N.D.; da Domingues, L.S.; Casagrande, C.R. Genetic control of iron concentration in Mesoamerican and Andean common bean seeds. Pesq. Agrop. Bras. 2015, 383–391. [Google Scholar] [CrossRef] [Green Version]
  14. Schmutz, J.; McClean, P.E.; Mamidi, S.; Wu, G.A.; Cannon, S.B.; Grimwood, J.; Jenkins, J.; Shu, S.; Song, Q.; Chavarro, C.; et al. A reference genome for common bean and genome-wide analysis of dual domestications. Nat. Genet. 2014, 46, 707–713. [Google Scholar] [CrossRef] [Green Version]
  15. Frahm, M.A.; Rosas, J.C.; Mayek-Pérez, N.; López-Salinas, E.; Acosta-Gallegos, J.A. Breeding beans for resistance to terminal drought in the lowland tropics. Euphytica 2004, 136, 223–232. [Google Scholar] [CrossRef]
  16. Miklas, P.N.; Singh, S.P. Common bean. In Genome Mapping and Molecular Breeding in Plants: Pulses, Sugar and Tuber Crops; Kole, C., Ed.; Springer: Berlin/Heidelberg, Germany, 2007; pp. 1–31. [Google Scholar]
  17. Singh, V.; Nguyen, C.T.; Yang, Z.; Chapmand, S.C.; van Oosteroma, E.J.; Hammera, G.L. Genotypic differences in effects of short episodes of high-temperature stress during reproductive development in sorghum. Crops Sci. 2016, 56, 1561–1572. [Google Scholar] [CrossRef]
  18. Porch, T.; Beaver, J.; Debouck, D.; Jackson, S.A.; Kelly, J.D.; Dempewolf, H. Use of wild relatives and closely related species to adapt common bean to climate change. Agronomy 2013, 3, 433–461. [Google Scholar] [CrossRef] [Green Version]
  19. Oliveira, M.R.C.; Corrêa, A.S.; de Souza, G.A.; Guedes, R.N.C.; Oliveira, L.O. Mesoamerican origin and pre- and post-columbian expansions of the ranges of Acanthoscelides obtectus Say, a cosmopolitan insect pest of the common bean. PLoS ONE 2013, 8, e70039. [Google Scholar] [CrossRef] [Green Version]
  20. Ramirez-Vallejo, P.; Kelly, J.D. Traits related to drought resistance in common bean. Euphytica 1998, 99, 127–136. [Google Scholar] [CrossRef]
  21. Rosales-Serna, R.; Kohashi-Shibata, J.; Acosta-Gallegos, J.A.; Trejo-Lópeza, C.; Ortiz-Cereceres, J.; Kelly, J.D. Biomass distribution, maturity acceleration and yield in drought-stressed common bean cultivars. Field Crop. Res. 2004, 85, 203–211. [Google Scholar] [CrossRef]
  22. Asfaw, A.; Ambachew, D.; Shah, T.; Blair, M.W. Trait associations in diversity panels of the two common bean (Phaseolus vulgaris L.) gene pools grown under well-watered and water-stress conditions. Front. Plant Sci. 2017, 8, 1–15. [Google Scholar] [CrossRef] [PubMed]
  23. Kavar, T.; Maras, M.; Kidrič, M.; Suštar-Vozlič, J.; Meglič, V. Identification of genes involved in the response of leaves of Phaseolus vulgaris to drought stress. Mol. Breed. 2008, 21, 159–172. [Google Scholar] [CrossRef]
  24. Asfaw, A.; Blair, M.W. Quantitative trait loci for rooting pattern traits of common beans grown under drought stress versus non-stress conditions. Mol. Breed. 2012, 30, 681–695. [Google Scholar] [CrossRef]
  25. Toaldo, D.; de Morais, P.P.P.; Battilana, J.; Coimbra, J.L.M.; Guidolin, A.F. Selection in early generations and the occurrence of heterosis for the character root distribution. Euphytica 2013, 190, 335–344. [Google Scholar] [CrossRef]
  26. Hoyos-Villegas, V.; Song, Q.; Wright, E.M.; Beebe, S.E.; Kelly, J.D. Joint linkage QTL mapping for yield and agronomic traits in a composite map of three common bean RIL populations. Crops Sci. 2016, 56, 2546–2563. [Google Scholar] [CrossRef]
  27. Schneider, K.A.; Rosales-Serna, R.; Ibarra-Perez, F.; Cazares-Enriquez, B.; Acosta-Gallegos, J.A.; Ramirez-Vallejo, P.; Wassimi, N.; Kelly, J.D. Improving common bean performance under drought stress. Crops Sci. 1997, 37, 43–50. [Google Scholar] [CrossRef]
  28. De Ron, A.M.; Kalavacharla, V.K.; Álvarez-García, S.; Casquero, P.A.; Rodríguez-González, A.; Suárez-Villanueva, S.; Rodiño, A.P.; Beaver, J.S.; Porch, T.; Galván, M.Z.; et al. Common Bean Genetics, Breeding, and Genomics for Adaptation to Changing to New Agri-environmental Conditions. In Genomic Designing of Climate-Smart Pulse Crops; Kole, C., Ed.; Springer: New York, NY, USA, 2019; pp. 1–106. [Google Scholar] [CrossRef]
  29. Blum, A. Drought resistance-is it a complex trait? Funct. Plant Biol. 2011, 38, 753–757. [Google Scholar] [CrossRef]
  30. Muñoz-Perea, C.G.; Terán, H.; Allen, R.G.; Wright, J.L.; Westermann, D.T.; Singh, S.P. Selection for drought resistance in dry bean landraces and cultivars. Crops Sci. 2006, 46, 2111–2120. [Google Scholar] [CrossRef] [Green Version]
  31. Miklas, P.N.; Kelly, J.D.; Beebe, S.E.; Blair, M.W. Common bean breeding for resistance against biotic and abiotic stresses: From classical to MAS breeding. Euphytica 2006, 147, 105–131. [Google Scholar] [CrossRef]
  32. Dramadri, I.O.; Nkalubo, S.T.; Kelly, J.D. Identification of QTL Associated with Drought Tolerance in Andean Common Bean. Crops Sci. 2019, 59, 1007–1020. [Google Scholar] [CrossRef]
  33. Beebe, S.; Ramirez, J.; Jarvis, A.; Rao, I.M.; Mosquera, G.; Bueno, J.M.; Blair, M.W. Genetic Improvement of Common Beans and the Challenges of Climate Change. In Crop Adaptation to Climate Changes; John Wiley & Sons, Ltd.: Chichester, UK, 2011; pp. 356–369. [Google Scholar]
  34. Polania, J.A.; Rao, I.M.; Cajiao, C.; Rivera, M.; Raatz, B.; Beebe, S.E. Physiological traits associated with drought resistance in Andean and Mesoamerican genotypes of common bean (Phaseolus vulgaris L.). Euphytica 2016, 210, 17–29. [Google Scholar] [CrossRef]
  35. Valdisser, P.A.M.R.; Müller, B.S.F.; Almeida Filho, J.E.; Morais Júnior, O.P.; Guimarães, C.M.; Borba, T.C.O.; Souza, I.P.; Zuc-chi, M.I.; Neves, L.G.; Coelho, A.S.G.; et al. Genome-wide association studies detect multiple QTLs for productivity in Mesoamerican Diversity Panel of common bean under drought stress. Front. Plant Sci. 2020, 11, 574674. [Google Scholar] [CrossRef]
  36. Song, Q.; Jia, G.; Hyten, D.L.; Jenkins, J.; Hwang, E.Y.; Schroeder, S.G.; Cregan, P.B. SNP assay development for linkage map construction, anchoring whole genome sequence and other genetic and genomic applications in common bean. G3 GenesGenomesGenet. 2015, 5, 2285–2290. [Google Scholar] [CrossRef] [Green Version]
  37. Yuste-Lisbona, F.J.; Santalla, M.; Capel, C.; García-Alcázar, M.; Fuente, M.D.L.; Capel, J.; Ron, A.M.D.; Lozano, R. Marker-based linkage map of Andean common bean (Phaseolus vulgaris L.) and mapping of QTLs underlying popping ability traits. BMC Plant. Biol. 2012, 12, 136. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Frei, A.; Blair, M.W.; Cardona, C.; Beebe, S.E.; Gu, H.; Dorn, S. QTL mapping of resistance to Thrips palmi karny in common bean. Crops Sci. 2005, 45, 379–387. [Google Scholar] [CrossRef]
  39. Assefa, T.; Beebe, S.E.; Rao, I.M.; Cuasquer, J.B.; Duque, M.C.; Rivera, M.; Battisti, A.; Lucchinc, L. Pod harvest index as a selection criterion to improve drought resistance in white pea bean. Field Crop. Res. 2013, 148, 24–33. [Google Scholar] [CrossRef]
  40. Beaver, J.S.; Osorno, J.M. Achievements and limitations of contemporary common bean breeding using conventional and molecular approaches. Euphytica 2009, 168, 145–175. [Google Scholar] [CrossRef]
  41. Blair, M.W.; Iriarte, G.; Beebe, S. QTL analysis of yield traits in an advanced backcross population derived from a cultivated Andean x wild common bean (Phaseolus vulgaris L.) cross. Theor. Appl. Genet. 2006, 112, 1149–1163. [Google Scholar] [CrossRef]
  42. Checa, O.E.; Blair, M.W. Mapping QTL for climbing ability and component traits in common bean (Phaseolus vulgaris L.). Mol. Breed. 2008, 22, 201–215. [Google Scholar] [CrossRef]
  43. Beattie, A.D.; Larsen, J.; Michaels, T.E.; Pauls, K.P. Mapping quantitative trait loci for a common bean (Phaseolus vulgaris L.) ideotype. Genome Natl. Res. Counc. Can. 2003, 46, 411–422. [Google Scholar] [CrossRef]
  44. Nodari, R.O.; Tsail, S.M.; Gilbertson, R.L.; Gepts, P. Towards an integrated linkage map of common bean 2. Development of an RFLP-based linkage map. Theor. Appl. Genet. 1993, 85, 513–520. [Google Scholar] [CrossRef]
  45. Koinage, E.M.K.; Singh, S.P.; Gepts, P. Genetic control of the domestication syndrome in common bean. Crops Sci. 1996, 36, 1037–1045. [Google Scholar] [CrossRef] [Green Version]
  46. Diaz, L.M.; Ricaurte, J.; Tovar, E.; Cajiao, C.; Teran, H.; Grajales, M.; Polania, J.; Rao, I.; Beebe, S.; Raatz, B. QTL analyses for tolerance to abiotic stresses in a common bean (Phaseolus vulgaris L.) population. PLoS ONE 2018, 13, e0202342. [Google Scholar] [CrossRef]
  47. Wright, E.M.; Kelly, J.D. Mapping QTL for seed yield and canning quality following processing of black bean (Phaseolus vulgaris L.). Euphytica 2011, 179, 471–484. [Google Scholar] [CrossRef]
  48. Tar’an, B.; Michaels, T.E.; Pauls, K.P. Genetic mapping of agronomic traits in common bean. Crops Sci. 2002, 42, 544–556. [Google Scholar] [CrossRef]
  49. Kwak, M.; Velasco, D.; Gepts, P. Mapping homologous sequences for determinacy and photoperiod sensitivity in common bean (Phaseolus vulgaris). Jour. Hered. 2008, 99, 283–291. [Google Scholar] [CrossRef]
  50. Moda-Cirino, V.; Oliari, L.; Lollato, M.A.; Fonseca, S. Iapar 81–Common bean. Crop. Breed. Appl. Biotechnol. 2001, 203–204. [Google Scholar] [CrossRef]
  51. Molina, J.; Moda-Cirino, V.; Júnior, N.S.F.; de Faria, R.T.; Destro, D. Response of common bean cultivars and lines to water stress. Crop. Breed. Appl. Biot. 2001, 363–372. [Google Scholar] [CrossRef]
  52. Resende, M.D.V. Selegen Reml/Blup–Sistema estatístico e seleção genética computadorizada. In Manual Complementar do Selegen-Reml/Blup; Embrapa Florestas: Colombo, Sri Lanka, 2014; pp. 1–17. [Google Scholar]
  53. Resende, M.D.V.; Duarte, J.B. Precisão e controle de qualidade em experimentos de avaliação de cultivares. Pesq. Agrop. Trop. 2007, 37, 182–194. [Google Scholar]
  54. Van Ooijen, J.W. JoinMap®4, Software for the Calculation of Genetic Linkage Maps in Experimental Populations; Kyazma, B.V., Ed.; Wageningen University: Wageningen, The Netherlands, 2006. [Google Scholar]
  55. Voorrips, R.E. MapChart: Software for the graphical presentation of linkage maps and QTLs. The J. Hered. 2002, 93, 77–78. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Wang, S.; Basten, C.J.; Zeng, Z.B. Windows QTL Cartographer 2.5. J. Infect. Dis. 2011, 204, 198–199. Available online: http://statgen.ncsu.edu/qtlcart/WQTLCart.htm (accessed on 1 June 2016).
  57. Miklas, P.N.; Porch, T. Guidelines for common bean QTL nomenclature. Annu. Rep. Bean Improv. Coop. 2010, 53, 202–204. [Google Scholar]
Figure 1. Soil water tension for drought stress (DS) and no-drought stress (NS) conditions and precipitation (mm) in the experiments conducted at the Technical Irrigation Center (CTI) in Maringá, PR, Brazil in 2014 and 2015. Dotted lines indicate the field capacity and wilt point according to the soil water tension values.
Figure 1. Soil water tension for drought stress (DS) and no-drought stress (NS) conditions and precipitation (mm) in the experiments conducted at the Technical Irrigation Center (CTI) in Maringá, PR, Brazil in 2014 and 2015. Dotted lines indicate the field capacity and wilt point according to the soil water tension values.
Plants 10 01568 g001
Figure 2. Distributions of the frequencies for the traits evaluated in the field: GY = grain yield; SYD = seed yield per day; NP = number of pods per plant; SW = 100-seed weight; PH = plant height; NDF = number of days for flowering; NDM = number of days to mature; and NDG = number of days to grain filling for the population of recombinant lines (RIL) from IAPAR 81 × LP97-28 cross under drought stress (DS) and no-drought stress (NS) conditions in Maringá, PR, Brazil in the years 2014 and 2015.
Figure 2. Distributions of the frequencies for the traits evaluated in the field: GY = grain yield; SYD = seed yield per day; NP = number of pods per plant; SW = 100-seed weight; PH = plant height; NDF = number of days for flowering; NDM = number of days to mature; and NDG = number of days to grain filling for the population of recombinant lines (RIL) from IAPAR 81 × LP97-28 cross under drought stress (DS) and no-drought stress (NS) conditions in Maringá, PR, Brazil in the years 2014 and 2015.
Plants 10 01568 g002
Figure 3. Genetic map for the common bean RIL population from IAPAR 81 × LP97-28 cross using 773 SNPs markers assigned to the 11 common bean chromosomes. QTLs mapped for each trait are colored as follows: GY = grain yield in dark blue; SYD = seed yield per day in green; NP = number of pods per plant in brown; SW = 100-seed weight in red; PH = plant height in black; NDF = number of days for flowering in pink; NDM = number of days to mature in light blue.
Figure 3. Genetic map for the common bean RIL population from IAPAR 81 × LP97-28 cross using 773 SNPs markers assigned to the 11 common bean chromosomes. QTLs mapped for each trait are colored as follows: GY = grain yield in dark blue; SYD = seed yield per day in green; NP = number of pods per plant in brown; SW = 100-seed weight in red; PH = plant height in black; NDF = number of days for flowering in pink; NDM = number of days to mature in light blue.
Plants 10 01568 g003aPlants 10 01568 g003b
Table 1. Drought intensity indices for IAPAR 81, LP97-28 and 160 recombinant inbred lines from this cross and controls for the agricultural years 2014 and 2015 in Maringá, PR, Brazil.
Table 1. Drought intensity indices for IAPAR 81, LP97-28 and 160 recombinant inbred lines from this cross and controls for the agricultural years 2014 and 2015 in Maringá, PR, Brazil.
TraitsDrought Intensity Indices
20142015
x _ g x _ L x _ P x _ P 1 x _ P 2 x _ L + P x _ g x _ L x _ P x _ P 1 x _ P 2 x _ L + P
GY0.590.590.700.720.650.650.300.290.380.210.500.34
SYD0.410.410.440.400.500.420.350.340.400.260.500.37
NP0.300.300.590.640.440.460.220.210.120.040.180.17
SW0.060.060.160.27-0.110.030.030.130.110.150.08
PH0.090.090.110.150.070.100.100.100.060.15-0.08
NDF--0.01-0.02-0.010.01--0.02-
NDM0.030.030.040.020.060.04----0.01-
NDG0.100.100.090.050.130.09------
x _ g : general average; x _ L : average of the 160 RILs; x _ P : average of the parents; x _ P 1 : average of IAPAR 81; x _ P 2 : average of LP97-28; x _ L + P : averages of RILs and parents. GY = grain yield; SYD = seed yield per day; NP = number of pods per plant; SW = 100-seed weight; PH = plant height; NDF = number of days for flowering; NDM = number of days to mature; and NDG = number of days to grain filling.
Table 2. Traits evaluated in the parents IAPAR 81 and LP97-28, and in 160 RILs derived from this cross for the experiments performed under drought stress (DS) and no-drought stress condition (NS).
Table 2. Traits evaluated in the parents IAPAR 81 and LP97-28, and in 160 RILs derived from this cross for the experiments performed under drought stress (DS) and no-drought stress condition (NS).
TraitsYearCon.ParentsRILs§PEV
IAPAR 81LP97-28AverageIntervalLRT (x2)
Yield and Its Components
Grain Yield
(kg ha−1)
2014NS3381.502048.172363.492020.68–3581.05169.22 ****31,674.07
DS980.23819.52897.38820.89–1068.3245.05 **5134.52
2015NS959.611207.77993.92744.72–1508.50297.85 ****3147.65
DS765.77751.01754.47521.96–1621.87319.35 **2369.41
Seed Yield per day (g m2)2014NS3.644.393.963.64–4.7418.73 ****0.079
DS2.632.152.372.155–2.9153.28 **0.037
2015NS2.202.742.271.71–3.55ns**0.017
DS1.651.631.601.11–3.45ns0.013
Number of pods per plant2014NS17.4910.2311.4110.12–18.13125.40 ****0.613
DS7.957.058.037.07–10.40214.06 **0.180
2015NS8.109.468.617.13–14.11118.32 ****0.464
DS7.448.887.235.59–12.67179.12 **0.323
Number of seeds per pod2014NS5.585.595.605.53–5.76nsns0.020
DS5.115.135.135.08–5.24ns0.019
2015NS5.105.115.004.91–5.17ns*0.040
DS4.584.554.514.36–4.80ns0.075
100-seed weight (g)2014NS25.5317.8119.3917.84–25.68374.43 ****0.268
DS18.9716.7417.7916.77–19.70150.26 **0.349
2015NS21.1319.9320.9119.60–23.7874.88 ****0.736
DS20.4519.1720.2818.97–22.6846.98 **1.110
Plant height (Cm)2014NS56.3955.4256.0054.28–59.8726.98 ****3.656
DS49.8050.1350.2349.31–52.118.26 **2.371
2015NS51.7549.8750.4948.71–53.4710.51 ****6.116
DS45.7847.3745.6544.02–49.098.00 **6.347
Phenology
Number of days for flowering2014NS48.0248.8748.1547.59–48.1592.33 ****0.168
DS48.0248.0348.2147.92–49.0234.46 **0.130
2015NS48.899.2249.0048.80–49.24ns**0.241
DS48.7148.7748.6448.39–48.895.87 *0.196
Number of days to maturity2014NS79.0080.9879.6778.93–80.9037.3 ****0.612
DS76.7078.1177.6776.25–80.01141.25 **0.641
2015NS87.9589.3188.2287.23–91.1611.10 ****1.640
DS94.4294.4094.4394.40–94.46ns0.153
Number of days to grain filling2014NS31.3632.1731.7231.33–32.4510.03 ****0.376
DS28.8429.7429.3928.29–31.5173.86 **0.620
2015NS39.2839.2739.0438.41–41.013.87 ***1.307
DS46.0346.0146.0546.01–46.09ns0.147
* and ** = significant by the chi-square test with one degree of freedom: at 5% (3.84) and 1% (6.63), respectively; ns = not significantly; PEV = variance of the error of prediction for genotypic values; Con = condition; RIL = recombinant inbred line; § Gê = interaction between genotype treatments and stress condition.
Table 3. Pearson’s correlation coefficients between the average values of the characteristics evaluated in the RIL population IAPAR 81 × LP97-28 cultivated under drought stress (DS) and no-drought stress (NS) conditions in Maringá, PR in the years 2014 and 2015. Values above the diagonal are the correlation coefficients between the traits under drought stress (DS); values below the diagonal are the correlation coefficients between traits no-drought stress (NS).
Table 3. Pearson’s correlation coefficients between the average values of the characteristics evaluated in the RIL population IAPAR 81 × LP97-28 cultivated under drought stress (DS) and no-drought stress (NS) conditions in Maringá, PR in the years 2014 and 2015. Values above the diagonal are the correlation coefficients between the traits under drought stress (DS); values below the diagonal are the correlation coefficients between traits no-drought stress (NS).
TraitsGYSYDNPSWPHNDFNDMNDG
GY -0.07 ns0.00 ns0.05 ns0.07 ns−0.04 ns0.00 ns−0.01 ns
SYD0.73 *** 0.72 ***0.46 ***0.14 ns−0.28 ***−0.34 **−0.27 ***
NP0.65 ***0.57 *** 0.21 **0.11 ns−0.11 ns−0.21 **−0.21 **
SW0.42 ***0.35 ***0.29 *** 0.18 *−0.25 **−0.21 **−0.13 ns
PH0.09 ns0.08 ns0.02 ns0.14 ns −0.07 ns−0.04 ns−0.01 ns
NDF−0.16 *−0.15 *−0.13 ns−0.21 **0.00 ns 0.04 ***0.07 ns
NDM−0.01 ns−0.07 ns0.03 ns−0.20 *0.09 ns0.37 *** 0.93 **
NDE0.08 ns0.04 ns0.05 ns−0.10 ns0.08 ns−0.14 ns0.84 ***
*, ** and *** = significance at 0.05, 0.01 and 0.001 probability levels, respectively; ns = not significant; GY = grain yield; SYD = seed yield per day; NP = number of pods per plant; SW = 100-seed weight; PH = plant height; NDF = number of days for flowering; NDM = number of days to mature; and NDG = number of days to grain filling.
Table 4. Number of SNP markers, chromosome size and average distance between markers for the 11 common bean chromosomes for the RIL population derived from the IAPAR 81 × LP97-28 cross genotyped with the BARCBean6K_3 BeadChip.
Table 4. Number of SNP markers, chromosome size and average distance between markers for the 11 common bean chromosomes for the RIL population derived from the IAPAR 81 × LP97-28 cross genotyped with the BARCBean6K_3 BeadChip.
ChromosomeNumber of SNP
Markers
Chromosome SizeAverage Distance between Markers
––––––––––––––– cM ––––––––––––––
Pv0143117.42.73
Pv02195101.40.52
Pv036184.11.38
Pv042313.10.57
Pv053171.12.29
Pv068495.31.13
Pv077291.31.27
Pv0813266.40.50
Pv094961.31.25
Pv106066.61.11
Pv112347.82.07
Total773815.8
Average 1.34
Table 5. Description of the quantitative trait loci (QTL) for drought tolerance identified in the common bean RIL population IAPAR 81 × LP97-28 (IL) cultivated under drought stress (DS) and no-drought stress (NS) conditions in Maringá, PR in the agricultural years of 2014 and 2015.
Table 5. Description of the quantitative trait loci (QTL) for drought tolerance identified in the common bean RIL population IAPAR 81 × LP97-28 (IL) cultivated under drought stress (DS) and no-drought stress (NS) conditions in Maringá, PR in the agricultural years of 2014 and 2015.
§ Traits/QTLYearCond. *ChrPeak of QTL PositionQTL
INTERVAL
Nearest
SNP Marker
# LOD ScoreLOD
Threshold
¥ Add R2
––––––––cM–––––––––––– %
Grain yield
GY9IL2014NS959.146.7–59.9ss7156403023.133.1−136.678.7
Yield per day
SYD7L2014DS752.1 49.8–57.1ss7156404873.83 0.6112.1
SYD8ILCombinedDS843.028.1–49.2ss7156482313.532.99−0.569.0
SYD9IL2014NS922.513.5–27.3ss7156493573.17 0.348.3
100-seed weight
SW7IL2014DS761.857.8–63.3ss7156393857.32.99−0.6115.4
2015NS761.050.6–64.1ss7156477284.54 −0.7211.2
2015DS761.057.3–66.1ss7156477283.13 −0.616.9
CombinedNS761.048.1–61.4ss7156477284.253.33−0.5410.2
CombinedDS761.058.0–61.4ss7156477286.42 −0.6213.8
SW8IL2014DS83.62.1–5.1ss7156496049.444.67−0.7320.6
2014DS86.45.1–9.5ss7156482209.254.02−0.7521.0
2015NS85.10.6–5.4ss7156485494.493.42−0.7510.5
2015NS87.15.4–11.9ss7156489294.873.6−0.7811.3
2015DS84.41.9–5.1ss7156482193.98 −0.749.4
2015DS811.16.5–16.5ss7156489296.333.92−0.9716.6
CombinedNS82.00.0–3.6ss7156507833.33 −0.58.5
CombinedNS87.15.4–12.1ss7156489293.7 −0.528.7
CombinedNS80.00.0–1.6ss7156480433.08 −0.137.5
CombinedDS84.42.7–5.1ss7156482198.654.87−0.7519.3
CombinedDS86.45.1–12.2ss7156482208.964.67−0.7820.5
CombinedDS815.314.3–16.4ss7156493786.89 −0.7117.2
Number of pods per plant
NP11IL2014DS1138.631.1–41.4ss7156474654.882.650.7113.0
2015DS1138.631.1–41.4ss7156474654.88 0.7113.0
Plant height
PH1IL2014DS167.061.0–79.7ss7156468683.15 −1.488.9
2015DS167.061.0–79.7ss7156468683.15 −1.488.9
CombinedDS168.061.0–78.3ss7156468684.342.93−1.0211.2
PH2.1ILCombinedNS219.418.5–29.8ss7156409413.48 −1.268.7
PH2.2IL2014DS296.793.8–98.8ss7156461443.953.69−1.539.8
2014DS299.298.8–100.2ss7156469293.513.24−1.468.8
2015DS296.793.8–98.8ss7156461443.953.6−1.539.8
2015DS299.298.8–100.2ss7156469293.513.24−1.468.8
CombinedNS298.190.4–100.2ss7156451863.45 −0.958.5
CombinedDS296.794.9–98.8ss7156461444.183.19−0.929.5
PH3ILCombined DS366.460.7–70.4ss7156505804.324.540.959.5
Number of days for flowering
NDF2IL2014NS239.729.0–42.8ss7156488195.553.070.3414.1
2014NS248.443.5–52.7ss7156475264.263.170.310.2
CombinedNS228.822.9–35.8ss7156476523.51 0.2913.1
NDF3IL2014NS348.442.4–50.5ss7156476895.02 −0.3212.3
NDF10ILCombinedDS100.00.0–7.0ss7156397833.382.880.178.6
Number of days to maturity
NDM7IL2014DS754.1150.6–57.2ss7156404873.74 0.6612.1
2015NS760.3155–65.3ss7156477344.863.380.9612.5
2015NS773.4171.2–74.5ss7156485533.86 0.8710.3
CombinedNS759.3150.8–64.2ss7156486194.733.340.5313.2
CombinedDS761.8160.3–65.9ss7156393853.17 0.538.4
* DS—drought stress condition; NS—no-drought stress condition; Chr—chromosome; # LOD, odds log. LOD thresholds were calculated by the performance of 1000 permutations with p = 0.05; ¥ effect of replacing a single allele from one parent to another. Positive values indicate allele of origin for IAPAR 81 and negative for LP97-28; portion of the phenotypic variance explained by the QTL LOD peak using Win Cartographer’s composite interval mapping; § common bean QTL nomenclature guidelines (http://bic.css.msu.edu/_pdf/Guidelines_QTL_Nomenclature.pdf (accessed on 1 December 2017)).
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Elias, J.C.F.; Gonçalves-Vidigal, M.C.; Vaz Bisneta, M.; Valentini, G.; Vidigal Filho, P.S.; Gilio, T.A.S.; Moda-Cirino, V.; Song, Q. Genetic Mapping for Agronomic Traits in IAPAR 81/LP97-28 Population of Common Bean (Phaseolus vulgaris L.) under Drought Conditions. Plants 2021, 10, 1568. https://doi.org/10.3390/plants10081568

AMA Style

Elias JCF, Gonçalves-Vidigal MC, Vaz Bisneta M, Valentini G, Vidigal Filho PS, Gilio TAS, Moda-Cirino V, Song Q. Genetic Mapping for Agronomic Traits in IAPAR 81/LP97-28 Population of Common Bean (Phaseolus vulgaris L.) under Drought Conditions. Plants. 2021; 10(8):1568. https://doi.org/10.3390/plants10081568

Chicago/Turabian Style

Elias, Júlio César Ferreira, Maria Celeste Gonçalves-Vidigal, Mariana Vaz Bisneta, Giseli Valentini, Pedro Soares Vidigal Filho, Thiago Alexandre Santana Gilio, Vânia Moda-Cirino, and Qijian Song. 2021. "Genetic Mapping for Agronomic Traits in IAPAR 81/LP97-28 Population of Common Bean (Phaseolus vulgaris L.) under Drought Conditions" Plants 10, no. 8: 1568. https://doi.org/10.3390/plants10081568

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop