<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" article-type="research-article" dtd-version="2.3" xml:lang="EN">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Plant Sci.</journal-id>
<journal-title>Frontiers in Plant Science</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Plant Sci.</abbrev-journal-title>
<issn pub-type="epub">1664-462X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fpls.2023.1236147</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Plant Science</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Variability in drought stress response in a panel of 100 faba bean genotypes</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Balko</surname>
<given-names>Christiane</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="author-notes" rid="fn001">
<sup>*</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2182230"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Torres</surname>
<given-names>Ana M.</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/31647"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Gutierrez</surname>
<given-names>Natalia</given-names>
</name>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2053012"/>
</contrib>
</contrib-group>
<aff id="aff1">
<sup>1</sup>
<institution>Julius K&#xfc;hn-Institut (JKI), Federal Research Centre for Cultivated Plants, Institute for Resistance Research and Stress Tolerance</institution>, <addr-line>Sanitz</addr-line>, <country>Germany</country>
</aff>
<aff id="aff2">
<sup>2</sup>
<institution>&#xc1;rea de Mejora Vegetal y Biotecnolog&#xed;a, Andalusian Institute of Agricultural and Fisheries Research and Training (IFAPA)</institution>, <addr-line>Centro Alameda del Obispo, C&#xf3;rdoba</addr-line>, <country>Spain</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>Edited by: Leif Skot, Aberystwyth University, United Kingdom</p>
</fn>
<fn fn-type="edited-by">
<p>Reviewed by: Elsayed Mansour, Zagazig University, Egypt; Odd Arne Rognli, Norwegian University of Life Sciences, Norway</p>
</fn>
<fn fn-type="corresp" id="fn001">
<p>*Correspondence: Christiane Balko, <email xlink:href="mailto:christiane.balko@julius-kuehn.de">christiane.balko@julius-kuehn.de</email>
</p>
</fn>
</author-notes>
<pub-date pub-type="epub">
<day>30</day>
<month>08</month>
<year>2023</year>
</pub-date>
<pub-date pub-type="collection">
<year>2023</year>
</pub-date>
<volume>14</volume>
<elocation-id>1236147</elocation-id>
<history>
<date date-type="received">
<day>07</day>
<month>06</month>
<year>2023</year>
</date>
<date date-type="accepted">
<day>14</day>
<month>08</month>
<year>2023</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2023 Balko, Torres and Gutierrez</copyright-statement>
<copyright-year>2023</copyright-year>
<copyright-holder>Balko, Torres and Gutierrez</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p>
</license>
</permissions>
<abstract>
<p>Faba bean is an important protein crop for food and feed worldwide and provides a range of advantages in crop rotations. Its limited use in modern agriculture is mainly due to the high fluctuations in yield. A well known limiting factor in most legumes, and particularly in faba bean, is the high sensitivity to water shortage, which is further aggravated by climate change. The present study was undertaken to exploit the genetic variation in drought stress response in a faba bean collection of 100 accessions with diverse origins and to assess selection criteria for identifying drought tolerant genotypes. Physiological, phenological and yield related traits evaluated under drought or water-sufficient conditions responded significantly to the end-terminated drought stress. Comparison of yield relations showed the advantage of using a stress tolerance index (STI) to identify genotypes combining high yield potential with high stress yield. With regard to physiological traits, SPAD (chlorophyll content) values were significantly related to yield as well as to STI, while the other traits also contributed to different extents to variation in yield formation. Among the yield related traits, seeds per plant proved to be the most important trait followed by pods per plant. Interestingly, the eight genotypes with the best STI performance use different strategies to cope with drought stress.</p>
</abstract>
<kwd-group>
<kwd>faba bean</kwd>
<kwd>drought stress</kwd>
<kwd>phenotyping</kwd>
<kwd>physiological traits</kwd>
<kwd>senescence</kwd>
<kwd>genetic diversity</kwd>
<kwd>heritability</kwd>
<kwd>yield</kwd>
</kwd-group>
<counts>
<fig-count count="10"/>
<table-count count="5"/>
<equation-count count="2"/>
<ref-count count="86"/>
<page-count count="15"/>
<word-count count="7536"/>
</counts>
<custom-meta-wrap>
<custom-meta>
<meta-name>section-in-acceptance</meta-name>
<meta-value>Plant Breeding</meta-value>
</custom-meta>
</custom-meta-wrap>
</article-meta>
</front>
<body>
<sec id="s1" sec-type="intro">
<title>Introduction</title>
<p>Grain legumes are a major source of plant protein for food and feed worldwide and provide multiple advantages in crop rotations regarding soil fertility, plant health and sustainability. Faba bean (<italic>Vicia faba</italic> L.) is one of the most widely distributed crops, being grown both as a grain (pulse) and green-manure legume (<xref ref-type="bibr" rid="B34">Jensen et&#xa0;al., 2010</xref>). With a production varying between 4.5 and 5.5 Mt in the last decade, faba bean ranks 6th in terms of world production of pulse crops (<xref ref-type="bibr" rid="B22">FAOSTAT, 2022</xref>). China is the largest producer followed by Ethiopia, the United Kingdom, Australia, France and Germany. In the last 70 years, there has been more than 50% decline in the global cropping area of the faba bean, which is mainly caused by poor yield stability owing to climate variability, diseases, weeds and other pests, which represent the major constraints of faba bean production.</p>
<p>Faba bean yield fluctuations are likely to increase with the predicted climate change in this century. Increasing atmospheric CO2 concentrations will change surface temperatures and precipitation patterns, causing more periods of extreme precipitation and drought (<xref ref-type="bibr" rid="B23">Farooq et&#xa0;al., 2014</xref>). Therefore, drought stress has become one of the most uncontrollable and unpredictable agricultural challenges and is today one of the major yield limiting factors in grain legume crops (<xref ref-type="bibr" rid="B77">Swann et&#xa0;al., 2016</xref>; <xref ref-type="bibr" rid="B55">Nadeem et&#xa0;al., 2019</xref>). Improving drought tolerance is a key strategy to enhance performance and stability of yield in these crops. Nevertheless, the complexity of the trait, the lack of efficient selection protocols and the mixed faba bean breeding system pose major challenges for effective implementation in plant breeding programs (<xref ref-type="bibr" rid="B51">Muktadir et&#xa0;al., 2020</xref>).</p>
<p>Legumes respond to drought with morphological, physiological and biochemical changes in roots, stems and leaves (<xref ref-type="bibr" rid="B71">Shakeel et&#xa0;al., 2011</xref>). Drought stress results in many interactive modifications including changes in the expression of drought-resistance genes, the synthesis of hormones, the overproduction of reactive oxygen species and the osmotic adjustments through active ions or organic compound such as proline and carbohydrates (<xref ref-type="bibr" rid="B24">Farooq et&#xa0;al., 2012</xref>; <xref ref-type="bibr" rid="B37">Kaur and Asthir, 2015</xref>; <xref ref-type="bibr" rid="B83">Wahab et&#xa0;al., 2022</xref>). Plants minimize water loss by closing stomata and reducing light absorbance, canopy leaf area and peaked water absorption.</p>
<p>The extent and type of responses will depend on the intensity and length of the stress, but in all cases grain yield is substantially reduced. A meta-analysis of legume yield responses to drought under field conditions from 1980 to 2014 revealed a yield reduction of 40% following a 65% reduction in water availability, with cultivar and environmental factors being important cofactors (<xref ref-type="bibr" rid="B20">Daryanto et&#xa0;al., 2015</xref>). The study confirmed, that legumes have a high yield potential, but respond in a sensitive way to water shortage (<xref ref-type="bibr" rid="B52">M&#xfc;ller et&#xa0;al., 1985</xref>; <xref ref-type="bibr" rid="B28">Grashoff, 1990</xref>; <xref ref-type="bibr" rid="B39">Khan et&#xa0;al., 2010</xref>; <xref ref-type="bibr" rid="B80">Torres et&#xa0;al., 2010</xref>; <xref ref-type="bibr" rid="B62">Redden et&#xa0;al., 2014</xref>; <xref ref-type="bibr" rid="B51">Muktadir et&#xa0;al., 2020</xref>). Faba bean is more sensitive to drought than other field crops and largely exceeds other grain legumes such as common bean, pea and chickpea (<xref ref-type="bibr" rid="B46">McDonald and Paulsen, 1997</xref>; <xref ref-type="bibr" rid="B8">Amede and Schubert, 2003</xref>; <xref ref-type="bibr" rid="B44">Manning et&#xa0;al., 2020</xref>).</p>
<p>Drought impairs grain legume yield at all growth stages. Germination, leaf area and photosynthetic activity are significantly reduced (<xref ref-type="bibr" rid="B55">Nadeem et&#xa0;al., 2019</xref>). However, the most sensitive period is the reproductive stage where drought stress leads to earlier flowering, and ultimately, reduced pod and grain set (<xref ref-type="bibr" rid="B59">Plies-Balzer et&#xa0;al., 1995</xref>; <xref ref-type="bibr" rid="B38">Khan et&#xa0;al., 2007</xref>; <xref ref-type="bibr" rid="B6">Alghamdi et&#xa0;al., 2015</xref>). Against this background, improvement of drought tolerance is an important faba bean breeding goal. Understanding of drought response patterns and associated traits are key factors for achieving higher yield stability under stress conditions. However, compared to other crops, improvement of productivity under drought stress has rarely been included in faba bean breeding programmes. Breeding progress has been relatively slow so far due to the limited number of genotypes included in the studies (<xref ref-type="bibr" rid="B1">Abdelmula et&#xa0;al., 1999</xref>), the low heritability of advantageous traits and the lack of efficient and reproducible screening methods (<xref ref-type="bibr" rid="B76">Stoddard et&#xa0;al., 2006</xref>; <xref ref-type="bibr" rid="B51">Muktadir et&#xa0;al., 2020</xref>). However, a wide genotypic variation in faba bean water stress response has been reported (<xref ref-type="bibr" rid="B42">Link et&#xa0;al., 1999</xref>; <xref ref-type="bibr" rid="B9">Amede et&#xa0;al., 2001</xref>; <xref ref-type="bibr" rid="B64">Ricciardi et&#xa0;al., 2001</xref>; <xref ref-type="bibr" rid="B38">Khan et&#xa0;al., 2007</xref>), indicating a high potential for breeding in drought-prone environments. To exploit this genetic variation, the application of an accurate and relevant phenotyping method plays a key role for the selection of drought-resilient genotypes and for the dissection and genetic analysis of a complex trait such as the adaptive response of crops to drought (<xref ref-type="bibr" rid="B82">Tuberosa, 2012</xref>).</p>
<p>Phenotypic attributes are the most frequently used criteria used for identifying drought-tolerant genotypes. In legumes, selection for drought tolerance based on highly heritable morphological and phenological traits, together with physiological attributes such as the accumulation of proline or soluble sugars, has proven highly successful for screening genotypes under limited water supply (<xref ref-type="bibr" rid="B41">Lafitte et&#xa0;al., 2003</xref>; <xref ref-type="bibr" rid="B65">Richards, 2006</xref>; <xref ref-type="bibr" rid="B76">Stoddard et&#xa0;al., 2006</xref>; <xref ref-type="bibr" rid="B11">Annicchiarico and Iannucci, 2008</xref>; <xref ref-type="bibr" rid="B5">Alderfasi and Alghamdi, 2010</xref>; <xref ref-type="bibr" rid="B10">Ammar et&#xa0;al., 2015</xref>; <xref ref-type="bibr" rid="B7">Ali et&#xa0;al., 2016</xref>). However, phenotyping under field conditions is time- and labour-consuming, and the reproducibility of water stress conditions is often poor, because of variations in the timing of onset, duration and severity of the drought (<xref ref-type="bibr" rid="B38">Khan et&#xa0;al., 2007</xref>). The recent use of a controlled water supply combined with rainout shelters reduces these uncertainties and generates results that are more accurate. These tools allow phenotyping of large field populations under conditions of adequate light intensity and quality, avoiding the effects of unpredictable rainfall patterns.</p>
<p>Direct selection for grain yield under drought conditions often proves inefficient to identify stress-tolerant genotypes. For this reason, several selection indices have been developed to evaluate yield stability based on grain yield under normal and stress conditions (reviewed in <xref ref-type="bibr" rid="B75">Sofi et&#xa0;al., 2018</xref>). Among these, the stress tolerance index (STI) defined by <xref ref-type="bibr" rid="B26">Fernandez (1992)</xref> has been consistently correlated with other indexes (<xref ref-type="bibr" rid="B25">Farshadfar et&#xa0;al., 2013</xref>; <xref ref-type="bibr" rid="B45">Mansour et&#xa0;al., 2021</xref>; <xref ref-type="bibr" rid="B72">Sharifi et&#xa0;al., 2021</xref>; <xref ref-type="bibr" rid="B47">Memari et&#xa0;al., 2022</xref>), indicating that it can be used as an alternative to select drought tolerant genotypes with high yield performance in stress and non-stress conditions.</p>
<p>Considering the scarce number of studies so far reported and the limited number of faba bean genotypes used, breeding progress for drought stress tolerance in faba bean has been relatively slow. The goal of the present study was to exploit the genetic variation of a wider collection of 100 faba bean accessions with diverse origins to assess the selection criteria for identifying drought tolerant genotypes and to dissect a complex trait such as the adaptive response to drought. Thus, our main aims were: (i) to evaluate the performance of this faba bean collection, one of the largest screenings reported so far, under water-sufficient and drought stress conditions and (ii) to investigate the relevance of morphological and phenological traits with respect to yield and yield stability in water stress conditions in an effort to identify and to improve selection efficiency under drought stress environments.</p>
</sec>
<sec id="s2" sec-type="materials|methods">
<title>Materials and methods</title>
<sec id="s2_1">
<title>Plant material and experimental design</title>
<p>The test set included 100 faba bean inbred lines of diverse origins, which were grown in 2019 and 2020 in Gro&#xdf; L&#xfc;sewitz (north-eastern Germany). The location is characterized by an annual mean temperature of 8.3&#xb0;C. Soil type is a slightly loamy sand with a pH of 5.7.</p>
<p>The worldwide diversity panel includes genetic stocks aiming at gathering a wide range of genetic diversity from diverse geographical origins: Africa (8 accessions), North and South America (2), Asia (27) and Europe (39). The remaining 24 accessions were provided by the International Center for Agricultural Research in the Dry Areas (ICARDA), but their exact origin is unknown. Europe, with 9 countries, is the most represented geographical area, followed by Asia, Africa and America (7, 4 and 2 countries, respectively). Spain accounts for the highest number of accessions (23). The diversity panel was built in collaboration with the public institutes ICARDA, the Andalusian Institute of Agricultural and Fisheries Research and Training (IFAPA) and the French National Institute for Agricultural Research (INRA). Prior to the genotyping analysis, all the Spanish lines had been self-pollinated for at least four generations. The remaining accessions were purified for two generations by single seed descent (SSD) in insect-proof cages. A detailed description of this collection is provided in <xref ref-type="supplementary-material" rid="ST1">
<bold>Supplementary Table S1</bold>
</xref>.</p>
<p>The experimental set up was a randomized complete-block design with 4 replications in both treatments. Drought stress treatment was performed in two mobile rainout shelters, with two blocks in each. The same design was applied in the open field control treatment. Plants were grown in single row plots of 1.2 m length with 14 plants each (one seed per hill) and the distance between plots was 0.5 m. Additional water was applied by drip irrigation (<xref ref-type="supplementary-material" rid="SF1">
<bold>Supplementary Figure S1</bold>
</xref>). Irrigation has been scheduled in the range of 60 - 70% of field capacity of the soil, determined over winter after excessive rainfall. Water content in the soil was assessed by time-domain reflectance (TDR) probes in about 0.40 m depth (<xref ref-type="fig" rid="f1">
<bold>Figure&#xa0;1</bold>
</xref>). Irrigation in the stress treatment was stopped when about 30% of the genotypes had started flowering. During rainfall the shelter covered the plots. An end-terminated drought stress was applied. The key data of the experiments are explained in <xref ref-type="supplementary-material" rid="ST2">
<bold>Supplementary Table S2</bold>
</xref>.</p>
<fig id="f1" position="float">
<label>Figure&#xa0;1</label>
<caption>
<p>Soil water content of the experimental areas measured by TDR probes in about 0.40 m depth for the first and second experimental year.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-14-1236147-g001.tif"/>
</fig>
</sec>
<sec id="s2_2">
<title>Weather conditions</title>
<p>The mean temperature, the Crop Heat Unit (CHU) based heat sums and the precipitation conditions during the vegetative period in the two years of the experiment are summarized in <xref ref-type="table" rid="T1">
<bold>Table&#xa0;1</bold>
</xref>. The experimental year 2019 was moderate with respect to rainfall but warm with an early heat period in June. The year 2020 was cooler with slightly more rainfall and a late heat period in August, when faba beans were partly harvested already. This resulted in a more intensive drought stress in 2019 compared to 2020.</p>
<table-wrap id="T1" position="float">
<label>Table&#xa0;1</label>
<caption>
<p>Weather conditions during the vegetative period in 2019 and 2020 at Gross L&#xfc;sewitz.</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="bottom" align="left"/>
<th valign="bottom" colspan="2" align="center">Mean Temperature (&#xb0;C)</th>
<th valign="bottom" colspan="2" align="center">CHU* heat sums (&#xb0;C)</th>
<th valign="bottom" colspan="2" align="center">Precipitation (mm)</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="bottom" align="left"/>
<td valign="middle" align="center">2019</td>
<td valign="middle" align="center">2020</td>
<td valign="middle" align="center">2019</td>
<td valign="middle" align="center">2020</td>
<td valign="middle" align="center">2019</td>
<td valign="middle" align="center">2020</td>
</tr>
<tr>
<td valign="bottom" align="left">
<bold>April</bold>
</td>
<td valign="bottom" align="center">9.2</td>
<td valign="bottom" align="center">8.4</td>
<td valign="bottom" align="center">259</td>
<td valign="bottom" align="center">188</td>
<td valign="bottom" align="center">23.9</td>
<td valign="bottom" align="center">10.9</td>
</tr>
<tr>
<td valign="bottom" align="left">
<bold>May</bold>
</td>
<td valign="bottom" align="center">10.9</td>
<td valign="bottom" align="center">10.9</td>
<td valign="bottom" align="center">311</td>
<td valign="bottom" align="center">279</td>
<td valign="bottom" align="center">33.1</td>
<td valign="bottom" align="center">20.5</td>
</tr>
<tr>
<td valign="bottom" align="left">
<bold>June</bold>
</td>
<td valign="bottom" align="center">19.1</td>
<td valign="bottom" align="center">16.6</td>
<td valign="bottom" align="center">665</td>
<td valign="bottom" align="center">579</td>
<td valign="bottom" align="center">87.4</td>
<td valign="bottom" align="center">131.2</td>
</tr>
<tr>
<td valign="bottom" align="left">
<bold>July</bold>
</td>
<td valign="bottom" align="center">17.7</td>
<td valign="bottom" align="center">16.0</td>
<td valign="bottom" align="center">644</td>
<td valign="bottom" align="center">571</td>
<td valign="bottom" align="center">79.7</td>
<td valign="bottom" align="center">62.1</td>
</tr>
<tr>
<td valign="bottom" align="left">
<bold>August</bold>
</td>
<td valign="bottom" align="center">18.5</td>
<td valign="bottom" align="center">19.4</td>
<td valign="bottom" align="center">685</td>
<td valign="bottom" align="center">708</td>
<td valign="bottom" align="center">55.7</td>
<td valign="bottom" align="center">82.2</td>
</tr>
<tr>
<td valign="middle" align="left">
<bold>Sum</bold>
</td>
<td valign="middle" colspan="2" align="center"/>
<td valign="middle" align="center">2564</td>
<td valign="middle" align="center">2325</td>
<td valign="middle" align="center">279.8</td>
<td valign="middle" align="center">306.9</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn>
<p>*Crop Heat Units.</p>
</fn>
</table-wrap-foot>
</table-wrap>
</sec>
<sec id="s2_3">
<title>Measurements of traits</title>
<p>Four physiological traits were assessed in the youngest fully-developed leaves from the main stem of the inner 10 plants of each plot. Two leaf samples mixed from 5 plants per plot each were taken about 4 weeks after onset of stress, according to the stress development, for determining the content of free proline (PRO) and total soluble sugars (TSS). Samples were immediately frozen in liquid nitrogen. PRO was analyzed as described by <xref ref-type="bibr" rid="B16">Bates et&#xa0;al. (1973)</xref> and TSS was estimated by the anthrone reagent method according to <xref ref-type="bibr" rid="B85">Yemm and Willis (1954)</xref>. SPAD measurements as an indirect parameter for chlorophyll content were carried out with a Chlorophyll Meter SPAD 502 Plus (Konica Minolta) at the beginning of stress treatment (SPAD1) and in the same time slot of TSS and PRO sampling (SPAD2). The mean value of 30 clippings per plot was used. The difference between these SPAD values, DiffSPAD (DiffSPAD = SPAD1 &#x2013; SPAD2), was calculated to standardize the possible differences in chlorophyll content among genotypes and to obtain a more accurate measure of chlorophyll degradation (senescence). Two phenological traits were analyzed: end of flowering (EF), recorded as the days after sowing when no open flowers were visible any more in the plot (BBCH 69); and maturity date (MAT), defined as the date when more than 90% of the pods have ripened (BBCH 89) following the extended BBCH Code (<xref ref-type="bibr" rid="B30">Hack et&#xa0;al., 1992</xref>). Moreover, five morphological traits were recorded: plant height in cm (PH); number of pods per plant (PP); number of seeds per plant (SP); hundred seed weight (HSW) in grams; and plot yield (PY) in kg. Plant height was measured close to maturity (BBCH 80-85) from the ground to the tip of 10 inner plants per plot. At the end of the vegetative period, 10 single plants per plot were harvested by hand and pods and seeds were counted. Dry matter content was determined at 105&#xb0;C. PY and HSW are given on a basis of 86% dry matter.</p>
</sec>
<sec id="s2_4">
<title>Data analyses</title>
<p>Statistical analyses were performed by SAS 9.4 (SAS Institute, USA). Analysis of variance (ANOVA) used a linear mixed model with year as a random effect. Means are shown as LSmeans. Minimum variance quadratic unbiased estimators (MIVQUEs) of variance components were used to calculate broad sense heritability h&#xb2; according to <xref ref-type="bibr" rid="B17">Becker (2011)</xref>. The response of the accessions to drought stress was calculated as relative yield of the drought stress treatment:</p>
<disp-formula>
<mml:math display="block" id="M1">
<mml:mrow>
<mml:mi>Y</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>l</mml:mi>
<mml:mo>=</mml:mo>
<mml:mi>P</mml:mi>
<mml:mi>Y</mml:mi>
<mml:mo>_</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mi>D</mml:mi>
<mml:mi>S</mml:mi>
<mml:mo>&#xd7;</mml:mo>
<mml:mn>100</mml:mn>
<mml:mo stretchy="false">/</mml:mo>
<mml:mi>P</mml:mi>
<mml:mi>Y</mml:mi>
<mml:mo>_</mml:mo>
<mml:mi>C</mml:mi>
</mml:mrow>
</mml:math>
</disp-formula>
<p>with PY_DS as plot yield of the drought stress and PY_C as plot yield of control treatment.</p>
<p>The Stress Tolerance Index (STI) was calculated for the yield measured in both treatments (control and drought) according to <xref ref-type="bibr" rid="B26">Fernandez (1992)</xref> as follows:</p>
<disp-formula>
<mml:math display="block" id="M2">
<mml:mrow>
<mml:mi>S</mml:mi>
<mml:mi>T</mml:mi>
<mml:mi>I</mml:mi>
<mml:mo>=</mml:mo>
<mml:mi>P</mml:mi>
<mml:mi>Y</mml:mi>
<mml:mo>_</mml:mo>
<mml:mi>C</mml:mi>
<mml:mo>&#xd7;</mml:mo>
<mml:mi>P</mml:mi>
<mml:mi>Y</mml:mi>
<mml:mo>_</mml:mo>
<mml:mi>D</mml:mi>
<mml:mi>S</mml:mi>
<mml:mo stretchy="false">/</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mo stretchy="false">(</mml:mo>
<mml:mover accent="true">
<mml:mi>&#x3d2;</mml:mi>
<mml:mo>&#xaf;</mml:mo>
</mml:mover>
<mml:mo>_</mml:mo>
<mml:mi>C</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:math>
</disp-formula>
<p>with PY _ C as plot yield of control, PY_DS as plot yield of the drought stress treatment and <inline-formula>
<mml:math display="inline" id="im1">
<mml:mrow>
<mml:mover accent="true">
<mml:mi>&#x3d2;</mml:mi>
<mml:mo>&#xaf;</mml:mo>
</mml:mover>
<mml:mo>_</mml:mo>
<mml:mi>C</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>as mean yield of the control environment.</p>
<p>Pearson&#x2019;s correlations between STI, PY_C and PY_DS data were calculated in order to check for possible dependencies between the performance of the accessions under control conditions and the strength of response to drought stress. A principal component analysis (PCA) was calculated based on correlations using JMP Genomics 9 (SAS Institute, USA). The results were represented using a biplot that combines both the principal component scores and the loading vectors in a single plot. In addition, a cluster analysis of observations was performed to find subgroups within phenotypic data in control and drought conditions by means of the K-means clustering method (<xref ref-type="bibr" rid="B32">Hartigan and Wong, 1979</xref>). The data were scaled to make variables comparable. The average silhouette approach was used to determine the optimal number of clusters. K-means clustering was performed using the built-in k-means libraries of R suite and similarities were measured by squared Euclidean distance to classify all the accessions into groups (<xref ref-type="bibr" rid="B61">R Core Team, 2022</xref>). The results were visualised by means of the &#x2018;fviz_cluster&#x2019; function in the &#x2018;factoextra&#x2019; R-package (<xref ref-type="bibr" rid="B40">Kassambara and Mundt, 2020</xref>).</p>
</sec>
</sec>
<sec id="s3" sec-type="results">
<title>Results</title>
<sec id="s3_1">
<title>Statistics of the data</title>
<p>Descriptive statistics and analysis of variance (ANOVA) of the traits recorded over two years is shown in <xref ref-type="table" rid="T2">
<bold>Table&#xa0;2</bold>
</xref>. The table includes the mean, range and Least Significant Difference (LSD) values together with the heritability (h&#xb2;) of the traits in both conditions. LSDs of all the traits studied in the 100 faba bean genotypes in control and drought stress conditions are provided in <xref ref-type="supplementary-material" rid="ST3">
<bold>Supplementary Table S3</bold>
</xref>. Most of the parameters displayed a wide range of differences among the treatments. A significant effect of genotype and drought stress treatment as well as an interaction between them could be observed for all traits.</p>
<table-wrap id="T2" position="float">
<label>Table&#xa0;2</label>
<caption>
<p>Descriptive statistics and analysis of variance (ANOVA) of the traits over both years.</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="top" rowspan="2" colspan="2" align="left">Trait&#xa0;<break/>&#xa0;</th>
<th valign="top" rowspan="2" align="center">Control</th>
<th valign="top" rowspan="2" align="center">Stress</th>
<th valign="top" rowspan="2" align="center">Unit<break/>&#xa0;</th>
<th valign="top" colspan="4" align="center">F values</th>
<th valign="top" colspan="2" align="center">Heritability (h&#xb2;)</th>
</tr>
<tr>
<th valign="top" align="center">Genotype</th>
<th valign="top" align="center">Treatment</th>
<th valign="top" align="center">G x T</th>
<th valign="top" align="center">Year</th>
<th valign="top" align="center">Control</th>
<th valign="top" align="center">Stress</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">Soluble sugars<break/>(TSS)</td>
<td valign="top" align="left">
<break/>Range<break/>LSD</td>
<td valign="top" align="center">1117.6<break/>769.5 - 1709.4<break/>587.1</td>
<td valign="top" align="center">1387.9<break/>998.2 - 2242.4<break/>659.9</td>
<td valign="top" align="center">&#xb5;mol/ g DM</td>
<td valign="top" align="center">7.44 ***</td>
<td valign="top" align="center">339.79 ***</td>
<td valign="top" align="center">1.86 ***</td>
<td valign="top" align="center">196.73 ***<break/>2019&gt;2020</td>
<td valign="top" align="center">0.56</td>
<td valign="top" align="center">0.66</td>
</tr>
<tr>
<td valign="top" align="left">Free proline<break/>(PRO)</td>
<td valign="top" align="left">
<break/>Range<break/>LSD</td>
<td valign="top" align="center">2.29<break/>1.50 - 6.55<break/>1.20</td>
<td valign="top" align="center">8.45<break/>1.79 - 57.02<break/>26.75</td>
<td valign="top" align="center">&#xb5;mol/ g DM</td>
<td valign="top" align="center">2.93 ***</td>
<td valign="top" align="center">171.09 ***</td>
<td valign="top" align="center">2.57 ***</td>
<td valign="top" align="center">128.78 ***<break/>2019&gt;2020<break/>stress only</td>
<td valign="top" align="center">0.63</td>
<td valign="top" align="center">0.3</td>
</tr>
<tr>
<td valign="top" align="left">SPAD 2</td>
<td valign="top" align="left">
<break/>Range<break/>LSD</td>
<td valign="top" align="center">42.8<break/>29.2 - 57.4<break/>8.5</td>
<td valign="top" align="center">24.0<break/>11.5 - 51.3<break/>9.8</td>
<td valign="top" align="center">SPAD Units</td>
<td valign="top" align="center">22.57 ***</td>
<td valign="top" align="center">6780.40 ***</td>
<td valign="top" align="center">9.30 ***</td>
<td valign="top" align="center">208.32 ***<break/>2020&gt;2019</td>
<td valign="top" align="center">0.73</td>
<td valign="top" align="center">0.89</td>
</tr>
<tr>
<td valign="top" align="left">Difference SPAD<break/>(DiffSPAD)</td>
<td valign="top" align="left">=Senescence<break/>Range<break/>LSD</td>
<td valign="top" align="center">-5.0<break/>-15.1 - 5.2<break/>10.1</td>
<td valign="top" align="center">15.5<break/>-6.8 - 25.2<break/>11.0</td>
<td valign="top" align="center">SPAD Units</td>
<td valign="top" align="center">11..56 ***</td>
<td valign="top" align="center">6100.08 ***</td>
<td valign="top" align="center">8.30 ***</td>
<td valign="top" align="center">54.74 ***<break/>2019&gt;2020</td>
<td valign="top" align="center">0.38</td>
<td valign="top" align="center">0.86</td>
</tr>
<tr>
<td valign="top" align="left">End of flowering<break/>(EF)</td>
<td valign="top" align="left">
<break/>Range<break/>LSD</td>
<td valign="top" align="center">97,0<break/>85.5-118.6<break/>15.7</td>
<td valign="top" align="center">86.3<break/>82.0-100.3<break/>11.2</td>
<td valign="top" align="center">DAS</td>
<td valign="top" align="center">8.13 ***</td>
<td valign="top" align="center">1169.23***</td>
<td valign="top" align="center">2.33***</td>
<td valign="top" align="center">230.27***<break/>2019&gt;2020</td>
<td valign="top" align="center">0.73</td>
<td valign="top" align="center">0.44</td>
</tr>
<tr>
<td valign="top" align="left">Maturity<break/>(MAT)</td>
<td valign="top" align="left">
<break/>Range<break/>LSD</td>
<td valign="top" align="center">127.4<break/>118.0-143.0<break/>10.1</td>
<td valign="top" align="center">115.1<break/>105.9 - 126.5<break/>9.2</td>
<td valign="top" align="center">DAS</td>
<td valign="top" align="center">9.28 ***</td>
<td valign="top" align="center">2404.57 ***</td>
<td valign="top" align="center">1.51 **</td>
<td valign="top" align="center">54.41 ***<break/>2019&gt;2020 C<break/>2020&gt;2019 S</td>
<td valign="top" align="center">0.75</td>
<td valign="top" align="center">0.68</td>
</tr>
<tr>
<td valign="top" align="left">Plant height<break/>(PH)</td>
<td valign="top" align="left">
<break/>Range<break/>LSD</td>
<td valign="top" align="center">64.6<break/>42.1 - 102.3<break/>18.5</td>
<td valign="top" align="center">56.5<break/>37.1 - 71.6<break/>12.2</td>
<td valign="top" align="center">cm</td>
<td valign="top" align="center">17.19 ***</td>
<td valign="top" align="center">359.77 ***</td>
<td valign="top" align="center">3.19 ***</td>
<td valign="top" align="center">170.93 ***<break/>2020&gt;2019<break/>control only</td>
<td valign="top" align="center">0.87</td>
<td valign="top" align="center">0.82</td>
</tr>
<tr>
<td valign="top" align="left">Pods per plant<break/>(PP)<break/>&#xa0;</td>
<td valign="top" align="left">
<break/>Range<break/>LSD</td>
<td valign="top" align="center">10.8<break/>4.5-27.1<break/>6.9</td>
<td valign="top" align="center">6.4<break/>3.5 - 12.0<break/>3.2</td>
<td valign="top" align="center">number</td>
<td valign="top" align="center">19.11 ***</td>
<td valign="top" align="center">1124.49***</td>
<td valign="top" align="center">5.22***</td>
<td valign="top" align="center">1.59***<break/>2020&gt;2019<break/>control only</td>
<td valign="top" align="center">0.89</td>
<td valign="top" align="center">0.88</td>
</tr>
<tr>
<td valign="top" align="left">Seeds per plant<break/>(SP)<break/>&#xa0;</td>
<td valign="top" align="left">
<break/>Range<break/>LSD</td>
<td valign="top" align="center">24.9<break/>8.9 - 77.9<break/>17.3</td>
<td valign="top" align="center">14.2<break/>5.8 - 27.0<break/>7.5</td>
<td valign="top" align="center">number</td>
<td valign="top" align="center">21.3 ***</td>
<td valign="top" align="center">1039.57 ***</td>
<td valign="top" align="center">6.86 ***</td>
<td valign="top" align="center">1.60 ***<break/>2020&gt;2019</td>
<td valign="top" align="center">0.9</td>
<td valign="top" align="center">0.85</td>
</tr>
<tr>
<td valign="top" align="left">Hundred seed weight<break/>(HSW)</td>
<td valign="top" align="left">
<break/>Range<break/>LSD</td>
<td valign="top" align="center">66.93<break/>30.27 - 115.38<break/>23.67</td>
<td valign="top" align="center">60.44<break/>26.07 - 103.24<break/>16.28</td>
<td valign="top" align="center">g/100 seeds</td>
<td valign="top" align="center">28.6 ***</td>
<td valign="top" align="center">133.13 ***</td>
<td valign="top" align="center">1.39 **</td>
<td valign="top" align="center">420.52 ***<break/>2019&gt;2020</td>
<td valign="top" align="center">0.83</td>
<td valign="top" align="center">0.91</td>
</tr>
<tr>
<td valign="top" align="left">Plot yield<break/>(PY)<break/>&#xa0;</td>
<td valign="top" align="left">
<break/>Range<break/>LSD</td>
<td valign="top" align="center">0.155<break/>0.047 - 0.321<break/>0.119</td>
<td valign="top" align="center">0.083<break/>0.031 - 0.117<break/>0.054</td>
<td valign="top" align="center">kg/plot</td>
<td valign="top" align="center">9.19 ***</td>
<td valign="top" align="center">1003.33 ***</td>
<td valign="top" align="center">3.26 ***</td>
<td valign="top" align="center">5.0 *<break/>2019&gt;2020</td>
<td valign="top" align="center">0.83</td>
<td valign="top" align="center">0.65</td>
</tr>
<tr>
<td valign="top" align="left">Stress tolerance index<break/>(STI)</td>
<td valign="top" align="left">
<break/>Range<break/>LSD</td>
<td valign="top" colspan="2" align="center">0.564<break/>0.074 - 1.378<break/>0.048</td>
<td valign="top" align="center"/>
<td valign="top" align="center">4.94 ***</td>
<td valign="top" align="center"/>
<td valign="top" align="center"/>
<td valign="top" align="center">1.38 n.s.</td>
<td valign="top" colspan="2" align="center">0.80</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn>
<p>***, ** and * are significant for &#x3b1; &#x2264; 0.1, 1 and 5%, respectively, n.s., not significant; LSD, Least Significant Difference; DAS, Days after sowing.</p>
</fn>
</table-wrap-foot>
</table-wrap>
<p>The results obtained for the physiological traits showed clear differences between treatments. TSS and PRO were measured when drought stress was fully developed and senescence started in the stress treatment. As expected, both parameters increased significantly (24% and 269% respectively) under drought stress (<xref ref-type="fig" rid="f2">
<bold>Figure&#xa0;2</bold>
</xref>). The significant interaction between genotype x treatment indicates a different response of genotypes to drought stress (<xref ref-type="table" rid="T2">
<bold>Table&#xa0;2</bold>
</xref>). The leaf chlorophyll concentration, assessed as SPAD values, showed a marked decrease (44%) after 4 weeks of drought stress compared to the control (SPAD2 and DiffSPAD, <xref ref-type="fig" rid="f3">
<bold>Figure&#xa0;3</bold>
</xref>), which was clearly visible (<xref ref-type="fig" rid="f4">
<bold>Figure&#xa0;4</bold>
</xref>). This was confirmed by the difference between the two SPAD values (DiffSPAD = SPAD1 &#x2013; SPAD2) which minimizes the effects of variability in chlorophyll content among genotypes and stands for the degree of senescence. Whereas in the control treatment values were increasing until full flowering (negative difference), chlorophyll content was already decreasing in the stress treatment (positive difference).</p>
<fig id="f2" position="float">
<label>Figure&#xa0;2</label>
<caption>
<p>Boxplots for total content of soluble sugars (TSS) and free proline (PRO) for control (C) and drought stress (DS) treatment.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-14-1236147-g002.tif"/>
</fig>
<fig id="f3" position="float">
<label>Figure&#xa0;3</label>
<caption>
<p>Chlorophyll content by SPAD Meter about 4 weeks after onset of drought stress (SPAD 2) and difference to SPAD 1 values (DiffSPAD) for control (C) and drought stress (DS) treatment.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-14-1236147-g003.tif"/>
</fig>
<fig id="f4" position="float">
<label>Figure&#xa0;4</label>
<caption>
<p>Visible signs of early senescence under drought stress in rain out-shelter (left side) compared to the control (right side) at the beginning of July 2019.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-14-1236147-g004.tif"/>
</fig>
<p>For all the phenological and morphological traits, the mean values also decreased under drought stress conditions with values ranging from 9.7% for MAT to 11.3% for EF (<xref ref-type="table" rid="T2">
<bold>Table&#xa0;2</bold>
</xref>). PH decreased by 12.5% and the taller forms were more affected than the shorter types (data not shown). In contrast, strict declines were observed for the remaining yield related traits PP, SP and PY. Mean PY value under drought stress was significantly reduced by 46.4%, while PP and SP were reduced by 40.8% and 42.8%, respectively. Similarly, HSW values were lower, compared to the control, but the mean decrease of 9.7% was not as high as that of the other yield components (<xref ref-type="fig" rid="f5">
<bold>Figure&#xa0;5</bold>
</xref>). A significant interaction between genotypes and drought stress treatment was detected.</p>
<fig id="f5" position="float">
<label>Figure&#xa0;5</label>
<caption>
<p>Boxplots for seeds per plant (SP), hundred seed weight (HSW) and plot yield (PY) for control (C) and drought stress (DS) treatment.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-14-1236147-g005.tif"/>
</fig>
<p>Heritability computed across the two water regimes (control vs. drought) ranged from 0.3 (PRO) to 0.91 (HSW), both in stress and non-stress conditions (<xref ref-type="table" rid="T2">
<bold>Table&#xa0;2</bold>
</xref>). Heritability was similarly high (&gt; 0.82) under both conditions for the yield related traits PH, PP, SP and HSW. Moderate to high heritability was recorded for the physiological traits TSS (0.56-0.66), MAT (0.75-0.68), SPAD2 (0.73-0.89), DiffSPAD (0.38-0.86) and PY (0.83-0.65) in both conditions. The most contrasting heritability values between water regimes were found in EF (0.73-0.44) and PRO (0.63-0.3).</p>
</sec>
<sec id="s3_2">
<title>Yield relationships and genotypes</title>
<p>For a better understanding of the relationship between yield under control and stress conditions, a regression analysis was conducted. This relationship provides a benchmark for identifying the most drought tolerant genotypes (<xref ref-type="fig" rid="f6">
<bold>Figure&#xa0;6</bold>
</xref>). A positive correlation was detected between the two parameters (r = 0.75), whereby the high yielding genotypes were more affected by drought stress than the low yielding ones.The non-linear regression line shows that stress yield may not exceed a distinct level predefined by the stress environment. A correlation coefficient of -0.72 (<xref ref-type="fig" rid="f7">
<bold>Figure&#xa0;7</bold>
</xref>) revealed the negative and relatively close correlation between the yield under controlled conditions (standing for yield potential) and the relative yield (standing for yield stability). The figure clearly depicts that selection of genotypes with a high relative yield in most cases leads to rather low yielding genotypes. For this reason, a stress tolerance index (STI) based on seed yield (PY) was calculated for each genotype in response to severe drought and well-watered conditions, in order to identify genotypes producing high yields under both conditions. The variability of this index is shown in <xref ref-type="fig" rid="f8">
<bold>Figure&#xa0;8</bold>
</xref>. By using STI as target parameter, it was possible to identify genotypes that combine a high yield potential with a reasonable stress yield in this faba bean assortment. Accession ID and origin of the eight most drought-tolerant genotypes (5, 26, 35, 49, 56, 64, 75 and 85), according to STI are given in <xref ref-type="table" rid="T3">
<bold>Table&#xa0;3</bold>
</xref>.</p>
<fig id="f6" position="float">
<label>Figure&#xa0;6</label>
<caption>
<p>Relation between yield under drought stress (PY_DS) and yield in the control (PY_C); the solid black line stands for 100% yield stability.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-14-1236147-g006.tif"/>
</fig>
<fig id="f7" position="float">
<label>Figure&#xa0;7</label>
<caption>
<p>Relationship between relative yield and yield in the control (PY_C).</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-14-1236147-g007.tif"/>
</fig>
<fig id="f8" position="float">
<label>Figure&#xa0;8</label>
<caption>
<p>Genotypes of the test set ordered according to their STI value. PY_C: control yield. PY_DS: drought stress yield.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-14-1236147-g008.tif"/>
</fig>
<table-wrap id="T3" position="float">
<label>Table&#xa0;3</label>
<caption>
<p>STI value, ID and origin of the eight best genotypes.</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="middle" align="center">STI value</th>
<th valign="middle" align="center">Genotype</th>
<th valign="middle" align="center">Accession ID</th>
<th valign="middle" align="center">Accession</th>
<th valign="middle" align="center">Origin</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="middle" align="center">1.371</td>
<td valign="middle" align="center">35</td>
<td valign="middle" align="center">EUC_VF_077</td>
<td valign="middle" align="center">1272-1</td>
<td valign="middle" align="center">RUS</td>
</tr>
<tr>
<td valign="middle" align="center">1.330</td>
<td valign="middle" align="center">49</td>
<td valign="middle" align="center">EUC_VF_130</td>
<td valign="middle" align="center">INRA 2394</td>
<td valign="middle" align="center">DEU</td>
</tr>
<tr>
<td valign="middle" align="center">1.323</td>
<td valign="middle" align="center">64</td>
<td valign="middle" align="center">EUC_VF_194</td>
<td valign="middle" align="center">INRA 612</td>
<td valign="middle" align="center">EGY</td>
</tr>
<tr>
<td valign="middle" align="center">1.176</td>
<td valign="middle" align="center">5</td>
<td valign="middle" align="center">EUC_VF_ 009</td>
<td valign="middle" align="center">359-1</td>
<td valign="middle" align="center">TUR</td>
</tr>
<tr>
<td valign="middle" align="center">1.116</td>
<td valign="middle" align="center">75</td>
<td valign="middle" align="center">EUC_VF_ 303</td>
<td valign="middle" align="center">IG 72242</td>
<td valign="middle" align="center">CHN</td>
</tr>
<tr>
<td valign="middle" align="center">1.068</td>
<td valign="middle" align="center">85</td>
<td valign="middle" align="center">EUC_VF_336</td>
<td valign="middle" align="center">Misr 3</td>
<td valign="middle" align="center">EGY</td>
</tr>
<tr>
<td valign="middle" align="center">1.001</td>
<td valign="middle" align="center">56</td>
<td valign="middle" align="center">EUC_VF_174</td>
<td valign="middle" align="center">1248-3</td>
<td valign="middle" align="center">IRN</td>
</tr>
<tr>
<td valign="middle" align="center">0.987</td>
<td valign="middle" align="center">26</td>
<td valign="middle" align="center">EUC_VF_064</td>
<td valign="middle" align="center">1157-1</td>
<td valign="middle" align="center">ESP</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>The K-means cluster analysis using the morphological, phenological and physiological faba bean data and the STI index in control and drought conditions classified the accessions into groups based on their dissimilarities (<xref ref-type="fig" rid="f9">
<bold>Figure&#xa0;9</bold>
</xref>). For control conditions, two clusters were identified: cluster 1 (15 genotypes) showed highest STI value and highest average values for all traits (except for HSW), containing six of the eight best performing STI genotypes (5, 35, 56, 64, 75 and 85), while cluster 2 (84 genotypes) had the maximum average value for HSW and contained the remaining genotypes 26 and 49 (<xref ref-type="table" rid="T4">
<bold>Table&#xa0;4</bold>
</xref>). Under drought conditions, the genotypes were classified in three clusters: cluster 1 (56 genotypes), cluster 2 (13 genotypes) and cluster 3 (30 genotypes). Cluster 2 showed the best average value for PRO, SPAD2 (and correspondingly low DiffSPAD), MAT, PH, PP, SP, PY and STI, while cluster 1 had a similar yield but an intermediate value for STI and cluster 3 had the lower yield and STI mean values but the highest ones for TSS, EF and HWS (<xref ref-type="table" rid="T4">
<bold>Table&#xa0;4</bold>
</xref>). Again, six of the best STI genotypes (accessions 5, 35, 49, 56, 64 and 85) were included in cluster 2 while cluster 1 contained the other two accessions (26 and 75). According to cluster analysis, genotypes from cluster 2 were superior regarding the STI index and yield related traits, so cultivation of these genotypes would be recommended under drought stress conditions.</p>
<fig id="f9" position="float">
<label>Figure&#xa0;9</label>
<caption>
<p>Optimum number of clusters and clustering of the 100 faba bean accessions under control <bold>(A)</bold> and drought <bold>(B)</bold> conditions.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-14-1236147-g009.tif"/>
</fig>
<table-wrap id="T4" position="float">
<label>Table&#xa0;4</label>
<caption>
<p>Average value based on physiological, phenological and morphological and the STI index under control (C) and drought stress conditions (DS) after K-means cluster analysis.</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="top" align="left">Cluster</th>
<th valign="top" align="center">TSS_C</th>
<th valign="top" align="center">PRO_C</th>
<th valign="top" align="center">SPAD2_C</th>
<th valign="top" align="center">DiffSPAD_C</th>
<th valign="top" align="center">EF_C</th>
<th valign="top" align="center">MAT_C</th>
<th valign="top" align="center">PH_C</th>
<th valign="top" align="center">PP_C</th>
<th valign="top" align="center">SP_C</th>
<th valign="top" align="center">HSW_C</th>
<th valign="top" align="center">PY_C</th>
<th valign="top" align="center">STI</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="top" align="left">
<bold>1 (15)*</bold>
</td>
<td valign="top" align="center">1212.53</td>
<td valign="top" align="center">2.95</td>
<td valign="top" align="center">44.66</td>
<td valign="top" align="center">-5.51</td>
<td valign="top" align="center">102.27</td>
<td valign="top" align="center">132.62</td>
<td valign="top" align="center">84.89</td>
<td valign="top" align="center">17.85</td>
<td valign="top" align="center">44.33</td>
<td valign="top" align="center">51.04</td>
<td valign="top" align="center">0.21</td>
<td valign="top" align="center">0.82</td>
</tr>
<tr>
<td valign="top" align="left">
<bold>2 (84)</bold>
</td>
<td valign="top" align="center">1099.55</td>
<td valign="top" align="center">2.13</td>
<td valign="top" align="center">42.45</td>
<td valign="top" align="center">-4.90</td>
<td valign="top" align="center">96.01</td>
<td valign="top" align="center">126.39</td>
<td valign="top" align="center">60.48</td>
<td valign="top" align="center">9.46</td>
<td valign="top" align="center">21.28</td>
<td valign="top" align="center">69.97</td>
<td valign="top" align="center">0.14</td>
<td valign="top" align="center">0.51</td>
</tr>
<tr>
<th valign="top" align="left">Cluster</th>
<th valign="top" align="center">TSS_DS</th>
<th valign="top" align="center">PRO_DS</th>
<th valign="top" align="center">SPAD2_DS</th>
<th valign="top" align="center">DiffSPAD_DS</th>
<th valign="top" align="center">EF_DS</th>
<th valign="top" align="center">MAT_DS</th>
<th valign="top" align="center">PH_DS</th>
<th valign="top" align="center">PP_DS</th>
<th valign="top" align="center">SP_DS</th>
<th valign="top" align="center">HSW_DS</th>
<th valign="top" align="center">PY_DS</th>
<th valign="top" align="center">STI</th>
</tr>

<tr>
<td valign="top" align="left">
<bold>1 (56)</bold>
</td>
<td valign="top" align="center">1276.416</td>
<td valign="top" align="center">6.671</td>
<td valign="top" align="center">21.529</td>
<td valign="top" align="center">18.183</td>
<td valign="top" align="center">84.221</td>
<td valign="top" align="center">112.875</td>
<td valign="top" align="center">54.782</td>
<td valign="top" align="center">6.405</td>
<td valign="top" align="center">14.869</td>
<td valign="top" align="center">62.408</td>
<td valign="top" align="center">0.091</td>
<td valign="top" align="center">0.612</td>
</tr>
<tr>
<td valign="top" align="left">
<bold>2 (13)</bold>
</td>
<td valign="top" align="center">1310.345</td>
<td valign="top" align="center">17.398</td>
<td valign="top" align="center">37.205</td>
<td valign="top" align="center">2.385</td>
<td valign="top" align="center">88.485</td>
<td valign="top" align="center">119.037</td>
<td valign="top" align="center">65.862</td>
<td valign="top" align="center">9.148</td>
<td valign="top" align="center">20.596</td>
<td valign="top" align="center">45.726</td>
<td valign="top" align="center">0.092</td>
<td valign="top" align="center">0.899</td>
</tr>
<tr>
<td valign="top" align="left">
<bold>3 (30)</bold>
</td>
<td valign="top" align="center">1634.117</td>
<td valign="top" align="center">6.995</td>
<td valign="top" align="center">22.449</td>
<td valign="top" align="center">16.267</td>
<td valign="top" align="center">88.961</td>
<td valign="top" align="center">117.319</td>
<td valign="top" align="center">55.341</td>
<td valign="top" align="center">5.157</td>
<td valign="top" align="center">10.324</td>
<td valign="top" align="center">63.628</td>
<td valign="top" align="center">0.065</td>
<td valign="top" align="center">0.321</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn>
<p>* Number of genotypes in each cluster.</p>
</fn>
</table-wrap-foot>
</table-wrap>
</sec>
<sec id="s3_3">
<title>Correlations and principal component analysis (PCA)</title>
<p>A correlation analysis was performed to evaluate the relationships among morphological, phenological and physiological traits as well as their association with PY and STI. As mentioned above, the correlations of plot yield under control (PY_C) and drought stress conditions (PY_DS) with STI are similar and highly significant, with PY_C having the greater impact on STI (<xref ref-type="table" rid="T5">
<bold>Table&#xa0;5</bold>
</xref>). The most yield influencing factors were SP and PP, although physiological traits such as chlorophyll content (SPAD values), PRO and TSS together with phenological traits such as EF and MAT, also contributed significantly to the respective yield and to the STI values.</p>
<table-wrap id="T5" position="float">
<label>Table&#xa0;5</label>
<caption>
<p>Correlation matrix with Pearson&#x2019;s correlation coefficient r for all relevant traits under control (above) and drought stress conditions (down).</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="middle" align="left">Trait</th>
<th valign="middle" align="center">TSS _C</th>
<th valign="middle" align="center">PRO_C</th>
<th valign="middle" align="center">SPAD2_C</th>
<th valign="middle" align="center">DiffSPAD_C</th>
<th valign="middle" align="center">EF_C</th>
<th valign="middle" align="center">MAT_C</th>
<th valign="middle" align="center">PH_C</th>
<th valign="middle" align="center">PP_C</th>
<th valign="middle" align="center">SP_C</th>
<th valign="middle" align="center">HSW_C</th>
<th valign="middle" align="center">PY_C</th>
<th valign="middle" align="center">STI</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="middle" align="left">
<bold>TSS_C</bold>
</td>
<td valign="middle" align="center">1</td>
<td valign="middle" align="center">-0.097</td>
<td valign="middle" align="center">
<bold>-0.290</bold>
</td>
<td valign="middle" align="center">
<bold>0.482</bold>
</td>
<td valign="middle" align="center">
<bold>0.385</bold>
</td>
<td valign="middle" align="center">
<bold>0.276</bold>
</td>
<td valign="middle" align="center">
<bold>0.284</bold>
</td>
<td valign="middle" align="center">0.118</td>
<td valign="middle" align="center">0.026</td>
<td valign="middle" align="center">-0.219</td>
<td valign="middle" align="center">-0.210</td>
<td valign="middle" align="center">
<bold>-0.256</bold>
</td>
</tr>
<tr>
<td valign="middle" align="left">
<bold>PRO_C</bold>
</td>
<td valign="middle" align="center"/>
<td valign="middle" align="center">1</td>
<td valign="middle" align="center">
<bold>0.260</bold>
</td>
<td valign="middle" align="center">-0.164</td>
<td valign="middle" align="center">0.185</td>
<td valign="middle" align="center">
<bold>0.357</bold>
</td>
<td valign="middle" align="center">
<bold>0.560</bold>
</td>
<td valign="middle" align="center">
<bold>0.562</bold>
</td>
<td valign="middle" align="center">
<bold>0.616</bold>
</td>
<td valign="middle" align="center">
<bold>-0.358</bold>
</td>
<td valign="middle" align="center">
<bold>0.421</bold>
</td>
<td valign="middle" align="center">
<bold>0.352</bold>
</td>
</tr>
<tr>
<td valign="middle" align="left">
<bold>SPAD2_C</bold>
</td>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center">1</td>
<td valign="middle" align="center">
<bold>-0.812</bold>
</td>
<td valign="middle" align="center">0.035</td>
<td valign="middle" align="center">
<bold>0.372</bold>
</td>
<td valign="middle" align="center">
<bold>0.263</bold>
</td>
<td valign="middle" align="center">
<bold>0.310</bold>
</td>
<td valign="middle" align="center">
<bold>0.345</bold>
</td>
<td valign="middle" align="center">0.129</td>
<td valign="middle" align="center">
<bold>0.509</bold>
</td>
<td valign="middle" align="center">
<bold>0.474</bold>
</td>
</tr>
<tr>
<td valign="middle" align="left">
<bold>DiffSPAD_C</bold>
</td>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center">1</td>
<td valign="middle" align="center">-0.004</td>
<td valign="middle" align="center">
<bold>-0.259</bold>
</td>
<td valign="middle" align="center">-0.163</td>
<td valign="middle" align="center">
<bold>-0.261</bold>
</td>
<td valign="middle" align="center">
<bold>-0.267</bold>
</td>
<td valign="middle" align="center">-0.200</td>
<td valign="middle" align="center">
<bold>-0.474</bold>
</td>
<td valign="middle" align="center">
<bold>-0.448</bold>
</td>
</tr>
<tr>
<td valign="middle" align="left">
<bold>EF_C</bold>
</td>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center">1</td>
<td valign="middle" align="center">
<bold>0.694</bold>
</td>
<td valign="middle" align="center">
<bold>0.408</bold>
</td>
<td valign="middle" align="center">
<bold>0.255</bold>
</td>
<td valign="middle" align="center">0.132</td>
<td valign="middle" align="center">-0.106</td>
<td valign="middle" align="center">-0.017</td>
<td valign="middle" align="center">-0.105</td>
</tr>
<tr>
<td valign="middle" align="left">
<bold>MAT_C</bold>
</td>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center">1</td>
<td valign="middle" align="center">
<bold>0.619</bold>
</td>
<td valign="middle" align="center">
<bold>0.444</bold>
</td>
<td valign="middle" align="center">
<bold>0.377</bold>
</td>
<td valign="middle" align="center">-0.007</td>
<td valign="middle" align="center">
<bold>0.374</bold>
</td>
<td valign="middle" align="center">
<bold>0.275</bold>
</td>
</tr>
<tr>
<td valign="middle" align="left">
<bold>PH_C</bold>
</td>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center">1</td>
<td valign="middle" align="center">
<bold>0.783</bold>
</td>
<td valign="middle" align="center">
<bold>0.767</bold>
</td>
<td valign="middle" align="center">
<bold>-0.340</bold>
</td>
<td valign="middle" align="center">
<bold>0.590</bold>
</td>
<td valign="middle" align="center">
<bold>0.484</bold>
</td>
</tr>
<tr>
<td valign="middle" align="left">
<bold>PP_C</bold>
</td>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center">1</td>
<td valign="middle" align="center">
<bold>0.925</bold>
</td>
<td valign="middle" align="center">
<bold>-0.515</bold>
</td>
<td valign="middle" align="center">
<bold>0.663</bold>
</td>
<td valign="middle" align="center">
<bold>0.566</bold>
</td>
</tr>
<tr>
<td valign="middle" align="left">
<bold>SP_C</bold>
</td>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center">1</td>
<td valign="middle" align="center">
<bold>-0.499</bold>
</td>
<td valign="middle" align="center">
<bold>0.737</bold>
</td>
<td valign="middle" align="center">
<bold>0.650</bold>
</td>
</tr>
<tr>
<td valign="middle" align="left">
<bold>HSW_C</bold>
</td>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center">1</td>
<td valign="middle" align="center">0.153</td>
<td valign="middle" align="center">0.201</td>
</tr>
<tr>
<td valign="middle" align="left">
<bold>PY_C</bold>
</td>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center">1</td>
<td valign="middle" align="center">
<bold>0.961</bold>
</td>
</tr>
<tr>
<td valign="middle" align="left">
<bold>Trait</bold>
</td>
<td valign="middle" align="center">
<bold>TSS _DS</bold>
</td>
<td valign="middle" align="center">
<bold>PRO_DS</bold>
</td>
<td valign="middle" align="center">
<bold>SPAD2_DS</bold>
</td>
<td valign="middle" align="center">
<bold>DiffSPAD_DS</bold>
</td>
<td valign="middle" align="center">
<bold>EF_DS</bold>
</td>
<td valign="middle" align="center">
<bold>MAT_DS</bold>
</td>
<td valign="middle" align="center">
<bold>PH_DS</bold>
</td>
<td valign="middle" align="center">
<bold>PP_DS</bold>
</td>
<td valign="middle" align="center">
<bold>SP_DS</bold>
</td>
<td valign="middle" align="center">
<bold>HSW_DS</bold>
</td>
<td valign="middle" align="center">
<bold>PY_DS</bold>
</td>
<td valign="middle" align="center">
<bold>STI</bold>
</td>
</tr>
<tr>
<td valign="middle" align="left">
<bold>TSS_DS</bold>
</td>
<td valign="middle" align="center">1</td>
<td valign="middle" align="center">0.217</td>
<td valign="middle" align="center">0.197</td>
<td valign="middle" align="center">0.187</td>
<td valign="middle" align="center">
<bold>0.267</bold>
</td>
<td valign="middle" align="center">
<bold>0.333</bold>
</td>
<td valign="middle" align="center">0.035</td>
<td valign="middle" align="center">
<bold>-0.293</bold>
</td>
<td valign="middle" align="center">
<bold>-0.414</bold>
</td>
<td valign="middle" align="center">0.123</td>
<td valign="middle" align="center">
<bold>-0.441</bold>
</td>
<td valign="middle" align="center">
<bold>-0.410</bold>
</td>
</tr>
<tr>
<td valign="middle" align="left">
<bold>PRO_DS</bold>
</td>
<td valign="middle" align="center"/>
<td valign="middle" align="center">1</td>
<td valign="middle" align="center">
<bold>0.510</bold>
</td>
<td valign="middle" align="center">
<bold>-0.459</bold>
</td>
<td valign="middle" align="center">0.100</td>
<td valign="middle" align="center">
<bold>0.343</bold>
</td>
<td valign="middle" align="center">
<bold>0.279</bold>
</td>
<td valign="middle" align="center">0.171</td>
<td valign="middle" align="center">
<bold>0.254</bold>
</td>
<td valign="middle" align="center">-0.052</td>
<td valign="middle" align="center">0.187</td>
<td valign="middle" align="center">
<bold>0.369</bold>
</td>
</tr>
<tr>
<td valign="middle" align="left">
<bold>SPAD2_DS</bold>
</td>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center">1</td>
<td valign="middle" align="center">
<bold>-0.905</bold>
</td>
<td valign="middle" align="center">
<bold>0.282</bold>
</td>
<td valign="middle" align="center">
<bold>0.558</bold>
</td>
<td valign="middle" align="center">
<bold>0.442</bold>
</td>
<td valign="middle" align="center">
<bold>0.423</bold>
</td>
<td valign="middle" align="center">
<bold>0.478</bold>
</td>
<td valign="middle" align="center">
<bold>-0.325</bold>
</td>
<td valign="middle" align="center">0.137</td>
<td valign="middle" align="center">
<bold>0.373</bold>
</td>
</tr>
<tr>
<td valign="middle" align="left">
<bold>DiffSPAD_DS</bold>
</td>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center">1</td>
<td valign="middle" align="center">
<bold>-0.269</bold>
</td>
<td valign="middle" align="center">
<bold>-0.504</bold>
</td>
<td valign="middle" align="center">
<bold>-0.392</bold>
</td>
<td valign="middle" align="center">
<bold>-0.420</bold>
</td>
<td valign="middle" align="center">
<bold>-0.432</bold>
</td>
<td valign="middle" align="center">
<bold>0.327</bold>
</td>
<td valign="middle" align="center">-0.068</td>
<td valign="middle" align="center">
<bold>-0.317</bold>
</td>
</tr>
<tr>
<td valign="middle" align="left">
<bold>EF_DS</bold>
</td>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center">1</td>
<td valign="middle" align="center">
<bold>0.680</bold>
</td>
<td valign="middle" align="center">
<bold>0.253</bold>
</td>
<td valign="middle" align="center">-0.014</td>
<td valign="middle" align="center">-0.172</td>
<td valign="middle" align="center">-0.029</td>
<td valign="middle" align="center">
<bold>-0.405</bold>
</td>
<td valign="middle" align="center">
<bold>-0.250</bold>
</td>
</tr>
<tr>
<td valign="middle" align="left">
<bold>MAT_DS</bold>
</td>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center">1</td>
<td valign="middle" align="center">
<bold>0.393</bold>
</td>
<td valign="middle" align="center">0.076</td>
<td valign="middle" align="center">-0.090</td>
<td valign="middle" align="center">0.032</td>
<td valign="middle" align="center">-0.182</td>
<td valign="middle" align="center">0.023</td>
</tr>
<tr>
<td valign="middle" align="left">
<bold>PH_DS</bold>
</td>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center">1</td>
<td valign="middle" align="center">
<bold>0.500</bold>
</td>
<td valign="middle" align="center">
<bold>0.398</bold>
</td>
<td valign="middle" align="center">-0.195</td>
<td valign="middle" align="center">0.206</td>
<td valign="middle" align="center">
<bold>0.381</bold>
</td>
</tr>
<tr>
<td valign="middle" align="left">
<bold>PP_DS</bold>
</td>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center">1</td>
<td valign="middle" align="center">
<bold>0.819</bold>
</td>
<td valign="middle" align="center">
<bold>-0.626</bold>
</td>
<td valign="middle" align="center">
<bold>0.332</bold>
</td>
<td valign="middle" align="center">
<bold>0.429</bold>
</td>
</tr>
<tr>
<td valign="middle" align="left">
<bold>SP_DS</bold>
</td>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center">1</td>
<td valign="middle" align="center">
<bold>-0.651</bold>
</td>
<td valign="middle" align="center">
<bold>0.499</bold>
</td>
<td valign="middle" align="center">
<bold>0.523</bold>
</td>
</tr>
<tr>
<td valign="middle" align="left">
<bold>HSW_DS</bold>
</td>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center">1</td>
<td valign="middle" align="center">
<bold>0.268</bold>
</td>
<td valign="middle" align="center">0.129</td>
</tr>
<tr>
<td valign="middle" align="left">
<bold>PY_DS</bold>
</td>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center"/>
<td valign="middle" align="center">1</td>
<td valign="middle" align="center">
<bold>0.836</bold>
</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn>
<p>Numbers in bold are significant for &#x3b1; &#x2264; 5%. PRO, content of free proline; TSS, total content of soluble sugars; EF, end of flowering; MAT, maturity; PH, plant height; PP, pods per plant; SP, seeds per plant; HSW, hundred seeds weight; PY, plot yield; STI, stress tolerance index; C, control; DS, drought stress.</p>
</fn>
</table-wrap-foot>
</table-wrap>
<p>Under control conditions, PH was significantly correlated with most traits, particularly PP and SP, highlighting their contribution to yield. Strong positive correlations were also observed for PP and SP with STI, PY and PRO. Interestingly, HSW showed a strong negative correlation with PP and SP. Under drought stress conditions, DiffSPAD showed a significant negative correlation with all traits except HSW. SPAD2 showed strong positive correlation with MAT and PRO and lower but still significant correlation with PH, PP and SP, which showed positive inter-correlations. In control conditions, all measured yield-related traits (PY, PP, and SP) were associated positively except HSW, indicating that they could simultaneously improve drought tolerance at the expense of seed weight. Correlations were stronger in the control treatment in nearly all cases with the exception of EF and TSS (<xref ref-type="table" rid="T5">
<bold>Table&#xa0;5</bold>
</xref>).</p>
<p>In order to determine the most contributing traits to yield under both conditions, a principal component analysis (PCA) was performed (<xref ref-type="fig" rid="f10">
<bold>Figure&#xa0;10</bold>
</xref>). Results showed that the first and second component (PC1 and PC2) represented 63.9% of the variance observed in control conditions and 59.3% of the variance in the drought stress treatment. Biplots for control and drought conditions were constructed from PC1 (42.6% and 36% variation explained, respectively) and PC2 (21.3% and 23.3% variation explained, respectively) with the distribution of the 100 genotypes, the morphological, phenological and physiological traits and STI. The PCA matrix revealed that the genotypes assayed followed different strategies with respect to seed set, seed size or physiological parameters (<xref ref-type="fig" rid="f10">
<bold>Figure&#xa0;10</bold>
</xref>). In the control condition, the PC1 axis was primarily associated with SP and PP, followed by PH and PRO, while TSS and differences in chlorophyll content (DiffSPAD) were the most contributing factors for the dispersion of the genotypes along PC2. In the drought stress treatment, chlorophyll content (SPAD2), followed by SP and PP had the largest contribution to PC1, while EF and MAT mainly contributed to PC2. The nearly 180&#xb0; angle between SPAD2 and DiffSPAD indicated a strong negative association, as expected. PCA biplot analysis further demonstrates that the eight superior genotypes according to STI follow different strategies with respect to phenological, physiological and yield components. Thus, in drought conditions the genotypes 26, 35, 64 and 75 showed a higher number of pods (PP_DS), seeds per plant (SP_DS) and yield (PY_DS) as well as a higher STI value, but low total soluble sugar content (TSS_DS). On the other hand, genotypes 5, 49, 56 and 85 showed higher values for plant height (PH_DS) and the physiological traits SPAD2 and PRO, together with a low hundred seed weight (HSW_DS).</p>
<fig id="f10" position="float">
<label>Figure&#xa0;10</label>
<caption>
<p>PCA based on all Least Squares (LS) means of the 100 genotypes for the respective traits assessed under control (above) and drought stress conditions (down). Points with numbers are the eight best genotypes according to STI (see ). The lines indicate eigenvectors representing the strength (length of the vector) and the direction of the parameter correlation relative to the first two principal componente (PC1 and PC2). TSS, total content of soluble sugars; PRO, free proline content; SPAD, measures of chlorophyll content 4 weeks after onset of stress (SPAD 2) and differences in chlorophyll content (DiffSPAD); EF, end of flowering; MAT, maturity; PH, plant height; PP, pods per plant; SP, seeds per plant; HSW, hundred seed weight; PY, plot yield; STI, stress tolerance index; C, control, DS, drought stress.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fpls-14-1236147-g010.tif"/>
</fig>
</sec>
</sec>
<sec id="s4" sec-type="discussion">
<title>Discussion</title>
<p>In this study, a faba bean collection of 100 accessions was used to assess the relevance of morphological, phenological and physiological traits with respect to yield under different water regimes. One determining factor is the chlorophyll content, measured indirectly by SPAD metre. A decrease in chlorophyll under drought stress has been widely reported (<xref ref-type="bibr" rid="B10">Ammar et&#xa0;al., 2015</xref>; <xref ref-type="bibr" rid="B74">Siddiqui et&#xa0;al., 2015</xref>; <xref ref-type="bibr" rid="B2">Abid et&#xa0;al., 2017</xref>). Drought-induced early senescence leading to early chlorophyll degradation was also observed in the present study with genotypes responding differently to the stress factor. Late chlorophyll degradation in connection with prolonged photosynthetic activity is a desirable trait, especially under drought stress,</p>
<p>Accumulation of free proline and of carbohydrates as free sugars contributes to osmotic adjustment in plants (<xref ref-type="bibr" rid="B86">Yoshiba et&#xa0;al., 1997</xref>; <xref ref-type="bibr" rid="B33">Hekneby et&#xa0;al., 2006</xref>). An increase in free proline between 42% and 202%, under 30% of field capacity was reported in a faba bean greenhouse trial (<xref ref-type="bibr" rid="B2">Abid et&#xa0;al., 2017</xref>), These results evidenced that accumulation of proline has a role in the faba bean tolerance to drought. Our study further supports this idea, as the amount of proline was four times higher under drought stress than in controlled conditions. <xref ref-type="bibr" rid="B74">Siddiqui et&#xa0;al. (2015)</xref> reported that drought-tolerant faba bean genotypes accumulated more proline than sensitive genotypes. In contrast, <xref ref-type="bibr" rid="B48">Migdadi et&#xa0;al. (2016)</xref> found a negative correlation between proline and yield over all treatments, which was mainly due to the high accumulation of proline under drought stress conditions rather than to differences between genotypes within the treatment. This is in accordance with our study, in which there was no close correlation between proline accumulation and other traits scored in the drought stress treatment. The significant correlation in the control treatment (PY_C) could be due to limitations in the measurements, as the content of free proline and therefore the genotypic differences are very low in the absence of osmotic stress (<xref ref-type="bibr" rid="B15">Balko, 2005</xref>; <xref ref-type="bibr" rid="B76">Stoddard et&#xa0;al., 2006</xref>; <xref ref-type="bibr" rid="B35">Kabbadj et&#xa0;al., 2017</xref>)</p>
<p>Another adaptation to water stress is the accumulation of carbohydrates, one of the main components of osmoregulation in many plant species (<xref ref-type="bibr" rid="B50">Morgan, 1984</xref>). In the forage legume <italic>Lotus japonicus</italic> L., accumulation of sugars such as fructose, galactose, glucose and maltose occurred when affected by drought (<xref ref-type="bibr" rid="B67">Sanchez et&#xa0;al., 2012</xref>). In soybean leaves, about 30% of the metabolites identified were soluble sugars and sugar alcohols, which are important for plant adaptation to stresses, especially during water deficit (<xref ref-type="bibr" rid="B18">Benkeblia et&#xa0;al., 2007</xref>). These findings are supported by the present study where a 24% increase of TSS under drought conditions was observed.</p>
<p>Plant height in faba bean is strongly determined by the growth type. Indeterminate types reach a greater height under favorable conditions while determinate (topless) types are lower. In this study, the indeterminate types showed a broad variability in plant height. Drought stress reduced plant height and this effect was more pronounced in the tall forms. Importantly, the yield-deciding parameter is not plant height itself but the number of reproductive nodes. Drought during the reproductive phase may reduce plant height and therefore the number of reproductive nodes (<xref ref-type="bibr" rid="B27">Gnanasambandam et&#xa0;al., 2012</xref>). In our field trials, plant height in both treatments was most significantly correlated to the number of pods per plant as well as to seeds per plant.</p>
<p>Pods per plant, seeds per pod and seed weight are among the main yield components (<xref ref-type="bibr" rid="B14">Ayaz et&#xa0;al., 2004</xref>; <xref ref-type="bibr" rid="B78">Tadesse et&#xa0;al., 2011</xref>; <xref ref-type="bibr" rid="B57">Ouji et&#xa0;al., 2017</xref>). Under drought stress, the most susceptible developmental stage is the reproductive phase including flowering, early podding and pod setting (<xref ref-type="bibr" rid="B51">Muktadir et&#xa0;al., 2020</xref>). The early podding was the most sensitive stage reported by <xref ref-type="bibr" rid="B54">Mwanamwenge et&#xa0;al. (1999)</xref>, causing a reduction in the number of pods and seeds due to high rates of abscission (<xref ref-type="bibr" rid="B69">Sekara et&#xa0;al., 2001</xref>). In our study, seeds per plant showed the closest relation to seed yield in both the control and the drought stress treatment, followed by pods per plant. Similar outcomes were reported by <xref ref-type="bibr" rid="B58">Pilbeam et&#xa0;al. (1992)</xref> and <xref ref-type="bibr" rid="B79">Tofiq et&#xa0;al. (2016)</xref>. <xref ref-type="bibr" rid="B43">Lopez et&#xa0;al. (1996)</xref> reported that the number of pods per plant was most affected by drought stress during flowering, with yield reductions up to 70%.</p>
<sec id="s4_1">
<title>Yield relations and genotypes</title>
<p>Previous studies pointed out that high yield potential under optimal conditions does not necessarily correlate with yield under limited water irrigation (<xref ref-type="bibr" rid="B60">Ramirez-Vallejo and Kelly, 1998</xref>; <xref ref-type="bibr" rid="B1">Abdelmula et&#xa0;al., 1999</xref>; <xref ref-type="bibr" rid="B21">El-Hendawy et&#xa0;al., 2017</xref>). This is in agreement with the results of our study where the performance of the faba bean lines was not consistent across the two treatments. The negative correlation between yield potential (yield under non-stress) and yield stability (yield under severe stress) suggests that the use of yield stability as a target parameter for drought tolerance could lead to the selection of low yielding genotypes. Nevertheless, varieties with high yield potential will often have an advantage over varieties with lower yield potential under moderate drought stress. Whether genotypes with medium yield potential and high stress yield might have an advantage for regions with regularly occurring drought stress remains a matter of debate. For regions with occasional drought stress, the relatively close correlation between stress yield and control yield (<xref ref-type="fig" rid="f6">
<bold>Figure&#xa0;6</bold>
</xref>) implies that both yield under well watered and drought stress should be considered, as shown by the STI value (<xref ref-type="bibr" rid="B26">Fernandez, 1992</xref>). Several studies point towards STI as the most suitable index for selecting the best yielding genotypes under contrasting conditions (<xref ref-type="bibr" rid="B26">Fernandez, 1992</xref>; <xref ref-type="bibr" rid="B21">El-Hendawy et&#xa0;al., 2017</xref>; <xref ref-type="bibr" rid="B47">Memari et&#xa0;al., 2022</xref>). According to our results, STI exhibits a strong correlation with yield (PY) and yield components (PP and SP) both under control conditions and drought stress, and appears to be an effective index for the selection of genotypes with high yield potential under both environmental conditions.</p>
</sec>
<sec id="s4_2">
<title>Correlations, cluster analysis and principal component analysis</title>
<p>Genotypes with higher yields in favorable as well as stress environments are the ones preferably selected in breeding programs. Nevertheless, the presence of genotype by environment (GxE) interactions is a major concern, since it may reduce correlation between genotypic and phenotypic values and slow down the selection progress (<xref ref-type="bibr" rid="B66">Romagosa et&#xa0;al., 2009</xref>).</p>
<p>Our results revealed a high GxE interaction for all the traits studied. Nevertheless, a clear association between the yield related traits, PP and SP and the physiological trait SPAD2 was observed through their significant correlation in both control and drought stress conditions. A moderate correlation was also found between PRO and SP in both environments. The high heritability estimates and the positive correlations of SPAD2 with PP and SP in stress conditions indicates that SPAD2 may be a suitable selection criterion for drought tolerance in faba bean.</p>
<p>A positive correlation between maturity and yield was found in the control condition, but not in the stress treatment. <xref ref-type="bibr" rid="B36">Katsoulieri et&#xa0;al. (2020)</xref> described that later ripening genotypes showed greater susceptibility to drought stress during early pod set compared to genotypes with early maturity. This was supported by <xref ref-type="bibr" rid="B39">Khan et&#xa0;al. (2010)</xref> as well as <xref ref-type="bibr" rid="B44">Manning et&#xa0;al. (2020)</xref>, where early flowering is considered an important adaptation trait for cool-season legumes growing under semi-arid conditions.</p>
<p>PCA is a powerful strategy to uncover genetic factors that contribute to complex traits. It has been considered a useful complementary tool in drought tolerance screenings and in the selection of the most tolerant genotypes (<xref ref-type="bibr" rid="B49">Mohamamdi et&#xa0;al., 2017</xref>; <xref ref-type="bibr" rid="B68">Santos et&#xa0;al., 2020</xref>; <xref ref-type="bibr" rid="B63">Reyes et&#xa0;al., 2022</xref>; <xref ref-type="bibr" rid="B84">Weng et&#xa0;al., 2022</xref>). This approach has been successfully adopted in drought stress tolerance evaluations in different crop legumes such as soybean (<xref ref-type="bibr" rid="B19">Bouslama and Schapaugh, 1984</xref>; <xref ref-type="bibr" rid="B81">Toum et&#xa0;al., 2022</xref>), lentil (<xref ref-type="bibr" rid="B73">Siahsar et&#xa0;al., 2010</xref>), chickpea (<xref ref-type="bibr" rid="B70">Shah et&#xa0;al., 2020</xref>; <xref ref-type="bibr" rid="B12">Arif et&#xa0;al., 2021</xref>) and faba bean (<xref ref-type="bibr" rid="B45">Mansour et&#xa0;al., 2021</xref>; <xref ref-type="bibr" rid="B4">Afzal et&#xa0;al., 2022</xref>). Here we applied the same strategy to integrate multiple parameters for grouping the faba bean accessions according to their drought tolerance. Our results show that in controlled conditions, the morphological traits SP, PP and PH were the most contributing factors for the dispersion of the genotypes along the PC1 axis, while PC2 was positively loaded with the TSS in the leaf, which is often increased during grain filling and may point towards late senescence.</p>
<p>Under drought stress conditions, there was a shift towards traits related to a longer photosynthesis period in PC1 (as late chlorophyll content or SPAD2), resulting in a higher number of seeds per plant. On the other hand, the phenological traits EF and MAT, which are strongly related to the length of the vegetation period, were the main factors contributing to PC2. These findings were supported by investigations in <italic>Miscanthus</italic>, where a PCA performed with relative values and the physiological traits chlorophyll content and chlorophyll fluorescence revealed their important role in drought stress response (<xref ref-type="bibr" rid="B84">Weng et&#xa0;al., 2022</xref>). The relevance of both physiological traits for final yield under drought stress was also reported in durum wheat (<xref ref-type="bibr" rid="B49">Mohamamdi et&#xa0;al., 2017</xref>; <xref ref-type="bibr" rid="B63">Reyes et&#xa0;al., 2022</xref>), further evidencing that a delay of chlorophyll degradation during senescence can improve performance under drought conditions. These observations reinforce the idea that screening for chlorophyll content under drought conditions could significantly contribute to efficient selection in breeding programs. The PCA was also helpful to group the faba bean genotypes according to their phenotypic response (<xref ref-type="bibr" rid="B3">Abou-Zaitoun et&#xa0;al., 2018</xref>; <xref ref-type="bibr" rid="B56">Or-Rashid et&#xa0;al., 2021</xref>). It further evidences that the lines tested follow different strategies to reach high yields under drought stress.</p>
<p>Cluster analysis is a multivariate analysis that categorizes the genotypes into several subgroups according to their similarities, based on the capability and performance of the traits examined. This approach has been used to predict the best-performing genotypes under drought stress in other crops (<xref ref-type="bibr" rid="B31">Hannok et&#xa0;al., 2021</xref>; <xref ref-type="bibr" rid="B53">Mutlu et&#xa0;al., 2022</xref>; <xref ref-type="bibr" rid="B13">Aslam et&#xa0;al., 2023</xref>). Grouping breeding materials using multivariate data under a given condition allows plant breeders to select breeding lines more efficiently. Here we applied the same strategy to integrate phenological, physiological as well as yield components for grouping the faba bean accessions according to their drought tolerance and the clustering method was effective to differentiate the best performing genotypes under drought stress.</p>
<p>PCA biplot and cluster analysis demonstrate that the eight superior genotypes according to STI follow different strategies with respect to phenological, physiological and yield components. Thus, in drought conditions the genotypes of cluster 2, that included six of the best performing STI genotypes (5, 35, 64, 56, 64 and 85), were the best group in terms of performance. These genotypes are characterized not only by the highest yield (PY_DS) and STI values but also by the highest mean values for proline (PRO_DS), chlorophyll content (SPAD2_DS), plant height (PH_DS), number of pods (PP_DS), seeds per plant (SP_DS) and later maturity days (MAT_DS), thus being the most promising traits for selection. On the other hand, these genotypes displayed lower total soluble sugar content (TSS_DS) and hundred seed weight values (HSW).</p>
<p>In summary, our results revealed that the Stress Tolerance Index (STI) is a useful criterion for selecting drought tolerant and high yielding faba bean genotypes. Using STI and biplot analysis, eight accessions with relatively high yield under both normal and drought stress conditions were identified and six of them further validated by K-means cluster analyses. Thus, genotypes 5, 35, 64, 56, 64 and 85 could potentially be used as genetic resources in faba bean breeding programs to develop varieties with enhanced drought resistance traits. Our study also suggests that SPAD2, followed by SP and PP can be reliable indicators for selection. To bridge the gap between traditional and molecular breeding, this faba bean collection has recently been genotyped with a high-density SNP genotyping array to conduct a genome-wide association (GWAS) study for drought resistance (<xref ref-type="bibr" rid="B29">Gutierrez et&#xa0;al., 2023</xref>). After validation, the candidate genes identified can be used for marker-assisted selection to accelerate and improve faba bean yield in agricultural areas where long water deficit periods are expected.</p>
</sec>
</sec>
<sec id="s5" sec-type="data-availability">
<title>Data availability statement</title>
<p>The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.</p>
</sec>
<sec id="s6" sec-type="author-contributions">
<title>Author contributions</title>
<p>AMT and NG selected and multiplied the seeds for the experiments. CB planned and carried out the drought evaluations and processed the experimental data. CB and AT wrote the manuscript and NG contributed in the review and editing. CB and AMT acquired the financial support for the project. All authors contributed to the article and approved the submitted version.</p>
</sec>
</body>
<back>
<sec id="s7" sec-type="funding-information">
<title>Funding</title>
<p>This project has received funding from the European Union&#x2019;s Horizon 2020 Programme for Research &amp; Innovation under grant agreement n&#xb0; 727312 (EUCLEG Project) and PID2020-114952RR-I00, co-financed by ERDF.</p>
</sec>
<ack>
<title>Acknowledgments</title>
<p>Many thanks to the seed companies and institutes which provided us with faba bean seeds of their accessions: INRA, Dijon (Jean-Bernard Magnin-Robert), ICARDA (Fouad Maalouf), and IFAPA.</p>
</ack>
<sec id="s8" sec-type="COI-statement">
<title>Conflict of interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
<p>The handling editor LS declared a past collaboration with the author AT.</p>
</sec>
<sec id="s9" sec-type="disclaimer">
<title>Publisher&#x2019;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<sec id="s10" sec-type="supplementary-material">
<title>Supplementary material</title>
<p>The Supplementary Material for this article can be found online at: <ext-link ext-link-type="uri" xlink:href="https://www.frontiersin.org/articles/10.3389/fpls.2023.1236147/full#supplementary-material">https://www.frontiersin.org/articles/10.3389/fpls.2023.1236147/full#supplementary-material</ext-link>
</p>
<supplementary-material xlink:href="Table_1.docx" id="ST1" mimetype="application/vnd.openxmlformats-officedocument.wordprocessingml.document"/>
<supplementary-material xlink:href="Table_2.docx" id="ST2" mimetype="application/vnd.openxmlformats-officedocument.wordprocessingml.document"/>
<supplementary-material xlink:href="Table_3.docx" id="ST3" mimetype="application/vnd.openxmlformats-officedocument.wordprocessingml.document"/>
<supplementary-material xlink:href="Image_1.jpeg" id="SF1" mimetype="image/jpeg"/>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Abdelmula</surname> <given-names>A. A.</given-names>
</name>
<name>
<surname>Link</surname> <given-names>W.</given-names>
</name>
<name>
<surname>von Kittlitz</surname> <given-names>E.</given-names>
</name>
<name>
<surname>Stelling</surname> <given-names>D.</given-names>
</name>
</person-group> (<year>1999</year>). <article-title>Heterosis and inheritance of drought tolerance in faba bean, <italic>Vicia faba</italic> L</article-title>. <source>Plant Breed.</source> <volume>118</volume> (<issue>6</issue>), <fpage>485</fpage>&#x2013;<lpage>490</lpage>. doi: <pub-id pub-id-type="doi">10.1046/j.1439-0523.1999.00411.x</pub-id>
</citation>
</ref>
<ref id="B2">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Abid</surname> <given-names>G.</given-names>
</name>
<name>
<surname>Hessini</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Aouida</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Baudoin</surname> <given-names>J. P.</given-names>
</name>
<name>
<surname>Aroua</surname> <given-names>I.</given-names>
</name>
<name>
<surname>Baudoin</surname> <given-names>J.-P.</given-names>
</name>
<etal/>
</person-group>. (<year>2017</year>). <article-title>Agro-physiological and biochemical responses of faba bean (<italic>Vicia faba</italic> L. var &#x2018;minor&#x2019;) genotypes to water deficit stress</article-title>. <source>Biotechnol. Agron. Soc Environ.</source> <volume>21</volume> (<issue>2</issue>), <fpage>146</fpage>&#x2013;<lpage>159</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.25518/1780-4507.13579</pub-id>
</citation>
</ref>
<ref id="B3">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Abou-Zaitoun</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Chandresakar</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Assili</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Shtaya</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Jamous</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Mallah</surname> <given-names>O.</given-names>
</name>
<etal/>
</person-group>. (<year>2018</year>). <article-title>Unlocking the genetic diversity within a middle-east panel of Durum wheat landraces for adaptation to semi-arid climate</article-title>. <source>Agronomy</source> <volume>8</volume>, <elocation-id>233</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.3390/agronomy8100233</pub-id>
</citation>
</ref>
<ref id="B4">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Afzal</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Alghamdi</surname> <given-names>S. S.</given-names>
</name>
<name>
<surname>Migadi</surname> <given-names>H. H.</given-names>
</name>
<name>
<surname>El-Harty</surname> <given-names>E.</given-names>
</name>
<name>
<surname>Al-Faifi</surname> <given-names>S. A.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Agronomical and physiological responses of faba bean genotypes to salt stress</article-title>. <source>Agriculture</source> <volume>12</volume>, <elocation-id>235</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.3390/agriculture12020235</pub-id>
</citation>
</ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Alderfasi</surname> <given-names>A. A.</given-names>
</name>
<name>
<surname>Alghamdi</surname> <given-names>S. S.</given-names>
</name>
</person-group> (<year>2010</year>). <article-title>Integrated water supply with nutrient requirements on growth, photosynthesis productivity, chemical status and seed yield of faba bean</article-title>. <source>Am.-Eurasian J. Agron.</source> <volume>3</volume>, <fpage>8</fpage>&#x2013;<lpage>17</lpage>.</citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Alghamdi</surname> <given-names>S. S.</given-names>
</name>
<name>
<surname>Al-Schameri</surname> <given-names>A. M.</given-names>
</name>
<name>
<surname>Migdadi</surname> <given-names>H. M.</given-names>
</name>
<name>
<surname>Ammar</surname> <given-names>M. H.</given-names>
</name>
<name>
<surname>El-Harty</surname> <given-names>E. H.</given-names>
</name>
<name>
<surname>Khan</surname> <given-names>M. A.</given-names>
</name>
<etal/>
</person-group>. (<year>2015</year>). <article-title>Physiological and molecular characterization of faba bean (<italic>Vicia faba</italic> L.) genotypes for adaptation to drought stress</article-title>. <source>J. Agron. Crop Sci.</source> <volume>201</volume>, <fpage>401</fpage>&#x2013;<lpage>409</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1111/jac.12110</pub-id>
</citation>
</ref>
<ref id="B7">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ali</surname> <given-names>M. B. M.</given-names>
</name>
<name>
<surname>Welna</surname> <given-names>G. C.</given-names>
</name>
<name>
<surname>Sallam</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Martsch</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Balko</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Gebser</surname> <given-names>B.</given-names>
</name>
<etal/>
</person-group>. (<year>2016</year>). <article-title>Association analyses to genetically improve drought and freezing tolerance of faba bean (<italic>Vicia faba</italic> L.)</article-title>. <source>Crop Sci.</source> <volume>56</volume>, <fpage>1036</fpage>&#x2013;<lpage>1048</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.2135/cropsci2015.08.0503</pub-id>
</citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Amede</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Schubert</surname> <given-names>S.</given-names>
</name>
</person-group> (<year>2003</year>). <article-title>Mechanisms of drought resistance in grain legumes. I. Osmotic adjustment</article-title>. <source>Ethiop J. Sci.</source> <volume>26</volume>, <fpage>37</fpage>&#x2013;<lpage>46</lpage>.</citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Amede</surname> <given-names>T.</given-names>
</name>
<name>
<surname>von Kittlitz</surname> <given-names>E.</given-names>
</name>
<name>
<surname>Schubert</surname> <given-names>S.</given-names>
</name>
</person-group> (<year>2001</year>). <article-title>Differential drought responses of faba bean (<italic>Vicia faba</italic> L.) inbred lines</article-title>. <source>J. Agron. Crop Sci.</source> <volume>183</volume>, <fpage>35</fpage>&#x2013;<lpage>45</lpage>. doi: <pub-id pub-id-type="doi">10.1046/j.1439-037x.1999.00310.x</pub-id>
</citation>
</ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ammar</surname> <given-names>M. H.</given-names>
</name>
<name>
<surname>Anwar</surname> <given-names>F.</given-names>
</name>
<name>
<surname>El-Harty</surname> <given-names>E. H.</given-names>
</name>
<name>
<surname>Migdadi</surname> <given-names>H. M.</given-names>
</name>
<name>
<surname>Abdel-Khalik</surname> <given-names>S. M.</given-names>
</name>
<name>
<surname>Al-Faifi</surname> <given-names>S. A.</given-names>
</name>
<etal/>
</person-group>. (<year>2015</year>). <article-title>Physiological and Yield Responses of Faba bean (<italic>Vicia faba</italic> L.) to Drought Stress in Managed and Open Field Environments</article-title>. <source>J. Agron. Crop Sci.</source> <volume>201</volume> (<issue>4</issue>), <fpage>280</fpage>&#x2013;<lpage>287</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1111/jac.12112</pub-id>
</citation>
</ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Annicchiarico</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Iannucci</surname> <given-names>A.</given-names>
</name>
</person-group> (<year>2008</year>). <article-title>Adaptation strategy, germplasm type and adaptive traits for field pea improvement in Italy based on variety responses across climatically contrasting environments</article-title>. <source>Field Crops Res.</source> <volume>108</volume> (<issue>2</issue>), <fpage>133</fpage>&#x2013;<lpage>142</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.fcr.2008.04.004</pub-id>
</citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Arif</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Parveen</surname> <given-names>N.</given-names>
</name>
<name>
<surname>Waheed</surname> <given-names>M. Q.</given-names>
</name>
<name>
<surname>Atif</surname> <given-names>R. M.</given-names>
</name>
<name>
<surname>Waqar</surname> <given-names>I.</given-names>
</name>
<name>
<surname>Shah</surname> <given-names>T. M.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>A comparative study for assessing the drought-tolerance of chickpea under varying natural growth environments</article-title>. <source>Front. Plant Sci.</source> <volume>11</volume>. doi:&#xa0;<pub-id pub-id-type="doi">10.3389/fpls.2020.607869</pub-id>
</citation>
</ref>
<ref id="B13">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Aslam</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Hussain</surname> <given-names>S. B.</given-names>
</name>
<name>
<surname>Baber</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Shaheen</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Aslam</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Waheed</surname> <given-names>R.</given-names>
</name>
<etal/>
</person-group>. (<year>2023</year>). <article-title>Estimation of drought tolerance indices in upland cotton under water deficit conditions</article-title>. <source>Agronomy</source> <volume>13</volume>, <elocation-id>984</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.3390/agronomy13040984</pub-id>
</citation>
</ref>
<ref id="B14">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ayaz</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Mckenzie</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Hill</surname> <given-names>G.</given-names>
</name>
<name>
<surname>Mcneil</surname> <given-names>D.</given-names>
</name>
</person-group> (<year>2004</year>). <article-title>Variability in yield of four grain legume species in a subhumid temperate environment. II. Yield components</article-title>. <source>J. Agric. Sci.</source> <volume>142</volume>, <fpage>21</fpage>&#x2013;<lpage>28</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1017/S0021859604004113</pub-id>
</citation>
</ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Balko</surname> <given-names>C.</given-names>
</name>
</person-group> (<year>2005</year>). <article-title>Drought tolerance in potatoes and faba beans. Through breeding research to suitable crop plants.Forschungs-Report</article-title>. <source>Ernahrung Landwirtschaft Forsten</source> <volume>1)</volume>, <fpage>10</fpage>&#x2013;<lpage>13</lpage>.</citation>
</ref>
<ref id="B16">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bates</surname> <given-names>C. J.</given-names>
</name>
<name>
<surname>Waldren</surname> <given-names>R. P.</given-names>
</name>
<name>
<surname>Teare</surname> <given-names>I. D.</given-names>
</name>
</person-group> (<year>1973</year>). <article-title>Rapid determination of free proline for water - stress studies</article-title>. <source>Plant Soil</source> <volume>39</volume>, <fpage>205</fpage>&#x2013;<lpage>207</lpage>. doi: <pub-id pub-id-type="doi">10.1007/BF00018060</pub-id>
</citation>
</ref>
<ref id="B17">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Becker</surname> <given-names>H.</given-names>
</name>
</person-group> (<year>2011</year>). <source>Pflanzenz&#xfc;chtung</source> (<publisher-loc>Stuttgart</publisher-loc>: <publisher-name>UTB GmbH</publisher-name>).</citation>
</ref>
<ref id="B18">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Benkeblia</surname> <given-names>N.</given-names>
</name>
<name>
<surname>Shinano</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Osaki</surname> <given-names>M.</given-names>
</name>
</person-group> (<year>2007</year>). <article-title>Metabolite profiling and assessment of metabolome compartmentation of soybean leaves using non-aqueous fractionation and GC-MS analysis</article-title>. <source>Metabolomics</source> <volume>3</volume>, <fpage>297</fpage>&#x2013;<lpage>305</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s11306-007-0078-y</pub-id>
</citation>
</ref>
<ref id="B19">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Bouslama</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Schapaugh</surname> <given-names>W. T.</given-names>
</name>
</person-group> (<year>1984</year>). <article-title>Stress tolerance in soybeans. I. Evaluation of three screening techniques for heat and drought tolerance</article-title>. <source>Crop Sci.</source> <volume>24</volume>, <fpage>933</fpage>&#x2013;<lpage>937</lpage>. doi: <pub-id pub-id-type="doi">10.2135/cropsci1984.0011183X002400050026x</pub-id>
</citation>
</ref>
<ref id="B20">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Daryanto</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Jacinthe</surname> <given-names>P. A.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Global synthesis of drought effects on food legume production</article-title>. <source>PloS One</source> <volume>10</volume>, <elocation-id>e0127401</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1371/journal.pone.0127401</pub-id>
</citation>
</ref>
<ref id="B21">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>El-Hendawy</surname> <given-names>S. E.</given-names>
</name>
<name>
<surname>Hassan</surname> <given-names>W. M.</given-names>
</name>
<name>
<surname>Al-Suhaibani</surname> <given-names>N. A.</given-names>
</name>
<name>
<surname>Schmidhalter</surname> <given-names>U.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Spectral assessment of drought tolerance indices and grain yield in advanced spring wheat lines grown under full and limited water irrigation</article-title>. <source>Agric. Water Manage.</source> <volume>182</volume> (<issue>7</issue>), <fpage>1</fpage>&#x2013;<lpage>12</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.agwat.2016.12.003</pub-id>
</citation>
</ref>
<ref id="B22">
<citation citation-type="other">
<collab>FAOSTAT</collab>
(<year>2022</year>). <article-title>FAOSTAT (FAOSTAT)</article-title>. Available at: <uri xlink:href="http://www.fao.org/faostat/en/#data">http://www.fao.org/faostat/en/#data</uri> (Accessed October, 2021).</citation>
</ref>
<ref id="B23">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Farooq</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Hussain</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Siddique</surname> <given-names>K. H. M.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>Drought stress in wheat during flowering and grain-filling periods</article-title>. <source>Crit. Rev. Plant Sci.</source> <volume>33</volume>, <fpage>331</fpage>&#x2013;<lpage>349</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1080/07352689.2014.875291</pub-id>
</citation>
</ref>
<ref id="B24">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Farooq</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Hussain</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Wahid</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Siddique</surname> <given-names>K. H. M.</given-names>
</name>
</person-group> (<year>2012</year>). &#x201c;<article-title>Drought stress in plants: An overview</article-title>,&#x201d; in <source>Plant responses to drought stress: From morphological to molecular features</source>, <edition>1st edn</edition>. Ed. <person-group person-group-type="editor">
<name>
<surname>Aroca</surname> <given-names>R.</given-names>
</name>
</person-group> (<publisher-loc>Berlin Heidelberg</publisher-loc>: <publisher-name>Springer</publisher-name>), <fpage>1</fpage>&#x2013;<lpage>33</lpage>.</citation>
</ref>
<ref id="B25">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Farshadfar</surname> <given-names>E.</given-names>
</name>
<name>
<surname>Mohammadi</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Farshadfar</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Dabiri</surname> <given-names>S.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>Relationships and repeatability of drought tolerance indices in wheat-rye disomic addition lines</article-title>. <source>Aust. J. Crop Sci.</source> <volume>7</volume> (<issue>1</issue>), <fpage>130</fpage>&#x2013;<lpage>138</lpage>.</citation>
</ref>
<ref id="B26">
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Fernandez</surname> <given-names>G. C.</given-names>
</name>
</person-group> (<year>1992</year>). &#x201c;<article-title>Effective selection criteria for assessing plant stress tolerance</article-title>,&#x201d; in <conf-name>Proceedings of the International Symposium on Adaptation of Vegetables and Other Food Crops in Temperature and Water Stress</conf-name>, Adaptation of food crops to temperature and water stress. (<publisher-loc>Shanhua, Taiwan</publisher-loc>: <publisher-name>AVRDC</publisher-name>) Kuo, C. George Book editor. Vol. <volume>18</volume>. <fpage>257</fpage>&#x2013;<lpage>270</lpage>.</citation>
</ref>
<ref id="B27">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gnanasambandam</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Paull</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Torres</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Kaur</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Leonforte</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Li</surname> <given-names>H.</given-names>
</name>
<etal/>
</person-group>. (<year>2012</year>). <article-title>Impact of molecular technologies on faba bean (<italic>Vicia faba</italic> L.) breeding strategies</article-title>. <source>Agronomy</source> <volume>2</volume>, <fpage>132</fpage>&#x2013;<lpage>166</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.3390/agronomy2030132</pub-id>
</citation>
</ref>
<ref id="B28">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Grashoff</surname> <given-names>C.</given-names>
</name>
</person-group> (<year>1990</year>). <article-title>Effect of pattern of water supply on <italic>Vicia faba</italic> L. 1. Dry matter partitioning and yield variability</article-title>. <source>Netherlands J. Agric. Sci.</source> <volume>38</volume> (<issue>1</issue>), <fpage>21</fpage>&#x2013;<lpage>44</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.18174/njas.v38i1.16608</pub-id>
</citation>
</ref>
<ref id="B29">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gutierrez</surname> <given-names>N.</given-names>
</name>
<name>
<surname>P&#xe9;gard</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Balko</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Torres</surname> <given-names>A. M.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Genome-wide association analysis for drought tolerance and associated traits in faba bean (<italic>Vicia faba</italic> L.)</article-title>. <source>Front. Plant Sci.</source> <volume>14</volume>. doi:&#xa0;<pub-id pub-id-type="doi">10.3389/fpls.2023.1091875</pub-id>
</citation>
</ref>
<ref id="B30">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hack</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Bleiholder</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Buhr</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Meier</surname> <given-names>U.</given-names>
</name>
<name>
<surname>Schnock-Fricke</surname> <given-names>U.</given-names>
</name>
<name>
<surname>Weber</surname> <given-names>E.</given-names>
</name>
<etal/>
</person-group>. (<year>1992</year>). <article-title>Einheitliche Codierung der ph&#xe4;nologischen Entwicklungsstadien mono- und dikotyler Pflanzen &#x2013; Erweiterte BBCH-Skala, Allgemein</article-title>. <source>Nachr.bl. Dtsch. Pflanzenschutzd</source> <volume>44</volume> (<issue>12</issue>), <fpage>265</fpage>&#x2013;<lpage>270</lpage>.</citation>
</ref>
<ref id="B31">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hannok</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Kaewunta</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Khunyota</surname> <given-names>W.</given-names>
</name>
<name>
<surname>Kaewnut</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Pachimkanthong</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Tangthavonkarn</surname> <given-names>C.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Exploration of phenotypic dissimilarity for drought tolerance in maize inbred line collection</article-title>. <source>Curr. Appl. Sci. Technol.</source> <volume>23</volume>, <fpage>(1)</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.55003/cast.2022.01.23.001</pub-id>
</citation>
</ref>
<ref id="B32">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hartigan</surname> <given-names>J. A.</given-names>
</name>
<name>
<surname>Wong</surname> <given-names>M. A.</given-names>
</name>
</person-group> (<year>1979</year>). <article-title>Algorithm AS 136: A K-means clustering algorithm</article-title>. <source>Appl. Stat.</source> <volume>28</volume>, <fpage>100</fpage>&#x2013;<lpage>108</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.2307/2346830</pub-id>
</citation>
</ref>
<ref id="B33">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hekneby</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Antolin</surname> <given-names>M. C.</given-names>
</name>
<name>
<surname>Sanchez-Diaz</surname> <given-names>M.</given-names>
</name>
</person-group> (<year>2006</year>). <article-title>Frost resistance and biochemical changes during cold acclimation in different annual legumes</article-title>. <source>Environ. Exp. Bot.</source> <volume>55</volume> (<issue>3</issue>), <fpage>305</fpage>&#x2013;<lpage>314</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.envexpbot.2004.11.010</pub-id>
</citation>
</ref>
<ref id="B34">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jensen</surname> <given-names>E. S.</given-names>
</name>
<name>
<surname>Peoples</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Hauggaard-Nielsen</surname> <given-names>H.</given-names>
</name>
</person-group> (<year>2010</year>). <article-title>Review: Faba beans in cropping systems</article-title>. <source>Field Crop Res.</source> <volume>115</volume>, <fpage>203</fpage>&#x2013;<lpage>216</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.fcr.2009.10.008</pub-id>
</citation>
</ref>
<ref id="B35">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kabbadj</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Makoudi</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Mouradi</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Pauly</surname> <given-names>N.</given-names>
</name>
<name>
<surname>Frendo</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Ghoulam</surname> <given-names>C.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Physiological and biochemical responses involved in water deficit tolerance of nitrogen-fixing <italic>Vicia faba</italic>
</article-title>. <source>PloS One</source> <volume>12</volume> (<issue>12</issue>), <elocation-id>e0190284</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1371/journal</pub-id>
</citation>
</ref>
<ref id="B36">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Katsoulieri</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Papastylianou</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Travlos</surname> <given-names>I.</given-names>
</name>
<name>
<surname>Vlachostergios</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Tigka</surname> <given-names>E.</given-names>
</name>
<name>
<surname>Kargiotidou</surname> <given-names>A.</given-names>
</name>
<etal/>
</person-group>. (<year>2020</year>). <article-title>Yield performance of faba bean cultivars under different environmental conditions in Greece</article-title>. <source>Bull. UASVM Horticult.</source> <volume>77</volume> (<issue>2</issue>), <fpage>30</fpage>&#x2013;<lpage>36</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.15835/buasvmcn-hort:2020.0016</pub-id>
</citation>
</ref>
<ref id="B37">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kaur</surname> <given-names>G.</given-names>
</name>
<name>
<surname>Asthir</surname> <given-names>B.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Proline: a key player in plant abiotic stress tolerance</article-title>. <source>Biol. Plant</source> <volume>59</volume>, <fpage>609</fpage>&#x2013;<lpage>619</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s10535-015-0549-3</pub-id>
</citation>
</ref>
<ref id="B38">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Khan</surname> <given-names>H. R.</given-names>
</name>
<name>
<surname>Link</surname> <given-names>W.</given-names>
</name>
<name>
<surname>Hocking</surname> <given-names>T. J.</given-names>
</name>
<name>
<surname>Stoddard</surname> <given-names>F. L.</given-names>
</name>
</person-group> (<year>2007</year>). <article-title>Evaluation of physiological traits for improving drought tolerance in faba bean (<italic>Vicia faba</italic> L.)</article-title>. <source>Plant Soil</source> <volume>292</volume>, <fpage>205</fpage>&#x2013;<lpage>217</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s11104-007-9217-5</pub-id>
</citation>
</ref>
<ref id="B39">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Khan</surname> <given-names>H. R.</given-names>
</name>
<name>
<surname>Paull</surname> <given-names>J. G.</given-names>
</name>
<name>
<surname>Siddique</surname> <given-names>K. H. M.</given-names>
</name>
<name>
<surname>Stoddard</surname> <given-names>F. L.</given-names>
</name>
</person-group> (<year>2010</year>). <article-title>Faba bean breeding for drought-affected environments: A physiological and agronomic perspective</article-title>. <source>Field Crops Res.</source> <volume>115</volume>, <fpage>279</fpage>&#x2013;<lpage>286</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.fcr.2009.09.003</pub-id>
</citation>
</ref>
<ref id="B40">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kassambara</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Mundt</surname> <given-names>F.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Extract and visualize the results of multivariate data analyses [R package factoextra] (CRAN)</article-title>.</citation>
</ref>
<ref id="B41">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Lafitte</surname> <given-names>H. R.</given-names>
</name>
<name>
<surname>Blum</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Atlin</surname> <given-names>G.</given-names>
</name>
</person-group> (<year>2003</year>). &#x201c;<article-title>Using Secondary Traits to Help Identify Drought-Tolerant Genotypes</article-title>,&#x201d; in <source>Breeding Rice for Drought-Prone Environments</source>. Eds. <person-group person-group-type="editor">
<name>
<surname>Fischer</surname> <given-names>K. S.</given-names>
</name>
<name>
<surname>Lafitte</surname> <given-names>R. H.</given-names>
</name>
<name>
<surname>Fukai</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Atlin</surname> <given-names>G.</given-names>
</name>
<name>
<surname>Hardy</surname> <given-names>B.</given-names>
</name>
</person-group> (<publisher-loc>Manila</publisher-loc>: <publisher-name>International Rice Research Institute</publisher-name>), <fpage>37</fpage>&#x2013;<lpage>48</lpage>.</citation>
</ref>
<ref id="B42">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Link</surname> <given-names>W.</given-names>
</name>
<name>
<surname>Abdelmula</surname> <given-names>A. A.</given-names>
</name>
<name>
<surname>Kittlitz</surname> <given-names>E. V.</given-names>
</name>
<name>
<surname>Bruns</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Riemer</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Stelling</surname> <given-names>D.</given-names>
</name>
</person-group> (<year>1999</year>). <article-title>Genotypic variation for drought tolerance in <italic>Vicia faba</italic>
</article-title>. <source>Plant Breed</source> <volume>118</volume>, <fpage>477</fpage>&#x2013;<lpage>483</lpage>. doi: <pub-id pub-id-type="doi">10.1046/j.1439-0523.1999.00412.x</pub-id>
</citation>
</ref>
<ref id="B43">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lopez</surname> <given-names>F. B.</given-names>
</name>
<name>
<surname>Johansen</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Chauhan</surname> <given-names>Y. S.</given-names>
</name>
</person-group> (<year>1996</year>). <article-title>Effects of timing of drought stress on phenology, yield and yield components of short-duration pigeonpea</article-title>. <source>J. Agron. Crop Sci.</source> <volume>177</volume> (<issue>5</issue>), <fpage>311</fpage>&#x2013;<lpage>320</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1111/j.1439-037X.1996.tb00251.x</pub-id>
</citation>
</ref>
<ref id="B44">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Manning</surname> <given-names>B. K.</given-names>
</name>
<name>
<surname>Adhikari</surname> <given-names>K. N.</given-names>
</name>
<name>
<surname>Trethowan</surname> <given-names>R.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Impact of sowing time, genotype, environment and maturity on biomass and yield components in faba bean (<italic>Vicia faba</italic>)</article-title>. <source>Crop Pasture Sci.</source> <volume>71</volume>, <fpage>147</fpage>&#x2013;<lpage>154</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1071/CP19214</pub-id>
</citation>
</ref>
<ref id="B45">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Mansour</surname> <given-names>E.</given-names>
</name>
<name>
<surname>Mahgoub</surname> <given-names>H. A. M.</given-names>
</name>
<name>
<surname>Mahgoub</surname> <given-names>S. A.</given-names>
</name>
<name>
<surname>El-Sobky</surname> <given-names>E.-S.E.A.</given-names>
</name>
<name>
<surname>Abdul-Hamid</surname> <given-names>M. I.</given-names>
</name>
<name>
<surname>Kamara</surname> <given-names>M. M.</given-names>
</name>
<etal/>
</person-group>. (<year>2021</year>). <article-title>Enhancement of drought tolerance in diverse Vicia faba cultivars by inoculation with plant growth-promoting rhizobacteria under newly reclaimed soil conditions</article-title>. <source>Sci. Rep.</source> <volume>11</volume> (<issue>1</issue>), <fpage>24142</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/s41598-021-02847-2</pub-id>
</citation>
</ref>
<ref id="B46">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>McDonald</surname> <given-names>G. K.</given-names>
</name>
<name>
<surname>Paulsen</surname> <given-names>G. M.</given-names>
</name>
</person-group> (<year>1997</year>). <article-title>High temperature effects on photosynthesis and water relations of grain legumes</article-title>. <source>Plant Soil</source> <volume>196</volume>, <fpage>47</fpage>&#x2013;<lpage>58</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1023/A:1004249200050</pub-id>
</citation>
</ref>
<ref id="B47">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Memari</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Dadashi</surname> <given-names>M. R.</given-names>
</name>
<name>
<surname>Sheikh</surname> <given-names>F.</given-names>
</name>
<name>
<surname>Ajamnoruzi</surname> <given-names>H.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Evaluation of drought stress tolerance in some faba bean genotypes using drought tolerance indices</article-title>. <source>Romanian Agric. Res</source>. n.39. <volume>11</volume>, <fpage>139</fpage>&#x2013;<lpage>151</lpage> Print ISSN 1222-4227; online ISSN 2067-5720. doi: <pub-id pub-id-type="doi">10.59665/rar3904</pub-id>
</citation>
</ref>
<ref id="B48">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Migdadi</surname> <given-names>H. M.</given-names>
</name>
<name>
<surname>El-Harty</surname> <given-names>E. H.</given-names>
</name>
<name>
<surname>Salamh</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Khan</surname> <given-names>M. A.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Yield and proline content of faba bean genotypes under water stress treatments</article-title>. <source>J. Anim. Plant Sci.</source> <volume>26</volume> (<issue>6</issue>), <fpage>1772</fpage>&#x2013;<lpage>1779</lpage>.</citation>
</ref>
<ref id="B49">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Mohamamdi</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Etminan</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Shooshtari</surname> <given-names>L.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Statistical methods for identifying traits associated with high yield potential in durum wheat under drought conditions</article-title>. <source>Indian J. Genet. Plant Breed.</source> <volume>77</volume> (<issue>4</issue>), <fpage>513</fpage>&#x2013;<lpage>523</lpage>. doi: <pub-id pub-id-type="doi">10.5958/0975-6906.2017.00068.2</pub-id>
</citation>
</ref>
<ref id="B50">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Morgan</surname> <given-names>J. M.</given-names>
</name>
</person-group> (<year>1984</year>). <article-title>Osmoregulation and water stress in higher plants</article-title>. <source>Ann. Rev. Plant Physiol.</source> <volume>35</volume>, <fpage>299</fpage>&#x2013;<lpage>319</lpage>. doi: <pub-id pub-id-type="doi">10.1146/annurev.pp.35.060184.001503</pub-id>
</citation>
</ref>
<ref id="B51">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Muktadir</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Adhikari</surname> <given-names>K. M.</given-names>
</name>
<name>
<surname>Merchant</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Belachew</surname> <given-names>K. Y.</given-names>
</name>
<name>
<surname>Vandenberg</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Stoddard</surname> <given-names>F. L.</given-names>
</name>
<etal/>
</person-group>. (<year>2020</year>). <article-title>Physiological and biochemical basis of faba bean breeding for drought adaptation &#x2013; a review</article-title>. <source>Agronomy</source> <volume>10</volume>, <elocation-id>1345</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.3390/agronomy10091345</pub-id>
</citation>
</ref>
<ref id="B52">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>M&#xfc;ller</surname> <given-names>U.</given-names>
</name>
<name>
<surname>Meyer</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Ehlers</surname> <given-names>W.</given-names>
</name>
<name>
<surname>Bohm</surname> <given-names>W.</given-names>
</name>
</person-group> (<year>1985</year>). <article-title>Water-uptake and water-use of field beans and oats grown on a loess-derived grey-brown podzolic soil (EUTROBORALF)</article-title>. <source>Z. Pflanzenern&#xe4;hr. Bodenkd.</source> <volume>4</volume>, <fpage>389</fpage>&#x2013;<lpage>400</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1002/jpln.19851480404</pub-id>
</citation>
</ref>
<ref id="B53">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Mutlu</surname> <given-names>S. S.</given-names>
</name>
<name>
<surname>Mutlu</surname> <given-names>N.</given-names>
</name>
<name>
<surname>Karaguzel</surname> <given-names>O.</given-names>
</name>
<name>
<surname>Sari</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Shearman</surname> <given-names>R. C.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Drought resistance of Bermudagrass accessions collected from Eastern Mediterranean</article-title>. <source>Eur.J.Hortic.Sci</source> <volume>87</volume> (<issue>3</issue>), <fpage>1</fpage>&#x2013;<lpage>12</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.17660/eJHS.2022/034</pub-id>
</citation>
</ref>
<ref id="B54">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Mwanamwenge</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Loss</surname> <given-names>S. P.</given-names>
</name>
<name>
<surname>Siddique</surname> <given-names>K. H. M.</given-names>
</name>
<name>
<surname>Cocks</surname> <given-names>P. S.</given-names>
</name>
</person-group> (<year>1999</year>). <article-title>Effect of water stress during floral initiation, flowering and podding on the growth and yield of faba bean (Vicia faba L.)</article-title>. <source>Eur. J. Agron.</source> <volume>11</volume> (<issue>1</issue>), <fpage>1</fpage>&#x2013;<lpage>11</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/S1161-0301(99)00003-9</pub-id>
</citation>
</ref>
<ref id="B55">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Nadeem</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Li</surname> <given-names>J. J.</given-names>
</name>
<name>
<surname>Yahya</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Sher</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Ma</surname> <given-names>C. X.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>X. B.</given-names>
</name>
<etal/>
</person-group>. (<year>2019</year>). <article-title>Research progress and perspective on drought stress in legumes: A review</article-title>. <source>Int. J. Mol. Sci.</source> <volume>20</volume> (<issue>10</issue>), <elocation-id>2541</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.3390/ijms20102541</pub-id>
</citation>
</ref>
<ref id="B56">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Or-Rashid</surname> <given-names>M. H.</given-names>
</name>
<name>
<surname>Miah</surname> <given-names>M. A. B.</given-names>
</name>
<name>
<surname>Islam</surname> <given-names>S. M. S.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>
<italic>In vitro</italic> morphophysiological screening of drought tolerant potato genotypes</article-title>. <source>In Vitro Cell. Dev. Biol.</source> <volume>57</volume>, <fpage>519</fpage>&#x2013;<lpage>528</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s11627-020-10155-z</pub-id>
</citation>
</ref>
<ref id="B57">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ouji</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Naouari</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Mouelhi</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Ben Younes</surname> <given-names>M.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Yield and yield components of faba bean (Vicia faba L.) as influenced by supplemental irrigation under semi-arid region of Tunisia</article-title>. <source>World J. Agric. Res. World J. Agric. Res.</source> <volume>5</volume> (<issue>1</issue>), <fpage>52</fpage>&#x2013;<lpage>57</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.12691/wjar-5-1-7</pub-id>
</citation>
</ref>
<ref id="B58">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Pilbeam</surname> <given-names>C. J.</given-names>
</name>
<name>
<surname>Akatse</surname> <given-names>J. K.</given-names>
</name>
<name>
<surname>Hebblethwaite</surname> <given-names>P. D.</given-names>
</name>
<name>
<surname>Wright</surname> <given-names>C. D.</given-names>
</name>
</person-group> (<year>1992</year>). <article-title>Yield production in two contrasting forms of spring-sown faba beans in relation to water supply</article-title>. <source>Field Crops Res.</source> <volume>29</volume>, <fpage>273</fpage>&#x2013;<lpage>287</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/0378-4290(92)90030-D</pub-id>
</citation>
</ref>
<ref id="B59">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Plies-Balzer</surname> <given-names>E.</given-names>
</name>
<name>
<surname>Kong</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Schubert</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Mengel</surname> <given-names>K.</given-names>
</name>
</person-group> (<year>1995</year>). <article-title>Effect of water stress on plant growth, nitrogenase activity and nitrogen economy of 4 different cultivars of <italic>Vicia faba L</italic>
</article-title>. <source>Eur. J. Agron.</source> <volume>4</volume> (<issue>2</issue>), <fpage>167</fpage>&#x2013;<lpage>173</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/S1161-0301(14)80043-9</pub-id>
</citation>
</ref>
<ref id="B60">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ramirez-Vallejo</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Kelly</surname> <given-names>J. D.</given-names>
</name>
</person-group> (<year>1998</year>). <article-title>Traits related to drought resistance in common bean</article-title>. <source>Euphytica</source> <volume>99</volume> (<issue>2</issue>), <fpage>127</fpage>&#x2013;<lpage>136</lpage>. doi: <pub-id pub-id-type="doi">10.1023/A:1018353200015</pub-id>
</citation>
</ref>
<ref id="B61">
<citation citation-type="book">
<person-group person-group-type="author">
<collab>R Core Team</collab>
</person-group> (<year>2022</year>). <source>A Language and Environment for Statistical Computing</source> (<publisher-loc>Vienna, Austria</publisher-loc>: <publisher-name>R Foundation for Statistical Computing</publisher-name>). URL <uri xlink:href="https://www.R-project.org/">https://www.R-project.org/</uri>.</citation>
</ref>
<ref id="B62">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Redden</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Paull</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Zong</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Sass</surname> <given-names>O.</given-names>
</name>
<name>
<surname>Yang</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Lintg</surname> <given-names>L.</given-names>
</name>
</person-group> (<year>2014</year>). &#x201c;<article-title>Faba bean</article-title>,&#x201d; in <source>Broadening the genetic base of grain legumes</source>. Eds. <person-group person-group-type="editor">
<name>
<surname>Singh</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Bisht</surname> <given-names>I. S.</given-names>
</name>
<name>
<surname>Dutta</surname> <given-names>M.</given-names>
</name>
</person-group> (<publisher-loc>New Delhi</publisher-loc>: <publisher-name>Springer</publisher-name>), <fpage>75</fpage>&#x2013;<lpage>93</lpage>.</citation>
</ref>
<ref id="B63">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Reyes</surname> <given-names>J. A. O.</given-names>
</name>
<name>
<surname>Delfin</surname> <given-names>E. F.</given-names>
</name>
<name>
<surname>Reyes</surname> <given-names>M. E. C.</given-names>
</name>
<name>
<surname>Ocampo</surname> <given-names>E. T. M.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Physiological and yield responses of selected mungbean genotypes to terminal drought</article-title>. <source>Philipp J. Sci.</source> <volume>151</volume> (<issue>2</issue>), <fpage>739</fpage>&#x2013;<lpage>750</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.56899/151.02.16</pub-id>
</citation>
</ref>
<ref id="B64">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ricciardi</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Polignano</surname> <given-names>G. B.</given-names>
</name>
<name>
<surname>De Giovanni</surname> <given-names>C.</given-names>
</name>
</person-group> (<year>2001</year>). <article-title>Genotypic response of faba bean to water stress</article-title>. <source>Euphytica</source> <volume>118</volume>, <fpage>39</fpage>&#x2013;<lpage>46</lpage>. doi: <pub-id pub-id-type="doi">10.1023/A:1004078017159</pub-id>
</citation>
</ref>
<ref id="B65">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Richards</surname> <given-names>R. A.</given-names>
</name>
</person-group> (<year>2006</year>). <article-title>Physiological traits used in the breeding of new cultivars for water-scarce environments</article-title>. <source>Agric. Water Manag</source> <volume>80</volume> (<issue>1-3</issue>), <fpage>197</fpage>&#x2013;<lpage>211</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.agwat.2005.07.013</pub-id>
</citation>
</ref>
<ref id="B66">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Romagosa</surname> <given-names>I.</given-names>
</name>
<name>
<surname>van Eeuwijk</surname> <given-names>F. A.</given-names>
</name>
<name>
<surname>Thomas</surname> <given-names>W. T.</given-names>
</name>
</person-group> (<year>2009</year>). &#x201c;<article-title>Statistical analyses of genotype by environment data</article-title>,&#x201d; in <source>Cereals</source> (<publisher-loc>Dordrecht, The Netherlands</publisher-loc>: <publisher-name>Springer</publisher-name>), <fpage>291</fpage>&#x2013;<lpage>331</lpage>.</citation>
</ref>
<ref id="B67">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sanchez</surname> <given-names>D. H.</given-names>
</name>
<name>
<surname>Schwabe</surname> <given-names>F.</given-names>
</name>
<name>
<surname>Erban</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Udvardi</surname> <given-names>M. K.</given-names>
</name>
<name>
<surname>Kopka</surname> <given-names>J.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>Comparative metabolomics of drought acclimation in model and forage legumes</article-title>. <source>Plant Cell Environ.</source> <volume>35</volume>, <fpage>136</fpage>&#x2013;<lpage>149</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1111/j.1365-3040.2011.02423.x</pub-id>
</citation>
</ref>
<ref id="B68">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Santos</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Carvalho</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Rosa</surname> <given-names>E.</given-names>
</name>
<name>
<surname>Carnicle</surname> <given-names>V.</given-names>
</name>
<name>
<surname>Castro</surname> <given-names>I.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Root and agro-morphological traits performance in cowpea under drought stress</article-title>. <source>Agronomy</source> <volume>10</volume> (<issue>10</issue>), <elocation-id>1604</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.3390/agronomy10101604</pub-id>
</citation>
</ref>
<ref id="B69">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sekara</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Poniedzialek</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Ciura</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Jedrszczyk</surname> <given-names>E.</given-names>
</name>
</person-group> (<year>2001</year>). <article-title>The effect of meteorological factors upon flowering and pod setting of faba bean (Vicia faba L.) at different sowing times</article-title>. <source>Veg. Crop Res. Bull.</source> <volume>54</volume> (<issue>1</issue>), <fpage>65</fpage>&#x2013;<lpage>68</lpage>.</citation>
</ref>
<ref id="B70">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Shah</surname> <given-names>T. M.</given-names>
</name>
<name>
<surname>Imran</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Atta</surname> <given-names>B. M.</given-names>
</name>
<name>
<surname>Ashraf</surname> <given-names>M. Y.</given-names>
</name>
<name>
<surname>Hameed</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Waqar</surname> <given-names>I.</given-names>
</name>
<etal/>
</person-group>. (<year>2020</year>). <article-title>Selection and screening of drought tolerant high yielding chickpea genotypes based on physio-biochemical indices and multi-environmental yield trials</article-title>. <source>BMC Plant Biol.</source> <volume>20</volume>, <fpage>1</fpage>&#x2013;<lpage>16</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1186/s12870-020-02381-9</pub-id>
</citation>
</ref>
<ref id="B71">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Shakeel</surname> <given-names>A. A.</given-names>
</name>
<name>
<surname>Xiao-yu</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Long-chang</surname> <given-names>W.</given-names>
</name>
<name>
<surname>Muhammad</surname> <given-names>F. S.</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>L.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>Morphological, physiological and biochemical responses of plants to drought stress</article-title>. <source>Afr. J. Agric. Res.</source> <volume>6</volume> (<issue>9</issue>), <fpage>2026</fpage>&#x2013;<lpage>2032</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.5897/AJAR10.027</pub-id>
</citation>
</ref>
<ref id="B72">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sharifi</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Astereki</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Sheikh</surname> <given-names>F.</given-names>
</name>
<name>
<surname>Khorasanizadeh</surname> <given-names>M.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Evaluation of faba bean genotypes in normal and drought stress conditions by tolerance and susceptibility indices</article-title>. <source>CAJPSI</source> <volume>1</volume> (<issue>4</issue>), <fpage>176</fpage>&#x2013;<lpage>191</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.22034/CAJPSI.2021.04.0</pub-id>
</citation>
</ref>
<ref id="B73">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Siahsar</surname> <given-names>B. A.</given-names>
</name>
<name>
<surname>Ganjali</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Allahdoo</surname> <given-names>M.</given-names>
</name>
</person-group> (<year>2010</year>). <article-title>Evaluation of drought tolerance indices and their relationship with grain yield of lentil lines in drought-stressed and irrigated environments</article-title>. <source>Aust. J. basic Appl. Sci.</source> <volume>4</volume> (<issue>9</issue>), <fpage>4336</fpage>&#x2013;<lpage>4346</lpage>.</citation>
</ref>
<ref id="B74">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Siddiqui</surname> <given-names>M. H.</given-names>
</name>
<name>
<surname>Al-Khaishany</surname> <given-names>M. Y.</given-names>
</name>
<name>
<surname>Al-Qutami</surname> <given-names>M. A.</given-names>
</name>
<name>
<surname>Al-Whaibi</surname> <given-names>M. H.</given-names>
</name>
<name>
<surname>Grover</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Ali</surname> <given-names>H. M.</given-names>
</name>
<etal/>
</person-group>. (<year>2015</year>). <article-title>Response of different genotypes of faba bean plant to drought stress</article-title>. <source>Int. J. Mol. Sci.</source> <volume>16</volume>, <fpage>10214</fpage>&#x2013;<lpage>10227</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.3390/ijms160510214</pub-id>
</citation>
</ref>
<ref id="B75">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sofi</surname> <given-names>P. A.</given-names>
</name>
<name>
<surname>Djanaguiraman</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Siddique</surname> <given-names>K. H. M.</given-names>
</name>
<name>
<surname>Prasad</surname> <given-names>P. V. V.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Reproductive fitness in common bean (<italic>Phaseolus vulgaris</italic> L.) under drought stress is associated with root length and volume</article-title>. <source>Ind. J. Plant Physiol.</source> <volume>23</volume>, <fpage>796</fpage>&#x2013;<lpage>809</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s40502-018-0429-x</pub-id>
</citation>
</ref>
<ref id="B76">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Stoddard</surname> <given-names>F. L.</given-names>
</name>
<name>
<surname>Balko</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Erskine</surname> <given-names>W.</given-names>
</name>
<name>
<surname>Khan</surname> <given-names>H. R.</given-names>
</name>
<name>
<surname>Link</surname> <given-names>W.</given-names>
</name>
<name>
<surname>Sarker</surname> <given-names>A.</given-names>
</name>
</person-group> (<year>2006</year>). <article-title>Screening techniques and sources of resistance to abiotic stresses in cool-season food legumes</article-title>. <source>Euphytica</source> <volume>147</volume> (<issue>1-2</issue>), <fpage>167</fpage>&#x2013;<lpage>186</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s10681-006-4723-8</pub-id>
</citation>
</ref>
<ref id="B77">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Swann</surname> <given-names>A. L.</given-names>
</name>
<name>
<surname>Hoffman</surname> <given-names>F. M.</given-names>
</name>
<name>
<surname>Koven</surname> <given-names>C. D.</given-names>
</name>
<name>
<surname>Randerson</surname> <given-names>J. T.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Plant responses to increasing CO2 reduce estimates of climate impacts on drought severity</article-title>. <source>Proc. Natl. Acad. Sci. U S A.</source> <volume>113</volume> (<issue>36</issue>), <fpage>10019</fpage>&#x2013;<lpage>10024</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1073/pnas.1604581113</pub-id>
</citation>
</ref>
<ref id="B78">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Tadesse</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Fikere</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Legesse</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Parven</surname> <given-names>A.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>Correlation and path coefficient analysis of yield and its components in faba beans</article-title>. <source>Int. J. Biodivers. Conserv.</source> <volume>3</volume>, <fpage>376</fpage>&#x2013;<lpage>382</lpage>.</citation>
</ref>
<ref id="B79">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Tofiq</surname> <given-names>S. E.</given-names>
</name>
<name>
<surname>Aziz</surname> <given-names>O. K.</given-names>
</name>
<name>
<surname>Salih</surname> <given-names>S. H.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Correlation and path coefficient analysis of seed yield and yield components in some faba bean genotypes in Sulaimani region</article-title>. <source>ARO-The Sci. J. Koya Univ.</source> <volume>4</volume> (<issue>2</issue>), <fpage>1</fpage>&#x2013;<lpage>6</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.14500/aro.10081</pub-id>
</citation>
</ref>
<ref id="B80">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Torres</surname> <given-names>A. M.</given-names>
</name>
<name>
<surname>Avila</surname> <given-names>C. M.</given-names>
</name>
<name>
<surname>Gutierrez</surname> <given-names>N.</given-names>
</name>
<name>
<surname>Palomino</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Moreno</surname> <given-names>M. T.</given-names>
</name>
<name>
<surname>Cubero</surname> <given-names>J. I.</given-names>
</name>
</person-group> (<year>2010</year>). <article-title>Marker-assisted selection in faba bean (Vicia faba L.)</article-title>. <source>Field Crops Res.</source> <volume>115</volume> (<issue>3</issue>), <fpage>243</fpage>&#x2013;<lpage>252</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.fcr.2008.12.002</pub-id>
</citation>
</ref>
<ref id="B81">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Toum</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Perez-Borroto</surname> <given-names>L. S.</given-names>
</name>
<name>
<surname>Pe&#xf1;a-Malavera</surname> <given-names>A. N.</given-names>
</name>
<name>
<surname>Luque</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Welin</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Berenstein</surname> <given-names>A.</given-names>
</name>
<etal/>
</person-group>. (<year>2022</year>). <article-title>Selecting putative drought-tolerance markers in two contrasting soybeans</article-title>. <source>Sci. Rep.</source> <volume>12</volume>, <fpage>10872</fpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1038/s41598-022-14334-3</pub-id>
</citation>
</ref>
<ref id="B82">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Tuberosa</surname> <given-names>R.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>Phenotyping for drought tolerance of crops in the genomics era</article-title>. <source>Front. Physiol.</source> <volume>3</volume>. doi:&#xa0;<pub-id pub-id-type="doi">10.3389/fphys.2012.00347</pub-id>
</citation>
</ref>
<ref id="B83">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wahab</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Abdi</surname> <given-names>G.</given-names>
</name>
<name>
<surname>Saleem</surname> <given-names>M. H.</given-names>
</name>
<name>
<surname>Ali</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Ullah</surname> <given-names>S.</given-names>
</name>
<name>
<surname>Shah</surname> <given-names>W.</given-names>
</name>
<etal/>
</person-group>. (<year>2022</year>). <article-title>Plants physio-biochemical and phyto-hormonal responses to alleviate the adverse effects of drought stress: A comprehensive review</article-title>. <source>Plants-Basel</source> <volume>11</volume>, <elocation-id>1620</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.3390/plants11131620</pub-id>
</citation>
</ref>
<ref id="B84">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Weng</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Nakashima</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Villanueva-Morales</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Stewart</surname> <given-names>J. R.</given-names>
</name>
<name>
<surname>Sacks</surname> <given-names>E.</given-names>
</name>
<name>
<surname>Yamada</surname> <given-names>T.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Assessment of drought tolerance of Miscanthus genotypes through dry-down treatment and fixed-soil-moisture-content techniques</article-title>. <source>Agriculture</source> <volume>12</volume>, <elocation-id>6</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.3390/agriculture12010006</pub-id>
</citation>
</ref>
<ref id="B85">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yemm</surname> <given-names>E. W.</given-names>
</name>
<name>
<surname>Willis</surname> <given-names>A. J.</given-names>
</name>
</person-group> (<year>1954</year>). <article-title>The estimation of carbohydrates in plant extracts by anthrone</article-title>. <source>Biochem. J.</source> <volume>57</volume>, <fpage>508</fpage>&#x2013;<lpage>514</lpage>. doi: <pub-id pub-id-type="doi">10.1042/bj0570508</pub-id>
</citation>
</ref>
<ref id="B86">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yoshiba</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Kiyosue</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Nakashima</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Yamaguchi-Shinozaki</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Shinozaki</surname> <given-names>K.</given-names>
</name>
</person-group> (<year>1997</year>). <article-title>Regulation of levels of proline as an osmolyte in plants under water stress</article-title>. <source>Plant Cell Physiol.</source> <volume>38</volume> (<issue>10</issue>), <fpage>1095</fpage>&#x2013;<lpage>1102</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1093/oxfordjournals.pcp.a029093</pub-id>
</citation>
</ref>
</ref-list>
</back>
</article>