doi: 10.15389/agrobiology.2012.3.47eng

УДК 633.111.1:631.523.4:575.116:581.1:631.811.1

IDENTIFICATION AND MAPPING OF PHYSIOLOGICAL-AGRONOMIC DETERMINANTS OF SPRING SOFT WHEAT (Triticum aestivum L.) IN DOSE GRADIENT OF NITRIC NUTRITION

Yu.V. Chesnokov1, E.A. Goncharova1, N.V. Pochepnya1, M.N. Sitnikov1, N.V. Kocherina1, U. Lohwasser2, A. Borner2

In spring soft wheat (Triticum aestivum L.) with the use of ITMI mapping population, which was obtained as result of crossing of spring wheat the Opata 85 variety with W7984 synthetic hexaploid, isolated during hybridization of Triticum tauschii (Coss.) Schmahlh. [syn. Aegilops tauschii Coss., Aegilopssquarrosa auct. non L.] (accession CIGM86.940, DD) and tetraploid wheat the Altar 84 (AABB) on different backgrounds on N (without of nitrogen fertilizer, with half and overall dose of nitrogen fertilizer) the authors identified and located on chromosomes 122 QTL (quantitative trait loci). The part of phenotypic variability was established, which determined by each of revealed QTL. The authors detected the significant correlation between established loci and polymorphism on that or other determinant, which was characterized on basis of threshold value of logarithm of odds of LOD-score. The complex estimation of compared average quantities on nitrogen doses was made by variance analysis with calculation of parameters of determinant variations. The coefficients of correlation permit to determine a character adjoined correlation between determinants and gradient of doses of nitrogen nutrition. The obtained data suggests that the combining of specific QTL, revealed in soils with low content of mineral nitrogen and which became in different ecological conditions, can be used for isolation of variants with stable crop capacity.

Keywords: physiological and agronomic traits, mapping, quantitative trait loci, soft spring wheat, gradient doses of nitrogen nutrition.

 

Intensification of nitrogen nutrition of plants (especially – cereals) is quite important task aimed at raising the productivity of agricultural crops. Improved plants with efficient absorption and assimilation of nitrogen fertilizers have better potential to fair yield of grain with a high-quality protein. Today, creation of such plants is mediated by identification of relevant genes involved in these processes and subsequent introduction of them into recipient forms, or using direct methods such as marker assisted selection (MAS).
NO3- and NH4+ use to be a source of nitrogen for plant growth and development. During assimilation these ions share the same metabolic pathways under certain conditions, and they also have a number of similar features: being actively absorbed by root cells at low external concentrations; the presence of two high-affinity transport systems (HATS) for NO3- and one – for NH4+ -  established by measuring their inflow; the absorption of both ions is permanently controlled process affecting nitrogen status of a plant. Molecular studies suggest that activity of certain genes encoding NO3- -transporters is induced by NO3-, while expression of some other such genes is controlled by repression with glutamine. Along with it, there are some differences in absorption of these ions. Besides, different types of plants have peculiarities in utilization of these sources of nitrogen. In this regard, identification of genetic determinants involved in absorption and assimilation of mineral nitrogen is of particular importance for revealing their physiological mechanisms.
Mapping gene loci for quantitative traits (quantitative trait loci, QTLs) is one of major modern tools in studies of their genetic variability. One of such complex commercially valuable traits is nitrogen use efficiency (NUE). A consistency between QTLs for agronomic features and QTLs for physiological traits associated with NUE allows the development of complex economically valuable traits featured by physiological importance, ecological adaptivity and genetic value. Moreover, if gene mapping reveals QTL location coinciding with that of genes involved in nitrogen assimilation, it suggests so-called candidate genes whose allelic variants are responsible for manifestation of observed variability. In such cases, a detection of new candidate genes is possible as well. In the case of clearly identified candidate genes, a desired allele can be transferred into a recipient genotype (eg, carrier of an undesirable allele) to test it for the expected effect. QTLs for adaptation to environmental and physiological stresses, eg. drought tolerance (1, 2), phosphorus deficit (3) or insufficient nitrogen nutrition (4, 5), have been already mapped in some cereals, but only a few works on identification and mapping QTLs for NUE in wheat are available (6, 7). Though, these reports review mainly the detection of genetic determinants for genotype-environment interactions (8-10) without identification and mapping of genetic determinants of mineral nutrition.
Heterogeneous bioavailability of soil NO3- and NH4+  is one of key factors provoking violation of mineral nitrogen nutrition. In addition, seasonal and diurnal dynamics of plant growth and other physiological processes cause gradient requirements for nutrients assimilated from soil. This makes quite important to reveal all peculiarities of nitrogen nutrition in plants of different ecological zones (11). Criteria for early diagnostics of plants’ responsiveness to increasing doses of mineral fertilizers were developed in earlier authors’ studies, along with revealing the yield structure in agricultural crops under extreme environmental conditions including limited soil nutrition (12-15). The obtained data suggest that genotype’s responsiveness to mineral nutrition varies from highly reliable strong to almost unreliable one in various combinations of hydrothermal environmental factors. It was found that optimum doses of mineral fertilizers introduced in soil stimulate the activity of attractive system of plants, while excessive doses – suppress it (13). In the case of increased doses, the enhanced intake of assimilates by developing grains was revealed by calculated parameters of attractivity and micro-distribution of plastic substances in the ear (11, 14). These data show genotypic selectivity of plants in respect to soil nutrition elements, such as the gradient of mineral nitrogen supply (15, 16). Experimental simulation and impact on the system genotype-environment can be used as the effective method for rapid detection of genotypes valuable in terms of selection (11). However, a reliable estimate of all costs on target plant products and recoupment of applied soils fertilizers should consider general physiological and specific genetically determined requirements of plants.
The purpose of this study was identification and mapping of chromosomal loci responsible for the expression of important physiological and agronomic quantitative traits, whose manifestation depends on nitrogen use efficiency under different conditions of mineral soil nitrogen in the North-West region of Russia.
Technique. The object of study was spring soft wheat Triticum aestivum L. The mapped population ITMI was obtained through a cross of spring wheat cv Opata 85 with synthetic hexaploid W7984 selected at hybridization of Triticum tauschii (Coss.) Schmahlh. [Syn. Aegilops tauschii Coss., Aegilops squarrosa auct. non L.] (sample CIGM86.940, DD) and tetraploid wheat Altar 84 (AABB). The interspecific hybridization was conducted by Dr. A. Mujeeb-Kazi (Centro Internacional de Mejoramiento de Maíz y Trigo - CIMMYT, International Centre for Maize and Wheat Improvement, El Batan, Mexico) (8). 150 recombinant inbred lines were obtained through a mono-seeded selection up to F8 or F9 (Cornell University, Ithaca, USA). Seeds (10 pcs.) from 114 randomly selected abovementioned lines were kindly provided by Dr. P. Leroy (Institut National de la Rech-erche Agronomique - INRA, the National Institute of Agronomic Research, Clermont-Ferrand, France). Four lines with least viable seeds were excluded and 110 remaining lines were used in this study. Analysis of traits in plants was performed in experimental fields of VIR Pushkin branch  (Vserossiiskii Institut Rastenievodstva – VIR, N.I. Vavilov All-Russia Research and Development Institute of Plant Industry, St. Petersburg-Pushkin, Russia).
In the experimental field there were arranged three linear ditches with the length of about 120 m, width 100-105 cm and depth 40-45 cm with gaps of 70 cm. From the inside the ditches were laid with plastic film and above it filled with nitrogen-depleted soil (pH = 7,1; contents of phosphorus oxide, potassium oxide and nitric oxide – respectively, 6,5; 8,2 and 1,9 mg/100 g soil). Uniformly developed plants were obtained at watering once every 2-3 days providing 100-150 ml water per 5 kg soil. A source of mineral nitrogen, potassium and phosphorus were, respectively, NH4NO3, K2SO4 and Ca(H2PO4) x H2O introduced considering a calculated physiological norm for cereals per 1 kg dry soil: N - 0,15; P – 0,10 and K - 0,10 g of the active ingredient (according to Z.I. Zhurbitsky) (17). The experiment was set in three variants: 0N - without the introduction of nitrogen fertilizer, 1/2N and 1N –  respectively, with half and full dose of N fertilizer (1N was equal to the physiological norm of nitrogen; P and K were provided at such quantities in all variants). The fertilizers were applied during tillering using a common practice (17). In all variants there were sown 25 seeds of each line in two replicates and 20 different quantitative traits were investigated throughout the growing season. The analysis of traits was carried out by methods adopted in VIR (18) in three replications.
QTL-analysis was performed in MapMaker/QTL (19). This computer program is based on calculations of Haldane’s map function (20), so the data of gene maps from GrainGenes database (gopher: http://www.greenge-nes.cit.cornell.edu) were used to convert map distances in MapMaker/EXP 3.0 (21). The obtained data on phenotypic analysis were integrated into the base map of population ITMI (22, 23). QTLs’ location was mapped considering only the markers consistent with Kosambi’s mapping function (24). The gene map with localized QTLs was compared with the known map of chromosomes using QGENE computer program (22).
Reliability of correlations between the identified loci and polymorphism of a particular trait was evaluated upon a threshold value of the likelihood ratio of LOD-score (logarithm of odds) (25-27). Each trait was individually assessed by QTL-analysis accounting only the loci with LOD ≥ 3,0 (p < 0,001), 3,0 > LOD ≥ 2,0 (p < 0,01) and 2,0 > LOD ≥ 1,0 (p < 0,05) (21). The nature of correlations between traits and the gradient nitrogen doses was determined by a calculated correlation coefficients rxy. The ratio rxy to its error was a criterion in reliability test (t-Student test) (28). To perform a complex assessment, mean value of each trait was associated with varying doses of nitrogen using the analysis of variance with calculation of performance variation by F-Fisher's exact test and reliability of the results (29). P<0,05 was adopted as acceptable level of statistical reliability including the probability of error of 5%. Results with p <0,01 were regarded as statistically reliable, p <0,005 and p <0,001 – highly reliable. All the data were calculated in the computer program Statistica 6.0.
Results. During the experiment there were performed in total 1800 measurements of 20 quantitative traits. 122 QTLs were identified, including 43 ones with LOD ≥ 3; 65 — LOD ≥ 2 and 14 — LOD ≥ 1 (Table 1 – see descriptions of traits, location of graphic QTL peaks corresponding to a maximum LOD-score for this QTL), and 19 QTLs with LOD > 4. The loci with LOD-score of 1-2, 2-3 and over 3 were respectively determined as minor, strong and basic. LOD-score and the proportion of phenotypic variation (R2,%) provided by the identified QTLs are shown for each trait separately.
Vegetative growth traits. QTL-analysis of traits related to vegetative growth and development revealed the location of most of them on chromosomes 5A and 7B. Two such QTLs were identified on chromosome 5D. It should be noted that QTLs on 5A chromosome were inherited from the maternal form Opata 85, while QTLs for the same traits inherited from the paternal form W7984 were detected on 7B and 5D. All QTLs located on 5A were mapped in one region of the chromosome regardless of nitrogen nutrition; probably, they belong to one chromosomal locus; a similar fact was found in QTLs on 7B. The proportion of phenotypic variation determined by these loci was quite high: from 11,15% (5A) to 31,87% (7B).
Growth and development of plants was accompanied by changes in distribution of genetic determinants responsible for physiological dynamics of these processes. In early ontogeny, the traits were determined by maternal QTLs located on 5A, later the paternal loci on 7B chromosome “woke up”. Such modification of genetic determinants caused consistent changes in the degree of variation of these traits (from early to later ones), and at the phase sprouts-ripening QTLs on 5D chromosome joined in expression along with QTLs on 7B chromosome.
It is known that genetic determinants for flowering and fruiting traits expressed in higher plants at later developmental stages are often linked with M-V genes affecting growth and viability in early ontogeny. Such systems of genes present in homologous chromosomal segments of related species within the genera Gossypium, Lycopersicon, Triticum, Phaseolus, etc. (30). Several genes determining one trait can be linked to form a block, or they can be located in different chromosomes or in different arms, and their activation is controlled by one gene-coordinator (31). That’s why chromosomal loci shouldn’t be seen as only physical linkage of genes, but as a certain degree of their organic order in a group of functionally related genes, i.e. a block co-adapted genes (32). Apparently, QTLs for traits of vegetative growth and development whose manifestation is controlled by the loci mapped on 5A and 7B chromosomes belong to this type of genes.
Morphophysiological traits. QTL for these traits were dispersed over different chromosomes depending on varying conditions of nitrogen nutrition. So, all QTLs for plant height were detected on 4B and 5A chromosomes, and all of them were inherited from the paternal form. However, increasing the doses of mineral nitrogen caused activation of QTLs located on 2A, 3B and 6A chromosomes which joined the abovementioned QTLs on 4B and 5A (QTL on 2A – inherited from maternal form, 3B and 6A – from paternal one). The degree of phenotypic variability for this trait ranged from 10,27 to 19,27%.
Location of QTLs for the length of top (ear-bearing) internode and stem node size wasn’t stable as well, it was influenced by introduced doses of the fertilizer. Several interesting facts were revealed. Firstly, alleles determining the length of top (ear-bearing) internode were found to be of paternal origin, while stem node size – maternal. Secondly, the length of top (ear-bearing) internode in 1N variant was determined by the same loci with plant height in ½N and 1N variants, which shows the presence of genetic relationship between expression of these important physiological and agronomic traits.
Length and width of flag leaf are the traits important for large-scale production of wheat. QTLs for these traits were identified on different chromosomes (most of them –  on 6B and 5B), though exact location of their peaks wasn’t constant as it changed depending on introduced doses of the fertilizer. All the alleles (except the only one on 3B) were introduced by the maternal form, which suggest the possible presence of blocks of co-adapted genes determining these traits.

 

1. Characterization of major physiological and agronomical traits, their symbolic notations and designations of corresponding QTLs identified in plants from the mapped population ITMI of soft spring wheat Triticum aestivum L. in respect to the gradient of nitrogen nutrition (St. Petersburg – Pushkin)

Trait

Symbolic notation

0N

1/2N

1N

Total

location

LOD

R2

location

LOD

R2

location

LOD

R2

Duration of a period:

 

 

 

 

 

 

 

 

 

 

 

sprouts—tillering

VST

5A (74,4)a

2,57

11,15

5A (74,4)a

2,57

11,15

5A (74,4)a

2,57

11,15

0b + 3c + 0d

sprouts—booting

VSB

5A (89,0)a

3,98

16,90

5A (89,0)a
7B (341,3)

3,99
3,00

16,93
19,72

5A (89,0)a
7B (341,3)

4,24
2,90

17,89
19,13

4b + 1c + 0d

sprouts—heading

VSH

7B (341,3)
5A (89,0)a

4,84
2,55

29,78
11,18

7B (341,3)
5A (89,0)a

5,25
2,64

31,87
11,55

7B (341,3)
5A (89,0)a

5,16
2,88

31,41
12,56

3b + 3c + 0d

sprouts—flowering

VSF

7B (341,3)

3,55

22,87

7B (341,3)

4,27

26,83

7B (341,3)

3,93

24,95

3b + 0c + 0d

sprouts—ripening

VSM

7B (341,3)

4,28

28,88

7B (341,3)
5D (72,1)

4,70
2,77

29,06
18,04

7B (341,3)
5D (130,8)

4,55
2,93

28,29
12,49

3b + 2c + 0d

Plant height

PH

5A (172,1)
4B (129,3)
2A (233,2)a

3,17
2,16
2,11

19,09
13,40
12,95

5A (172,1)
3B (256,9)
4B (79,8)

2,48
2,32
2,26

11,09
11,54
10,27

6A (99,6)
4B (129,3)
5A (172,1)

3,02
2,53
2,40

19,27
11,08
10,66

2b + 7c + 0d

Length of top (ear-bearing)
internode 

StLuI

3B (256,9)
6D (160,1)a

3,50
2,17

21,69
19,53

3B (256,9)
7A (228,2)

3,97
2,85

18,95
12,90

4B (79,8)
5A (172,1)

4,64
2,94

19,95
13,04

3b + 3c + 0d

Size of stem node

StNS

6B (180,8)a
4A (237,9)a

3,36
3,05

22,11
12,87

7A (121,8)a
2D (250,3)a

2,72
1,65

17,30
11,34

2D (61,0)a
2A (198,2)a

2,97
2,48

22,36
18,48

2b + 3c + 1d

Flag leaf:

 

 

 

 

 

 

 

 

 

 

 

length

LFL

1D (42,6)a
2A (233,2)a
2D (287,8)

2,61
2,45
2,26

16,63
10,47
14,58

6B (91,3)a
1B (174,7)a
5B (154,9)a

3,14
2,71
2,56

13,47
18,21
16,37

3D (97,1)a
3B (18,1)
5B (293,0)a

1,81
1,65
1,60

9,43
13,84
12,34

1b + 5c + 3d

width

LFW

4A (248,6)a
5A (190,8)a
6B (180,8)a

2,69
2,63
2,57

18,92
20,42
17,36

5D (130,8)
6B (48,5)a
1B (169,2)a

3,48
2,77
2,21

14,81
17,57
16,62

4A (61,8)a
6B (175,2)a
2D (127,0)a

3,13
3,09
2,72

30,26
22,41
13,28

3b + 6c + 0d

Ear length

SpL

5B (96,3)
4A (258,1)a

2,79
2,20

28,10
20,12

4A (258,1)a

3,37

21,83

4A (258,1)a

2,98

20,15

1b + 3c + 0d

Number:

 

 

 

 

 

 

 

 

 

 

 

spikelets per ear

NSpt

4A (120,7)a
7A (65,8)a

3,63
2,65

21,54
28,12

7A (65,8)a

2,64

20,85

4A (120,7)a
7A (65,8)a

3,16
2,58

13,93
20,76

2b + 3c + 0d

grains per ear

NSeSp

7A (11,2)a
2D (285,7)

2,78
2,72

25,71
24,80

7B (341,3)
2B (153,8)

2,38
1,70

16,44
13,97

4D (21,2)a
4A (149,9)a

2,23
2,06

15,03
13,95

0b + 5c + 1d

Weight:

 

 

 

 

 

 

 

 

 

 

 

grains per ear

GMSp

7A (11,2)a
2D (285,7)

2,84
2,44

26,22
22,55

2B (151,0)

1,47

10,37

5A (141,8)

1,05

7,38

0b + 2c + 2d

1000 grains

TGW

5B (242,0)
6D (21,9)a
7D (164,0)

1,65
1,49
1,79

15,83
14,10
16,73

6A (101,9)
6D (15,8)a
7D (151,9)

2,22
1,60
1,81

10,00
7,23
15,63

6A (83,4)
6B (91,3)a
7D (111,1)a

2,63
2,21
3,10

11,85
16,62
20,59

1b + 3c + 5d

Wax coating:

 

 

 

 

 

 

 

 

 

 

 

on inside leaf surface

LWBi

2D (295,8)a
5D (193,0)a

3,45
2,90

21,66
19,10

2D (295,8)a
5D (105,8)

2,29
2,13

9,83
13,80

2D (300,0)a

2,56

11,33

1b + 4c + 0d

on outside leaf surface

LWBo

2D (300,0)a

4,34

26,12

2D (271,3)a
2D (300,0)a
7D (151,9)a

5,30
6,01
2,58

53,34
24,62
20,79

2D (295,8)a
7D (151,9)a
4A (120,7)

3,06
2,71
2,55

12,92
21,74
10,88

4b + 3c + 0d

on stem

StWB

2D (300,0)a
4B (31,8)a

11,72
2,36

42,36
16,10

2D (295,8)a
7D (151,9)a

13,35
2,47

45,27
20,01

2D (295,8)a
2A (230,5)

10,10
3,29

36,61
22,96

4b + 2c + 0d

on ear

SpWB

2D (300,0)a
1D (236,2)a
5D (194,3)a

6,54
3,56
2,15

26,46
15,00
14,15

2D (295,8)a
7D (151,9)a
5D (191,9)a
1D (238,4)a

6,92
2,83
2,83
2,69

26,83
22,53
18,44
11,42

2D (285,7)a
2A (230,5)
1D (238,4)a

5,30
3,53
2,62

35,83
27,77
12,28

5b + 5c + 0d

Difficulty of threshing

DifThC

2D (250,3)
7D (286,9)
2A (45,0)

2,17
1,95
1,83

20,70
18,11
17,07

2D (250,3)
6A (101,9)
5B (60,5)

3,98
3,27
2,50

26,35
14,39
16,70

2D (263,1)
7A (239,6)
1B (0,0)

4,42
3,43
2,10

21,52
29,07
14,43

4b + 3c + 2d

Total

 

 

14b + 22c + 5d

 

 

15b + 23c + 5d

 

 

14b + 20c + 4d

 

43b + 65c + 14d = 122

Note. The mapped population ITMI was derived through a cross of spring wheat cv Opata 85 with synthetic hexaploid W7984 selected at hybridization Triticum tauschii (Coss.) Schmahlh. [syn. Aegilops tauschii Coss., Aegilops squarrosa auct. non L.] (sample CIGM86.940, DD) and tetraploid wheat Altar 84 (AABB). 0N, 1/2N and 1N — variants of experiment, respectively, without the addition of mineral nitrogen fertilizer, with its half and full dose; LOD – values of logarithm of odds, R2 — proportion of phenotypic variability (%) determined by a corresponding QTL. Referent maps were presented in the report of M.W. Ganal et al. (23). In rounded brackets following a chromosome number with genome lettering is shown the distance for QTL inherited from the maternal form Opata 85 (upper index a) and parental form Synthetic (without upper index). For major QTL (b) LOD ≥ 3, for strong QTL (c) 3 > LOD ≥ 2, for minor QTL (d) 2 > LOD ≥ 1.

 

Despite low LOD-score in these cases (sometimes a little above 1,60), the proportion of phenotypic variability was almost the same as for QTLs with LOD = 3,14, which allows considering the location of these minor QTLs on chromosomes as quite reliable. Morphological and physiological traits with similar presentation were described in an earlier authors’ work on barley (15).
Distribution of QTLs for wax coating on the outer and inner side of leaves, as well as on stem and ear, wasn’t much diverse. In all variants these traits were determined primarily by major QTLs located in D-genome (except a few single QTLs on 2A, 4A and 4B chromosomes) and inherited from Opata 85. High LOD-score for wax coating trait – from 2,13 to 13,35 – confirms high stability of its manifestation. At the same time, the gradient of nitrogen supply influenced the number of loci encoding local peculiarities of wax coating.
Difficulty of thrashing was described much similar to the previous trait. Known location of these genes in D-genome was confirmed by the authors’ data. Its major QTLs were mapped on 2D chromosome, along with a few loci located on other chromosomes of A, B genomes. Gradient of nitrogen dose caused variation in number of additional QTLs for this trait, whose LOD-score  ranged from 1,83 to 4,42, and the proportion of phenotypic variation - from 17,07 to 29,07%.
Grain yield traits. Ear length trait was developed mainly under the control of major QTLs located on 4A chromosome (location peak 258,1 cM), all of which were inherited from cv Opata 85. One additional QTL for this trait was detected on 5B chromosome. The number of spikelets per spike was similarly determined by major QTLs located on 4B chromosome (location peak 120,7 cM). A number of spikelets per ear was also influenced by one additional QTL on 7A. Redistribution of these loci for both traits occurred only under the change in dose of nitrogen. QTLs for both these traits were inherited from maternal cv Opata 85, and three of them showed LOD-score over 3.
Number of grains per ear and grain weight per ear had almost similar number of identified and mapped QTLs. In the absence of nitrogen they were located on 7A and 2D chromosomes. In ½N  variant QTLs for both traits were detected on 2B, and one additional QTL for number of grains per ear was mapped on 7B. In 1N variant, location of QTLs for both traits changed – they were detected on 4A, 4D and 5A chromosomes.
Weight of 1000 grains also showed unstable location of corresponding QTLs, which though had a certain pattern of changes in the gradient of nitrogen nutrition. Thus, in all three variants there was detected QTL on 7D chromosome. At the same time, in both 0N and ½N variants one additional QTL was identified on 6D, while in 1N – on 6B.  Similarly, another QTL for weight of 1000 grains was detected on 6A chromosome in variants ½N and 1N,  but in 0N its location changed to 5B. Even when location of these loci coincided, their peaks were recorded on different regions of chromosomes. Despite a relatively low LOD-score for this trait (1,49-3,10), the proportion of phenotypic variation provided by the identified QTLs was, respectively, from 14,10 to 20,59%.
Stability of identified QTLs and reliability of their localization. Molecular genetic mapping of QTLs is one of major modern techniques allowing identification and a controlled transfer of chromosomal loci determining manifestation of economically important traits (33, 34). Specific genetic material not peculiar to a certain region can be quite valuable for breeders improving local varieties and breeding lines. The mapped population of spring wheat ITMI was derived from parental forms quite unusual and distinctive to climatic conditions of the Russian Federation (8, 35, 36), and it is of an undoubted breeding value. At the same time, only reproducible results of QTL-analysis can be useful for plant breeders. A number of QTLs described in the authors’ works satisfy it. Furthermore, the data on QTLs activity depending on environmental conditions can be used in precision agriculture providing necessary conditions for growing and best performance of desired traits.
All the identified QTLs can be divided into two groups - dependent and independent on the gradient of nitrogen dozes, i.e., on a target environmental impact.
Earlier, genetic variability for the trait of grain yield at low doses of nitrogen was revealed in experiments of other researchers who showed a significant interaction between genotype and this environmental impact (37). The absorption of nitrogen was shown as a factor providing a large share of diversity in NUE at low nitrogen content and very important for grain yield. Direct selection of wheat was found to be most effective at low doses of nitrogen fertilizer (38). The efficiency of direct selection was also noted in works performed on double haploid lines grown at low doses of nitrogen fertilizer (39) and a gradient availability of soil nitrogen (40).
Applied importance of grain yield-specific QTLs revealed at low levels of nitrogen fertilizer was shown as a promising feature useful for development of new forms stably yielding even at a low content of mineral nitrogen in soil under different environmental conditions (41). At the same time, it should be noted that in these experiments identification and analysis of QTLs was performed in field conditions typical to a certain ecological point (42), or plants were grown in pots in a greenhouse (6), which conditions were quite different from the authors’ experiments. In the latter case (6), statement of the experiment was not clear (eg., the number of plants of one line per replication grown in one vegetation pot).

 

2. Coefficients of correlation between traits and the gradient of nitrogen nutrition  in plants from the mapped population ITMI of soft spring wheat Triticum aestivum L. (St. Petersburg – Pushkin)

Symbolic
 notation

Correlation
coefficient rxy

 t-statistics

p

VST

0,00

0,00

1,000

VSB

0,07

1,17

0,241

VSH

0,04

0,72

0,473

VSF

0,09

1,52

0,128

VSM

0,63

14,11

0,000

PH

-0,03

-0,45

0,649

StLuI

-0,01

-0,03

0,974

StNS

0,15

2,50

0,013

LFL

0,19

3,12

0,002

LFW

0,15

2,40

0,017

SpL

-0,09

-1,19

0,236

NSpt

-0,06

-0,80

0,419

NSeSp

0,12

1,63

0,105

GMSp

0,07

1,00

0,317

TGM

-0,15

-2,09

0,038

LWBi

-0,07

-1,21

0,226

LWBo

-0,05

-0,87

0,386

StWB

0,01

0,19

0,846

SpWB

0,01

0,10

0,916

DifThC

-0,05

-0,64

0,523

Note. Description of traits for symbolic notations – see Table 1.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

In this study, identification and mapping of QTLs was carried out using the same set of lines of the mapped population subject to gradient doses of nitrogen in vegetation experiment performed in field conditions. Such combination of field and pot experiments provided the accuracy of pot tests preventing loss and washing off the fertilizers to adjacent soil along with exact monitoring of vegetation and development of experimental plants, while being quite close to real field conditions important for practical breeding. In most cases, changes in contents of available mineral nitrogen contributed to changes in location of the identified QTLs, or activation of additional QTLs affecting a particular trait. This effect is most likely associated with nature of the interaction genotype-environment. For example, location of QTLs for wax coating traits was stable in all variants of this experiment, while QTLs for grain yield traits changed location depending on the dose of fertilizers. This fact displays realization of adaptive potentials of plants under certain environmental conditions (30, 32, 43-45). Possibly, complex impact of environmental factors significantly defines the features of evolutionarily developed co-adapted gene blocks in each plant species, such as wheat, as well as regularities of its genetic co-adaptation system. The same basis is enabled by evolutionary and developmental "memory" of genetic systems F and R specific to each type of plant (30-32, 46, 47).
At the same time, doses of nitrogen fertilizer were reliably and positively correlated with duration of period sprouts-ripening (correlation coefficient 0,63; without the addition of nitrogen this period was shorter than in two other variants), and weakly positive though reliably correlated (rxy = 0,15-0,19) with flag leaf length and width, as well as stem nod size (Table 2). Besides, there was also a weak negative reliable correlation (rxy= – 0,15) between the weight of 1000 grains and nitrogen supply: in ½N and 1N variants the trait value was lower than in 0N variant. These conclusions are reliably confirmed by data shown in Table 3 (especially F-test and its high significance). The analysis of variance revealed that plant height was mainly influenced by the gradient doses of nitrogen, which caused a non-smooth diversity of this trait over the variants (in particular, plants grown with ½N were on average smaller than in two other variants). The length of top (ear-bearing) internode had a similar regularity, which fact indicates interrelations between these traits.

3. Results of one-factor variance analysis of manifested physiological and agronomical traits in plants from the mapped population ITMI of soft spring wheat Triticum aestivum L. in respect to the gradient of nitrogen nutrition (St. Petersburg – Pushkin)

Symbolic notation

A

B

C

p

Residual variation (error)

D

E

VST

2

0,0

0,00

1,000

306

4,3

VSB

2

8,1

0,88

0,414

306

9,2

VSH

2

9,3

0,33

0,722

306

28,5

VSF

2

29,0

1,49

0,228

306

20,0

VSM

2

7321,0

167,53

0,000

306

44,0

PH

2

597,0

5,83

0,003

201

102,0

StLuI

2

233,2

6,22

0,002

198

37,5

StNS

2

0,0200

4,90

0,008

267

0,0040

LFL

2

37,4

5,02

0,007

264

7,4

LFW

2

0,0014

2,89

0,057

264

0,0005

SpL

2

4222,7

0,98

0,375

186

4285,8

NSpt

2

127,5

0,70

0,494

186

180,4

NSeSp

2

214,3

1,57

0,211

186

136,7

GMSp

2

8,7

0,59

0,553

186

14,6

TGM

2

306,9

4,11

0,018

186

74,7

LWBi

2

3,6

1,35

0,259

306

2,4

LWBo

2

1,1

0,42

0,657

306

2,6

StWB

2

5,3

0,65

0,523

306

8,2

SpWB

2

1,5

0,19

0,827

267

8,0

DifThC

2

2,3

0,60

0,546

186

3,8

Note. See Table 2. A – number of degrees of freedom, kА; B — mean square of deviations, , sА2; C — variance ratio F = sA2/se2; D — number of degrees of freedom for residual variation, ke; E — mean square of deviations for residual variation, se2.


Another interesting finding of this research was relatively low LOD-score for some of identified QTLs (Table 1). In the literature, LOD-score below 3 is often referred to a low reliability level owing to multiple testing in QTL-analysis (21). In fact, LOD-score is the common logarithm for probability that null hypothesis is false (null hypothesis: two classes of recombinant lines carrying paternal (AA) and maternal (aa) alleles have no reliable phenotypic differences). Therefore, LOD=2 means that the hypothesis alternative to the null hypothesis is more probable in 102 times, LOD = 3 – 103 times, etc. (27). Earlier, it was reported about a comparison of c2-test and  1/A = 1/10-LOD-value where A – limit of a function of type I error (25). It was shown (48) that at high LOD-score 1/A is close to a type I error and, on the contrary, at low LOD-score the error is consistently smaller than 1/A, which suggests LOD-score as relatively conservative estimate. In this case the critical LOD-score = 3 corresponds to the maximum value of a type I error (p <0,001). At the same time, in the case of choosing a very high particular (individual) type I error, eg. 5%, high level of linkage can be reliably found at random (27). In this case, the data of various authors (8, 9) suggest that both major and minor QTLs can be located in the same positions in different experiments even in different years of research, therefore, in this study LOD-score below 3 is also considerable.
Thus, the performed study have demonstrated that activity of genetic determinants providing manifestation of physiological, morphological and agronomic quantitative traits depends on doses of introduced mineral nitrogen. The authors have revealed genomic regions involved in control of mineral nitrogen metabolism, including formation of traits related to vegetative growth and grain yield of spring wheat, which data can be used in breeding programs using identification and practical transfer of allelic variants of genetic determinants for physiological and economically important traits.

REFERENCES

1. Agrama H.A.S. and Moussa M.E., Mapping QTL in Breeding for Drought Tolerance in Maize (Zea mays L.), Euphitica, 1996, vol. 91, pp. 89-97.
2. Tuberosa R., Sanguineti M.C., Landi P., Salvi S., Casarani E. and Conti S., RFLP Mapping of Quantitative Trait Loci Controlling Abscisic Acid Concentration in Leaves of Drought-Stressed Maize (Zea mays L.), Theor. Appl. Genet., 1998, vol. 97, pp. 744-755.
3. Reiter H.S., Coors J.G., Sussman M.R. and Gabelman W.H., Genetic Analysis of Tolerance to Low Phosphorus Stress in Maize Using RFLPs, Theor. Appl. Genet., 1991, vol. 82, pp. 561-568.
4. Bertin P. and Gallais A., Physiological and Genetic Basis of Nitrogen Use Efficiency in Maize. II. QTL Detection and Coincidences, Maydica, 2001, vol. 46, pp. 53-68.
5. Hirel B., Bertin P., Quillere I. et al., Towards a Better Understanding of the Genetic and Physiological Basis of Nitrogen Use Efficiency in Maize, Plant Physiol., 2001, vol. 125, pp. 1258-1270.
6. Habash D.Z., Bernard S., Schondelmaier J., Weyen J. and Quarrie S.A., The Genetics of Nitrogen Use in Hexaploid Wheat: N Utilization, Development and Yield, Theor. Appl. Genet., 2007, vol. 114, pp. 403-419
7. Fontaine J.X., Ravel C., Pageau K., Heumez E., Dubois F., Hirel B. and Le Gouis J., A Quantitative Genetic Study for Elucidating the Contribution of Glutamine Synthetase, Glutamate Dehydrogenase and Other Nitrogen-Related Physiological Traits to the Agronomic Performance of Common Wheat, Theor. Appl. Genet., 2009, vol. 119, pp. 645-662.
8. Börner A., Schumann E., Fürste A., Göster H., Leithold B., Röder M.S. and Weber W.E., Mapping of Quantitative Trait Loci Determining Agronomic Important Characters in Hexaploid Wheat (Triticum aestivum L.), Theor. Appl. Genet., 2002, vol. 105, pp. 921-936.
9. Chesnokov Yu.V., Pochepnya N.V., Berner A., Lovasser U. Goncharova E.A. and Dragavtsev V.A., Ecological-Genetic Organization of Plant Quantitative  Traits and Mapping of the Loci Determining Agronomically Important Traits in Soft Wheat, Doklady AN, 2008, vol. 408, pp. 693-696.
10. Zheng B.S., Le Gouis J., Leflon M., Rong W.Y., Laperche A. and Brancourt-Hulmel M., Using Probe Genotypes to Dissect QTL ½ Environment Interactions for Grain Yield Components in Winter Wheat, Theor. Appl. Genet., 2010, vol. 121, pp. 1501-1517.
11. Dragavtsev V.A., Udovenko G.V., Preobrazhensky M.G., Nikulin P.L. and Stepanova N.A., Interaction Genotype-Environment at the Impact of Different Environmental Factors, Doklady RASKhN,  1994, vol. 6, pp. 3-5.
12. Udovenko G.V. and Goncharova E.A., Vliyanie ekstreml’nykh uslovii sredy na strukturu urozhaya sel’skokhozyaistvennykh rastenii (Influence of Extreme Environmental Conditions on Yield Components of Agricultural Crops), Leningrad, 1982.
13. Udovenko G.V. and Gusakova S.P., Reaction of Different Barley Cultivars to Changes in Levels of Mineral Soil Nutrition and Soil Moisture, Agrokhimiya, 1994, vol. 12, pp. 15-22.
14. Udovenko G.V. and Volkova A.M. Rannyaya diagnostika otzyvchivosti rastenii na uvelichenie dozy mineral’nukh udobrenii (Early Diagnostics of Plants’ Responsiveness to Increased Doses of Mineral Fertilizers: Methodological Guidelines), St. Petersburg, 1993.
15. Goncharova E.A., Ontogenesis of Barley and Formation of Plants’ Productivity as Response  to the Deficit of Soil Nutrition, Tr. po prikladnoi botanike, genetike i selektsii:Geneticheskie resursy rzhy, yachmenya i ovsa, St. Petersburg, 2009, issue 165, pp. 105-110.
16. Pochepnya N.V., Goncharova E.A. and Dragavtsev V.A., Genetic and Physiological Basis of Productivity of Spring Barley in Gradient Nitrogen Nutrition, in Mat. Mezhd. nauch.-prakt. konf., posvyaschennoi 45-letiyu osnovaniya Instituta genetiki i tsitologii NAN Belarusi (Papers of Sci.-Pract. Int. Congress to the 45th Anniversary of Foundation of the Institute of Genetics and Cytology of Belorussian National Academy of Sciences), Minsk, 2010, p. 71.
17. Zhurbitsky Z.I., Teoriya i praktika vegetatsionnogo metoda (Theory and Practice of Pot Trials), Moscow, 1968.
18. Filatenko A.A. and Shitova I.P., Shirokii universal’nyi klassificator SEV roda Triticum L. (COMECON Extended Unified Classification of the Genus Triticum L.), Leningrad, 1989.
19. Paterson A.H., Lander E.S., Hewitt J.D., Peterson S., Lincoln S.E. and Tanksley S.D.,Resolution of Quantitative Traits into Mendelian Factors by Using a Complete Linkage Map of Restriction Fragment Length Polymorphisms, Nature, 1988, vol. 335, pp. 721-726.
20. Haldane J.B.S., The Recombination of Linkage Values and the Calculation of Distance Between the Loci of Linkage Factors, J. Genet., 1919, vol. 8, pp. 299-309.
21. Lander E.S. and Botstein  D.,Mapping Mendelian Factors Underlying Quantitative Traits Using RFLP Linkage Maps, Genetics, 1989, vol. 129, pp. 185-199.
22. Nelson J.C., QGENE: Software for Mapping-Based Genomic Analysis and Breeding, Mol. Breed., 1997, no. 3, pp. 239-245.
23. Ganal M.W. and Röder M.S., Microsattelite and SNP Markers in Wheat Breeding. Genomics Assisted Crop Improvement: Genomics Applications in Crops, Varshney R.K. and Tuberosa R., Eds., New York-Berlin: Springer, 2007, vol. 2, pp. 1-24.
24. Kosambi D.D.,The Estimation of Map Distances from Recombination Values, Ann. Eugen., 1944, vol. 12, pp. 172-175.
25. Morton N.E., Sequential Test for the Detection of Linkage, Am. Journal Hum. Genetics, 1955, vol. 7, pp. 277-318.
26. Liu B.H., Statistical Genomics: Linkage, Mapping, and QTL Analysis, London-New York-Washington: CRC Press, 1998.
27. Kochergina N.V., Artem’eva A.M. and Chesnokov Yu.V.,  Using LOD-Score in Mapping Quantitative Trait Loci of Plants, Doklady RASKhN, 2011, vol. 3, pp. 14-17.
28. Lakin G.F., Biometriya (Biometrics), Moscow, 1990.
29. Fisher R.A., Statisticheskie metody dlya issledovatelei (Statistical Methods for Research Workers), Moscow, 1958.
30. Zhuchenko A.A., Adaptivnoe rastenievodstvo (ecologo-geneticheskie osnovy). Teoriyaipraktika. V 3 tomakh(Adaptive Plant Production (Ecological and Genetic Fundamentals). Theory and Practice), Moscow, 2008, vol. 1-3.
31. Zhuchenko A.A., Adaptivnoe rastenievodstvo (ecologo-geneticheskie osnovy)(Adaptive Plant Production (Ecological and Genetic Fundamentals), Kishinev, 1990.
32. Zhuchenko A.A., Adaptivnyi potentsial kul’turnykh rastenii: ekologo-geneticheskie osnovy (Adaptive Potential of Cultivated Crops: Ecological and Genetic Fundamentals), Kishinev, 1988.
33. Tanksley S.D.,  Mapping Polygenes, Annu. Rev. Genet., 1993, vol. 27, pp. 205-233.
34. Chesnokov Yu.V., Kartirovanie lokusov kolichestvennykh priznakov u rastenii (Mapping of Quantitative Trait Loci of Plants), St. Petersburg, 2009.
35. Nelson J.C., Van Deynze A.E., Autrique E., Sorrels M.E., Lu Y.H., Negre S., Bernard M. and Leroy P., Molecular Mapping of Wheat. Homoeologous Group 3, Genome, 1995, vol. 38, pp. 525-533.
36. Van Beem J., Mohler V., Lukman R., Van Ginkel M., William M., Crossa J. and Worland A.J., Analysis of Genetic Factors Influencing the Developmental Rate of Globally Important CIMMYT Wheat Cultivars, Crop Sci., 2005, vol. 45, pp. 2113-2119.
37. Le Gouis J., Beghin D., Heumez E. and Pluchard P., Genetic Differences for Nitrogen Uptake and Nitrogen Utilization Efficiencies in Wither Wheat, Eur. J. Argonomy, 2000, vol. 12, pp. 163-173.
38. Brancourt-Hulmel M., Heumez E., Pluchard P., Beghin D., Depatureaux C., Giraud  A. and Le Gouis J., Indirect Versus Direct Selection of Winter Wheat for Low Input or High Input Levels, Crop Sci., 2005, vol. 45, pp. 1427-1431.
39. Laperche A., Brabcourt-Hulmel M., Heumez E., Gardet O. and LeGouis J.,Estimation of Genetic Parameters of a DH Wheat Population Grown at Different N Stress Levels Characterized by Probe Genotypes, Theor. Appl. Genet., 2006, vol. 112, pp. 797-807.
40. An D., Su J., Liu Q., Zhu Y., Tong Y., Li J., Jing R. and Li Z., Mapping QTLs for Nitrogen Uptake in Relation to Early Growth of Wheat (Triticum aestivum L.), Plant and Soil, 2006, vol. 284, pp. 73-84.
41. Quarrie S.A., Steed A., Calestani C. et al., High Density Genetic Map of Hexaploid Wheat (Triticum aestivum L.) from the Cross Chinese Spring ½ SQ1 and Its Use to Compare QTLs from Grain Yield across a Range of Environments, Theor. Appl. Genet., 2005, vol. 110, pp. 865-880.
42. Laperche A., Brabcourt-Hulmel M., Heumez E., Gardet O., Hanocq E., Devienn-Barret F. and Le Gouis J., Using Genotype ½ Nitrogen Interaction Variables to Evaluate the QTL Involved in Wheat Tolerance to Nitrogen Constrains, Theor. Appl. Genet., 2007, vol. 115, pp. 399-415.
43. Worland A.J., The Influence of Flowering Time Genes on Environmental Adaptability in European Wheats, Euphytica, 1996, vol. 89, pp. 49-57.
44. Worland A.J., Börner A., Korzun V., Li W.M., Petrovic S. and Sayers E.J.,The Influence of Photoperiod Genes on the Adaptability of European Winter Wheats, Euphytica, 1998, vol. 100, pp. 385-394.
45. Iwaki K., Haruna S., Niwa T. and Kato K., Adaptation and Ecological Differentiation in Wheat with Special Reference to Geographical Variation of Growth Habit and Vrn Genotype, Plant Breed., 2001, vol. 120, pp. 107-114.
46. Zhuchenko A.A., Ekologicheskaya genetika kul’turnykh rastenii (adaptatsiya, rekombinogenez, agrobiotsenoz) (Ecological Genetics of Cultivated Crops (Adaptation, Recombinogenesis, Agrobiocenosis), Kishinev, 1980.
47. Zhuchenko A.A. and Korol’ A.B., Rekombinatsiya v evolutsii i selektsii (Recombination in Evolution and Selection), Moscow, 1985.
48. Gerber S. and Rodolphe F., Estimation and Test for Linkage between Markers: a Comparison of LOD Score and χ2 Test in a Linkage Study of Maritime Pine (Pinus pinaster Ait.), Theor. Appl. Genet., 1994, vol. 88, pp. 293-297

1N.I. Vavilov All-Russia Research and Development Institute of Plant Growing, RAAS,
St. Petersburg 190000, Russia

e-mail: yu.chesnokov@vir.nw.ru;
2Leibniz Institute of Plant Genetics and Crop
Plant Research (IPK),

Corrensstr. 3, 06466 Gatersleben, Germany

Received January 11, 2012

 

back