The Development and Characterization of Near-Isogenic and Pyramided Lines Carrying Resistance Genes to Brown Planthopper with the Genetic Background of Japonica Rice (Oryza sativa L.)

The brown planthopper (BPH: Nilaparvata lugens Stål.) is a major pest of rice, Oryza sativa, in Asia. Host plant resistance has tremendous potential to reduce the damage caused to rice by the planthopper. However, the effectiveness of resistance genes varies spatially and temporally according to BPH virulence. Understanding patterns in BPH virulence against resistance genes is necessary to efficiently and sustainably deploy resistant rice varieties. To survey BPH virulence patterns, seven near-isogenic lines (NILs), each with a single BPH resistance gene (BPH2-NIL, BPH3-NIL, BPH17-NIL, BPH20-NIL, BPH21-NIL, BPH32-NIL and BPH17-ptb-NIL) and fifteen pyramided lines (PYLs) carrying multiple resistance genes were developed with the genetic background of the japonica rice variety, Taichung 65 (T65), and assessed for resistance levels against two BPH populations (Hadano-66 and Koshi-2013 collected in Japan in 1966 and 2013, respectively). Many of the NILs and PYLs were resistant against the Hadano-66 population but were less effective against the Koshi-2013 population. Among PYLs, BPH20+BPH32-PYL and BPH2+BPH3+BPH17-PYL granted relatively high BPH resistance against Koshi-2013. The NILs and PYLs developed in this research will be useful to monitor BPH virulence prior to deploying resistant rice varieties and improve rice’s resistance to BPH in the context of regionally increasing levels of virulence.


Introduction
The brown planthopper (BPH: Nilaparvata lugens Stål.) is a major pest of rice (Oryza sativa L.) in tropical and subtropical Asia [1]. BPH damages rice by sucking phloem from the plants (mechanical damage) or by transmitting viruses such as rice grassy stunt virus (RGSV), rice ragged stunt phytoreovirus (RRSV) and rice wilted stunt virus (RWSV) [2,3]. At high planthopper densities, rice crops display patches of desiccated rice known as 'hopperburn.' Insecticides have been widely used to reduce BPH populations [4]. However, insecticides are damaging to human health and the environment, and are increasingly recognized as contributing to BPH outbreaks through physiological and ecological pest resurgence mechanisms [1]. Host plant resistance is considered a potentially effective alternative to harmful insecticides, that reduces BPH damage without detrimental effects on the natural enemies of BPH [5].
Monogenic resistance is vulnerable to rapid adaptation by BPH populations. Research indicates that BPH populations have sufficient genetic variability to enable them to overcome specific resistance genes when selected on a resistant host over multiple generations [4,42,43]. In the late 1970s, BPH populations adapted to varieties carrying the BPH1 and/or BPH2 genes after these were widely deployed in rice varieties across Asia [4,42,44]. A recent multi-national study has indicated that BPH populations across much of Asia have adapted to feed on rice carrying the BPH1, BPH2, BPH5, BPH7, BPH8, BPH9, BPH10 and BPH18 genes [45]. Under laboratory conditions, BPH populations continually reared for between seven to 15 generations on resistant rice varieties were capable of adapting to resistance from a range of genes, including BPH1, BPH2, BPH3, BPH8, BPH9, BPH10 and BPH32 [43,[46][47][48][49]. Through adaptation to resistance genes, BPH acquires stronger virulence against resistance genes and BPH virulence remains stable for several decades [4,50]. Therefore, it is important to preserve the effects of resistance genes by preventing BPH adaptation.
To prevent further adaptation by BPH populations to available resistance genes, a strategy for deploying resistance based on insect virulence is necessary [4]. However, BPH virulence varies under different environments depending on the predominant rice cultivars, BPH migration routes, and the length of population exposure to different resistance genes [4]. Without exposing resistance genes to BPH populations under controlled conditions prior to deployment, the potential effectiveness of the resistance genes for target regions is difficult to predict. In previous studies, the virulence of BPH has been characterized using resistant varieties [45,48,[51][52][53]. However, many BPH-resistant varieties have multiple resistance genes, such that the effects of any single resistance gene cannot be assessed using these varieties. In contrast, the effects of any single resistance gene may be revealed in detail by using near-isogenic lines (NILs) that carry the gene on the genetic background of a susceptible variety. Recently, more than 16 NILs with BPH resistance genes (BPH3, BPH4, BPH6, BPH9, BPH10, BPH12, BPH14, BPH15, BPH17, BPH18, BPH20, BPH21, BPH25, BPH26, BPH30 and BPH32) have been developed on the genetic backgrounds of several indica and japonica-susceptible varieties. These NILs have been evaluated against different BPH populations from China, the Philippines and Japan [14,[54][55][56][57][58].
Because the resistance of rice varieties carrying single genes is weaker and less durable (i.e., allowing rapid BPH adaptation) to BPH than varieties with multiple resistance genes, several researchers have proposed the pyramiding of two or more genes to enhance resistance levels and thereby avoid pest adaptation [59]. Combinations of multiple BPH resistance genes have been reported to increase levels of plant resistance to BPH. For example, a pyramided line (PYL) with BPH14 and BPH15 enhanced resistance against BPH from China compared to monogenic NILs with either BPH14 or BPH15 alone [60]. Similarly, the pyramided lines BPH6 + BPH12 PYL and BPH3 + BPH27 PYL exhibited greater resistance levels in bulk seedling tests than monogenic lines with each of the genes present alone [14,55], and BPH17 + BPH21 PYL had greater resistance against BPH in the Philippines than lines with either gene alone [57]. Pyramiding the BPH25 and BPH26 genes into a single rice line was reported to have positive epistatic effects against BPH populations collected in Vietnam, the Philippines and Japan [51,61]. Therefore, the development of rice varieties carrying multiple BPH resistance genes might be an effective way to enhance BPH resistance.
In this study, seven NILs with BPH resistance genes (BPH2, BPH3, BPH17, BPH20, BPH21, BPH32 and BPH17-ptb) and a japonica rice genetic background were developed to evaluate the effects of different resistance genes on BPH populations. Based on the NILs developed, 15 pyramided lines (PYLs) carrying two or three resistance genes were developed to enhance levels of resistance against BPH. Additionally, using the NILs and PYLs we developed, the study compared resistance against two BPH populations collected in Japan: the first was collected in 1966 (before resistant varieties were widely released) and the second was collected in 2013 (recently migrated from China to Japan). Comparisons of the reactions by BPH from each population to the NILs and PYLs indicates the utility of resistance genes and their different combinations (some with epistatic effects) against modern BPH populations.

Development of Seven NILs for BPH Resistance
Seven NILs with BPH resistance genes from three donor parents on the genetic background of Taichung 65 (T65) were developed through marker-assisted selection (MAS) and backcrossing (Tables 1  and 2). For three donor parents, IR71033-121-15 has BPH20 and BPH21; Rathu Heenati contains BPH3 and BPH17; and PTB33 carries BPH2, BPH17-ptb and BPH32 based on previous studies. For PTB33, there has been no previous report of a BPH resistance gene on chromosome 4S. However, amino acid sequences for the BPH17 locus in PTB33 were identical to those of Rathu Heenati [8]. Thus, we assume that PTB33 contains a gene for BPH resistance on chromosome 4S and tentatively named this as BPH17-ptb. The substituted chromosomal segments of the NILs were detected by polymorphic simple sequence repeat (SSR) markers that were equally distributed across the whole genome (Table 3; Figure 1). The genetic background of BPH2-NIL was analyzed using 203 polymorphic SSR markers. The ratio of substituted segments from PTB33 on BPH2-NIL was 9.1-14.8% (total 33.9-55.0 Mbp). One substituted segment with a size of 21.3-25.4 Mbp encompassing BPH2 was detected between RM247 and RM5479 on chromosome 12. The other three segments were detected between RM5426 and RM248 on chromosome 7 with a size of 3.4-4.2 Mbp, between RM5688 and RM444 on chromosome 9 with a size of 4.2-9.2 Mbp and between RM7492 and RM216 on chromosome 10 with a size of 5.0-16.2 Mbp.
The genetic background of BPH3-NIL was confirmed using 195 polymorphic SSR markers, and the ratio of substituted segment from Rathu Heenati was 1.0-3.0% (total 3.8-11.3 Mbp). One segment with a size of 1.6-1.8 Mbp including BPH3 was detected between MSSR1 and RM1369 on the short arm of chromosome 6. The other substituted segments were detected between RM1359 and RM1155 on chromosome 4 with a size of 0.5-4.2 Mbp and between RM1345 and RM3155 on chromosome 8 with a size of 1.8-4.9 Mbp.
The genetic background of BPH17-NIL was surveyed using 173 polymorphic SSR markers. The ratio of substituted segments was 1.0-4.8% (total 3.8-17.6 Mbp) containing one segment located between RM8213 and B40 on chromosome 4, including the BPH17 region.
The genetic background of BPH20-NIL was confirmed using 237 polymorphic SSR markers and the ratio of substituted segments of IR71033-121-15 was 5.6-9.6% (total 20.6-35.5 Mbp). One segment with a size of 13.8-19.9 Mb containing BPH20 was detected between RM335 and RM5900 on chromosome 4. Two other substituted segments were detected between RM224 and RM5926 on chromosome 11 (1.5-7.7 Mbp) and between RM7315 and RM3103 on chromosome 12 (5.3-7.9 Mbp).
The genetic background of BPH32-NIL was confirmed using 233 polymorphic SSR markers. The ratio of substituted segments of PTB33 on BPH32-NIL was 1.9-4.1% (total 7.1-15.1 Mbp). One segment with a size of 1.6-3.2 Mbp containing BPH32 was detected between RM6775 and RM190 on chromosome 6. Three other segments from the donor parent were detected between RM5755 and RM3280 on chromosome 3 with a size of 4.9-8.1 Mbp, between RM1306 and RM248 on chromosome 7 with a size of 0.4-3.2 Mbp and between RM5349 and RM5961 on chromosome 11 with a size of 0.2-0.6 Mbp.  Table 3. Background survey analysis of seven near-isogenic lines using SSR polymorphic markers.

Figure 2.
Damage scores (A) and fresh biomass reduction rates (B) of near-isogenic lines and pyramided lines infested with the Hadano-1966 Nilaparvata lugens population using the modified seedbox screening test at the seedling stage. The lower damage scores and fresh biomass reduction rates indicate higher resistance levels.

Entry
Average of Agronomic Trait (AVE ± SD)

Discussion
The seven NILs we developed carried BPH resistance genes on the short arm of chromosome 4 (BPH17-NIL, BPH20-NIL and BPH17-ptb-NIL), on the short arm of chromosome 6 (BPH3-NIL and BPH32-NIL) and on the long arm of chromosome 12 (BPH2-NIL and BPH21-NIL). One of the resistance genes on chromosome 12, BPH2, was originally identified from ASD7 which was used as a donor parent for many modern resistant varieties (e.g., IR36, IR42 and so on) [31,67]. BPH2 from ASD7 is identical to BPH26 in DNA sequence and resistance level [10]. BPH2 from ASD7 was resistant against the Hatano-66 population (synonym of Hadano-66) but susceptible to Nishigoshi-05, a BPH population collected in Koshi, Kumamoto Prefecture in 2005 [51]. PTB33 was reported to carry one dominant and one recessive gene [40] that were confirmed to be BPH2 and BPH3 using conventional genetic analysis [41]. However, there was no report of the exact location of BPH2 from PTB33. In our study, BPH2-NIL had similar resistance patterns to BPH2 in ASD7: BPH2-NIL was highly resistant (ADM of 68.9%) against the Hadano-66 population but less effective (ADM of 4.0%) against the recently collected population, Koshi-2013. Moreover, BPH2-NIL (PTB33) and BPH26-NIL had similar resistance levels against both Hadano-66 and Koshi-2013, suggesting that PTB33, ADR52 and ASD7 might harbor the same resistance gene. Further sequence analysis for BPH2 from PTB33 is necessary to understand its genetic basis. Another gene on chromosome 12, BPH21, was originally identified from IR71033-121-15, an introgression line derived from O. minuta and estimated to be located between two markers, S12094A and B122, on the long arm of chromosome 12 [17]. Recently, BPH21 has been reported to be allelic to BPH26 [6] and BPH18 [68] based on amino acid sequences. Both BPH18 and BPH26 were isolated and located at 22.9 Mbp on chromosome 12 [9,10]. Therefore, we estimated that the location of BPH21 was around 22.9 Mbp on chromosome 12, and the region carrying BPH21 from IR71033-121-15 was selected using RM1246 (19.2 Mbp) and RM28493 (23.3 Mbp) in this study.
The BPH17 locus on chromosome 4S from Rathu Heenati has been reported by Sun et al. (2005) [16]. BPH17 was mapped between two markers, RHD9 (6.2 Mbp) and RHC10 (7.0 Mbp), on chromosome 4S and isolated by Liu et al. (2014) [8]. The amino acid sequence and chromosomal location of BPH17 from Rathu Heenati were the same as those of BPH17-ptb from PTB33 [8]. In this study, resistance of BPH17-NIL and BPH17-ptb-NIL against the Hadano-66 population differed; however, both NILs had similar effects on the Koshi-2013 population. The different resistant levels might be because the loci were derived from different accessions or varieties of rice. Therefore, the amino acid sequences of PTB33 and Rathu Heenati used in this study on the BPH17 locus should be determined for future research. Additionally, BPH20 was detected between two markers, B42 (8.7 Mbp) and B44 (8.9 Mbp) on chromosome 4 [17]. Two NILs for BPH17 and BPH20 on the genetic background of 9311 varieties developed by   [51] showed different resistance levels against a BPH population from China [68]. In our study, the resistance levels of BPH17 and BPH20 were different in both MSST and antibiosis tests against the Hadano-66 population and against the Koshi-2013 population, which corresponds well with previous research by   [51]. Therefore, the genes on chromosome 4S of IR71033-121-15, PTB33 and Rathu Heenati might be different. To confirm this, further sequence analyses are needed for the three loci BPH17, BPH17-ptb and BPH20.
Among six genes/QTLs that have been identified on the short arm of chromosome 6 of O. sativa and its wild relatives [3,12], BPH3 and BPH32 have been widely introduced to elite rice cultivars to improve BPH resistance and were related to durable and broad-spectrum resistance in PTB33 and Rathu Heenati [67,69]. In previous research, BPH3 was mapped onto chromosome 6 between two markers, RM19291 (1.2 Mbp) and RM8072 (1.4 Mbp) [69]. BPH32 from PTB33 was identified at the same location as BPH3 from Rathu Heenati, but the amino acid sequence of BPH3 was not identical to that of BPH32 [12]. In our study, the resistance levels of the BPH3-NIL were slightly different from those of the BPH32-NIL, suggesting that BPH3 might be different from BPH32. A comparison of amino acid sequences between BPH3 and BPH32 would be necessary to confirm whether these resistance genes are different.
Among the developed NILs, the BPH3-NIL, BPH17-NIL, BPH17-ptb-NIL and BPH32-NIL had around 97.0% of their chromosomal segments from the recurrent parent. This proportion coincides with the theoretical ratio for substituting chromosomal segments from recurrent parents by backcrossing four times. The other NILs had fewer chromosomal segments from T65 than the theoretical rate. The substituted chromosomal segments from the donor parents might be related to undesirable traits such as the suppression of the associated BPH resistance gene. Additionally, due to the low density of available polymorphic SSR markers between T65 and donor DNA around the target genes, the BPH resistance genes on the NILs were selected using two flanking markers that were relatively far apart. Furthermore, BPH17, BPH20 and BPH17-ptb on the NILs were selected by flanking markers with longer intervals because of the low density of polymorphic markers between donor parents and T65 around the chromosome 4S region. The intervals between each of the flanking marker pairs for BPH17 and BPH17-ptb were 3.7 Mbp, and that for BPH20 was 5.7 Mbp. Similarly, the interval for each of the two flanking markers for BPH2 and BPH21 on chromosome 12 was 4.1 Mbp because the exact locations of genes had not been identified before we started to develop the NILs by MAS. Therefore, many of the NILs had relatively long chromosome segments derived from the donor parents and there is a possibility that the remaining chromosomal segments from donor DNA around the target genes included linkage drag associated with susceptibility. In further research, ensuring that flanking markers are tightly linked to target genes will avoid linkage drag from donors through MAS and backcrossing.
In a previous study, virulence of a BPH population collected during 2005 in Japan had increased compared with the virulence of a population collected in 1966 [51]. Through antibiosis tests, we evaluated BPH resistance against the populations collected in 1966 (Hadano-66) and in 2013 (Koshi-2013). Both represented BPH arriving as migrants to Japan. The Hadano-66 population was virulent to T65 (with no resistance gene) but avirulent to all plants with resistance genes, including Mudgo (BPH1), ASD7 (BPH2), Rathu Heenai (BPH3 and BPH17), Babawee (BPH4), Chin Saba (BPH8), Balamawee (BPH9) and two NILs, BPH25-NIL and BPH26-NIL [51,71]. In the present study, most of the NILs, all of the PYLs and the donor parents were still effective against the Hadano-66 population. In contrast, all of the NILs and most of the PYLs were susceptible to the Koshi-2013 population, suggesting that BPH recently arriving to Japan from China has greater virulence than was evident about 50 years ago (i.e., 1966). Among the PYLs, two PYLs, BPH20 + BPH32-PYL and BPH2 + BPH3 + BPH17-PYL, had relatively high resistance, suggesting that PYLs with combinations of these genes are likely to provide good resistance against the current BPH populations that arrive to Japan (Koshi-2013). Finding new sources of resistance genes will be necessary to further improve resistance against contemporary BPH populations as they gain virulence.
In comparison to the corresponding NILs and PYLs, the resistance levels of PTB33, Rathu Heenati and IR71033-121-15 were higher. This suggests that PTB33, Rathu Heenati and IR71033-121-15 might also contain other BPH resistance gene(s). The other genetic factor(s) for BPH resistance can be revealed by analyzing the segregating populations derived from crosses between the developed PYLs and their donor parents in future studies. Additionally, Rathu Heenati had QBPH4.1 (5.8-7.8 Mbp) and QBPH4.2 (15.2-17.2 Mbp) on chromosome 4S rather than BPH3 and BPH17 [13]. Therefore, the NILs and PYLs carrying QBPH4.1 and QBPH4.2 should be developed and evaluated in further analyses. On the other hand, the lower resistance levels of the NILs and PYLs might be related to the relatively high ratio of substituted chromosomal segments from donors in the NILs (from 3.8 to 55.0 Mbp) and PYLs. There is a possibility that the retained donor chromosomal segments in the genetic background of the NILs and PYLs might be linked to the suppression of BPH resistance. To gain further knowledge of BPH resistance controlled by multiple genes, it will be essential to reduce the donor parent chromosomal segments on the NILs and PYLs by further backcrossing and MAS.

Plant Materials
To develop NILs with BPH resistance genes, a japonica rice variety, T65, that is susceptible to BPH, was used as a recurrent parent, and three rice varieties resistant to BPH were donor parents. The donor lines were IR71033-121-15, PTB33 and Rathu Heenati. IR71033-121-15 contains two BPH resistance genes, BPH20 and BPH21, from the wild rice species O. minuta (Accession number: IRGC101141) [17]. PTB33 (Accession number: IRGC19325) that originated from India contains BPH2, BPH17-ptb and BPH32. Rathu Heenati (Acc. no. IRGC 11730), that originated from Sri Lanka, carries BPH3 and BPH17 [16,39]. T65 was crossed with these donor parents and F 1 plants were backcrossed four times with T65 to generate BC 4 F 1 plants (Figure 3). At each generation of backcrossing, plants carrying BPH resistance genes from the donor parents were selected by MAS using flanking SSR markers of the target BPH resistance genes ( Table 1). The selected BC 4 F 1 plants were self-pollinated to produce BC 4 F 3 , BC 4 F 4 and BC 4 F 5 plants with BPH resistance genes. Finally, seven NILs with either BPH2, BPH3, BPH17, BPH20, BPH21, BPH32 or BPH17-ptb were developed. The NILs were used to survey the genetic background and evaluate BPH resistance levels as well as agronomic traits. Two additional NILs, BPH25-NIL and BPH26-NIL were used in the development of the PYLs [54].

The Development of PYLs with BPH Resistance Genes
All the PYLs for two or three BPH resistance genes were developed using the NILs descended from the BC4F1 generation, except BPH20 + BPH21-PYL and BPH32 + BPH17-ptb-PYL that were descended from the BC3F1 generation. The F1 plants derived from crosses between NILs were selfpollinated to produce F2 plants. From 96 F2 plants, plants that were homozygous for two or three BPH

The MAS for BPH Resistance Genes
To conduct MAS, approximately 2 cm of leaves from two-week old seedlings were collected and dried in a freeze drier for 48 h, and total DNA was extracted using the potassium acetate method [72]. The genotypes of SSR markers on plants in each generation were determined by polymerase chain reaction (PCR) and electrophoresis. The PCR amplification mix (8 µL) contained 3 µL of 1X GoTaq ® Green Master Mix (pH 8.5), 0.25 µM of primer and 4 µL of 20 times-diluted DNA. Each PCR amplification included one cycle at 96 • C for 5 min, 35 cycles at 96 • C for 30 s, 55 • C for 30 s and 72 • C for 30 s, followed by one extension cycle at 25 • C for 1 min. PCR products were analyzed by electrophoresis at 200 V using 4% agarose gel with 0.5 µg/mL ethidium bromide in 0.5X TBE buffer for 1 h and photographed under ultraviolet light. During MAS for resistance genes on chromosome 4S, the plants with BPH17 and BPH17-ptb, were selected using two markers, RM8213 and MS10, and the plants with BPH20 were selected using MS10 and RM5900 ( Table 1). The plants with BPH3 and BPH32 on the short arm of chromosome 6 were selected using two flanking markers, RM508 and RM588. The plants carrying BPH2 and BPH21 located on the long arm of chromosome 12 were screened using RM1246 and RM28493. The plants with BPH25 were selected using S00310 and MSSR1, and the plants with BPH26 were selected using RM309, RM28438, InD14, RM28466, RM28481 and MSSR2.

The Genetic Background Survey of the NILs
In the genetic background survey of the NILs, the bulk DNA from five plants was used. A total of 384 SSR markers distributed on 12 rice chromosomes were used during polymorphism tests with T65 and the donor parents [62]. Among the 384 SSR markers, 254 SSR markers with polymorphisms between IR71033-121-15 and T65 were utilized to identify substituted chromosomal segments from IR71033-121-15 on BPH20-NIL and BPH21-NIL. Additionally, 244 of 384 SSR markers with polymorphisms between PTB33 and T65 were used to detect substituted chromosomal segments from PTB33 on BPH2-NIL, BPH32-NIL and BPH17-ptb-NIL. To identify substituted chromosomal segments from Rathu Heenati on BPH3-NIL and BPH17-NIL, 204 of 384 SSR markers with polymorphisms between Rathu Heenati and T65 were used. The whole genome compositions of the developed NILs were graphically displayed following the concept of the graphical genotype proposed by Young and Tanksley (1989) using GGT software version 2.0 [73]. To evaluate resistance, an adaptation of the modified seedbox screening test (MSST) [45,74] was applied at 25 • C using the Hadano-66 strain. To conduct the test, 30 seeds of each of the NILs, PYLs and parent lines were sown to single rows in a plastic tray (23.0 × 30.0 × 2.5 cm) with 2.5 cm between successive rows of seedlings. Two sets of trays-one tray infested by BPH and the other without infestation (the control tray), were used to measure the effects of BPH on plant biomass. One row of Rathu Heenati was added as a resistant control, while three rows of T65 were sown at the center and the two edges as a susceptible control. At seven days after sowing (DAS), the plants in the trays were thinned to 20 plants per row. One tray was infested by the second and third instar nymphs at a density of around 20 BPHs per plant. The experiment was replicated three times. When all the plants of T65 were completely desiccated due to BPH feeding, the DSs of all lines were graded following the standard evaluation system for rice of the International Rice Research Institute [75]. The plants from each row in the two trays were cut above the soil surface and weighed. The fresh biomass reduction rate (FBRR) was calculated using the following formula: Fresh biomass reduction rate (FBRR)(%) = 1− Infested plant weight (g) Non-infested plant weight (g) ×100. (1)

Antibiosis Tests
Antibiosis tests were conducted at 25 • C following the method described by   [51]. Five plants of each NIL, PYL and parent line were individually sown in 200 mL plastic cups. At four weeks after sowing, the plants were trimmed to 15 cm height and covered with a plastic cage with insect screen windows for ventilation. Each cage was infested with five thin-abdomen brachypterous female BPHs. At five days after infestation, the ADM was recorded (i.e., the number of dead females).

Characterization of NILs and PYLs for Agronomic Traits
The NILs and PYLs were grown in a paddy field at Saga University (Saga, Japan) in 2018 and characterized for their agronomic traits compared to those of T65. Seedlings were transplanted at 28 DAS as one plant per hill, with 20 cm between hills and 25 cm between rows. Each entry was planted as at least three rows (12 plants per a row). Six agronomic traits: DTH, CL, PL, LL, LW and PN were measured for five plants in the same row. DTH was the days from sowing until 50% of panicles flowered. CL was measured from the soil surface to the panicle neck. PL is the length from tip to panicle neck of the longest panicle. The flag leaf width and length were measured from the largest and longest flag leaf of each sampled plant. Panicle number is the number of reproductive panicles of each plant at maturity.

Statistical Analysis
Mean values of BPH resistance (DS, FBRR and ADM) for the NILs and PYLs and agronomic traits were compared using one-way ANOVA. Dunnett's test and Tukey Kramer's test were conducted for multiple comparisons of BPH resistance and agronomic traits, respectively, using R software version 3.5.2. for developing materials and taking care of the plants. We also wish to thank the Government of Vietnam for the doctoral fellowship granted to C.D.N.

Conflicts of Interest:
The authors declare no conflict of interest.