Next Article in Journal
Effect of Statin Lipophilicity on the Proliferation of Hepatocellular Carcinoma Cells
Next Article in Special Issue
Combining Historical and Molecular Data to Study Nearly Extinct Native Italian Grey Partridge (Perdix perdix) at the Turn of the Twentieth Century
Previous Article in Journal
Exposure to Waterpipe Smoke Disrupts Erythrocyte Homeostasis of BALB/c Mice
Previous Article in Special Issue
Ovotestis Isolation and Cryopreservation of Nesiohelix samarangae (Oriental Snail) as a Snail Model for Conserving Other Endangered Snail Species
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genetic and Phenotypic Evaluation of European Maize Landraces as a Tool for Conservation and Valorization of Agrobiodiversity

by
Carlotta Balconi
1,*,†,
Agustin Galaretto
2,†,
Rosa Ana Malvar
3,
Stéphane D. Nicolas
2,
Rita Redaelli
1,
Violeta Andjelkovic
4,
Pedro Revilla
3,
Cyril Bauland
2,
Brigitte Gouesnard
5,
Ana Butron
3,
Alessio Torri
1,
Ana Maria Barata
6,
Natalija Kravic
4,
Valérie Combes
2,
Pedro Mendes-Moreira
7,8,
Danela Murariu
9,
Hrvoje Šarčević
10,
Beate Schierscher-Viret
11,
Morgane Vincent
5,
Anne Zanetto
5,
Bettina Kessel
12,
Delphine Madur
2,
Tristan Mary-Huard
2,13,
André Pereira
7,8,
Domnica Daniela Placinta
9,
Alexandre Strigens
14,
Alain Charcosset
2 and
Sandra Goritschnig
15,*
add Show full author list remove Hide full author list
1
CREA—Council for Agricultural Research and Economics, Research Centre for Cereal and Industrial Crops, via Stezzano 24, 24126 Bergamo, Italy
2
INRAE, CNRS, AgroParisTech, GQE—Le Moulon, Université Paris-Saclay, 12 route 128, 91190 Gif-sur-Yvette, France
3
Misión Biológica de Galicia Consejo Superior de Investigaciones Científicas, Pazo de Salcedo Carballeira, 8 Salcedo, 36143 Pontevedra, Spain
4
Maize Research Institute Zemun Polje, 11000 Belgrade, Serbia
5
UMR AGAP Institut, CIRAD, INRAE, Institut Agro, University Montpellier, F-34398 Montpellier, France
6
Banco Português de Germoplasma Vegetal, Quinta de S. José, S.Pedro de Merelim, 4700-859 Braga, Portugal
7
Coimbra School of Agriculture, Polytechnic University of Coimbra (ESAC-IPC), 3045-093 Coimbra, Portugal
8
CERNAS—Research Centre for Natural Resources, Environment and Society, Bencanta, 3045-601 Coimbra, Portugal
9
Suceava Genebank, B-Dul. 1 Mai 17, 720224 Suceava, Romania
10
Faculty of Agriculture, University of Zagreb, Svetošimunska Cesta 25, 10000 Zagreb, Croatia
11
Agroscope, Route de Duillier 60, 1260 Nyon, Switzerland
12
KWS SAAT SE & Co. KGaA, Grimsehlstr. 31, 37574 Einbeck, Germany
13
INRAE, UMR MIA Paris-Saclay, Université Paris-Saclay, AgroParisTech, 91120 Paris, France
14
DSP—Delley Semences et Plantes SA, Route de Portalban 40, 1567 Delley, Switzerland
15
ECPGR, Alliance of Bioversity International and CIAT, Via di San Domenico 1, 00153 Rome, Italy
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Biology 2024, 13(6), 454; https://doi.org/10.3390/biology13060454
Submission received: 15 April 2024 / Revised: 30 May 2024 / Accepted: 4 June 2024 / Published: 19 June 2024

Abstract

:

Simple Summary

Maize is one of the major crops of the world for feed, food, and industrial uses. It originated in Central America and was first introduced into Europe at the end of the 15th century. Due to its adaptability, farmers and breeders across Europe have developed a wide diversity of local maize varieties with different characteristics over the past centuries. Many of these are conserved in genebanks’ seed collections, but little is known about their specific characteristics. Here, we present results obtained by the European Evaluation Network for Maize, a private–public partnership with partners from nine countries aimed at promoting the valorization of maize genetic resources in breeding programs. The work describes the selection and the genetic and phenotypic evaluation of a collection of 626 maize landraces preserved in European genebanks, providing evidence for historic introductions and geographic adaptation. In a world where climate change, rising food prices, and other issues are affecting food security and the environment, the conservation and use of crop diversity is becoming increasingly important. The results of our study will facilitate the use of maize genetic resources in breeding for resilience to climate change, for sustainable agriculture, food security, and food quality.

Abstract

The ECPGR European Evaluation Network (EVA) for Maize involves genebanks, research institutions, and private breeding companies from nine countries focusing on the valorization of maize genetic resources across Europe. This study describes a diverse collection of 626 local landraces and traditional varieties of maize (Zea mays L.) from nine European genebanks, including criteria for selection of the collection and its genetic and phenotypic diversity. High-throughput pool genotyping grouped the landraces into nine genetic groups with a threshold of 0.6 admixture, while 277 accessions were designated admixed and likely to have resulted from previous breeding activities. The grouping correlated well with the geographic origins of the collection, also reflecting the various pathways of introduction of maize to Europe. Phenotypic evaluations of 588 accessions for flowering time and plant architecture in multilocation trials over three years confirmed the great diversity within the collection, although phenotypic clusters only partially correlated with the genetic grouping. The EVA approach promotes conservation of genetic resources and opens an opportunity to increase genetic variability for developing improved varieties and populations for farmers, with better adaptation to specific environments and greater tolerance to various stresses. As such, the EVA maize collection provides valuable sources of diversity for facing climate change due to the varieties’ local adaptation.

Graphical Abstract

1. Introduction

Maize (Zea mays L.) was brought to Europe for the first time in 1494, when Columbus arrived at Seville in Southern Spain [1]. Successive introductions of maize were made in Southwestern Europe. Genotyping analyses [2] have confirmed that maize was first introduced in Southern Spain from the Caribbean area, and shortly after from North America. The most likely hypotheses for this second introduction are that it arrived through the European Atlantic coast [1] and French explorations conducted by Giovanni da Verrazzano. From these two introductions originated two main European maize pools, the Mediterranean and European flints, which were crossed afterwards and selected for adaptation to diverse environments and uses.
Plant genetic resources for food and agriculture (PGRFA) are the backbone of the world’s food security as they provide the necessary genetic diversity for researchers, breeders, and farmers to continuously develop new varieties that can withstand challenges associated with the climate crisis and support sustainable agricultural production. In view of the changing climate, enhanced efforts are needed to conserve existing locally adapted PGRFA diversity and make it accessible for breeding programs through characterization.
The European Cooperative Programme for Plant Genetic Resources (ECPGR) is a pan-European network active since 1980, connecting national programs on PGRFA through expert working groups and other initiatives with the aim to support conservation and sustainable utilization of PGRFA, in line with the Plant Genetic Resources Strategy for Europe [3]. An ECPGR initiative launched in 2019, the European Evaluation Network (EVA) brings together national genebanks and other stakeholders from the public and private sectors into public–private partnerships (PPPs) jointly working on creating and sharing knowledge on PGRFA conserved in European genebanks [4]. Using standardized protocols, crop-specific EVA networks evaluate genebank materials, which are often understudied, in multiple locations across Europe, thereby identifying climate-resilient breeding material and enriching genebank inventories, ultimately promoting their sustainable use. The EVA Maize Network, launched in 2020, currently has partners from nine countries, including eleven genebanks and research institutes and eight breeding companies (https://www.ecpgr.cgiar.org/eva/eva-networks/maize, accessed on 14 April 2024). It builds on existing national networks, connecting members of the ECPGR Maize Working Group with breeding companies involved in the French PPP initiative ProMaïs (http://pro-mais.org/, accessed on 14 April 2024). Network activities are focused on the valorization of maize genetic resources across Europe, with a focus on landraces and traditional varieties [5], which are valuable sources of genetic diversity for facing climate change, due to their local adaptation.
The EVA Maize Network also links to a previous initiative, the European program RESGEN88 (1997–2001, https://www.agap-ge2pop.org/eu-genres-088/, accessed on 14 April 2024). This project aimed to obtain a better knowledge of the maize landraces grown in Europe in the past centuries and now maintained in genebanks, and to define the best conditions for their valorization and exploitation. The RESGEN88 program involved seven countries: the Netherlands, Germany, Greece, France, Portugal, Italy, and Spain and shared a set of 394 local maize populations. Each country defined a core set of landraces representative of their national maize collection. A preliminary evaluation through agronomical trials, collection of morphological descriptors, and molecular analyses via RFLP [6] resulted in the constitution of a European Maize Landraces Core Collection (EUMLCC) of 96 populations, which synthesized the genetic variability present in the original set of genotypes [1,7]. Evaluations of the EUMLCC by partners explored grain composition through near-infrared (NIR) spectroscopy, revealing large variability especially of carotenoid compounds [8], as well as resistance and tolerance to corn borer attack [9] and tolerance to cold conditions.
In this paper, we present the concept applying PPP to the multilocation evaluation of European maize genetic resources, including part of the EUMLCC. We describe the contributing maize collections and criteria for selection of interesting landraces and present results of diversity studies based on genotypic data and phenotypic traits. We show how the EVA approach provides a model for effective collective generation of knowledge on PGRFA, making them more accessible for use in research and breeding.

2. Materials and Methods

2.1. Plant Material

A total of 626 different European maize landrace populations provided by nine genebank partners in the project were selected for inclusion in the EVA collection. Criteria considered for selection are described in detail in the results. Accessions were regenerated by holding institutes following an established protocol [10] and distributed with SMTA under the terms of the Multilateral System of the International Treaty for Plant Genetic Resources for Food and Agriculture (ITPGRFA) to partners for evaluations.
A set of five common registered check hybrids from company partners representing FAO maturity groups (210-300-400-500-600) was included in evaluations to enable cross-comparison of the different trials. Table 1 summarizes the EVA Maize collection and Table S1 lists the passport data of all accessions.

2.2. Genotyping

For genotyping of the 626 landraces in the EVA maize collection, each landrace population was represented through randomly sampling 15 individual plants from which equal amounts of leaves were mixed together prior to DNA extraction, to obtain one DNA pool per landrace as described previously [26]. Plants were germinated directly from genebank seeds, with no prior multiplication.
The DNA pools of the 626 landraces were genotyped using the 50K Illumina Infinium HD array [27] according to manufacturer’s instructions (Illumina Inc., San Diego, CA, USA). SNPs were filtered based on their suitability for diversity analysis and their quality for predicting allelic frequency in DNA pools, following a procedure described previously [26], resulting in a final set of 23,412 SNPs for further analysis.
Fluorescence data for alleles A and B of the 23,412 SNPs were extracted for each landrace with GenomeStudio’s Genotyping Module software (v2010.2, Illumina Inc., San Diego, CA, USA) and used for predicting allelic frequency using a freely available script (https://doi.org/10.15454/GANJ7J, accessed on 14 April 2024), following the two-step procedure described in [26] based on the fluorescence intensity ratio (FIR) of alleles A and B for each SNP. First, we tested whether SNPs were monomorphic or polymorphic for each landrace. For SNPs that were polymorphic within each landrace, we then estimated the allelic frequency of the B allele using a generalized linear model calibrated on FIR data from 1000 SNPs from two series of controlled pools (see [26] for more detail and Equation (2) for the model). This two-step approach led to a global mean absolute error of 3% and was more conservative for SNPs that were fixed or close to fixation than for SNPs with balanced allelic frequency [26]. The threshold to reject the hypothesis that landraces were monomorphic was set to 5%, indicating that 5% of landraces were expected to be declared polymorphic while actually being monomorphic (false positive).

2.3. Diversity Analysis

2.3.1. Estimation of Genetic Diversity Parameters

For each landrace, we averaged across 23,412 SNPs the mean allele number (A), the minor allele frequency (MAF), and the expected heterozygosity (H) [28,29]. Considering that all landraces were represented by 15 different plants (30 gametes), we did not apply a correction for the number of individuals when estimating these parameters because it would have resulted in only a small increase for diversity parameters (3.4% according to [30]). We also estimated genetic differentiation (Gst) between individual landraces using the 23,412 SNPs, according to [28]. We then extended the diversity measures at the complete panel level and for specific genetic groups (described below), using locus average allelic frequency. For genetic groups, only highly assigned landraces (>0.6) were considered.

2.3.2. Genetic Structure and Relationship between Landraces

We estimated the genetic distance between all landraces using modified Roger’s distance (MRD) [31] with ad hoc scripts, and fixation index (FST) based on the allelic frequencies of 23,412 SNPs using R package HFst [32]. To decipher the structure of genetic diversity within our panel of 626 landraces and their relations with worldwide diversity previously described in [33], we used three approaches:
(1)
A distance-based approach using MRDs between the 626 landraces to perform principal coordinate analysis (PCoA) [34], as well as hierarchical clustering. We applied either Ward or neighbor-joining algorithms using the “hc” and “bionj” functions of “ape” R package v 5.0 [35], respectively;
(2)
A Bayesian multi-locus approach, implemented in the ADMIXTURE v1.3.0 software, to assign probabilistically each landrace to K ancestral populations assumed to be in Hardy–Weinberg equilibrium [36]. Different methods were used to identify the most appropriate number of ancestral populations (K): cross-validation error or difference between successive cross-validations [36] and graphical methods [37]. Since ADMIXTURE requires multi-locus genotypes of individual plants, we simulated the genotype of five individuals for each population for a subset of 2500 independent SNPs to avoid artifacts of linkage disequilibrium (See Method S5 in [33] for more details).
(3)
A linear penalized regression approach, to quantitatively assign each of the 626 landraces to seven genetic groups established for 156 landraces representing European and American diversity [33]. Allelic frequencies at 23,412 SNPs of each of the seven genetic groups were estimated considering landraces with a membership superior to 0.6. Admixture coefficients associated with the new landraces were then obtained through a penalized regression approach via fitting, for each of the 626 landraces, a linear regression model using landrace allelic frequencies as response variables and group frequencies as explanatory variables, using the R package “quadprog” [38] v1.5-8. The coefficients of each regression were constrained to be positive and sum to one and were considered as estimates of the admixture coefficients.
We used R package “ggtree” [39] to draw the dendrograms derived from hierarchical clustering. Barplots from admixture analysis and penalized regression as well as the color scales for Hs, DS, and PH were added using “ggtreeExtra” R package [40,41].
Finally, to investigate the relationship between admixture and geographical origins, we mapped individual pie plots representing the admixture coefficients of accessions to their collection sites, using “ggplot2” [42] and “scatterpie” [43] R packages.

2.4. Phenotypic Evaluations

Among the 626 landraces that were genotyped, 38 could not be phenotyped due to limited seed availability. The 588 other landraces were evaluated by project partners in eleven locations in three sets of 175–218 accessions per set, over three years. Each set was evaluated on at least two locations in one growing season, yielding a total of 22 environments (Table 2). Each accession set covered maturity ratings from FAO100 to 800 and was distributed across locations to match local capacities and growing conditions, resulting in 2–9 trial locations per accession (Table S2). Evaluations followed a standard experimental protocol, using standard agronomic practices. The protocol consisted of evaluating a set of landraces along with a set of check hybrids, but the number of repetitions and the number of check hybrids varied in the different trials.
Data were collected for phenology traits, including days to tasseling [DT] and days to silking [DS] and for plant architecture traits, including plant height [PH] and ear height [EH], following standard protocols and using standardized scoring scales (Table 3) and were stored in a project specific database. Derived traits including anthesis–silking interval (ASI = DS-DT) and relative ear height (EPHR) were calculated from the raw data.

2.5. Statistical Analysis of Phenotypic Data

A combined analysis of variance was performed for the phenotypic traits studied using the PROC MIXED procedure in SAS software v. 9.4 (SAS Institute Inc.: Cary, NC, USA). The best linear unbiased estimators (BLUEs) for each variety were calculated based on pooled data for all environments. Varieties were considered as fixed effects while replication (trials) was considered as random.
A principal component analysis was carried out on the BLUEs of the 588 landraces and 5 check hybrids using the six evaluated traits (Table 3). The data were standardized (mean = 0 and standard deviation = 1) to avoid scaling problems. Principal components with an eigenvalue > 1 became new variables, with which cluster analysis was carried out. For cluster analysis, we first used the FASTCLUS procedure of SAS to perform preliminary clustering using Euclidean distance (the recommended method for more than 100 varieties), choosing 15 preliminary clusters based on the pseudo F statistic. Subsequently, we performed the CLUSTER procedure using average linkage as a clustering method. Finally, the results of the hierarchical clustering analysis were displayed as a dendrogram using the TREE procedure (SAS Institute Inc.: Cary, NC, USA).

3. Results

3.1. Creating the EVA Maize Collection

The EVA maize collection included materials provided by genebank partners from Croatia, France, Italy, Portugal, Romania, Serbia, Spain, and Switzerland, reaching a total of 626 accessions, originating from 26 countries, and focusing on locally adapted landraces of various FAO maturity groups (Table S1). These constituted key accessions for each holding institute (see Supplementary Materials S1 for details about each genebank). Four genebanks previously involved in the RESGEN88 project contributed 75 accessions belonging to the EUMLCC core collection. The selection processes adopted by each genebank to include accessions into the EVA maize collection are described below and summarized in Table 1.
Croatia—University of Zagreb, Faculty of Agriculture (HRV041): The subset of 50 landraces from the collection of the University of Zagreb Faculty of Agriculture (Croatian genebank) included in the EVA maize collection was selected to represent the collection based on its wide geographical distribution as well as on previous characterization data [11]. It included accessions from the continental and Mediterranean regions of Croatia with a wide range of variation in phenological and morphological traits such as maturity group, plant height, ear length, and ear width.
France—INRAE—Montpellier (FRA015): The CRB GAMéT genebank made available for the EVA maize collection the French representative collection (RMNC, maintained with the Promaïs partnership), initially defined for the RESGEN88 project. A set of 80 landraces (16 of them inserted in EUMLCC) were chosen based on morphological traits using Mstrat software (v 4.1) to maximize the number of phenotypic classes with at least one population. The morphological traits considered for selection were measured and analyzed in the classification of French maize populations [7,12].
Italy—CREA—Cereal and Industrial Crops, Bergamo (ITA386): The set of 65 landraces shared by CREA for the EVA maize collection included 19 varieties of the Italian collection belonging to the EUMLCC selected in the frame of RESGEN88 [8]. Additionally, accessions originally collected around 1950 from different Italian areas were chosen. The Italian landraces selected for the EVA maize collection have also recently been included in different programs devoted to valorization of local genetic resources and their possible utilization as sources of stress-tolerance traits in pre-breeding activities [13].
Portugal—INIAV, Braga (PRT001): The set of 42 Portuguese landraces shared by INIAV for the EVA project included 17 varieties of the Portuguese collection belonging to the EUMLCC. Twenty-five additional accessions were chosen from within the frame of the national collection of 70 varieties based on morphological, agronomic, and molecular characterization, as well as on previous knowledge from national breeding programs [14,15].
Portugal—ESAC-IPC, Coimbra (PRT053): The six Portuguese landraces shared by the Escola Superior Agrária de Coimbra (ESAC, Portugal) in the EVA maize collection included accessions collected during a collecting mission focused on maize varieties for traditional bread production, and others were tested in drought trials [16]. Additionally, some accessions were selected based on previous studies, in which yield evaluations were conducted at different altitudes. The ‘Bilhó’ OPV population (white flint maize) was selected at higher-altitude eco-geographic areas, generally with a lack of water, in contrast with ‘Sangalhos’ from low altitudes.
Romania—Suceava Genebank (ROM007): The maize collection in the Suceava Genebank represents 32% of the entire Romanian national collection. Most of the 6143 maize accessions in the national collection are local landraces from Romania (60%), especially from mountainous and sub-mountainous areas (up to 400 m altitude). The 51 maize landraces included in the EVA maize collection were selected from among early and semi-early accessions (FAO groups 200–300), collected between 20 and 50 years ago (1973–2000), with resistance to low temperatures [17,18].
Serbia—Maize Research Institute Zemun Polje (SRB001): The criteria followed for selection of 91 accessions by MRIZP, Serbia, included drought tolerance observed in a two-year screening program under a managed stress environment in Egypt, stability and high grain yield in a two-year trial under rain-fed conditions and with higher plant density applied under temperate conditions at Zemum Polje, and good performance regarding traits important for breeding [19,20], representing valuable donors of genes responsible for higher grain yield in crossings with germplasm from different heterotic groups. The selection included mainly accessions from the Western Balkans but also some exotic materials that may be a valuable source for drought tolerance.
Spain—Misión Biológica de Galicia, Pontevedra (MBG-CSIC)—(ESP004, ESP007, ESP009, ESP016, ESP019): Among the 137 Spanish landraces included in the EVA maize collection, 23 landraces were also contained within the EUMLCC. As Spanish germplasm is very variable, all Spanish landraces with enough seed availability were included. Previous knowledge on corn borer resistance, drought and cold tolerance or suitability for making bread was available for many of them [1,9,21]. Pontevedra (ESP009) coordinated the distribution of the accessions from other Spanish holding institutes.
Switzerland—Agroscope, Nyon (CHE001): The 86 accessions included in the EVA maize collection are Swiss and Liechtenstein maize landraces selected to represent different collecting phases as well as different regions of Switzerland. The accessions have phenological and morphological data such as plant height, ear length, kernel color, and others [22,23,24,25].
The EVA maize collection of 626 European landraces was distributed to EVA Maize Network partners for field evaluations and genotyping, with the aim to assess the genetic diversity available in the collection and to evaluate the landraces in multilocation trials for phenotypic characteristics important for further use in breeding and research.

3.2. Genotypic Diversity of the Collection

Out of 23,412 useful SNPs, 4 were monomorphic for the whole EVA panel. Minor allelic frequency (MAF) was 0.25 (± 0.14) (Table 4). At group level the average number of alleles per locus was 2.00 (± 0.01), while within landraces it was 1.71 (± 0.13) and ranged from 1.23 (SRB0772, Zuti Poluzuban) to 1.92 (HRV0770, Domaci Osmak) (see Table S3). The mean total expected heterozygosity among (Ht) and within landraces (Hs) in the whole panel were 0.33 (± 0.15) and 0.21 (± 0.05), respectively (Table 4), leading to a mean genetic differentiation between landraces (Gst) of 0.35 (± 0.14) across the whole panel. The expected heterozygosity within landraces varied from 0.03 (SRB0772, Zuti poluzuban) to 0.30 (ITA0083, Giallo Tosoratti) (Table S3, Figure S1). Only 13 landraces had an Hs below 0.10, suggesting that these landraces were either not well maintained in the genebank or that they had been under strong bottleneck during their initial collection.
We analyzed the genetic structure of 626 landraces using the ADMIXTURE program [36]. Likelihood analysis of ADMIXTURE results indicated that the optimal numbers of genetic groups were K = 2, K = 4, K = 6, and K = 9 (Figure S2). We compared this new genetic structuration for successive K values with the assignment to 7 genetic groups determined in [33], using a penalized regression approach (Table S3, Figure S3). At K = 2, “Northern Flint” landraces split from other landraces, which is highly consistent with the first axis of the PCoA analysis (Figure 1) and the first split in the dendrogram (Figure 2). At K = 4, genetic groups corresponded to those identified in Europe [33]: Northern Flint, Pyrenean–Galician Flint, Italian Flint, and Corn Belt Dent. At K = 6, the groups of Northern Flint and Italian Flint were conserved, while the Pyrenean–Galician group was split into two equal genetic groups corresponding to collections from humid Northern Spain (landraces from Galicia) and French/Spanish landraces from the Pyrenees. Finally, a small group of popcorn landraces from Spain, Romania, and France was separated from the group composed of Corn Belt Dent landraces. At K = 9, the groups Northern Flint (K8), Corn Belt Dent (K4), popcorn (K1), and Pyrenean (K9) were maintained. On the other hand, Portuguese landraces, which were considered an admixture between Pyrenean–Galician, Andean, and Northern Flint Groups [45], became a new genetic group (K6), closely related to the genetic group of Galician landraces (K7). Similarly, the Italian Flint group [45] was split into two equal genetic groups corresponding to landraces from Northern Italy (K5) and South–Central Italy (K3). Interestingly, landraces collected along the Adriatic coast of Croatia and Montenegro were strongly assigned to the South–Central Italy (K3) genetic group. In addition, Southern dry Spanish landraces, which penalized regression showed to be an admixture between Caribbean and Italian Flint groups [45], became a specific genetic group (K2).
The new genetic structuration with K = 9 was highly consistent with our previous analysis based on 156 worldwide landraces [33] and with the history of both introductions of maize in Europe, the first from the Caribbean into Southern Spain and the second from North America into Northern Europe. In addition, it allowed us to go deeper in the understanding of spatial organization of the genetic diversity along the European continent. For further analysis, we focused on this last structure of ancestry.
The size of these nine genetic groups varied from 14 (K1, K6) to 75 (K8) landraces (Table 4). We compared the diversity parameters of landraces assigned to nine genetic groups, considering only landraces with a membership above 0.6 (Table S3, Table 4). The Ht in the nine groups varied from 0.22 (K8, Northern Flint) to 0.33 (K4, Corn Belt Dent), with the latter value equivalent to the whole panel. The mean within-gene diversity in each group (Hs) varied from 0.14 (K8, Northern Flint) to 0.24 (K4, Corn Belt Dent). Average genetic distance (MRD) between landraces in each group ranged from 0.16 (K6, Portuguese and K7, Galician) to 0.23 (K2, Southern Spain). The number of monomorphic markers in each group varied from 956 in group K1 (popcorn) to 81 in group K4 (Corn Belt Dent). Genetic differentiation (Gst) within genetic groups varied from 0.17 (K6, Portuguese) to 0.36 (K2, Southern Spain and K8, Northern Flint). The estimation of fixation index (FST) between genetic groups added information on their relationship (Table S4). The largest FST (0.32) was found between the group of K2 (Southern Spain) and K8 (Northern Flint) while the least differentiated genetic groups were K7 (Galician) and K9 (Pyrenean)—0.05. In addition, the genetic group K6 (Portuguese) was also weakly differentiated from these two groups. Little differentiation was found between landraces from groups K5 (Northern Italian) and K3 (South–Central Italian; 0.08) or between them and group K2 (Southern Spain; 0.10 and 0.09, respectively). Finally, K4 (Corn Belt Dent) and K6 (Portuguese) groups were weakly differentiated (0.10).
To better understand the relationships between landraces and the nine genetic groups, we performed PCoA analysis based on MRD distance (Figure 1) and clustering analysis using both the Ward algorithm based on MRD distance (Figure 2) and neighbor joining based on FST (Figure S4).
In Figure 1, the first axis (PC1, 13.9% of total variation) discriminated landraces from Northern Europe (Switzerland and some landraces from Northeastern France) corresponding to Northern Flint landraces (K8) against landraces from Southern Spain (K2). This was consistent with the two groups identified with the ADMIXTURE analysis at K = 2 and corresponded to the first split in the dendrogram tree (Figure 2). The second axis (PC2, 7.27% of total variation) separated French, Spanish, and Portuguese landraces collected in the regions of the Pyrenees and Galicia (K6, K7, K9) from Corn Belt Dent and popcorn-like landraces from France, Romania, and Southern Spain (K1 and K4 genetic groups). Interestingly, flint landraces from Northern and South–Central Italy were intermediate between Southern Spanish landraces and Pyrenean–Galician–Portuguese Flint landrace groups (K6, K7, K9) on the first two axes of the PCoA. The third axis (PC3, 4.72% of total variation) separated Spanish and Galician Flint landraces from Northern and South–Central Italian Flint landraces (K3 and K5). The fourth PCA axis (PC4, 1.9% of total variation) separated popcorn landraces (K1) from the rest of the panel (Figure S5).
Dendrogram trees based on Ward clustering (Figure 2) or neighbor joining (Figure S4) provided a finer view of the relationship between the clustering and ancestry of the landraces. At a higher level of the dendrogram tree, landraces from Switzerland and Northern France (K8) assigned to the Northern Flint Arca group [33] were split from all other landraces/genetic groups, following the PCoA first axis and ADMIXTURE groups at K = 2. Accordingly, the Northern Flint group displayed the highest differentiation from landraces from five other genetic groups (K1, K2, K3, K4, K5 with FST ranging from 0.19 to 0.32 in Table S4). Northern Flint landraces were closer to K6 (Portuguese, FST = 0.14), K7 (Galician, FST = 0.21), and K9 (Pyrenean, FST = 0.17). This was expected since the Pyrenean–Galician varieties resulted from admixture events between Northern Flint from Northern Europe and Southern Spanish landraces from Caribbean origins [46]. The landraces were then divided into two distinct groups in the tree: (i) landraces from the Pyrenees, Galicia, and Portugal (K6, K7, K9) corresponding to the Pyrenean–Galician Arca group [33], which were the least differentiated groups (FST < 0.07 in Table S4), and (ii) landraces from Italy, Spain, and Southeastern Europe (K2, K3, K4) corresponding to Italian Flint and Corn Belt Dent Arca groups, respectively. At a lower level, the Pyrenean–Galician–Portuguese landraces split in different subgroups: (i) a first cluster of pure French–Pyrenean landraces (K9) with admixed French landraces with Northern Flint or Corn Belt Dent and admixed Spanish–Galician landraces; and (ii) a second cluster of Portuguese (K6) and Galician landraces (K7) with some admixed landraces between these two groups, originating from Northern Portugal. The Pyrenean–Galician–Portugal cluster was therefore mainly composed of three genetic groups: K6 (mainly Portuguese landraces), K7 (Spanish–Galician landraces), and K9 (mainly landraces collected in the French Pyrenees). Interestingly, some subgroups were admixed almost exclusively between these or with K4 (Corn Belt Dent) and K8 (Northern Flint).
On the other branches of the tree, the landraces were divided into two major groups: (i) landraces highly or partially assigned to the Italian Flint Arca group (or admixed with it), and (ii) landraces mainly from Southeastern Europe assigned in great part to the Corn Belt Dent Arca group. The Italian Flint Arca group was split in three subclusters: (i) one subcluster of landraces from Northern Italy and Southern Switzerland, that were purely assigned to Italian Flint (K5); (ii) a new subcluster of landraces from South–Central Italy and from Croatia along the Adriatic coast (K3), and (iii) a third subcluster of landraces from Southern Spain, cultivated in dry areas (K2), that appeared as admixed between Italian Flint and the Caribbean Arca group. A small subcluster presented a level of admixture with the K1 genetic group and was composed of popcorn varieties. The second subcluster, partially assigned to Italian Flint, included highly admixed landraces, also not clearly assigned to any ancestral population. It was mainly composed of landraces collected in Spain, Italy, France, and Croatia.
Finally, the cluster highly assigned to Corn Belt Dent was made up of three clearly separated subclusters. The first one was a small group including popcorn varieties mainly from Spain and identified as a separate genetic group through the ADMIXTURE analysis (K1). The second one was composed of accessions highly assigned to the Corn Belt Dent ancestral group and K4 genetic group. This group included landraces mainly collected in Croatia, Serbia, and Bosnia. The third subcluster was composed of admixed landraces mainly from Romania, Croatia, and Bulgaria. Landraces in this group were partially assigned to a mixture of most ancestral groups, except for Caribbean. In terms of ADMIXTURE results, most landraces showed different levels of admixture between K4, K8, and K9.
The organization of the neighbor-joining dendrogram based on FST was almost identical, except that the Corn Belt Dent (K4) cluster was closer to the Northern Flint genetic group (K8) (Figure S4). In the same direction, a set of subclusters of mainly Romanian landraces were placed between Northern Flint and Corn Belt Dent landraces following a gradual decrease in Northern Flint ancestry. Specifically, the subcluster integrated with admixed Pyrenean–Northern Flint landraces, placed close to Pyrenean in the Ward dendrogram, was moved next to the Northern Flint landraces.
To analyze the spatial structure of these genetic groups, we mapped the admixture coefficients of the nine genetic groups for 483 landraces for which GPS coordinates of collection sites were available (Figure 3). Landraces belonging to the nine genetic groups clustered according to their geographic origins except for popcorn landraces that spread from Southern Spain to Southwestern France. The genetic structure matched very well with eco-geographical regions corresponding to various environments encountered in Europe, suggesting that landraces of these different genetic groups corresponded to adaptation to specific environments. For instance, Corn Belt Dent landraces (K4) were mapped mainly in continental climate areas in Southeastern Europe with hot summers, while Northern Flint landraces (K8) were cultivated in the North of Europe and well adapted to cold temperatures during spring sowings. Landraces from Spain were split into two groups corresponding to two contrasted environments: humid oceanic climate in Galicia vs. hot and dry Mediterranean climate in Southern Spain. Contrasting environments could also explain the shift into different genetic groups between Pyrenean and Galician (transition from the Mediterranean to Atlantic/Lusitanian region), as well as between Galician and Portuguese landraces (transition from the Atlantic to Mediterranean region). Landraces from Southern Spain appeared as the closest relatives to landraces from the Canary Islands that were assigned to the Caribbean Arca group (Figure S4). In addition, cultural or historical processes may have contributed to the geographic distribution of genetic groups. This is likely to have been the case for a group of Swiss landraces highly assigned to the K5 group that were collected in the Ticino region, known for being in close relationship with Italian culture. The geographic structure was also evident when looking at different K values detected in the ADMIXTURE analysis (Figure S6). Interestingly, at K = 4 and K = 6, landraces from the regions of Southern Spain and Southern Italy mainly belonged to the same genetic groups (K1 and K3, respectively).

3.3. Phenotypic Diversity of the Collection

Based on seed availability and field capacity for evaluations, 588 out of the 626 landraces were evaluated in field trials. Five commercial hybrids of different FAO cycles (from 210 to 600) were included in the trials as checks to find out which landraces were best suited to specific environments. To assess the phenotypic diversity of the collection, we analyzed phenological traits including days to silking (DS), days to tasseling/anthesis (DT), and anthesis–silking interval (ASI) as well as plant architecture traits including plant height (PH), ear height (EH), and relative ear height (EPHR); mean data were based on data recorded in 2–9 locations per accession (Table S2). These traits were chosen for preliminary phenotypic description of the European landraces because they provide key information for adaptation to different environments. BLUEs for each of the traits studied and for all landraces and checks can be found in Table S2.
The great amplitude of range of variation showed that there was large variability among landraces for all traits studied (Table 5). Regarding male and female flowering, about 25% of the landraces were earlier than the earliest FAO210 hybrid check. Also, there were landraces that flowered later than the latest hybrid FAO600. In terms of plant and ear height, 75% and 50% of the landraces, respectively, were shorter than any check hybrid but some landraces were also taller than the checks. In terms of ASI, the wide range of variation of the landraces contrasted with that of the hybrids as did the ratio between plant and ear heights.
The wide range of variation in the EVA collection for female flowering and plant height is shown in Figure 4 and Figure S7. Early landraces were assigned in general to genetic groups K7 and K9, with K7 landraces being shorter. As for late flowering landraces, they were assigned mainly to K1, K2, and K4, landraces from K4 being taller than the others. Admixed landraces were also found among the extreme phenotypes for flowering or plant height. In general, all genetic groups showed large variability in both flowering time and plant height. Interestingly, landraces assigned to K3 (South–Central Italian) were shorter and earlier than landraces from K5 (Northern Italian), which could be associated with specific uses or management practices. In addition, ASI was similar for all the genetic groups except for K8 (Northern Flint), which showed the highest ASI.

3.4. Principal Component Analysis

The highest correlations among the traits studied were between flowering and plant (0.97) and ear height (0.89) while the lowest correlations were observed between ASI and the other traits, ranging from 0.17 for EPHR to 0.47 for DS.
Principal component analysis was carried out to summarize the information relating to the six traits evaluated. Only the first principal component (PRIN1) had an eigenvalue greater than 1, and this explained 72% of the variability. The second principal component (PRIN2), with an eigenvalue close to 1, explained 16% of the total variance. In total, 95% of the total variance was explained through the first three principal components (Table 6).
Days to flowering and plant architecture were the characteristics with the highest (positive) influence on PRIN1, while ASI presented the highest (positive) coefficient for the second principal component. This means that tall and late varieties showed high PRIN1 values while high PRIN2 values corresponded to high ASI (Table 6).
In the graphic representation of the PCA, landraces were positioned in all four quadrants, corresponding to substantial phenotypic variation in all genetic groups (Figure 5). The hybrid checks were mainly located in quadrant IV (positive PRIN1 and negative PRIN2 values) as they were tall and late flowering and presented small anthesis–silking intervals. For Axis 1, landraces assigned to K3 and K7 were placed mainly in the quadrants II and III, while landraces assigned to K1, K2, K4, and K5 were positioned in the quadrants I and IV. In reference to Axis 2, most varieties from Switzerland assigned to K1 were positioned in the two upper parts of the PCA chart (quadrants I and II).

3.5. Cluster Analysis

Cluster analysis based on the phenotypic data of 588 landraces and five hybrid checks grouped them into 15 subclusters, which, in turn, were grouped into five main clusters A–E, with clusters B and E containing most landraces (Figure 6, Table 7 and Table S2).
Cluster A was the smallest cluster, formed of only subcluster 1, containing two landraces characterized as very late and tall. These Spanish popcorn populations (ESP0200, Sangonera La Seca and ESP0115, Alcantarilla) were embedded in the K2 genetic group, but were considered admixed, as they shared ancestry with other genetic groups, including K1 (popcorn).
Cluster B was formed of 308 landraces with trait values slightly higher than the average landrace. Cluster B was formed of six subclusters (2, 3, 5, 7, 11, and 14). The closest subclusters 3 and 11 grouped 86 and 66 landraces, respectively, and contained the earliest flowering landraces of cluster B. In subcluster 3, more than half of the landraces were admixed (<0.6 probability of assignment to a particular genetic group). More than half of these admixed landraces had K4 as their nearest group (Corn Belt Dent germplasm), and ten landraces from Eastern Europe were also assigned to K4. Likewise, almost two thirds of the landraces in subcluster 11 were admixed, although many of them were related to K4. The remaining landraces in this subcluster were distributed across genetic groups. None of the landraces belonged to the ‘Caribbean’ group [33], while few of them were included in the Northern Flint group described in the same study. Subclusters 7 and 5, with 45 and 48 landraces, respectively, grouped the latest flowering landraces of cluster B. In both subclusters, more than half of the landraces belonged or were related to the K4 (Corn Belt Dent germplasm) genetic group. Subcluster 14 consisted of 19 intermediate flowering landraces. Genetically, most of the landraces were included in or related to groups K1 and K2 and came from Southeastern Spain, including popcorn landraces. Subcluster 2 consisted of 44 landraces, most of them from Switzerland, characterized by large ASI and ascribed, in a high proportion, to the K8 genetic group (Northern Flint germplasm). In summary, more than 60% of the populations belonging to groups K1, K2, K4 and K5 were included in phenotypic cluster B, while less than 15% of the populations belonging to group K3 were in this cluster. So, this phenotypic cluster encompassed Caribbean, Northern Italian, and Corn Belt Dent germplasms.
Cluster C consisted of 44 landraces; 29 belonging to subcluster 8 and 15 to subcluster 4. Cluster C was a group of late landraces with high insertion of the main ear. In subgroup 4, 14 out of the 15 landraces originated from Eastern or Southern Spain or the Canary Islands. Subcluster 4 included some popcorn landraces and was highly related to genetic groups K1 and K2. Subcluster 8 was formed of landraces from different origins that were highly related to the genetic group K4 (Corn Belt Dent).
Cluster D, representing only subcluster 6, was made up of eight early landraces, most of them being part of the K9 (Pyrenean) or K7 (Galician) genetic groups. They stood out, on average, for shedding silks before anthers and, consequently, presenting negative ASI values.
Cluster E was made up of 227 early and short varieties (226 landraces and the earliest hybrid check F210) grouped into five subclusters: 13 (32 landraces), 15 (25 landraces), 12 (46 landraces), 10 (76 landraces), and 9 (47 landraces). Out of these 226 landraces, 114 were admixed, 15 belonged to K3 and 5 to K5 genetic groups (Italian landraces), 29 to K7 (Galician genetic group), 19 to K9 (Pyrenean genetic group), and 38 (mainly in sub-cluster 13) to K8 (Northern Flint group). Sub-clusters 12 and 10 were joined first, followed by sub-clusters 9, 15, and 13.
Clusters B and C joined to form a new group of late and tall landraces with large anthesis–silking intervals. Genetically, they belonged or were mostly related to the K4 genetic group, but also to the K1, K2, and K5 groups. They represented the European landraces with the greatest contribution from the Corn Belt Dent and Caribbean germplasms. Clusters D and E united to form a new group of early and short landraces with a high presence of landraces from the Northern Flint (K8), Galician (K7), and Pyrenean (K9) genetic groups. Clusters B–C and D–E finally joined cluster A, which was a small cluster formed by two popcorn landraces that stood out for being late and presenting remarkably high values for ear height and relative ear height.

4. Discussion

The EVA Maize Network enabled for the first time the systematic evaluation of a large collection of European maize germplasm in multilocation trials, which, together with the SNP genotyping of all accessions, provided an overview of the diversity of maize available in European genebanks. Previous studies on the European maize core collection (EUMLCC) only included accessions from six countries and thus captured only a limited range of European germplasm. In addition, characterization of the EUMLCC was carried out for specific traits in specific environments and was therefore not as extensive as that carried out in the EVA network.
The EVA maize collection presented in this study captured a large proportion of European maize diversity. Some landraces had been analyzed in previous studies using SNP [33], RFLP [6,45,47], SSR [48,49], or isozymes [1]. Seventy-five accessions in the EVA collection were also included in the EUMLCC and evaluated for various phenotypes in the RESGEN88 project. Despite not including German and Greek landraces from the original EUMLCC panel, the EVA collection extended the diversity to populations from Switzerland, Croatia, Romania, and Serbia, and included additional materials from France, Italy, Spain, and Portugal. This new study reveals genetic diversity that had not previously been represented in the EUMLCC and could be used to review the composition of the core collection. For example, we identified new, unique clusters in the genotypic grouping, which did not contain any accessions previously included in the EUMLCC. Genotyping data highlighted a clear clustering of accessions into nine genetic groups. The EVA collection contained 349 accessions that could be assigned to a prevalent group and 277 admixed accessions not assigned to any genetic group (membership <0.6 to any genetic group). The mixed genetic background of the latter is likely to have resulted from breeding activities since the introduction of maize to Europe [50]. The main lines of our results appear consistent with historical records and findings from previous molecular studies. Gauthier et al. [6] classified maize populations from Europe into three main clusters consistent with geographic origins: a Northeastern cluster, a Southeastern cluster with semi-dent populations from Southeastern Europe, and finally a Northwestern cluster with early flint populations from Northwestern Europe. They also detected a possible direction of gene flow indicating that Northern Flint and Caribbean populations were introduced to Northern and Southern Europe, respectively [51,52,53]. The large representation of diverse germplasm from several countries in the present study highlighted finer-scale patterns.
The K1 molecular group, highlighted in the present research, is a special type of Mediterranean germplasm, composed mainly of popcorn accessions. These were probably genetically isolated from the rest of the Mediterranean populations due to the presence of the gametophytic factor in most popcorn populations [54]. Due to the long growth cycle of these materials, some of them were included in phenotypic cluster B because of their large plant size, while other accessions were in cluster C because of their late cycle and high ear insertion. Following this clustering, two populations from the FRA015 collection were re-investigated and validated as popcorn (M. Vincent, personal communication), showing how the present study can contribute to genetic resources characterization.
Group K2 could be explained through the first introductions of maize from the Caribbean area into Seville in 1494 [1] and subsequent introductions throughout the South of Spain. Though these first introductions were probably not very successful, their singularity is still present in that area. These varieties were included in phenotypic cluster A because they are late and tall, as would be expected from the photoperiod sensitivity related to their Caribbean origin [55]. This may also explain why some of these populations were included in phenotypic clusters B and C, with large, late flowering plants. Similarly, groups K3 (South-Central Italy) and K5 (Northern Italy) represent the diverse introductions of maize made throughout the Mediterranean area. Maize was introduced in Spain and Italy almost simultaneously from Central American sources, giving rise to singular types of germplasm that remained in both countries. Some of the ancient introductions also had large plants because of their tropical origins [56,57]. The proximity of populations from Italy to populations from Spain has been shown to represent one of five genetic groups in 67 European and American inbred lines [46]. Interestingly, populations from Southern Italy, which are genetically close to Spanish populations from Southeastern Spain (K2), were significantly earlier than those populations and also those of Northern Italy. This “switch” in phenological traits may be explained by the fact that farmers in Southern Italy selected for materials with shorter cycles to avoid drought and/or high summer temperatures. Additionally, short-cycle maize landraces provided optimal results in terms of crop rotation for second cropping after wheat cultivation typical in these areas [58]. Finally, a small group of landraces from the Adriatic coast of Croatia were strongly assigned to this genetic group. These Croatian landraces were also phenotypically distinct from the remaining landraces collected in the same region, matching the phenotypic characteristics of the K3 group. This may suggest direct introduction from Italy and isolation, possibly related to the same cultivation practices applied on the other side of the Adriatic Sea.
Northern Flints (K8) were introduced during the early 16th century [2]. The genetic groups K6, K7, and K9 are likely to have derived from admixture between these Northern Flints and tropical introductions. This is consistent with the origin of populations in Southwest France [59]. Some of these accessions were included in the phenotypic cluster D because they were probably selected for adaptation to short growth cycles. However, a larger number of them were included in cluster E, representing the introduction of flint from North America along the Atlantic coast.
Finally, in the late 19th and 20th century with the release of improved varieties, maize was introduced continuously from America, mainly from the USA. Therefore, the presence of Corn Belt Dent (K4) was established later than other groups, mainly in Eastern European countries where late varieties with high yield were very successful. Accordingly, many populations from Eastern Europe belonged to the Corn Belt Dent genetic group (K4) and to phenotypic clusters B and C and showed larger plant development than accessions from other clusters.
The analysis of phenotypic data for flowering time and plant or ear height explained only partially the separation patterns of the genotypic clustering. This indicates that the genetic groups exhibit large phenotypic diversity for these traits. The genotypic study shows the phylogenetic relationships due to the origins of the populations, while the similarities between varieties in plant height or flowering dates reflect selection for adaptation to diverse European environments into which the genetic groups were introduced. Similar observations were also reported in other studies [60,61]. Other traits may be better correlated to the genetic groups, and analysis of other evaluations of the EVA collection will provide information on this hypothesis. Another possible explanation for the missing separation into groups is the use of BLUEs across multiple trial locations for the analysis of phenotypic diversity. Perhaps disaggregating data according to location would show clearer patterning in the phenotypic PCA correlating with genetic groups.

5. Conclusions

The EVA Maize collection described in this study highlights the large genotypic diversity available across large ecogeographic regions in Europe. This collection provides good material for breeding programs to tackle climate change, and some landraces have also been registered as conservation varieties, making them available for growers seeking unique characteristics for specialty food products [5]. Others are being used in participatory breeding programs, which adapt a holistic transdisciplinary approach aiming at valorization along the entire value chain and promoting cultural legacy and identity [62].
Although not all European regions are represented in this collection, the goal of the ECPGR EVA initiative is to further explore genebank holdings through extension of the EVA Maize Network to countries in Central and Eastern Europe. Additional evaluations for agronomic traits and stress tolerance will add further knowledge of the collection, enabling identification of new variations to be used in the development of pre-breeding materials [63,64,65,66,67]. Furthermore, connecting the results with other ongoing (Dromamed, MineLandDiv) and future projects will ensure effective exploitation of the genetic resources conserved in European genebanks.
The EVA Network connects genebank collections with research institutes and public and private sector breeders. With its unique approach and much of the joint work provided as in-kind contributions by partners, these networks have been set up as long-term initiatives to continue exploring genetic diversity conserved in genebanks or on farms and make material and information available to end users, breeders, and growers.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/biology13060454/s1, Supplementary information 1: Description of Genebanks involved in the EVA Maize network. Figure S1: Genetic diversity (Hs) frequency for the complete panel and genetic groups from ADMIXTURE (K = 9). Figure S2: Determination of the K value for ADMIXTURE using deltaK. Figure S3: Ward dendrogram of genetic diversity based on Modified Rogers Distance of 626 landraces, for different K values. Figure S4: Neighbour joining dendrogram based on Fst between 626 landraces. Figure S5: Principal Coordinate Analysis based on Modified Roger’s Distance estimated with 23,412 SNP, representing (i) axis 1 and axis 3, (ii) axis 1 and axis 4. Figure S6: Map of origins of 483 EVA Maize accessions colored according to their belonging to genetic groups. Figure S7: Relationship for female flowering (days to silking, DS) and plant height (PH) among 588 maize landraces of the EVA collection. Table S1: Passport information of 626 landraces in the EVA Maize panel, including GPS location of collection sites and information about membership in the European Maize landraces core collection (EUMLCC) defined by the RES GEN 88 project. Table S2: Phenotypic data and statistical analysis for 588 maize accessions in the EVA panel and five check hybrids. Table S3: Genetic structure analyses of 626 EVA Maize landraces: diversity parameters, penalized regression to 7 ancestral groups of Arca et al. (2023), Admixture analysis (K = 9) and Principal Coordinate analysis (PCoA). Table S4: Fixation index (FST) within and across each genetic group (K = 9) of the EVA Maize panel.

Author Contributions

Conceptualization, C.B. (Carlotta Balconi) and S.G.; formal analysis, A.G., R.A.M., S.D.N. and T.M.-H.; investigation, V.A., C.B. (Carlotta Balconi), C.B. (Cyril Bauland), V.C., B.G., B.K., N.K., D.M. (Delphine Madur), P.M.-M., D.M. (Danela Murariu), A.P., D.D.P., R.R., P.R., H.Š., A.S. and A.T.; resources, V.A., A.M.B., A.B., N.K., R.A.M., P.M.-M., D.M. (Danela Murariu), H.Š., B.S.-V., A.T., M.V. and A.Z.; data curation, A.G, S.G., R.A.M. and S.D.N.; writing—original draft preparation, V.A., C.B. (Carlotta Balconi), A.M.B., A.B., A.G, B.G., S.G., R.A.M., P.M.-M., D.M. (Danela Murariu), S.D.N., R.R., H.Š., B.S.-V., A.T., M.V. and A.Z.; writing—review and editing, C.B. (Carlotta Balconi), A.C., V.C., A.G., B.G., S.G., R.A.M., S.D.N., P.R. and R.R.; supervision, C.B. (Carlotta Balconi) and S.G.; project administration, S.G.; funding acquisition, S.G. All authors have read and agreed to the published version of the manuscript.

Funding

The authors are grateful for the financial support for this research provided by the German Federal Ministry of Food and Agriculture, grant GenRes 2019-2 to ECPGR, which allowed the implementation of the EVA networks.

Institutional Review Board Statement

Not applicable.

Data Availability Statement

Data generated within the EVA project are subject to an embargo period under the EVA cooperation agreement and will thereafter be made publicly available according to FAIR principles. Reasonable requests for access to raw data can be addressed to the corresponding authors.

Acknowledgments

We want to thank all partners in the EVA Maize Network for the collaboration and their comments on the manuscript and Nora Capozio for help with designing the graphical abstract. The authors pay tribute to Jacques Dallard, coordinator of the European RES GEN 88 project, who died in 2023.

Conflicts of Interest

The authors declare that no conflicts of interest exist. Dr. Bettina Kessel is employed by the company KWS SAAT SE & Co, Germany; Dr. Alexandre Strigens is employed by the company DSP—Delley Semences et Plantes SA, Switzerland. All 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.

References

  1. Revilla, P.; Soengas, P.; Cartea, M.E.; Malvar, R.A.; Ordás, A. Isozyme variability among European maize populations and the introduction of maize on Europe. Maydica 2003, 48, 141–152. [Google Scholar]
  2. Rebourg, C.; Chastanet, M.; Gouesnard, B.; Welcker, C.; Dubreuil, P.; Charcosset, A. Maize introduction into Europe: The history reviewed in the light of molecular data. Theor. Appl. Genet. 2003, 106, 895–903. [Google Scholar] [CrossRef]
  3. ECPGR Plant Genetic Resources Strategy for Europe. European Cooperative Programme for Plant Genetic Resources, Rome, Italy. Available online: https://www.ecpgr.org/resources/ecpgr-publications/publication/plant-genetic-resources-strategy-for-europe-2021 (accessed on 14 April 2024).
  4. ECPGR: European Evaluation Network (EVA). Available online: https://www.ecpgr.cgiar.org/eva (accessed on 1 May 2024).
  5. Revilla, P.; Alves, M.L.; Andjelkovic, V.; Balconi, C.; Dinis, I.; Reis Mendes-Moreira, P.M.; Redaelli, R.; Ruiz De Galarreta, J.I.; Vaz Patto, M.C.; Slađana, Ž.; et al. Traditional foods from maize (Zea mays L.) in Europe. Front. Nutr. 2022, 8, 683399. [Google Scholar] [CrossRef]
  6. Gauthier, P.; Gouesnard, B.; Dallard, J.; Redaelli, R.; Rebourg, C.; Charcosset, A.; Boyat, A. RFLP diversity and relationships among European maize populations. Theor. Appl. Genet. 2002, 105, 91–99. [Google Scholar] [CrossRef]
  7. Gouesnard, B.; Dallard, J.; Bertin, P.; Boyat, A.; Charcosset, A. European maize landraces: Genetic diversity, core collection definition and methodology of use. Maydica 2005, 50, 225–234. [Google Scholar]
  8. Berardo, N.; Mazzinelli, G.; Valoti, P.; Redaelli, R. Characterization of maize germplasm for the chemical composition of the grain. J. Agric. Food Chem. 2009, 57, 2378–2384. [Google Scholar] [CrossRef]
  9. Malvar, R.A.; Butrón, A.; Álvarez, A.; Ordás, B.; Soengas, P.; Revilla, P.; Ordás, A. Evaluation of the European Union maize landrace core collection for resistancee to Sesamia nonagrioides (Lepidoptera: Noctuidae) and Ostrinia nubilalis (Lepidoptera: Crambidae). J. Econ. Entomol. 2004, 97, 628–634. [Google Scholar] [CrossRef]
  10. Taba, S.; Twumasi-Afriyie, S. Regeneration guidelines: Maize. In Crop Specific Regeneration Guidelines [CD-ROM]; Dulloo, M.E., Thormann, I., Jorge, M.A., Hanson, J., Eds.; CGIAR System-Wide Genetic Resource Programme: Rome, Italy, 2008; 10p, Available online: https://cropgenebank.sgrp.cgiar.org/images/file/maize/Maize_ENG.pdf (accessed on 14 April 2024).
  11. Simeonovska, E.; Gadžo, D.; Jovović, Z.; Murariu, D.; Kondic, D.; Mandic, D.; Fetahu, S.; Šarčević, H.; Elezi, F.; Prodanović, S.; et al. Collecting local landraces of maize and cereals in South Eastern Europe during 2009 and 2010. Rom. Agric. Res. 2013, 30, 37–43. [Google Scholar]
  12. Gouesnard, B.; Dallard, J.; Panouille, A.; Boyat, A. Classification of French maize populations based on morphological traits. Agronomie 1997, 17, 491–498. [Google Scholar] [CrossRef]
  13. Torri, A.; Lanzanova, C.; Locatelli, S.; Valoti, P.; Balconi, C. Screening of local Italian maize varieties for resistance to Fusarium verticillioides. Maydica 2015, 60, 1–8. [Google Scholar]
  14. Rocha, F.; Bettencourt, E.; Carlos Gaspar, C. Genetic erosion assessment through the re-collecting of crop germplasm. Counties of Arcos de Valdevez, Melgaço, Montalegre, Ponte da Barca and Terras de Bouro (Portugal). Plant Gen. Res. Newsl. 2008, 154, 6–13. [Google Scholar]
  15. Rocha, F.; Gaspar, C.; Barata, A.M. The legacy of collecting missions to the valorisation of agro-biodiversity. Agric. For. 2017, 63, 25–38. [Google Scholar] [CrossRef]
  16. Leitão, S.T.; Ferreira, E.; Bicho, M.C.; Alves, M.L.; Pintado, D.; Santos, D.; Vaz Patto, M.C. Maize open-pollinated populations physiological improvement: Validating tools for drought response participatory selection. Sustainability 2019, 11, 6081. [Google Scholar] [CrossRef]
  17. Murariu, D.; Murariu, M.; Placinta, D.D. Field and laboratory screening of Romanian maize landraces very resistant to low temperatures. Maydica 2015, 60, 11. [Google Scholar]
  18. Murariu, D.; Placinta, D.D.; Simioniuc, D. Assessing genetic diversity in Romanian maize landraces, using molecular markers. Rom. Agric. Res. 2019, 36, 3–9. [Google Scholar] [CrossRef]
  19. Babic, V.; Vancetovic, J.; Prodanovic, S.; Kravić, N.; Babic, M.; Andjelkovic, V. Numerical Classification of Western Balkan Drought Tolerant Maize (Zea mays L.). Landraces. J. Agric. Sci. Technol. 2015, 17, 455–468. [Google Scholar]
  20. Popović, A.; Kravić, N.; Prodanović, S.; Filipović, M.; Sečanski, M.; Babić, V.; Miriţescu, M. Characterisation and evaluation towards selection of maize landraces with the best per se performances. Rom. Agric. Res. 2020, 37, 49–58. [Google Scholar] [CrossRef]
  21. Malvar, R.A.; Butrón, A.; Alvarez, A.; Padilla, G.; Cartea, M.E.; Revilla, P.; Ordás, A. Performance of the European Union Maize Landrace Core Collection for yield under multiple corn borer infestations. Crop Protect. 2007, 26, 775–781. [Google Scholar] [CrossRef]
  22. Schilperoord, P. Kulturpflanzen in der Schweiz—Mais; Verein für alpine Kulturpflanzen: Alveneu, Switzerland, 2017. [Google Scholar] [CrossRef]
  23. Eschholz, T.W.; Peter, R.; Stamp, P.; Hund, A. Swiss maize landraces—Their diversity and genetic relationships. Acta Agron. Hung. 2006, 54, 321–328. [Google Scholar] [CrossRef]
  24. Freitag, N.M. The Swiss Core Set of Maize Landraces and Its Genetic and Phenotypic Diversity. Ph.D. Thesis, ETH, Zurich, Switzerland, 2011. [Google Scholar]
  25. Peter, R.; Eschholz, T.W.; Stamp, P.; Liedgens, M. Swiss maize landraces—Early vigour adaptation to cool conditions. Acta Agron. Hung. 2006, 54, 329–336. [Google Scholar] [CrossRef]
  26. Arca, M.; Mary-Huard, T.; Gouesnard, B.; Berard, A.; Bauland, C.; Combes, V.; Madur, D.; Charcosset, A.; Nicolas, S.D. Deciphering the genetic diversity of landraces with high throughput SNP genotyping of DNA bulks: Methodology and application to the maize 50k array. Front. Plant Sci. 2021, 11, 568699. [Google Scholar] [CrossRef]
  27. Ganal, M.W.; Durstewitz, G.; Polley, A.; Berard, A.; Buckler, E.S.; Charcosset, A.; Clarke, J.D.; Graner, E.-M.; Hansen, M.; Joets, J.; et al. A large maize (Zea mays L.) SNP genotyping array: Development and germplasm genotyping, and genetic mapping to compare with the B73 Reference genome. PLoS ONE 2011, 6, e28334. [Google Scholar] [CrossRef]
  28. Nei, M. Analysis of gene diversity in subdivided populations. Proc. Natl. Acad. Sci. USA 1973, 70, 3321–3323. [Google Scholar] [CrossRef]
  29. Nei, M. F-statistics and analysis of gene diversity in subdivided populations. Ann. Hum. Genet. 1977, 41, 225–233. [Google Scholar] [CrossRef]
  30. Nei, M.; Chesser, R.K. Estimation of fixation indices and gene diversities. Ann. Hum. Genet. 1983, 47, 253–259. [Google Scholar] [CrossRef]
  31. Rogers, J.S. Measures of Genetic Similarity and Genetic Distance. In Studies in Genetics VII; University of Texas Publication 7213; University of Texas: Austin, TX, USA, 1972; pp. 145–153. [Google Scholar]
  32. Mary-Huard, T.; Balding, D. Fast and Accurate Joint Inference of Coancestry Parameters for Populations and/or Individuals. PLoS Genet. 2023, 19, e1010054. [Google Scholar] [CrossRef]
  33. Arca, M.; Gouesnard, B.; Mary-Huard, T.; Le Paslier, M.C.; Bauland, C.; Combes, V.; Madur, D.; Charcosset, A.; Nicolas, S.D. Genotyping of DNA pools identifies untapped landraces and genomic regions to develop next-generation varieties. Plant Biotechnol. J. 2023, 21, 1123–1139. [Google Scholar] [CrossRef]
  34. Gower, J.C. Some distance properties of latent root and vector methods used in multivariate analysis. Biometrika 1966, 53, 325–338. [Google Scholar] [CrossRef]
  35. Paradis, E.; Schliep, K. ape 5.0: An environment for modern phylogenetics and evolutionary analyses in R. Bioinformatics 2019, 35, 526–528. [Google Scholar] [CrossRef]
  36. Alexander, D.H.; Novembre, J.; Lange, K. Fast model-based estimation of ancestry in unrelated individuals. Genome Res. 2009, 19, 1655–1664. [Google Scholar] [CrossRef]
  37. Evanno, G.; Regnaut, S.; Goudet, J. Detecting the number of clusters of individuals using the software Structure: A simulation study. Mol. Ecol. 2005, 14, 2611–2620. [Google Scholar] [CrossRef]
  38. Turlach, B.A. quadprog: Functions to Solve Quadratic Programming Problems. R Package Version 1.5-8. 2019. Available online: https://CRAN.R-project.org/package=quadprog (accessed on 14 April 2024).
  39. Yu, G.; Smith, D.K.; Zhu, H.; Guan, Y.; Lam, T.T. ggtree: An r Package for Visualization and Annotation of Phylogenetic Trees with Their Covariates and Other Associated Data. Methods Ecol. Evol. 2017, 8, 28–36. [Google Scholar] [CrossRef]
  40. Xu, S.; Dai, Z.; Guo, P.; Fu, X.; Liu, S.; Zhou, L.; Tang, W.; Feng, T.; Chen, M.; Zhan, L.; et al. ggtreeExtra: Compact Visualization of Richly Annotated Phylogenetic Data. Mol. Biol. Evol. 2021, 38, 4039–4042. [Google Scholar] [CrossRef]
  41. Xu, S.; Li, L.; Luo, X.; Chen, M.; Tang, W.; Zhan, L.; Dai, Z.; Lam, T.T.; Guan, Y.; Yu, G. Ggtree: A serialized data object for visualization of a phylogenetic tree and annotation data. iMeta 2022, 1, e56. [Google Scholar] [CrossRef]
  42. Wickham, H. ggplot2: Elegant Graphics for Data Analysis; Springer: New York, NY, USA, 2016. [Google Scholar]
  43. Yu, G. scatterpie: Scatter Pie Plot. R Package Version 0.2.1. 2023. Available online: https://CRAN.R-project.org/package=scatterpie (accessed on 14 April 2024).
  44. IBPGR. Descriptors for Maize; International Maize and Wheat Improvement Center, Mexico City/International Board for Plant Genetic Resources: Rome, Italy, 1991. [Google Scholar]
  45. Dubreuil, P.; Charcosset, A. Genetic diversity within and among maize populations: A comparison between isozyme and nuclear RFLP loci. Theor. Appl. Genet. 1998, 96, 577–587. [Google Scholar] [CrossRef]
  46. Brandenburg, J.-T.; Mary-Huard, T.; Rigaill, G.; Hearne, S.J.; Corti, H.; Joets, J.; Vitte, C.; Charcosset, A.; Nicolas, S.D.; Tenaillon, M.I. Independent Introductions and Admixtures Have Contributed to Adaptation of European Maize and Its American Counterparts. PLoS Genet. 2017, 13, e1006666. [Google Scholar] [CrossRef]
  47. Rebourg, C.; Gouesnard, B.; Charcosset, A. Large scale molecular analysis of traditional European maize populations. Relationships with morphological variation. Heredity 2001, 86, 574–587. [Google Scholar] [CrossRef]
  48. Camus-Kulandaivelu, L.; Veyrieras, J.B.; Madur, D.; Combes, V.; Fourmann, M.; Barraud, S.; Dubreuil, P.; Gouesnard, B.; Manicacci, D.; Charcosset, A. Maize adaptation to temperate climate: Relationship between population structure and polymorphism in the dwarf8 gene. Genetics 2006, 172, 2449–2463. [Google Scholar] [CrossRef]
  49. Mir, C.; Zerjal, T.; Combes, V.; Dumas, F.; Madur, D.; Bedoya, C.; Dreisigacker, S.; Franco, J.; Grudloyma, P.; Hao, P.X.; et al. Out of America: Tracing the genetic footprints of the global diffusion of maize. Theor. Appl. Genet. 2013, 126, 2671–2682. [Google Scholar] [CrossRef]
  50. Holker, A.C.; Mayer, M.; Presterl, T.; Bolduan, T.; Bauer, E.; Ordas, B.; Brauner, P.C.; Ouzunova; Melchinger, A.; Schön, C. European maize landraces made accessible for plant breeding and genome-based studies. Theor. Appl. Genet. 2019, 132, 3333–3345. [Google Scholar] [CrossRef]
  51. Bedoya, C.A.; Dreisigacker, S.; Hearne, S.; Franco, J.; Mir, C.; Prasanna, B.M.; Taba, S.; Charcosset, A.; Warburton, M.L. Genetic diversity and population structure of native maize populations in Latin America and the Caribbean. PLoS ONE 2017, 12, e0173488. [Google Scholar] [CrossRef] [PubMed]
  52. Rojas-Barrera, I.C.; Wegier, A.; De Jesús Sánchez González, J.; Owens, G.L.; Rieseberg, L.H.; Piñero, D. Contemporary evolution of maize landraces and their wild relatives influenced by gene flow with modern maize varieties. Proc. Natl. Acad. Sci. USA 2019, 116, 21302–21311. [Google Scholar] [CrossRef]
  53. McLean-Rodríguez, F.D.; Costich, D.E.; Camacho-Villa, T.C.; Pè, M.E.; Dell’Acqua, M. Genetic diversity and selection signatures in maize landraces compared across 50 years of in situ and ex situ conservation. Heredity 2021, 126, 913–928. [Google Scholar] [CrossRef]
  54. Zhang, H.; Liu, X.; Zhang, Y.; Jiang, C.; Cui, D.; Liu, H.; Li, D.; Wang, L.; Chen, T.; Lihua, N. Genetic analysis and fine mapping of the Ga1-S gene region conferring cross-incompatibility in maize. Theor. Appl. Genet. 2012, 124, 459–465. [Google Scholar] [CrossRef]
  55. Gouesnard, B.; Rebourg, C.; Welcker, C.; Charcosset, A. Analysis of photoperiod sensitivity within a collection of tropical maize populations. Genet. Resour. Crop Evol. 2002, 49, 471–481. [Google Scholar] [CrossRef]
  56. Rashid, Z.; Sofi, M.; Harlapur, S.I.; Kachapur, R.M.; Dar, Z.A.; Singh, P.K.; Zaidi, P.H.; Vivek, B.S.; Nair, S.K. Genome-wide association studies in tropical maize germplasm reveal novel and known genomic regions for resistance to Northern corn leaf blight. Sci. Rep. 2020, 10, 21949. [Google Scholar] [CrossRef]
  57. Barbosa, P.A.; Fritsche-Neto, R.; Andrade, M.C.; Petroli, C.D.; Burgueño, J.; Galli, G.; Willcox, M.C.; Sonder, K.; Vidal-Martínez, V.A.; Sifuentes-Ibarra, E.; et al. Introgression of maize diversity for drought tolerance: Subtropical maize landraces as source of new positive variants. Front. Plant Sci. 2021, 12, 691211. [Google Scholar] [CrossRef]
  58. Brandolini, A.; Brandolini, A. Maize introduction, evolution and diffusion in Italy. Maydica 2009, 54, 233–242. [Google Scholar]
  59. Diaw, Y.; Tollon-Cordet, C.; Charcosset, A.; Nicolas, S.; Madur, D.; Ronfort, J.; David, J.; Gouesnard, B. Genetic diversity of maize landraces from the South-West of France. PLoS ONE 2021, 16, e0238334. [Google Scholar] [CrossRef]
  60. Romero Navarro, J.A.; Willcox, M.; Burgueno, J.; Romay, C.; Swarts, K.; Trachsel, S.; Preciado, E.; Terron, A.; Vallejo Delgado, H.; Vidal, V.; et al. A study of allelic diversity underlying flowering-time adaptation in maize landraces. Nat. Genet. 2017, 49, 476–480. [Google Scholar] [CrossRef]
  61. Wang, L.; Josephs, E.B.; Lee, K.M.; Roberts, L.M.; Rellán-Álvarez, R.; Ross-Ibarra, J.; Hufford, M.B. Molecular parallelism underlies convergent highland adaptation of maize landraces. Mol. Biol. Evol. 2021, 38, 3567–3580. [Google Scholar] [CrossRef]
  62. Mendes-Moreira, P.; Satovic, Z.; Mendes-Moreira, J.; Santos, J.P.; Nina Santos, J.P.; Pego, S.; Vaz Patto, M.C. Maize participatory breeding in Portugal: Comparison of farmer’s and breeder’s on-farm selection. Plant Breed. 2017, 136, 861–871. [Google Scholar] [CrossRef]
  63. Wambugu, P.W.; Ndjiondjop, M.N.; Henry, R.J. Role of genomics in promoting the utilization of plant genetic resources in genebanks. Brief. Funct. Genom. 2018, 17, 198–206. [Google Scholar] [CrossRef]
  64. Westengen, O.T.; Skarbø, K.; Mulesa, T.H.; Berg, T. Access to genes: Linkages between genebanks and farmers’ seed systems. Food Secur. 2018, 10, 9–25. [Google Scholar] [CrossRef]
  65. Mascher, M.; Schreiber, M.; Scholz, U.; Graner, A.; Reif, J.C.; Stein, N. Genebank genomics bridges the gap between the conservation of crop diversity and plant breeding. Nat. Genet. 2019, 51, 1076–1081. [Google Scholar] [CrossRef]
  66. Wang, B.; Lin, Z.; Li, X.; Zhao, Y.; Zhao, B.; Wu, G.; Ma, X.; Wang, H.; Xie, Y.; Li, Q.; et al. Genome-wide selection and genetic improvement during modern maize breeding. Nat. Genet. 2020, 52, 565–571. [Google Scholar] [CrossRef]
  67. Mayer, M.; Hölker, A.C.; González-Segovia, E.; Eva Bauer, E.; Presterl, T.; Ouzunova, M.; Melchinger, A.E.; Schön, C.C. Discovery of beneficial haplotypes for complex traits in maize landraces. Nat. Commun. 2020, 11, 4954. [Google Scholar] [CrossRef]
Figure 1. Principal coordinate analysis of 626 maize landraces based on modified Roger’s distance estimated with 23,412 SNP. Different colors represent landraces assigned (>0.6) to one of 9 genetic groups of ADMIXTURE. Symbols represent the country of origin of each landrace (BIH: Bosnia and Herzegovina, CHE: Switzerland, ESP: Spain, FRA: France, HRV: Croatia, ITA: Italy; PRT: Portugal, ROU: Romania, SRB: Serbia).
Figure 1. Principal coordinate analysis of 626 maize landraces based on modified Roger’s distance estimated with 23,412 SNP. Different colors represent landraces assigned (>0.6) to one of 9 genetic groups of ADMIXTURE. Symbols represent the country of origin of each landrace (BIH: Bosnia and Herzegovina, CHE: Switzerland, ESP: Spain, FRA: France, HRV: Croatia, ITA: Italy; PRT: Portugal, ROU: Romania, SRB: Serbia).
Biology 13 00454 g001
Figure 2. Dendrogram of genetic diversity based on modified Roger’s distance (MRD) of 626 landraces including additional data on (from internal to external layer): (i) expected heterozygosity (Hs), (ii) female flowering time/days to silking (DS), and (iii) plant height (PH), with (iv) barplots from penalized regression to seven ancestral groups [33] and (v) barplots from ADMIXTURE analysis (K = 9). Individual accessions are labeled with the landrace code and common name, with label color corresponding to assignment to genetic groups (see also Table S1).
Figure 2. Dendrogram of genetic diversity based on modified Roger’s distance (MRD) of 626 landraces including additional data on (from internal to external layer): (i) expected heterozygosity (Hs), (ii) female flowering time/days to silking (DS), and (iii) plant height (PH), with (iv) barplots from penalized regression to seven ancestral groups [33] and (v) barplots from ADMIXTURE analysis (K = 9). Individual accessions are labeled with the landrace code and common name, with label color corresponding to assignment to genetic groups (see also Table S1).
Biology 13 00454 g002
Figure 3. Map of collection sites of 483 EVA Maize accessions colored according to their membership of nine genetic groups. Each accession is represented by a pie chart where the color composition reflects the proportion of genome ancestry from nine genetic groups in the ADMIXTURE analysis (see legend).
Figure 3. Map of collection sites of 483 EVA Maize accessions colored according to their membership of nine genetic groups. Each accession is represented by a pie chart where the color composition reflects the proportion of genome ancestry from nine genetic groups in the ADMIXTURE analysis (see legend).
Biology 13 00454 g003
Figure 4. Boxplots of phenotypic BLUEs correlated to genetic grouping (K = 9), showing days to tasseling (DT), days to silking (DS), anthesis–silking interval (ASI), plant height in cm (PH), ear height in cm (EH), and relative ear height in % (EPHR) for landraces belonging to each genetic group K1 to K9 (>0.6) or admixed (<0.6), using the same color coding as other figures. Whiskers extend to minimum and maximum values (Q0, Q4), if not considered outliers.
Figure 4. Boxplots of phenotypic BLUEs correlated to genetic grouping (K = 9), showing days to tasseling (DT), days to silking (DS), anthesis–silking interval (ASI), plant height in cm (PH), ear height in cm (EH), and relative ear height in % (EPHR) for landraces belonging to each genetic group K1 to K9 (>0.6) or admixed (<0.6), using the same color coding as other figures. Whiskers extend to minimum and maximum values (Q0, Q4), if not considered outliers.
Biology 13 00454 g004
Figure 5. Representation of the first two components of the principal component analysis for phenotypic data collected from 588 maize landraces of the EVA collection. Different symbols represent country of origin ISO codes (BIH: Bosnia and Herzegovina, CHE: Switzerland, ESP: Spain, FRA: France, HRV: Croatia, ITA: Italy; PRT: Portugal, ROU: Romania, SRB: Serbia), while color represents the genetic structure of landraces (see legend). Hybrid checks are represented by both different symbols and colors. The first principal component (PRIN1) was positively related to plant architecture and days to flowering and the second principal component (PRIN2) was positively related to the anthesis–silking interval (ASI). The proportion of variance explained through each component is also shown.
Figure 5. Representation of the first two components of the principal component analysis for phenotypic data collected from 588 maize landraces of the EVA collection. Different symbols represent country of origin ISO codes (BIH: Bosnia and Herzegovina, CHE: Switzerland, ESP: Spain, FRA: France, HRV: Croatia, ITA: Italy; PRT: Portugal, ROU: Romania, SRB: Serbia), while color represents the genetic structure of landraces (see legend). Hybrid checks are represented by both different symbols and colors. The first principal component (PRIN1) was positively related to plant architecture and days to flowering and the second principal component (PRIN2) was positively related to the anthesis–silking interval (ASI). The proportion of variance explained through each component is also shown.
Biology 13 00454 g005
Figure 6. Dendrogram of pre-cluster and cluster analysis of 588 landraces and 5 hybrid checks based on phenotypic data. Euclidean distance and average linkage were used as a clustering method. The numbers on the left are the sub-clusters (see also Table 7).
Figure 6. Dendrogram of pre-cluster and cluster analysis of 588 landraces and 5 hybrid checks based on phenotypic data. Euclidean distance and average linkage were used as a clustering method. The numbers on the left are the sub-clusters (see also Table 7).
Biology 13 00454 g006
Table 1. Overview of the EVA Maize collection, including criteria for selection. Detailed passport data for all accessions are available in Table S1.
Table 1. Overview of the EVA Maize collection, including criteria for selection. Detailed passport data for all accessions are available in Table S1.
CountryEVA Maize GenebankHolding Institute Code (* WIEWS)Number of AccessionsPart of ** EUMLCCCriteria for Selection of MaterialsReferences
CroatiaUniversity of Zagreb, Faculty of AgricultureHRV04150naRepresentative of national collection based on a wide geographical distribution;
phenological and morphological data.
[11]
FranceINRAE—MontpellierFRA0158016Morphological traits to maximize number of classes;
representative of national collection.
[7,12]
ItalyCREA-Cereal and Industrial Crops BergamoITA3866519Representative of different geographical areas in Italy;
promising sources of stress-tolerance traits.
[8,13]
PortugalINIAV, BragaPRT0014217Morphological, agronomical, and molecular characterization;
previous knowledge from national breeding programmes.
[14,15]
ESAC-IPC, CoimbraPRT0536naDrought tolerance;
yield evaluation at different altitudes.
[16]
RomaniaSuceava GenebankROM00751naEarly and semi-early (FAO groups 200–300);
collected between 1973 and 2000;
cold tolerance.
[17,18]
SerbiaMaize Research Institute Zemun PoljeSRB00191naDrought tolerance;
stability and high grain yield;
good performance in breeding programmes.
[19,20]
SpainMisión Biológica de Galicia, Pontevedra (MBG-CSIC)ESP0042naRepresentative of available diversity;
corn borer resistance;
good for breadmaking;
cold and drought tolerance.
[1,9,21]
ESP0078na
ESP00913723
ESP0164na
ESP119 4na
SwitzerlandAgroscope, NyonCHE00186naRepresentative of the collection from different regions in Switzerland;
phenological and morphological data.
[22,23,24,25]
Total62675
* WIEWS, World Information and Early Warning System on Plant Genetic Resources for Food and Agriculture; ** EUMLCC, European Maize Landraces Core Collection; na, not applicable.
Table 2. Trial locations for phenotypic evaluations in the EVA Maize network 2021–2023.
Table 2. Trial locations for phenotypic evaluations in the EVA Maize network 2021–2023.
Country Location EVA PartnerRange of FAO Maturity Rating# Accessions EvaluatedYears of Evaluation
CroatiaZagrebUniversity of Zagreb100–6001872022–2023
FrancePloudanielINRAE100–8005982021–2023
FranceAlzonneKWS300–5003262021–2023
GermanyBernburg (Saale)KWS100–3001612021–2023
ItalyBergamoCREA—Cereal and Industrial Crops200–5001782021–2023
ItalyMonseliceKWS500–8002182021–2023
PortugalCoimbraESAC400–600532022–2023
RomaniaSuceavaSuceava genebank100–3001082021–2022
SerbiaBelgradeMRIZP300–8002992021–2023
SpainPontevedraMBG-CSIC100–8002162021
SwitzerlandDelleyDelley Semences et Plantes100–250842021–2022
Table 3. Description of traits evaluated in EVA Maize phenotypic trials, following published descriptors [44].
Table 3. Description of traits evaluated in EVA Maize phenotypic trials, following published descriptors [44].
Trait NameTrait AcronymTrait DescriptionUnitCrop Ontology
Days to tasseling (anthesis, male flowering)DTIPGRI descriptor 4.1.1: number of days from sowing to when 50% of the plants have shed pollen.dCO_322:0000030
Days to silking (female flowering)DSIPGRI descriptor 4.1.2: number of days from sowing to when silks have emerged on 50% of the plants.dCO_322:0000031
Anthesis–silking intervalASITime difference [in days] between anthesis and silking, calculated as (ASI = DS − DT).dCO_322:0000001
Plant heightPHIPGRI descriptor 4.1.4: from ground level to the base of the tassel after milk stage, measured in cm.; observed value recorded from an average of 10 plants per plot.cmCO_322:0000994
Ear heightEHIPGRI descriptor 4.1.5: from ground level to the node bearing the uppermost ear after milk stage, measured in cm.; observed value recorded from an average of 10 plants per plot.cmCO_322:0000017
Relative ear heightEPHREar to plant height ratio calculated as (EH/PH) × 100.%
Table 4. Genetic diversity parameters within the EVA panel and each genetic group (K = 9) based on highly assigned landraces (>0.6) from admixture analysis.
Table 4. Genetic diversity parameters within the EVA panel and each genetic group (K = 9) based on highly assigned landraces (>0.6) from admixture analysis.
Whole PanelK1K2K3K4K5K6K7K8K9
Mean ± S.D.Mean ± S.D.Mean ± S.D.Mean ± S.D.Mean ± S.D.Mean ± S.D.Mean ± S.D.Mean ± S.D.Mean ± S.D.Mean ± S.D.
Size of the panel/group626142931403214437539
Number of monomorphic SNPs49567319281161573122229590
Allele number (A) at group level2.00 ± 0.011.96 ± 0.22.00 ± 0.061.99 ± 0.092.00 ± 0.061.99 ± 0.081.98 ± 0.151.99 ± 0.071.99 ± 0.11.97 ± 0.16
Allele number (A) average within landraces1.71 ± 0.131.57 ± 0.131.67 ± 0.091.69 ± 0.061.79 ± 0.081.65 ± 0.111.75 ± 0.051.72 ± 0.071.52 ± 0.111.67 ± 0.09
Minor allele frequency (MAF) at group level0.25 ± 0.140.19 ± 0.150.21 ± 0.150.19 ± 0.150.25 ± 0.140.2 ± 0.150.21 ± 0.150.19 ± 0.150.16 ± 0.150.20 ± 0.16
Minor allele frequency (MAF) average within landraces0.15 ± 0.030.13 ± 0.030.13 ± 0.020.14 ± 0.020.18 ± 0.020.14 ± 0.030.17 ± 0.010.16 ± 0.020.10 ± 0.020.15 ± 0.03
Total expected heterozygosity (Ht) at group level0.33 ± 0.150.27 ± 0.170.29 ± 0.160.26 ± 0.170.33 ± 0.150.27 ± 0.170.28 ± 0.170.27 ± 0.170.22 ± 0.170.27 ± 0.18
Expected heterozygosity (Hs) average within landraces0.21 ± 0.050.17 ± 0.040.19 ± 0.030.19 ± 0.020.24 ± 0.030.19 ± 0.040.23 ± 0.020.22 ± 0.030.14 ± 0.030.20 ± 0.03
Differentiation (Gst) between landraces (within group Ht)0.35 ± 0.140.35 ± 0.140.36 ± 0.100.25 ± 0.090.26 ± 0.100.30 ± 0.130.17 ± 0.070.18 ± 0.100.36 ± 0.150.25 ± 0.13
Modified Roger’s distance (MRD) between landraces0.24 ± 0.040.22 ± 0.040.23 ± 0.030.18 ± 0.030.21 ± 0.030.20 ± 0.030.16 ± 0.020.16 ± 0.020.20 ± 0.040.18 ± 0.04
Table 5. Mean, range, and quartiles of 588 landraces and five check hybrids for field evaluation parameters.
Table 5. Mean, range, and quartiles of 588 landraces and five check hybrids for field evaluation parameters.
* DT* DS* ASI* PH* EH* EPHR
Landraces
Mean72.374.82.91697644.3
Minimum54.251.8−2.21012122.9
Q2567.569.62.01525939.2
Median (Q50)72.275.12.81697444.2
Q7576.279.53.91859048.9
Maximum94.695.67.623815872.1
Check hybrids
Mean78.077.80.32079144.0
Check1 FAO21068.268.60.91997437.4
Check2 FAO30075.475.50.52129845.9
Check3 FAO40078.678.40.32059245.0
Check4 FAO50082.482.30.32099444.8
Check5 FAO60085.484.3−0.72099846.8
* Days to tasseling/anthesis (DT), days to silking (DS), anthesis–silking interval (ASI), plant height in cm (PH), ear height in cm (EH), and relative ear height in % (EPHR).
Table 6. Coefficients of the different traits, Eigenvalues and explained variability for the first (PRIN1), second (PRIN2), and third (PRIN3) principal components.
Table 6. Coefficients of the different traits, Eigenvalues and explained variability for the first (PRIN1), second (PRIN2), and third (PRIN3) principal components.
* TraitsPRIN1PRIN2PRIN3
DT0.450−0.017−0.445
DS0.4520.187−0.370
ASI0.1800.9240.176
PH0.412−0.1240.730
EH0.458−0.2110.254
EPHR0.426−0.224−0.192
Eigenvalue4.3310.97230.3976
% variability explained72167
* Days to tasseling/anthesis (DT), days to silking (DS), anthesis–silking interval (ASI), plant height (PH), ear height (EH), and relative ear height (EPHR).
Table 7. Means of clusters and sub-clusters of 588 landraces and five hybrid checks for field evaluation parameters.
Table 7. Means of clusters and sub-clusters of 588 landraces and five hybrid checks for field evaluation parameters.
ClusterSub- ClusterN* DT* DS* ASI* PH* EH* EPHR
A 292.493.91.9521015371.08
1292.493.91.9521015371.08
B 31274.977.83.341798547.3
24474.178.95.211637043.2
38673.675.13.281828647.0
116672.173.82.111697644.8
141974.977.83.341798547.3
54877.281.54.771929750.1
74977.679.52.311959950.6
C 4484.086.93.2820411756.9
41585.487.72.6517710558.4
82986.990.94.4422713657.8
D 857.155.4−1.321213431.2
6857.155.4−1.321213431.2
E 22766.768.72.501505637.8
94766.467.00.951566240.0
107668.270.82.951616540.5
124666.668.62.461394936.4
152662.363.01.071304334.1
133266.971.34.861474633.3
Total 59372.374.82.941697644.3
* Days to tasseling/anthesis (DT), days to silking (DS), anthesis–silking interval (ASI), plant height in cm (PH), ear height in cm (EH), and relative ear height as % (EPHR). N, number of accessions in each cluster/subcluster.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Balconi, C.; Galaretto, A.; Malvar, R.A.; Nicolas, S.D.; Redaelli, R.; Andjelkovic, V.; Revilla, P.; Bauland, C.; Gouesnard, B.; Butron, A.; et al. Genetic and Phenotypic Evaluation of European Maize Landraces as a Tool for Conservation and Valorization of Agrobiodiversity. Biology 2024, 13, 454. https://doi.org/10.3390/biology13060454

AMA Style

Balconi C, Galaretto A, Malvar RA, Nicolas SD, Redaelli R, Andjelkovic V, Revilla P, Bauland C, Gouesnard B, Butron A, et al. Genetic and Phenotypic Evaluation of European Maize Landraces as a Tool for Conservation and Valorization of Agrobiodiversity. Biology. 2024; 13(6):454. https://doi.org/10.3390/biology13060454

Chicago/Turabian Style

Balconi, Carlotta, Agustin Galaretto, Rosa Ana Malvar, Stéphane D. Nicolas, Rita Redaelli, Violeta Andjelkovic, Pedro Revilla, Cyril Bauland, Brigitte Gouesnard, Ana Butron, and et al. 2024. "Genetic and Phenotypic Evaluation of European Maize Landraces as a Tool for Conservation and Valorization of Agrobiodiversity" Biology 13, no. 6: 454. https://doi.org/10.3390/biology13060454

APA Style

Balconi, C., Galaretto, A., Malvar, R. A., Nicolas, S. D., Redaelli, R., Andjelkovic, V., Revilla, P., Bauland, C., Gouesnard, B., Butron, A., Torri, A., Barata, A. M., Kravic, N., Combes, V., Mendes-Moreira, P., Murariu, D., Šarčević, H., Schierscher-Viret, B., Vincent, M., ... Goritschnig, S. (2024). Genetic and Phenotypic Evaluation of European Maize Landraces as a Tool for Conservation and Valorization of Agrobiodiversity. Biology, 13(6), 454. https://doi.org/10.3390/biology13060454

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