BUSCADOR RECOLECTA

Dataset are structured following well-established data formats Darwing Core. Three files are provided. The first file (Don_liz-gec_event_20221201) contains the information of the project, the institution and the description each event (time of occurrence, geographical coordinates, habitat type, etc…); the second file (Don_liz-gec_occ_20221201) contains the information of the occurrences of species recorded in each transect, taxonomic classification, geographoical coordinates of its observation, etc…; and the third file (Don_liz-gec_mof_20221201) provide information of the description of other variables measured during counts., The monitoring of lizards and geckos’ community in Doñana was initiated in 2005 as part of the monitoring program of natural resources and processes. One of the aims of this project was to obtain a temporal and continuous series of data of the presence and abundance of these species to detect changes and trends in their wild populations within the protected area. Lizard and gecko counts are collected annually by members of the monitoring team three times per year in the periods when reptile activity is high (two samplings in March-June and one sampling in September-October), with good environmental conditions (temperature between 17 and 25 ºC and absence of rain or strong wind conditions). The method used to record the presence and abundance of lizard and gecko species (kilometric index: number of individuals/km) are the transect censuses. Individuals are visually searched along seven transects (linear or circular) that are carried out by one trained person on foot. Each transect have an established length, but the length surveyed is different (averaged 1858 ± 51.86 meters). Each transect is located in dunes or mediterranean vegetation habitats, representative of Doñana ecosystems. Two transects runs through wooden footpaths (within the Natural Park and five transects are placed on sand-trails (4 within the National Park and one in Natural Park which contains a small part of wooden footpath). Linear transects have been surveyed in the outward track direction and after waiting 15 minutes it has been again performed on the way back. This information is valid to account for species presence in the area. However, we suggest to choose the first survey performed for statistical analyses that require independence of samples. Eight species can be potentially observed during the samplings: Mediterranean house gecko (Hemidactylus turcicus), Common Wall Gecko (Tarentola mauritanica), Fringe-fingered Lizard (Acanthodactylus erythrurus), Algerian Psammodromus (Psammodromus algirus), Western Psammodromus (Psammodromus occidentalis), Carbonell's wall lizard (Podarcis carbonelli), Andalusian wall lizard (Podarcis vaucheri), Ocellated lizard (Timon lepidus). Other reptile species present in Doñana have not been included in this study as the detection with this method is very low. For instance, to detect species with a fossorial behaviour (the Mediterranean Worm Lizard Blanus cinereus) or those cryptic as adders (Vipera latastei), it is required larger investment of survey which consist in looking under vegetation, burrows or logs (there are not stones in Doñana). Data recorded during the surveys include weather description (cloud cover, temperature, rain, or wind speed), species identification, number of individuals, sex and life stage of the reptiles when possible, as well as time and georeferenced data of the observation. Between 2005-2007 data was registered in Excel file and since 2008 data is recorded with the app CyberTracker (see protocol). The protocol used has been supervised by herpetological researchers and the data have been validated by the members who performed the transects., National Parks Autonomous Agency (OAPN) between 2002-2007; Singular Scientific and Technical Infrastructures from the Spanish Science and Innovation Ministry (ICTS-MICINN); Ministry of Agriculture, Livestock, Fisheries and Sustainable Development from the Regional Government of Andalusia (CAGPDES-JA) since 2007; and Doñana Biological Station from the Spanish National Research Council (EBD-CSIC) since all the study period (2005)., 1. Don_liz-gec_event_20230524: eventID, institutionCode, institutionID, datasetName, continent, country, countryCode, stateProvince, county, locality, eventDate, eventTime, decimalLatitude, decimalLongitude, verbatimCoordinate, habitat, sampleSizeValue, sampleSizeUnit, samplingEffort, samplingProtocol 2. Don_liz-gec_occ_20230524: recordedBy, eventID, occurrenceID, dynamicPropeerties, decimalLatitude, decimalLongitude, basisOfRecord, individualCount, occurrenceStatus, lifeStage, sex, occurrenceRemarks, behavior, kingdom, class, family, scientificName, genus, specificEpithet, scientificNameAuthorship, taxonRank 3. Don_liz-gec_mof_20230524: eventID, measurementID, measurementType, measurementValue, measurementUnit, measurementMethod., Peer reviewed

Dataset are structured following well-established data formats. Three files are provided. The first file (Meta-data) contains the information of each event (time of occurrence, geographical coordinates, Ecosystem, Sampling mehtod, etc…); the second file (Fish) contains the information of the occurrences of fish species recorded in each station, taxonomic classification, etc…; and the third file (Macroinvertebrates) provide information of the occurrences of macroinvertebrates recorded in each station, taxonomic classification, abundance clases, etc…, The monitoring of the macroinvertebrates and fish community in Doñana wetlands was initiated in 2004 as part of the Monitoring Program of Natural Resources and Processes. The aim was to obtain a temporal and continuous series of data in the abundance and distribution of macroinvertebrates and fish species to analyze the evolution of their numbers and estimates biodiversity values. Data were recorded annually between 2004-2019 by more than 2 members of the monitoring team which performed samplings in different locations twice per year in winter-spring and summer seasons when the study sites are flooded. The macroinvertebrates and fishes were sampled at the 140 stations classified according to their location (on either aeolian sands or marshland). Funnel traps were used as a sampling method. Between 5-9 funnel traps were randomly distributed (until 50 cm of depth) in each location, depending of the flooded area and depth. The traps were left for 24 hours and emptied the content into white sorting pans. Individuals were counted and identified until the maximun taxonomic level in the field and realease. During samplings, it was identified 66 and 16 families, of macroinvertebrates and fishes respectively. The most abundances were Notonectidae and Corixidae in macroinvertebrates, and Poecilidae and Cyprinidae in fishes. Data recorded during the surveys included species identification, number of individuals, sex and life stage (pupa, larvae, inmature, adult) of the organisms when possible, as well as the time and georreferenced data of the observation. Between 2004-2007 data was registered in Excel file and since 2008 data was recorded in CyberTracker sequence). The protocol used has been supervised by researchers and the data have been validated by the members who performed the sampling., We acknowledge financial support from National Parks Autonomous Agency (OAPN) between 2002-2007; Singular Scientific and Technical Infrastructures from the Spanish Science and Innovation Ministry (ICTS-MICINN); Ministry of Agriculture, Livestock, Fisheries and Sustainable Development from the Regional Government of Andalusia (CAGPDES-JA) since 2007; and Doñana Biological Station from the Spanish National Research Council (EBD-CSIC) since all the study period (2005)., 1. Metadata: Taxa group, Site ID, Site name, Country, y coordinate, x coordinate, Ecosystem River/lake name, Sampling method, Starting year, Ending year, 1st Name, 1st Mail, 2nd Name, 2nd Mail, 3rd Name, 3rd Mail.-- 2. Fish: Site ID, Sample ID, Sampling date, Taxon name, Taxon ID, Definition of abundance class, Abundance class.-- 3. Macroinvertebrates: Site ID, Sample ID, Sampling date, Taxon name, Taxon ID, 0+, 1+, Adult, All., Peer reviewed

Dataset are structured following well-established data formats. Two files are provided. The first file (icts-rbd-dungBe_event_20221107) contains the information of each event (time of occurrence, geographical coordinates, habitat, sampling effort, etc…); the second file (icts-rbd-dungBe_occ_20221107) contains the information of the occurrences of dung-beetles species recorded in each site, numbers of individual recorded and taxonomic classification., The monitoring of the roller dung-beetles (Scarabaeinidae) in Doñana, southwestern Spain, was initiated in 2004 as part of the Monitoring Program of Natural Resources and Processes. The aim was to obtain a temporal and continuous series of data in the abundance and distribution of two species (Scarabeus sacer and S. cicatricosus) present in the area. Data were recorded annually from 2004 to 2012 by members of the monitoring team which performed one sampling (between May and August) in different habitats (sand dunes, mediterranean schrublands, flooplain meadows, and marshlands). Dung-baited pitfalls traps were used as a method to obtain samples of individuals of the two species. These traps were plastic cilinder of 30 cm diameter x 20 cm high buried on the ground. A baited grill of 2cm x 2 cm mesh rested on top of the trap. Bait was fresh horse or cow feaces (250 g) collected around the area early in the morning the day before of trapping. Five pitfall traps were established at each site separated 15 m each other during 24 hours. Two checking were conducted every 12 hours after baiting to avoid the mortality of individuals. Individual of each species were counted and release after it. Data recorded during the surveys included species identification and number of individuals. Between 2004-2008 data was registered in Excel file and since 2008 data was recorded in CyberTracker sequence. The protocol used has been supervised by researchers and the data have been validated by the members who performed the sampling., We acknowledge financial support from National Parks Autonomous Agency (OAPN) between 2004-2007; Singular Scientific and Technical Infrastructures from the Spanish Science and Innovation Ministry (ICTS-MICINN); Ministry of Agriculture, Livestock, Fisheries and Sustainable Development from the Regional Government of Andalusia (CAGPDES-JA) since 2007; and Doñana Biological Station from the Spanish National Research Council (EBD-CSIC) since all the study period (2004-2012)., 1.icts-rbd-dungBeetles_event_20221107: intitutionCode, institutionID, datasetName, eventID, eventDate, eventTime, continent, country, countryCode, stateProvince, locality, decimalLatitude, decimalLongitude, habitat, eventRemarks, sampleSizeValue, sampleSizeUnit, sampleEffort, dynamicPropertiesEvents, recordyBy and scientificName 2.icts-rbd-dungBeetles_occ_20221107: eventID, OccurrenceID, basisOfRecords, individualCount, kingdom, class, family, scientificName, genus, specificEpithet and scientificNameAuthorship., Peer reviewed

Dataset are structured following well-established data formats. Three files are provided. The first file (Don_macroinv_ev_20221222) contains the information of each event (eventID, event date, geographical coordinates, sample effort, etc…); the second file (Don_macroinv_occ_20221222) contains the information of the occurrences of individuals recorded in each station and its taxonomic classification; and the third file (Don_macroinv_mof_20221222) provide information of biometric variable (weithg) of macroinvertebrates samples recorded in each occurrence., The monitoring of the macroinvertebrates community in Doñana wetlands was initiated in 2004 as part of the Monitoring Program of Natural Resources and Processes. The aim was to obtain a temporal and continuous series of data in the abundance and distribution of macroinvertebrates species to analyze the evolution of their numbers and estimates biodiversity values. Data were recorded annually between 2004-2019 by more than 2 members of the monitoring team which performed samplings in different locations twice per year in winter-spring and summer seasons when the study sites are flooded. The macroinvertebrates were sampled at the 139 stations classified according to their location (on either aeolian sands or marshland). Funnel traps were used as a sampling method. Between 5-9 funnel traps were randomly distributed (until 50 cm of depth) in each location, depending of the flooded area and depth. The traps were left for 24 hours and emptied the content into white sorting pans. Individuals were counted and identified until the maximun taxonomic level in the field and realease. During samplings, it was identified 65 families. The most abundances were Notonectidae and Corixidae. Data recorded during the surveys included species identification, number of individuals, sex and life stage (pupa, larvae, juvenile, adult) of the organisms when possible, as well as the time and georreferenced data of the observation. Between 2004-2007 data was registered in Excel file and since 2008 data was recorded in CyberTracker sequence). The protocol used has been supervised by researchers and the data have been validated by the members who performed the sampling., We acknowledge financial support from National Parks Autonomous Agency (OAPN) between 2002-2007; Singular Scientific and Technical Infrastructures from the Spanish Science and Innovation Ministry (ICTS-MICINN); Ministry of Agriculture, Livestock, Fisheries and Sustainable Development from the Regional Government of Andalusia (CAGPDES-JA) since 2007; and Doñana Biological Station from the Spanish National Research Council (EBD-CSIC) since all the study period (2005)., 1. Don_macroinv_ev_20221222: eventID, intitutionCode, datasetName, eventDate, year, month, day, continent, country, stateProvince, location, localityID, locality, sampleSizeUnit, sampleSizeEffort, DynamicPropiertiesEvent, eventRemarks, recordedBy.-- 2. Don_macroinv_occ_20221222: eventID, occurrenceID, basisOfRecord, individualCount, sex, lifeStage, kingdom, phylum, order, family, genus, specificEpithet, scientificName.-- 3. Don_macroinv_mof_20221222: ocurrenceID, measurementID, measurementValue, measurementUnit, measurementType, measurementAccuracy, measurementMethod., Peer reviewed

Dataset are structured following well-established data formats. Three files are provided. The first file (Don_fish_ev_20221222) contains the information of each event (eventID, event date, geographical coordinates, sample effort, etc…); the second file (Don_fish_occ_20221222) contains the information of the occurrences of fish species recorded in each station, taxonomic classification; and the third file (Don_fish_mof_20221222) provide information of the biometric variable (weight) of fish sample in each occurrence., The monitoring of the fish community in Doñana wetlands was initiated in 2004 as part of the Monitoring Program of Natural Resources and Processes. The aim was to obtain a temporal and continuous series of data in the abundance and distribution of fish species to analyze the evolution of their numbers and estimates biodiversity values. Data were recorded annually between 2004-2019 by more than 2 members of the monitoring team which performed samplings in different locations twice per year in winter-spring and summer seasons when the study sites are flooded. The fishes were sampled at the 139 stations classified according to their location (on either aeolian sands or marshland). Funnel traps were used as a sampling method. Between 5-9 funnel traps were randomly distributed (until 50 cm of depth) in each location, depending of the flooded area and depth. The traps were left for 24 hours and emptied the content into white sorting pans. Individuals were counted and identified until the maximun taxonomic level in the field and realease. During samplings, it was identified 15 families. The most abundances were Poecilidae and Cyprinidae. Data recorded during the surveys included species identification, number of individuals, sex and life stage (pupa, larvae, inmature, mature) of the organisms when possible, as well as the time and georreferenced data of the observation. Between 2004-2007 data was registered in Excel file and since 2008 data was recorded in CyberTracker sequence). The protocol used has been supervised by researchers and the data have been validated by the members who performed the sampling., We acknowledge financial support from National Parks Autonomous Agency (OAPN) between 2002-2007; Singular Scientific and Technical Infrastructures from the Spanish Science and Innovation Ministry (ICTS-MICINN); Ministry of Agriculture, Livestock, Fisheries and Sustainable Development from the Regional Government of Andalusia (CAGPDES-JA) since 2007; and Doñana Biological Station from the Spanish National Research Council (EBD-CSIC) since all the study period (2005)., 1. Don_fish_ev_20221222: eventID, intitutionCode, institutionID, datasetName, eventDate, year, month, day, country, stateProvince, location, localityID, locality, decimalLatitude, decimalLongitude, habitat, sampleSizeUnit, sampleSizeEffort, DynamicPropiertiesEvent, eventRemarks, recordeBy.-- 2. Don_fish_occ_20221222: eventID, occurrenceID, individualCount, sex, lifeStage, kingdom, phylum, order, family, genus, specificEpithet, scientificName.-- 3. Don_fish_mof_20221222: OccurrenceID, measurementID, measurementType, measurementValue, measurementUnit, measurementMethod., Peer reviewed

Dataset are structured following well-established data formats. Three files are provided., The first file (Don_biom_red-swamp-crayfish_ev_20230111) contains the information of each event (eventID, event date, geographical coordinates, sample effort, etc…); the second file (Don_ biom_red-swamp-crayfish _occ_20230111) contains the information of the occurrences of individuals recorded in each station and its taxonomic classification; and the third file (Don_ biom_red-swamp-crayfish _mof_20230111) provide information of biometric variables of individual crayfish recorded (total body length, cephalothorax length, cephalothorax width and weight)., The monitoring of biometric parameters (total body length, cephalothorax length, cephalothorax width and weight) of red-swamp crayfish (Procambarus clarkii) in Doñana wetlands was initiated in 2004 as part of the Monitoring Program of Natural Resources and Processes. The aim was to obtain a temporal and continuous series of data in the biometry of the species to analyze the evolution in these variables. Data were recorded annually between 2004-2022 by more than 2 members of the monitoring team which performed samplings in different locations twice per year in winter-spring and summer seasons when the study sites are flooded. The individuals were sampled at the 64 stations classified according to their location (on either aeolian sands or marshland). Modified commercial traps, similar to funnel traps, with 4 mm mesh size were used as sampling method, following the method by Bravo et. al., 1994. Between 5-9 traps were randomly distributed (until 50 cm of depth) in each location, depending of the flooded area and depth. The traps were left for 24 hours and emptied the content into white sorting pans. Individuals were counted and identified until the maximun taxonomic level in the field and realease. Each individual captured was weigthed using an electronic balance within 0.2 g accuracy. Cephalothorax length (from the tip of the rostrum to the tip of the cephalatorax) and cephalothorax width (maximum width of the cephalothorax) were measured using a slide caliper within 0.1 mm accuracy. Total body length (from the tip of the rostrum to the tip of the telson) were recorded using a ruler within 0.5 mm accuracy. Data recorded during the surveys also included sex and maturation stage (inmature, mature) of the organisms when possible, as well as the time and georreferenced data of the capture. Between 2004-2007 data was registered in Excel file and since 2008 data was recorded in CyberTracker sequence. The protocol used has been supervised by researchers and the data have been validated by the members who performed the sampling. Bravo, M. A.; Duarte, C. M. & Montes, C. (1994) Environmental factors controlling the life history of Procambarus clarkii (Decapoda, Cambaridae) in a temporary marsh of the Doñana National Park (SW Spain). January 1994. Limnology 25(4):2450-2453. Romaire, R. P. & Lutz, C. G., 1989 Population dynamics of Procambarus clarkii (Girard) and Procambarus acutus acutus (Girard) (Decapoda: Cambaridae) in commercial ponds. Aquaculture 81(3-4):253-274., National Parks Autonomous Agency (OAPN) between 2002-2007; Singular Scientific and Technical Infrastructures from the Spanish Science and Innovation Ministry (ICTS-MICINN); Ministry of Agriculture, Livestock, Fisheries and Sustainable Development from the Regional Government of Andalusia (CAGPDES-JA) since 2007; and Doñana Biological Station from the Spanish National Research Council (EBD-CSIC) since all the study period (2004)., 1. Don_biom_red-swamp-crayfish_ev_20230111: eventID, intitutionCode, institutionID, datasetName, eventDate, eventTime, continent, country, stateProvince, location, localityID, locality, decimalLatitude, decimalLongitude, habitat, sampleSizeValue, sampleSizeUnit, sampleSizeEffort. 2. Don_ biom_red-swamp-crayfish _occ_20230111: eventID, occurrenceID, occurenceTime, individualCount, basisOfRecord, sex lifestage, kingdom, phylum, class, family, genus, specificEpithet, scientificName, scientificNameAuthorship. 3. Don_ biom_red-swamp-crayfish _mof_20230111: ocurrenceID, measurementID, measurementValue, measurementUnit, measurementType, measurementAccuracy, measurementMethod., Peer reviewed

Dataset are structured following well-established data formats Darwing Core. Three files are provided. The first file (icts_small_mam_ev_20230515) contains the information of each event (eventID, event date, geographical coordinates, sample effort, etc…). The eventID code has been built with a abbrevation of the plot, subplot, trap (eventRemarks) and eventDate. The information of the treatments applied in the subplots is included in eventRermarks. The second file (Don_small_mam_occ_20230515) contains the information of the occurrences of individuals recorded in each trap, occurrence remarks, ear tag codes of individual captured and its taxonomic classification. The occurrenceID has built with the eventID code plus a number regarding the number of observations recorded for this event.; and the third file (Don_small_mam_mof_20230515) provide information of individual weight captured in the traps. The measurementID has built with the occurrenceID code plus a number regarding the number of measurements recorded for this occurrence., All the data provided in the dataset by the authors are own data., To record the abundance, density and distribution of small-mammals in Doñana Natural Area, a protocol of capture-mark and recapture monitoring of individuals is applied using Sherman’s live-traps (https://www.semice.org/es/el-proyecto/; Flowerdew et al. 2004, Torre el al., 2018). Seven localities of Doñana (MANECORRO, MATASGORDAS, MEDIANA, MORO, LORO AND SABINAR) of different mediterranean habitats (grassland, shrubland, pine forest, juniper woodlands and oak forest) have been sampled twice per year (between May-Dec) since 2011. Each locality has been divided in plots and subplots respectively: MANECORRO3 (A, B, C, D), MATASGORDAS3 (B, C, D), MATASGORDAS4 (A), MATASGORDAS5 (A), MEDIANA2 (B, C, D), MEDIANA3 (A, C, D), MORO1 (B, D), MORO2 (A, C), LORO1 (A, B, D), LORO2 (A, B, C), LORO4 (B, C, D), LORO6 (A, C, D), SABINAR1 (A), SABINAR2 (A) and SABINAR3 (A). In each subplot has been placed between 36 traps spaced 15 meters apart covering an area of 0.56 ha. The small mammals species captured have been: white-toothed shrew Crocidura russula, black rat Rattus rattus, garden dormouse Eliomys quercinus, algerian mouse Mus spretus and wood mouse Apodemus sylvaticus. Traps were installed during the morning and remained active for three consecutive nights. A daily check was made in the morning (a night check is unnecessary since Sherman traps include food, synthetic cotton inside and lined with film to withstand inclement weather). Once captured, individual was identified until species level and sexed when it was posible. Additionally, body mass of indviduals were taken using a pesola with 0.5 g accuracy. To control the recaptures individuals of R. rattus, A. sylvaticus, M. spretus, and E. quercinus were marked with a ear tag. On the other hand, individuals of E. quercinus also was marked with transponder. Individuals of C. russula were marked with hair cut or due to its high metabolism and risk of death was released without performing the procedures at moment. All animals were released in the same place of capture. The trapping method was consistent throughout the study period. The study is divided in two periods (2011-2016 and 2017-2021) with different sampling effort and different treatments applied in the study plots, since the study belong two different projects (ICTS-RBD 2011-2021, LIFE-ADAPTAMED 2017-2021). LORO1 and LORO2 was affected by a wildfire in June of 2017. The different treatments applied in the study area are the followings: 1.- Grazing exclusion using fences. 2.- Planting with protection of heap of branches. 3.- Planting without protection of heap of branches. 4.- Plant with and without protection of heap of branches. 5.- Cutting 60% of vegetation. 6.- Not treatment., We acknowledge financial support from Singular Scientific and Technical Infrastructures from the Spanish Science and Innovation Ministry (ICTS-MICINN); the Ministry of Agriculture, Livestock, Fisheries and Sustainable Development from the Regional Government of Andalusia (CAGPDES-JA) since 2011; the Doñana Biological Station from the Spanish National Research Council (EBD-CSIC) since all the study period (2011); the Ministry of Environmetal sustainability and blue economy from the Regional Goverment of Andalusia since 2017 with the LIFE-ADAPTAMED project,; Ministry of Science and Innovation (Recovery, Transformation and Resilence Plan); and the European Comision with the Long-term Ecosystem Research in Europe (eLTER) (a HORIZON funding coordination of the European funding programme for research and innovation) and NextGenerationEU funding., 1.icts-rbd_smallma_ev_20230515: institutionCode, institutionID, datasetName, eventID, year, month, day, eventDate, continent, country, stateProvince, county, locality, decimalLatitude, decimalLongitude, habitat, eventRemarks, sampleSizeValue, samplingSizeUnit, samplingEffort, samplingProtocol.-- 2. icts-rbd_smallma_occ_20230515: eventID, occurrenceID, occurrenceRemarks, basisOfRecord, individualCount, sex, lifeStage, kingdom, class, family, genus, specificEpithet, scientificName, scientificNameAuthorship, dynamicProperties, recordedBy.-- 3. icts-rbd_smallma_mof_20230515: occurrenceID, measurementID, measurementValue, measurementUnit, measurementType, measurementMethod., Peer reviewed

Table S6 Candidate genes mapped within the genomic regions associated with the significantly associated SNPs detected in response to U. pisi. Chromosomal linkage disequilibrium (LD) decay was considered to limit the genomic regions were to look for candidate genes., Uromyces pisi ([Pers.] D.C.) Wint. is an important foliar biotrophic pathogen infecting grass pea (Lathyrus sativus L.), compromising their yield stability. To date, few efforts have been made to assess the natural variation in grass pea resistance and to identify the resistance loci operating against this pathogen, limiting its efficient breeding exploitation. To overcome this knowledge gap, the genetic architecture of grass pea resistance to U. pisi was investigated using a worldwide collection of 182 accessions through a genome-wide association approach. The response of the grass pea collection to rust infection under controlled conditions and at the seedling stage did not reveal any hypersensitive response but a continuous variation for disease severity, with the identification of promising sources of partial resistance. A panel of 5,651 high-quality single-nucleotide polymorphism (SNP) markers previously generated was used to test for SNP-trait associations, based on a mixed linear model accounting for population structure. We detected seven SNP markers significantly associated with U. pisi disease severity, suggesting that partial resistance is oligogenic. Six of the associated SNP markers were located in chromosomes 4 and 6, while the remaining SNP markers had no known chromosomal position. Through comparative mapping with the pea reference genome, a total of 19 candidate genes were proposed, encoding for leucine-rich repeat, NB-ARC domain, and TGA transcription factor family, among others. Results presented in this study provided information on the availability of partial resistance in grass pea germplasm and advanced our understanding of the molecular mechanisms of quantitative resistance to rust in grass pea. Moreover, the detected associated SNP markers constitute promising genomic targets for the development of molecular tools to assist disease resistance precision breeding., Peer reviewed

Table S5. Inflation factors observed for the models not accounting (Naïve) or accounting for population structure (Eigen), and familiar relatedness (Kinship) tested in the GWAS related to the inoculation experiments with U. pisi., Uromyces pisi ([Pers.] D.C.) Wint. is an important foliar biotrophic pathogen infecting grass pea (Lathyrus sativus L.), compromising their yield stability. To date, few efforts have been made to assess the natural variation in grass pea resistance and to identify the resistance loci operating against this pathogen, limiting its efficient breeding exploitation. To overcome this knowledge gap, the genetic architecture of grass pea resistance to U. pisi was investigated using a worldwide collection of 182 accessions through a genome-wide association approach. The response of the grass pea collection to rust infection under controlled conditions and at the seedling stage did not reveal any hypersensitive response but a continuous variation for disease severity, with the identification of promising sources of partial resistance. A panel of 5,651 high-quality single-nucleotide polymorphism (SNP) markers previously generated was used to test for SNP-trait associations, based on a mixed linear model accounting for population structure. We detected seven SNP markers significantly associated with U. pisi disease severity, suggesting that partial resistance is oligogenic. Six of the associated SNP markers were located in chromosomes 4 and 6, while the remaining SNP markers had no known chromosomal position. Through comparative mapping with the pea reference genome, a total of 19 candidate genes were proposed, encoding for leucine-rich repeat, NB-ARC domain, and TGA transcription factor family, among others. Results presented in this study provided information on the availability of partial resistance in grass pea germplasm and advanced our understanding of the molecular mechanisms of quantitative resistance to rust in grass pea. Moreover, the detected associated SNP markers constitute promising genomic targets for the development of molecular tools to assist disease resistance precision breeding., Peer reviewed

Table S4. Averaged DS scores (±standard error of the mean) measured from inoculation experiments with U. pisi, in accessions grouped based on their geographical origin, seed color, and seed size. The small letter represents significant differences (P-value ≤ 0.05) among accessions grouped by origin, seed size or color., Uromyces pisi ([Pers.] D.C.) Wint. is an important foliar biotrophic pathogen infecting grass pea (Lathyrus sativus L.), compromising their yield stability. To date, few efforts have been made to assess the natural variation in grass pea resistance and to identify the resistance loci operating against this pathogen, limiting its efficient breeding exploitation. To overcome this knowledge gap, the genetic architecture of grass pea resistance to U. pisi was investigated using a worldwide collection of 182 accessions through a genome-wide association approach. The response of the grass pea collection to rust infection under controlled conditions and at the seedling stage did not reveal any hypersensitive response but a continuous variation for disease severity, with the identification of promising sources of partial resistance. A panel of 5,651 high-quality single-nucleotide polymorphism (SNP) markers previously generated was used to test for SNP-trait associations, based on a mixed linear model accounting for population structure. We detected seven SNP markers significantly associated with U. pisi disease severity, suggesting that partial resistance is oligogenic. Six of the associated SNP markers were located in chromosomes 4 and 6, while the remaining SNP markers had no known chromosomal position. Through comparative mapping with the pea reference genome, a total of 19 candidate genes were proposed, encoding for leucine-rich repeat, NB-ARC domain, and TGA transcription factor family, among others. Results presented in this study provided information on the availability of partial resistance in grass pea germplasm and advanced our understanding of the molecular mechanisms of quantitative resistance to rust in grass pea. Moreover, the detected associated SNP markers constitute promising genomic targets for the development of molecular tools to assist disease resistance precision breeding., Peer reviewed