Genetic analyses of lodging resistance and yield provide insights into post-Green-Revolution breeding in rice

Genetic Dissection of Major Rice QTLs for Strong Culms and Fine Mapping of qWS5 for Breeding Application in Transplanted System

BACKGROUND

Rice is one of the major staples that feeds about one half of the global populations, and it is important to identify the genetic loci for the traits related to yield improvement. Lodging will cause severe yield loss when it happens, and stem diameter has been characterized as an important trait for lodging resistance. However, most QTLs for stem diameter have not been finely dissected due to their sensitivity to environmental fluctuation.

RESULT

In this study, we performed QTL analysis for stem diameter using populations derived from Nipponbare (NIP) and strong culm variety YYP1, and confirmed the single and combined effect of three major QTLs by recombinant inbred lines (RILs). Based on the QTL location, we found that qWS5 is a novel QTL not well characterized before. To finely dissect the novel locus, several recombinant heterogeneous inbred families (HIFs) were selected from the RILs for linkage analysis and their derived nearly isogenic lines (NILs) were subjected to detailed trait investigation throughout different years. The HIF-NILs strategy confined the QTL to about 380 kb region supported by repeated genotype and phenotype data, and it lays the foundation for QTL cloning in the future. In addition, introgression of the QTL to an elite japonica variety SD785 was performed by successive backcrossing, and it confirmed the value of qWS5 in increasing stem diameter and other agronomic traits during rice breeding.

CONCLUSIONS

We prove that qWS5 is a novel QTL with relatively stable effect for stem diameter and the QTL can be finely mapped to small region by the HIF-NILs strategy. The result will facilitate the improvement of rice lodging resistance by molecular marker assisted selection breeding.

The OsEIL1-OsWOX11 transcription factor module controls rice crown root development in response to soil compaction

Optimizing the root architecture of crops is an effective strategy for improving crop yields. Soil compaction is a serious global problem that limits crop productivity by restricting root growth, but the underlying molecular mechanisms are largely unclear. Here, we show that ethylene stimulates rice (Oryza sativa) crown root development in response to soil compaction. First, we demonstrate that compacted soil promotes ethylene production and the accumulation of ETHYLENE INSENSITIVE 3-LIKE 1 (OsEIL1) in rice roots, stimulating crown root primordia initiation and development, thereby increasing crown root number in lower stem nodes. Through transcriptome profiling and molecular analyses, we reveal that OsEIL1 directly activates the expression of WUSCHEL-RELATED HOMEOBOX 11 (OsWOX11), an activator of crown root emergence and growth, and that OsWOX11 mutations delay crown root development, thus impairing the plant's response to ethylene and soil compaction. Genetic analysis demonstrates that OsWOX11 functions downstream of OsEIL1. In summary, our results demonstrate that the OsEIL1-OsWOX11 module regulates ethylene action during crown root development in response to soil compaction, providing a strategy for the genetic modification of crop root architecture and grain agronomic traits.

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.

Genetic analyses and prediction for lodging‑related traits in a diverse Iranian hexaploid wheat collection

Lodging is one of the most important limiting environmental factors for achieving the maximum yield and quality of grains in cereals, including wheat. However, little is known about the genetic foundation underlying lodging resistance (LR) in wheat. In this study, 208 landraces and 90 cultivars were phenotyped in two cropping seasons (2018-2019 and 2019-2020) for 19 LR-related traits. A genome-wide association study (GWAS) and genomics prediction were carried out to dissect the genomic regions of LR. The number of significant marker pairs (MPs) was highest for genome B in both landraces (427,017) and cultivars (37,359). The strongest linkage disequilibrium (LD) between marker pairs was found on chromosome 4A (0.318). For stem lodging-related traits, 465, 497, and 478 marker-trait associations (MTAs) and 45 candidate genes were identified in year 1, year 2, and pooled. Gene ontology exhibited genomic region on Chr. 2B, 6B, and 7B control lodging. Most of these genes have key roles in defense response, calcium ion transmembrane transport, carbohydrate metabolic process, nitrogen compound metabolic process, and some genes harbor unknown functions that, all together may respond to lodging as a complex network. The module associated with starch and sucrose biosynthesis was highlighted. Regarding genomic prediction, the GBLUP model performed better than BRR and RRBLUP. This suggests that GBLUP would be a good tool for wheat genome selection. As a result of these findings, it has been possible to identify pivotal QTLs and genes that could be used to improve stem lodging resistance in Triticum aestivum L.

OsTCP19 coordinates inhibition of lignin biosynthesis and promotion of cellulose biosynthesis to modify lodging resistance in rice

Lignin and cellulose are two essential elements of plant secondary cell walls that shape the mechanical characteristics of the culm to prevent lodging. However, how the regulation of the lignin and cellulose composition is combined to achieve optimal mechanical characteristics is unclear. Here, we show that increasing OsTCP19 expression in rice coordinately repressed lignin biosynthesis and promoted cellulose biosynthesis, resulting in enhanced lodging resistance. In contrast, repression of OsTCP19 coordinately promoted lignin biosynthesis and inhibited cellulose biosynthesis, leading to greater susceptibility to lodging. We found that OsTCP19 binds to the promoters of both MYB108 and MYB103L to increase their expression, with the former being responsible for repressing lignin biosynthesis and the latter for promoting cellulose biosynthesis. Moreover, up-regulation of OsTCP19 in fibers improved grain yield and lodging resistance. Thus, our results identify the OsTCP19-OsMYB108/OsMYB103L module as a key regulator of lignin and cellulose production in rice, and open up the possibility for precisely manipulating lignin-cellulose composition to improve culm mechanical properties for lodging resistance.

Dissection of the Genetic Basis of Resistance to Stem Rot in Cultivated Peanuts (Arachis hypogaea L.) through Genome-Wide Association Study

Peanut (Arachis hypogaea) is an important oilseed and cash crop worldwide, contributing an important source of edible oil and protein for human nutrition. However, the incidence of stem rot disease caused by Athelia rolfsii poses a major challenge to peanut cultivation, resulting in significant yield losses. In this study, a panel of 202 peanut accessions was evaluated for their resistance to stem rot by inoculating plants in the field with A. rolfsii-infested oat grains in three environments. The mean disease index value of each environment for accessions in subsp. fasitigiate and subsp. hypogaea showed no significant difference. Accessions from southern China displayed the lowest disease index value compared to those from other ecological regions. We used whole-genome resequencing to analyze the genotypes of the accessions and to identify significant SNPs associated with stem rot resistance through genome-wide association study (GWAS). A total of 121 significant SNPs associated with stem rot resistance in peanut were identified, with phenotypic variation explained (PVE) ranging from 12.23% to 15.51%. A total of 27 candidate genes within 100 kb upstream and downstream of 23 significant SNPs were annotated, which have functions related to recognition, signal transduction, and defense response. These significant SNPs and candidate genes provide valuable information for further validation and molecular breeding to improve stem rot resistance in peanut.

SoybeanGDB: A comprehensive genomic and bioinformatic platform for soybean genetics and genomics

Soybean (Glycine max (L.) Merr.) is a globally significant crop, widely cultivated for oilseed production and animal feeds. In recent years, the rapid growth of multi-omics data from thousands of soybean accessions has provided unprecedented opportunities for researchers to explore genomes, genetic variations, and gene functions. To facilitate the utilization of these abundant data for soybean breeding and genetic improvement, the SoybeanGDB database (https://venyao.xyz/SoybeanGDB/) was developed as a comprehensive platform. SoybeanGDB integrates high-quality de novo assemblies of 39 soybean genomes and genomic variations among thousands of soybean accessions. Genomic information and variations in user-specified genomic regions can be searched and downloaded from SoybeanGDB, in a user-friendly manner. To facilitate research on genetic resources and elucidate the biological significance of genes, SoybeanGDB also incorporates a variety of bioinformatics analysis modules with graphical interfaces, such as linkage disequilibrium analysis, nucleotide diversity analysis, allele frequency analysis, gene expression analysis, primer design, gene set enrichment analysis, etc. In summary, SoybeanGDB is an essential and valuable resource that provides an open and free platform to accelerate global soybean research.

SEMI-ROLLED LEAF 10 stabilizes catalase isozyme B to regulate leaf morphology and thermotolerance in rice (Oryza sativa L.)

Plant architecture and stress tolerance play important roles in rice breeding. Specific leaf morphologies and ideal plant architecture can effectively improve both abiotic stress resistance and rice grain yield. However, the mechanism by which plants simultaneously regulate leaf morphogenesis and stress resistance remains elusive. Here, we report that SRL10, which encodes a double-stranded RNA-binding protein, regulates leaf morphology and thermotolerance in rice through alteration of microRNA biogenesis. The srl10 mutant had a semi-rolled leaf phenotype and elevated sensitivity to high temperature. SRL10 directly interacted with catalase isozyme B (CATB), and the two proteins mutually increased one other's stability to enhance hydrogen peroxide (H₂ O₂ ) scavenging, thereby contributing to thermotolerance. The natural Hap3 (AGC) type of SRL10 allele was found to be present in the majority of aus rice accessions, and was identified as a thermotolerant allele under high temperature stress in both the field and the growth chamber. Moreover, the seed-setting rate was 3.19 times higher and grain yield per plant was 1.68 times higher in near-isogenic line (NIL) carrying Hap3 allele compared to plants carrying Hap1 allele under heat stress. Collectively, these results reveal a new locus of interest and define a novel SRL10-CATB based regulatory mechanism for developing cultivars with high temperature tolerance and stable yield. Furthermore, our findings provide a theoretical basis for simultaneous breeding for plant architecture and stress resistance.

Transcriptome analysis provides insights into lignin synthesis and MAPK signaling pathway that strengthen the resistance of bitter gourd (Momordica charantia) to Fusarium wilt

Fusarium wilt is a typical soil-borne disease caused by Fusarium oxysporum f. sp. momordicae (FOM) in bitter gourd. In this study, by comparing sequencing data at multiple time points and considering the difference between resistant (R) and susceptible (S) varieties, differentially expressed genes were screened out. Short time-series expression miner analysis revealed the upregulated expression trend of genes, which were enriched in phenylpropanoid biosynthesis, plant-pathogen interaction, and mitogen-activated protein kinase signaling pathway. Further, observation of the microstructure revealed that the R variety may form tyloses earlier than the S variety to prevent mycelium diffusion from the xylem vessel. After Fusarium wilt infection, the enzymatic activities of superoxide dismutase, peroxidase, phenylalanine ammonia lyase, and catalaseas well as levels of superoxide anion and malondialdehyde were increased in the R variety higher than those in the S variety. This study provides a reference to elucidate the disease resistance mechanism of bitter gourd.

Identification of QTL underlying the main stem related traits in a doubled haploid barley population

Lodging reduces grain yield in cereal crops. The height, diameter and strength of stem are crucial for lodging resistance, grain yield, and photosynthate transport in barley. Understanding the genetic basis of stem benefits barley breeding. Here, we evaluated 13 stem related traits after 28 days of heading in a barley DH population in two consecutive years. Significant phenotypic correlations between lodging index (LI) and other stem traits were observed. Three mapping methods using the experimental data and the BLUP data, detected 27 stable and major QTLs, and 22 QTL clustered regions. Many QTLs were consistent with previously reported traits for grain filling rate, internodes, panicle and lodging resistance. Further, candidate genes were predicted for stable and major QTLs and were associated with plant development and adverse stress in the transition from vegetative stage to reproductive stage. This study provided potential genetic basis and new information for exploring barley stem morphology, and laid a foundation for map-based cloning and further fine mapping of these QTLs.

RNA-Seq and Genome-Wide Association Studies Reveal Potential Genes for Rice Seed Shattering

The loss of the shattering ability is one of the key events in rice domestication. The strength of the seed shattering ability is closely related to the harvest yield and the adaptability of modern mechanical harvesting methods. In this study, using a population of 587 natural rice cultivars, quantitative trait loci associated with seed shattering were detected by genome-wide association studies (GWASs). We consider the quantitative trait loci (QTLs) qBTS1 and qBTS3 to be the key loci for seed shattering in rice. Additionally, the abscission zone (AZ) and nonabscission zone (NAZ) of materials with a loss of shattering (DZ129) and easy shattering (W517) were subjected to RNA-Seq, and high-quality differential expression profiles were obtained. The AZ-specific differentially expressed genes (DEGs) of W517 were significantly enriched in plant hormone signal transduction, while the AZ-specific DEGs of DZ129 were enriched in phenylpropanoid biosynthesis. We identified candidate genes for the lignin-associated laccase precursor protein (LOC_Os01g63180) and the glycoside hydrolase family (LOC_Os03g14210) in the QTLs qBTS1 (chromosome 1) and qBTS3 (chromosome 3), respectively. In summary, our findings lay the foundation for the further cloning of qBTS1 and qBTS3, which would provide new insights into seed shattering in rice.

Shading decreases lodging resistance of wheat under different planting densities by altering lignin monomer composition of stems

To clarify the influences of shading stress and planting density on the lignin monomer composition of wheat stems and their relationship with lodging resistance, Lodging resistant variety Shannong 23 (SN23) and lodging sensitive variety Shannong 16 (SN16) were grown during 2018-2019 and 2019-2020 growing seasons. The planting densities were 150 × 10⁴ plants ha^-1 (D1), 225 × 10⁴ plants ha^-1 (D2) and 300 × 10⁴ plants ha^-1 (D3). At the jointing stage, an artificial shading shed was used to simulate shading stress. Then the effects of shading on stem morphological characteristics, lignin monomer composition and lodging resistance of wheat under different planting densities were studied. Results indicate that shading at the jointing stage increased the length of basal internodes and the plant height and moved the height of center of gravity (CG) upward. Moreover, the stem diameter and the wall thickness decreased by 0.10-0.53 mm and 0.18-0.40 mm, respectively. The stem filling degree was reduced accordingly. As indicated by the correlation analysis and the stepwise regression analysis, shading-induced lodging mainly resulted from changes in the stem morphological characteristics and lignin accumulation. The influential magnitude of these factors was ordered as follows: stem filling degree, wall thickness, lignin content, contents and proportions of monomers S and H, and length of the second internode. The expression abundance of TaPAL, TaCOMT, TaCCR, and TaCAD declined in response to shading stress and high planting density. As a result, the distribution ratios of photosynthetic carbon sources to lignin monomers S, G and H were changed. The lignin content of stems on the day 42 after the jointing stage decreased by 18.48%. The monomer S content decreased, while the content and proportion of monomer H increased, thus weakening the breaking strength of wheat stems.

Genome-wide association study reveals a GLYCOGEN SYNTHASE KINASE 3 gene regulating plant height in Brassica napus

Rapeseed (Brassica napus) is an allotetraploid crop that is the main source of edible oils and feed proteins in the world. The ideal plant architecture breeding is a major objective of rapeseed breeding and determining the appropriate plant height is a key element of the ideal plant architecture. Therefore, this study aims to improve the understanding of the genetic controls underlying plant height. The plant heights of 230 rapeseed accessions collected worldwide were investigated in field experiments over two consecutive years in Wuhan, China. Whole-genome resequencing of these accessions yielded a total of 1,707,194 informative single nucleotide polymorphisms (SNPs) that were used for genome-wide association analysis (GWAS). GWAS and haplotype analysis showed that BnaA01g09530D, which encodes BRASSINOSTEROID-INSENSITIVE 2 and belongs to the GLYCOGEN SYNTHASE KINASE 3 (GSK3) family, was significantly associated with plant height in B. napus. Moreover, a total of 31 BnGSK3s with complete domains were identified from B. napus genome and clustered into four groups according to phylogenetic analysis, gene structure, and motif distribution. The expression patterns showed that BnGSK3s exhibited significant differences in 13 developmental tissues in B. napus, suggesting that BnGSK3s may be involved in tissue-specific development. Sixteen BnGSK3 genes were highly expressed the in shoot apical meristem, which may be related to plant height or architecture development. These results are important for providing new haplotypes of plant height in B. napus and for extending valuable genetic information for rapeseed genetic improvement of plant architecture.

A natural uORF variant confers phosphorus acquisition diversity in soybean

Phosphorus (P) is an essential element for all organisms. Because P fertilizers are a non-renewable resource and high fixation in soils, sustainable agriculture requires researchers to improve crop P acquisition efficiency. Here, we report a strong association signal at a locus of CPU1 (component of phosphorus uptake 1), from a genome-wide association study of P acquisition efficiency in a soybean core collection grown in the field. A SEC12-like gene, GmPHF1, is identified as the causal gene for CPU1. GmPHF1 facilitates the ER (endoplasmic reticulum) exit of the phosphate transporter, GmPT4, to the plasma membrane of root epidermal cells. A common SNP in an upstream open reading frame (uORF) of GmPHF1, which alters the abundance of GmPHF1 in a tissue-specific manner, contributes to P acquisition diversity in soybean. A natural genetic variation conditions diversity in soybean P acquisition, which can be used to develop P-efficient soybean genotypes.

Combined BSA-Seq Based Mapping and RNA-Seq Profiling Reveal Candidate Genes Associated with Plant Architecture in Brassica napus

Plant architecture involves important agronomic traits affecting crop yield, resistance to lodging, and fitness for mechanical harvesting in Brassica napus. Breeding high-yield varieties with plant architecture suitable for mechanical harvesting is the main goal of rapeseed breeders. Here, we report an accession of B. napus (4942C-5), which has a dwarf and compact plant architecture in contrast to cultivated varieties. A BC₈ population was constructed by crossing a normal plant architecture line, 8008, with the recurrent parent 4942C-5. To investigate the molecular mechanisms underlying plant architecture, we performed phytohormone profiling, bulk segregant analysis sequencing (BSA-Seq), and RNA sequencing (RNA-Seq) in BC₈ plants with contrasting plant architecture. Genetic analysis indicated the plant architecture traits of 4942C-5 were recessive traits controlled by multiple genes. The content of auxin (IAA), gibberellin (GA), and abscisic acid (ABA) differed significantly between plants with contrasting plant architecture in the BC₈ population. Based on BSA-Seq analysis, we identified five candidate intervals on chromosome A01, namely those of 0 to 6.33 Mb, 6.45 to 6.48 Mb, 6.51 to 6.53 Mb, 6.77 to 6.79 Mb, and 7 to 7.01 Mb regions. The RNA-Seq analysis revealed a total of 4378 differentially expressed genes (DEGs), of which 2801 were up-regulated and 1577 were down-regulated. There, further analysis showed that genes involved in plant hormone biosynthesis and signal transduction, cell structure, and the phenylpropanoid pathway might play a pivotal role in the morphogenesis of plant architecture. Association analysis of BSA-Seq and RNA-Seq suggested that seven DEGs involved in plant hormone signal transduction and a WUSCHEL-related homeobox (WOX) gene (BnaA01g01910D) might be candidate genes responsible for the dwarf and compact phenotype in 4942C-5. These findings provide a foundation for elucidating the mechanisms underlying rapeseed plant architecture and should contribute to breed new varieties suitable for mechanization.

Identification and Allele Combination Analysis of Rice Grain Shape-Related Genes by Genome-Wide Association Study

Grain shape is an important agronomic character of rice, which affects the appearance, processing, and the edible quality. Screening and identifying more new genes associated with grain shape is beneficial to further understanding the genetic basis of grain shape and provides more gene resources for genetic breeding. This study has a natural population containing 623 indica rice cultivars. Genome-wide association studies/GWAS of several traits related to grain shape (grain length/GL, grain width/GW, grain length to width ratio/GLWR, grain circumferences/GC, and grain size/grain area/GS) were conducted by combining phenotypic data from four environments and the second-generation resequencing data, which have identified 39 important Quantitative trait locus/QTLs. We analyzed the 39 QTLs using three methods: gene-based association analysis, haplotype analysis, and functional annotation and identified three cloned genes (GS3, GW5, OsDER1) and seven new candidate genes in the candidate interval. At the same time, to effectively utilize the genes in the grain shape-related gene bank, we have also analyzed the allelic combinations of the three cloned genes. Finally, the extreme allele combination corresponding to each trait was found through statistical analysis. This study’s novel candidate genes and allele combinations will provide a valuable reference for future breeding work.

Genome-Wide Association Study of the Genetic Basis of Effective Tiller Number in Rice

BACKGROUND

Effective tiller number (ETN) has a pivotal role in determination of rice (Oryza sativa L.) grain yield. ETN is a complex quantitative trait regulated by both genetic and environmental factors. Despite multiple tillering-related genes have been cloned previously, few of them have been utilized in practical breeding programs.

RESULTS

In this study, we conducted a genome-wide association study (GWAS) for ETN using a panel of 490 rice accessions derived from the 3 K rice genomes project. Thirty eight ETN-associated QTLs were identified, interestingly, four of which colocalized with the OsAAP1, DWL2, NAL1, and OsWRKY74 gene previously reported to be involved in rice tillering regulation. Haplotype (Hap) analysis revealed that Hap5 of OsAAP1, Hap3 and 6 of DWL2, Hap2 of NAL1, and Hap3 and 4 of OsWRKY74 are favorable alleles for ETN. Pyramiding favorable alleles of all these four genes had more enhancement in ETN than accessions harboring the favorable allele of only one gene. Moreover, we identified 25 novel candidate genes which might also affect ETN, and the positive association between expression levels of the OsPILS6b gene and ETN was validated by RT-qPCR. Furthermore, transcriptome analysis on data released on public database revealed that most ETN-associated genes showed a relatively high expression from 21 days after transplanting (DAT) to 49 DAT and decreased since then. This unique expression pattern of ETN-associated genes may contribute to the transition from vegetative to reproductive growth of tillers.

CONCLUSIONS

Our results revealed that GWAS is a feasible way to mine ETN-associated genes. The candidate genes and favorable alleles identified in this study have the potential application value in rice molecular breeding for high ETN and grain yield.

Genome-Wide Association Study Identified Novel Candidate Loci/Genes Affecting Lodging Resistance in Rice

Lodging reduces rice yield, but increasing lodging resistance (LR) usually limits yield potential. Stem strength and leaf type are major traits related to LR and yield, respectively. Hence, understanding the genetic basis of stem strength and leaf type is of help to reduce lodging and increase yield in LR breeding. Here, we carried out an association analysis to identify quantitative trait locus (QTLs) affecting stem strength-related traits (internode length/IL, stem wall thickness/SWT, stem outer diameter/SOD, and stem inner diameter/SID) and leaf type-associated traits (Flag leaf length/FLL, Flag leaf angle/FLA, Flag leaf width/FLW, leaf-rolling/LFR and SPAD/Soil, and plant analyzer development) using a diverse panel of 550 accessions and evaluated over two years. Genome-wide association study (GWAS) using 4,076,837 high-quality single-nucleotide polymorphisms (SNPs) identified 89 QTLs for the nine traits. Next, through "gene-based association analysis, haplotype analysis, and functional annotation", the scope was narrowed down step by step. Finally, we identified 21 candidate genes in 9 important QTLs that included four reported genes (TUT1, OsCCC1, CFL1, and ACL-D), and seventeen novel candidate genes. Introgression of alleles, which are beneficial for both stem strength and leaf type, or pyramiding stem strength alleles and leaf type alleles, can be employed for LR breeding. All in all, the experimental data and the identified candidate genes in this study provide a useful reference for the genetic improvement of rice LR.