Next Article in Journal
Genetic Diversity of Stratiotes aloides L. (Hydrocharitaceae) Stands across Europe
Next Article in Special Issue
Evaluation and Selection of Interspecific Lines of Groundnut (Arachis hypogaea L.) for Resistance to Leaf Spot Disease and for Yield Improvement
Previous Article in Journal
The Interplay between Light Quality and Biostimulant Application Affects the Antioxidant Capacity and Photosynthetic Traits of Soybean (Glycine max L. Merrill)
Previous Article in Special Issue
Exploring the Loci Responsible for Awn Development in Rice through Comparative Analysis of All AA Genome Species
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Plants
Volume 10
Issue 5
10.3390/plants10050862
plants-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Phenotypic Diversity and Productivity of Medicago sativa Subspecies from Drought-Prone Environments in Mediterranean Type Climates

by
Luis Inostroza
1,*,
Soledad Espinoza
2,
Viviana Barahona
2,
Macarena Gerding
3,
Alan Humphries
4,
Alejandro del Pozo
5 and
Carlos Ovalle
6,*
1
CRI-Quilamapu, Instituto de Investigaciones Agropecuaria, Chillán 3780000, Chile
2
CRI-Raihuen Instituto de Investigaciones Agropecuaria, Cauquenes 3690000, Chile
3
Facultad de Agronomía, Universidad de Concepción, Chillán 3780000, Chile
4
South Australian Research and Development Institute (SARDI), Adelaide, SA 5000, Australia
5
Centro de Mejoramiento Genético y Fenómica Vegetal, Facultad de Ciencias Agrarias, Universidad de Talca, Talca 3460000, Chile
6
CRI-La Cruz, Instituto de Investigaciones Agropecuaria, La Cruz 228000, Chile
*
Authors to whom correspondence should be addressed.
Plants 2021, 10(5), 862; https://doi.org/10.3390/plants10050862
Submission received: 25 March 2021 / Revised: 15 April 2021 / Accepted: 19 April 2021 / Published: 24 April 2021
(This article belongs to the Special Issue Pre-breeding towards the Effective Utilization of Plant Genetic Resources)
Browse Figures
Versions Notes

Abstract

:
The phenotypic diversity and productivity of a diverse alfalfa (M. sativa subspp.) panel of cultivars, landraces and wild relatives with putative drought tolerance were evaluated in two Mediterranean environments (central Chile and Southern Australia). In Chile, 70 accessions were evaluated in rainfed conditions and in Australia 30 accessions under rainfed and irrigated conditions, during three growing seasons. Large phenotypic variation was observed among and within subspecies for NDVI, stem length, intercepted PAR and forage yield. Principal component analysis indicated that the first two principal components (PC) accounted for 84.2% of total variance; fall dormancy, taxa, and breeding status were closely related to the agronomical performance of alfalfa accessions. Forage yield varied largely among accessions across years and locations. A linear relationship was found between annual forage yield and annual water added to the experiments (R2 = 0.60, p < 0.001). The GxE analysis for forage yield allowed the detection of the highest yielding accessions for each of the two mega-environments identified. The accessions CTA002 and CTA003 showed greater forage yield in both Chile and Australia environments. It is concluded that new breeding lines derived from crosses between cultivated alfalfa (M. sativa subsp. sativa) and wild relatives belonging to the primary (M. sativa subsp. falcata) and tertiary (M. arborea) gene pool, achieve outstanding agronomical performance in drought-prone environments.

1. Introduction

Mediterranean drought-prone environments are one the most threatened by climate change [1,2]. The increase in temperature along with the decline and larger interannual variability of rainfall, [3,4] are affecting dryland farming systems and their profitability, due to lower water availability for forage yield and therefore livestock production [5].
Alfalfa (Medicago sativa) is a perennial legume of high protein content and yield potential in Mediterranean type climates [6,7]. Its deep root system allows extending the growing season into early summer and autumn period and consequently increasing forage yield in rainfed environments [8].
There is evidence that the genetic diversity of Medicago sativa has been reduced during the domestication process; indeed, cultivated alfalfa has lost about 30% of its genetic diversity compared to wild populations [9]. It has been proposed that landraces and wild relatives of alfalfa encompass greater genotypic and phenotypic diversity than commercial cultivars, as well as adaptive traits to survive in extreme drought stress conditions [10,11,12,13].
Alfalfa belongs to the Medicago sativa-falcata complex, which includes subspecies that are partially sympatric and potentially freely hybridizing in nature [11]. The taxa encompassed in the complex contain both diploids and tetraploids, and gene flow can occur within and between ploidy levels. The most up-to-date list of taxa in the M. sativa-falcata complex include: M. sativa subsp. caerulea (2n), M. sativa subsp. falcata (2n and 4n), M. sativa subsp. x hemicycla (2n), M. sativa subsp. glutinosa (2n), M. sativa subsp. sativa (4n), M. sativa subsp. x varia (4n), and M. sativa subsp. glomerata (4n) [14]. These subspecies represent the primary gene pool for alfalfa improvement [10,11,14]. Recently, Humphries et al. [10] pointed out that primary gene pool offers great potential for alfalfa drought tolerance improvement. Indeed, most of the subspecies are originally from arid zones or drought-prone environments of Eurasia. For instance, M. sativa subsp. caerulea germplasm collected from the Caspian shore is reported to be late flowering, high yielding, adapted to arid conditions, and over-grazing. M. sativa subsp. falcata hybridizes readily with wild populations of alfalfa (M. sativa subsp. sativa) at the tetraploid level. These two taxa have hybridized extensively, giving rise to stabilized introgressant populations, which have been assigned to the collective taxon M. sativa subsp. x varia. Many (by no means all) populations of subsp. falcata and varia occur in arid areas, and are potential sources of drought-resistant germplasm [10]. Some wild M. sativa subsp. falcata populations have shown lower stomatal conductance [15] and higher leaf chlorophyll content under severe drought compared to domesticated alfalfa [16], indicating ability to delay leaf senescence and retain photosynthesis.
Medicago arborea L. is a leguminous shrub well-adapted to the Mediterranean area. It has many traits of potential use in alfalfa breeding; can reach up to 4 m tall; is remarkably drought-resistant and is the longest-lived Medicago species. In addition, M. arborea also has large seeds, disease resistance and morphological traits that could restructure alfalfa [17]. Recently, researchers have attempted the introgression of favorable traits from M. arborea genome to M. sativa subsp. sativa such as drought tolerance [10,17,18]. Hybrids between M. sativa subsp. sativa and M. arborea, which have been given the common name Alborea, have been produced in the USA and Australia [17,18]. A study by Tani et al. [18] showed that at the seedling stage, Alborea plants exhibited higher stem elongation rate, photosynthetic rate, stomatal conductance, and transpiration rate than their parent M. sativa subsp. sativa and M. arborea accessions.
The genetic diversity of crop wild relatives conserved in world’s germplasm banks has not been sufficiently used for crop improvement [19]. According to Dempewolf et al. [20], the use of crop wild relatives in cultivar development have been limited due to: (i) the existence of sufficient variation in the cultivated germplasm of some crops; (ii) differences in ploidy and other hybridization barriers between wild species and cultivated germplasm; (iii) inferiority of wild relatives of crops for desired traits; (iv) insufficient phenotypic and genotypic data on wild relative accessions; (v) inadequate human or financial resources to carry out the necessary research and development.
In this article we studied the phenotypic diversity and productivity of a diverse alfalfa panel composed by cultivars, landraces and wild relatives, all of them with putative drought tolerance [10], in two Mediterranean environments (central Chile and Southern Australia). The objectives were (a) to compare morpho-physiological traits of four M. sativa subspecies (M. sativa subsp. sativa, M. sativa subsp. caerulea, M. sativa subsp. varia and M. alborea); (b) to evaluate the broad sense heritability (H2) of traits; (c) to evaluate their productivity in dryland Mediterranean conditions. The ultimate goal is to select candidate materials to be incorporated into a selection program for drought tolerant varieties for the Mediterranean zone of Chile and Australia. The accessions and lines described in this paper are available from the Australian Pastures Genebank, (https://apg.pir.sa.gov.au/gringlobal/search.aspx, accessed on 21 April 2021) and further information on this project is available at www.projectwebsite, accessed on 21 April 2021.

2. Results

2.1. Plant Survival, Canopy Traits and Forage Yield

The alfalfa diverse panel included four M. sativa subspecies (M. sativa subsp. sativa, M. sativa hybr. (Alborea), M. sativa subsp. caerulea, and M. sativa subsp. varia; hereafter sativa, alborea, caerulea, and varia, respectively; Supplementary Table S1). Within the sativa pool, 20 out of 51 phenotypically characterized accessions were cultivars originated in China, Kazakhstan, Australia, and the USA. Two cultivars belonging to varia were also included in the diversity panel. The remaining accessions were landraces (36) and pre-bred lines (11) developed by crossing cultivars of sativa with alfalfa wild relatives (Supplementary Table S1).
The phenotypic characterization was performed over three years in South Australia (2016–2019) and Chile (2017–2020). The first six months was considered the establishment year and therefore was not included in the analysis. In terms of growing seasons, data from the first and second growing season after the establishment year were considered in both South Australia and Chile. Plant survival was similar among the four subspecies, but there was important variability within subspecies (Table 1). Phenotypic variation was observed among and within subspecies for normalized difference vegetation index (NDVI), plant height, fraction of intercepted photosynthetic active radiation (FIPAR) and annual forage yield; alborea showed the highest NDVI, SL, FIPAR and forage yield values (Table 1). In all phenotypic traits evaluated, excepting plant survival, more than 50% of the total phenotypic variance was attributed to the alfalfa accessions component (Table 1); broad sense heritability ranged between 0.53 and 0.85. Plant survival was highly affected by the environmental condition; it showed a H2 value of 0.10.
Principal component analysis indicated that the first two principal components (PC) accounted for 84.2% of total variance; PC1 and PC2 explained 71.3 and 12.9% of the variability, respectively (Figure 1). Each variable, except forage yield in the first growing season (2018–2019) and plant survival, contributed 13% in average to PC1. Whereas the traits with greater contribution to PC2 were plant survival (74.1%) and forage yield during the first growing season (12.6%). The forage yield in the first and second growing seasons were significantly correlated with all phenotypic traits evaluated (p < 0.05) with r-values ranging between 0.27 (plant survival) and 0.85 (FIPAR) (Figure 1B).
Fall dormancy, taxon, and breeding status determined the agronomic performance of alfalfa accessions. To observe the fall dormancy effects, the alfalfa accessions were assigned to three categories of fall dormancy: dormant type (fall dormancy 0–3), semidormant type (fall dormancy 4–7), and nondormant type (fall dormancy 8–10). The dormant type showed the lowest PC1 values and the highest PC2 (Figure 1A), which were associated with lower and higher values of forage yield and plant survival, respectively (Figure 1A). The nondormant type exhibited the highest PC1 values and the lowest PC2. In both categories, there were some outstanding alfalfa accessions. For instance, accessions APG45675 and CTA008 (dormant type), reached middle-high values of PC1. Accessions CTA003, APG45669 and APG 58574 (nondormant type) reached middle-high values of PC2, which is associated with high annual forage yield and high plant survival. Semidormant accessions exhibited PC1 and PC2 values that were between those of dormant and nondormant types.
As mentioned before, PC1 mostly described agronomic performance and PC2 plant survival. The caerulea accessions showed the lowest agronomic performance in comparison to other subspecies. Sativa accessions were distributed across the entire PC1, which reflects a high level of genetic diversity for the evaluated traits. The varia and alborea accessions reached PC1 values within the sativa pool (Supplementary Figure S2A). In regard to breeding status, some landraces (APG6567, APG45669, APG40234, AltaSierra12) and pre-bred lines (CTA003, CTA004, CTA008) exhibited similar and even better agronomic performance than some cultivars (Supplementary Figure S2B).

2.2. Morphological Traits

Alfalfa subspecies exhibited similar morphological expression, excluding caerulea accessions, which showed the lowest values of StemDM, LeavesDM, LSratio, SLA, Leaf Size, StemD and StemL (Table 2). The sativa pool showed the broadest range of variation in all morphological traits evaluated. Excluding caerulea, alborea accessions showed between 10 and 15% higher values of StemDM, LeavesDM, LSratio, SLA, Leaf Size, StemD and StemL than the other subspecies (Table 2). Excluding StemD and StemL, morphological traits exhibited low genetic control; broad sense heritability varied between 0.10 (LeavesDM) and 0.47 (StemDM).
Two-way cluster analysis (heatmap) identified two groups for morphological traits. One of them included only LSratio. All other morphological traits were clustered together. Three clusters were found for alfalfa accessions (Figure 2). At accession level, Cluster I, included caerulea accessions and two sativa landraces (APG45671 and APG16453). This cluster was characterized for showing the lowest values of LeavesDM, StemDM, StemL, StemD, Leaf Size and SLA, but high values of LSratio. On the other hand, Cluster III was characterized for including three subspecies (sativa, varia and alborea) and accessions with different breeding status. Cluster III included cultivars (WL903HQ, Sardi10, Darkhan90, and Genesis), pre-bred lines (CTA004 and CTA003) and landraces (APG45680, AltaSierra5, Mediterranea, and APG58574). Alfalfa accessions included in Cluster III were characterized for exhibiting morphological characteristics totally opposite to Cluster I (Figure 2). Morphological traits were all significantly correlated among them. The LSratio was the unique morphological trait negatively correlated with each of the other traits (Figure 3).

2.3. Forage Yield in Drought-Prone Environments in Chile and Australia

Annual forage yield data corresponds to the first and second growing seasons after the establishment year, in both Australia and Chile. In Australia, during the first growing season (2017–2018), forage yield under irrigated condition was close to 35% higher than rainfed condition. However, during the second growing season (2018–2019), the forage yield in irrigated condition was almost three times higher than that observed in rainfed condition. In Chile, the average forage yield was 7.2 and 8.9 Mg ha−1 during the first and second growing season, respectively (Figure 4A). In all environments, except in Au18Rf, the forage yield showed broad range of variation among alfalfa accessions (Figure 4A). The GGE analysis accounted for 86.1% for the total variance in forage yield and the GGE biplot analyzed the Accession x Environment (AxE) interaction (Figure 4B). Environments were clustered in two polygons or mega-environments (ME). In ME including Au18Rf, Au19Ir, Au18Ir, Au19Rf, and Cl19Rf pre-bred lines CTA002 and CTA003 showed greater forage yield. In the ME including Cl18Rf, the landrace APG6567 expressed the highest forage yield.
The total amount of water added to each experimental site in every growing season was estimated (rainfall + irrigation). A linear relationship was found between annual forage yield and annual water added to the experiments (R2 = 0.60, p < 0.001; Figure 5).

3. Discussion

This work confirms the existence of a broad phenotypic diversity among accessions of the M. sativa-falcata complex [11,12,14,21,22,23,24], but it also adds new evidence of the genetic contribution of alfalfa wild relatives for improving agronomic performance of cultivated alfalfa in Mediterranean drought-prone environments. Furthermore, new insights about the outstanding genetic contribution of Medicago arborea to the M. sativa-falcata complex were observed at field condition in pre-bred lines.
Fall dormancy is a critical trait for successful cultivation of alfalfa across the world [25,26] since it contributes to alfalfa adaptation/survival in harsh winter environments. In this work, nondormant accessions exhibited higher forage yield than dormant (Figure 1A), which is coincident with results reported in other works [25,26]. In environments with mild temperatures during autumn–winter, nondormant alfalfa genotypes produce more herbage in autumn, resume shoot growth earlier in spring, and initiate shoot regrowth quickly after harvest in summer [27]. As expected, the dormant accessions showed higher plant survival than the nondormant accessions. However, there were some nondormant (CTA003) and semidormant (APG6567) accessions with similar plant survival to dormant alfalfa accessions (Figure 1A). It is important to highlight the dormant accessions APG45675 and CTA008, which expressed similar agronomic performance to nondormant alfalfa (Figure 1A). Accession APG45675 is a sativa landrace grown by farmers in the Bolivian Altiplano. Whereas, accession CTA008 is a varia accession bred from the cross between a falcata wild material originally from Kazakhstan and the Australian cultivar Sardi-Grazer [10]. These two dormant accessions could be used for breeding high-yielding alfalfa cultivars in cold environments [26].

3.1. Phenotypic Diversity

A large body of literature supports the abundant genetic diversity within the M. sativa-falcata complex [10,11,13,14,22,23,24,26,28]. However, genetic diversity used for breeding has been mostly based on cultivated populations. The use of wild populations in breeding programs has been fairly limited [29]. In this work, accessions CTA008 and CTA004, which are varia pre bred-lines, represent the great contribution of the primary gene pool on sativa genome, because it exhibited one of the highest plant survivals with suitable agronomic performance in Mediterranean drought-prone environments (Supplementary Figure S2B).
In the last decades, great efforts have been undertaken in order to produce hybrids between M. sativa and M. arborea, known as the alborea hybrids, which are expected to perform better in drought-prone environments [18]. In this agronomic study, the pre-bred line CTA003, which corresponds to an alborea hybrid derived from the cross between M. sativa subsp. sativa cv. Genesis and M. arborea [10], showed one of the highest DM yields and plant survival (Figure 1). Additionally, the accession CTA003 exhibited similar DM yield to cv. Genesis but with higher plant survival in Mediterranean drought-prone environment. This insight could be considered as empirical evidence of the drought tolerance improvement transferred from M. arborea to M. sativa. Tani et al. [18], showed that at seedling stage, alborea hybrids exhibited higher stem elongation rate, photosynthetic rate, stomatal conductance, and transpiration rate (E) than sativa and arborea accessions under drought stress, indicating that the hybrids correspond to a drought tolerance ideotype [30,31,32,33,34,35].
Morphological characterization revealed that leaf and stem traits are highly correlated within the M. sativa-falcata complex (Figure 2 and Figure 3). Additionally, ShootDM (an estimator of forage yield) forage yield was significantly associated with SLA, leaf size and StemD, which coincide with results of other works [36]. From an agronomic point of view, the alfalfa LSratio has been used as a positive estimator of forage quality, owing to the greater quality of leaves relative to stems [32,37,38]. The LSratio performed as a differentiator trait of the agronomic categories identified by the cluster analysis. For instance, Cluster I, which grouped low-yielding accessions (Figure 2) and small-sized plants (lowest values of LeavesDM, StemDM, ShootDM, LeafSize, StemD and StemL), showed the highest values of LSratio. In contrast, Cluster III, which grouped high-yielding materials and large-sized plants, showed lower values of LSratio. Clusters I and III showed LSratio values of 0.99 and 2.0 on average, respectively. In others works, high forage quality alfalfa cultivars or populations had shown LSratio near to one [37,38,39]. Thus, pre-bred lines (CTA004 and CTA003) and landraces (APG58574, APG45680, Mediterranea, and Altasierra5) grouped in Cluster III could putatively be considered as high forage quality materials.
The relationship between LSratio and forage yield has been scarcely studied. In this work, the LSratio was negatively correlated with all measured morphological traits including ShootDM (Figure 3). Based on this association, a negative trade-off between forage yield and LSratio could be expected. Recently, Annicchiarico [40] pointed out that selection for high-yielding alfalfa populations in Italy did not have any impact on LSratio. However, the range of variation in LSratio of those materials was between 0.89 and 1.0.
Alfalfa cultivars grown in optimal conditions have shown SLA values ranging between 250 and 460 cm2 g−1 [41], which is in agreement with the range of variation observed in this work, excluding the values of the diploid materials (caerulea; Table 2). SLA is a morphological trait strongly associated with alfalfa drought tolerance [36,42]. Erice et al. [36] reported genotypic differences in SLA expression; drought tolerant accessions reduced their SLA under drought conditions, whereas the sensitive accession did not modify it. Mickky et al. [42] also observed the reduction in SLA under drought condition, which was associated with increased leaf thickness and reduced phloem and xylem area. All these modifications could be considered as water saving physiological mechanisms. In this work, accessions belonging to varia and alborea exhibited higher values of SLA, which would support the putative drought tolerance associated with these two subspecies [10,14,17,18]. Accessions with higher values of SLA (344 cm2 g−1, in average) were landraces APG58574 (varia), APG40234 (sativa), and Altasierra5 (sativa); pre-bred lines CTA011, CTA008, CTA004 (varia), and CTA002 (alborea); cultivars Zhungeer (varia) and Genesis (sativa).

3.2. Forage Yield and Accession by Environment Interaction (AxE)

Forage yield in Mediterranean drought-prone environments varied broadly among accessions and across years and locations (Figure 4A). In Australia, irrigation had a large impact on annual forage production, particularly in the third growing season (Figure 4A and Figure 5), where annual forage yield was three-fold higher in irrigated relative to rainfed condition. This could be explained by the rainfall conditions during the 2019 growing season, since South Australia recorded the driest year in more than a decade, with 170 mm less than the long-term average (528.3 mm). Similar response to irrigation was reported by Li and Su [43] in alfalfa grown under a gradient of water treatments in Northern China. Additionally, forage yield throughout locations and years was significantly related to the amount of rainfall or water added to the experiments (Figure 5), which is the typical agronomic responses of crops grown in Mediterranean environments [8,44,45]. The GGE analysis accounted for AxE interaction. Two mega-environments were identified (Figure 4B). However, one of them included all environments, except Cl18Rf. This information is relevant because the Chilean and Australian environments are largely different in their edaphic condition. Therefore, our results suggest that climate has a greater weight on phenotypic expression than soil. This supports the idea that cultivars developed with specific adaptability to Mediterranean environments of Chile or Australia, will perform well in other Mediterranean regions of the world. Recently, Annicchiarico [40] demonstrated the importance of AxE on alfalfa cultivars development. He pointed out that genetic gains on forage yield are significantly higher when cultivars are developed with specific adaptability to target environments. Our results agree with this statement, since the accessions characterized in this work were selected based on their history of drought tolerance or specific adaptability to drought-prone environments [10]. For instance, the landrace APG6567 (varia), was top yielding in Chile because it was collected in arid zones of Spain [10]. In the same way, the alborea hybrids (CTA002 and CTA003), which exhibited the highest forage yield in Australia (2018 and 2019 growing seasons) and Chile (2019 growing season).
From a plant breeding point of view, the inheritance of morphological and agronomic traits evaluated in this alfalfa diverse panel was similar to that reported in other works. For instance, forage yield heritability values around 0.5 has been also reported by Julier et al. [46] and Ray et al. [47]. The LSratio heritability values were similar to values reported by Guines et al. [48]. However, we found a broad range of plant height heritability values in literature; values of 0.09, 0.29, and 0.45 have been reported [46,47,48]. Differences in heritability among different studies may well arise from differences in the population type and size and/or differences in experiment error size. Furthermore, all heritability values used to compare our results were estimated in alfalfa breeding populations.
In conclusion, the wild and pre-bred alfalfa germplasm studied in this work increased the morphological diversity in the M. sativa-falcata complex and allowed to find candidate materials to be incorporated into a selection program for drought tolerant varieties for the Mediterranean zone of Chile and Australia. Additionally, this work revealed that new breeding lines derived from crosses between cultivated alfalfa and wild relatives belonging to the primary (M. sativa-falcata complex) and tertiary (M. arborea) gene pool achieve outstanding agronomic performance in drought-prone environments.

4. Materials and Methods

4.1. Plant Material

A set of 70 alfalfa accessions with putative drought and salinity tolerance was identified using curator knowledge and acquired through the CWR-alfalfa research project supported by Crop Trust (https://www.cwrdiversity.org/partnership/alfalfa-pre-breeding-project-2/, accessed on 21 April 2021). The set included landraces, cultivars and advanced genetic lines (here after accessions) originated from Kazakhstan, Azerbaijan, Spain, Australia, USA, and Chile (Table S1; [10]).

4.2. Experimental Sites and Plants Establishment

In Chile, seeds of 70 accessions (Supplementary Table S1) were sown in 200-holes germination trays containing peat moss as substrate (Kekkila, Finland). Substrate was daily irrigated and periodically fertilized with a solution of 1.1 g L−1 Phostrogen (Bayer, Cambridge, UK). Seedlings were inoculated with a suspension of Sinorhizobium meliloti (strain WSM2141) and grown for two months under greenhouse condition. One week before transplanting, seedlings were moved to a shelter for hardening. On August 2017, seedlings were transplanted at the Cauquenes Research Station of INIA-Chile (35°57′ S; 72°19′ W). Seedlings were arranged in plots of five rows of 2.5 m length separated by 40 cm among them (125 plants per plot). Before transplanting, the soil was prepared using a chisel plow and disc harrows. Plots received 200 kg ha−1 of triple phosphate (45% P2O5), 2000 kg ha−1 of CaCO3, 100 kg ha−1 of potassium sulphate (50% K2O and 54% SO4) and 20 kg ha−1 of boronatrocalcite (11% B). The experiment was arranged in an α-lattice design with three replicates. Each replicate had five incomplete blocks including 14 accessions. The soil is classified as Ultic Palexeralfs, with a pH of 5.7 (1:2.5 in water, 0–20 cm), organic matter content of 2.7% and available N, P, and K content at the top 20 cm of 772 18, 17, and 250 mg kg−1, respectively. The experiment was managed under rainfed condition with 5 or 6 months of negative water balance (Evapotranspiration (ET) > rainfall) during a growing season (Supplementary Figure S1).
In Australia, 30 accessions (Supplementary Table S1) were sown on 26 August 2016 in 1.2 × 2.8 m plots (3.36 m2) at sowing rate of 7 kg ha−1. Seeds were inoculated with Group AL rhizobium prior to sowing. The experiment was established with two blocks at the Waite Research Institute of University of Adelaide (34°58′ S, 138°38′ E). The blocks differed in their water regime, with one managed under rainfed conditions whilst the other block was irrigated. Similar to Chile, the rainfed experiment included 5 or 6 months with negative water balance (Evapotranspiration (ET) > rainfall) during a growing season (Supplementary Figure S1).
Plots were arranged in an augmented (unreplicated) row and column design. The experiments included 12 rows with eight accessions each (columns). The sativa cultivar SARDI 7 Series 2 was grown once in every eight accessions as a highly repeated check to remove spatial variation. The soil was fine sandy loam of the non-sodic Urrbrae series. Soil pH (in CaCl2) was 6.0 with negligible calcium carbonate [49].
The irrigated experiment had subsurface drip irrigation, with two lines running 50 cm apart, 20 cm beneath each plot, and drip intervals of 50 cm. Irrigation was applied monthly at 30 mm between November and May, increasing to twice per week when temperature was above 35 °C.

4.3. Phenotypic Characterization

In Chile, several morphological, physiological and agronomic traits were evaluated, which are described in the following sections. In Australia, only forage yield and plant height were determined.

4.4. Plant Survival and Forage Yield

Plant survival was evaluated by counting the number of surviving plants per plot at the beginning of the fourth growing season (spring 2020) and expressing the result as a percentage of the establishment density (125 plants per plot).
Forage yield was evaluated by cutting 1 m2 of each plot at 5 cm aboveground. Samples were oven-dried with forced air ventilation at 70 °C until reaching constant weight for dry matter determination. Two and three biomass cuts were performed during 2018 and 2019 growing seasons, respectively. The dates of cuts were 13 September 2018, 21 November 2018, 3 September 2019, 24 October 2019, and 19 December 2019. Annual forage yield was determined adding the cuts per years. Plant height (PH) was evaluated in every plot one day before forage yield determination with a 1-m ruler.

4.5. Morphological Traits

Five randomly selected stems per plot were collected. Stem length was measured with a ruler and stem diameter was determined in the middle section of the stem with a digital caliper. Five fully expanded trifoliated leaves were selected for leaf area (LA), petiole length (PL) and dry weight of leaves determination. LA and PL were measured with leaf area meter device (MK2, Delta-T Devices, Cambridge, UK) and digital caliper, respectively. Then, leaves were collected in paper bags for drying. Specific leaf area (SLA) was determined as the relationship between LA and leaves dry weight. Finally, leaves (LeavesDM), stem (StemDM), and shoot dry matter (ShootDM = StemDM + LeavesDM), were determined. Additionally, the LeavesDM:StemDM ratio was calculated. Dry matter contents of leaves and stems were determined after drying of samples at 70 °C in forced air ventilation oven until they reached a constant weight. Morphological characterization was performed at the beginning of the fourth growing season, when plants expressed their higher vigor (28 October 2020).

4.6. Canopy Characterization

The fraction of intercepted photosynthetic active radiation (FIPAR) was evaluated in each plot with a ceptometer, which includes a one meter long probe with 64 PAR sensors and a BF5 reference PAR sensor (SunScan canopy analyzer); three determinations at the bottom of each plot were taken, placing the ceptometer in parallel to the crop row. FIPAR measurements were taken one day before the forage yield determination between 11:00 and 16:00 h, in 2018 and 2019 growing seasons. The NDVI was determined simultaneously with FIPAR measurements, with a handheld spectroradiometer (GreenSeeker, Trimble, Sunnyvale, CA, USA). Measurements were performed on clear days before the sampling for forage yield, passing the sensor 60 cm above the top of the canopy. The NDVI corresponds to the differences between the reflectance (R) in the near-infrared (760 nm) and red (660 nm) band, and is calculated as: (R760 − R660)/(R760 + R660).

4.7. Statistical Analysis

A phenotypic linear mixed model including correlated error for accounting fine-scale spatial variation among experimental units was implemented in the Chilean and Australian data sets. Both experiments were arranged in rows and columns for modeling spatial variation. The model was implemented to estimate the variance components and the best linear unbiased prediction (BLUP) using the restricted maximum likelihood method within the ASReml-R package [50] in R software (https://www.r-project.org/, accessed on 21 April 2021) using the following equation:
Y i k = μ + r k + g i + ε i k
where Yik is the phenotypic value of ith alfalfa accession (g) in the kth replicate (r), μ is the overall population mean, r is the fixed effect of replicate, g is the random effect of the alfalfa accession ~IDD(0,σ_g^2), and ε the random experimental error ~IDD(0,σ_ε^2). Year effect and its interactions with other terms were not considered in the model. The first-order autoregressive anisotropic covariance structure (AR1 x AR1) was used to model spatial variation in row and column directions. This model is reportedly adequate to account for spatial variation in yield trials of cultivars and breeding lines [51,52]. For the unreplicated Australian experiments, the standard check information (Sardi 7) was used for calculating variance components and BLUPs in accordance with Piepho and Williams [53]. In brief, the structure of the model is similar to the model described above, but the alfalfa accession effect (g) was simultaneously modeled as fixed and random effects when standard check and other accessions were adjusted, respectively.
The variance components of the Chilean data sets were used to estimate the broad-sense heritability (H2) on a plot-mean basis, which was calculated as follows:
H 2 = σ g 2 σ g 2 + σ ε 2 r
Predicted means were used to perform multivariate analyses. The relationship among phenotypic traits and alfalfa accessions and the effect of fall dormancy, taxon, and status of breeding on agronomic performance were studied with principal component analysis (PCA) using the packages FactoMinR [54] and factoextra [55] in R software. The morphological diversity expressed by the Medicago subspecies was explored with a two-way clustering analysis using heatmaply package in R software [56].
The annual forage yields of alfalfa accessions evaluated in Chile and Australia were used to explore the accession by environment interaction (AxE). Annual forage yield was the total throughout a growing season. In Chile and Australia, the growing season was extended from spring (September) to autumn (March). Two growing seasons were considered in this analysis, the first and second after the establishment year. Full combinations of locations and growing seasons were considered as single environment. Therefore, the analysis was performed with six environments: rainfed experiment in Australia in 2017–2018 (Au18Rf) and 2018–2019 (Au19Rf) growing seasons, irrigated experiment in Australia in 2017–2018 (Au18Ir) and 2018–2019 (Au19Ir) growing seasons, and rainfed experiment in Chile in 2018–2019 and 2019–2020 growing seasons (Cl18Rf and Cl19Rf). The AxE interaction was explored with the GGE biplot method using the GGEbiplot package in R software [57,58].

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/plants10050862/s1, Table S1: Passport information of the alfalfa diversity panel and sites where accessions were evaluated. Figure S1: Water balance in Chile (A) and Australia (B) of growing seasons in which the alfalfa diverse panel was cultivated. Reference evapotranspiration (ET) was calculated according to FAO standards; monthly accumulated rainfall (Rainfall). Meteorological data in Chile and Australia were obtained from Agromet-INIA (https://agrometeorologia.cl, accessed on 21 April 2021) and the Australian Bureau of Meteorology (http://www.bom.gov.au/?ref=hdr, accessed on 21 April 2021), respectively. Figure S2: Biplot of the first two principal components (PC1 and PC2) for the principal component analysis of nine traits evaluated in the alfalfa diversity panel in Chile. The traits are forage yield (FY), in first (FY2018) and second (FY2019) growing seasons, plant survival (Survival), Normalized Difference Vegetation Index (NDVI), fraction of intercepted photosynthetic active radiation (FIPAR), and plant height (PH) evaluated in 2018 and 2019 growing seasons. Biplots are highlighting the effect of taxon (A) and status of breeding (B).

Author Contributions

C.O., A.H. and A.d.P. planed the experimental work and obtained the funding for this research. S.E., V.B. and M.G. did most of the field work and experimental management in Chile. S.E., L.I. and A.d.P. performed the phenotyping in Chile. L.I. performed the data analysis and prepared the first draft. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Crop Trust GS19001and FONDECYT 1180821.

Institutional Review Board Statement

This work did not include animals and humans. Please exclude this section.

Informed Consent Statement

This work did not include animals and humans. Please exclude this section.

Data Availability Statement

We don’t have a link. However, if somebody needs our data set, they can contact to corresponding authors.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Giorgi, F.; Lionello, P. Climate change projections for the Mediterranean region. Glob. Planet. Chang. 2008, 63, 90–104. [Google Scholar] [CrossRef]
  2. Del Pozo, A.; Brunel-Saldias, N.; Engler, A.; Ortega-Farias, S.; Acevedo-Opazo, C.; Lobos, G.A.; Jara-Rojas, R.; Molina-Montenegro, M.A. Climate Change Impacts and Adaptation Strategies of Agriculture in Mediterranean-Climate Regions (MCRs). Sustainability 2019, 11, 276. [Google Scholar] [CrossRef] [Green Version]
  3. Montecinos, A.; Aceituno, P. Seasonality of the ENSO-Related Rainfall Variability in Central Chile and Associated Circulation Anomalies. J. Clim. 2003, 16, 281–296. [Google Scholar] [CrossRef] [Green Version]
  4. Garreaud, R.D.; Boisier, J.P.; Rondanelli, R.; Montecinos, A.; Sepúlveda, H.H.; Veloso-Aguila, D. The Central Chile Mega Drought (2010–2018): A climate dynamics perspective. Int. J. Clim. 2020, 40, 421–439. [Google Scholar] [CrossRef]
  5. Moore, A.D.; Ghahramani, A. Climate change and broadacre livestock production across southern Australia. 1. Impacts of climate change on pasture and livestock productivity, and on sustainable levels of profitability. Glob. Chang. Biol. 2013, 19, 1440–1455. [Google Scholar] [CrossRef]
  6. Dear, B.S.; Moore, G.A.; Hughes, S.J. Adaptation and potential contribution of temperate perennial legumes to the southern Australian wheatbelt: A review. Aust. J. Exp. Agric. 2003, 43, 1–18. [Google Scholar] [CrossRef]
  7. Porqueddu, C.; Ates, S.; Louhaichi, M.; Kyriazopoulos, A.P.; Moreno, G.; del Pozo, A.; Ovalle, C.; Ewing, M.A.; Nichols, P.G.H. Grasslands in ‘Old World’ and ‘New World’ Mediterranean-climate zones: Past trends, current status and future research priorities. Grass Forage Sci. 2016, 71, 1–35. [Google Scholar] [CrossRef]
  8. Del Pozo, A.; Ovalle, C.; Espinoza, S.; Barahona, V.; Gerding, M.; Humphries, A. Water relations and use-efficiency, plant survival and productivity of nine alfalfa (Medicago sativa L.) cultivars in dryland Mediterranean conditions. Eur. J. Agron. 2017, 84, 16–22. [Google Scholar] [CrossRef]
  9. Muller, M.-H.; Poncet, C.; Prosperi, J.M.; Santoni, S.; Ronfort, J. Domestication history in the Medicago sativa species complex: Inferences from nuclear sequence polymorphism. Mol. Ecol. 2006, 15, 1589–1602. [Google Scholar] [CrossRef]
  10. Humphries, A.; Ovalle, C.; Hughes, S.; del Pozo, A.; Inostroza, L.; Barahona, V.; Yu, L.; Yerzhanova, S.; Rowe, T.; Hill, J.; et al. Characterization, preliminary evaluation and pre-breeding of diverse alfalfa crop wild relatives originating from drought-stressed environments. Crop Sci. 2020, 61, 69–88. [Google Scholar] [CrossRef]
  11. İlhan, D.; Li, X.; Brummer, E.C.; Şakiroğlu, M. Genetic Diversity and Population Structure of Tetraploid Accessions of the Medicago sativa–falcata Complex. Crop Sci. 2016, 56, 1146–1156. [Google Scholar] [CrossRef]
  12. Sakiroglu, M.; Moore, K.J.; Brummer, E.C. Variation in Biomass Yield, Cell Wall Components, and Agronomic Traits in a Broad Range of Diploid Alfalfa Accessions. Crop Sci. 2011, 51, 1956–1964. [Google Scholar] [CrossRef] [Green Version]
  13. Basbag, M.; Aydin, A.; Sakiroglu, M. Evaluating Agronomic Performance and Investigating Molecular Structure of Drought and Heat Tolerant Wild Alfalfa (Medicago sativa L.) Collection from the Southeastern Turkey. Biochem. Genet. 2017, 55, 63–76. [Google Scholar] [CrossRef] [PubMed]
  14. Şakiroğlu, M.; İlhan, D. Medicago sativa species complex: Re-visiting the century old problem in the light of molecular tools. Crop Sci. 2020, 61, 827–838. [Google Scholar] [CrossRef]
  15. Ray, I.M.; Segovia-Lerma, A.; Murray, L.W. Diallel analysis of carbon isotope discrimination and its association with forage yield among nine historically recognized alfalfa germplasms. Crop Sci. 2004, 44, 1970–1975. [Google Scholar] [CrossRef]
  16. Hanson, A.; Xu, L.; Boe, A.; Johnson, P.S.; Gates, R.N.; Wu, Y. Identification and characterization of drought tolerant alfalfa (Medicago sativa subsp. falcata) germplasm. Proc. S. D. Acad. Sci. 2016, 94, 263–272. [Google Scholar]
  17. Bingham, E.; Armour, D.; Irwin, J. The Hybridization Barrier between Herbaceous Medicago sativa and Woody M. arborea Is Weakened by Selection of Seed Parents. Plants 2013, 2, 343–353. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Tani, E.; Chronopoulou, E.G.; Labrou, N.E.; Sarri, E.; Goufa, Μ.; Vaharidi, X.; Tornesaki, A.; Psychogiou, M.; Bebeli, P.J.; Abraham, Ε.M. Growth, Physiological, Biochemical, and Transcriptional Responses to Drought Stress in Seedlings of Medicago sativa L., Medicago arborea L. and Their Hybrid (Alborea). Agronomy 2019, 9, 38. [Google Scholar] [CrossRef] [Green Version]
  19. Byrne, P.F.; Volk, G.M.; Gardner, C.; Gore, M.A.; Simon, P.W.; Smith, S. Sustaining the Future of Plant Breeding: The Critical Role of the USDA-ARS National Plant Germplasm System. Crop Sci. 2018, 58, 451–468. [Google Scholar] [CrossRef] [Green Version]
  20. Dempewolf, H.; Baute, G.; Anderson, J.; Kilian, B.; Smith, C.; Guarino, L. Past and Future Use of Wild Relatives in Crop Breeding. Crop Sci. 2017, 57, 1070–1082. [Google Scholar] [CrossRef]
  21. Sakiroglu, M.; Charles Brummer, E. Presence of phylogeographic structure among wild diploid alfalfa accessions (Medicago sativa L. subsp. microcarpa Urb.) with evidence of the center of origin. Genet. Resour. Crop Evol. 2013, 60, 23–31. [Google Scholar] [CrossRef]
  22. Zhang, T.; Kesoju, S.; Greene, S.L.; Fransen, S.; Hu, J.; Yu, L.-X. Genetic diversity and phenotypic variation for drought resistance in alfalfa (Medicago sativa L.) germplasm collected for drought tolerance. Genet. Resour. Crop Evol. 2018, 65, 471–484. [Google Scholar] [CrossRef]
  23. Qiang, H.; Chen, Z.; Zhang, Z.; Wang, X.; Gao, H.; Wang, Z. Molecular Diversity and Population Structure of a Worldwide Collection of Cultivated Tetraploid Alfalfa (Medicago sativa subsp. sativa L.) Germplasm as Revealed by Microsatellite Markers. PLoS ONE 2015, 10, e0124592. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Şakiroğlu, M.; Doyle, J.J.; Charles Brummer, E. Inferring population structure and genetic diversity of broad range of wild diploid alfalfa (Medicago sativa L.) accessions using SSR markers. Appl. Genet. 2010, 121, 403–415. [Google Scholar] [CrossRef] [PubMed]
  25. Zhang, S.; Shi, Y.; Cheng, N.; Du, H.; Fan, W.; Wang, C. De Novo Characterization of Fall Dormant and Nondormant Alfalfa (Medicago sativa L.) Leaf Transcriptome and Identification of Candidate Genes Related to Fall Dormancy. PLoS ONE 2015, 10, e0122170. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Ariss, J.J.; Vandemark, G.J. Assessment of Genetic Diversity among Nondormant and Semidormant Alfalfa Populations Using Sequence-Related Amplified Polymorphisms. Crop Sci. 2007, 47, 2274–2284. [Google Scholar] [CrossRef]
  27. Cunningham, S.M.; Gana, J.A.; Volenec, J.J.; Teuber, L.R. Winter Hardiness, Root Physiology, and Gene Expression in Successive Fall Dormancy Selections from ‘Mesilla’ and ‘CUF 101’ Alfalfa. Crop Sci. 2001, 41, 1091–1098. [Google Scholar] [CrossRef] [Green Version]
  28. Riday, H.; Brummer, E.C. Heterosis in a Broad Range of Alfalfa Germplasm. Crop Sci. 2005, 45, 8–17. [Google Scholar] [CrossRef]
  29. Annicchiarico, P.; Barrett, B.; Brummer, E.C.; Julier, B.; Marshall, A.H. Achievements and Challenges in Improving Temperate Perennial Forage Legumes. Crit. Rev. Plant Sci. 2015, 34, 327–380. [Google Scholar] [CrossRef]
  30. Inostroza, L.; Acuña, H.; Méndez, J. Multi-physiological-trait selection indices to identify Lotus tenuis genotypes with high dry matter production under drought conditions. Crop Pasture Sci. 2015, 66, 90–99. [Google Scholar] [CrossRef]
  31. Inostroza, L.; Acuña, H. Water use efficiency and associated physiological traits of nine naturalized white clover populations in Chile. Plant Breed. 2010, 129, 700–706. [Google Scholar] [CrossRef]
  32. Inostroza, L.; Ortega-Klose, F.; Vásquez, C.; Wilckens, R. Changes in Root Architecture and Aboveground Traits of Red Clover Cultivars Driven by Breeding to Improve Persistence. Agronomy 2020, 10, 1896. [Google Scholar] [CrossRef]
  33. Inostroza, L.; Acuña, H.; Méndez, J.; Bhakta, M.; Gezan, S.A. Using genome conservation between Lotus japonicus and agronomically important Lotus species for discovering drought tolerance QTLs. Euphytica 2019, 215, 153. [Google Scholar] [CrossRef]
  34. Acuña, H.; Inostroza, L.; Sánchez, M.P.; Tapia, G. Drought-tolerant naturalized populations of Lotus tenuis for constrained environments. Acta Agric. Scand. Sect. B Soil Plant Sci. 2010, 60, 174–181. [Google Scholar] [CrossRef]
  35. Inostroza, L.; Acuña, H.; Tapia, G. Relationships between phenotypic variation in osmotic adjustment, water-use efficiency, and drought tolerance of seven cultivars of Lotus corniculatus L. Chil. J. Agric. Res. 2015, 75, 3–12. [Google Scholar] [CrossRef] [Green Version]
  36. Erice, G.; Louahlia, S.; Irigoyen, J.J.; Sanchez-Diaz, M.; Avice, J.-C. Biomass partitioning, morphology and water status of four alfalfa genotypes submitted to progressive drought and subsequent recovery. J. Plant Physiol. 2010, 167, 114–120. [Google Scholar] [CrossRef] [PubMed]
  37. Biazzi, E.; Nazzicari, N.; Pecetti, L.; Brummer, E.C.; Palmonari, A.; Tava, A.; Annicchiarico, P. Genome-Wide Association Mapping and Genomic Selection for Alfalfa (Medicago sativa) Forage Quality Traits. PLoS ONE 2017, 12, e0169234. [Google Scholar] [CrossRef]
  38. Annicchiarico, P. Alfalfa forage yield and leaf/stem ratio: Narrow-sense heritability, genetic correlation, and parent selection procedures. Euphytica 2015, 205, 409–420. [Google Scholar] [CrossRef]
  39. Grev, A.M.; Wells, M.S.; Catalano, D.N.; Martinson, K.L.; Jungers, J.M.; Sheaffer, C.C. Stem and leaf forage nutritive value and morphology of reduced lignin alfalfa. Agron. J. 2020, 112, 406–417. [Google Scholar] [CrossRef]
  40. Annicchiarico, P. Breeding gain from exploitation of regional adaptation: An Alfalfa case study. Crop Sci. 2020. [Google Scholar] [CrossRef]
  41. Malik, W.; Boote, K.J.; Hoogenboom, G.; Cavero, J.; Dechmi, F. Adapting the CROPGRO Model to Simulate Alfalfa Growth and Yield. Agron. J. 2018, 110, 1777–1790. [Google Scholar] [CrossRef] [Green Version]
  42. Mickky, B.M.; Abbas, M.A.; El-Shhaby, O.A. Alterations in photosynthetic capacity and morpho-histological features of leaf in alfalfa plants subjected to water deficit-stress in different soil types. Indian J. Plant Physiol. 2018, 23, 426–443. [Google Scholar] [CrossRef]
  43. Li, Y.; Su, D. Alfalfa Water Use and Yield under Different Sprinkler Irrigation Regimes in North Arid Regions of China. Sustainability 2017, 9, 1380. [Google Scholar] [CrossRef] [Green Version]
  44. Richards, R.A.; Rebetzke, G.J.; Condon, A.G.; van Herwaarden, A.F. Breeding Opportunities for Increasing the Efficiency of Water Use and Crop Yield in Temperate Cereals. Crop Sci. 2002, 42, 111–121. [Google Scholar] [CrossRef]
  45. Condon, A.G.; Richards, R.A.; Rebetzke, G.J.; Farquhar, G.D. Breeding for high water-use efficiency. J. Exp. Bot. 2004, 55, 2447–2460. [Google Scholar] [CrossRef] [Green Version]
  46. Julier, B.; Huyghe, C.; Ecalle, C. Within- and Among-Cultivar Genetic Variation in Alfalfa: Forage Quality, Morphology, and Yield. Crop Sci. 2000, 40, 365–369. [Google Scholar] [CrossRef] [Green Version]
  47. Ray, I.M.; Townsend, M.S.; Muncy, C.M. Heritabilities and Interrelationships of Water-Use Efficiency and Agronomic Traits in Irrigated Alfalfa. Crop Sci. 1999, 39, 1088–1092. [Google Scholar] [CrossRef]
  48. Guines, F.; Julier, B.; Ecalle, C.; Huyghe, C. Genetic control of quality traits of lucerne (Medicago sativa L.). Aust. J. Agric. Res. 2002, 53, 401–407. [Google Scholar] [CrossRef]
  49. Grace, P.R.; Oades, J.M.; Keith, H.; Hancock, T.W. Trends in wheat yields and soil organic carbon in the Permanent Rotation Trial at the Waite Agricultural Research Institute, South Australia. Aust. J. Exp. Agric. 1995, 35, 857–864. [Google Scholar] [CrossRef]
  50. Butler, D. asreml: Fits the Linear Mixed Model; R Package Version 4.1.0.130; VSNI: Brisbane, Australia, 2020. [Google Scholar]
  51. Sripathi, R.; Conaghan, P.; Grogan, D.; Casler, M.D. Spatial Variability Effects on Precision and Power of Forage Yield Estimation. Crop Sci. 2017, 57, 1383–1393. [Google Scholar] [CrossRef]
  52. Andrade, M.H.M.L.; Fernandes Filho, C.C.; Fernandes, M.O.; Bastos, A.J.R.; Guedes, M.L.; de Marçal, T.S.; Gonçalves, F.M.A.; Pinto, C.A.B.P.; Zotarelli, L. Accounting for spatial trends to increase the selection efficiency in potato breeding. Crop Sci. 2020, 60, 2354–2372. [Google Scholar] [CrossRef]
  53. Piepho, H.-P.; Williams, E.R. Augmented Row–Column Designs for a Small Number of Checks. Agron. J. 2016, 108, 2256–2262. [Google Scholar] [CrossRef]
  54. Lê, S.; Josse, J.; Husson, F. FactoMineR: An R Package for Multivariate Analysis. J. Stat. Softw. 2008, 25, 1–18. [Google Scholar] [CrossRef] [Green Version]
  55. Kassambara, A.; Mundt, F. Factoextra: Extract and Visualize the Results of Multivariate Data Analyses; R Package Version 1.0.7.; 2020; Available online: https://cran.r-project.org/web/packages/factoextra/index.html (accessed on 21 April 2021).
  56. Galili, T.; O’Callaghan, A.; Sidi, J.; Sievert, C. heatmaply: An R package for creating interactive cluster heatmaps for online publishing. Bioinformatics 2018, 34, 1600–1602. [Google Scholar] [CrossRef] [PubMed]
  57. Yan, W.; Kang, M.S.; Ma, B.; Woods, S.; Cornelius, P.L. GGE Biplot vs. AMMI Analysis of Genotype-by-Environment Data. Crop Sci. 2007, 47, 643–653. [Google Scholar] [CrossRef]
  58. Dumble, S. GGEBiplots: GGE Biplots with “ggplot2”; R Package Version 0.1.1. 2017. Available online: https://cran.r-project.org/web/packages/GGEBiplots/index.html (accessed on 21 April 2021).
Plants 10 00862 g001 550
Figure 1. Biplot (A) and correlogram (B) of the first two principal components (PC1 and PC2) for the principal component analysis of nine traits evaluated in the alfalfa diversity panel in Chile. The traits are forage yield (FY), in the first (FY2018) and second (FY2019) growing seasons, plant survival (Survival), Normalized Difference Vegetation Index (NDVI), fraction of intercepted photosynthetic active radiation (FIPAR), and plant height (PH) evaluated in 2018 and 2019 growing seasons. In the correlogram * and *** refer to p < 0.05 and p < 0.001, respectively.
Figure 1. Biplot (A) and correlogram (B) of the first two principal components (PC1 and PC2) for the principal component analysis of nine traits evaluated in the alfalfa diversity panel in Chile. The traits are forage yield (FY), in the first (FY2018) and second (FY2019) growing seasons, plant survival (Survival), Normalized Difference Vegetation Index (NDVI), fraction of intercepted photosynthetic active radiation (FIPAR), and plant height (PH) evaluated in 2018 and 2019 growing seasons. In the correlogram * and *** refer to p < 0.05 and p < 0.001, respectively.
Plants 10 00862 g001
Plants 10 00862 g002 550
Figure 2. Heatmap visualization of two-way cluster analyses for eight morphological traits (column tree) and 70 alfalfa accessions (row tree) evaluated in the second growing season (2019–2020) in Mediterranean drought-prone environment of Chile. Row annotation indicates taxon and fall dormancy. The greater values in phenotypic expression bar indicate greater phenotypic values in a given trait. Morphological traits: shoot dry matter (ShootDM), stem dry matter (StemDM), leaves dry matter (LeavesDM), LeavesDM to StemDM ratio (LSratio), specific leaf area (SLA), stem length (StemL), stem diameter (StemD), and leaf size (LeafSize).
Figure 2. Heatmap visualization of two-way cluster analyses for eight morphological traits (column tree) and 70 alfalfa accessions (row tree) evaluated in the second growing season (2019–2020) in Mediterranean drought-prone environment of Chile. Row annotation indicates taxon and fall dormancy. The greater values in phenotypic expression bar indicate greater phenotypic values in a given trait. Morphological traits: shoot dry matter (ShootDM), stem dry matter (StemDM), leaves dry matter (LeavesDM), LeavesDM to StemDM ratio (LSratio), specific leaf area (SLA), stem length (StemL), stem diameter (StemD), and leaf size (LeafSize).
Plants 10 00862 g002
Plants 10 00862 g003 550
Figure 3. Correlogram of eight morphological traits evaluated in 70 alfalfa accessions in the second growing season (2019–2020) in a Mediterranean drought-prone environment of Chile. In the correlogram * and *** refer to p < 0.05 and p < 0.001, respectively. Morphological traits: shoot dry matter (ShootDM), stem dry matter (StemDM), leaves dry matter (LeavesDM), LeavesDM to StemDM ratio (LSratio), specific leaf area (SLA), stem length (StemL), stem diameter (StemD), and leaf size (LeafSize).
Figure 3. Correlogram of eight morphological traits evaluated in 70 alfalfa accessions in the second growing season (2019–2020) in a Mediterranean drought-prone environment of Chile. In the correlogram * and *** refer to p < 0.05 and p < 0.001, respectively. Morphological traits: shoot dry matter (ShootDM), stem dry matter (StemDM), leaves dry matter (LeavesDM), LeavesDM to StemDM ratio (LSratio), specific leaf area (SLA), stem length (StemL), stem diameter (StemD), and leaf size (LeafSize).
Plants 10 00862 g003
Plants 10 00862 g004 550
Figure 4. Frequency distribution (box plot; A) and GGE biplot (B) of annual forage yield evaluated in 30 alfalfa accessions in drought-prone environments of Chile and Australia during the first and second growing season after the establishment year. GGE biplot corresponds to the first two principal components (PC1 and PC2) for the principal component analysis of annual forage yield evaluated in six environments: rainfed experiment in Australia in 2017–2018 (Au18Rf) and 2018–2019 (Au19Rf) growing seasons, irrigated experiment in Australia in 2017–2018 (Au18Ir) and 2018–2019 (Au19Ir) growing seasons, and rainfed experiment in Chile in 2018–2019 and 2019–2020 growing seasons (Cl18Rf and Cl19Rf).
Figure 4. Frequency distribution (box plot; A) and GGE biplot (B) of annual forage yield evaluated in 30 alfalfa accessions in drought-prone environments of Chile and Australia during the first and second growing season after the establishment year. GGE biplot corresponds to the first two principal components (PC1 and PC2) for the principal component analysis of annual forage yield evaluated in six environments: rainfed experiment in Australia in 2017–2018 (Au18Rf) and 2018–2019 (Au19Rf) growing seasons, irrigated experiment in Australia in 2017–2018 (Au18Ir) and 2018–2019 (Au19Ir) growing seasons, and rainfed experiment in Chile in 2018–2019 and 2019–2020 growing seasons (Cl18Rf and Cl19Rf).
Plants 10 00862 g004
Plants 10 00862 g005 550
Figure 5. Relationship between annual forage yield evaluated in Mediterranean environments of Chile and Australia and amount of water (rainfall or rainfall + irrigation) added to the experimental plots. Environments were rainfed experiment in Australia in 2018 (Au18Rf) and 2019 (Au19Rf) growing seasons, irrigated experiment in Australia in 2018 (Au18Ir) and 2019 (Au19Ir) growing seasons, and rainfed experiment in Chile in 2018 and 2019 growing season (Cl18Rf and Cl19Rf).
Figure 5. Relationship between annual forage yield evaluated in Mediterranean environments of Chile and Australia and amount of water (rainfall or rainfall + irrigation) added to the experimental plots. Environments were rainfed experiment in Australia in 2018 (Au18Rf) and 2019 (Au19Rf) growing seasons, irrigated experiment in Australia in 2018 (Au18Ir) and 2019 (Au19Ir) growing seasons, and rainfed experiment in Chile in 2018 and 2019 growing season (Cl18Rf and Cl19Rf).
Plants 10 00862 g005
Table
Table 1. Mean ± standard deviation and broad sense heritability (H2 ± SE) of plant survival, normalized difference vegetation index (NDVI), plant height, fraction of intercepted photosynthetic active radiation (FIPAR) and dry matter production, evaluated in an alfalfa diversity panel including accessions belonging to four subspecies M. sativa subsp. sativa, M. sativa subsp. caerulea, M. sativa subsp. varia and M. sativa subsp. alborea. Phenotypic characterization was performed during two growing seasons (2018–2019 and 2019–2020) in a Mediterranean drought-prone environment in Chile.
Table 1. Mean ± standard deviation and broad sense heritability (H2 ± SE) of plant survival, normalized difference vegetation index (NDVI), plant height, fraction of intercepted photosynthetic active radiation (FIPAR) and dry matter production, evaluated in an alfalfa diversity panel including accessions belonging to four subspecies M. sativa subsp. sativa, M. sativa subsp. caerulea, M. sativa subsp. varia and M. sativa subsp. alborea. Phenotypic characterization was performed during two growing seasons (2018–2019 and 2019–2020) in a Mediterranean drought-prone environment in Chile.
TaxonNDVIPlant Height (cm)FIPARForage Yield (Mg ha−1)Plant Survival
(%)
2018/192019/202018/192019/202018/192019/202018/192019/20
M.s. hybr. (alborea)
(* n = 2)		0.62 ± 0.110.80 ± 0.0752.9 ± 0.143.3 ± 6.80.88 ± 0.010.87 ± 0.048.88 ± 0.111.65 ± 2.884.2 ± 12.8
M.s. caerulea
(n = 2)		0.35 ± 0.050.32 ± 0.019.3 ± 4.28.3 ± 5.40.035 ± 0.010.13 ± 0.094.10 ± 0.041.25 ± 0.279.3 ± 9.2
M.s. sativa
(n = 50)		0.59 ± 0.120.79 ± 0.0839.8 ± 10.836.3 ± 10.30.69 ± 0.210.79 ± 0.186.67 ± 1.88.83 ± 2.685.2 ± 5.9
M.s. varia
(n = 14) 		0.56 ± 0.100.76 ± 0.1037.6 ± 9.534.0 ± 10.30.68 ± 0.170.72 ± 0.237.15 ± 0.98.75 ± 2.287.4 ± 4.5
H20.85 ± 0.020.67 ± 0.08 0.65 ± 0.070.53 ± 0.090.10 ± 0.01
* n is the number of accessions in each taxon.
Table
Table 2. Mean ± standard deviation and broad sense heritability (H2 ± SE) of eight morphological traits evaluated in an alfalfa diversity panel including accessions belonging to four subspecies M. sativa subsp. sativa, M. sativa subsp. caerulea, M. sativa subsp. varia and M. sativa subsp. alborea. Phenotypic characterization was performed during the second growing season (2019–2020) in Mediterranean drought-prone environment of Chile.
Table 2. Mean ± standard deviation and broad sense heritability (H2 ± SE) of eight morphological traits evaluated in an alfalfa diversity panel including accessions belonging to four subspecies M. sativa subsp. sativa, M. sativa subsp. caerulea, M. sativa subsp. varia and M. sativa subsp. alborea. Phenotypic characterization was performed during the second growing season (2019–2020) in Mediterranean drought-prone environment of Chile.
TaxonShootDM
(g)		StemDM
(g)		LeavesDM
(g)		LSratio
(g g−1)		SLA
(g cm−2)		StemL
(cm)		StemD
(mm)		LeafSize (cm2)
M.s. hybr. (alborea)
(* n = 2)		3.70 ± 0.571.78 ± 0.25 1.93 ± 0.32 1.12 ± 0.06 319.7 ± 70.8 50.4 ± 5.73.05 ± 0.1096.07 ± 9.8
M.s. caerulea
(n = 2)		1.23 ± 0.600.45 ± 0.28 0.78± 0.322.31 ± 0.81 86.3 ± 16.925.1 ± 3.21.69 ± 0.2110.81 ± 6.8
M.s. sativa
(n = 50)		3.06 ± 0.861.35 ± 0.51 1.69 ± 0.39 1.40 ± 0.43 242.5 ± 50.542.3 ± 7.92.71 ± 0.3967.62 ± 22.8
M.s. varia
(n = 14) 		3.18 ± 0.711.36 ± 0.471.82 ± 0.30 1.58 ± 0.58 269.4 ± 64.943.1 ± 7.02.63 ± 0.3179.10 ± 25.4
H20.28 ± 0.100.47 ± 0.120.10 ± 0.010.30 ± 0.010.23 ± 0.070.62 ± 0.070.52 ± 0.060.26 ± 0.02
Morphological traits: shoot dry matter (ShootDM), stem dry matter (StemDM), leaves dry matter (LeavesDM), LeavesDM to StemDM ratio (LSratio), specific leaf area (SLA), stem length (StemL), stem diameter (StemD), and leaf size (LeafSize).
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Inostroza, L.; Espinoza, S.; Barahona, V.; Gerding, M.; Humphries, A.; del Pozo, A.; Ovalle, C. Phenotypic Diversity and Productivity of Medicago sativa Subspecies from Drought-Prone Environments in Mediterranean Type Climates. Plants 2021, 10, 862. https://doi.org/10.3390/plants10050862

AMA Style

Inostroza L, Espinoza S, Barahona V, Gerding M, Humphries A, del Pozo A, Ovalle C. Phenotypic Diversity and Productivity of Medicago sativa Subspecies from Drought-Prone Environments in Mediterranean Type Climates. Plants. 2021; 10(5):862. https://doi.org/10.3390/plants10050862

Chicago/Turabian Style

Inostroza, Luis, Soledad Espinoza, Viviana Barahona, Macarena Gerding, Alan Humphries, Alejandro del Pozo, and Carlos Ovalle. 2021. "Phenotypic Diversity and Productivity of Medicago sativa Subspecies from Drought-Prone Environments in Mediterranean Type Climates" Plants 10, no. 5: 862. https://doi.org/10.3390/plants10050862

APA Style

Inostroza, L., Espinoza, S., Barahona, V., Gerding, M., Humphries, A., del Pozo, A., & Ovalle, C. (2021). Phenotypic Diversity and Productivity of Medicago sativa Subspecies from Drought-Prone Environments in Mediterranean Type Climates. Plants, 10(5), 862. https://doi.org/10.3390/plants10050862

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

Article Metrics

Back to TopTop