Email updates

Keep up to date with the latest news and content from BMC Genomics and BioMed Central.

Open Access Research article

DArT markers for the rye genome - genetic diversity and mapping

Hanna Bolibok-Brągoszewska1*, Katarzyna Heller-Uszyńska2, Peter Wenzl2, Grzegorz Uszyński2, Andrzej Kilian2 and Monika Rakoczy-Trojanowska1

Author Affiliations

1 Department of Plant Genetics, Breeding and Biotechnology, Warsaw University of Life Sciences (SGGW), 159 Nowoursynowska Str, 02-776 Warsaw, Poland

2 Diversity Arrays Technology P/L, 1 Wilf Crane Crescent, Yarralumla, Canberra, ACT, 2600, Australia

For all author emails, please log on.

BMC Genomics 2009, 10:578  doi:10.1186/1471-2164-10-578

The electronic version of this article is the complete one and can be found online at:

Received:27 November 2008
Accepted:3 December 2009
Published:3 December 2009

© 2009 Bolibok-Brągoszewska et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.



Implementation of molecular breeding in rye (Secale cereale L.) improvement programs depends on the availability of high-density molecular linkage maps. However, the number of sequence-specific PCR-based markers available for the species is limited. Diversity Arrays Technology (DArT) is a microarray-based method allowing for detection of DNA polymorphism at several thousand loci in a single assay without relying on DNA sequence information. The objective of this study was the development and application of Diversity Arrays technology for rye.


Using the PstI/TaqI method of complexity reduction we created a rye diversity panel from DNA of 16 rye varieties and 15 rye inbred lines, including parents of a mapping population consisting of 82 recombinant inbred lines. The usefulness of a wheat diversity panel for identification of DArT markers for rye was also demonstrated. We identified 1022 clones that were polymorphic in the genotyped ILs and varieties and 1965 clones that differentiated the parental lines L318 and L9 and segregated in the mapping population. Hierarchical clustering and ordination analysis were performed based on the 1022 DArT markers to reveal genetic relationships between the rye varieties and inbred lines included in the study. Chromosomal location of 1872 DArT markers was determined using wheat-rye addition lines and 1818 DArT markers (among them 1181 unique, non-cosegregating) were placed on a genetic linkage map of the cross L318 × L9, providing an average density of one unique marker every 2.68 cM. This is the most saturated rye linkage map based solely on transferable markers available at the moment, providing rye breeders and researches with a better choice of markers and a higher probability of finding polymorphic markers in the region of interest.


The Diversity Arrays Technology can be efficiently and effectively used for rye genome analyses - assessment of genetic similarity and linkage mapping. The 11520-clone rye genotyping panel with several thousand markers with determined chromosomal location and accessible through an inexpensive genotyping service is a valuable resource for studies on rye genome organization and in molecular breeding of the species.


Winter rye (Secale cereale L.) is an important cereal crop in Central, Eastern and Northern Europe, mainly due to its exceptional ability to thrive and to produce high yields in adverse environmental conditions. Nutrient efficiency and tolerance of diseases exhibited by rye, allowing for a reduced usage of pesticides and fertilizers during production, increase its attractiveness for farmers and consumers. The recently better recognized dietary value of rye contributed to the noticeable growth of consumers' interest in rye products. Breeding progress in rye, however, is rather slow, as the traditional breeding is hampered by the outcrossing nature of this crop, self-incompatibility and occurrence of inbreeding depression.

Implementation of molecular breeding in rye improvement programs depends on the availability of high-density molecular linkage maps. Several genetic maps of rye have been published so far [1-4], but the possibilities of their practical application are rather limited, mainly due to an insufficient saturation. The average density of the most saturated rye linkage map published to date is 2.9 cM. However, there are seven gaps greater than 20 cM on the map [5]. The average map interval length on remaining rye maps exceeds 4.0 cM [2], in the majority of maps it is greater than 7 cM, i.e. 7.4 [6], 7.9 [4], 9.9 [7], 12.2 [8]. On the other hand, as the number of sequence-specific, PCR-based markers (SSRs and SCARs) available for the species was until recently below 300 [4,6,9-11], rye genetic maps were constructed predominantly with the use of usually poorly-transferable AFLP and RAPD markers, or labor intensive, time consuming and low-throughput RFLP markers. The highest number of sequence-specific, PCR-based markers on a rye linkage map equals 58 [4]. Lately, with the creation of BAC library specific for 1RS [12] and the development of 74 SSR markers for 1RS [13], a significant advance took place in rye genomics. However, the knowledge about the remaining part of the rye genome remains very limited.

The Diversity Arrays Technology (DArT) - a microarray, hybridization-based platform - has a capacity to deliver several thousand of sequence-specific markers without relying on sequence information [14]. To date, the performance of the method was validated in several species including cereals such as barley (Hordeum vulgare L.) [15], wheat (Triticum aestivum L.) [16] and sorghum (Sorghum bicolor (L.) Moench) [17]. The current list of species for which DArT arrays are available commercially as service is at webcite and the most comprehensive review of the technology was presented by Kilian et al [18]. The technology finds increasing use in creating high density genetic maps [19] and association studies [20], but it is also expanding into the area of plant and animal biodiversity and population genetics [21,22]. The objective of this study was (i) to develop DArT markers for rye, (ii) to test their usefulness for assessing genetic similarity in rye inbred lines and varieties and (iii) to use DArT markers for constructing a high-density linkage map of the cross L318 × L9 [7].

Results and Discussion

Applicability of the PstI/TaqI complexity reduction method for the rye genome

Complexity reduction is a process, which reproducibly selects a defined fraction of genomic fragments and is a critical step in DArT. The PstI/TaqI method is one of routinely used methods of genome complexity reduction in DArT assays and was shown to be superior to other method tested in wheat [16] and barley [15]. Consequently, this method was chosen for preparing targets form rye DNA in the preliminary experiments with a wheat diversity panel. The performance of PstI/TaqI generated rye targets during hybridizations was satisfactory, as demonstrated by good values of quality parameters, and for this reason the PstI/TaqI method of genome complexity reduction was used in all subsequent DArT analyses performed. Additionally, the preliminary hybridizations revealed wheat DArT markers that differentiated rye ILs. In total 768 candidate markers were identified on the wheat DArT array as polymorphic among the rye accessions tested, when the quality criteria of Q >77%, call rate > 80% and discordance <0.01 were applied. These clones ('wheat rearrayed for rye array' - two 384-well plates) were included in microarrays used for genotyping in genetic diversity and mapping analyses along with the clones from rye libraries. The utility of wheat DArT markers for rye genotyping is consistent with the previous data on cross transferability of molecular markers between wheat and rye - it was proven that wheat RFLPs [3] and SSRs [4] can be successfully used for rye genome analyses.

The experiments conducted with the rye genotyping array confirmed the usefulness of the PstI/TaqI method for the analysis of the very large (ca. 8 Gbp), highly methylated and composed predominantly of repetitive sequences rye genome. It is clearly demonstrated by the proportion of clones displaying polymorphism in the materials genotyped. From the total of 1022 clones that were polymorphic in the genotyped ILs and varieties, 747 originated from the rye libraries, which corresponds to 12.2 polymorphism rate. The proportion of polymorphic clones in mapping experiments was slightly higher - 16.2 (from the total of 1965 clones that differentiated the parental lines L318 and L9 and segregated in the mapping population, 1742 originated from the 10752-clone rye library). However, those two numbers are not directly comparable as different marker selection criteria were used in diversity and mapping experiments. Moreover, the array extended with 4608 DNA clones from genetically diverse parental components of rye mapping populations was used for mapping. Those clones (mapping array 2.0) constituted 40% of the rye genotyping array 2.0 and 14.9% of them (686 markers) were polymorphic. Another 40% of the rye genotyping array 2.0 originated form DNA of parental lines of the population used in this study (mapping array 1.0) and 18.5% of them (851 markers) segregated in the RILs. The percentage of polymorphic markers observed for the diversity array during mapping experiments was 13.3 (205 markers). In the case of genetic diversity experiments the percentages of polymorphic markers were 13.1 (605 markers) and 9.2 (142 markers) for the arrays mapping 1.0 and diversity, respectively. The low number of polymorphic clones in the diversity array can be attributed to the low genetic diversity represented by the rye varieties used for library creation revealed in subsequent analyses. In general, the proportion of polymorphic markers in rye is higher than those reported for other cereals analyzed using DArT - wheat [16] and barley [15] - 9.4 and 10.4 respectively.

From 768 wheat DArT markers identified in the preliminary experiments, 275 differentiated ILs and varieties (35.8%), 223 were useful for mapping (29.0%) at P > 80%, Q>80% and a call rate of at least 90%.

Genetic Diversity Analysis

An UPGMA dendrogram based on the Jaccard's similarity matrix data obtained with 1022 DArT markers is shown in Figure 1. IL 541 separated from the remaining lines and varieties at about 38% similarity level. IL 541 is very sensitive to preharvest sprouting and originated from the cross involving the rye variety Smolickie and lines selected from the variety Kazimierskie [23]. While Kazimierskie was included in the study, to our knowledge no other IL analysed is related to Smolickie. In our opinion those are the two main factors that contributed to the observed outcome of the cluster analysis.

thumbnailFigure 1. UPGMA dendrogram based on the Jaccard's similarity matrix data obtained with 1022 DArT markers. Bootstrap support values are shown if greater than 50%.

All the varieties analyzed were clustered together at ca. 94% similarity level. Very similar results were obtained by Ma et al [24], who analyzed spring and winter rye varieties using RAPD markers and reported clustering of winter varieties at similarity level slightly lower than 90%. Rye varieties are highly heterozygous [25]. Since DArT markers were scored in the described experiments as dominant markers (presence - score '1' - vs. absence - score '0' - of a marker in genomic representation hybridized to an array), it was not possible to distinguish between markers contained in heterozygous state (one copy of a target sequence per genomic representation) and homozygous markers (two copies per genomic representation). In addition, the genotypic and allelic frequencies are likely to vary among the varieties for a number of markers, and because DArT assays measure the differences in allelic frequencies well [26], the binarisation of signal intensity data by DArTsoft was likely seriously compromised, resulting in reduced frequency of markers discriminating varieties.

Genetic similarity (GS) coefficients for all 600 possible pairs of genotypes calculated on the basis of the Jaccard's coefficient ranged from 0.26 in the pair of ILs Ot 1-3 and 541 to 0.98 in each pair of varieties Dańkowskie Nowe and Dańkowskie Złote, Dańkowskie Nowe and Tetra, Dańkowskie Nowe and Włoszanowskie, Dańkowskie Nowe and Petkus. The ILs included in the experiment were analyzed previously with respect to genetic similarity based on data from 486 PCR-based markers [27]. The result of the 2-way Mantel test used for the comparison of similarity matrices from both studies was not statistically significant. This result is not surprising since DArT markers exploit a different source of polymorphism then the PCR- based markers used earlier. A poor correlation between estimates of similarity based on data derived using different markers systems was observed also in other studies [28]. At this point, it is not possible to judge, which markers - DArT or PCR-based - reflect the genetic relationships between rye genotypes more accurately, since only a small number of genotypes was common for both studies and for some of them a detailed prodigy information is not available. Moreover, pedigree information can be an insufficient criterion to conclude about genetic relationships between rye lines obtained using selfing. Since rye varieties are highly heterozygous, it can be expected that inbred lines selected from the same variety or cross will differ in a number of loci [29]. To resolve this issue an additional study of a larger set of rye genotypes involving an analysis of phenotypic features is needed.

The average GS coefficient values were 0.47 and 0.49 for DArT and PCR-based markers, respectively. Hence, in case of this study, the discriminating power of DArT markers didn't differ significantly from that of PCR-based markers. DArT markers, however, offer substantial advantages over PCR-based markers such as independence form gel electrophoresis and sequence information, transferability, reproducibility, automated scoring and high throughput, which allow obtaining more reliable results in a more cost effective way. These advantages make DArT a method of choice over serially produced low-plex markers for diversity studies [21].

The principal correspondence plot generated from DArT data is shown in Figure 2. The first principal axis, explaining 25. 9% of the variation, differentiates rye varieties from ILs and resolves the differences within ILs, the second principal axis, explaining 11.2% of the variation, differentiates IL 541 from the remaining ILs and varieties.

thumbnailFigure 2. Principal correspondence plot based on DArT data. a. Patterns of relationships among 26 winter rye genotypes. b. Close up of the PCO plot section containing rye varieties. The first and the second principal axis explain 25.9% and 11.2% of the variation, respectively.

Chromosomal location of DArT markers

Genotyping of wheat-rye addition lines on the rye genotyping array 2.0 revealed 1872 (16.3%) DArT markers that were present in genomic representation of only one of the addition lines (while marker quality criteria Q>80% and a call rate of at least 90% were applied), and therefore their chromosomal location could be determined. Proportions of DArT markers localized on individual chromosomes were not uniform and varied from 9.4% for the chromosome 1R (175 DArT markers) to 18.0% for the chromosome 6R (337 DArT markers). The result of the Spearman rank correlation analysis of these proportions with physical lengths of rye chromosomes [30] (Figure 3) was not statistically significant. From the DArT marker with determined chromosomal location, 367 (19.6%) displayed polymorphism in parental lines of mapping population and could be used as anchor loci during map construction.

thumbnailFigure 3. Relation between DArT markers with determined chromosomal location and physical length of rye chromosomes. Relation between the proportion of DArT markers localized on individual chromosomes using wheat-rye addition lines (blue bars) and the proportion of the physical length of the rye genome represented by individual chromosomes [30] (green bars).

Genetic mapping of DArT markers

The Mendelian nature of DArT markers was demonstrated previously by Jacccoud at al. [14], Wenzl et al. [15] and Akbari et al. [16]. DArTs turned out to be highly efficient also in genetic mapping of rye. In total 1965 (17.1%) DArT markers differentiated ILs L9 and L318 and segregated in the RILs (under marker selection criteria described earlier). RECORD separated them into 43 linkage groups with marker number ranging from 2 (5 groups) to 182 (2 groups). Based on the presence of anchor loci (367 DArT and 20 SSR) the linkage groups were assembled using JoinMap 4 into 7 larger linkage groups representing rye chromosomes (Figure 4).

thumbnailFigure 4. Genetic linkage map of the cross L318xL9 containing 1818 loci. Loci with distorted segregations at p < 0.01 are indicated by asterisks and shown in red. DArT markers with chromosomal location determined using addition lines are indicated by the symbol of the respective chromosome added after underscore, their map locations are shown in blue.

The resulting map contains 1818 markers. Among them 1181 loci (65%) were unique, non co-segregating with other markers. Number of markers per chromosome ranged from 178 (115 unique) for 3R to 380 (248 unique) for 4R with an average of 260 markers (169 unique) per chromosome. The map spans 3144.6 cM, providing an average density of one unique marker every 2.7 cM (Table 1). Markers were not evenly distributed along the chromosomes - clustering of markers was apparent in certain chromosome regions, on the other hand, 4 intervals that are longer than 20 cM long can be found on the map on the chromosomes 2R (20.5 cM), 4R (21.5 cM), 6R (20.1 cM) and 7R (20.7 cM). A consensus map of rye based on data from five mapping populations was published recently by Gustafson et al. [31]. The existence of five gaps of 10 cM or longer in the map was reported: four at the telomeric ends of the chromosomes - 1R (two gaps), 4R and 5R - and one in the middle of the short arm of the chromosome 6R. The complete lack of common markers between the L318 × L9 and the consensus map makes the comparison of gap locations impossible. We expect, however, that the ongoing creation of a DArT-based integrated map of several rye mapping population, including the DS2 × RXL10 cross (used also by Gustafson et al. [31]), will deliver the information necessary to answer, whether the particular gaps are population specific, result from limitation of the marker system used, or are conserved in rye.

Table 1. Characteristics of the L318 × L9 genetic map

The length of the newly constructed L318 × L9 linkage map exceeded more than twice the length of the most saturated rye linkage map published to date, created by Bednarek et al. [5], based on analysis of a F2 population, which spanned 1386 cM and comprised 480 loci. This can be explained by map expansion due to the use of RILs, resulting from the multiple rounds of meiosis undergone [32,33]. A similar length relation between F2 and RIL-based genetic maps was observed earlier in maize [34]. As shown by Knox and Ellis [35], another factor that contributes to genetic map expansion is excess heterozygosity. An overrepresentation of heterozygots in L318 × L9 RILs was revealed in the earlier study of the population employing a larger number of SSR markers [7].

Of 1985 polymorphic markers, 729 (36.7%) deviated from the expected segregation ratio at the 1% level. Consistently with the results of previous mapping experiments conducted in rye on the same [7] and a different mapping population [8], the highest proportion of distorted unique markers was found on the chromosome 7R (62.8%). The percentage of distorted markers on remaining chromosomes varied from 24.1 (6R) to 47.8 (1R). Hackauf and Wehling [8] attributed the high percentage of distorted markers on the rye chromosome 7R maps to the existence of a locus governing zygotic selection on that chromosome. A high proportion of markers with segregation distortion observed on the chromosome 1R can be explained in a similar way as the self-incompability locus S is located there [36]. The segregation distortion is very likely to occur in a population of an out crossing crop suffering from inbreeding depression due to several cycles of selfing. Segregation distortion is a common phenomenon in rye [2,3,6,37] and in other plant species [34,38,39]. Importantly, it was shown not to influence the quality of mapping results both with simulated [40] and experimental data [34,41], even in the case of a RIL-based map with a higher proportion of distorted markers that reported in this study [34].

The DArT based L318 × L9 linkage map is the most saturated of the rye maps based solely on transferable markers available at the moment. In our opinion, the potential for its practical application is increased by the kind of the population used. RILs allow for repeated measurements and therefore for identifying loci influencing various traits of interest for breeders and researchers, but are relatively rarely used in rye genetics. The extent and uniformity of the genome coverage are still not sufficient for map based cloning and alignment of the (future) physical map. However, using the same array as used in this study for genotyping several other mapping populations, we identified over 5,000 markers segregating in at least one of these populations (data not reported). Integrating the information from several populations typed on the commercially available DArT array resulted in the construction of a high density consensus map of barley [19] and a similar strategy will be deployed to construct a 5,000+ DArT marker map of the rye genome. Having such a resource combined with a BAC anchoring strategy using DArT arrays (which was reported recently for wheat by Paux et al [42]), there is an opportunity to rapidly integrate genetic and physical mapping information for rye. Importantly, the relatively low level of marker redundancy observed in the mapping data points clearly to a potential for further array expansion using more candidate clones from additional varieties representing different germplasm pools. As some of the apparent marker redundancy in mapping data is likely to be due to genetic linkage (especially with relatively small RIL population size and non-random distribution along the chromosomes) rather than sequence redundancy, we can anticipate doubling the number of markers on the next generation array for the PstI/TaqI representation. An even larger increase in DArT marker density can be achieved by deploying other methods of complexity reduction opening opportunities for rapid identification of tightly linked markers and providing a platform for positional cloning of rye genes.


The presented study has demonstrated that the Diversity Arrays Technology can be efficiently and effectively used for the rye genome analyses - assessment of genetic similarity and linkage mapping. The 11520-clone rye genotyping panel, including 1872 DArT markers with chromosomal location determined using wheat-rye addition lines, and the DArT based genetic linkage map of the population L318 × L9 can be a valuable resource for studies on rye genome organization and in molecular breeding of the species. The array developed in the course of this study has been already applied with success to type a number of rye mapping populations.


Plant material and DNA isolation

DNA used in the experiments originated from 15 rye inbred lines (ILs): L9, L318, Ot 1-3, 541, L29, DW28, K3, L329, M2, M10, LM2020, RXL10, DS2, S76, S120, 16 rye varieties: Dańkowskie Złote, Kazimierskie, Turbo, Pancerne, Tetra Lubelskie, Dańkowskie Nowe, Włoszanowskie Nowe, Wojcieszyckie, Puławskie, Garczyńskie, Vjatka Moskovskaja, Imperial, Forrajero Klein, Harlan, Petkus, Pekka, 82 recombinant inbred lines (RILs) from the cross L318 × L9 and Chinese Spring-Imperial wheat - rye addition lines. DNA was isolated according to Murray and Thompson [43] from ca. 1 g of bulked leaf tissue of 4 weeks old greenhouse grown plants. Each IL and RIL was represented by 6 plants, each variety - by 16.

DArT Analyses

Complexity reduction

Genomic representations for microarray preparation and for genotyping were prepared using the same complexity reduction method according to the protocol described by Yang et al. [44]. Briefly, ca. 100 ng of DNA were digested with restriction enzymes PstI and TaqI (NEB). Simultaneously PstI adapter was ligated. One μl of restriction/ligation reaction was used as a template in 50 μl amplification reaction using PstI+0 primer. Sequences of adaptor and primer and cycling conditions are given in Yang et al. [44].

Microarray preparation

Preparation of Diversity Arrays Technology genomic libraries (diversity panels) and amplification of inserts from bacterial clones were carried out as described in Xia et al. [45].

A 1536-clone DArT genomic library ('diversity array'- four 384-well plates) was created from a mixture of DNA samples of 10 rye IL (without LM2020, RXL10, DS2, S76, S120) and all rye varieties included in the study. A 4608-clone library ('mapping array 1.0' - twelve 384-well plates) was prepared from a mixture of DNA of ILs L9 and L318 to increase the number of DArTs displaying polymorphism in RILs. Another 4608-clone library ('mapping array 2.0') was prepared from mixture of DNA samples of rye ILs 541, Ot1-3, LM2020, RXL10 and DS2.

Amplified inserts were precipitated, resuspended in spotting buffer DArT Spotter 2 (Wenzl et al., in preparation) and then printed on polylysine-coated slides using a MicroGrid II arrayer (Genomic Solutions, Lincoln, Neb.). For the first experimental trial, comprising genotyping of 10 rye ILs, rye varieties, and RILs from the cross L318 × L9, clones from arrays: diversity, wheat rearrayed for rye and mapping 1.0 (constituting together rye genotyping array 1.0) were printed in duplicate. For the second experimental trial, comprising genotyping of wheat-rye addition lines and RILs from the cross L318 × L9, clones from arrays: diversity, wheat rearrayed for rye, mapping 1.0 and mapping 2.0 (constituting together rye genotyping array 2.0) were printed without replication. After printing slides were put in water bath for 2 min at 95°C for DNA denaturation, then briefly immersed in 0.1 mM DTT, 0.1 mM EDTA solution and dried by centrifugation at 500 × g for 7 min.

Genotyping of DNA samples

Genomic representations of individual rye ILs, varieties, RILs and addition lines, prepared as described above were concentrated tenfold by precipitation with one volume of isopropanol, denatured for 3 min at 95°C and labeled with 1 μl 500 μM Cy3-labeled random decamers using the exo- Klenow fragment of E. coli DNA polymerase I (NEB).

Each labeled representation (target) was added to 50 μl of hybridization buffer (2 mM EDTA pH 8.0 solution containing 50 parts of ExpressHyb buffer (Clontech), 5 parts of 10 g × l-1 herring sperm DNA (Promega) and one part of the FAM-labeled polylinker fragment of the plasmid used for library preparation).

After denaturing targets were hybridized to microarrays overnight at 65°C. Subsequently slides were washed in 0.1 mM DTT, 1 × SSC, 0.1%SDS for 5 min, in 0.1 mM DTT, 1 × SSC for 5 min, in 0.1 mM DTT, 0.2 × SSC for 2 min and in 0.1 mM DTT, 0.02 × SSC for 30s. Immediately after washing slides were dried by centrifugation at 500 × g for 7 min at 30°C and then scanned with a fluorescent microarray scanner (Tecan LS300 scanner).

Batches of TIFF images pairs of individual slides were automatically analyzed by DArTsoft, a purpose-built software package developed at DArT P/L, to identify and score polymorphic markers as '1' - being present or '0' - being absent in the representation hybridized to a slide. Additionally, several parameters were computed by DArTsoft for evaluating the quality of markers, i.e. parameter Q (which measures the fraction of the of the relative target hybridization intensity as a percentage of the total variance), call rate (the percentage of DNA samples with defined '0' or '1' allele calls) and discordance (the fraction of concordant calls for replicate assays).

Testing the performance of the PstI/TaqI method of genome complexity reduction

For testing the applicability of PstI/TaqI method of complexity reduction for the DArT analysis of the rye genome, 12288 clones from a diversity panel constructed from wheat DNA [[16], webcite] were printed onto microarrays. Targets prepared from DNA samples of 10 rye ILs and 16 varieties were hybridized to those microarrays in the manner described above.

Genetic Diversity Analysis

The DArTsoft generated 0-1 scores were used to calculate a pair-wise genetic similarity matrix on the basis of the Jaccard's coefficient [46] using the NTSYS-pc, Version 2.1 [47] and then a cluster analysis was done to construct an UPGMA dendrogram. Cophenetic values were computed for the phenogram resulting in construction of cophenetic matrix. The goodness of fit of the phenogram was assessed by comparing the cophenetic matrix with the similarity matrix using the Mantel's statistic Z [48]. Bootstraping and calculation of the Jaccard's coefficient matrices from the resulting multiple datasets was performed with the PhylTools software [49]. A consensus tree was constructed with the Consense module of the PHYLIP V3.6 software package [50]. The Jaccard's similarity matrix generated from marker scores was also used for Principal Coordinate Analysis (PCO) to visualize the genetic relationships between the rye genotypes. Only clones with Q > 80% and a call rate of at least 90% were included in the analyses.

Map Construction

Marker scores obtained after hybridizations of RILs and the parental lines L9 and L318 to the rye genotyping array 1.0 and the rye genotyping array 2.0 were combined and DArT markers with P values from both experiments greater than 80% and additionally markers with P value from one experiment above 85% were chosen for map construction. The scores of polymorphic DArT markers were converted into genotype codes according to the scores of the parents. RECORD [51] was used to determine linkage groups and the order of markers within linkage groups. Linkage groups were assembled into chromosomes based on the presence of anchor loci - SSR (Table 2) and DArT. The segregations of SSR anchor loci were determined in an earlier study of the same mapping population [7]. Order of markers within chromosomes was established applying the maximum likelihood algorithm of JoinMap 4.0 [52] with the use of marker orders determined by RECORD. Graphical genotypes were inspected visually to verify the ordering of the markers. Genetic maps of the chromosomes were calculated for the established marker orders with the help of MAPMAKER/EXP 3.0 [53]. Recombination fractions were converted to cM with the Kosambi mapping function. The graphical representation of the map was drawn using MapChart software [54].

Table 2. SSR anchor loci used during map construction

Data deposition

All DArT marker scores reported in this paper were deposited in the GrainGenes database webcite and they can be retrieved using the link:;query=DArT* webcite

Competing interests

Employees of DArT P/L co-authoring this paper (Andrzej Kilian, Katarzyna Heller-Uszyńska, Peter Wenzl and Grzegorz Uszyński) may benefit financially from this work.

Authors' contributions

HBB participated in molecular analyses, genetic mapping and the design of the study, performed hierarchical clustering and ordination analysis and drafted the manuscript, KHU participated in molecular and DArTsoft analyses, PW participated in genetic mapping, GU participated in DArTsoft data analyses and downstream processing, AK participated in the design and coordination of the study, guided array development and data production and helped to draft the manuscript, MRT conceived of the study, participated in its design and coordination and helped to draft the manuscript. All authors read and approved the final manuscript.


We are grateful to Prof. A.J. Łukaszewski (University of California, Riverside, USA) for seeds of wheat-rye addition lines, to M. Niedzielski (Botanical Garden - Center for Biological Diversity Conservation of the Polish Academy of Sciences, Warsaw, Poland) for seeds of rye varieties, to Dr. L. Madej (Plant Breeding and Acclimatization Institute, Radzikow, Poland) for seeds of ILs M2 and M10 and to Dr. B. Myśków, Dr. P. Milczarski and Dr. S. Stojałowski (West Pomeranian University of Technology, Szczecin, Poland) for DNA of ILs 541, Ot 1-3, S120, S76, L2020, RXL10 and DS2. Help from our colleagues at Diversity Arrays Technology P/L, where the DArT procedures were performed is appreciated.

This research project was financially supported by the Ministry of Education and Science projects No 2P06A01929 and No NN302361833.


  1. Devos KM, Atkinson MD, Chinoy CN, Francis HA, Harcourt RL, Koebner RMD, Liu CJ, Masojæ P, Xie DX, Gale MD: Chromosomal rearrangements in the rye genome relative to that of wheat.

    Theor Appl Genet 1993, 85:673-680. Publisher Full Text OpenURL

  2. Ma XF, Wanous MK, Houchins K, Rodriguez Milla MA, Goicoechea PG, Wang Z, Xie M, Gustafson JP: Molecular linkage mapping in rye (Secale cereale L).

    Theor Appl Genet 2001, 102:517-523. Publisher Full Text OpenURL

  3. Korzun V, Malyshev S, Voylokov AV, Börner A: A genetic map of rye (Secale cereale L.) combining RFLP, isozyme, protein, microsatellite and gene loci.

    Theor Appl Genet 2001, 102:709-717. Publisher Full Text OpenURL

  4. Khlestkina EK, Than MHM, Pestsova EG, Röder MS, Malyshev SV, Korzun V, Börner A: Mapping of 99 new microsatellite-derived loci in rye (Secale cereale L.) including 39 expressed sequence tags.

    Theor Appl Genet 2004, 109:725-732. PubMed Abstract | Publisher Full Text OpenURL

  5. Bednarek PT, Masojæ P, Lewandowska R, Myśków B: Saturating rye genomic map with amplified fragment length plymorphism (AFLP) and random amplified polymorphic DNA (RAPD) markers.

    J Appl Genet 2003, 44:21-33. PubMed Abstract | Publisher Full Text OpenURL

  6. Saal B, Wricke G: Development of simple sequence repeat markers in rye (Secale cereale L.).

    Genome 1999, 42:964-972. PubMed Abstract | Publisher Full Text OpenURL

  7. Bolibok H, Gruszczyńska A, Hromada-Judycka A, Rakoczy-Trojanowska M: Identification of QTLs associated with the in vitro response of rye (Secale cereale L.).

    Cel Mol Biol Lett 2007, 12:523-535. Publisher Full Text OpenURL

  8. Hackauf B, Wehling P: Development of microsatellite markers in rye: map construction.

    Plant Breed Seed Science 2003, 48:143-151. OpenURL

  9. Hackauf B, Wehling P: Identification of microsatellite polymorphisms in expressed portion of the rye genome.

    Plant Breed 2002, 121:17-25. Publisher Full Text OpenURL

  10. Chebotar S, Röder MS, Korzun V, Saal B, Weber WE, Börner A: Molecular studies on genetic integrity of open-pollinating species rye (Secale cereale L.) after long-term genebank maintenance.

    Theor Appl Genet 2003, 107:1469-1476. PubMed Abstract | Publisher Full Text OpenURL

  11. Matos M, Camacho MV, Perez-Flores V, Pernaute B, Pinto-Carnide O, Benito C: A new aluminum tolerance gene located on rye chromosome arm 7RS.

    Theor Appl Genet 2005, 111:360-369. PubMed Abstract | Publisher Full Text OpenURL

  12. Simkova H, Safar J, Suchankova P, Kovarova P, Bartos J, Kubalakova M, Janda J, Cihalikova J, Mago R, Lelley T, Dolezel J: A novel resource for genomic of Triticeae; BAC library specific for the short arm of rye (Secale cereale L.) chromosome 1R (1RS).

    BMC Genomics 2008, 9:237. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  13. Kofler R, Bartos J, Gong L, Stift G, Suchankova P, Simkova H, Berenyi M, Burg K, Dolezel J, Lelley T: Development of microsatellite markers specific for the short arm of rye (Secale cereale L.) chromosome 1.

    Theor Appl Genet 2008, 117:915-926. PubMed Abstract | Publisher Full Text OpenURL

  14. Jaccoud D, Peng K, Feinstein D, Kilian A: Diversity arrays: a solid state technology for sequence independent genotyping.

    Nucleic Acids Res 2001, 29:e25. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  15. Wenzl P, Carling J, Kudrna D, Jaccoud D, Huttner E, Kleinhofs A, Kilian A: Diversity arrays technology (DArT) for whole-genome profiling of barley.

    Proc Natl Acad Sci USA 2004, 101:9915-9920. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  16. Akbari M, Wenzl P, Caig V, Carling J, Xia L, Yang S, Uszynski G, Mohler V, Lehmensiek A, Kuchel H, Hayden MJ, Howes N, Sharp P, Vaughan P, Rathmell B, Huttner E, Kilian A: Diversity arrays technology (DArT) for high-throughput profiling of the hexaploid wheat genome.

    Theor Appl Genet 2006, 113:1409-1420. PubMed Abstract | Publisher Full Text OpenURL

  17. Mace ES, Xia L, Jordan DR, Halloran K, Parh DK, Huttner E, Wenzl P, Kilian A: DArT markers: diversity analyses and mapping in Sorghum bicolor.

    BMC Genomics 2008, 9:26. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  18. Kilian A, Huttner E, Wenzl P, Jaccoud D, Carling J, Caig V, Evers M, Heller-Uszynska K, Cayla C, Patarapuwadol S, Xia L, Yang S, Thomson B: The fast and the cheap: SNP and DArT-based whole genome profiling for crop improvement. In Proceedings of the International Congress In the Wake of the Double Helix: From the Green Revolution to the Gene Revolution May 27-31, 2003, Bologna, Italy. Edited by Tuberosa R, Phillips RL, Gale M. Avenue Media; 2005:443-461. OpenURL

  19. Wenzl P, Li H, Carling J, Zhou M, Raman H, Paul E, Hearnden P, Maier C, Xia L, Caig V, Ovesná J, Cakir M, Poulsen D, Wang J, Raman R, Smith KP, Muehlbauer GJ, Chalmers KJ, Kleinhofs A, Huttner E, Kilian A: A high-density consensus map of barley linking DArT markers to SSR, RFLP and STS loci and phenotypic traits.

    BMC Genomics 2006, 7:206. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  20. Crossa J, Burgueno J, Dreisigacker S, Vargas M, Herrera-Foessel SA, Lillemo M, Singh RP, Trethowan R, Warburton M, Franco J, Reynolds M: Association analysis of historical bread wheat germplasm using additive genetic covariance of relatives and population structure.

    Genetics 2007, 177:1889-1913. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  21. James KE, Schneider H, Ansell SW, Evers M, Robba L, Uszynski G, Pedersen N, Newton AE, Russell SJ, Vogel JC, Kilian A: Diversity arrays technology (DArT) for pan-genomic evolutionary studies of non-model organisms.

    PLoS ONE 2008, 3(2):e1682. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  22. Bonin A, Paris M, Despres L, Tetreau G, David JP, Kilian A: A MITE-based genotyping method to reveal hundreds of DNA polymorphisms in an animal genome after a few generations of artificial selection.

    BMC Genomics 2008, 9(1):459. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  23. Masojc P, Banek-Tabor A, Milczarski P, Twardowska M: QTLs for resistance to preharvest sprouting in rye (Secale cereale L.).

    J Appl Genet 2007, 48(3):211-217. PubMed Abstract | Publisher Full Text OpenURL

  24. Ma R, Yli-Mattila T, Pulli S: Phylogenetic relationships among genotypes of worldwide collection of spring and winter ryes (Secale cereale L.) determined by RAPD-PCR markers.

    Hereditas 2004, 140:210-221. PubMed Abstract | Publisher Full Text OpenURL

  25. Wilde K, Geiger HH, Miedaner T: Significance of host complexity and diversity for race-specific leaf-rust resistance in self-fertile synthetic rye populations.

    Plant Breeding 2006, 125:225-230. Publisher Full Text OpenURL

  26. Wenzl P, Raman H, Wang J, Zhou M, Huttner E, Kilian A: A DArT platform for quantitative bulked segregant analysis.

    BMC Genomics 2007, 8:196. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  27. Bolibok H, Rakoczy-Trojanowska M, Hromada A, Pietrzykowski R: Effciency of different PCR-based marker systems in assessing genetic diversity among winter rye (Secale sereale L.) inbred lines.

    Euphytica 2005, 146:109-216. Publisher Full Text OpenURL

  28. Sigh A, Chaudhury A, Srivastava PS, Lakshmikumaran M: Comparison of AFLP and SAMPL markers for assessment of intra-population genetic variation in Azadirachta indica A. Juss.

    Plant Sci 2002, 162:17-25. Publisher Full Text OpenURL

  29. Myśków B, Masojæ P, Banek-Tabor A, Szołkowski A: Genetic diversity of inbred rye lines evaluated by RAPD analysis.

    J Appl Genet 2001, 42(1):1-14. PubMed Abstract OpenURL

  30. Schlegel R, Melz G, Nestrowicz R: A universal reference karyotype in rye Secale cerale L.

    Theor Appl Genet 1987, 74:820-826. Publisher Full Text OpenURL

  31. Gustafson JP, Ma X-F, Korzun V, Snape JW: A consensus map of rye integrating mapping data from five mapping populations.

    Theor Appl Genet 2009, 118:793-800. PubMed Abstract | Publisher Full Text OpenURL

  32. Haldane JBS, Waddington CH: Inbreeding and linkage.

    Genetics 1931, 16:357-374. PubMed Abstract | PubMed Central Full Text OpenURL

  33. Teuscher F, Guiard V, Rudolph PE, Brockmann GA: The map expansion obtained with recombinant inbred strains and intermated recombinant inbred populations for finite generation designs.

    Genetics 2005, 170:875-879. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  34. Sharopova N, McMullen MD, Schultz L, Schroeder S, Sanchez-Villeda H, Gardiner J, Bergstrom D, Houchins K, Melia-Hancock S, Musket T, Duru N, Polacco M, Edwards K, Ruff T, Register JC, Brouwer C, Thompson R, Velasco R, Chin E, Lee M, Woodman-Clikeman W, Long MJ, Liscum E, Cone K, Davis G, Coe EH Jr: Development and mapping of SSR markers for maize.

    Plant Mol Biol 2002, 48:463-481. PubMed Abstract | Publisher Full Text OpenURL

  35. Knox MR, Ellis THN: Excess heterozygosity contributes to genetic map expansion in pea recombinant inbred populations.

    Genetics 2002, 162:861-873. PubMed Abstract | PubMed Central Full Text OpenURL

  36. Hackauf B, Wehling P: Approaching the self-incompatibility locus Z in rye (Secale cereale L.) via comparative genetics.

    Theor Appl Genet 2005, 110:832-845. PubMed Abstract | Publisher Full Text OpenURL

  37. Senft P, Wricke G: An extended genetic map of rye (Secale cereale L.).

    Plant Breeding 1996, 115:508-510. Publisher Full Text OpenURL

  38. Blair MW, Pedraza F, Buendia HF, Gaitan-Solis E, Beebe SE, Gepts P, Tohme J: Development of a genome-wide anchored microsatellite map for common bean (Phaseolus vulgaris L.).

    Theor Appl Genet 2003, 107:1362-1374. PubMed Abstract | Publisher Full Text OpenURL

  39. Li JZ, Sjakste TG, Röder MS, Röder MS, Ganal MW: Development and genetic mapping of 127 new microsatelllite markers in barley.

    Theor Appl Genet 2003, 107:1021-1027. PubMed Abstract | Publisher Full Text OpenURL

  40. Hackett CA, Broadfoot LB: Effects of genotyping errors, missing values and segregation distortion in molecular marker data on the construction of linkage maps.

    Heredity 2003, 90:33-38. PubMed Abstract | Publisher Full Text OpenURL

  41. Liu S-S, Kowalsky SP, Lan T-H, Feldmann KA, Peterson AH: Genome-wide high-resolution mapping by recurrent intermating using Arabidopsis thaliana as a model.

    Genetics 1996, 142:247-258. PubMed Abstract | PubMed Central Full Text OpenURL

  42. Paux E, Sourdille P, Salse J, Saintenac C, Choulet F, Leroy P, Korol P, Michalak M, Kianian S, Spielmeyer W, Lagudah E, Somers D, Kilian A, Alaux M, Vautrin S, Bergès H, Eversole K, Appels R, Safar J, Simkova H, Dolezel J, Bernard M, Feuillet C: A physical map of the 1-Gigabase bread wheat chromosome 3B.

    Science 2008, 322:101-104. PubMed Abstract | Publisher Full Text OpenURL

  43. Murray MG, Thompson WF: Rapid isolation of high molecular weight plant DNA.

    Nucleic Acid Res 1980, 8:4321-4325. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  44. Yang S, Pang W, Ash G, Harper J, Carling J, Wenzl P, Huttner E, Zong X, Kilian A: Low level of genetic diversity in cultivated Pigeonpea compared to its wild relatives is revealed by Diversity Arrays Technology.

    Theor Appl Genet 2006, 113:585-595. PubMed Abstract | Publisher Full Text OpenURL

  45. Xia L, Peng K, Yang S, Wenzl P, Carmen de Vincente M, Fregene M, Kilian A: DArT for high-throughput genotyping of Cassava (Manihot esculenta) and its wild relatives.

    Theor Appl Genet 2005, 110:1092-1098. PubMed Abstract | Publisher Full Text OpenURL

  46. Jaccard P: Nouvelles recherches sur la distribution florale.

    Bull Soc Vaudoise Sc Nat 1908, 44:223-270. OpenURL

  47. Rohlf FJ: NTSYS-pc: numerical taxonomy and multivariate system. Version 2.1. Exeter Software, Setauket, NY; 2000.

  48. Mantel NA: The detection of disease clustering and a generalized regression approach.

    Cancer Res 1967, 27:209-220. PubMed Abstract | Publisher Full Text OpenURL

  49. Buntjer JB: Phylogenetic Computer Tools v. 1.3. Wageningen University; 1997.

  50. Felsenstein J: PHYLIP - phylogeny inference package.

    Cladistics 1989, 5:164-166. OpenURL

  51. Van Os H, Stam P, Visser RGF, van Eck HJ: RECORD: a novel method for ordering loci on a genetic linkage map.

    Theor Appl Genet 2005, 112:30-40. PubMed Abstract | Publisher Full Text OpenURL

  52. Van Ooijen JW: JoinMap® 4, Software for the calculation of genetic linkage maps in experimental populations. Edited by Kyazma BV. Wageningen, Netherlands; 2006.

  53. Lander ES, Green P, Abrahamson J, Barlow A, Daly MJ, Lincoln SE, Newberg L: MAPMAKER: an interactive computer package for constructing primary genetic linkage maps of experimental and natural populations.

    Genomics 1987, 1:174-181. PubMed Abstract | Publisher Full Text OpenURL

  54. Voorrips RE: MapChart: Software for the graphical presentation of linkage maps and QTLs.

    J Hered 2002, 93:77-78. PubMed Abstract | Publisher Full Text OpenURL