CC BY
19UDC 577.2:631[633.111:664.64.016.8] http://dx.doi.org/10.21498/2518-1017.3(32).2016.75973
Associations of alleles of microsatellite markers with agronomical traits of modern bread winter wheat varieties in Southern Ukraine
O. O. Kolesnyk*, O. M. Khokhlov, S. V. Chebotar
The Plant Breeding and Genetics Institute - National Center of Seed and Cultivar Investigations, 3, Ovidiopolska doroha st., Odesa, 65036, Ukraine, *e-mail: emerald-olga@ukr.net
Purpose. Defining marker-trait associations of microsatellite markers with specific regions of the genome that control important agronomical traits in the investigated varieties originated in the Plant Breeding and Genetics Institute - National Center of Seed and Cultivar Investigations and entered into the State register of plant varieties suitable for dissemination in Ukraine during different years. Methods. Molecular genetic methods (extraction of genomic DNA, polymerase chain reactions (PCR), electrophoresis of amplification products in polyacrylamide gel), field methods (phenological observations of heading date and analysis of plant height, visual assessment of the colour and length of wheat ear and awns), statistical methods (evaluation of trait means by descriptive statistic instruments of EXCEL package, ANOVA method performed by GLM instrument from AGROBASE 21 package). Results. During four growing years (2010/11, 2011/12, 2012/13, 2013/14), 47 bread winter wheat varieties were phenotypically measured and analyzed with 17 microsatellite loci. 35 marker-trait associations (MTA) for heading date, 39 for plant height, 33 for awn size, 20 for ear colour and 8 for ear size were found to be stable and significant during two-four different growing years. Conclusions. Microsatellite markers that showed substantial and stable during different growing years associations with agronomical traits can be useful and suitable for marker-assisted selection (MAS) in Ukrainian wheat breeding programs.
Keywords: marker-trait associations, Triticum aestivum L., bread winter wheat varieties.
Introduction
Bread wheat (Triticum aestivum L.) is one of the most important and widely grown food crop in many countries, and therefore the main task of breeding programs all over the world is to increase grain yield (GYLD) of wheat. Basic components of wheat yield include such breeding traits as ear number per plant, grain number per ear and grain weight [1]. Grain size and (or) weight is the target for breeding not merely as a component of grain yield, but also because of its impact on flour yield and protein content, that also affect milling and baking quality [2]. Some other agronomic traits, including anthesis date, plant height (PH) and grain volume weight, also have an effect on GYLD [3]. The efforts of breeders to increase yield are of great significance as they have to feed the expanding mass of human population taking into account the constantly changing market requirements. Such efforts focused on the selection of de-
Olga Kolesnyk
http://orcid.org/0000-0003-0603-5311 Alexander Khokhlov
http://orcid.org/0000-0003-4459-195X Sabina Chebotar
http://orcid.org/0000-002-9130-7272
sired genotypes to develop high-yielding varieties are implemented in Ukraine where wheat breeding began more than 100 years ago. For more than one century, the South of Ukraine has been regarded as the main bread winter wheat producing region of the country due to its rich black soils and a relatively dry climate. The special approach to breeding that was applied at the Plant Breeding and Genetics Institute (PBGI) [4] included traditional and non-traditional methods such as anther culture, biochemical and molecular markers, resulted in the release of new winter wheat varieties with high yield potentials and well expressed adaptation traits. Nowadays, a great progress in a conventional breeding is enhanced by in-trogression of new molecular marker technology especially microsatellite analysis that facilitates the traditional breeding process [5]. In many studies a number of QTLs for agro-nomically important traits have been identified in wheat by microsatellite analysis including plant height (PH), heading date (HD), GYLD, tiller number per plant, ear number per unit area, ear length, spikelet number, grain number, compactness of ear and grain weight [6-12].
The purpose of the research is to find marker-trait associations in the investigated varieties.
Materials and methods
The analyzed material consists of 47 bread winter wheat varieties (Triticum aestivum L.) originated at the Plant Breeding and Genetics Institute - National Center of Seed and Cultivar Investigations (PBGI - NCSCI) and entered into the State Register of plant varieties suitable for dissemination in Ukraine during different years: 'Hospodynia' (2007), 'Scarbnytsia' (2007), 'Kosovytsia' (2008), 'Anto-nivka' (2008), 'Zamozhnist' (2008), 'Blaho-darka odeska' (2009), 'Misiia odeska' (2009), 'Dalnytska' (2005), 'Yednist' (2008), 'Kiriia' (2004), 'Liona' (2005), 'Kuialnyk' (2003), 'Poshana' (2004), 'Zaporuka' (2008), 'Bunchuk'
(2009), 'Podiaka' (2008), 'Oksana' (2007), 'Zah-rava odeska' (2010), 'Epokha odeska' (2010), 'Lytanivka' (2008), 'Sluzhnytsia odeska' (2009), 'Hoduvalnytsia odeska' (2009), 'Istyna odeska'
(2010), 'Zmina' (2007), 'Dovira' (2009), 'Krasen' (2009), 'Otaman' (2008), 'Borvii' (2010), 'Tu-runchuk' (2008), 'Diuk' (2008), 'Nebokrai'
(2011), 'Khyst' (2013), 'Pylypivka' (2011), 'Zo-repad' (2011), 'Zhaivir' (2010), 'Uzhynok' (2010), 'Hurt' (2013), 'Dobrochyn' (2013), 'Vatazhok' (2011), 'Poliovyk' (2009), 'Holubka odeska' (2011), 'Kniahynia Olha' (2011), 'Lebidka odeska' (2 samples, 2011), 'Zhuravka odeska' (2011), 'Bezmezhna' (2008), 'Lastivka odeska' (2011). Such variety as 'Albatros odeskyi' (1990) and two collection samples of 'Bezosta 1' variety (1955) were taken into research as reference samples according to the recommendation of the Ukrainian Institute for Plant Variety Examination of the Ministry of Agriculture and Food of Ukraine (UIPVE of MAFU) for distinctness, uniformity and stability (DUS) testing of new varieties of bread wheat (Triti-cum aestivum L.) for the purpose of granting the Breeders' Right [13]. 'Albatros odeskyi' variety was created and recommended for breeding practices and production by Academician N. A. Lytvynenko as a source that ensures genetic productivity, frost resistance, drought tolerance, heat resistance and has excellent grain quality [14], so we used this variety as a reference crop in field experiments.
The mentioned genotypes of 47 bread winter wheat varieties were characterized using 17 SSR markers abundantly described in the literature [1, 5], the part of which was given in the guidelines on the differentiation and identification of wheat varieties [15], and additionally they were used in our previous researches [16, 17].
Genomic DNA was extracted from seedlings using modified CTAB method [18]. Polymerase
chain reactions (PCR) was performed on a Ter-tsyk thermocycler (DNA Technology, Russia) according to Röder et al. [19], with 35 instead of 45 cycles of denaturation for 1 min at 94 °C, annealing for 1 min at 50 °C (55 or 60 °C depending on the primer) and extension for 1 min at 72 °C followed by a final extension step for 5 min at 72 °C. PCR was carried out in a final reaction volume of 25 pL containing: 60 ng of DNA; 0.25 pM of each primer; 10x PCR-buffer (40 mM Tris-HCl pH 8,4; 25 mM KCl; 1.5 mM of MgCl2; 0,01% Tween-20); 0,2 mM of each dNTP; 1 unit of Taq-polymerase. The amplification products (10-pL aliquot of the PCR mixture) were separated in 7% polyacrylamide gel in 1 X TBE using Höfer scientific instruments (USA). Visualization of PCR products was performed by staining of gels in AgNO3 according to "Silver sequence TM DNA Sequencing System Technical Manual".
Image Master VDS video system (Amer-sham Pharmacia Biotech, USA) was used to assess the fragment size of the alleles at each microsatellite locus. The pUC19 DNA/ MspI and 100 bp DNA Ladder were used as standard ladders. Statistical processing of the results obtained was carried out by standard methods [20].
The seeds provided by the different research departments of PBGI - NCSCI have been sown on separate plots in field experiment conditions of 2010/11, 2011/12, 2012/13, 2013/14 growing years by the scientists of the laboratory of Variety Investigation and Breeding Process Modeling of PBGI located in Odesa, Ukraine (46° 27' 3", 30° 39' 18"). 5 seedlings of each variety were taken in a random way for microsatellite analysis (MS-analysis), among which one was used for the further growing under ear-row scheme. Each season the reference microsatellite analysis was performed to determine the purity of seeds and prevent possible contamination. The obtained grain material was analyzed for a number of important agronomic traits: heading date (HD), plant height (PH), awn size (AS), ear colour (EC), ear size (ES). HD was estimated during three growing years (2010/11, 2011/12, 2012/13) as days since the 1st of May when approximately half of ears in a plot have half emerged from the flag leaf. PH (cm) is made up of a main stem length including the length of the ear and was evaluated during four growing years (2010/11, 2011/12, 2012/13, 2013/14). Awn size was determined according to the scale of 1 up to 7 as it was recommended by [13]: 1 - absolutely without awns (as 'Bezosta 1' va-
riety), 2 - short awns, 3 - average awns, 5 -long awns, 7 - very long awns. It was estimated during three growing years (2010/11, 2011/12, 2012/13). Among the studied varieties, there was found none that would have long or very long awns. EC and ES were evaluated during two growing years (2011/12, 2012/13).
EC was determined as the degree of colour intensity with the help of the scale of 1 to 5 (1 - white, 3 - straw yellow, 5 - red). Among tested varieties we have not revealed any possessing smoky grey or black ears which were mentioned to exist according to [13]. For determining the ear size, a 9-point scale (1 -very short, 3 - short, 5 - medium, 7 - long, 9 - very long) has been used.
'Albatros odeskyi' was included in the trial as a reference variety. It was planted once every 8 studied varieties, which were sown in two rows, for more accurate registration of the whole set of features according to the method provided by [13]. Scoring indicators for ear size, awn size and ear colour were made as it was suggested by UPOV and UIEPV of MAFU [13]. Thus, 'Albatros odeskyi' variety is characterized by 2 points for AS, 3 for EC and 5 for ES.
All data were evaluated by means of descriptive statistic instruments of EXCEL package. Averaged values and their errors were calculated for every group carrying al-lele of interest, for each season separately. Significance of between-group differences was estimated by t-criterion. For illustration of differences between alternative pairs of alleles (as shown in Figures 2 and 3), along with investigation of the structure of total variation of HD and PH traits, ANOVA method was used with alleles and years as main factors; the calculations were performed by GLM instrument from AGROBASE 21 package (Agronomix Software Inc., Canada) [21].
Results and discussion
Mean values of evaluated traits in all years analyzed showed considerable levels of inter-variety diversity (Table 1). Early HD was recorded first in 'Zaporuka' variety, while 'Ok-sana' had the latest date of HD. The highest value of PH was revealed in 'Polovyk', while the smallest one 'Zaporuka' variety had. The shortest AS was in 'Zmina' variety, whereas 'Hospodynia', 'Misiia odeska', 'Kiriia', Liona', 'Kuialnyk', 'Poshana', 'Podiaka', 'Zah-rava odeska', 'Dovira' varieties had average AS. Both samples of 'Lebidka odeska' variety
had the smallest value of EC, while 'Krasen' variety showed the highest one (5, red). 'Ok-sana' variety had the smallest value of ES, whereas 'Epokha odeska' variety was characterized by the highest one among the studied varieties.
Table 1
Levels of inter-variety diversity for heading date (HD), plant height (PH), awn size (AS), ear colour (EC), ear size (ES)*
Data HD PH AS EC ES
Mean 15,36 79,57 2,64 2,21 5,13
SD 1,25 5,82 0,72 0,74 0,44
CV 8,15 7,32 27,12 33,56 8,62
min 12,33 65,00 1,00 1,00 4,00
max 18,33 95,00 3,00 5,00 6,00
* HD was estimated as days since the 1st of May; PH is given in cm; AS, EC and ES are given in points evaluated by 7-point, 5-point and 9-point scales, respectively; Mean -average values; SD - standard deviation; CV - coefficient of variation, %.
Significant marker-trait associations for HD, PH, AS, EC, ES during two-four analyzed years are denoted by asterisks (*p < 0.05 and **p < 0.01) in Tables 2, 3 (supplementary material). The highest number of associations was revealed for HD and PH. The SSR markers Xgwm186-5A and Xgwm3-3D to a great degree were associated with all analyzed traits and shown stability during three years for PH. Markers Xtaglgap-1B, Xgwm155-3A showed stable associations with HD, PH, AS during all years analyzed and significant associations with EC.
Three traits were considerably associated with three markers (Xgwm408-5B, Xgwm325-6D, Xgwm577-7B), namely Xgwm408-5B showed association for HD, PH and EC; Xgwm325-6D with HD, PH and AS; marker Xgwm577-7B with HD, EC and ES. Seven SSR markers were associated with two studied traits, whereas three markers showed associations with only one trait, namely for HD (Xgwm389-3B), AS (Xgwm190-5D), ES (Xgwm18-1B).
Table 2
The number of alleles according to their associations with heading date (HD), plant height (PH), awn size (AS), ear colour (EC), ear size (ES)
Characteristics HD* PH AS EC** ES
Number of alleles increasing
the value of a studied trait 20 17 15 10 3
Number of alleles decreasing
the value of a studied trait 15 22 17 10 5
Alternate 5 2 1 3 1
Rest 69 68 76 86 100
* Increasing values mean the earlier date of HD; ** Increasing values mean the lighter EC.
Supplementary material
Table 3
Significant associations of SSR alleles with heading date (HD), plant height (PH), awn size (AS), ear colour (EC), ear size (ES), represented as signs, revealed in field experiments during 2-4 years (the years of harvesting are shown in the header)
Locus Allele HD PH AS EC ES
2011 2012 2013 2011 2012 2013 2014 2011 2012 2013 2012 2013 2012 2013
Xgwm357 (1A) 116 i** i* i** na i* na T** T* T* T*
119
121
123 T** T* T** T* i* i* i**
125 T** T** T** T** T** T**
128 T* T* T** T* T** T**
134 T* T** T* i** i** i**
Xgwm18 (1B) 186 i*
188
192 T* T*
196 T* i* i*
Taglgap (1B) 215 T** T** T** T* i* T** T** T** i** i* T*
218 i** i** i** i* i* i* i** T** T** T** T** T* i*
235 i*
238 T** T** T** T* T* T* T** i** i** i**
265
Xgwm095 (2A) 110 na na T** i*
120 i** i* i** i* i** i* i* i*
122 T* T** T* T** i* i* T*
128 T* T* T*
130 i* i** na i** i* i*
Xgwm3 (2D) 75
77 T* T** T* i** i** i* i* i* T*
79 T* T* i** i** i** i* T*
81 i* T* T* T* T* i*
83 T* T** na T* T* T* T* i*
86 T** T* i** i* i* T* T** T* i**
Xgwm3 (2D) 88 i** i** i* i* i* T** T** T** T** T** i** i*
Xgwm155 (3A) 129 na na i*
135
137 na na i* na na na i*
139 T* i* i* T* T** T** T** T* i*
141 T* T** T* T* T** T** i* i** i** i** T*
143 i* i* i** i* T** T** T** T** i*
145 i* i* T** T** i* i* i*
147 T* i** i* i* i** i* T** T** T** i* i* T**
149 T* T* i* i* i* i** T** T** T** i** i** T**
152 i* i* na i* i** na T** T** T** i**
Xgwm389 (3B) 117 T** T* T+ i* i-
119 T* T** i* i*
134 T* T* i*
136 i** i* i+ T* i*
Xgwm389 (3B) 138 T** T*
142 T* na T* na T* T+ T*
145 i*
Xgwm165/1 (4D) 185 T*
189 na na i* i*
191 i* T* T** T** i** T*
193 i* i** i** i** T** i*
195 T* T** T** T** T*
Locus Allele HD PH AS EC ES
2011 2012 2013 2011 2012 2013 2014 2011 2012 2013 2012 2013 2012 2013
Xgwm186 (5A) 102 T** T** T** T** T**
107 T* T* na T** T* na i* i* i** i* T**
113 i** i* T** T**
115 T** i** T** i** i**
125 i** i** i* i** T* T* T** T* T** T+
129 i* i* i** i* ** *
135 i** i** i* i*
139 T** na na i* i* i* i*
Xgwm186 (5A) 142 i* T* T* T* i* i* T*
Xgwm408 (5B) 148
162 na na na i*
178 T*
185
188 i** i** i** i* i** T** T*
192 T** T** T** T* T** i* i** i* i*
Xgwm190 (5D) 204 T** T**
208 i** i**
210
Xgwm325 (6D) 115 i** i** i** na i** i* i** na T** T**
120 na i** i* na i** i** na T** T**
128 T* T** T** na T** T**
134 na na T* na na na i*
138 na i** na i* i* na T** T**
Xgwm325 (6D) 142
144 i** i* T* i* T** T** T** T**
146 i** i** i* i* T** T**
148 T** T** T* T** i** i** i** i*
150 T** T** T** na T** T** T** na i** i** i**
Xgwm577 (7B) 137 T* T* i** i** T**
152 T** T**
171 T* i* T+ T** i** i*
173 T* i** i**
175 i* i* T** i* T* T**
Xgwm437 (7D) 105
107 i* i* i* T** T** T**
109 T* T* T* i** i** i**
113 i*
Xgwm44 (7D) 176 i** T* i*
178
Xgwm44 (7D) 180 i* T* i*
183 i* T** i* i* T*
185 T** T* T** i*
187 na i** i* T*
Barc126 (7D) 138 T* na T** T**
142
146 na na i**
152 na i** na na na i**
156 T** i** i* T+ T** T** T**
158
160
162 na i** i**
164 T** i** T** i- i** i** i**
166 T** T** T**
Wmc405 (7D) 210
212
Locus Allele HD PH AS EC ES
2011 2012 2013 2011 2012 2013 2014 2011 2012 2013 2012 2013 2012 2013
Wmc405 (7D) 216 na i* T*
218 T* i** i* T* i*
220 T* i* i* i* i*
222 T* T* T* T** T* i* T*
* Significant at p < 0.05.
** Significant at p < 0.01; T (+) and i (-) symbols indicate increasing or reducing, respectively, of the studied values depending on alleles associated with these values; shading area reveals the most significant and stable marker-trait associations; na - not available.
For HD 13, marker-trait associations (MTA) were shown to be stable during three growing years, especially alleles Xgwm357123, Xgwm357128, Xtaglgap215 (Fig. 1), Xtaglgap238 (Fig. 1), Xgwm325128, Xgwm325148, Xgwm325150 and Xwmc405222 were found to be significantly associated with the later date of HD while alleles Xgwm357116, Xtaglgap218, Xgwm155
143'
Xgwm325115, Xgwm325146 were strongly associated with the earlier date of HD. Additionally, during two growing years 12 alleles were significantly associated with the later date of HD while 10 alleles showed association with the earlier date of HD. The HD values associated with presence of contrast alleles pairs in seven different loci are shown in Fig. 2.
1 2 3 4 5 6 7 8 9 10 M 11 12 13 14 15 M 16 17 1819 20
Fig. 1. Electrophoregram of DNA amplification products by the Xtaglgap-1B locus of 'Dalnytska' (1-5), 'Yednist' (6-10), 'Bezosta 1' (l1-15), 'Oksana' (16-20) varieties.
(M - molecular weight marker pUC 19/Msp I. The arrows indicate sizes in bp of fragment amplification)
The largest quantity of MTAs (39) during several growing years was revealed for PH. 17 MTA were found to be stable in three-four growing years, and 22 MTA were proved to be significant in two growing years. Overall 17 SSR markers were associated with the higher value of PH, while 22 alleles were associated with its smaller value. Fig. 3 shows PH values associated with presenting of contrast alleles pairs in seven different loci.
33 MTAs for AS were found to be stable and significant during two-three growing years. The shorter value of AS was greatly associated with alleles Xtaglgap238, Xgwm379, Xgwm155141, Xgwm325148, Xgwm437109, Xgwm357123, Xgwm357148, Xgwm095120, Xgwm095122, Xgwm3 77, Xgwm186107, Xgwm186^ Xgwm19 0 208, Xgwm325150, Barc126164, while the larger value of AS showed associations with alleles Xtaglgap218, Xgwm388, Xgwm155139, Xgwm155143, Xgwm155147,
Xgwm155149, Xgwm325144, Xgwm437107, Xgwm357116, Xgwm357125, Xgwm095128, Xgwm383, Xgwm155145, Xgwm1551^2, Xgwm186125, Xgwm19 0 204, Barc126156. Allele Xtaglgap 215 appeared to be significantly associated with both smaller and higher values of AS during three growing seasons.
20 MTAs for EC and 8 MTAs for ES were proved to be significant in two growing years. The analysis revealed that 10 alleles were associated with the darker EC and 10 alleles with the lighter EC; 5 alleles were associated with the smaller value of ES while alleles Xgwm18192, Xgwm186102, Xgwm577175 were significantly associated with the higher value of ES.
Analysis of variance (ANOVA) for HD values associated with 6 alleles in Xgwm325-6D locus and also for PH values associated with 5 alleles in Xgwm325-6D locus revealed significant affection on HD and PH, respectively, despite the fact that total variation was most-
ly attributed (88-95% of the total sum of squares) to the years.
Thus, we have detected that marker Xtagl-gap-lB showed stable association with HD, PH, AS during all years analyzed and significant associations with EC. In studies carried out by Landjeva et al. [22] it was noted that
the presence of allele Xtaglgap244 was strictly correlated with the red colour of the ear. Our investigation proved allele Xtaglgap215 to be largely associated with the darker EC while allele Xtaglgap218 to a great extent was associated with the lighter EC. In general, it was shown that the dominant gene Rg1, which
Fig. 2. HD values associated with presence of contrast alleles pairs in seven different loci, averaged values for 3 years; Y-axis shows days since the 1st of May when approximately half of ears in a plot have half emerged from the flag leaf; X-axis shows alleles - fragment pairs (bp) in seven different loci: C - Xgwm357-1A; E - Xgwm3-3D;
F - Xgwm165/1-4A; G - Xgwm155-3A; H - Xgwm389-3B; J - Xgwm325-6D; P - Xtaglgap-IB. Bars span range Mean ± (Sd/2) and whiskers - range Mean ± (LSD05/2). Solid black bars indicate pair with inter-allelic differences
to be stable over the years, but with certainty below p < 0.05.
Fig. 3. PH values associated with contrast alleles pairs in seven different loci, averaged values for 4 years; Y-axis shows PH (cm) given as a main stems length including the length of the ear; X-axis shows alleles -
fragment pairs (bp) in seven different loci: B - Xgwm095-2A; C - Xgwm357-1A; E - Xgwm3-3D; F - Xgwm165/1-4A; H - Xgwm389-3B; J - Xgwm325-6D; P - Xtaglgap-IB. Bars span range Mean ± (Sd/2) and whiskers - range Mean ± (LSD05/2). Solid black bars indicate pair with inter-allelic differences to be stable
over the years, but with certainty below p < 0.05
determines the red colour of glumes is closely linked with Gli-B1 region on chromosome IBS [23, 24]. We detected that the SSR marker Xgwm408-5B showed significant associations with HD, PH and EC. According to Leonova et al. [25], the SSR marker Xgwm408-5B appeared to be in close linkage to the vernalization response gene Vrn-B1 on 5BL. In our study the SSR marker Xgwm577-7B showed significant associations with HD, EC and ES. It was previously established that the SSR marker Xgwm577-7B is closely associated with powdery mildew resistance and the Stb8 resistance gene. Two loci on lB chromosome near Xgwm18 were found to be closely linked to a Yr (yellow rust) gene which was placed 1.9 cM distally to these SSR markers. In our study allele Xgwm18 proved to be strongly associated with the smaller value of ES while allele Xgwm18192 was considerably associated with the higher value of ES. Brbaklic et al. [5] identified that the QTL located near microsatellite marker Xg-wm18-1B explained about 23% of the total variability for PH in three growing years. The indicated marker showed association with HD and flowering date for two years. In our study, the same SSR marker did not have any association with HD or PH in tested varieties under climatic conditions of Southern Ukraine. Zanke et al. [12] reported the detection of 34 SSR alleles located on different wheat chromosomes which were associated with PH. 16 SSR alleles were associated with reducing the PH while 18 SSR alleles were associated with its enhancing. We have detected 13 SSR that were significantly associated with increasing PH while 16 SSR alleles - with its reducing. Additionally, 2 SSR alleles showed alternative effect on the value of PH. Zanke et al. [12] identified that allele Xgwm357 was significantly associa-tiated with PH and proved to be one of the most PH reducing alleles. In our study we have found that SSR allele Xgwm357128 was proved to be significantly associated with the later date of HD in three growing years under Ukrainian climatic conditions. It was pointed out by Chaturvedi et al. [26] that there is the direct positive effects of PH on grain yield, so direct identification of ge-nomic regions which affect expression of PH is important for the breeding process. In our study, 12 SSR markers associated with PH were determined on chromosomes 2A, 3A, 5A, 1B, 3D, 4D, 6D and 7D. According to our investigations we have not been analyzing reduced plant height (Rht) genes that have a
significant effect on plant height. Confirmation of marker-trait associations in our material opens new opportunities for using the mentioned microsatellite markers for MAS.
Conclusions
In order to facilitate and improve breeding process, extensive studies, focused on finding traits of interest associated with microsatellite markers and different genomic regions of wheat, are required and could be assisted by molecular marker technologies. In conclusion, 35 marker-trait associations (MTA) for heading date, 39 for plant height, 33 for awn size, 20 for ear colour and 8 for size of ears were proved to be stable and significant during two-four different growing years. 20 SSR alleles were found to be significantly associated with the later date of HD while 15 SSR alleles - with its earlier date. 5 SSR al-leles showed significant alternative effect on the value of HT. 17 SSR alleles for PH, 15 for AS and 3 for ES were found to be significantly associated with the increase of these values, respectively, while 22 SSR alleles for PH, 17 for AS and 5 for ES were associated with their reduction. As for EC, 10 SSR alleles were found to be significantly associated with the lighter EC and 10 SSR alleles - with the darker EC; 3 SSR alleles showed alternative effect on the value of EC.
Acknowledgements
We greatly appreciate V. M. Tsevma for his generous help with currying out of field experiments and for his assistance in the data processing procedures.
Bibliography
1. Röder M. S. Fine mapping of the region on wheat chromosome 7D controlling grain weight / M. S. Röder, X.-Q. Huang, A. B örner // Funct Integr Genomics. - 2008. - Vol. 8, Iss. 1. - P. 79-86. doi: 10.1007/s10142-007-0053-8
2. Gupta P. K. Genetic and molecular basis of grain size and grain number and its relevance to grain productivity in higher plants / P. K. Gupta, S. Rustgi, N. Kumar // Genome. - 2006. - Vol. 49, No. 6. - P. 565-571. doi: 10.1139/g06-063
3. High-density mapping and comparative analysis of agronomi-cally important traits on wheat chromosome 3A / M. Dilbirligi, M. Eraymana, B. T. Campbell [et al.] // Genomics. - 2006. -Vol. 88, No. 1. - P. 74-87. doi: 10.1016/j.ygeno.2006.02.001
4. Litvinenko N. A. Breeding intensive winter bread wheat varieties for Southern Ukraine / N. A. Litvinenko // Euphy-tica. - 1998. - Vol. 100, No. 1-3. - P. 7-14.
doi: 10.1023/A:1018383317939
5. Detection of QTLs for important agronomical traits in hexa-ploid wheat using association analysis / L. Brbaklic, D. Trkulja, A. Kondic-Spika [et al.] // Czech J. Genet. Plant Breed. - 2013. - Vol. 49, No. 1. - P. 1-8.
6. Mapping of quantitative trait loci determining agronomic important characters in hexaploid wheat (Triticum aestivum L.) / A. Börner, E. Schumann, A. Fürste [et al.] // Theor. Appl.
Genet. - 2002. - Vol. 105, Iss. 6-7. - P. 921-936. doi: 10.1007/s00122-002-0994-1
7. Chromosomal locations and genetic relationships of tiller and spike characters in wheat / W. L. Li, J. C. Nelson, C. Y. Chu [et al.] // Euphytica. - 2002. - Vol. 125, Iss. 3. - P. 357-366. doi: 10.1023/A:1016069809977
8. Identification and mapping of genetic loci affecting the free-threshing habit and spike compactness in wheat (Triticum aes-tivum L.) / C. Jantasuriyarat, M. I. Vales, C. J. W. Watson, 0. Riera-Lizarazu // Theor. Appl. Genet. - 2004. - Vol. 108, Iss. 2. - P. 261-273. doi: 10.1007/s00122-003-1432-8
9. Markers associated with a QTL for grain yield in wheat under drought / F. M. Kirigwi, M. V. Ginkel, G. Brown-Guedira [et al.] // Mol. Breeding. - 2007. - Vol. 20, Iss. 4. - P. 401-413. doi: 10.1007/s11032-007-9100-3
10. QTL mapping for yield and yield contributing traits in two mapping populations of bread wheat / N. Kumar, P. L. Kulwal, H. S. Balyan, P. K. Gupta // Mol. Breeding. - 2007. - Vol. 19, Iss. 2.
- P. 163-177. doi: 10.1007/s11032-006-9056-8
11. Mapping QTL for grain yield, yield components, and spike features in a doubled haploid population of bread wheat / B. Heidari, B. E. Sayed-Tabatabaei, G. Saeidi, K. Sue-naga // Genome. - 2011. - Vol. 54, Iss. 6. - P. 517-527. doi: 10.1139/g11-017
12. Genetic architecture of main effect QTL for heading date in European winter wheat (Triticum aestivum L.) [Електронний ресурс] / C. Zanke, J. Ling, J. Plieske [et al.] // PLoS ONE. - 2014.
- Vol. 9, Iss. 11. e113287. - Режим доступу : http://journals. plos.org/plosone/article?id=10.1371/journal.pone.0113287. doi: 10.1371/journal.pone.0113287
13. Методика проведення експертизи сорлв пшениц м'яко'1 (Triticum aestivum L. emend. Fiori et Paol.) на впдмпнш'сть, однор1'дтсть i стаб"|'льтсть [Електронний ресурс]. - Режим доступу : http://sops.gov.ua/uploads/files/documents/Meto-diki/12.pdf.
14. Литвиненко М. А. Теоретичт основи та методи селекц'' ози-моТ м'яко'1 пшениц на п'двищення адаптивного потенц'алу для умов Степу УкраТни : автореф. дис. ... д-ра с.-г. наук : спец. 06.01.05. «Селекц'я рослин» / М. А. Литвиненко ; 1н-т землеробства УААН. - К., 2001. - 46 с.
15. Диференц'йна здаттсть методiв iдентифiкацil сорлв пшениц' за допомогою м1'кросател1'тного анал1'зу та комп'ютер-ного визначення морфометричних параметр1'в зерна / 0. 0. Колесник, С. В. Чеботар, 0. М. Хохлов, Ю. М. Сиволап // Вкник 0НУ. Сер1'я : Бтолоп'я. - Одеса, 2009. - Т. 14, Вип. 14. - С. 27-42.
16. Диференц'аф'я сучасних сорлв озимо'' м'яко'1 пшениц' п'в-дня УкраТни за алельним складом мкросателлних локус1'в / 0. 0. Колесник, С. В. Чеботар, 0. М. Хохлов, Ю. М. Сиволап // Зб. наук. праць СГ1 - НЦНС. - 0деса : СГ1 - НЦНС, 2012. -Вип. 19. - С. 47-59.
17. Iдентифiкацiя i реестрац'я генотип'в м'яко'1 пшениф (Triticum aestivum L.), ячменю (Hordeum vulgare L.), кукурудзи (Zea mays L.), соняшника (Helianthus annuus L.) за допомогою анал1'зу мкросател^них локус1'в : метод. рекоменд. / Ю. М. Сиволап, В. В. Волкодав, М. С. Бальвинська [та 1'н.]. - 0деса : Т0В «Зовт'шрекламсервк», 2004. - 14 с.
18. Использование ПЦР-анализа в генетико-селекционных исследованиях : науч.-метод. руководство / под ред. Ю. М. Сиволапа. - К. : Аграрна наука, 1998. - С. 34-40.
19. A microsatellite map of wheat / M. S. Roder, V. Korzun, K. Wendehake [et al.] // Genetics. - 1998. - Vol. 149, No. 4. -P. 2007-2023. doi: 10.1016/B0-12-227620-5/00113-0
20. Рокицкий П. Ф. Биологическая статистика / П. Ф. Рокицкий.
- 3-е изд., испр. - М. : Колос, 1973. - 320 с.
21. Stroup W. W. Removing spatial variation from wheat yield trials: a comparison of methods / W. W. Stroup, P. S. Baenziger, D. K. Mulitze // Crop Science. - 1994. - Vol. 34, No. 1. - P. 6266. doi: 10.2135/cropsci1994.0011183X003400010011x
22. Landjeva S. Evaluation of genetic diversity among Bulgarian winter wheat (Triticum aestivum L.) varieties during the pe-
riod 1925-2003 using microsatellites / S. Landjeva, V. Korzun, G. Ganeva // Gen. Res. Crop Evol. - 2006. - Vol. 53, Iss. 8. -P. 1605-1614. doi: 10.1007/s10722-005-8718-4
23. Linkage mapping of four gene loci, Glu-B1, Gli-B1, Rg1 and Yr10 on chromosome 1B of bread wheat / P. I. Payne, L. M. Holt, R. Johnson, J. W. Snape // Genetica Agraria. - 1986. - Vol. 40. - P. 231-242.
24. Khlestkina E. K. Genes determining the coloration of different organs in wheat / E. K. Khlestkina // Russ. J of Genetics: Appl Research. - 2013. - Vol. 3, Iss. 1. - P. 54-65. doi: 10.1134/S2079059713010085
25. Mapping of the Vrn-B1 gene in Triticum aestivum L. using microsatellite markers / I. Leonova, E. Pestsova, E. Salina [et al.] // Plant Breeding. - 2003. - Vol. 122, No. 3. - P. 209-212. doi: 10.1046/j.1439-0523.2003.00818.x
26. Chaturvedi B. K. Selection parameters for some grain and quality attributes in spring wheat (Triticum aestivum L.) / B. K. Chaturvedi, R. R. Gupta // Agricultural Science Digest Kernal. -1995. - Vol. 15, No. 4. - P. 186-190.
References
1. Röder, M. S., Huang, X.-Q., & Börner, A. (2008). Fine mapping of the region on wheat chromosome 7D controlling grain weight. Funct Integr Genomics, 8(1), 79-86. doi: 10.1007/s10142-007-0053-8
2. Gupta, P. K., Rustgi, S., & Kumar, N. (2006). Genetic and molecular basis of grain size and grain number and its relevance to grain productivity in higher plants. Genome, 49(6), 565-571. doi: 10.1139/g06-063
3. Dilbirligi, M., Eraymana, M., Campbell, B. T., Randhawa, H. S., Baenziger, P. S., Dweikat, I., & Gill, K. S. (2006). High-density mapping and comparative analysis of agronomically important traits on wheat chromosome 3A. Genomics, 88(1), 74- 87. doi: 10.1016/j.ygeno.2006.02.001
4. Litvinenko N. A. (1998). Breeding intensive winter bread wheat varieties for Southern Ukraine. Euphytica, 100(1-3), 7-14. doi: 10.1023/A:1018383317939
5. Brbaklic, L., Trkulja, D., Kondic-Spika, A., Treskic, S., & Kobiljski, B. (2013). Detection of QTLs for important agronomical traits in hexaploid wheat using association analysis. Czech J. Genet. Plant Breed, 49(1), 1-8.
6. Börner, A., Schumann, E., Fürste, A., Cöster, H., Leithold, B., Röder, M., & Weber, W. (2002). Mapping of quantitative trait loci determining agronomic important characters in hexaploid wheat (Triticum aestivum L.). Theor. Appl. Genet, 105(6-7), 921936. doi: 10.1007/s00122-002-0994-1
7. Li, W. L., Nelson, J. C., Chu, C. Y., Shi, L. H., Huang, S. H., & Liu, D. J. (2002). Chromosomal locations and genetic relationships of tiller and spike characters in wheat. Euphytica, 125(3), 357-366. doi: 10.1023/A:1016069809977
8. Jantasuriyarat, C., Vales, M. I., Watson, C. J. W., & Riera-Lizara-zu, O. (2004). Identification and mapping of genetic loci affecting the free-threshing habit and spike compactness in wheat (Triticum aestivum L.). Theor. Appl. Genet, 108(2), 261273. doi: 10.1007/s00122-003-1432-8
9. Kirigwi, F. M., Ginkel, M. V., Brown-Guedira, G., Gill, B. S., Paulsen, G. M., & Fritz, A. K. (2007). Markers associated with a QTL for grain yield in wheat under drought. Mol. Breeding, 20(4), 401413. doi: 10.1007/s11032-007-9100-3
10. Kumar, N., Kulwal, P. L., Balyan, H. S., & Gupta, P. K. (2007). QTL mapping for yield and yield contributing traits in two mapping populations of bread wheat. Mol. Breeding, 19(2), 163-177. doi: 10.1007/s11032-006-9056-8
11. Heidari, B., Sayed-Tabatabaei, B. E., Saeidi, G., & Suenaga, K. (2011). Mapping QTL for grain yield, yield components, and spike features in a doubled haploid population of bread wheat. Genome, 54(6), 517-527. doi: 10.1139/g11-017
12. Zanke, C., Ling, J., Plieske, J. Kollers, S., Ebmeyer, E., Korzun, V., & Röder, M. S. (2014). Genetic architecture of main effect QTL for heading date in European winter wheat (Triticum aestivum L.). PLoS ONE, 9(11), e113287. Retrieved from http://journals.
pLos.org/pLosone/articLe?id=10.1371/journaL.pone.0113287. doi: 10.1371/journaL.pone.0113287
13. Metodyka provedennia ekspertyzy sortiv pshenytsi miakoi (Triticum aestivum L. emend. Fiori et PaoL.) na vidminnist, odnorid-nist i stabiLnist [ReguLations on the procedure and the conduct of tests for distinctness, uniformity and stability (DUS) of new bread wheat (Triticum aestivum L. emend. Fiori et PaoL.) varieties for the purpose of granting the Breeders' Right]. (n.d.). Retrieved from http://sops.gov.ua/upLoads/fiLes/documents/ Metodiki/12.pdf. [in Ukrainian]
14. Lytvynenko, M. A. (2001). Teoretychni osnovy ta metody selektsii ozymoi miakoi pshenytsi na pidvyshchennia adaptyvnoho poten-tsialu dlia umovStepu Ukrainy [TheoreticaL basis and methods of breeding bread winter wheat of universaL type for increasing its adaptive abiLity under the conditions of the Steppe zone of Ukraine]. (Extended Abstract of Dr. Agric. Sci. Diss.). Institute of AgricuLture, Kyiv, Ukraine. [in Ukrainian]
15. KoLesnik, O. O., Chebotar, S. V., KhokhLov, O. M., & SivoLap, Yu. M. (2009). Differentiating abiLity of wheat variety identification methods using microsateLLite anaLysis and computer-based determination of the grain morphometric parameters. Visnyk Odeskoho natsionalnoho universytetu. Biolohiia [Odesa NationaL University HeraLd. BioLogy], 34(14), 27-42. [in Ukrainian]
16. KoLesnyk, O. O., Chebotar, S. V., KhokhLov, O. M., & SivoLap, Yu. M. (2012). Differentiation of the modern soft winter wheat varieties from the south of Ukraine by aLLeLe composition of microsateLLite Loci. Zbirnyk naukovykh prats SHI-NTsNS [CoLLected scientific articLes of PBGI-NCSCI], 39(59), 47-59. [in Ukrainian]
17. SyvoLap, Yu. M., VoLkodav, V. V., & BaLvynska, M. S. (2004). Iden-tyfikatsiia i reiestratsiia henotypiv miakoi pshenytsi (Triticum aestivum L.), yachmeniu (Hordeum vulgare L.), kukurudzy (Zea mays L.), soniashnyka (Helianthus annuus L.) za dopomohoiu analizu mikrosatelitnykh lokusiv [Identification and registration of genotypes of common wheat (Triticum aestivum L.),
barley (Hordeum vulgare L.), maize (Zea mays L.) and sunflower (Helianthus annuus L.) by microsatellite locus analysis]. Odesa: TOV "Zovnishreklamservis". [in Ukrainian]
18. Syvolap, Yu. M. (Ed.). (1998). Ispol'zovanie PTsR-analiza v ge-netiko-selektsionnykh issledovaniyakh [The Use of PCR-analysis in Genetic and Breeding Research]. Kyiv: Ahrarna nauka. [in Russian]
19. Röder, M. S., Korzun, V., Wendehake, K., Plaschke, J., Tixier, M. H., Leroy, P., & Ganal, M. W. (1998). A microsatellite map of wheat. Genetics, 149(4), 2007-2023. doi: 10.1016/B0-12-227620-5/00113-0
20. Rokitskiy, P. F. (1973). Biologicheskaya statistika [Biological Statistics]. Moscow: Kolos. [in Russian]
21. Stroup, W. W., Baenziger, P. S., & Mulitze, D. K. (1994). Removing spatial variation from wheat yield trials: a comparison of methods. Crop Science, 34(1), 62-66. doi: 10.2135/cropsci1994.0011183X003400010011x
22. Landjeva, S., Korzun, V., & Ganeva, G. (2006). Evaluation of genetic diversity among Bulgarian winter wheat (Triticum aestivum L.) varieties during the period 1925-2003 using microsatellites. Gen. Res. Crop Evol, 53(8), 1605-1614. doi: 10.1007/s10722-005-8718-4
23. Payne, P. I., Holt, L. M., Johnson, R., & J. W. Snape, (1986). Linkage mapping of four gene loci, Glu-B1, Gli-B1, Rg1 and Yr10 on chromosome 1B of bread wheat. Genetica Agraria, 40, 231-242.
24. Khlestkina, E. K. (2013). Genes determining the coloration of different organs in wheat. Russ. J. of Genetics: Appl Research, 3(1), 54-65. doi: 10.1134/S2079059713010085
25. Leonova, I., Pestsova, E., Salina, E. Efremova T., Röder, M. S., Börner, A., & Fischbeck, G. (2003). Mapping of the Vrn-B1 gene in Triticum aestivum L. using microsatellite markers. Plant Breeding, 122(3), 209-212. doi: 10.1046/j.1439-0523.2003.00818.x
26. Chaturvedi, B. K., & Gupta, R. R. (1995). Selection parameters for some grain and quality attributes in spring wheat (Triticum aestivum L.). Agric. Sci. Digest. Kernal, 15(4), 186-190.
УДК 577.2:631[633.111:664.64.016.8]
Колесник О. 0.*, Хохлов 0. M., Чеботар С. В. АсоцтацГп алел1-в мжросател1-тних маркер1-в з агроном1-чними ознаками сучасних сорлв пшениц' м'яко'' озимо'' П'вдня Укра'ни // Сортовивчення та охорона прав на сорти рослин. - 2016. - № 3. - С. 19-29. http://dx.doi.org/10.21498/2518-1017.3(32).2016.75973
Селекщйно-генетичний iнститут - Нащональний центр нааннсзнавства та сортовивчення, вул. 0в1дюпольська дорога, 3, м. Одеса, 65036, Украина, * e-mail: emerald-olga@ukr.net
Мета. Визначення асоц1'ац1'й м1кросател1тних марке-ртв 31 специф1'чними областями геному, як' контролюють важлив1' агроном1'чн1' ознаки, на виборц1' сорлв пшени-Ц1 м'яко'1 озимо'1 селекцИ' СГ1 - НЦНС, зареестрованих в Державному реестр1' сорлв, придатних для поширення в Укра'т, протягом р1'зних рок'в. Методи. Молекулярно-генетичт методи (вид1'лення геномно'1 ДНК, пол1мераз-на ланцюгова реакц1'я (ПЛР), електрофорез продуклв ампл1ф1кацИ' в пол1'акрилам1'дному гел1'), польов1' методи (фенолопчт спостереження дати колос1'ння та анал1'з ознак «висота рослин», «довжина колосу» та «довжи-на остюк1'в», в1'зуальн1 - оц1'нка ознаки «кол1'р колосу»), статистичт методи (визначення середн1'х значень агро-ном1'чних ознак за допомогою описово'1 статистики з ви-користанням програми ExceL 1'з пакету Microsoft Office, дисперс'йний анал1'з виконано процедурою GLM (гене-
ральна л1'н1'йна модель) з1 спец1'ал1'зованого програмно-го комплексу AGROBASE 21). Результати. Протягом чо-тирьох рок1в (2010/11, 2011/12, 2012/13, 2013/14) 47 сорлв пшениц' м'яко'1 озимо'1, були вивчет за феноти-повими ознаками та проанал1зован1 за 17 мжросателгг-ними локусами. Виявлено 35 асоц1'ац1'й маркер1'в з датою колос'ння, 39 - з висотою рослин, 33 - з розм1'ром остей, 20 - з кольором колосу i 8 - з розм1'ром колосу; асо-цтацт! були стаб1льними та достов1рними протягом 2-4 рок'в достдження. Висновки. Мжросателггт маркери, як1 виявили значт та стаб1льн1 в р1зн1 роки асоц'ацИ' з агроном1чними ознаками, можуть бути корисними й при-датними для добору за допомогою маркер1в (MAS) у рамках програм селекцИ' пшениц' Укра'ни.
Ключов1 слова: асоц'ацй маркер-ознака, Triticum aestivum L., сорти пшениц м'якоï озимо)'.
УДК 577.2:631[633.111:664.64.016.8]
Колесник О. А.*, Хохлов А. Н., Чеботарь С. В. Ассоциации аллелей микросателлитных маркеров с агрономическими признаками современных сортов пшеницы мягкой озимой Юга Украины // Сортовивчення та охоро-на прав на сорти рослин. - 2016. - № 3. - С. 19-29. Ийр:/^х.сЫ.огд/10.21498/2518-1017.3(32).2016.75973
Селекционно-генетический институт - Национальный Центр семеноведения и сортоизучения, ул. Овидиопольская дорога, 3, г. Одесса, 65036, Украина, "е-таН: emerald-olga@ukr.net
Цель. Определение ассоциаций микросателлитных маркеров со специфическими областями генома, которые контролируют важные агрономические признаки, на выборке сортов пшеницы мягкой озимой селекции СГИ - НЦНС, зарегистрированных в Государственном реестре сортов растений, пригодных для выращивания в Украине, на протяжении разных лет. Методы. Моле-кулярно-генетические методы (выделение геномной ДНК, полимеразная цепная реакция (ПЦР), электрофорез продуктов амплификации в полиакриламидном геле), полевые методы (фенологические наблюдения даты колошения и анализ признаков «высота растений», «длина колоса» и «длина остей», визуальные - оценка признака «цвет колоса»), статистические методы (определение средних значений агрономических признаков с помощью описательной статистики с использованием программы Excel из пакета Microsoft Office, дисперсионный анализ выполнен процедурой GLM (генеральная
линейная модель) из специализированного программного комплекса AGROBASE 21). Результаты. В течение четырех лет (2010/11, 2011/12, 2012/13, 2013/14) 47 сортов пшеницы мягкой озимой были изучены по фе-нотипическим признакам и проанализированы с помощью 17 микросателлитных локусов. Выявлено 35 ассоциаций маркеров с датой колошения, 39 - с высотой растений, 33 - с размером остей, 20 - с цветом колоса и 8 - с размером колоса; ассоциации были стабильными и достоверными на протяжении 2-4 лет исследования. Выводы. Микросателитные маркеры, которые обнаружили значительные и стабильные в разные годы ассоциации с агрономическими признаками, могут быть полезными и пригодными для отбора с помощью маркеров (MAS) в рамках программ селекции пшеницы Украины.
Ключевые слова: ассоциации маркер-признак, Triticum aestivum L., сорта пшеницы мягкой озимой.
Hадiйшла 16.02.2016
