Plant materials
CSSLs with chromosomal segments of “Nona Bokra” with the genetic background of “Koshihikari”, developed by Takai et al. (2007), were used for QTL analysis from 2016 to 2018. Seeds of 44 CSSLs and “Koshihikari” were sown in a greenhouse in late April 2016 and 2017, and four CSSLs with the chromosomal segment on chromosome 4 (line nos. 511–514) were sown in 2018. In late May, 10 plants were transplanted, with a single plant per hill spaced at 30 × 15 cm, into a paddy field in Utsunomiya, Japan (latitude 36°N, longitude 139°E), and grown under natural conditions. For QTL characterization, CSSL512 (CSSL-PRL4) and “Koshihikari” were grown in each plot of 10 rows with 12 hills per row in 2017 and 2018. These seeds were sown in late April and transplanted to a single plant per hill, spaced at 30 × 15 cm, in late May. In 2019 and 2020, the NIL with PRL4 (NIL-PRL4) and “Koshihikari” was grown under the same cultivation conditions as in 2017 and 2018. Additionally, NIL-LPW4 was grown in each plot of two rows, with 12 hills per row, in 2019 and 2020. N, P2O5, and K2O fertilizers were applied as basal dressing at 6 g m− 2 before transplanting, respectively, in 2016–2020.
QTL analysis
QTLs for pushing resistance of the lower part per plant or tiller, tiller number, and days to heading were determined by comparing the phenotypes and the genotypes of 140 simple sequence repeats and five expressed sequence tag markers (Takai et al. 2007). QTL detection was based on a t-test, with the threshold at a probability level of 0.05, of the difference between each CSSL and “Koshihikari”, according to Madoka et al. (2008). The additive effect of a QTL from the “Nona Bokra” allele was expressed as positive (plus) or negative (minus) relative percentage compared with the “Koshihikari” phenotype.
Selection of CSSL and NIL with PRL4
From 44 CSSLs, only CSSL512, carrying 10.3 Mb of the “Nona Bokra” homozygous allele, was selected to characterize the QTL effect for the PRL4 position on chromosome 4. After the backcross of CSSL512 with “Koshihikari”, the genotypes of 288 BC1F2 plants were analyzed using eight SSR markers (International Rice Genome Sequencing Project, 2005). NIL-PRL4 was selected from 288 BC1F3 plants by genotyping using 13 SSR markers. NIL-PRL4 had 6.6 Mb of the “Nona Bokra” homozygous allele, localized nearest PRL4 marker RM5749 (Fig. 3a). In addition, because the locus for low panicle weight, which was named LPW4 in this study, was identified at the nearest marker RM3308 by Ujiie et al. (2012), NIL-LPW4, carrying 3.2 Mb of “Nona Bokra” homozygous allele between RM7113 and RM3308, was selected from the BC1F3 plants (Supplementary Fig. 1a). DNA extraction was conducted using TPS buffer according to Monna et al. (2002). PCR was carried out using KAPA Taq Extra (Kapa Biosystems, USA), and the PCR profile was 94°C for 2 min, 35 cycles of 94°C for 30 sec, 55°C for 15 s, 72°C for 30 s, and a final extension at 72°C for 5 min. The resultant PCR products were visualized using 4% agarose gel stained with Midori Green Advance (Nippon Genetics, Tokyo, Japan).
Lodging resistance traits
At seven weeks after heading (7WAH), the pushing resistance of the lower part (PRL) was analyzed according to Kashiwagi and Ishimaru (2004). The upper part of each plant was removed at a height of 40 cm, and a prostrate tester (Daiki Rika Kogyou Co., Tokyo, Japan) was placed vertically to the culm at a height of 20 cm above the soil surface. The pushing resistance was measured until the plants had inclined to an angle of 45° vertically. The pushing resistance of the lower part per tiller (PRL/TN) was calculated after measurement of tiller number. Physical properties of a basal internode were analyzed using universal table-top testing instruments (EZ-SX, Shimadzu, Kyoto, Japan). An average-sized basal internode was compressed with a compression jig (1 mm radius × 80 mm width) at a constant velocity of 1 mm min− 1, and compression resistance, strain at the maximum compression resistance, and culm stiffness were measured. The angle of each plant to the ground was measured using a protractor under field conditions. All lodging resistance traits were determined using six plants in each line.
Morphological traits, SPAD value of leaves, and yield components
Plant length and tiller number were measured before the harvest season, and the lengths and diameters of six average basal internodes per plant were measured at 7WAH. In 2018, basal internodes of “Koshihikari” and CSSL-PRL4 were measured to determine the thicknesses of the culm wall and cortical fiber tissue. Sampled basal culms were dried at 80°C for 3 d, and the dried weight per unit length was calculated. SPAD values of the top three leaves were measured using a chlorophyll meter (SPAD-502, KONICA MINOLTA Inc., Tokyo, Japan) immediately after heading and at 7WAH. Measurements of morphological traits and SPAD values were conducted in six replicates. The panicles of six plants from each line were randomly sampled at 7WAH, and the yield components were measured at the individual plant level.
Carbohydrate content of basal culms
Dried basal culms were ground to a powder using a Wonder Blender (Osaka Chemical Co., Osaka, Japan). Holocellulose content was measured in accordance with the sodium chlorite method of Wise et al. (1946). A powdered sample of 0.3 g was degreased in ethanol–benzene (1:2 v/v) with a Soxhlet extraction for 6 h and then air-dried. After drying at 80°C for 1 h, the sample was added to hot distilled water at 80°C and delignified by the addition of sodium chlorite and glacial acetic acid at 80°C for 1 h. This process was repeated four times. The sample solution was filtrated through a glass filter, and the residue was washed in cool distilled water and acetone. After drying at 105°C for 1 h, the residue was cooled to room temperature in a desiccator and the holocellulose content was weighed. Lignin content was measured in accordance with Suzuki et al. (2009). A powdered sample of approximately 15 mg was transferred to a microcentrifuge tube with a screw cap and dried at 60°C for 1 h. The sample was extracted with distilled water and centrifuged for 10 min at 16,200 × g. The pellet was extracted in methanol at 60°C for 20 min and then centrifuged for 10 min at 16,200 × g. This process was then repeated. After drying with a rotary evaporator (CC-105, TOMY, Tokyo, Japan), 3M HCl and thioglycolic acid were added to the pellet. The sample was heated at 80°C for 3 h and centrifuged for 10 min at 16,200 × g. The pellet was vortexed in distilled water and resuspended in 1 M NaOH for 16 h with a culture rotator (RT-50, TAITEC, Saitama, Japan). After centrifugation, the supernatant was acidified with concentrated HCl and chilled at 4°C for 4 h. The sample was then centrifuged and the pellet was dissolved and diluted with 1 M NaOH. The absorbance at 280 nm of diluted TGAL solution was recorded using a V-730BIO spectrophotometer (JASCO, Tokyo, Japan). NSC content was measured in accordance with Ishimaru et al. (2001). A powdered sample of approximately 20 mg was extracted twice with 80% (v/v) ethanol at room temperature, once with 80% ethanol at 80°C, and then centrifuged at 2,400 × g for 5 min. The supernatant was collected and dried in a vacuum, and then the sample was dissolved in distilled water, centrifuged at 14,000 × g for 5 min, and used for the determination of sucrose and hexoses by an enzymatic method (Bergmeyer and Bernt 1974). To determine starch content, the pellet was boiled in distilled water for 1 h and then digested with amyloglucosidase for 20 min at 55°C. The resultant hexoses were determined by the enzymatic method. Six plants from each line were used to determine carbohydrate content.
Statistical analysis
Comparisons between “Koshihikari” and the CSSLs or NIL-PRL4 were carried out using a t-test at a probability of P < 0.05 only when the F-test showed significance at 0.05 probability levels with Microsoft Excel 2018 for Mac.