Skip to main content
Download PDF

Selective sweeps identification in distinct groups of cultivated rye (Secale cereale L.) germplasm provides potential candidate genes for crop improvement

BMC Plant Biology volume 23, Article number: 323 (2023) Cite this article

Abstract

Background

During domestication and subsequent improvement plants were subjected to intensive positive selection for desirable traits. Identification of selection targets is important with respect to the future targeted broadening of diversity in breeding programmes. Rye (Secale cereale L.) is a cereal that is closely related to wheat, and it is an important crop in Central, Eastern and Northern Europe. The aim of the study was (i) to identify diverse groups of rye accessions based on high-density, genome-wide analysis of genetic diversity within a set of 478 rye accessions, covering a full spectrum of diversity within the genus, from wild accessions to inbred lines used in hybrid breeding, and (ii) to identify selective sweeps in the established groups of cultivated rye germplasm and putative candidate genes targeted by selection.

Results

Population structure and genetic diversity analyses based on high-quality SNP (DArTseq) markers revealed the presence of three complexes in the Secale genus: S. sylvestre, S. strictum and S. cereale/vavilovii, a relatively narrow diversity of S. sylvestre, very high diversity of S. strictum, and signatures of strong positive selection in S. vavilovii. Within cultivated ryes we detected the presence of genetic clusters and the influence of improvement status on the clustering. Rye landraces represent a reservoir of variation for breeding, and especially a distinct group of landraces from Turkey should be of special interest as a source of untapped variation. Selective sweep detection in cultivated accessions identified 133 outlier positions within 13 sweep regions and 170 putative candidate genes related, among others, to response to various environmental stimuli (such as pathogens, drought, cold), plant fertility and reproduction (pollen sperm cell differentiation, pollen maturation, pollen tube growth), and plant growth and biomass production.

Conclusions

Our study provides valuable information for efficient management of rye germplasm collections, which can help to ensure proper safeguarding of their genetic potential and provides numerous novel candidate genes targeted by selection in cultivated rye for further functional characterisation and allelic diversity studies.

Peer Review reports

Background

During domestication and subsequent diversification and improvement plants were subjected to intensive positive selection. Consequently, several key traits differentiate crop plant from their wild progenitors. In the case of cereal crops, these traits include: larger grain size, loss of natural seed dispersal mechanisms (causing seed retention until harvest), changes in the plant’s architecture (apical dominance), and in plant physiology (changes related to seed dormancy, photoperiodic sensitivity, vernalization requirements) [1, 2]. A number of genes responsible for domestication traits had been already identified and characterized in major crops, such as maize, rice or wheat, for example Q (controlling inflorescence structure in wheat), teosinte branched1 (tb1, controlling shoot architecture in maize), Shattering1 (Sh1, causing the loss of seed shattering), Btr1 and Btr2 (required for the disarticulation of rachis) [1, 3, 4]. Diversification genes, targeted by selection after domestication, are responsible for inter-varietal differences and are typically related to yield, biotic and abiotic stress resistance, grain quality and adaptation [2]. Well know examples of such genes are: maize Y1 gene, related to high carotenoids levels and yellow kernels [5], wheat Rht gene, controlling reduced height, and rice Hd1, controlling flowering time [5,6,7].

At first, the QTL approach was predominately used to identify domestication/improvement loci [8,9,10]. More recently, various population genetic approaches were developed to detect selective sweeps based on genome-wide scans, including population differentiation and environmental association methods [11,12,13,14].

Early studies suggested that several, large effect loci underlie the phenotypic switch from a wild progenitor to a domesticate. More recently, studies based on genome-wide SNPs revealed hundreds of loci showing signatures of selection [11, 15, 16], providing a new insight into the influence of domestication and breeding on the genome and numerous potential candidate genes for crop improvement programs. Nevertheless, despite extensive research on the subject, the knowledge of the influence of domestication and improvement on the genome is still very incomplete in many crops.

Rye (Secale cereale L.) is a cereal closely related to wheat, and an important crop in Central, Eastern and Northern Europe. It is mainly used for the production of flour for bread making, as animal feed, and in distilleries to produce whiskey and vodka. Rye has the highest tolerance of abiotic and biotic stresses (including cold temperature, low soil fertility, and high soil acidity) among the small grain temperate cereals and is a widely used source of genetic variation for wheat improvement [17, 18]. Contrary to its closest crop relatives wheat and barley, cultivated rye is outcrossing, making rye improvement more challenging. For example, self-incompatibility and inbreeding depression occurring in many rye accessions hamper the development of inbred lines and recombinant inbred line populations [19, 20]. Recently genome sequences of two rye accessions were published [21, 22]. Rye genome size ranges from 7.68 to 8.03 Gbp, and the repetitive elements account for 85%-90% of the assemblies.

According to Germplasm Resource Information Network (GRIN), there are four species recognized in the Secale genus: S. cereale, S. strictum, S. vavilovii, S. sylvestre. Many molecular studies indicate, however, that S. vavilovii is a part of S. cereale complex, and postulate a revision of Secale classification [23,24,25,26]. A possible explanation for the discrepancies regarding classification of S. vavilovii was provided by Zohary et al. [27], who proposed that the four complexes within Secale are S. cereale, S. strictum, S. iranicum, and S. sylvestre. “True” S. vavilovii forms belong to S. cereale complex, which is supported be extensive molecular data mentioned above. S. iranicum (Kolbylansky) is poorly know, and was at a point of time erroneously described as S. vavilovii and sent to several germplasm collection under this description causing confusion, with some researches working on the ‘true’ S. vavilovii, and others on S. iranicum, only mistakenly described as S. vavilovii. Thus, the matter of Secale classification is not fully resolved yet.

Rye domestication happened approximately four thousand years ago [27, 28], much later than the domestication of wheat or barley (ca. 10 thousand years ago [29]). Prior to that, rye occurred as a weed in wheat and barley plantations [30]. For these reasons, rye is referred to as a secondary domesticate [28]. There is no consensus regarding the immediate wild progenitor of cultivated rye (S. cereale subsp. cereale). Both S. vavilovii and S. strictum, among others, had been suggested as likely candidates [28, 31]. Central and Eastern Turkey and adjacent regions are reported to be the main centre of diversity of rye wild species [30]. Recent genetic diversity scans indicate, that there is considerable diversity within rye genetic resources and that the current breeding pool is genetically relatively narrow and distant from accessions representing genebank collections. Additionally, no clear correspondence of genetic diversity patterns with geographic origins was observed [23, 32,33,34,35]. An intense germplasm exchange between rye breeding programs in different parts of the world is usually indicated as a possible cause [19, 36,37,38].

The aim of this study was to: i) assess the genetic diversity structure in a diverse collection of 478 rye accessions representing different geographic origins and improvement status based on genome-wide, high-quality GBS markers, ii) identify selective sweeps in established germplasm clusters and iii) indicate potential candidate genes targeted by selection in rye.

Results

GBS (DArTseq) genotyping

In total 79 877 SNP markers (Dataset-1) differentiating 478 rye accessions (Table S1) were identified and 49 977 (62.65%) of them could be aligned to the Lo7 rye reference genome sequence [21], with the markers fairly evenly distributed among chromosomes. Specifically, the percentage of SNPs mapped to individual chromosomes ranged from 11.3 for 1R to 15.6 for 2R. After quality filtering, 12 846 high quality (HQ) SNP markers (Dataset-2) were identified and used for population structure and phylogenetic analyses (Table S2). Of those, 10 607 (82.6%) aligned to the Lo7 genome sequence and spanned 99.7% of the assembly (6.72 Gb). Percentage of HQ SNP markers mapped to individual chromosomes varied between 11.6 for 1R to 17.4 for 2R (Table S3). The average MAF and PIC values for all 12 846 HQ SNPs and 10,607 HQ SNPs mapped to genome sequence were 0.15 and 0.17, respectively. Distribution of MAF and PIC values of 12,846 HQ SNPs is shown in Figure S1.

Population structure

Assignment tests

Our analysis showed that K = 2 explained best the population structure (Figure S2). Using the cut-off value of Q ≥ 70%, 430 and 34 accessions were assigned to populations 1 and 2, respectively, while 13 accessions were classified as admixtures (Figure S3, Table S1). Population 1 comprised all accessions representing S. vavilovii and Secale cereale ssp. included in the study (only one S. c. subsp. dighoricum accession was assigned to population 2), and 14 accessions described as unknown in genebank records. Five S. s. subsp. strictum and three S. s. subsp. anatolicum accessions were also assigned to population 1—ca. 24% and 38% of accessions representing these taxa included in the study, respectively. Population 2 included all ten S. sylvestre accessions and accessions representing S. strictum subspecies: most of S. s. ssp. kuprijanovii accessions (11 accessions, 79%) and also accessions of S. s. ssp. strictum (nine accessions, 43%) and S. s. ssp. anatolicum (three accessions, 38%). The remaining three S. s. ssp. kuprijanovii accessions were classified as admixtures, together with seven S. s. ssp. strictum, two S. s. subsp. anatolicum, and two S. s. subsp. africanum accessions.

We further performed two additional STRUCTURE analyses for the germplasm groups indicated in the initial STRUCTURE analysis. In both cases K = 2 was indicated as the most probable value explaining population structure in these subsets (Figure S4). The first additional STRUCTURE analysis was run on 48 accessions classified initially as population 2 (34 accessions) and admixtures (14 accessions). At Q ≥ 0.7 eleven accessions (ten S. sylvestre accession and one S. s. subsp. strictum accession) were assigned to subpopulation p2_a and the remaining 37 accessions (36 accessions of S. strictum subsp. and one S. c. subsp. dighoricum accession) were assigned to subpopulation p2_b. The second additional STRUCTURE analysis was run on 430 accessions classified initially as population 1. At Q ≥ 0.7 as many as 229 accessions were classified as admixtures, 192 accessions were assigned to subpopulation p1_a, and 9 accessions were assigned to subpopulation p1_b. We found the result of STRUCTURE analysis of subpopulation 1 not conclusive and analysed the genetic structure of the whole set further using PCoA and NJ clustering. The results of both additional STRUCTURE runs are shown in Figure S5 and the subpopulation assignments of the accessions are indicated in Table S1.

Principal Coordinates Analysis

PCoA and STRUCTURE results were in a very good agreement. PCoA clustered the accessions into three groups (Fig. 1). Accessions assigned to population 2 by STRUCTURE were divided into two groups in the PCoA plot – a group containing S. sylvestre accessions and a group containing S. strictum accessions. These groups corresponded, respectively, to subpopulations 2_a and 2_b identified in additional STRUCTURE runs. The third group indicated by PCoA, occupying a relatively small area of diversity space, corresponded to population 1 indicated by STRUCTURE. As expected, landraces were dispersed across a larger plot area than modern and historical cultivars, and thus turned out to be more diverse.

Fig. 1
figure 1

Principal Coordinates Analysis plot showing relationships between 478 rye accessions genotyped with 12 846 SNPs. The same plot labelled differently according to A STRUCTURE-based population assignments and B taxonomy

Full size image

NJ clustering

Three major clusters could be distinguished in the NJ tree showing phylogenetic relationships between accessions: A1, A2, and A3 (Fig. 2, Table S1). Cluster A3 could be further subdivided into four subclusters: A3.1- A3.4. The clustering was in very good agreement with the STRUCTURE and PCoA results. Accessions from population 2 were grouped in clusters A1 and A2 in the NJ tree, corresponding to two smaller groups of accessions visible in the PCoA plot (Figure S6). Accessions assigned to population 1 formed cluster A3. Admixtures were placed in the outer region of cluster A2, adjacent to cluster A3.1 (Fig. 2A).

Fig. 2
figure 2

Neighbor-joining tree based on 12 846 SNP markers showing relationships between 478 rye accessions. The same plot labelled differently wth respect to A STRUCTURE-based population assignments, B taxonomy, and C improvement status

Full size image

Clustering corresponded largely with the taxonomy, too (Fig. 2B). The most divergent cluster (A1) was composed of S. sylvestre accessions, and cluster A2 contained predominantly S. strictum accessions, specifically all S. strictum subsp. kuprijanovii analysed. Cluster A3 contained all cultivated rye S. cereale subsp. cereale and S. vavilovii accessions and almost all wild/weedy S. cereale accessions. Wild/weedy S. cereale accessions and S. vavilovii accessions were dispersed in all A3 subclusters, with exception of subcluster A3.4.

Improvement status of S. cereale subsp. cereale influenced the clustering (Fig. 2C). Inbred lines from a hybrid breeding program formed a separate group (subcluster A3.4), which was divided into two parts corresponding to heterotic pools Carsten and Petkus (A3.4a and A3.4b, respectively). All modern varieties and almost all historical varieties (63 varieties, 87.5%) were located in cluster A3.2. Historical varieties found in other clusters (A3.1 and A3.3) originated mostly from North America. Apart from varieties, cluster A3.2 contained also 94 landraces (58.4% of landraces analysed in this study). Landraces occurred also in clusters A3.1 and A.3.3 (24.8 and 16.8% of landraces analysed, respectively). The clustering of landraces did not correspond strongly with geographic origin. Cluster A3.2 contained the majority of European (including all landraces from the Balkan region and Southern Europe) and Asian landraces analysed and 31% of landraces from the Middle East. Landraces from the Middle East were the largest regional germplasm set included in the study (55 accessions from Turkey and Iran, obtained predominantly from the NSGC genebank (36 accessions), but also from PGRC, PAS BG and IPK – 9, 6, and 4 accessions, respectively). A subset of these landraces (ca. 62% of landraces from the Middle East, mostly Turkish, representing all four genebanks mentioned above) was clearly divergent from the rest and constituted the majority of cluster A3.1. Cluster A3.3 included landraces of various geographic origins: Northern Europe (Finnish and Norwegian), Eastern Europe (Russian), Asia (Afganistan), Western Europe (Germany, Austria and Switzerland) and the Middle East. These landraces were acquired almost exclusively from IPK and NordGen genebanks.

The highest numbers of accessions for this study were obtained from the following four genebanks: IPK (113), NSGC (97), NordGen (56), and PAS BG (52). Accessions obtained from the IPK genebank covered a broad spectrum of diversity within S. cereale/S. vavilovii group (cluster A3) and were dispersed in subclusters A3.1, A3.2 and A3.3, similarly to the accessions from NSGC (Figure S7, Table S1). Accessions from PAS BG were not represented in cluster A3.3 (with exception of a single S. c. subsp. ancestrale accession), while accessions from NordGen were absent from cluster A3.1. Taken together the accessions obtained for the study provided a good representation of diversity within the Secale genus.

AMOVA

A very high degree of differentiation was found between the two subpopulations indicated by STRUCTURE (FST = 0.468, 53% percent of total molecular variance attributed to variation within populations, P < 0.001). A very high degree of differentiation was also found for the three accessions groups indicated by PCoA and NJ clustering (corresponding to the A1, A2, and A3 clusters in the NJ tree), with the proportion of molecular variance explained by the differences among populations equal 55% and pairwise population FST values ranging from 0.428 between accessions groups A2 and A3 to 0.734 between accession groups A1 and A3 (P < 0.001). AMOVA analysis of the six germplasm groups defined based on NJ clustering: A1, A2, A3.1, A3.2, A3.3, and A3.4 attributed 31% of the total molecular variance to the differences among populations and 69% to the differences within populations. There was a very high degree of differentiation between groups A1 and A2 and the remaining germplasm groups (Table S4). For the pairwise comparisons involving A1 the FST values ranged from 0.785 (between A1 and A3.2) to 0.603 (between A1 and A2). For the pairwise comparisons involving A2 the FST values ranged from 0.603 (between A2 and A1) to 0.318 (between A1 and A3.4b). In the remaining germplasm group pairs the degree of differentiation was moderate – pairwise FST values ranged from 0.054 (between groups A3.1 and A3.2) to 0.205 (between groups A3.2 and A3.4). The only exception was the population pair A3.1 and A3.3, where the differentiation was low – pairwise FST = 0.035 (P < 0.001, Table S4). There was a moderate degree of differentiation (FST = 0.092, P < 0.001) between the two groups of lines from hybrid breeding program (A3.4a and A3.4b).

Diversity indices

Rye accessions were assigned to groups based on improvement status, taxonomy, results of population structure and phylogenetics analyses. Summary of SNP marker numbers, He (expected heterozygosity) and Ho (observed heterozygosity) values, and physical map length by germplasm group is given in Table 1, while the information on accessions’ membership in the defined groups is given in Table S1. Table 1 contains also information on germplasm groups that were defined for selective sweep detection (g1 to g4, see below, section ‘Selective sweeps in clusters of cultivated germplasm’). The numbers of mapped HQ SNPs differentiating the defined groups varied from 10 364 (97.7% of the HQ, mapped SNPs polymorphic in the whole set) in the S. strictum group to only 803 markers (7.6%) in the S. sylvestre group. A similar pattern of the chromosomal distribution of SNPs could be observed in most germplasm groups, with the highest proportion of SNPs on chromosomes 2R and 5R and the lowest on chromosomes 1R, 4R and 6R (Figure S8). When accessions were grouped according to taxonomy, a deviation from this pattern was noticeable in the S. sylvestre group, where a relatively low proportion of SNPs originated from 5R. In S. sylvestre and also in S. strictum group a relatively high proportion of SNPs was mapped on 6R. When accessions were grouped according to improvement status, a relatively low percentage of SNPs was observed on 5R and relatively high – on 6R in the wild/weedy group.

Table 1 Summary of SNP marker numbers, He and Ho values, and physical map length by germplasm group
Full size table

He and Ho values were 0.211 and 0.228 for the whole set. Within established germplasm groups He ranged from 0.196 for the A1 group established based on NJ clustering to 0.333 for the taxonomic group S. vavilovii (Table 1). Values of He above 0.3 were also obtained for the group A2 (0.308) and in the sweep detection set g2. The highest Ho values occurred in groups of historic cultivars (0.466) and S. vavilovii (0.463). The lowest Ho values, in the range 0.145–0.033, occurred in the germplasm groups containing inbred lines from hybrid breeding program. A low Ho value of 0.183 was obtained for the group of S. sylvestre.

Values of the diversity indices PS (proportion of polymorphic sites), ϴ (theta), π (nucleotide diversity), and Tajima D’s values computed for the whole set of 478 accessions and for the established germplasm groups are given in Table S5. When the accessions were grouped according to improvement status, the values of diversity indices were the highest in the wild/weedy group, and the lowest in modern cultivars. In taxonomical groups the highest values of diversity indices were observed in S. strictum and the lowest in S. sylvestre. Average Tajima’s D values were negative in each germplasm group indicating the occurrence of positive selection. The weakest negative Tajima’s D values were recorded in breeding lines, suggesting weak selection in this group. In the cultivated germplasm the strongest negative Tajima’s D value was obtained for modern cultivars, followed by historic cultivars, implying a strong selection in this groups. Among taxonomic groups, the strongest negative Tajima’s D value was noted for S. vavilovii.

Selective sweeps in clusters of cultivated germplasm

NJ analysis revealed the presence of genetics clusters in cultivated rye germplasm and this result was supported by the pairwise FST indicating moderate differentiation between these clusters. Thus, our selection inference is based on the following groups of cultivated accessions: g1 (historical and modern cultivars and related landraces from NJ cluster A3.2), g2 (divergent Turkish landraces from cluster A3.1), g3 (diverse landraces from cluster A3.3), and g4 (inbred lines from hybrid breeding program from cluster A3.4). Selective sweeps detection was also performed separately for groups g4a (inbred lines from Carsten heterotic pool) and g4b (inbred lines from Petkus heterotic pool). The information on accessions’ membership in the sweep detection groups is given in Table S1 and shown in Figure S9. At the adopted settings the algorithms used (SweeD, OmegaPlus, and RAiSD) identified 133 outlier positions in the rye genome in common (Table 2, Table S6). The outliers were located within 13 sweep regions ranging in size from 0.84 Mb to 11.76 Mb (Table S6). The number of sweeps per germplasm group ranged from one (in groups g4a and g4b) to four in group g3. The largest number of outliers (55) was identified in group g4. There were no selective sweep regions in common between the analysed germplasm clusters. The putative selective sweeps were dispersed across the rye genome, with the largest number of sweeps (four) detected in chromosome 7R. No sweeps were detected in chromosome 4R (Table 2).

Table 2 Chromosomal location of putative sweep regions and numbers of outliers and candidate genes identified in groups of cultivated rye accessions
Full size table

Candidate genes in the putative selective sweep regions

In total, for all the sweep detection groups, 170 putative candidate genes were found in the Lo7 genome in the vicinity of the outlier positions within the identified sweeps (Table S7). The number of candidate genes per chromosome ranged from 63 on chromosome 7R to nine on chromosome R2. Within germplasm groups the highest number of candidate genes was identified in g4 (inbred lines from hybrid breeding program) – 51 genes, followed by g3 (divergent Turkish landraces) – 45 genes. Within the candidate genes, we identified the ones located in the vicinity of the outlier/s with the highest values of the statistics computed by the sweep detection algorithms. These candidate genes are listed in Table 3.

Table 3 Candidate genes located in the vicinity of the outliers with the highest values of the statistics computed by the sweep detection algorithms
Full size table

Gene ontology (GO) enrichment analysis, performed on the identified putative candidate genes, revealed that 34 GO terms were significantly overrepresented (Table S8), including glycerolipid biosynthetic process, diacylglycerol O-acyltransferase activity, phosphorelay signal transduction system, polygalacturonase activity, embryo sac morphogenesis, and pollen sperm cell differentiation.

Correspondence to previously identified domestication genes

Based on the literature we compiled a list of ten cereal domestication/improvement genes and using BLAST [39] determined the location of their putative homologues in the Lo7 genome (Table S9). We found that these ten genes were located outside of selective sweeps identified in this study. Previously, selective sweep detection was performed in rye based on the Weining rye genome sequence by Li et al. [22]. We attempted to compare the results obtained by Li et al. [22] with our results, which we had obtained using the genomic sequence of Lo7 rye inbred line. For this purpose we identified in the genome sequence of Lo7 sequences homologous to putative candidate genes identified by Li et al. [22]. Several of these homologous sequences were located in the vicinity (less than 5 Mb) of outliers indicated in our study. This coinciding location was found for the following five candidate genes identified by Li et al. [22]: ScWN1R01G158700_LOC_Os01g08320 (auxin and brassinosteroid hormone responses and plant morphogenesis), ScWN2R01G091200_LOC_Os07g47670 (hypoxia signalling, Pi uptake and accumulation), ScWN2R01G169300_LOC_Os10g25130 (regulation of starch storage in endosperm, internode elongation, domestication traits), ScWN5R01G313900_LOC_Os08g41880 (phosphate deficiency adaptations), ScWN7R01G263700_LOC_Os08g44400 (disease resistance, stress response). We also compared our results to those obtained recently in rye by Sun et al. [40]. Many of the putative candidate genes identified in our study represented the same gene families as the putative candidate genes targeted by selection reported by Sun et al. [40], for example: lipase, gibberellin receptor GID1a, pentatricopeptide repeat-containing protein, and leucine-rich repeat protein kinase family protein. However, their genomic locations did not overlap.

Discussion

Genetic diversity within Secale genus and within cultivated ryes

The first aim of this study was a detailed analysis of genetic diversity and genetic structure in a broad collection of diverse rye germplasm and the identification of germplasm groups suitable for detection of selective sweeps. Several studies on rye genetic diversity were carried out to date, using SSR, array-based (DArT), GBS, and, recently, whole genome resequencing data [23, 24, 26, 33, 40, 41]. In the previous genome-wide studies deploying high-density genotyping to analyse rye genetic diversity [21, 24, 26, 40] up to 143 accessions were used. The present study involved the largest number of rye accessions to date (478). We have assembled possibly the most diverse, yet balanced germplasm set. The improvement status of the accessions ranged from wild species, random mating populations to hybrids to inbred lines used in hybrid breeding. Wild accessions, landraces, and cultivars/breeding lines each represented ca. one third of the set. To ensure the best possible representation of genetic diversity the accessions were obtained from multiple genebanks and breeding companies and selected to cover a broad spectrum of geographic origins (Table S1). The accessions derived from genebanks IPK and NSGC turned out to largely overlap with respect to their diversity, while some areas of rye diversity space were not represented by accessions derived from PAS BG and NordGen. However, the number of accessions derived from these genebanks and sampled in this study is too small to justify a suggestion that there are gaps in their rye germplasm collections.

For the detection of SNP variation, the DArTseq genotyping-by-sequencing method was used. This method was previously shown to efficiently target low copy regions of the very large – ca. 8 Gb [21], and highly repetitive (> 90% [22]) rye genome [42]. Similarly, most of DArTseq markers were found to align to intragenic regions in wheat [43]. DArTseq genotyping proved to be a suitable tool for high-density genome-wide genetic diversity studies and for the detection of selection signals [43,44,45].

The HQ SNPs identified in this study for the analysis of genetic diversity and population structure provided good coverage of the rye genome, spanning ca. 99.7% of the reference genome assembly, with similar proportion of markers originating from individual chromosomes (between ca. 11% and 17%, Table S3). Similarly like in other cereal species, such us barley [46], and wheat [43], a larger proportion of HQ SNP markers segregated in the wild than in cultivated rye (S. cereale subsp. cereale) accessions (over 99.9% vs. 75.7%, Table 1). It is consistent with the assumption that domestication and improvement resulted in a decrease in diversity, and that crop wild relatives are a treasure trove of untapped and potentially valuable variation for crop improvement [47, 48].

STRUCTURE analysis suggested the presence of two subpopulation in the analysed collection (K = 2), dividing the set in to two groups—the first consisting of S. sylvestre and S. strictum accessions and the second containing S. vavilovii and S. cereale (both cultivated and weedy) accessions. Only one S. cereale accession (S. c. subsp. dighoricum from the IPK genebank, accession number R 1722) was assigned to population 2, together with S. sylvestre and S. strictum accessions. We suspect, that this accession was misclassified or mislabelled during ex situ conservation. A verification of this hypothesis would require phenotypical evaluation, since in the genebank records there is no additional information available.

Previous studies had shown, that the implementation of Delta K method often indicates K = 2 (suggesting the presence of two subpopulations) as the highest level of hierarchical structure within the analysed germplasm set, even if the structure is more complex. Hence, the use of other methods in conjunction with Delta K is recommended [49,50,51]. Therefore, we further examined relationships between accessions using additional STRUCTURE runs and PCoA and NJ clustering. The results of these analyses were in very good agreement and indicated a more complex structure within Secale genus – three main complexes corresponding to the taxonomy: the most divergent S. sylvestre complex, the S. strictum complex, and the S. vavilovii/S. cereale complex. This outcome agrees largely with results of previous studies on rye genetic diversity [23, 24, 26], and is also well supported by the outcome of AMOVA analysis and pairwise FST values, indicating a very high degree of differentiation between these three germplasm groups.

A large number of S. strictum subsp. samples included in the study (45 accessions, which originated from different genebanks) allowed us to gather novel information concerning its genetic diversity. We revealed a considerable genetic diversity of this genus, as demonstrated by high values of genetic diversity indices and the results of STRUCTURE analysis, with its representatives present in both populations indicated and also classified as admixtures. NJ clustering indicated differences between S. strictum subspecies analysed. S. strictum subsp. kuprijanovii (14 accessions) turned out to be the most homogenous group among S. strictum subspecies and were located exclusively in cluster A2. Accessions of S. strictum subsp. anatolicum and S. strictum subsp. strictum (eight and 21 accessions, respectively) occurred in three clusters of the NJ tree: in cluster A2 (together with S. strictum subsp. kuprijanovii), and also in clusters A3.1 and A3.2. In clusters A3.1 and A3.2 S. strictum subsp. anatolicum and S. strictum subsp. strictum accessions were intermixed with cultivated and weedy accessions of Secale cereale. However, it cannot be excluded that some of this S. strictum samples placed within S. cereale complex are only erroneously described as S. strictum, since a morphological description was not performed within this study. We find that evaluation of morphological characters would be advisable as a part of future molecular studies on rye taxonomy and phylogeny to exclude possible misclassification of some accessions. On the other hand, the diversity of ten S. sylvestre accessions analysed in this study turned out to be very small, indicating a need for a follow-up, more detailed examination, to ensure proper safeguarding of the genetic potential of this wild relative of cultivated rye. In consistence with the earlier molecular reports examining taxonomic relationships within the genus Secale [23, 26, 40], all the 14 S. vavilovii accessions analysed were intermixed with the S. cereale accessions in clusters A3.1, A3.2 and A3.3 of the NJ tree. Referring to the explanation of Zohary et al. [27] on the controversies regarding the taxonomic position of S. vavilovii, this result would suggest that the samples analysed in this study were “true” S. vavilovii forms. Surprisingly, we observed a strong negative Tajima’s D value in this germplasm group which could be a result of an unintentional selection during genebank conservation. Such selection is possible and could be caused by a low germination rate during regeneration. A coanalysis with S. vavilovii samples from a recent collection would be needed to verify this hypothesis. Based on the analysis of our large germplasm set, comprising 58 wild/weedy Secale cereale accessions (S. c. subsp. afghanicum, S. c. subsp. ancestrale, S. c. subsp. dighoricum, S. c. subsp. rigidum, and S. c. subsp. segetale) we were able to confirm the lack of separation between weedy and cultivated forms of Secale cereale, indicated by earlier studies [24, 26, 40] and suggesting a strong gene flow between these two groups.

We detected the presence of genetic clusters within cultivated ryes and the influence of improvement status on the clustering. NJ clustering indicated a genetic distinctiveness of inbreds used in rye hybrid breeding from the remaining cultivated rye accessions, with the two heterotic pools forming separate clusters in the NJ tree. However, the degree of genetic differentiation between heterotic pools, measured by FST value, was moderate. Further, the remaining cultivated accessions formed three major clusters: one containing the majority of cultivars and related landraces, and two clusters comprising mostly rye landraces. Thus, the study confirmed the indication from earlier works [23, 32, 33], that a large portion of the genetic diversity of rye landraces is not represented in rye cultivars, especially in the modern ones. Similar patterns of genetic diversity distribution between landraces and cultivars were reported also for other crops [43, 52]. These patterns reflect strong selection for several key traits, which occurs during breeding and/or initial, historic choices of germplasm for breeding programs. In this study we identified a very distinct group of rye landraces (cluster A3.1) mostly from Turkey. This geographic region is important for rye evolution as the probable area of origin of cultivated rye and the main centre of diversity [27, 30]. Therefore, we postulate that this germplasm group should be of special interest in conservation efforts, and also future allele-mining projects aimed at identification of novel variation for rye improvement [32].

Detection of genome regions targeted by selection

Detection of selection targets was attempted in rye for the first time by Bauer et al. [53], who used the FST outlier approach and the XTX statistic in pairwise comparisons between three germplasm groups. The germplasm groups used consisted of 46 individuals representing rye genetic resources, 38 inbred lines used in rye hybrid breeding, representing the seed parent pool, and 46 inbred lines representing the pollen parent pool. The analyses were done based on genotypic data obtained with the use of Rye600k array. In each comparison numerous outlier markers were identified, which clustered in a few distinct genome regions, and in total 27 putative selection targets – rye orthologues of cloned and functionally characterised rice genes – were found in these regions. Functions of these putative orthologues were related, among others, to plant height, grain size and number, pollen germination ability, other plant development and morphology functions, abiotic and biotic stress, and regulation of various physiological processes. Subsequently, Li et al. [22] identified loci potentially involved in the domestication of rye based on GBS data of 101 accessions reported by Schreiber et al. [24]. Specifically, SNPs differentiating 81 cultivated rye and five S. vavilovii accessions and three selective sweep detection methods (reduction of diversity (ROD), genome-wide scan of fixation index (FST) and cross-population composite likelihood ratio (CP-CLR)) were used. As a result, 11 selective sweeps with the total of three candidate rye genes, related to brassinosteroid signalling, the transition from vegetative to floral development, and also to tillering and grain yield regulation, were detected in common by the three approaches. The number of candidate genes detected by at least one tool ranged from 10 to 21. Recently, based on resequencing data, Sun et al. [40] performed identification of genes targeted by selection during domestication using FST, XP CLR and ROD approaches on cultivated and weedy ryes from a worldwide set of 116 accessions. As a result selective sweeps with 279 candidate genes were identified by at least two methods and included genes related to plant height, disease resistance, tiller number and grain yield, and also genes with shattering-related functions.

In this study, we used a different approach to detect genome regions targeted by positive selection in rye. Each of the defined rye germplasm groups was analysed separately with three algorithms detecting various distinct signatures of selective pressure: SweeD [54], OmegaPlus [55], and RAiSD [56]. Previously, these algorithms were used successfully to identify selective sweeps signals in various plant germplasm sets, such as African rice, maize inbred lines adapted to African highlands, Canadian spring wheat cultivars, and wild strawberries [44, 49, 57, 58]. We performed sweep detection in cultivated accessions from each of the identified clusters within S. cereale complex. To reduce the length of the manuscript and to focus on the most probable selection targets, we reported only the sweeps that were identified in common by the three algorithms used and listed the potential candidate genes located in the vicinity of the detected outliers.

Genome-wide studies on the influence of selection on the crop plant genome demonstrated that numerous loci, scattered across the genome, are targeted by selection pressure [11, 15]. In accordance with those findings, we detected a total of 133 outliers that were located in 13 putative sweeps, dispersed in the rye genome. The lengths of selective sweeps regions identified in this study ranged from 0.84 to 11.76 Mb and were comparable to those reported in rye by Li et al. [22] (2 – 37 Mb) and Sun et al. [40] (0.11–11 Mb).

We were not able to find a clear correspondence between positions of outliers located in this study and the location of known cereal domestication genes orthologs in the rye genome. However, this result is not altogether surprising since the methods used in this study are dedicated to detection of recent and strong positive selection [12]. The overlap between candidate genes for selection reported previously in rye [22, 40] and those identified in this study was very small. Lack of consistency in the results of selective sweep detection studies is frequently encountered in literature [44, 58], and is associated with various factors, including differences in the type of molecular markers used, marker density and genome coverage, size and the diversity of the germplasm set analysed, and the sweep detection method applied [58, 59]. To circumvent some of these limitations and to improve the reliability of sweep detection, it has been suggested to use different sweep detection methods in parallel [60]. We followed this route and used three tools detecting different signatures of selective pressure and reported the overlap of the results from the three tools, providing numerous novel candidate genes targeted by selection in cultivated rye for further studies.

Potential functions of identified candidate genes and implication for rye genetic improvement

Candidate genes under selection identified so far in crop plants often represent functions related to response to environmental stimuli, such as biotic and abiotic stress resistance, plant architecture, seed size and composition [11, 22, 40, 57, 58]. Similarly, many of the putative candidate genes identified in this study are related to the aforementioned traits.

To identify the key candidate genes we used the values of the statistics computed by the algorithms SweeD, OmegaPlus and RAiSD, and we also performed GO term enrichment analysis. One of the significantly enriched GO terms was phosphorelay signal transduction system, which is involved in turning on and off cellular responses to environmental stimuli [61, 62], such as light, cold, or drought. GO enrichment analysis also indicated the GO terms glycerolipid biosynthetic process and diacylglycerol O-acyltransferase activity. Among candidate genes, these two GO terms were represented by several O-acyltransferases WSD1. Three of them were located in the vicinity of an outlier with the highest statistic score, which was identified in group g4a (Carsten heterotic pool). O-acyltransferases WSD1 are involved in cuticular wax biosynthesis, which acts in plants as a protective barrier against biotic and abiotic stresses, including drought [63,64,65]. Genes representing GO term polygalacturonase activity, were also found to be significantly enriched among putative candidates for selection. This group consisted of several pectin lyase-like genes. Multiple biological functions are being attributed to pectin lyases, such as roles as extracellular virulence agents and roles in plant growth and development, including pollen maturation and pollen tube growth [66]. We also observed significant enrichment of GO terms related to plant fertility and reproduction: embryo sac morphogenesis, and pollen sperm cell differentiation. The candidate genes located in the vicinity of most statistically significant outliers included further genes with obvious connections to traits relevant for plant breeding, such as, among others, genes encoding gibberellin-receptor GID1a, GDSL esterases/lipases, pentatricopeptide repeat-containing proteins, TIR-NBS-LRR class disease resistance protein. Gibberellin-receptors GID1 are key elements in gibberellin signal transduction in plants, and therefore play a role in the control of various aspects of plant growth and fertility, including seed germination and biomass production [67, 68]. GDSL esterases/lipases are active during seed germination and play important roles in plant metabolism, growth and development, including seed development. Pentatricopeptide repeat containing proteins are a large family of proteins regulating gene expression at the RNA level [69]. Most of the known Restorer of Fertility (Rf) genes are members of this family [70]. TIR-NBS-LRR genes constitute one of two groups of the nucleotide-binding leucine-rich repeat (NB-LRR) family, comprising most of the plant pathogen resistance genes [71].

Key targets in rye breeding are disease resistance, resilience to drought and heat stress, and yield improvement [72]. Also, broadening of the genetic base had been identified as one of the most important goals in rye hybrid breeding [38, 72]. Our study identified many putative candidate genes, with functions directly connected to the above mentioned target traits. Presently, it is not possible to determine if the indicated candidate genes indeed underlie the desired phenotypes the breeders select for during rye improvement programs. The identified genes present a starting point for further analyses, aimed at elucidating their specific functions in rye. Extensive allelic diversity studies of the putative candidate genes in broad collections of diverse germplasm might provide further verification for the assumption that these genes were strongly affected by selection and display reduced diversity in improved germplasm. Such approaches are also a powerful tool of identifying novel allelic variants of candidate genes [32, 73,74,75], for targeted broadening of diversity in breeding programs.

Conclusions

Based on a genome-wide detailed analysis of genetic diversity structure in a broad collection of diverse rye germplasm, we identified three complexes within the Secale genus: S. sylvestre, S. strictum and S. cereale/vavilovii. We revealed a relatively narrow diversity of S. sylvestre, very high diversity of S. strictum, and signatures of strong positive selection in S. vavilovii. Within cultivated ryes we detected the presence of genetic clusters and the influence of improvement status on the clustering. Rye landraces represent a reservoir of variation for breeding, since a large portion of their diversity is not represented in rye cultivars, and especially a distinct group of landraces from Turkey should be of special interest as a source of untapped variation. Using three sweep detection algorithms we found that signatures of selection are dispersed in the rye genome and identified 170 putative candidate genes targeted by selection in cultivated rye related, among others, to response to various environmental stimuli (such as pathogens, drought, cold), plant fertility and reproduction (pollen sperm cell differentiation, pollen maturation, pollen tube growth), plant growth and biomass production. Our study provides useful information for efficient management of rye germplasm collections, that can help to ensure proper safeguarding of their genetic potential and provides numerous novel candidate genes targeted by selection in cultivated rye for further functional characterisation and allelic diversity studies.

Methods

Plant material

In total 478 rye (Secale sp.) accessions representing different geographic origins and improvement status, belonging to 17 taxonomic units were analysed: 134 wild or weedy accessions, 161 landraces, 75 historic cultivars, 36 modern cultivars, 58 breeding lines and 14 accessions of an unknown improvement status. The germplasm collection used was composed of three sets: (i) germplasm set 1 consisted of 340 accessions obtained as seed samples from several genebanks and breeding companies, (ii) germplasm set 2 consisted of 54 rye inbred lines representing heterotic pools used in rye hybrid breeding at HYBRO Saatzucht [76], and (iii) germplasm set 3 consisted of 84 accessions used by Al-Beyroutiová et al. [26] in a study on phylogenetic relationships in the genus Secale. The accessions from germplasm set 3 were obtained from eight genebanks and consisted of 15 landraces and 69 wild/weedy accessions representing 13 taxa. Further details can be found in [26]. Information on accessions used in the study is listed in Table S1.

DNA isolation and genotyping

For accessions from germplasm set 1, tissue was collected and DNA isolated as described in [32]. Information on DNA isolation for samples from germplasm sets 2 and 3 can be found, respectively, in [76] and [26]. GBS genotyping (DArTseq) was performed at Diversity Array Technology Pty Ltd., Bruce ACT, Australia (http://diversityarrays.com) as described in [77]. For selection of high quality (HQ) SNPs for population structure and genetic diversity analyses, the following thresholds were adopted: reproducibility > 95%, minor allele frequency (MAF) > 0.01, and maximum missing data < 10%. SNP filtering was done with R package dartR [78]. LD filtering was not done prior to population structure and genetic diversity analyses. PIC values were calculated using the formula by Botstein et al. [79].

Population structure and genetic diversity analyses

The number of clusters (K) capturing the major structure in data was identified with STRUCTURE 2.3.4 software [80] using the admixture model and correlated allele frequencies and the following settings: length of the burn-in period 100,000, number of MCMC replications after burn-in: 10,000. For each number of K tested, ranging from 1 to 14, five independent iterations were performed. The Evanno method [51] implemented in the Structure Harvester [81] was used to identify the number of subpopulations (K) explaining the best population structure in the set. Computation of distance matrices was conducted using R packages adegenet [82, 83], stats and dartR [78]. MEGA11 software [84] was used to construct a Neighbor Joining dendrogram. Principal Coordinate Analysis was done using NTSYSpc 2.2 [85]. AMOVA was performed with GenAlEx 6.503 [86, 87]. He and Ho values were calculated for defined germplasm groups using dartR. A variant call format (vcf) file with markers, which aligned to the reference genome and fulfilled the above specified quality criteria, and the corresponding marker scores was converted using TASSEL v5.2.86 [88] to PHILIP interleaved format and used as input for MEGA11 [84] to calculate the number of segregating sites (S), the proportion of polymorphic sites (Ps), Theta (θ), nucleotide diversity (π), and Tajima’s D [89] for each defined germplasm group.

Selective sweep detection and GO term enrichment analysis

Three algorithms: i) SweeD [54], ii) OmegaPlus [55], and iii) RAiSD [56] were used to detect selective sweeps in defined groups of cultivated rye germplasm. SweeD is based on the Site Frequency Spectrum (SFS), and detects the increase of high- and low-frequency derived variants that is expected in the proximity of a beneficial mutation according to the theory of selective sweeps [90]. The OmegaPlus algorithm relies on Linkage Disequilibrium (LD) and detects a pattern of LD levels where high LD levels are found at each side of the beneficial mutation and drop dramatically for loci across the beneficial mutation. RAiSD detects selective sweeps using multiple signatures of a selective sweep: changes in the amount of genetic diversity, the shift of the SFS towards high- and low-frequency derived variants, and the specific pattern of LD levels, while relying on SNP vectors. The analyses were run as described earlier [44, 57, 91] on the vcf input file of Dataset-1. Only physically mapped SNPs located outside centromere regions were included in the analyses. The top 2% of the highest scores detected by all three algorithms were declared putative selective sweeps. Two or more adjacent, partially overlapping sweep regions were assumed to be one selective sweep. The identified outlier positions within each putative sweep (plus/ minus 500 kb) were used to search for candidate genes in the reference genome of the rye inbred line Lo7 and the accompanying annotation file [21]. To determine if any functional classes were over-represented among the candidate genes a GO enrichment analysis was performed using FUNC-E (github.com/SystemsGenetics/FUNC-E, accessed on 17 Nov. 2022) with the P-value criterion of < 0.01.

Availability of data and materials

GBS genotyping scores for 478 rye accessions analysed in this study were submitted to DRYAD data base, DOI: doi:10.5061/dryad.866t1g1vx.

References

  1. Doebley JF, Gaut BS, Smith BD. The molecular genetics of crop domestication. Cell. 2006;127:1309–21.

    Article  CAS  PubMed  Google Scholar 

  2. Abbo S, Pinhasi R, Gopher A, Saranga Y, Ofner I, Peleg Z. Plant domestication versus crop evolution : a conceptual framework for cereals and grain legumes. Trends Plant Sci. 2014;19:351–60.

    Article  CAS  PubMed  Google Scholar 

  3. Lin Z, Li X, Shannon LM, Yeh CT, Wang ML, Bai G, et al. Parallel domestication of the Shattering1 genes in cereals. Nat Genet. 2012;44:720–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Pourkheirandish M, Hensel G, Kilian B, Senthil N, Chen G, Sameri M, et al. Evolution of the grain dispersal system in barley. Cell. 2015;162:527–39.

    Article  CAS  PubMed  Google Scholar 

  5. Palaisa K, Morgante M, Tingey S, Rafalski A. Long-range patterns of diversity and linkage disequilibrium surrounding the maize Y1 gene are indicative of an asymmetric selective sweep. Proc Natl Acad Sci U S A. 2004;101:9885–90.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Yano M, Katayose Y, Ashikari M, Yamanouchi U, Monna L, Fuse T, et al. Hd1, a major photoperiod sensitivity quantitative trait locus in rice, is closely related to the Arabidopsis flowering time gene CONSTANS. Plant Cell. 2000;12:2473–83.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Pearce S, Saville R, Vaughan SP, Chandler PM, Wilhelm EP, Sparks CA, et al. Molecular characterization of Rht-1 dwarfing genes in hexaploid wheat. Plant Physiol. 2011;157:1820–31.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Doebley J, Stec A, Gustus C. teosinte branched1 and the origin of maize: Evidence for epistasis and the evolution of dominance. Genetics. 1995;141:333–46.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Frary A, Nesbitt TC, Frary A, Grandillo S, Van Der Knaap E, Cong B, et al. fw22: A quantitative trait locus key to the evolution of tomato fruit size. Science (80- ). 2000;289:85–8.

    Article  CAS  Google Scholar 

  10. Li C, Zhou A, Sang T. Rice domestication by reducing shattering. Science. 2006;311:1936–9.

    Article  CAS  PubMed  Google Scholar 

  11. Maccaferri M, Harris NS, Twardziok SO, Pasam RK, Gundlach H, Spannagl M, et al. Durum wheat genome highlights past domestication signatures and future improvement targets. Nat Genet. 2019;51:885–95.

    Article  CAS  PubMed  Google Scholar 

  12. Pavlidis P, Alachiotis N. A survey of methods and tools to detect recent and strong positive selection. J Biol Res. 2017;24:7.

  13. De Mita S, Thuillet AC, Gay L, Ahmadi N, Manel S, Ronfort J, et al. Detecting selection along environmental gradients: Analysis of eight methods and their effectiveness for outbreeding and selfing populations. Mol Ecol. 2013;22:1383–99.

    Article  PubMed  Google Scholar 

  14. Ayalew H, Sorrells ME, Carver BF, Baenziger PS, Ma XF. Selection signatures across seven decades of hard winter wheat breeding in the Great Plains of the United States. Plant Genome. 2020;13:1–10.

    Article  Google Scholar 

  15. Stetter MG, Gates DJ, Mei W, Ross-Ibarra J. How to make a domesticate. Curr Biol. 2017;27:R896-900.

    Article  CAS  PubMed  Google Scholar 

  16. Purugganan MD, Fuller DQ. The nature of selection during plant domestication. Nature. 2009;457:843–8.

    Article  CAS  PubMed  Google Scholar 

  17. Rakoczy-Trojanowska M, Bolibok-Brągoszewska H, Myśków B, Dzięgielewska M, Stojałowski S, Grądzielewska A, et al. Genetics and genomics of stress tolerance. In: Rabanus-Wallace MT, Stein N, editors., et al., The Rye Genome, Compendium of Plant Genomes. Springer: Cham; 2021. p. 213–36.

    Google Scholar 

  18. Crespo-Herrera LA, Garkava-Gustavsson L, Åhman I. A systematic review of rye (Secale cereale L.) as a source of resistance to pathogens and pests in wheat (Triticum aestivum L.). Hereditas. 2017;154:14.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Bolibok-Brągoszewska H, Rakoczy-Trojanowska M. Molecular marker based assessment of genetic diversity in rye. In: Ahuja M, Jain S, editors. Genetic diversity and erosion in plants vol 1: Indicators and prevention. Springer International Publishing Switzerland; 2015.

  20. Lundqvist A. Self-incompatibility rye I Genetic control in the diploid. Hereditas. 1956;42:293–348.

    Article  Google Scholar 

  21. Rabanus-Wallace TM, Hackauf B, Mascher M, Lux T, Wicker T, Gundlach H, et al. Chromosome-scale genome assembly provides insights into rye biology, evolution and agronomic potential. Nat Genet. 2021;53:564–73.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Li G, Wang L, Yang J, He H, Jin H, Li X, et al. A high-quality genome assembly highlights rye genomic characteristics and agronomically important genes. Nat Genet. 2021;53:574–84.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Bolibok-Bragoszewska H, Targonska M, Bolibok L, Kilian A, Rakoczy-Trojanowska M. Genome-wide characterization of genetic diversity and population structure in Secale. BMC Plant Biol. 2014;14:184.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Schreiber M, Himmelbach A, Börner A, Mascher M. Genetic diversity and relationship between domesticated rye and its wild relatives as revealed through genotyping-by-sequencing. Evol Appl. 2018:12:66–77.

  25. Targońska-Karasek M, Bolibok-Brągoszewska H, Oleniecki T, Sharifova S, Kopania M, Rakoczy-Trojanowska M. Verification of taxonomic relationships within the genus Secale (Poaceae: Pooideae: Triticeae) based on multiple molecular methods. Phytotaxa. 2018;383:128.

    Article  Google Scholar 

  26. Al-Beyroutiova M, Sabo M, Sleziak P, Dusinsky R, Bircak E, Hauptvogel P, et al. Evolutionary relationships in the genus Secale revealed by DArTseq DNA polymorphism. Plant Syst Evol. 2016;302:1083–91.

    Article  CAS  Google Scholar 

  27. Zohary D, Hopf M, Weiss E. Domestication of Plants in the Old World: The origin and spread of domesticated plants in Southwest Asia, Europe, and the Mediterranean Basin. Oxford University Press; 2012.

  28. Preece C, Livarda A, Christin PA, Wallace M, Martin G, Charles M, et al. How did the domestication of Fertile Crescent grain crops increase their yields? Funct Ecol. 2017;31:387–97.

    Article  PubMed  Google Scholar 

  29. Fuller DQ. Contrasting patterns in crop domestication and domestication rates: Recent archaeobotanical insights from the old world. Ann Bot. 2007;100:903–24.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Behre KE. The history of rye cultivation in Europe. Veg Hist Archaeobot. 1992;1:141–56.

    Article  Google Scholar 

  31. Tang ZX, Ross K, Ren ZL, Yang ZJ, Zhang HY, T C, et al. Secale. In: C K, editor. Wild Crop Relatives: Genomic and Breeding Resources, Cereals. Berlin Heidelberg: Springer; 2011. p. 367–96.

  32. Hawliczek A, Bolibok L, Tofil K, Borzęcka E, Jankowicz-Cieślak J, Gawroński P, et al. Deep sampling and pooled amplicon sequencing reveals hidden genic variation in heterogeneous rye accessions. BMC Genomics. 2020;21:845.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Targońska M, Bolibok-Brągoszewska H, Rakoczy-Trojanowska M. Assessment of genetic diversity in Secale cereale based on SSR markers. Plant Mol Biol Rep. 2016;34:37–51.

    Article  Google Scholar 

  34. Monteiro F, Vidigal P, Barros AB, Monteiro A, Oliveira HR, Viegas W. Genetic distinctiveness of rye in situ accessions from Portugal unveils a new hotspot of unexplored genetic resources. Front Plant Sci. 2016;7:1334.

    Article  PubMed  PubMed Central  Google Scholar 

  35. Sidhu JS, Ramakrishnan SM, Ali S, Bernardo A, Bai G, Abdullah S, et al. Assessing the genetic diversity and characterizing genomic regions conferring Tan Spot resistance in cultivated rye. PLoS ONE. 2019;14:1–22.

    Article  Google Scholar 

  36. Isik Z, Parmaksiz I, Coruh C, Geylan-Su YS, Cebeci O, Beecher B, et al. Organellar genome analysis of rye (Secale cereale) representing diverse geographic regions. Genome. 2007;50:724–34.

    Article  CAS  PubMed  Google Scholar 

  37. Chikmawati T, Skovmand B, Gustafson JP. Phylogenetic relationships among Secale species revealed by amplified fragment length polymorphisms. Genome. 2005;48:792–801.

    Article  CAS  PubMed  Google Scholar 

  38. Fischer S, Melchinger AE, Korzun V, Wilde P, Schmiedchen B, Möhring J, et al. Molecular marker assisted broadening of the Central European heterotic groups in rye with Eastern European germplasm. Theor Appl Genet. 2010;120:291–9.

    Article  PubMed  Google Scholar 

  39. Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;215:403–10.

    Article  CAS  PubMed  Google Scholar 

  40. Sun Y, Shen E, Hu Y, Wu D, Feng Y, Lao S, et al. Population genomic analysis reveals domestication of cultivated rye from weedy rye. Mol Plant. 2022;15:552–61.

    Article  CAS  PubMed  Google Scholar 

  41. Hagenblad J, Oliveira HR, Forsberg NEG, Leino MW. Geographical distribution of genetic diversity in Secale landrace and wild accessions. BMC Plant Biol. 2016;16:1–20.

    Article  Google Scholar 

  42. Borzęcka E, Hawliczek-Strulak A, Bolibok L, Gawroński P, Tofil K, Milczarski P, et al. Effective BAC clone anchoring with genotyping-by-sequencing and Diversity Arrays Technology in a large genome cereal rye. Sci Rep. 2018;8:8428.

    Article  PubMed  PubMed Central  Google Scholar 

  43. Sansaloni C, Franco J, Santos B, Percival-Alwyn L, Singh S, Petroli C, et al. Diversity analysis of 80,000 wheat accessions reveals consequences and opportunities of selection footprints. Nat Commun. 2020;11:1–12.

    Article  Google Scholar 

  44. Ndjiondjop MN, Alachiotis N, Pavlidis P, Goungoulou A, Kpeki SB, Zhao D, et al. Comparisons of molecular diversity indices, selective sweeps and population structure of African rice with its wild progenitor and Asian rice. Theor Appl Genet. 2019;132:1145–58.

    Article  CAS  PubMed  Google Scholar 

  45. Allan V, Vetriventhan M, Senthil R, Geetha S, Deshpande S, Rathore A, et al. Genome-Wide DArTSeq genotyping and phenotypic based assessment of within and among accessions diversity and effective sample size in the diverse sorghum, pearl millet, and pigeonpea landraces. Front Plant Sci. 2020;11:1–20.

    Article  Google Scholar 

  46. Milner SG, Jost M, Taketa S, Mazón ER, Himmelbach A, Oppermann M, et al. Genebank genomics highlights the diversity of a global barley collection. Nat Genet. 2019;51:319–26.

    Article  CAS  PubMed  Google Scholar 

  47. Zhang H, Mittal N, Leamy LJ, Barazani O, Song BH. Back into the wild—Apply untapped genetic diversity of wild relatives for crop improvement. Evol Appl. 2017;10:5–24.

    Article  PubMed  Google Scholar 

  48. Brozynska M, Furtado A, Henry RJ. Genomics of crop wild relatives: Expanding the gene pool for crop improvement. Plant Biotechnol J. 2016;14:1070–85.

    Article  CAS  PubMed  Google Scholar 

  49. Hu Y, Feng C, Yang L, Edger PP, Kang M. Genomic population structure and local adaptation of the wild strawberry Fragaria nilgerrensis. Hortic Res. 2022;9 October 2021.

  50. Janes JK, Miller JM, Dupuis JR, Malenfant RM, Gorrell JC, Cullingham CI, et al. The K = 2 conundrum. Mol Ecol. 2017;26:3594–602.

    Article  PubMed  Google Scholar 

  51. 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–20.

    Article  CAS  PubMed  Google Scholar 

  52. Müller T, Schierscher-Viret B, Fossati D, Brabant C, Schori A, Keller B, et al. Unlocking the diversity of genebanks: whole-genome marker analysis of Swiss bread wheat and spelt. Theor Appl Genet. 2018;131:407–16.

    Article  PubMed  Google Scholar 

  53. Bauer E, Schmutzer T, Barilar I, Mascher M, Gundlach H, Martis MM, et al. Towards a whole-genome sequence for rye (Secale cereale L.). Plant J. 2017;89:853–69.

    Article  CAS  PubMed  Google Scholar 

  54. Pavlidis P, Živković D, Stamatakis A, Alachiotis N. SweeD: Likelihood-based detection of selective sweeps in thousands of genomes. Mol Biol Evol. 2013;30:2224–34.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  55. Alachiotis N, Stamatakis A, Pavlidis P. OmegaPlus: a scalable tool for rapid detection of selective sweeps in whole-genome datasets. Bioinformatics. 2012;28:2274–5.

    Article  CAS  PubMed  Google Scholar 

  56. Alachiotis N, Pavlidis P. RAiSD detects positive selection based on multiple signatures of a selective sweep and SNP vectors. Commun Biol. 2018;1:79.

    Article  PubMed  PubMed Central  Google Scholar 

  57. Wegary D, Teklewold A, Prasanna BM, Ertiro BT, Alachiotis N, Negera D, et al. Molecular diversity and selective sweeps in maize inbred lines adapted to African highlands. Sci Rep. 2019;9:13490.

  58. Semagn K, Iqbal M, Alachiotis N, N’Diaye A, Pozniak C, Spaner D. Genetic diversity and selective sweeps in historical and modern Canadian spring wheat cultivars using the 90K SNP array. Sci Rep. 2021;11:1–16.

    Article  Google Scholar 

  59. Narum SR, Hess JE. Comparison of FST outlier tests for SNP loci under selection. Mol Ecol Resour. 2011;11(SUPPL. 1):184–94.

    Article  PubMed  Google Scholar 

  60. Vasemägi A, Nilsson J, Primmer CR. Expressed sequence tag-linked microsatellites as a source of gene-associated polymorphisms for detecting signatures of divergent selection in Atlantic salmon (Salmo salar L.). Mol Biol Evol. 2005;22:1067–76.

    Article  PubMed  Google Scholar 

  61. Mira-Rodado V. New insights into multistep-phosphorelay (Msp)/ two-component system (tcs) regulation: Are plants and bacteria that different? Plants. 2019;8:1–21.

    Article  Google Scholar 

  62. Oka A, Sakai H, Iwakoshi S. His-Asp phosphorelay signal transduction in higher plants: Receptors and response regulators for cytokinin signaling in Arabidopsis thaliana. Genes Genet Syst. 2002;77:383–91.

    Article  CAS  PubMed  Google Scholar 

  63. Li F, Wu X, Lam P, Bird D, Zheng H, Samuels L, et al. Identification of the wax ester synthase/acyl-coenzyme a:diacylglycerol acyltransferase WSD1 required for stem wax ester biosynthesis in Arabidopsis. Plant Physiol. 2008;148:97–107.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  64. Xia K, Ou X, Gao C, Tang H, Jia Y, Deng R, et al. OsWS1 involved in cuticular wax biosynthesis is regulated by osa-miR1848. Plant Cell Environ. 2015;38:2662–73.

    Article  CAS  PubMed  Google Scholar 

  65. Laskoś K, Myśków B, Dziurka M, Warchoł M, Dziurka K, Juzoń K, et al. Variation between glaucous and non ‑ glaucous near ‑ isogenic lines of rye ( Secale cereale L.) under drought stress. Sci Rep. 2022;12:22486.

  66. Cao J. The Pectin Lyases in Arabidopsis thaliana: Evolution Selection and Expression Profiles. PLoS One. 2012;7:e46944.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  67. Ge W, Steber CM. Positive and negative regulation of seed germination by the Arabidopsis GA hormone receptors, GID1a, b, and c. Plant Direct. 2018;2:e00083.

    Article  PubMed  PubMed Central  Google Scholar 

  68. Wang X, Li J, Ban L, Wu Y, Wu X, Wang Y, et al. Functional characterization of a gibberellin receptor and its application in alfalfa biomass improvement. Sci Rep. 2016;2017(7):1–12.

    Google Scholar 

  69. Manna S. An overview of pentatricopeptide repeat proteins and their applications. Biochimie. 2015;113:93–9.

    Article  CAS  PubMed  Google Scholar 

  70. Melonek J, Stone JD, Small I. Evolutionary plasticity of restorer-of-fertility-like proteins in rice. Sci Rep. 2016;6:1–12.

    Article  Google Scholar 

  71. Lee HA, Yeom SI. Plant NB-LRR proteins: Tightly regulated sensors in a complex manner. Brief Funct Genomics. 2015;14:233–42.

    Article  CAS  PubMed  Google Scholar 

  72. Miedaner T, Laidig F. Hybrid breeding in rye (Secale cereale L.). In: Al-Khayri J, Jain S, Johnson D, editors. Advances in Plant Breeding Strategies: Cereals. Cham: Springer, Cham; 2019. p. 343–72.

    Chapter  Google Scholar 

  73. Marroni F, Pinosio S, Di Centa E, Jurman I, Boerjan W, Felice N, et al. Large-scale detection of rare variants via pooled multiplexed next-generation sequencing: Towards next-generation Ecotilling. Plant J. 2011;67:736–45.

    Article  CAS  PubMed  Google Scholar 

  74. Duitama J, Kafuri L, Tello D, Leiva AM, Hofinger B, Datta S, et al. Deep assessment of genomic diversity in cassava for herbicide tolerance and starch biosynthesis. Comput Struct Biotechnol J. 2017;15:185–94.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  75. Kharabian-Masouleh A, Waters DLE, Reinke RF, Henry RJ. Discovery of polymorphisms in starch-related genes in rice germplasm by amplification of pooled DNA and deeply parallel sequencing. Plant Biotechnol J. 2011;9:1074–85.

    Article  CAS  PubMed  Google Scholar 

  76. Siekmann D, Jansen G, Zaar A, Kilian A, Fromme FJ, Hackauf B. A genome-wide association atudy pinpoints quantitative trait genes for plant height, heading date, grain quality, and yield in rye (Secale cereale L.). Front Plant Sci. 2021;12:1–23.

    Article  Google Scholar 

  77. Targońska-Karasek M, Bolibok-Brągoszewska H, Rakoczy-Trojanowska M. DArTseq genotyping reveals high genetic diversity of polish rye inbred lines. Crop Sci. 2017;57:1906–15.

    Article  Google Scholar 

  78. Gruber B, Unmack PJ, Berry OF, Georges A. dartR: An R package to facilitate analysis of SNP data generated from reduced representation genome sequencing. Mol Ecol Resour. 2018;18:691–9.

    Article  PubMed  Google Scholar 

  79. Botstein D, White RL, Skolnick M, Davis RW. Construction of a genetic linkage map in man using restriction fragment length polymorphisms. Am J Hum Gen. 1980;32:314–31.

    CAS  Google Scholar 

  80. Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics. 2000;155:945–59.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  81. Earl DA, vonHoldt BM. STRUCTURE HARVESTER: a website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv Genet Resour. 2011;4:359–61.

    Article  Google Scholar 

  82. Jombart T. Adegenet: A R package for the multivariate analysis of genetic markers. Bioinformatics. 2008;24:1403–5.

    Article  CAS  PubMed  Google Scholar 

  83. Jombart T. Ahmed I. adegenet 1.3–1: New tools for the analysis of genome-wide SNP data. Bioinformatics. 2011;27:3070–1.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  84. Tamura K, Stecher G, Kumar S. MEGA11: molecular evolutionary genetics analysis version 11. Mol Biol Evol. 2021;38:3022–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  85. Rohlf FJ. NTSYSpc: Numerical Taxonomy System, ver. 2.2. 2008.

    Google Scholar 

  86. Peakall R, Smouse PE. Genalex 6: genetic analysis in Excel. Population genetic software for teaching and research. Mol Ecol Notes. 2006;6:288–95.

    Article  Google Scholar 

  87. Peakall R, Smouse PE. GenAlEx 6,5: genetic analysis in Excel. Population genetic software for teaching and research–an update. Bioinformatics. 2012;28:2537–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  88. Bradbury PJ, Zhang Z, Kroon DE, Casstevens TM, Ramdoss Y, Buckler ES. TASSEL: Software for association mapping of complex traits in diverse samples. Bioinformatics. 2007;23:2633–5.

    Article  CAS  PubMed  Google Scholar 

  89. Tajima F. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics. 1989;123:585–95.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  90. Maynard Smith J, Haigh J. The hitch-hiking effect of a favourable gene. Genet Res. 2007;89:391–403.

    Article  Google Scholar 

  91. Ji F, Ma Q, Zhang W, Liu J, Feng Y, Zhao P, et al. A genome variation map provides insights into the genetics of walnut adaptation and agronomic traits. Genome Biol. 2021;22:1–22.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

Not applicable.

Funding

This research was funded by the Polish National Science Centre grant No. DEC-2014/14/E/NZ9/00285, by Slovak Research and Development Agency (No. APVV-20–0246), and by Slovak Grant Agency VEGA (No. 1/0180/22). The genotyping of rye inbred lines was funded by the German Federal Ministry of Food and Agriculture based on the decision of the Parliament of the Federal Republic of Germany through the Federal Office of Agriculture and Food (Grant No. 2814IP001). The funding bodies played no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.

Author information

Authors and Affiliations

  1. Department of Plant Genetics, Breeding and Biotechnology, Institute of Biology, Warsaw, University of Life Sciences-SGGW, Warsaw, Poland

    Anna Hawliczek, Ewa Borzęcka, Katarzyna Tofil, Piotr Gawroński & Hanna Bolibok-Brągoszewska

  2. Faculty of Electrical Engineering, Mathematics and Computer Science, University of Twente, Enschede, The Netherlands

    Nikolaos Alachiotis

  3. Department of Silviculture, Institute of Forest Sciences, Warsaw University of Life Sciences-SGGW, Warsaw, Poland

    Leszek Bolibok

  4. HYBRO Saatzucht GmbH & Co. KG, Schenkenberg, Germany

    Dörthe Siekmann

  5. Julius Kühn-Institut, Groß Lüsewitz, Germany

    Bernd Hackauf

  6. Department of Botany, Faculty of Natural Sciences, Comenius University in Bratislava, Bratislava, Slovakia

    Roman Dušinský & Miroslav Švec

Authors
  1. Anna Hawliczek
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Ewa Borzęcka
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Katarzyna Tofil
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Nikolaos Alachiotis
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Leszek Bolibok
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Piotr Gawroński
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Dörthe Siekmann
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Bernd Hackauf
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Roman Dušinský
    View author publications

    You can also search for this author in PubMed Google Scholar

  10. Miroslav Švec
    View author publications

    You can also search for this author in PubMed Google Scholar

  11. Hanna Bolibok-Brągoszewska
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

AH, EB, KT, HBB performed the experiments. AH, EB, NA, LB, PG, HBB analysed the data. NA, LB, BH, DS, RD, MS contributed data/analysis tools. HBB conceived, designed and supervised the experiments, wrote the manuscript. AH, NA, BH, DS, RD, MS, HBB reviewed and edited the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Hanna Bolibok-Brągoszewska.

Ethics declarations

Ethics approval and consent to participate

All methods were carried out in accordance with relevant guidelines and regulations.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Figure S1.

Summary of information on 12 846 HQ SNPs polymorphic in 478 diverse rye accessions. A. Histogram of polymorphic information content (PIC) values. B. Histogram of minor allele frequency (MAF) values. Figure S2. Plot of Delta K values for number of assumed subpopulations (K) ranging from 2 to 14. Figure S3. Population structure of 478 rye accessions at K=2 based on 12 846 SNPs. Each accessions is represented by a vertical stripe partitioned into coloured segments with lengths representing the membership fractions in the inferred clusters. The order of accessions in the plot is the same as in the Table S1. Figure S4. Plot of Delta K values for number of assumed subpopulations (K) ranging from 2 to 14. A. Plot based on STRUCTURE analysis of 48 accessions classified as population 2 and admixtures. B. Plot based on STRUCTURE analysis of 430 accessions classified as population 1. Figure S5. Population structure of rye accessions at K=2. Each accessions is represented by a vertical stripe partitioned into coloured segments with lengths representing the membership fractions in the inferred clusters. The order of accessions in the plot is indicated in the Table S1. A. Population structure of 48 rye accession classified as population 2 and admixtures in the initial STRUCTURE analysis. B. Population structure of 430 rye accession classified as population 1 in the initial STRUCTURE analysis. Figure S6. Principal Coordinates Analysis plot showing relationships between 478 rye accessions genotyped with 12846 SNPs with accessions labelled according to the outcome of NJ clustering. Figure S7. Neighbor-joining tree based on 12 846 SNP markers showing relationships between 478 rye accessions with accession labelled according to source. Figure S8. Chromosomal distribution of polymorphic SNPs by germplasm group. A. Accessions grouped according to improvement status. B. Accessions grouped according to taxonomy. C. Accessions grouped according to the outcome of NJ clustering. D. Accessions grouped according to their membership in sweep detection sets. Figure S9. Neighbor-joining tree based on 12 846 SNP markers showing relationships between 478 rye accessions with branch colour indicating cultivated rye accessions belonging to the respective sweep detection set.

Additional file 2: Table S1.

Information on rye accessions used in the study, including genebank accession number, name, source, taxon, improvement status, country of origin, geographic region, group memberships based on STRUCTURE anaysis and NJ clustering, and sweep detection set membership. Table S2. List of 12486 HQ DArTseq markers used in this study, including their sequences and position in the Lo7 reference genome. Table S3. Chromosomal distribution of SNPs by germplasm group. Table S4. Pairwise population FST values for accessions groups based on NJ clustering. Table S5. Values of genetic diversity indices for the collection of 478 rye accessions and each established germplasm group. Table S6. Information on common outliers and sweeps detected by all three methods (SweeD, OmegaPlus and RAiSD) in groups of cultivated rye accessions. Table S7. List of putative candidate genes targeted by selection in cultivated rye. Table S8. Enriched GO terms for putative candidate genes from the selective sweep regions. Table S9. Literature based list of known cereal domestication/improvement genes and locations of their putative homologues in the Lo7 genome.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Hawliczek, A., Borzęcka, E., Tofil, K. et al. Selective sweeps identification in distinct groups of cultivated rye (Secale cereale L.) germplasm provides potential candidate genes for crop improvement. BMC Plant Biol 23, 323 (2023). https://doi.org/10.1186/s12870-023-04337-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12870-023-04337-1

Keywords

BMC Plant Biology

ISSN: 1471-2229

Contact us