Genetic analysis of hybrid necrosis
Test crosses of WL711, Manitou and Pan555 (Ne2-carriers), as well as Spica and TA4152-60 (Ne1-carriers) (Table S1) showed that N9134 (or N0439) is a carrier of Ne1 allele, while ZH22 (or Z8425B) carries Ne2 allele and XN509 carries neither. Then, through systematic genetic analysis (Table S1) of winter wheat cultivars (lines) we confirmed that the hybrid necrosis in the winter wheat investigated in this study was also controlled by two complementary dominant alleles, and its genetic pattern was the same as that reported in spring wheat. The specific results are summarized as follows: The F1-generation plants derived from crosses of N9134 (or N0439)/ZH22 (or Z8425B) developed necrosis consistently; the F2 populations segregated at the ratio of 9:7(necrotic: normal) at the seedling stage; all of the BC1F1 populations showed a segregation ratio of 1:1 for necrosis compared with the normal phenotype; all plants in the F1 generations of XN509/ZH22 and XN509/N0439 exhibited the normal phenotype; and the populations from XN509//N0439/ZH22 and XN509/N0439//ZH22 crosses showed the segregation ratios of 1:3 and 1:1, respectively, for necrosis compared with the normal phenotype. The hybrid necrosis traits of plants at different growth stages and cultivation conditions are shown in Fig. 1.
The dosage effect still exists in the moderate and severe hybrid necrosis system.
The F1 plants derived from N9134 and ZH22 exhibited moderate or severe necrosis (Fig. 4c) according to the necrosis grade criteria defined by Hermsen (1963a). Through genetic analysis, we found there was one line had no living plants and only three lines exhibited all plant necrosis among the 140 F2:3 lines derived from randomly selected necrotic F2 plants (Table 1, Table S5). This means that there were no more than four double dominant homozygous lines. Even the ratio of double dominant homozygous lines to heterozygous lines was 4:136, it was not fit to, 1:8 (p = 7.83E− 04), the theoretical Mendelian ratios without dosage effect. Additionally, based on this observation combined with the theoretical (F2 = 9:7, F3 = 25:11) and practical separation (F2 = 8:7, F3 = 16:11) ratios of the F2 and F3 populations derived from necrotic plants(Table 1, Table S1), we concluded most plants with genotype Ne1Ne1Ne2Ne2 die before heading stage and are unable to produce offspring. That is, regardless of the necrosis grades, the Ne1ne1Ne2ne2-plants (double heterozygous, dosage 2) should be the weakest generally and also coincide with the corresponding F1-plants in terms of the necrotic intensity, while the plants with genotype Ne1Ne1Ne2Ne2 (double dominant homozygous, dosage 4) should always have the strongest degree of necrosis, and most will die before heading stage in the moderate and severe hybrid necrosis system.
Table 1
Genetic analysis of dosage effect associated with moderate and severe necrosis.
Years
|
Female Parent
|
Male Parent
|
Genetic population
|
Necrotic plants
|
Normal plants
|
Seeds
|
Investigation condition
|
Theoretical Mendelian ratios
(Necrotic: Normal)
|
p value
|
Theoretical Mendelian ratios with dosage effect
(Necrotic: Normal)
|
p value
|
2017-18
|
N0439
|
ZH22
|
F2
|
110
|
71
|
264
|
Seedling, Field
|
9:7
|
0.220
|
8:7
|
0.045 *
|
2012-13
|
ZH22
|
N0439
|
F2
|
137
|
132
|
380
|
Jointing, Field
|
9:7
|
0.079
|
8:7
|
0.429
|
2015-16
|
ZH22
|
N0439
|
F2
|
128
|
132
|
420
|
Heading, Field
|
9:7
|
0.023 *
|
8:7
|
0.185
|
2013-14
|
ZH22
|
N0439
|
F3
|
241
|
172
|
800 a
|
Heading, Field
|
25:11
|
9.93E− 07 **
|
16:11
|
0.708
|
2018-19 b
|
N9134
(ZH22)
|
ZH22
(N9134)
|
F3 (F2:3)
|
986
|
717
|
2637
|
Heading, Field
|
25:11
|
4.44E− 25 **
|
16:11
|
0.253
|
2018-19 b
|
N9134
(ZH22)
|
ZH22
(N9134)
|
F2:3
|
3 c
|
136 c
|
140 c
|
Heading, Field
|
1:8
|
7.83E− 04 **
|
0:8
|
——
|
Notes: a. The seeds were randomly selected from the mixed seeds of the F2 necrotic plants. b. The same population, investigated as individual plants and lines. The details of the F2:3 population are shown in Table S5. c. The number of the F2:3 lines having necrotic plants only; The number of F2:3 lines containing both necrotic and normal plants; The number of spikes selected from necrotic F2 plants. P -values obtained by chi-squared analysis of the ratios of the genetic population and theoretical Mendelian ratios. *P < 0.05; **P < 0.001. |
High-density genetic maps of the Ne1 and Ne2 in winter wheat.
Based on the linked markers of Ne1 (Chu et al. 2006) and Ne2 (Lr13 and LrLC10) (Chu et al. 2006; Qiu et al. 2020; Zhang et al. 2016) reported in spring wheat, we constructed two genetic maps of these genes in winter wheat using 269 and 264 BC1F1 plants, respectively (Fig. S2). Two high-density genetic linkage maps of these genes (Fig. 2) were then constructed mainly using the developed KASP markers and AS-PCR markers in the two BC1F1 populations consisting of 1,006 and 1,143 plants, respectively. The detailed genetic map of the Ne2 locus contained 25 molecular markers (Fig. 2a, Fig. 2b). Both the InDel marker Lseq102 (co-segregated with LrLC10) and SSR marker TC67744 co-segregated with Ne2. These markers were located at approximately 156.59 and 157.76 million bases (Mb), respectively, on chromosome arm 2BS in Chinese spring RefSeq v1.0. Furthermore, Ne2 was located between the two flanking tightly linked InDel markers, Lseq54 and Lseq22, with a genetic interval of 0.18 cM and a physical distance of about 4.45 Mb. The SNP markers Nwu_5B_4137 and Nwu_5B_5114 were identified as the two closest linked markers on either side of the Ne1 locus, with the genetic interval of 0.50 cM (0.20 and 0.30 cM). Both markers were found to be located separately at approximately 383.40 Mb and 388.01 Mb on chromosome arm 5BL in Chinese spring RefSeq v1.0, with the physical interval 4.61 Mb (Fig. 2c, Fig. 2d).
Distribution, proportion and genotype frequencies of Ne1 and Ne2 are discrete in China’s wheat regions.
Based on the planting area and yield of China’s 10 agro-ecological production zones (CHINA STATISTICAL YEARBOOK 2019, http://www.stats.gov.cn/tjsj/ndsj/2019/indexch.htm), we chose wheat cultivars (lines) for investigation according to the proportion grown in each area from the “Precise identification and innovative utilization of wheat germplasm resources” program germplasm nursery. The selected materials (consisting of 1,178 modern Chinese cultivars (MCC), 65 Chinese landraces (CL) and 121 introduced modern cultivars (IMC) ) were genotyped using molecular marker and/or hybrid test. Among them, 99 were detected by the both two methods, and 98 out of the 99 had the consistent genotyping for Ne1 or Ne2 (Table S2). In accordance with this, the two groups of genotype frequencies obtained separately using the two methods were also showed a highly positive correlation (r = 0.919, Fig. 3a). After de-redundancy, 26.2% of the cultivars (lines) were genotyped as Ne1Ne1ne2ne2, while 33.2% were genotyped as ne1ne1Ne2Ne2 (Fig. 3b). These results in Fig. 3c demonstrated discrete differences in the distribution and proportion of Ne genes from different wheat region, with highly variable correlation coefficients among the regions (Fig. 3d). Typically, Ne1Ne1ne2ne2 was significantly dominant in the region IV, whereas ne1ne1Ne2Ne2 was significantly dominant in the IMC population. Cultivars (lines) in the region III had the highest frequency of ne1ne1ne2ne2 (Fig. 3e).
Similarly, according to the CL, IMC and MCC classification (Fig. 3f), we found that the CL population had the highest frequency of both the ne1ne1ne2ne2 and Ne1Ne1ne2ne2 genotypes, while the IMC population showed the relatively lowest frequency of ne1ne1ne2ne2, but with the highest frequency of ne1ne1Ne2Ne2. Intriguingly, the Ne2 gene frequency in the MCC population was 2.3 times higher than that in the CL population due to the introduction of the IMC into wheat breeding program of China. Moreover, the proportions of the three genotypes distributed in the MCC population were highly correlated with the average of those in the CL and the IMC (r = 0.971; Fig. 3g). These findings indicated that the genotype frequencies of Ne genes in the MCC are formed by the interaction of the CL (contributing Ne1 allele) and the IMC (contributing Ne2 allele).
Ne1-nw is inherited from WE As846, while Ne2 in ZH22 is inherited from the IMC or hexaploid triticale
To clarify the origins of Ne1-nw and Ne2 in ZH22, we traced the pedigree of these genes separately. The pedigree of N9134 indicated that Ne1 is inherited from WE As846 (Fig. S1a). We found that chromosome 5B of N9134 can only be inherited from WE As846, which was consistent with a previous study of PmAS846 (Xue et al. 2012). In addition, when ZH22 was crossed with WE As846, the F1 plants showed necrosis similar to the F1-progeny of ZH22/N9134, while the offspring of ZH22/Abbondanza 5BN were normal. This result also confirmed that Ne1-nw is inherited from WE As846. However, based on these results, we could only infer that the Ne2 in ZH22 is inherited from one of the five ancestral parents, four introduced modern cultivars or one hexaploid triticale. The F1-plants obtained by crossing N9134 and Z8425B showed hybrid necrosis, while those of the other three crosses between N9134 and Zhoumai 12, Yumai 49, and Zhoumai 9 did not; therefore, we proposed that Ne2 in ZH22 was inherited from Z8425B. However, since Z8425B is derived from five ancestral cultivars (2 Italian cultivars (lines) named St2422/464 (Zhengyin 4) and St1472/506 (Zhengyin 1), 1 Russian named Прелгорная 2 (Erythrospermum-315H160/Wumang 1), 1 Mexican cultivar named Nainari60, and a hexaploid triticale line Guangmai 74) (Fig. S1b), the exact donor of Ne2 remains to be verified.