Combined genome-wide association study and epistasis analysis reveal multifaceted genetic architectures of plant height in Asian cultivated rice.

Plant height (PH) in rice (Oryza sativa) is an important trait for its adaptation and agricultural performance. Discovery of the semi-dwarf1 (SD1) mutation initiated the Green Revolution, boosting rice yield and fitness, but the underlying genetic regulation of PH in rice remains largely unknown. Here, we performed genome-wide association study (GWAS) and identified 12 non-repetitive QTL/genes regulating PH variation in 619 Asian cultivated rice accessions. One of these was an SD1 structural variant, not normally detected in standard GWAS analyses. Given the strong effect of SD1 on PH, we also divided 619 accessions into subgroups harbouring distinct SD1 haplotypes, and found a further 85 QTL/genes for PH, revealing genetic heterogeneity that may be missed by analysing a broad, diverse population. Moreover, we uncovered two epistatic interaction networks of PH-associated QTL/genes in the japonica (Geng)-dominant SD1NIP subgroup. In one of them, the hub QTL/gene qphSN1.4/GAMYB interacted with qphSN3.1/OsINO80, qphSN3.4/HD16/EL1, qphSN6.2/LOC_Os06g11130, and qphSN10.2/MADS56. Sequence variations in GAMYB and MADS56 were associated with their expression levels and PH variations, and MADS56 was shown to physically interact with MADS57 to coregulate expression of gibberellin (GA) metabolic genes OsGA2ox3 and Elongated Uppermost Internode1 (EUI1). Our study uncovered the multifaceted genetic architectures of rice PH, and provided novel and abundant genetic resources for breeding semi-dwarf rice and new candidates for further mechanistic studies on regulation of PH in rice.

Genetic architecture of seed glycerolipids in Asian cultivated rice.

Glycerolipids are essential for rice development and grain quality but its genetic regulation remains unknown. Here we report its genetic base using metabolite-based genome-wide association study and metabolite-based quantitative traits locus (QTL) analyses based on lipidomic profiles of seeds from 587 Asian cultivated rice accessions and 103 chromosomal segment substitution lines, respectively. We found that two genes encoding phosphatidylcholine (PC):diacylglycerol cholinephosphotransferase (OsLP1) and granule-bound starch synthase I (Waxy) contribute to variations in saturated triacylglycerol (TAG) and lyso-PC contents, respectively. We demonstrated that allelic variation in OsLP1 sequence between indica and japonica results in different enzymatic preference for substrate PC-16:0/16:0 and different saturated TAG levels. Further evidence demonstrated that OsLP1 also affects heading date, and that co-selection of OsLP1 and a flooding-tolerant QTL in Aus results in the abundance of saturated TAGs associated with flooding tolerance. Moreover, we revealed that the sequence polymorphisms in Waxy has pleiotropic effects on lyso-PC and amylose content. We proposed that rice seed glycerolipids have been unintentionally shaped during natural and artificial selection for adaptive or import seed quality traits. Collectively, our findings provide valuable genetic resources for rice improvement and evolutionary insights into seed glycerolipid variations in rice.

Genetic basis underlying tiller angle in rice (Oryza sativa L.) by genome-wide association study.

KEY MESSAGE: Novel alleles of two reported tiller angle genes and eleven candidate genes for rice tiller angle were identified by combining GWAS with transcriptomic, qRT-PCR and haplotype analysis. Rice tiller angle is a key agronomic trait determining rice grain yield. Several quantitative trait loci (QTLs) affecting rice tiller angle have been mapped in the past decades. Little is known about the genetic base of tiller angle in rice, because rice tiller angle is a complex polygenic trait. In this study, we performed genome-wide association study (GWAS) on tiller angle in rice using a population of 164 japonica varieties derived from the 3 K Rice Genomes Project (3 K RGP). We detected a total of 18 QTLs using 1135519 single-nucleotide polymorphisms (SNP) based on three GWAS models (GLM, FastLMM and FarmCPU). Among them, two identified QTLs, qTA8.3 and qTA8.4, overlapped with PAY1 and TIG1, respectively, and additional 16 QTLs were identified for the first time. Combined with haplotype and expression analyses, we further revealed that PAY1 harbors one non-synonymous variation at its coding region, likely leading to variable tiller angle in the population, and that nature variations in the promoter of TIG1 significantly affect its expression, closely correlating with tiller angle phenotypes observed. Similarly, using qRT-PCR and haplotype analysis, we identified 1 and 7 candidate genes in qTA6.1 and qTA8.1 that were commonly detected by two GWAS models, respectively. In addition, we identified 3 more candidate genes in the remaining 14 novel QTLs after filtering by transcriptome analysis and qRT-PCR. In summary, this study provides new insights into the genetic architecture of rice tiller angle and candidate genes for rice breeding.

Dissection of the Genetic Basis of Rice Panicle Architecture Using a Genome-wide Association Study.

Panicle architecture is one of the major factors influencing productivity of rice crops. The regulatory mechanisms underlying this complex trait are still unclear and genetic resources for rice breeders to improve panicle architecture are limited. Here, we have performed a genome-wide association study (GWAS) to analyze and identify genetic determinants underlying three panicle architecture traits. A population of 340 rice accessions from the 3000 Rice Genomes Project was phenotyped for panicle length, primary panicle number and secondary branch number over two years; GWAS was performed across the whole panel, and also across the japonica and indica sub-panels. A total of 153 quantitative trait loci (QTLs) were detected, of which 5 were associated with multiple traits, 8 were unique to either indica or japonica sub-panels, while 37 QTLs were stable across both years. Using haplotype and expression analysis, we reveal that genetic variations in the OsSPL18 promoter significantly affect gene expression and correlate with panicle length phenotypes. Three new candidate genes with putative roles in determining panicle length were also identified. Haplotype analysis of OsGRRP and LOC_Os03g03480 revealed high association with panicle length variation. Gene expression of DSM2, involved in abscisic acid biosynthesis, was up-regulated in long panicle accessions. Our results provide valuable information and resources for further unravelling the genetic basis determining rice panicle architecture. Identified candidate genes and molecular markers can be used in marker-assisted selection to improve rice panicle architecture through molecular breeding.

Gibberellins orchestrate panicle architecture mediated by DELLA-KNOX signalling in rice.

Panicle architecture is a key determinant of grain yield in cereals, but the mechanisms governing panicle morphogenesis and organ development remain elusive. Here, we have identified a quantitative trait locus (qPA1) associated with panicle architecture using chromosome segment substitution lines from parents Nipponbare and 9311. The panicle length, branch number and grain number of Nipponbare were significantly higher than CSSL-9. Through map-based cloning and complementation tests, we confirmed that qPA1 was identical to SD1 (Semi Dwarf1), which encodes a gibberellin 20-oxidase enzyme participating in gibberellic acid (GA) biosynthesis. Transcript analysis revealed that SD1 was widely expressed during early panicle development. Analysis of sd1/osga20ox2 and gnp1/ osga20ox1 single and double mutants revealed that the two paralogous enzymes have non-redundant functions during panicle development, likely due to differences in spatiotemporal expression; GNP1 expression under control of the SD1 promoter could rescue the sd1 phenotype. The DELLA protein SLR1, a component of the GA signalling pathway, accumulated more highly in sd1 plants. We have demonstrated that SLR1 physically interacts with the meristem identity class I KNOTTED1-LIKE HOMEOBOX (KNOX) protein OSH1 to repress OSH1-mediated activation of downstream genes related to panicle development, providing a mechanistic link between gibberellin and panicle architecture morphogenesis.