Landraces of temperate japonica rice have superior alleles for improving culm strength associated with lodging resistance

Strong culm: a crucial trait for developing next-generation climate-resilient rice lines

Lodging, a phenomenon characterized by the bending or breaking of rice plants, poses substantial constraints on productivity, particularly during the harvesting phase in regions susceptible to strong winds. The rice strong culm trait is influenced by the intricate interplay of genetic, physiological, epigenetic, and environmental factors. Stem architecture, encompassing morphological and anatomical attributes, alongside the composition of both structural and non-structural carbohydrates, emerges as a critical determinant of lodging resistance. The adaptive response of the rice culm to various biotic and abiotic environmental factors further modulates the propensity for lodging. Advancements in next-generation sequencing technologies have expedited the genetic dissection of lodging resistance, enabling the identification of pertinent genes, quantitative trait loci, and novel alleles. Concurrently, contemporary breeding strategies, ranging from biparental approaches to more sophisticated methods such as multi-parent-based breeding, gene pyramiding, genomic selection, genome-wide association studies, and haplotype-based breeding, offer perspectives on the genetic underpinnings of culm strength. This review comprehensively delves into physiological attributes, culm histology, epigenetic determinants, and gene expression profiles associated with lodging resistance, with a specialized focus on leveraging next-generation sequencing for candidate gene discovery.

Genome-wide association studies for a comprehensive understanding of the genetic architecture of culm strength and yield traits in rice

Lodging resistance in rice is a complex trait determined by culm morphological and culm physical strength traits, and these traits are a major determinant of yield. We made a detailed analysis of various component traits with the aim of deriving optimized parameters for measuring culm strength. Genotyping by sequencing (GBS)-based genome-wide association study (GWAS) was employed among 181 genotypes for dissecting the genetic control of culm strength traits. The VanRaden kinship algorithm using 6,822 filtered single-nucleotide polymorphisms (SNPs) revealed the presence of two sub-groups within the association panel with kinship values concentrated at<0.5 level, indicating greater diversity among the genotypes. A wide range of phenotypic variation and high heritability for culm strength and yield traits were observed over two seasons, as reflected in best linear unbiased prediction (BLUP) estimates. The multi-locus model for GWAS resulted in the identification of 15 highly significant associations (p< 0.0001) for culm strength traits. Two novel major effect marker-trait associations (MTAs) for section modulus and bending stress were identified on chromosomes 2 and 12 with a phenotypic variance of 21.87% and 10.14%, respectively. Other MTAs were also noted in the vicinity of previously reported putative candidate genes for lodging resistance, providing an opportunity for further research on the biochemical basis of culm strength. The quantitative trait locus (QTL) hotspot identified on chromosome 12 with the synergistic association for culm strength trait (section modulus, bending stress, and internode breaking weight) and grain number can be considered a novel genomic region that can serve a dual purpose of enhancing culm strength and grain yield. Elite donors in the indica background with beneficial alleles of the identified major QTLs could be a valuable resource with greater significance in practical plant breeding programs focusing on improving lodging resistance in rice.

Genome-wide association study of leaf photosynthesis using a high-throughput gas exchange system in rice

Enhancing leaf photosynthetic capacity is essential for improving the yield of rice (Oryza sativa L.). Although the exploitation of natural genetic resources is considered a promising approach to enhance photosynthetic capacity, genomic factors related to the genetic diversity of leaf photosynthetic capacity have yet to be fully elucidated due to the limitation of measurement efficiency. In this study, we aimed to identify novel genomic regions for the net CO2 assimilation rate (A) by combining genome-wide association study (GWAS) and the newly developed rapid closed gas exchange system MIC-100. Using three MIC-100 systems in the field at the vegetative stage, we measured A of 168 temperate japonica rice varieties with six replicates for three years. We found that the modern varieties exhibited higher A than the landraces, while there was no significant relationship between the release year and A among the modern varieties. Our GWAS scan revealed two major peaks located on chromosomes 4 and 8, which were repeatedly detected in the different experiments and in the generalized linear modelling approach. We suggest that high-throughput gas exchange measurements combined with GWAS is a reliable approach for understanding the genetic mechanisms underlying photosynthetic diversities in crop species.

Identification of Genomic Regions for Deep-Water Resistance in Rice for Efficient Weed Control with Reduced Herbicide Use

Deep-water (DW) management in rice fields is a promising technique for efficient control of paddy weeds with reduced herbicide use. Maintaining a water depth of 10-20 cm for several weeks can largely suppress the weed growth, though it also inhibits rice growth because the DW management is usually initiated immediately after transplanting. Improving the DW resistance of rice during the initial growth stage is essential to avoid suppressing growth. In this study, we demonstrate a large genetic variation in the above-ground biomass (AGB) after the end of DW management among 165 temperate japonica varieties developed in Japan. Because the AGB closely correlated with plant length (PL) and tiller number (TN) at the early growth stage, we analyzed genomic regions associated with PL and TN by conducting a genome-wide association study. For PL, a major peak was detected on chromosome 3 (qPL3), which includes a gene encoding gibberellin biosynthesis, OsGA20ox1. The rice varieties with increased PL had a higher expression level of OsGA20ox1 as reported previously. For TN, a major peak was detected on chromosome 4 (qTN4), which includes NAL1 gene associated with leaf morphological development and panicle number. Although there was less difference in the expression level of NAL1 between genotypes, our findings suggest that an amino acid substitution in the exon region is responsible for the phenotypic changes. We also found that the rice varieties having alternative alleles of qPL3 and qTN4 showed significantly higher AGB than the varieties with the reference alleles. Our results suggest that OsGA20ox1 and NAL1 are promising genes for improving DW resistance in rice.

Thickness-stiffness trade-off improves lodging resistance in rice

Lodging of cereal crops significantly reduces grain yield and quality, making lodging resistance a prime target for breeding programs. However, lodging resistance among different rice (Oryza sativa L.) cultivars in the field remains largely unknown, as is the relationship between the major properties of culms such as their morphological and mechanical properties. Here, we investigated the morphological and mechanical properties of 12 rice cultivars by considering different internodes within culms. We detected variation in these two traits among cultivars: one set of cultivars had thicker but softer culms (thickness-type), while the other set of cultivars showed stiffer but thinner culms (stiffness-type). We designate this variation as a thickness-stiffness trade-off. We then constructed a mechanical model to dissect the mechanical and/or morphological constraints of rice culms subjected to their own weight (self-weight load). Through modeling, we discovered that ear weight and the morphology of the highest internode were important for reducing deflection, which may be important factors to achieve higher lodging resistance. The mechanical theory devised in this study could be used to predict the deflection of rice culms and may open new avenues for novel mechanics-based breeding techniques.

Effective use of legacy data in a genome-wide association studies improves the credibility of quantitative trait loci detection in rice

Genome-wide association studies (GWASs) are used to detect quantitative trait loci (QTL) using genomic and phenotypic data as inputs. While genomic data are obtained with high throughput and low cost, obtaining phenotypic data requires a large amount of effort and time. In past breeding programs, researchers and breeders have conducted a large number of phenotypic surveys and accumulated results as legacy data. In this study, we conducted a GWAS using phenotypic data of temperate japonica rice (Oryza sativa) varieties from a public database. The GWAS using the legacy data detected several known agriculturally important genes, indicating reliability of the legacy data for GWAS. By comparing the GWAS using legacy data (L-GWAS) and a GWAS using phenotypic data that we measured (M-GWAS), we detected reliable QTL for agronomically important traits. These results suggest that an L-GWAS is a strong alternative to replicate tests to confirm the reproducibility of QTL detected by an M-GWAS. In addition, because legacy data have often been accumulated for many traits, it is possible to evaluate the pleiotropic effect of the QTL identified for the specific trait that we focused on with respect to various other traits. This study demonstrates the effectiveness of using legacy data for GWASs and proposes the use of legacy data to accelerate genomic breeding.

Identification of Novel Quantitative Trait Loci for Culm Thickness of Rice Derived from Strong-Culm Landrace in Japan, Omachi

Increasing the lodging resistance of rice through genetic improvement has been an important target in breeding. To further enhance the lodging resistance of high-yielding rice varieties amidst climate change, it is necessary to not only shorten culms but strengthen them as well. A landrace rice variety, Omachi, which was established more than 100 years ago, has the largest culm diameter and bending moment at breaking in the basal internodes among 135 temperate japonica accessions. Using unused alleles in such a landrace is an effective way to strengthen the culm. In this study, we performed quantitative trait locus (QTL) analysis to identify the genetic factors of culm strength of Omachi using recombinant inbred lines (RILs) derived from a cross between Omachi and Koshihikari, a standard variety in Japan. We identified three QTLs for the culm diameter of the 5th internode on chromosomes 3 (qCD3) and 7 (qCD7-1, qCD7-2). Among them, qCD7-2 was verified by QTL analysis using the F₂ population derived from a cross between one of the RILs and Koshihikari. RNA-seq analysis of shoot apex raised 10 candidate genes underlying the region of qCD7-2. The increase in culm strength by accumulating Omachi alleles of qCD3, qCD7-1 and qCD7-2 was 25.0% in 2020. These QTLs for culm diameter pleiotropically increased spikelet number per panicle but did not affect days to heading or culm length. These results suggest that the Omachi alleles of qCD3, qCD7-1 and qCD7-2 are useful for breeding to increase lodging resistance and yield.

Lodging resistance of rice plants studied from the perspective of culm mechanical properties, carbon framework, free volume, and chemical composition

In this study, four varieties of rice were cultivated on the same farmland under same conditions and for same duration. However, their lodging resistance was found to be obviously different from each other. Herein, three key factors that highly influenced the lodging resistance were identified. First, in terms of morphological trait, in contrast to the generally believed theory that the overall thickness of the fresh culm wall governs the strength, the thickness of the depressed region of the dried basal culm wall largely determined the mechanical properties by acting as the weak link. This depressed region represents the vulnerable part with high syneresis rate. Second, the culm and its carbon framework exhibited sufficient strength and rigidity for both support and stability of the rice stem. The constraint of high lodging resistance of rice plants is attributed to the culm flexibility. Furthermore, the results of the positron annihilation lifetime spectroscopy corroborate that the most amorphous part and the highest-fraction free volume in the culm carbon framework were found for samples that exhibited high lodging resistance. This result confirmed the significant influence of the culm flexibility on lodging resistance. Third, a higher level of nitrogen element content in the basal culm can benefit its growth and development, which may contribute to an increase in lodging resistance of rice plants.

Unveiling the genetic architecture for lodging resistance in rice (Oryza sativa. L) by genome-wide association analyses

Lodging is one of the major abiotic stresses, affecting the total crop yield and quality. The improved lodging resistance and its component traits potentially reduce the yield losses. The section modulus (SM), bending moment at breaking (M), pushing resistance (PR), and coefficient of lodging resistance (cLr) are the key elements to estimate the lodging resistance. Understanding the genetic architecture of lodging resistance-related traits will help to improve the culm strength and overall yield potential. In this study, a natural population of 795 globally diverse genotypes was further divided into two (indica and japonica) subpopulations and was used to evaluate the lodging resistance and culm strength-related traits. Significant diversity was observed among the studied traits. We carried out the genome-wide association evaluation of four lodging resistance traits with 3.3 million deep resolution single-nucleotide polymorphic (SNP) markers. The general linear model (GLM) and compressed mixed linear model (MLM) were used for the whole population and two subpopulation genome-wide association studies (GWAS), and a 1000-time permutation test was performed to remove the false positives. A total of 375 nonredundant QTLs were observed for four culm strength traits on 12 chromosomes of the rice genome. Then, 33 pleiotropic loci governing more than one trait were mined. A total of 4031 annotated genes were detected within the candidate genomic region of 33 pleiotropic loci. The functional annotations and metabolic pathway enrichment analysis showed cellular localization and transmembrane transport as the top gene ontological terms. The in silico and in vitro expression analyses were conducted to validate the three candidate genes in a pleiotropic QTL on chromosome 7. It validated OsFBA2 as a candidate gene to contribute to lodging resistance in rice. The haplotype analysis for the candidate gene revealed a significant functional variation in the promoter region. Validation and introgression of alleles that are beneficial to induce culm strength may be used in rice breeding for lodging resistance.

Identification of a unique allele in the quantitative trait locus for crown root number in japonica rice from Japan using genome-wide association studies

To explore the genetic resources that could be utilized to help improve root system architecture phenotypes in rice (Oryza sativa), we have conducted genome-wide association studies to investigate maximum root length and crown root number in 135 10-day-old Japanese rice accessions grown hydroponically. We identified a quantitative trait locus for crown root number at approximately 32.7 Mbp on chromosome 4 and designated it qNCR1 (quantitative trait locus for Number of Crown Root 1). A linkage disequilibrium map around qNCR1 suggested that three candidate genes are involved in crown root number: a cullin (LOC_Os04g55030), a gibberellin 20 oxidase 8 (LOC_Os04g55070), and a cyclic nucleotide-gated ion channel (LOC_Os04g55080). The combination of haplotypes for each gene was designated as a haploblock, and haploblocks 1, 2, and 3 were defined. Compared to haploblock 1, the accessions with haploblocks 2 and 3 had fewer crown roots; approximately 5% and 10% reductions in 10-day-old plants and 15% and 25% reductions in 42-day-old plants, respectively. A Japanese leading variety Koshihikari and its progenies harbored haploblock 3. Their crown root number could potentially be improved using haploblocks 1 and 2.

Potential of rice landraces with strong culms as genetic resources for improving lodging resistance against super typhoons

It is generally believed that rice landraces with long culms are susceptible to lodging, and have not been utilized for breeding to improve lodging resistance. However, little is known about the structural culm strength of landraces and their beneficial genetic loci. Therefore, in this study, genome-wide association studies (GWAS) were performed using a rice population panel including Japanese rice landraces to identify beneficial loci associated with strong culms. As a result, the landraces were found to have higher structural culm strength and greater diversity than the breeding varieties. Genetic loci associated with strong culms were identified, and it was demonstrated that haplotypes with positive effects of those loci were present in a high proportion of these landraces. These results indicated that the utilization of the strong culm-associated loci present in Japanese rice landraces may further improve the lodging resistance of modern breeding varieties that have relied on semi-dwarfism.